REVIEW 2 major objections 5 minor 1 cited by
Ultrathin hafnia's out-of-plane lattice expansion is a sign of suppressed polarization, not enhanced ferroelectricity.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-04 16:28 UTC pith:LVVYC26C
load-bearing objection A mostly convincing mechanistic explanation for the reverse size effect in hafnia, with a clean orientation-hyperferroelectricity argument; the main soft spot is the experimental comparison resting on fitted surface parameters from a free-standing slab. the 2 major comments →
Origin of Anomalous Size Effects in Ferroelectric Hafnia Thin Films
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For (111)-oriented HfO2 slabs under open-circuit conditions, the out-of-plane polar displacement drops to roughly one-third of its bulk value, yet the out-of-plane interplanar spacing expands by about 0.03 Å. This inverse correlation is the paper's central discovery: the expansion is driven by a negative effective longitudinal piezoelectric coefficient (computed as d33^eff = -1.02 pm/V) responding to the depolarization field, plus a thickness-dependent surface stress described by a surface elastic model with parameters fitted from first-principles slabs. The combined model gives η_OP(t) = -ν η_IP(t) + η_OP(E_d), matching experimental d111 data for Hf0.5Zr0.5O2 films. Separately, a free-energ
What carries the argument
The argument rests on two quantitative tools. First, a surface elastic model (total energy = bulk strain energy plus surface energy with surface stress σ and surface elastic constant S) yields the thickness-dependent in-plane strain η_IP(t) = -2σ/(Ct+2S), which enters the out-of-plane expansion through the Poisson effect. Second, a phenomenological free-energy model for a ferroelectric under open-circuit conditions adds a depolarization penalty proportional to P_s^2/[1+χ_e]^2, showing that the (111) orientation's reduced out-of-plane polarization component lowers this penalty enough to keep the ferroelectric state stable. The negative piezoelectric coefficient of HfO2, computed as d33^eff =
Load-bearing premise
The quantitative reproduction of the experimental d111 thickness dependence rests on surface-stress parameters fitted from a single four-layer slab and on treating a HfO2/graphene capacitor as representative of the partial screening in the actual Hf0.5Zr0.5O2 films; if those fitted parameters or that screening proxy are not transferable, the quantitative claim is not established.
What would settle it
Measure, on the same ultrathin hafnia film, both the out-of-plane lattice spacing and the switchable polarization as functions of thickness: if d111 increases while polarization also increases, or if d111 decreases when the depolarization field is screened out by improved electrodes, the paper's central claim would be contradicted.
If this is right
- If the model is correct, the measured increase in d111 with decreasing film thickness is direct evidence of a residual depolarization field and negative piezoelectricity, not of enhanced polarization.
- Orientation engineering along (111) provides a general strategy to eliminate the critical thickness limit in conventional perovskite ferroelectrics such as PbTiO3, extending beyond hafnia.
- Ultrathin capacitors using Ca2N electrodes should restore near-bulk polarization, avoiding the depolarization-driven suppression seen with weaker screening.
- The reverse size effect should vanish in perfectly screened capacitors, where only the surface-stress-induced Poisson expansion remains.
- For (001)-oriented hafnia, open-circuit conditions stabilize a nonpolar ground state, affirming the orientation-specific nature of the mechanism.
Where Pith is reading between the lines
- A direct experimental test could vary the screening environment (e.g., electrode material or capping layer) on identical ultrathin hafnia films; the model predicts that improved screening should decrease d111 while increasing switchable polarization.
- The orientation-induced hyperferroelectricity mechanism may apply to other fluorite-structure or non-perovskite ferroelectrics where the polarization vector has a component along the surface normal smaller than the total magnitude.
- The surface-stress parameters fitted from a single slab could be checked for transferability by computing them at several thicknesses and seeing whether Eq. 2 still reproduces the DFT strains; such a check would test the model's predictive power.
- If Ca2N behaves as calculated, it could serve as a practical electrode for nanoscale ferroelectric memory, potentially avoiding the dead-layer problem that plagues metal electrodes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper addresses the anomalous out-of-plane lattice expansion observed in ultrathin (111)-oriented HfO2 films. Using DFT on HfO2(111) slabs and capacitors with different electrodes, the authors argue that the expansion is not a signature of enhanced ferroelectricity but coexists with suppressed polar displacements. They identify two mechanisms: a negative longitudinal piezoelectric response to the residual depolarization field and a thickness-dependent surface stress. A surface elasticity model with parameters fitted to a DFT slab is used to compute the in-plane strain, and the total out-of-plane strain is decomposed into a Poisson contribution and a depolarization-field contribution. This model is claimed to quantitatively reproduce the experimental d111(t) data for Hf0.5Zr0.5O2 films. The paper also predicts orientation-induced hyperferroelectricity in HfO2(111) and PbTiO3(111) under open-circuit conditions and identifies Ca2N as a near-ideal electrode.
Significance. If the central claim holds, the paper offers a resolution to a long-standing puzzle in hafnia ferroelectrics and proposes a practical design principle: choosing a film orientation with a reduced out-of-plane polarization component can bypass the critical-thickness limit. The internal consistency of the DFT evidence is a notable strength. The DFPT value deff_33 = −1.02 pm/V agrees with the slope extracted from capacitor calculations (−1.28 pm/V), and the surface-elasticity model fitted at one slab thickness reproduces DFT free-standing slab strains at other thicknesses. This shows the model is not merely an interpolation. The LGD treatment provides a transparent mechanism for the orientation dependence. The main caveat, however, is that the quantitative match to experiments rests on unverified transferability assumptions: the fitted surface parameters come from a free-standing symmetric slab, and the partial-screening curve is obtained by imposing that slab's in-plane strain on a different capacitor geometry. The experimental comparison is also made to Hf0.5Zr0.5O2 films on SrTiO3, not to pure HfO2, and no HZO or supported-film calculation is provided.
major comments (2)
- [Surface elasticity model, Eqs. (2)-(3) and Fig. 3] The central quantitative claim in Fig. 3b rests on the in-plane strain η_IP(t) from Eq. (2), with σ and S fitted from a single four-layer free-standing OCBC HfO2(111) slab (Fig. S5). The model is then validated for other free-standing slab thicknesses (Fig. 3a). However, the partial-screening curve in Fig. 3b is computed by 'imposing the same in-plane strain' from the free-standing slab on a HfO2/graphene capacitor. This is an assumption, not a calculation. The equilibrium in-plane strain of the graphene-capped capacitor is governed by two HfO2/graphene interfaces, not two free surfaces, and no DFT minimization of η_IP(t) for the capacitor geometry is reported. Moreover, the experimental films are supported on SrTiO3(001), where substrate clamping and interfacial stress should alter the in-plane strain. Since the first term of Eq. (3), −ν η_IP(t), is linear in η_IP(t), any transferabilit
- [Fig. 3b and comparison to Refs. [39,40]] The experimental d111 data are for Hf0.5Zr0.5O2 films, while all DFT slabs and the LGD model are for pure HfO2. The partial-screening curve corresponds to a HfO2/graphene capacitor with a specific Ed = −10.4 MV/cm, which is not derived from the actual electrode/interface conditions of the measured films. Without an explicit calculation for Zr-alloyed HfO2, or at least a demonstration that alloying leaves the negative piezoelectric coefficient and the effective surface-stress parameters unchanged, the match in Fig. 3b is suggestive rather than a quantitative validation for the measured system. The manuscript should state this limitation clearly and, ideally, provide bounds on how composition and screening variations shift the predicted curve.
minor comments (5)
- [Title and abstract] The arXiv title reads 'Origin of Anomalous Size Effects', while the main-text title reads 'Origin of Reverse Size Effect'. Please harmonize the terminology.
- [Eq. (1) and reference [21]] Reference [21] (Lv et al., on perovskite nanocrystals) does not appear related to surface elasticity; the surface-stress model is better supported by Ref. [22]. Please check and correct this citation.
- [Eq. (2), elastic constant C] The value C = 4.44 eV/Å^3 is introduced without specifying how it is derived (bulk elastic constant? from the same DFPT calculation?). Please state the source.
- [Fig. 3b] The experimental data points from Refs. [39,40] are shown without error bars or film-thickness uncertainty. Adding these would help assess how 'quantitative' the agreement really is.
- [Fig. 4 and PbTiO3 slab] The generalization to PbTiO3 rests on one DFT slab calculation (Fig. S6). A brief statement about thickness dependence or the implications of the slab termination would strengthen the claim.
Circularity Check
No significant circularity: central claims rest on independent DFT calculations and are compared to external experimental data without fitting those data.
full rationale
The paper's derivation chain is self-contained. The two mechanisms proposed for the anomalous out-of-plane expansion—negative longitudinal piezoelectric response and surface stress—are each quantified from first principles: the effective piezoelectric coefficient comes from DFPT (d_eff^33 = -1.02 pm/V) and is cross-checked against the slope of the DFT-computed d_111-E_d relation; the surface stress parameters σ and S are fitted to one four-layer open-circuit slab's energy-strain curve and then used in the analytic surface-elasticity model. This is a parameterization of DFT data, not a fit to experimental data. The resulting thickness-dependent η_IP(t) is validated against DFT-computed strains at other thicknesses (Fig. 3a); the t=4 point is in-sample, but the thickness trend is not forced by the fit. The final quantitative comparison with experimental d_111 values in Fig. 3b uses DFT calculations of free-standing HfO2(111) slabs, HfO2/graphene capacitors, and bulk HfO2 with imposed in-plane strains—none of these curves are fitted to the experimental d_111 points. The LGD hyperferroelectricity analysis uses independently computed total energies, polarizations, and susceptibilities, and the PbTiO3(111) prediction is confirmed by direct DFT slab relaxation. Self-citations (e.g., Refs. [19], [21], [40]) are not load-bearing: the claims they support are either independently recomputed in this paper or corroborated by external references. The transferability of free-standing-slab surface-stress parameters to films on SrTiO3 substrates is a modeling assumption that affects predictive robustness, but it does not make any equation reduce to its own input by construction. Therefore no specific circular step can be exhibited.
Axiom & Free-Parameter Ledger
free parameters (2)
- surface stress sigma =
0.38 eV/A^2
- surface elastic constant S =
-4.78 eV/A^2
axioms (5)
- domain assumption DFT with PBE/PBEsol captures the relevant polar distortions, piezoelectric coefficients, and surface energetics of HfO2 and PbTiO3.
- domain assumption The surface elasticity model of Eq. 1 with bulk elastic energy and two-surface energy is valid for films down to 0.9 nm.
- domain assumption Bulk LGD quantities U(lambda), Ps(lambda), and chi_e(lambda) can be applied to thin films in Eq. 5.
- domain assumption For (111) orientation, only the [111] component of Ps contributes to the depolarization energy penalty.
- domain assumption Hf0.5Zr0.5O2 experimental films can be compared to HfO2 model results without composition effects altering the mechanisms.
read the original abstract
The persistence of ferroelectricity in ultrathin HfO$_2$ films challenges conventional theories, particularly given the paradoxical observation that the out-of-plane lattice spacing increases as the film thickness decreases, an anomalous size effect absent in perovskite ferroelectrics. Here, we resolve this puzzle by revealing that this lattice expansion is counterintuitively coupled to a suppressed out-of-plane polarization. First-principles calculations combined with analytical modeling identify two mechanisms behind this expansion: a negative longitudinal piezoelectric response to the residual depolarization field and a positive surface stress that becomes significant at reduced thickness. Their interplay quantitatively reproduces the experimentally observed lattice expansion. Furthermore, (111)-oriented HfO$_2$ films can support out-of-plane polarization even under open-circuit conditions, in contrast to (001) films that stabilize a nonpolar ground state. This behavior points to the emergence of orientation-induced hyperferroelectricity, an unrecognized mechanism that enables polarization persistence through orientation engineering without electrode screening. We further demonstrate that this principle generalizes to conventional perovskites such as PbTiO$_3$, offering a strategy to eliminate the critical thickness limit by choosing the appropriate film orientation. As a practical pathway to device integration, we also identify the two-dimensional electride Ca$_2$N as a near-ideal electrode that fully restores the ferroelectric properties of HfO$_2$ in ultrathin capacitors.
Figures
Forward citations
Cited by 1 Pith paper
-
Atomic-scale theory of robust out-of-plane ferroelectricity in ultrathin films
Atomic-scale model attributes persistent ultrathin ferroelectricity to self-polarization and switchable termination layers correlated with characteristic structure, plus electrode stabilization effects.
Reference graph
Works this paper leans on
-
[1]
S. S. Cheema, D. Kwon, N. Shanker, R. dos Reis, S.-L. Hsu, J. Xiao, H. Zhang, R. Wagner, A. Datar, M. R. McCarter, C. R. Serrao, A. K. Yadav, G. Karbasian, C.-H. Hsu, A. J. Tan, L.-C. Wang, V. Thakare, X. Zhang, A. Mehta, E. Karapetrova, R. V. Chopdekar, P. Shafer, E. Arenholz, C. Hu, R. Proksch, R. Ramesh, J. Ciston, and S. Salahuddin, Enhanced ferroelec...
2020
-
[2]
S. S. Cheema, N. Shanker, S.-L. Hsu, Y. Rho, C.-H. Hsu, V. A. Stoica, Z. Zhang, J. W. Free- land, P. Shafer, C. P. Grigoropoulos, J. Ciston, and S. Salahuddin, Emergent ferroelectricity in subnanometer binary oxide films on silicon, Science376, 648 (2022)
2022
-
[3]
Kim, I.-J
M.-K. Kim, I.-J. Kim, and J.-S. Lee, CMOS-compatible ferroelectric NAND flash memory for high-density, low-power, and high-speed three-dimensional memory, Sci. Adv.7, eabe1341 (2021)
2021
-
[4]
Schroeder, M
U. Schroeder, M. H. Park, T. Mikolajick, and C. S. Hwang, The fundamentals and applications of ferroelectric HfO 2, Nat. Rev. Mater.7, 653 (2022)
2022
-
[5]
J. P. B. Silva, R. Alcala, U. E. Avci, N. Barrett, L. B´ egon-Lours, M. Borg, S. Byun, S.-C. Chang, S.-W. Cheong, D.-H. Choe, J. Coignus, V. Deshpande, A. Dimoulas, C. Dubourdieu, I. Fina, H. Funakubo, L. Grenouillet, A. Gruverman, J. Heo, M. Hoffmann, H. A. Hsain, F.-T. Huang, C. S. Hwang, J. ´I˜ niguez, J. L. Jones, I. V. Karpov, A. Kersch, T. Kwon, S. ...
2023
-
[6]
Junquera and P
J. Junquera and P. Ghosez, Critical thickness for ferroelectricity in perovskite ultrathin films, Nature422, 506 (2003)
2003
-
[7]
Cao, Q.-L
R. Cao, Q.-L. Yang, H.-X. Deng, S.-H. Wei, J. Robertson, and J.-W. Luo, Softening of the optical phonon by reduced interatomic bonding strength without depolarization, Nature634, 1080 (2024)
2024
-
[8]
Shiraishi, K
T. Shiraishi, K. Katayama, T. Yokouchi, T. Shimizu, T. Oikawa, O. Sakata, H. Uchida, Y. Imai, T. Kiguchi, T. J. Konno, and H. Funakubo, Impact of mechanical stress on ferroelectricity in (Hf0.5Zr0.5)O2 thin films, Appl. Phys. Lett.108, 262904 (2016)
2016
-
[9]
Estand ´ ıa, N
S. Estand ´ ıa, N. Dix, M. F. Chisholm, I. Fina, and F. S´ anchez, Domain-matching epitaxy of ferro- electric Hf0.5Zr0.5O2(111) on La 2/3Sr1/3MnO3(001), Cryst. Growth Des.20, 3801 (2020). 14
2020
-
[10]
Fina and F
I. Fina and F. S´ anchez, Epitaxial ferroelectric HfO 2 films: Growth, properties, and devices, ACS Appl. Electron. Mater.3, 1530 (2021)
2021
-
[11]
H.-J. Lee, M. Lee, K. Lee, J. Jo, H. Yang, Y. Kim, S. C. Chae, U. Waghmare, and J. H. Lee, Scale-free ferroelectricity induced by flat phonon bands in HfO 2, Science369, 1343 (2020)
2020
-
[12]
S. Zhou, J. Zhang, and A. M. Rappe, Strain-induced antipolar phase in hafnia stabilizes robust thin-film ferroelectricity, Sci. Adv.8, eadd5953 (2022)
2022
-
[13]
Y. Wei, P. Nukala, M. Salverda, S. Matzen, H. J. Zhao, J. Momand, A. S. Everhardt, G. Ag- nus, G. R. Blake, P. Lecoeur, B. J. Kooi, J. ´I˜ niguez, B. Dkhil, and B. Noheda, A rhombohedral ferroelectric phase in epitaxially strained Hf 0.5Zr0.5O2 thin films, Nat. Mater.17, 1095 (2018)
2018
-
[14]
Y. Yun, P. Buragohain, M. Li, Z. Ahmadi, Y. Zhang, X. Li, H. Wang, J. Li, P. Lu, L. Tao, H. Wang, J. E. Shield, E. Y. Tsymbal, A. Gruverman, and X. Xu, Intrinsic ferroelectricity in Y-doped HfO 2 thin films, Nat. Mater.21, 903 (2022)
2022
-
[15]
P. W. Tasker, The stability of ionic crystal surfaces, J. Phys. C: Solid State Phys.12, 4977 (1979)
1979
-
[16]
L. Ma, J. Wu, T. Zhu, Y. Huang, Q. Lu, and S. Liu, Ultrahigh oxygen ion mobility in ferroelectric hafnia, Phys. Rev. Lett.131, 256801 (2023)
2023
-
[17]
Haastrup, M
S. Haastrup, M. Strange, M. Pandey, T. Deilmann, P. S. Schmidt, N. F. Hinsche, M. N. Gjerding, D. Torelli, P. M. Larsen, A. C. Riis-Jensen, J. Gath, K. W. Jacobsen, J. Jørgen Mortensen, T. Olsen, and K. S. Thygesen, The computational 2D materials database: High-throughput modeling and discovery of atomically thin crystals, 2D Mater.5, 042002 (2018)
2018
-
[18]
K. Lee, S. W. Kim, Y. Toda, S. Matsuishi, and H. Hosono, Dicalcium nitride as a two-dimensional electride with an anionic electron layer, Nature494, 336 (2013)
2013
-
[19]
J. Liu, S. Liu, J.-Y. Yang, and L. Liu, Electric auxetic effect in piezoelectrics, Phys. Rev. Lett. 125, 197601 (2020)
2020
-
[20]
Dutta, P
S. Dutta, P. Buragohain, S. Glinsek, C. Richter, H. Aramberri, H. Lu, U. Schroeder, E. Defay, A. Gruverman, and J. ´I˜ niguez, Piezoelectricity in hafnia, Nat. Commun.12, 7301 (2021)
2021
-
[21]
B. Lv, T. Zhu, Y. Tang, Y. Lv, C. Zhang, X. Wang, D. Shu, and M. Xiao, Probing permanent dipole moments and removing exciton fine structures in single perovskite nanocrystals by an electric field, Phys. Rev. Lett.126, 197403 (2021)
2021
-
[22]
Schmid, W
M. Schmid, W. Hofer, P. Varga, P. Stoltze, K. W. Jacobsen, and J. K. Nørskov, Surface stress, surface elasticity, and the size effect in surface segregation, Phys. Rev. B51, 10937 (1995). 15
1995
-
[23]
K. F. Garrity, K. M. Rabe, and D. Vanderbilt, Hyperferroelectrics: Proper ferroelectrics with persistent polarization, Phys. Rev. Lett.112, 127601 (2014)
2014
-
[24]
P. Li, X. Ren, G.-C. Guo, and L. He, The origin of hyperferroelectricity in LiBO 3 (B= V, Nb, Ta, Os), Sci. Rep.6, 34085 (2016)
2016
-
[25]
Adhikari and H
R. Adhikari and H. Fu, Hyperferroelectricity in ZnO: Evidence from analytic formulation and numerical calculations, Phys. Rev. B99, 104101 (2019)
2019
-
[26]
Acosta, J
A. Acosta, J. M. P. Martirez, N. Lim, J. P. Chang, and E. A. Carter, Effect of thickness and surface composition on the stability of polarization in ferroelectric Hf xZr1−xO2 thin films, Phys. Rev. Mater.7, 124401 (2023)
2023
-
[27]
S. Lv, Q. Hu, Y. Xue, C. Xue, Z. Wang, G. Teobaldi, T. Cao, and L.-M. Liu, Physical origin of depolarization effects in HfO 2 films, Phys. Rev. Mater.9, 084407 (2025)
2025
-
[28]
Kresse and J
G. Kresse and J. Furthm¨ uller, Efficient iterative schemes forab initiototal-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)
1996
-
[29]
P. E. Bl¨ ochl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)
1994
-
[30]
Kresse and D
G. Kresse and D. Joubert, From ultrasoft pseudopotentials to the projector augmented-wave method, Phys. Rev. B59, 1758 (1999)
1999
-
[31]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett.77, 3865 (1996)
1996
-
[32]
J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, Phys. Rev. Lett.100, 136406 (2008)
2008
-
[33]
H. J. Monkhorst and J. D. Pack, Special points for Brillouin-zone integrations, Phys. Rev. B13, 5188 (1976)
1976
-
[34]
Neugebauer and M
J. Neugebauer and M. Scheffler, Adsorbate-substrate and adsorbate-adsorbate interactions of Na and K adlayers on Al(111), Phys. Rev. B46, 16067 (1992)
1992
-
[35]
R. D. King-Smith and D. Vanderbilt, Theory of polarization of crystalline solids, Phys. Rev. B47, 1651 (1993)
1993
-
[36]
Vanderbilt and R
D. Vanderbilt and R. D. King-Smith, Electric polarization as a bulk quantity and its relation to surface charge, Phys. Rev. B48, 4442 (1993)
1993
-
[37]
X. Wu, D. Vanderbilt, and D. R. Hamann, Systematic treatment of displacements, strains, and electric fields in density-functional perturbation theory, Phys. Rev. B72, 035105 (2005). 16
2005
-
[38]
Momma and F
K. Momma and F. Izumi, VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data, J. Appl. Cryst.44, 1272 (2011)
2011
-
[39]
X. Lyu, F. Ali, T. Song, I. Fina, and F. S´ anchez, Cooperative effects of interface symmetry, redox conditions and low-thickness to improve polarization in ferroelectric Hf 0.5Zr0.5O2 films, ACS Appl. Mater. Interfaces17, 32596 (2025)
2025
-
[40]
J. Yang, J. Wu, J. Li, C. Zhou, Y. Sun, Z. Chen, and S. Liu, Theoretical lower limit of coercive field in ferroelectric hafnia, Phys. Rev. X15, 021042 (2025). 17
2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.