{"id":"567d8ec7-340c-421e-b54b-e2c8152ce674","arxiv_id":"2607.27752","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"An extra oxygen tilt in antiferroelectric PbZrO3 breaks the symmetry that cancels opposite dipoles, creating a predicted Pmc2_1 ferrielectric phase with net polarization.","lead":"By adding a particular oxygen octahedral tilt to antiferroelectric lead zirconate, the authors show the material can develop a ferrielectric state with small net electric polarization. The result combines symmetry analysis, density-functional calculations, and atom-resolved electron microscopy, offering a general route to ferrielectricity in antiferroelectrics.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite-temperature stabilization of Pmc2_1 over Ima2 rests on QHA phonon free energies with meV-scale errors; a small correction could reverse the ordering and weaken the claim of accessible ferrielectricity.","rationale":"The paper's symmetry argument is well posed and the DFT relaxations at 10 GPa show a genuinely polar Pmc2_1 phase with a plausible microscopic mechanism (M2+ tilt coupling to polarization). The experimental images provide suggestive but not definitive structural evidence. The weakest link is the finite-temperature ordering: the assertion that Pmc2_1 becomes thermodynamically preferred over Ima2 near room temperature is based on QHA phonon calculations where the dynamic instability of Pbam is handled by excluding imaginary modes, and where the relevant energy scales are only a few meV/f.u. This is precisely the kind of situation where small errors in phonon frequencies or anharmonic contributions can flip the relative stability of two competing phases. The reader's verdict of CONDITIONAL is appropriate; my concern aligns with the reader's primary weakest assumption, and the concrete test I propose would directly test whether the finite-temperature ordering is robust. I see no reason to change the verdict.","tokens_in":12609,"tokens_out":3584,"duration_ms":39681,"concrete_test":"Recompute the Gibbs free-energy difference between Pmc2_1 and Ima2 at 10 GPa and 300 K using an anharmonic method (e.g., self-consistent harmonic approximation or molecular dynamics with a machine-learned interatomic potential). If Ima2 remains lower than Pmc2_1, the central claim of room-temperature thermodynamic stabilization is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Pmc2_1 becomes thermodynamically more favorable than Ima2 near room temperature (Fig. 3b) rests on quasi-harmonic phonon free energies computed with the imaginary modes of Pbam excluded, and on DFT enthalpy differences of only ~3.6 meV/f.u. between competing phases. QHA errors from neglected anharmonicity and from the approximate treatment of dynamically unstable modes are easily of this magnitude; a phonon-DOS error of a few meV/f.u. at 300 K would reorder the Gibbs free energies. Since the experimental evidence is 'Pmc2_1-like' rather than a full structural refinement, the finite-temperature phase ordering is the main pillar connecting the symmetry mechanism to realizability. If the ordering reverses, the phase may still exist as a metastable state, but the claim that it 'gains increasing free-energy advantage' and is 'thermodynamically more favorable around room temperature' fails.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the antiferroelectric Pbam phase of PbZrO3 can be converted into a polar Pmc2_1 ferrielectric phase by condensing an additional in-phase octahedral tilt (M2+ mode). Group-theoretic analysis identifies symmetry-allowed couplings by which the tilt, together with the pre-existing Pbam distortions, induces the polar Gamma4- mode through quadratic-linear-linear and quadrilinear terms. First-principles DFT relaxations show that Pmc2_1 becomes stable under hydrostatic pressure above about 5 GPa and under biaxial compressive strain, and that at 10 GPa it is dynamically stable and 3.6 meV/f.u. lower in enthalpy than Pbam. Quasi-harmonic calculations are used to argue that Pmc2_1 overtakes the Ima2 ferrielectric phase near room temperature. STEM-HAADF and iDPC imaging are presented as evidence of Pmc2_1-like structures in thin films and single crystals. The paper also extends the mechanism to PbSnO3 and to 180-degree domain walls in NaNbO3.","tokens_in":1458,"tokens_out":1601,"duration_ms":48638,"significance":"If the central mechanism is correct, the paper identifies a new symmetry-based route to ferrielectricity that is qualitatively different from previously studied FiE states: the Pmc2_1 phase is connected to the parent Pbam phase via a direct group-subgroup relation and can be reached without reversing any Pb displacements. The competitive improper coupling picture, in which polarization is determined by a subtle balance of M2+ and S1 modes, is a valuable conceptual extension of hybrid-improper ferroelectricity, and the predicted polarization reversal within a fixed phase by tuning nonpolar mode amplitudes is a striking and falsifiable consequence. The paper also provides direct DFT relaxations, phonon stability checks, and a transparent Landau fit. If the finite-temperature and experimental evidence were as robust as the symmetry argument, this would be a significant advance. As it stands, the theoretical mechanism is credible but the paper overreaches in its kinetic and experimental claims.","major_comments":[{"comment":"The finite-temperature ordering of Pmc2_1 above Ima2 rests on QHA phonon free energies with imaginary modes of Pbam excluded and enthalpy differences of only a few meV/f.u. (3.6 meV/f.u. at 10 GPa). QHA errors from anharmonicity and from the treatment of dynamically unstable phases are easily of this magnitude; a small correction could reverse the room-temperature ordering. The paper acknowledges the approximation but still concludes that Pmc2_1 becomes thermodynamically more favorable around room temperature. This is load-bearing for the claim of accessible FiE. Please provide a quantitative error estimate, cross-check with another method (e.g., the deep-learning potential cited as Ref. [43]), or soften the thermodynamic claim to 'competitive metastable state'.","section":"Fig. 3(b) and QHA discussion"},{"comment":"The paper asserts that the Pbam-to-Pmc2_1 transition is 'barrier-free' and 'kinetically easily accessible' based on the direct group-subgroup relation and shared fourfold periodicity. A group-subgroup relation does not guarantee a barrier-free path; continuous symmetry-allowed transitions can still be first-order with a finite nucleation barrier. No nudged-elastic-band calculation, Landau energy profile along the transition path, or explicit kinetic model is provided. Since the argument that Pmc2_1 can compete with the Ima2 ground state relies on this kinetic accessibility, please either compute the barrier or replace 'barrier-free' with a more cautious statement.","section":"Fig. 3(a) inset and text on kinetic accessibility"},{"comment":"The experimental identification of Pmc2_1 is based on asymmetric oxygen-chain rippling measured from an iDPC image near a chemically intermixed PZO/SRO interface, and from regions with Pb vacancies in the single crystal. The rippling angles are reported without error bars, and no simulated iDPC images of alternative local distortions (e.g., Pbam with inhomogeneous strain, defect-induced tilt patterns, or interface polarity effects) are compared. The sentence 'provides definitive structural evidence supporting the presence of the Pmc2_1-like FiE state' is stronger than the data support. Please quantify the rippling amplitude with uncertainties, and test against plausible alternative structural models in the same local environment.","section":"Fig. 4(d-g) and iDPC evidence"}],"minor_comments":[{"comment":"Some arrow symbols in the abstract and body appear garbled (e.g., the polarization pattern symbols). Please ensure the typesetting of the 'up-up-down-down' pattern is consistent and correct.","section":"Throughout"},{"comment":"The notation 'Sym.δ Pb' and 'Asym.∆ O' is clear in the figure, but the main text would benefit from a one-sentence definition of δ and ∆ as displacement and ripple variables, respectively.","section":"Introduction and Fig. 1"},{"comment":"Equations (1)-(3) give the coupling terms symbolically, but the actual Landau polynomial is only in the Supplementary Material. A sample free-energy expression with fitted coefficient values in the main text would help readers assess the competition described in Fig. 2.","section":"Section on Landau model"},{"comment":"The label 'QΓ4−=0' is not defined in the caption; specify that it is the order-parameter amplitude of the polar mode.","section":"Fig. 2(b)"}],"recommendation":"major_revision","confidential_remarks":"The symmetry mechanism and DFT results are solid and likely correct, but the paper's central claims of finite-temperature stabilization and direct confirmation are overstated relative to the evidence. The QHA issue and the lack of a real barrier calculation are fixable but nontrivial; the experimental section needs more rigorous quantification. I would be willing to see a revised manuscript with those points addressed. The paper is potentially important and should not be rejected on the basis of the current overreach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's core claim is new and convincing: adding an M2+ in-phase tilt to the Pbam phase of PbZrO3 breaks the symmetry constraint on the antiparallel Pb displacements, yielding an uncompensated ↑↑↓↓ polar phase (Pmc2_1) that is reached without dipole reversal. That is a clean symmetry argument, and the DFT relaxations showing stabilization above ~5 GPa, the phonon stability at 10 GPa, and the Landau analysis with the dual induction pathways are all solid. The Landau coefficients are honestly fitted to DFT data and the phase stability itself comes from direct relaxation, so there is no circularity problem.\n\nThe soft spots are exactly where the reader's report and stress-test point. The QHA comparison in Fig. 3(b) excludes the imaginary modes of Pbam and rests on enthalpy differences of a few meV/f.u., so the claim that Pmc2_1 becomes more favorable than Ima2 near room temperature is not robust to systematic errors of that size. The paper even acknowledges the approximation, and the trend is plausible, but 'thermodynamically more favorable around room temperature' is stated too strongly. Similarly, the experimental evidence is 'Pmc2_1-like' — asymmetric oxygen-chain rippling without error bars or a full refinement, in regions near an intermixed PZO/SRO interface or with Pb vacancies. That is suggestive but not definitive, and the paper should say so more plainly. The 'barrier-free' kinetic picture is also inferred from the group-subgroup relation, not computed.\n\nNone of these kill the central mechanism. The symmetry argument and the direct DFT result stand on their own, and the generalization to PbSnO3 and NaNbO3 domain walls is a useful extension. The paper is for anyone working on antiferroelectrics, improper ferroelectricity, or perovskite phase competition. It deserves a serious referee: the novel phase and mechanism warrant careful vetting, and the finite-temperature and experimental claims can be revised or more tightly framed.","headline":"A genuinely new symmetry route to ferrielectricity in PbZrO3, with solid DFT support, but the finite-temperature ordering and the STEM identification are weaker than the central mechanism.","tokens_in":13383,"tokens_out":2259,"would_cite":true,"duration_ms":24658,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"An added oxygen tilt can turn antiferroelectric PbZrO3 into a polar ferrielectric phase.","keywords":["ferrielectricity","antiferroelectricity","PbZrO3","octahedral tilts","improper ferroelectricity","first-principles","STEM imaging","perovskite"],"falsifier":"A direct test would be to measure the oxygen-chain rippling and lead displacements in a clean, defect-poor PbZrO3 sample under controlled hydrostatic pressure around 5 GPa; if no fourfold-periodic uncompensated 'up-up-down-down' pattern with asymmetric ripples appears, the central claim fails. More cheaply, recompute the 10 GPa Gibbs free-energy difference between Pmc2_1 and Ima2 with a phonon method that does not discard imaginary modes; if Ima2 remains lower through 500 K, the room-temperature stabilization argument collapses.","tokens_in":12416,"feed_emoji":"⚛️","tokens_out":4178,"duration_ms":39874,"temperature":0.7,"pith_summary":"This paper argues that a hidden polar phase of PbZrO3 is created by adding an oxygen octahedral tilt to its known antiferroelectric structure. The tilt breaks a symmetry constraint that keeps the two antiparallel lead displacements equal in size, so the compensated 'up-up-down-down' pattern becomes an uncompensated 'up-up-down-down' pattern with a net polarization. First-principles calculations show this Pmc2_1 phase is dynamically stable under lattice contraction and gains free-energy advantage near room temperature. Atomic-resolution imaging in thin films and single crystals finds the asymmetric oxygen-chain rippling and uncompensated lead displacements that mark this phase. A sympathetic reader would care because it offers a kinetically easy, symmetry-governed route to ferrielectricity and a new kind of competitive improper ferroelectricity.","feed_headline":"Extra oxygen tilt makes PbZrO3 ferrielectric","feed_subtitle":"A symmetry-lifting tilt turns compensated antiparallel dipoles into a polar phase with net polarization.","key_machinery":"The M2+ in-phase octahedral tilt mode is the trigger. Condensing M2+ into the Pbam structure breaks inversion-related equivalence; it couples with inherited Pbam modes (Sigma2 and S2) to activate the S1 modulated-tilt mode and the Gamma4- polar mode via quadratic-linear-linear and quadrilinear couplings. The net polarization is set by the competition between M2+ and S1, which can push the polarization continuously through zero and reverse its sign without reversing the driving modes themselves.","core_discovery":"The central claim is that Pbam antiferroelectric PbZrO3 can be transformed into a polar Pmc2_1 ferrielectric phase by condensing an additional M2+ in-phase oxygen octahedral tilt. In Pbam, the 'up-up-down-down' lead displacements are related by a symmetry operation that forces the two antiparallel sublattices to have equal magnitudes, so the structure is nonpolar. The added tilt breaks that operation and makes the oxygen-chain ripples between the 'up-up' and 'down-down' pairs unequal, which removes the equality constraint and lets one antiparallel sublattice become larger. The result is an uncompensated 'up-up-down-down' pattern with a net polarization that belongs to space group Pmc2_1. The","pith_inferences":["If local Pb vacancies act as microscopic sources of lattice contraction, defect patterning could be used to imprint or erase Pmc2_1-like polar regions at the nanoscale in otherwise antiferroelectric films.","The ability to reverse polarization sign by tuning M2+/S1 amplitudes without reversing the driving modes suggests a continuous polarization-reversal route that might be observable as soft-mode behavior or as domains annihilating at zero net polarization.","The 'competitive' improper mechanism raises the possibility of multiple near-degenerate polar states with different polarization magnitudes; field cycling might show staircase-like or non-binary switching in energy-storage devices.","The same symmetry analysis could be applied to other Pbam antiferroelectrics such as PbHfO3 to predict pressure-induced ferrielectricity."],"forward_implications":["Above about 5 GPa hydrostatic pressure (approximately -1.2% linear strain), Pmc2_1 sits lower in enthalpy than Pbam and becomes dynamically stable; it remains a metastable competitor above 10 GPa.","The Pbam-to-Pmc2_1 path is group-subgroup and barrier-free in the dipole pattern, so the polar phase is kinetically reachable by adjusting dipole magnitudes rather than reversing dipoles, unlike the threefold Ima2 ferrielectric ground state.","At 10 GPa and near room temperature, QHA free-energy results put Pmc2_1 thermodynamically above Ima2, suggesting a heating path from Pmc2_1 directly to cubic without returning to Pbam.","Reciprocal couplings imply electric fields can excite M-point octahedral tilts, giving a way to control tilt order electrically.","Because the mechanism only needs polarization modulation plus in-phase tilts, it extends beyond PZO, as shown for PbSnO3 and 180-degree domain walls in the ferroelectric Q phase of NaNbO3."],"fun_headline_variants":["Oxygen tilt gives PbZrO3 a polar twist","Symmetry lift turns PbZrO3 ferrielectric","New tilt unbalances PbZrO3 dipoles","Tilt-induced polar phase in PbZrO3"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The finite-temperature conclusion depends on quasi-harmonic phonon free energies in which imaginary modes of Pbam are discarded, and the energy differences separating Pmc2_1 from Ima2 are only a few meV per formula unit; an error of that size could reverse the ordering in Fig. 3(b). The STEM evidence also relies on regions near a chemically intermixed interface or containing Pb vacancies, so alternative local distortions are not excluded.","fun_headline_variants_meta":{"raw":{"variants":["Oxygen tilt gives PbZrO3 a polar twist","Symmetry lift turns PbZrO3 ferrielectric","New tilt unbalances PbZrO3 dipoles","Tilt-induced polar phase in PbZrO3"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000699,"raw_usage":{"total_tokens":2982,"prompt_tokens":721,"completion_tokens":2261,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":465,"completion_tokens_details":{"reasoning_tokens":2205}},"tokens_in":465,"tokens_out":2261,"duration_ms":15887,"temperature":1.0,"reasoning_tokens":2205,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T01:55:52.059692+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to measure the oxygen-chain rippling and lead displacements in a clean, defect-poor PbZrO3 sample under controlled hydrostatic pressure around 5 GPa; if no fourfold-periodic uncompensated 'up-up-down-down' pattern with asymmetric ripples appears, the central claim fails. More cheaply, recompute the 10 GPa Gibbs free-energy difference between Pmc2_1 and Ima2 with a phonon method that does not discard imaginary modes; if Ima2 remains lower through 500 K, the room-temperature stabilization argument collapses.","supporting_citations":[],"review_version":1}