The true ground state and low-temperature phase diagram of PbZrO3: antiferroelectric or ferrielectric?

Lead zirconate (PbZrO3) is the archetypal antiferroelectric, yet first-principles calculations predict that a ferrielectric phase with an uncompensated up-up-down pattern of Pb displacements is lower in free energy than the canonical antipolar Pbam structure at low temperature, and that an 80-atom Pnam antiferroelectric phase should condense below roughly 115–140 K. The energy differences at stake are only a few meV per formula unit, comparable to the accuracy of the exchange-correlation functionals used to compute them. This problem asks for a settled determination of the equilibrium low-temperature phase diagram of PbZrO3 — the identity of the ground state, the transition temperatures, the experimental status of the predicted Pnam phase and of the controversial intermediate phase near 500 K — together with an effective-Hamiltonian or Ginzburg-Landau theory that captures the ferrielectric–antiferroelectric competition. Later work has realized the ferrielectric phase at room temperature in strained films and at translational boundaries, and machine-learned-potential simulations support the metastability of the observed antipolar phase, but as of 2026 the thermodynamic ground state and the effective theory remain unresolved.

Background

Lead zirconate, PbZrO3, is a perovskite oxide (ABO3 structure) and the archetype of antiferroelectricity: below about 500 K its Pb ions displace in an antiparallel up-up-down-down pattern along a pseudocubic axis, so that the 40-atom orthorhombic cell (space group $Pbam$, conventionally called AFE40) carries zero net polarization, yet an applied electric field switches the material reversibly into a polar phase, producing the double hysteresis loop that defines antiferroelectric (AFE) behavior and underpins energy-storage, electrocaloric, and electromechanical applications. A ferrielectric (FiE) is a structure in which antiparallel dipoles only partially compensate — for example an up-up-down arrangement — leaving a net polarization smaller than that of an ordinary ferroelectric (FE).

The accepted microscopic picture of AFE40 came from density-functional theory (DFT): the phase is driven not by the antipolar Pb displacements alone but by a trilinear coupling between antiphase rotations of the oxygen octahedra (the $R_4^+$ tilt mode at wave vector $\mathbf{q}=(\frac{1}{2},\frac{1}{2},\frac{1}{2})$), the antipolar $\Sigma_2$ mode at $\mathbf{q}=(\frac{1}{4},\frac{1}{4},0)$, and a small secondary distortion (Íñiguez et al. 2014). Effective Hamiltonians — atomistic models built from collective local-mode variables and parameterized against DFT, then simulated by Monte Carlo at finite temperature — were subsequently shown to reproduce this $Pbam$ ground state.

In 2021 two first-principles works unsettled this picture. Baker et al. found a dynamical instability of AFE40 that doubles the c axis and condenses into an 80-atom $Pnam$ structure (AFE80) of slightly lower energy. Aramberri, Cazorla, Stengel and Íñiguez instead proposed a ferrielectric phase: the simplest modification of AFE40 that keeps the $R_4^+$ tilts but rearranges the Pb displacements into the uncompensated up-up-down pattern, with a computed net polarization of about 0.11 C/m² (later work labels this phase $Ima2$). Within the LDA and PBEsol exchange-correlation functionals this FiE phase lies below AFE40 in free energy at low temperature, with a crossing estimated between roughly 255 and 460 K depending on the functional (about 370 K once thermal expansion is included in the quasi-harmonic approximation, which treats each phonon as a harmonic oscillator whose frequency shifts with the thermally expanded lattice volume). AFE80 is predicted, in most functionals, to appear as a phonon instability of AFE40 between 115 and 140 K. The difficulty is that all of these energy differences are only a few meV per formula unit — at or below the systematic differences between exchange-correlation functionals. Indeed PBE and SCAN, which the authors argue are unreliable for this compound, predict orderings incompatible with experiment (PBE even makes the FE phase the ground state up to 500 K).

