First-principles investigation of small polarons in rhombohedral NaNbO₃
Pith reviewed 2026-05-20 09:27 UTC · model grok-4.3
The pith
Small hole polarons trap stably on oxygen in rhombohedral NaNbO3 while excess electrons do not self-trap.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using density-functional theory corrected by a Hubbard U and the enforced-piecewise-linearity approach with finite-size scaling, the authors determine that a small hole polaron centered on an O-2p orbital possesses a trapping energy of −0.65 eV. Nudged-elastic-band calculations yield an adiabatic migration barrier of 0.32 eV for this polaron. Excess electrons do not form self-trapped states on Nb-4d orbitals, indicating weak electron-phonon coupling. The results identify oxygen as an intrinsic hole trap and show that hole polarons must be included in defect models of NaNbO3-based electroceramics.
What carries the argument
Enforced-piecewise-linearity DFT+U calculations with finite-size scaling for polaron trapping energies, paired with nudged-elastic-band paths for adiabatic migration barriers.
If this is right
- Oxygen atoms function as intrinsic hole traps within NaNbO3.
- Defect models for NaNbO3-based electroceramics must account for hole polarons to predict charge compensation correctly.
- The 0.32 eV migration barrier implies a thermally activated contribution to hole conductivity at operating temperatures.
- Weak electron-phonon coupling means excess electrons are compensated without forming small polarons on niobium sites.
Where Pith is reading between the lines
- These hole-trapping properties may help explain leakage currents and phase-transition behavior observed in NaNbO3 capacitors under bias.
- The same computational protocol could be used to compare polaron stability across related perovskite niobates.
- Targeted doping or oxygen-vacancy control might be tested to tune hole-trapping strength and raise electrical resistivity.
Load-bearing premise
The rhombohedral phase acts as a structurally stable model whose calculated polaron properties represent the intrinsic behavior relevant to capacitor performance.
What would settle it
Direct spectroscopic detection of oxygen-centered hole states or measured hole migration activation energies close to 0.32 eV in rhombohedral NaNbO3 would support the results; absence of localization or substantially different barriers would contradict them.
Figures
read the original abstract
Sodium niobate (NaNbO$_{3}$) is a perovskite oxide and a key component of emerging lead-free antiferroelectric capacitors for high-energy-density applications. However, its performance can be hindered by irreversible phase transitions and leakage currents associated with low electrical resistivity. Defect and doping engineering offers a potential way to overcome these problems, but its use requires a detailed understanding of electronic, ionic, and polaron charge-compensation mechanisms, where the role of polarons remains largely unexplored. Here, we investigate the stability of small hole and electron polarons in rhombohedral NaNbO$_{3}$, which is a structurally well-defined model system that avoids lattice-dynamical instabilities. Trapping energies are calculated using density-functional theory corrected by a Hubbard $U$, using the enforced-piecewise-linearity approach including finite-size scaling. For the small hole-polaron centered on O-2$p$ orbital, we find a trapping energy of $-$0.65 (eV) and an adiabatic migration barrier of 0.32 (eV) determined by nudged-elastic-band calculations. In contrast, we show that excess electrons do not self-trap on Nb-4$d$ orbitals, reflecting weak electron-phonon coupling in the conduction band manifold. These results identify oxygen as an intrinsic hole trap in NaNbO$_{3}$ and highlight the importance of including hole polarons in defect models of NaNbO$_{3}$-based electroceramics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates small polarons in rhombohedral NaNbO3 using DFT+U with enforced piecewise linearity and finite-size scaling. It reports a hole polaron trapping energy of -0.65 eV localized on O-2p orbitals together with a 0.32 eV adiabatic migration barrier obtained from nudged-elastic-band calculations, while demonstrating that excess electrons do not self-trap on Nb-4d states owing to weak electron-phonon coupling. The rhombohedral phase is adopted as a structurally stable model system, and the results are positioned as relevant to defect engineering and leakage-current mitigation in NaNbO3-based capacitors.
Significance. If the central claims hold, the work supplies quantitative, first-principles values for hole-polaron stability and mobility that can be directly incorporated into defect models for antiferroelectric capacitors. The explicit use of enforced piecewise linearity plus finite-size scaling for trapping energies is a methodological strength that reduces self-interaction artifacts and improves transferability. The reported absence of electron self-trapping provides a clear contrast that helps rationalize observed transport asymmetries.
major comments (2)
- [NEB calculations (results section on hole-polaron migration)] The trapping energy of -0.65 eV is obtained with DFT+U under enforced piecewise linearity and finite-size scaling. In contrast, the 0.32 eV adiabatic barrier is reported from separate NEB calculations; the manuscript does not state whether the same linearity-enforcement correction was applied to the NEB images. If standard DFT+U was used for the barrier, the self-interaction error is treated inconsistently, directly undermining the reliability of the hole-polaron transport claim.
- [Electron polaron results subsection] The conclusion that electrons do not self-trap on Nb-4d orbitals is based on the absence of a negative trapping energy. The manuscript should explicitly confirm that the same supercell sizes and finite-size scaling procedure used for holes were applied to the electron case, or provide the raw total-energy differences before scaling to allow independent assessment.
minor comments (2)
- [Methods] The value of the Hubbard U parameter and the precise procedure for enforcing piecewise linearity (e.g., which orbitals and how the potential is adjusted) should be stated in the methods section for reproducibility.
- [Figures] Figure captions for the polaron charge-density plots should include the isosurface value and the supercell size used.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our work and for the detailed comments that help improve the clarity and rigor of the manuscript. We address the major comments point by point below.
read point-by-point responses
-
Referee: [NEB calculations (results section on hole-polaron migration)] The trapping energy of -0.65 eV is obtained with DFT+U under enforced piecewise linearity and finite-size scaling. In contrast, the 0.32 eV adiabatic barrier is reported from separate NEB calculations; the manuscript does not state whether the same linearity-enforcement correction was applied to the NEB images. If standard DFT+U was used for the barrier, the self-interaction error is treated inconsistently, directly undermining the reliability of the hole-polaron transport claim.
Authors: We appreciate this important observation regarding methodological consistency. The NEB calculations for the migration barrier were indeed performed using the same DFT+U approach with enforced piecewise linearity applied to each image along the path, and finite-size scaling was used for the total energies. However, we acknowledge that this was not explicitly stated in the manuscript. In the revised version, we will add a clear statement in the methods and results sections confirming that the same corrections were applied to ensure consistency in the treatment of self-interaction errors. This will strengthen the reliability of the reported 0.32 eV barrier. revision: yes
-
Referee: [Electron polaron results subsection] The conclusion that electrons do not self-trap on Nb-4d orbitals is based on the absence of a negative trapping energy. The manuscript should explicitly confirm that the same supercell sizes and finite-size scaling procedure used for holes were applied to the electron case, or provide the raw total-energy differences before scaling to allow independent assessment.
Authors: We thank the referee for this suggestion to improve transparency. The same supercell sizes and finite-size scaling procedure were applied to the electron polaron calculations as for the holes. The trapping energy for electrons remained positive even after scaling, indicating no self-trapping. To allow for independent assessment, we will include the raw total-energy differences before and after scaling in the supplementary information of the revised manuscript. revision: yes
Circularity Check
No circularity: trapping energies and barriers are direct computational outputs
full rationale
The paper computes trapping energies via total-energy differences in supercell DFT+U calculations that incorporate the enforced-piecewise-linearity correction and finite-size scaling as methodological inputs. The adiabatic migration barrier is obtained separately from nudged-elastic-band calculations. Neither quantity is defined in terms of the other, nor does any reported value reduce by construction to a fitted parameter or self-citation chain. The central claims rest on explicit numerical results rather than tautological redefinitions, satisfying the criteria for a self-contained first-principles derivation.
Axiom & Free-Parameter Ledger
free parameters (1)
- Hubbard U
axioms (2)
- domain assumption The enforced-piecewise-linearity approach combined with finite-size scaling yields accurate polaron trapping energies.
- domain assumption Rhombohedral NaNbO3 avoids lattice-dynamical instabilities and serves as a well-defined model for polaron behavior.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Trapping energies are calculated using density-functional theory corrected by a Hubbard U, using the enforced-piecewise-linearity approach including finite-size scaling. For the small hole-polaron ... trapping energy of −0.65 (eV) and an adiabatic migration barrier of 0.32 (eV) determined by nudged-elastic-band calculations.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
We therefore enforceϵ p(Q) =ϵ p(0), withQ= +1 (−1) for a hole (electron) polaron as shown in Fig
Determination of theUparameter via enforced piecewise-linearity For an exact functional, the combination of piecewise linearity and Janak’s theorem [82] implies that the energy of the localized (polaron) levelϵ p(q) is constant between integer occupations. We therefore enforceϵ p(Q) =ϵ p(0), withQ= +1 (−1) for a hole (electron) polaron as shown in Fig. 2....
-
[2]
Trapping energy and stability of small hole polaron Fig. 3(a) and Fig. 3(b) display the density of states (DOS) for the localized and delocalized configurations, respectively. The delocalized state corresponds to an ex- tra hole being uniformly distributed across all ions, re- sulting in a metallic state where the Fermi energy crosses the valence band. In...
-
[3]
T. R. Shrout and S. Zhang, Journal of Electroceramics 19, 185 (2007)
work page 2007
-
[4]
J. R¨ odel, W. Jo, K. T. P. Seifert, E.-M. Anton, T. Granzow, and D. Damjanovic, Journal of the Amer- ican Ceramic Society92, 1153 (2009)
work page 2009
-
[5]
A. J. Bell and O. Deubzer, MRS Bulletin43, 581–587 (2018)
work page 2018
-
[6]
B. Xu, J. ´I˜ niguez, and L. Bellaiche, Nature Communi- cations8, 15682 (2017)
work page 2017
-
[7]
L. Yang, X. Kong, F. Li, H. Hao, Z. Cheng, H. Liu, J.-F. Li, and S. Zhang, Progress in Materials Science102, 72 (2019)
work page 2019
-
[8]
D. Yang, J. Gao, L. Shu, Y.-X. Liu, J. Yu, Y. Zhang, X. Wang, B.-P. Zhang, and J.-F. Li, Journal of Materials Chemistry A8, 23724 (2020)
work page 2020
- [9]
-
[10]
J. Wu, D. Xiao, and J. Zhu, Chem. Rev.115, 2559 (2015)
work page 2015
- [11]
- [12]
-
[14]
Z. Liu, T. Lu, J. Ye, G. Wang, X. Dong, R. Withers, and Y. Liu, Advanced Materials Technologies3, 1800111 (2018)
work page 2018
-
[15]
J.-F. Li, K. Wang, F.-Y. Zhu, L.-Q. Cheng, and F.-Z. Yao, Journal of the American Ceramic Society96, 3677 (2013)
work page 2013
- [16]
-
[17]
H. D. Megaw, Ferroelectrics7, 87 (1974)
work page 1974
-
[18]
S. Lanfredi, M. H. Lente, and J. A. Eiras, Applied Physics Letters80, 2731 (2002)
work page 2002
-
[19]
L. A. Reznitchenko, A. V. Turik, E. M. Kuznetsova, and V. P. Sakhnenko, Journal of Physics: Condensed Matter 13, 3875 (2001)
work page 2001
- [20]
- [21]
- [22]
-
[23]
J. Koruza, J. Tellier, B. Maliˇ c, V. Bobnar, and M. Kosec, Journal of Applied Physics108(2010), 10.1063/1.3512980
-
[24]
N. Luo, L. Ma, G. Luo, C. Xu, L. Rao, Z. Chen, Z. Cen, Q. Feng, X. Chen, F. Toyohisa, Y. Zhu, J. Hong, J.-F. Li, and S. Zhang, Nature Communications14, 1776 (2023)
work page 2023
-
[25]
K. E. Johnston, C. C. Tang, J. E. Parker, K. S. Knight, P. Lightfoot, and S. E. Ashbrook, Journal of the Amer- ican Chemical Society132, 8732 (2010)
work page 2010
-
[26]
C. S. Htet, S. Nayak, A. Manj´ on-Sanz, J. Liu, J. Kong, D. R. Sørensen, F. Marlton, M. R. V. Jørgensen, and A. Pramanick, Phys. Rev. B105, 174113 (2022)
work page 2022
-
[27]
H. Guo, H. Shimizu, Y. Mizuno, and C. A. Randall, Journal of Applied Physics118, 054102 (2015)
work page 2015
-
[28]
S. Egert, M.-H. Zhang, J. Koruza, P. B. Groszewicz, and G. Buntkowsky, The Journal of Physical Chemistry C 124, 23852 (2020)
work page 2020
- [29]
-
[30]
H. Shimizu, H. Guo, S. E. Reyes-Lillo, Y. Mizuno, K. M. Rabe, and C. A. Randall, Dalton Transactions44, 10763 (2015)
work page 2015
-
[31]
Z. Wang, D. Li, W. Liu, L. He, D. Xu, J. Liu, J. Ren, X. Wang, Y. Liu, G. He, J. Bao, Z. Fang, G. Yan, X. Liang, T. Zhou, W. Zhao, W. Liu, D. Wang, and D. Zhou, Nature Communications16, 2892 (2025)
work page 2025
- [32]
-
[33]
Y. Shiratori, A. Magrez, J. Dornseiffer, F.-H. Haegel, C. Pithan, and R. Waser, The Journal of Physical Chem- istry B109, 20122 (2005)
work page 2005
-
[34]
X. Tan, Z. Xu, X. Liu, and Z. Fan, Materials Research Letters6, 159 (2018)
work page 2018
-
[35]
S. Aso, H. Matsuo, and Y. Noguchi, Scientific Reports 13, 6771 (2023)
work page 2023
-
[36]
J. Ye, G. Wang, X. Chen, F. Cao, and X. Dong, Applied Physics Letters114, 122901 (2019)
work page 2019
-
[37]
F. Ye, H. Cheng, X. Jiang, L. Zhu, W. Shi, and J. Xue, Journal of Advanced Ceramics (2025), 10.26599/JAC.2025.9221119
-
[38]
H. Qi, A. Xie, J. Fu, and R. Zuo, Acta Materialia208, 116710 (2021)
work page 2021
-
[39]
B. Luo, H. Dong, D. Wang, and K. Jin, Journal of the American Ceramic Society101, 3460 (2018)
work page 2018
- [40]
-
[41]
W. Sakamoto, Y. ichi Hamazaki, H. Maiwa, M. Moriya, and T. Yogo, Thin Solid Films518, 4256 (2010)
work page 2010
- [42]
-
[43]
J. T. Devreese, Physics Today68, 11 (2015)
work page 2015
-
[44]
C. Franchini, M. Reticcioli, M. Setvin, and U. Diebold, Nature Reviews Materials6, 560 (2021)
work page 2021
-
[45]
M. I. Dykman and E. I. Rashba, Physics Today68, 10 (2015)
work page 2015
-
[46]
A. S. Alexandrov and J. T. Devreese,Advances in polaron physics, Vol. 159 (Springer, 2010)
work page 2010
-
[47]
Emin,Polarons(Cambridge University Press, 2012)
D. Emin,Polarons(Cambridge University Press, 2012)
work page 2012
- [48]
-
[49]
H. Fr¨ ohlich, H. Pelzer, and S. Z. and, The London, Edin- burgh, and Dublin Philosophical Magazine and Journal of Science41, 221 (1950)
work page 1950
- [50]
- [51]
- [52]
-
[53]
G.-Z. Wang, H. Chen, G. Wu, A.-L. Kuang, and H.-K. Yuan, ChemPhysChem17, 489 (2016)
work page 2016
-
[54]
N. Bein, B. Kmet, T. Rojac, A. B. c. v. Golob, B. Maliˇ c, J. Moxter, T. Schneider, L. Fulanovic, M. Azadeh, T. Fr¨ omling, S. Egert, H. Wang, P. van Aken, J. Schwarzkopf, and A. Klein, Phys. Rev. Mater. 6, 084404 (2022)
work page 2022
-
[55]
A. Klein, K. Albe, N. Bein, O. Clemens, K. A. Creutz, P. Erhart, M. Frericks, E. Ghorbani, J. P. Hofmann, B. Huang, B. Kaiser, U. Kolb, J. Koruza, C. K¨ ubel, K. N. S. Lohaus, J. R¨ odel, J. Rohrer, W. Rheinheimer, R. A. D. Souza, V. Streibel, A. Weidenkaff, M. Widen- meyer, B.-X. Xu, and H. Zhang, Journal of Electroce- ramics51, 147 (2023)
work page 2023
-
[56]
Y. Hermans, A. Klein, H. P. Sarker, M. N. Huda, H. Junge, T. Toupance, and W. Jaegermann, Advanced Functional Materials30, 1910432 (2020)
work page 2020
- [57]
-
[58]
T. Tadano and S. Tsuneyuki, Journal of the Ceramic Society of Japan127, 404 (2019)
work page 2019
-
[59]
I. Pallikara, P. Kayastha, J. M. Skelton, and L. D. Whal- ley, Electronic Structure4, 033002 (2022)
work page 2022
-
[60]
M. Kotiuga, S. Halilov, B. Kozinsky, M. Fornari, N. Marzari, and G. Pizzi, Phys. Rev. Res.4, L012042 (2022)
work page 2022
-
[61]
O. F. Schirmer, Journal of Physics: Condensed Matter 21, 123201 (2009)
work page 2009
-
[62]
F. Schmidt, A. L. Kozub, T. Biktagirov, C. Eigner, C. Sil- berhorn, and U. Gerstmann, Physical Review Research 2, 043002 (2020)
work page 2020
-
[63]
S. Sanna, in2015 Symposium on Piezoelectricity, Acous- tic Waves, and Device Applications (SPAWDA)(2015) pp. 523–527
work page 2015
- [64]
- [65]
-
[66]
M. Friedrich, W. G. Schmidt, A. Schindlmayr, and S. Sanna, Phys. Rev. Mater.1, 054406 (2017)
work page 2017
-
[67]
G. Kresse and J. Furthm¨ uller, Computational Materials Science6, 15 (1996)
work page 1996
-
[68]
P. E. Bl¨ ochl, Phys. Rev. B50, 17953 (1994)
work page 1994
- [69]
-
[70]
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)
work page 1996
-
[71]
J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Phys. Rev. Lett.100, 136406 (2008)
work page 2008
-
[72]
T. D. Pham and N. A. Deskins, Journal of Chemical The- ory and Computation16, 5264 (2020)
work page 2020
-
[73]
V. I. Anisimov, J. Zaanen, and O. K. Andersen, Phys. Rev. B44, 943 (1991)
work page 1991
-
[74]
S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, Phys. Rev. B57, 1505 (1998)
work page 1998
-
[75]
J. P. Allen and G. W. Watson, Physical Chemistry Chem- ical Physics16, 21016 (2014)
work page 2014
-
[76]
S. Falletta, J. Wiktor, and A. Pasquarello, Phys. Rev. B 102, 041115 (2020)
work page 2020
-
[77]
S. Falletta and A. Pasquarello, npj Computational Ma- terials8, 263 (2022)
work page 2022
-
[78]
Nudged elastic band method for finding minimum energy paths of transitions,
H. J´ onsson, G. Mills, and K. W. Jacobsen, “Nudged elastic band method for finding minimum energy paths of transitions,” inClassical and Quantum Dynamics in Condensed Phase Simulations, pp. 385–404
- [79]
-
[80]
S. K. Mishra, N. Choudhury, S. L. Chaplot, P. S. R. Krishna, and R. Mittal, Phys. Rev. B76, 024110 (2007)
work page 2007
-
[81]
H. H. Nahm and C. H. Park, Phys. Rev. B78, 184108 (2008)
work page 2008
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.