The authors reconcile their calculations with the universal experimental observation of AFE40 by a kinetic argument: the AFE40→FiE transformation is first-order (the two structures have different order parameters and cannot be connected by a continuous symmetry breaking), with an estimated energy barrier of order 100 meV per formula unit, so a sample that crystallizes from the cubic phase near 500 K may simply never reach the FiE minimum on cooling. Six issues were left open, which together constitute this problem: (i) the true FiE–AFE40 transition temperature, and whether FiE is thermodynamically stable at ambient conditions; (ii) whether ferrielectric polymorphs other than the single, simplest one explored exist with even lower (free) energy; (iii) the actual AFE40→FiE transformation pathway and barrier, and whether FiE can be realized reproducibly in bulk (its realization in films is now established); (iv) the experimental status of AFE80 — missed experimentally, or suppressed by lattice quantum fluctuations (zero-point motion of the atoms beyond quasi-harmonic corrections, known to suppress low-temperature transitions in quantum paraelectrics such as SrTiO3); (v) the structure of the experimentally reported intermediate phase between cubic and AFE40, a narrow window of order 25 K near 500 K that has been controversial since the earliest studies; and (vi) whether existing effective-Hamiltonian or Ginzburg-Landau theories can capture the FiE state, which is argued to require symmetry-allowed couplings spanning four wave vectors — $\mathbf{q}=(0,0,0)$, $\mathbf{q}=(\frac{1}{3},\frac{1}{3},0)$, $\mathbf{q}=(\frac{1}{2},\frac{1}{2},\frac{1}{2})$ and $\mathbf{q}=(\frac{1}{6},\frac{1}{6},\frac{1}{2})$ — including the fourth-order invariant $Q(\Gamma_4^-)(Q(\Sigma_2))^3$ between the zone-center polar mode $\Gamma_4^-$ and the antipolar $\Sigma_2$ mode, an interaction absent in AFE40.

Problem Statement

Determine the equilibrium low-temperature phase behavior of bulk PbZrO3 and the effective theory that describes it — that is, decide whether the archetypal antiferroelectric is genuinely antiferroelectric. The competing phases are the canonical antipolar $Pbam$ structure (AFE40, observed in all experiments below about 500 K), the 80-atom antiferroelectric $Pnam$ structure (AFE80, predicted to condense from a phonon instability of AFE40 between 115 and 140 K), the ferrielectric $Ima2$ structure (FiE, computed lower in free energy than AFE40 at low temperature with a net polarization of about 0.11 C/m²), and any further ferrielectric polymorphs not yet explored. The free-energy differences at stake are a few meV per formula unit, at or below the accuracy of standard exchange-correlation functionals. The following linked questions remain open.

1. Ground state and phase diagram of bulk PbZrO3 at ambient pressure. Determine, with quantified accuracy sufficient to resolve few-meV-per-formula-unit differences, which phase is the thermodynamic ground state at low temperature; the transition temperature between the FiE and AFE40 phases (and hence whether FiE is thermodynamically stable at ambient conditions); whether the AFE40→AFE80 transition predicted at 115–140 K actually occurs or is shifted or suppressed by lattice quantum fluctuations beyond zero-point corrections; and whether ferrielectric polymorphs other than the single simplest one exist with still lower (free) energy.

2. Status of the disputed phases. Establish whether the AFE80 $Pnam$ phase occurs in its predicted low-temperature stability window (roughly below 115–140 K) in equilibrated bulk samples — or is instead suppressed by lattice quantum fluctuations — and resolve the structure of the narrow (order 25 K, near 500 K) intermediate phase between the cubic and AFE40 phases.

3. Effective theory and transformation pathway. Construct or revise an effective-Hamiltonian or Ginzburg-Landau theory for PbZrO3 — including the symmetry-allowed couplings spanning $\mathbf{q}=(0,0,0)$, $\mathbf{q}=(\frac{1}{3},\frac{1}{3},0)$, $\mathbf{q}=(\frac{1}{2},\frac{1}{2},\frac{1}{2})$ and $\mathbf{q}=(\frac{1}{6},\frac{1}{6},\frac{1}{2})$, among them the fourth-order invariant $Q(\Gamma_4^-)(Q(\Sigma_2))^3$ — that reproduces the first-principles energy landscape of the competing phases and predicts the FiE–AFE40 competition and the finite-temperature phase diagram; and use it or an equivalent simulation to compute the AFE40→FiE transformation pathway and barrier, thereby determining whether the ferrielectric phase is realizable in equilibrated bulk samples (its stabilization in strained films is already established and is not part of what remains open).

An answer may proceed by any sound route — converged first-principles free-energy calculations with demonstrated robustness across functionals or against higher-accuracy references, path-integral or anharmonic treatments where quantum fluctuations matter, validated machine-learned interatomic potentials, or definitive experiments on equilibrated bulk samples — but it must adjudicate thermodynamic stability, not merely kinetic accessibility, and must be consistent with the experimental record: the universal observation of AFE40 below 500 K, the room-temperature ferrielectric structures found at translational boundaries and in strained films, the weak low-temperature polarity reported in single crystals, and the double-hysteresis switching behavior.

References

1. H. Aramberri, C. Cazorla, M. Stengel, J. Íñiguez, 'On the possibility that PbZrO3 not be antiferroelectric', npj Computational Materials 7, 196 (2021). DOI: 10.1038/s41524-021-00671-w, https://doi.org/10.1038/s41524-021-00671-w (source of the open question)
2. J. S. Baker, M. Paściak, J. K. Shenton, P. Vales-Castro, B. Xu, J. Hlinka, P. Márton, R. G. Burkovsky, G. Catalan, A. M. Glazer, D. R. Bowler, 'A re-examination of antiferroelectric PbZrO3 and PbHfO3: an 80-atom Pnam structure', arXiv:2102.08856 (2021). https://arxiv.org/abs/2102.08856
3. J. Íñiguez, M. Stengel, S. Prosandeev, L. Bellaiche, 'First-principles study of the multimode antiferroelectric transition in PbZrO3', Physical Review B 90, 220103(R) (2014). DOI: 10.1103/PhysRevB.90.220103, https://doi.org/10.1103/PhysRevB.90.220103
4. Y. Liu, R. Niu, A. Majchrowski, K. Roleder, K. Cordero-Edwards, J. M. Cairney, J. Arbiol, G. Catalán, 'Translational Boundaries as Incipient Ferrielectric Domains in Antiferroelectric PbZrO3', Physical Review Letters 130, 216801 (2023). DOI: 10.1103/PhysRevLett.130.216801, https://doi.org/10.1103/PhysRevLett.130.216801
5. Y. Yao, A. Naden, M. Tian, S. Lisenkov, Z. Beller, A. Kumar, J. Kacher, I. Ponomareva, N. Bassiri-Gharb, 'Ferrielectricity in the Archetypal Antiferroelectric, PbZrO3', Advanced Materials 35, 2206541 (2023). DOI: 10.1002/adma.202206541, https://doi.org/10.1002/adma.202206541
6. Z. Yu, N. Fan, Z. Fu, B. He, S. Yan, H. Cai, X. Chen, L. Zhang, Y. Zhang, B. Xu, G. Wang, F. Xu, 'Room-temperature stabilizing strongly competing ferrielectric and antiferroelectric phases in PbZrO3 by strain-mediated phase separation', Nature Communications 15, 3438 (2024). DOI: 10.1038/s41467-024-47776-6, https://doi.org/10.1038/s41467-024-47776-6
7. Y. Liu, R. Niu, R. Uriach, D. Pesquera, J. C. Roque, J. Santiso, J. Cairney, X. Liao, J. Arbiol, G. Catalán, 'Coexistence of ferroelectric and ferrielectric phases in ultrathin antiferroelectric PbZrO3 thin films', Microstructures 4, 2024012 (2024). DOI: 10.20517/microstructures.2024.12, https://doi.org/10.20517/microstructures.2024.12
8. H. Zhang, H.-C. Thong, L. Bastogne, C. Gui, X. He, P. Ghosez, 'Finite-temperature properties of the antiferroelectric perovskite PbZrO3 from a deep-learning interatomic potential', Physical Review B 110, 054109 (2024). DOI: 10.1103/PhysRevB.110.054109, https://doi.org/10.1103/PhysRevB.110.054109 (preprint: arXiv:2406.09011)
9. K. Patel, B. Xu, S. Prosandeev, R. Faye, B. Dkhil, P. Janolin, L. Bellaiche, 'Temperature-dependent properties of the antiferroelectric model PbZrO3: An effective Hamiltonian study', Physical Review B 106, 214108 (2022). DOI: 10.1103/PhysRevB.106.214108, https://doi.org/10.1103/PhysRevB.106.214108
10. K. Roleder, G. Catalan, A. M. Glazer, J. S. Baker, J.-H. Ko, F. H. Naqvi, S. B. Junaid, A. Majchrowski, Z. Trybuła, J. Zaręba, I. Lazar, D. Kajewski, J. Koperski, A. Soszyński, 'Weak low-temperature polarity in a PbZrO3 single crystal', Physical Review B 107, L140102 (2023). DOI: 10.1103/PhysRevB.107.L140102, https://doi.org/10.1103/PhysRevB.107.L140102
11. N. Maity, M. Haddad, N. Bassiri-Gharb, A. Kumar, L. Jones, S. Lisenkov, I. Ponomareva, 'Ferroelectricity at the extreme thickness limit in the archetypal antiferroelectric PbZrO3', npj Computational Materials 11, 48 (2025). DOI: 10.1038/s41524-025-01520-w, https://doi.org/10.1038/s41524-025-01520-w
12. K. Parmar, P. Dufour, E. Texier, C. Carrétéro, X. Li, F. Godel, J. Hlinka, B. Dkhil, D. Sando, H. Aramberri, J. Íñiguez-González, S. Fusil, A. Gloter, T. Maroutian, V. Garcia, 'Establishing a pure antiferroelectric PbZrO3 phase through tensile epitaxial strain', Nature Communications 16, 6536 (2025). DOI: 10.1038/s41467-025-61867-y, https://doi.org/10.1038/s41467-025-61867-y
13. M. Kniazeva, A. Ganzha, R. Gao, A. Dasgupta, A. V. Filimonov, R. Burkovsky, 'Revealing the phase transition scenario in antiferroelectric thin films by x-ray diffuse scattering', Journal of Applied Physics 135, 244102 (2024). DOI: 10.1063/5.0184724, https://doi.org/10.1063/5.0184724
14. G. Catalan, A. Gruverman, J. Íñiguez-González, D. Meier, M. Trassin, 'A modern perspective on antiferroelectrics', Nature Materials 25, 557–565 (2026). DOI: 10.1038/s41563-026-02483-z, https://doi.org/10.1038/s41563-026-02483-z
15. T. Zhang, Y. Si, 'Unresolved controversies in perovskite-based antiferroelectrics: Fundamentals and frontiers', Materials Science and Engineering: R: Reports 167, 101121 (2026). DOI: 10.1016/j.mser.2025.101121, https://doi.org/10.1016/j.mser.2025.101121
16. N. A. Arkhipov, P. D. Shulga, R. G. Burkovsky, 'Gibbs energy approach to the thermodynamical phase diagrams of PbZrO3 and PbHfO3', Physical Review B 111, L180102 (2025). DOI: 10.1103/PhysRevB.111.L180102, https://doi.org/10.1103/PhysRevB.111.L180102