pith. sign in

arxiv: 2605.18258 · v1 · pith:PHNR4JH5new · submitted 2026-05-18 · ❄️ cond-mat.mtrl-sci

Why hole polaron formation on oxygen is limiting the Fermi level in Fe acceptor doped BaTiO₃ under oxidizing conditions

Mohammad Amirabbasi , Emre Erdem , Denis Sudarikov , Jochen Rohrer , Andreas Klein , Karsten Albe This is my paper

Pith reviewed 2026-05-20 09:47 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords hole polaronFe-doped BaTiO3oxygen polaronFermi levelacceptor dopingligand holeDFT+Uperovskite
0
0 comments X

The pith

Oxidizing Fe-doped BaTiO3 favors oxygen-centered hole polarons over Fe4+ centers.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper establishes that in Fe acceptor-doped BaTiO3 under oxidizing conditions the oxidizing hole prefers to localize on oxygen rather than converting the Fe3+ acceptor into Fe4+. Density-functional calculations with a Hubbard correction on oxygen states show the resulting Fe3+-O- polaron complex lies lower in energy than the formal Fe4+ state. This resolves the mismatch between the expected oxidation product and the Fe3+-dominated EPR spectra that experiments actually record. A reader cares because the result shows how oxygen polarons can pin the Fermi level and shape the defect response in acceptor-doped ferroelectric perovskites.

Core claim

Using density-functional theory with occupation-matrix control and a piecewise-linearity-based Hubbard correction for O-2p states, we find that an oxygen-centered hole polaron forming a Fe3+-O- complex is lower in energy than the formal Fe4+ configuration. Our results identify ligand-hole formation as a favorable charge-compensation mechanism in oxidized Fe-doped BaTiO3 and provide an explanation for the predominance of Fe3+-based centers in spectroscopy. More broadly, they show how oxygen polarons can limit Fermi-level shifts and control the electronic response of acceptor-doped ferroelectric perovskites.

What carries the argument

The oxygen-centered hole polaron that forms the Fe3+-O- complex and carries the hole instead of oxidizing iron.

If this is right

  • Fermi-level movement is restricted once the oxygen polaron becomes the stable compensation defect.
  • Spectroscopic methods register mainly Fe3+ signatures because Fe4+ is not the ground-state configuration.
  • Charge compensation in oxidized acceptor-doped perovskites proceeds through ligand holes on oxygen.
  • The electronic and defect properties of the material are governed by the stability of these oxygen polarons.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same ligand-hole preference may appear in other acceptor-doped oxide perovskites and limit their Fermi-level range.
  • Engineering oxygen stoichiometry or vacancy concentration could offer a route to adjust the pinned Fermi position.
  • Similar calculations on related compounds would test whether oxygen polarons are a general feature of oxidized acceptor-doped ferroelectrics.

Load-bearing premise

The piecewise-linearity-based Hubbard correction applied to O-2p states with occupation-matrix control correctly identifies the polaron ground state without artificial stabilization or destabilization of the localized hole.

What would settle it

Direct detection of dominant Fe4+ centers by EPR or other spectroscopy in strongly oxidized Fe-doped BaTiO3, or a calculation showing Fe4+ lower in energy than the Fe3+-O- complex, would falsify the claimed preference.

Figures

Figures reproduced from arXiv: 2605.18258 by Andreas Klein, Denis Sudarikov, Emre Erdem, Jochen Rohrer, Karsten Albe, Mohammad Amirabbasi.

Figure 1
Figure 1. Figure 1: FIG. 1. (Color online) Two charge-compensation configu [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (Color online) Density of states (DOS) for (a) a lo [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (Color online) (a) Total energies for the Fe [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
read the original abstract

Oxidizing Fe-doped BaTiO$_3$ is commonly expected to convert substitutional Fe$^{3+}$ acceptors into formal Fe$^{4+}$ centers. Yet, the experimentally accessible picture based on electron-paramagnetic resonance (EPR) is dominated by Fe$^{3+}$-related signatures, while Fe$^{4+}$ is not a straightforward observable. Here we show that this apparent discrepancy reflects the preferred location of the oxidizing hole: not on Fe, but on oxygen. Using density-functional theory with with occupation-matrix control and a piecewise-linearity-based Hubbard correction (DFT+$U$) for O-2$p$ states, we find that an oxygen-centered hole polaron is forming a Fe$^{3+}$-O$^{-}$ complex that is lower in energy than the formal Fe$^{4+}$ configuration. Our results identify ligand-hole formation as a favorable charge-compensation mechanism in oxidized Fe-doped BaTiO$_3$ and provide an explanation for the predominance of Fe$^{3+}$-based centers in spectroscopy. More broadly, they show how oxygen polarons can limit Fermi-level shifts and control the electronic response of acceptor-doped ferroelectric perovskites.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript examines the preferred location of the oxidizing hole in Fe acceptor doped BaTiO3 under oxidizing conditions. Using DFT+U with occupation-matrix control and a piecewise-linearity-based Hubbard U for O-2p states, it concludes that an oxygen-centered hole polaron forming a Fe3+-O- complex is lower in energy than the formal Fe4+ configuration. This explains the predominance of Fe3+ signatures in EPR spectroscopy and how oxygen polarons limit Fermi level shifts in such materials.

Significance. If the central result holds, the work provides a microscopic rationale for experimental spectroscopic observations in acceptor-doped BaTiO3 and more generally in ferroelectric perovskites. It highlights ligand-hole formation as a key charge-compensation mechanism. The approach of deriving U from piecewise linearity is noted as a positive aspect for reducing empiricism.

major comments (2)
  1. [Abstract] The central claim that the oxygen-centered hole polaron is lower in energy than the formal Fe4+ is presented without any reported numerical energy differences, convergence tests, or quantitative comparison to experiment, which is essential to evaluate the robustness of the energy ordering.
  2. [Methods/Results on DFT+U] The piecewise-linearity-based Hubbard correction applied specifically to O-2p states combined with occupation-matrix control may preferentially stabilize the localized polaron configuration. The manuscript should demonstrate that the reported energy preference is not an artifact by showing results for different U values or without the O-specific correction, as the Fe4+ configuration may not receive equivalent treatment.
minor comments (1)
  1. [Abstract] There is a typographical error: 'density-functional theory with with occupation-matrix control' should be corrected to remove the duplicate 'with'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and the constructive major comments. We address each point below and will revise the manuscript to strengthen the presentation of quantitative results and methodological robustness.

read point-by-point responses

  1. Referee: [Abstract] The central claim that the oxygen-centered hole polaron is lower in energy than the formal Fe4+ is presented without any reported numerical energy differences, convergence tests, or quantitative comparison to experiment, which is essential to evaluate the robustness of the energy ordering.

    Authors: We agree that explicit numerical values are needed to allow readers to assess the energy ordering. In the revised manuscript we will report the calculated total-energy difference between the oxygen-centered hole polaron (Fe^{3+}-O^{-}) configuration and the formal Fe^{4+} state, together with the supercell-size and k-point convergence tests that were performed. Direct quantitative comparison with experiment is not straightforward because EPR reports the dominant paramagnetic center rather than formation energies; we will expand the discussion to clarify this point and to relate the computed preference to the observed predominance of Fe^{3+} signatures. revision: yes

  2. Referee: [Methods/Results on DFT+U] The piecewise-linearity-based Hubbard correction applied specifically to O-2p states combined with occupation-matrix control may preferentially stabilize the localized polaron configuration. The manuscript should demonstrate that the reported energy preference is not an artifact by showing results for different U values or without the O-specific correction, as the Fe4+ configuration may not receive equivalent treatment.

    Authors: The piecewise-linearity condition is applied separately to the O-2p manifold to remove self-interaction error for the hole, while Fe-3d states receive their own U derived from the same procedure; both configurations are therefore treated within the same DFT+U framework. Occupation-matrix control is used only to enforce the desired orbital occupancy and does not alter the relative energetics once the self-consistent solution is reached. To address the referee’s concern directly, the revised manuscript will include additional calculations performed with U_O = 0 eV and with a range of U_O values around the piecewise-linear value, confirming that the oxygen-polaron configuration remains lower in energy. We will also discuss results obtained without occupation-matrix control to demonstrate that the ordering is robust. revision: yes

Circularity Check

0 steps flagged

DFT+U with piecewise-linearity U on O-2p yields independent energy ordering

full rationale

The central result is an energy comparison between the oxygen-centered hole polaron (Fe3+-O-) and the formal Fe4+ configuration, obtained from DFT+U calculations that apply a Hubbard correction derived from the piecewise-linearity condition on O-2p states together with occupation-matrix control. This procedure is a standard first-principles route to determine U and localize the hole; the resulting ordering is not forced by construction because the U value is fixed by the linearity requirement rather than by the target energy difference or by any spectroscopic observable. No self-definitional step, fitted-input prediction, or load-bearing self-citation chain appears in the derivation. The method remains externally falsifiable against EPR data and is therefore scored as minor (non-circular) methodological choice.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the accuracy of the chosen Hubbard correction for oxygen p-states and the assumption that the supercell model captures the relevant polaron localization without significant finite-size artifacts.

free parameters (1)
  • Hubbard U value for O-2p states
    Piecewise-linearity-based U is introduced to stabilize the oxygen hole polaron; its specific numerical value is a free parameter chosen to enforce linearity.
axioms (1)
  • domain assumption Occupation-matrix control in DFT+U can reliably locate the ground-state polaron configuration on oxygen ligands.
    The method is invoked to compare energies of Fe3+-O- versus Fe4+ states.

pith-pipeline@v0.9.0 · 5765 in / 1267 out tokens · 43310 ms · 2026-05-20T09:47:39.686923+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

42 extracted references · 42 canonical work pages

  1. [1]

    Iqbal, A

    N. Iqbal, A. Dixit, P. S. Dobal, R. S. Katiyar, and A. S. Bhalla, Journal of Materials Science: Materials in Elec- tronics36, 1008 (2025)

    work page 2025

  2. [2]

    Acosta, N

    M. Acosta, N. Novak, V. Rojas, S. Patel, R. Vaish, J. Ko- ruza, G. A. Rossetti, and J. R¨ odel, Applied Physics Re- views4, 041305 (2017)

    work page 2017

  3. [3]

    Y. B. Adediji, A. M. Adeyinka, D. I. Yahya, and O. V. Mbelu, Energy, Ecology and Environment8, 401 (2023)

    work page 2023

  4. [4]

    G. Y. Yang, G. D. Lian, E. C. Dickey, C. A. Randall, D. E. Barber, P. Pinceloup, M. A. Henderson, R. A. Hill, J. J. Beeson, and D. J. Skamser, Journal of Applied Physics96, 7500 (2004)

    work page 2004

  5. [5]

    Suzuki, L

    I. Suzuki, L. Gura, and A. Klein, Physical Chemistry Chemical Physics21, 6238 (2019)

    work page 2019

  6. [6]

    R. N. Schwartz and B. A. Wechsler, Physical Review B 48, 7057 (1993)

    work page 1993

  7. [7]

    Possenriede, P

    E. Possenriede, P. Jacobs, and O. Schirmer, Journal of Physics: Condensed Matter4, 4719 (1992)

    work page 1992

  8. [8]

    V. V. Laguta, A. M. Slipenyuk, I. P. Bykov, M. D. Glinchuk, M. Maglione, D. Michau, J. Rosa, and L. Jas- trabik, Applied Physics Letters87, 022903 (2005)

    work page 2005

  9. [9]

    R. A. Theerthan, A. Artemenko, and M. Maglione, Jour- 5 nal of Physics: Condensed Matter24, 405901 (2012)

    work page 2012

  10. [10]

    O. F. Schirmer, Journal of Physics: Condensed Matter 18, R667 (2006)

    work page 2006

  11. [11]

    O. F. Schirmer, Journal of Physics: Condensed Matter 23, 334218 (2011)

    work page 2011

  12. [12]

    Franchini, M

    C. Franchini, M. Reticcioli, M. Setvin, and U. Diebold, Nature Reviews Materials6, 560 (2021)

    work page 2021

  13. [13]

    Varnhorst, O

    T. Varnhorst, O. F. Schirmer, H. Kr¨ ose, R. Scharfschw- erdt, and T. W. Kool, Phys. Rev. B53, 116 (1996)

    work page 1996

  14. [14]

    Klein, W

    A. Klein, W. Jaegermann, K. Aicher, T. Berger,et al., Journal of Electroceramics47, 405 (2023)

    work page 2023

  15. [15]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Computational Materials Science6, 15 (1996)

    work page 1996

  16. [16]

    Kresse and D

    G. Kresse and D. Joubert, Phys. Rev. B59, 1758 (1999)

    work page 1999

  17. [17]

    P. E. Bl¨ ochl, Phys. Rev. B50, 17953 (1994)

    work page 1994

  18. [18]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

    work page 1996

  19. [19]

    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

  20. [20]

    BaTiO 3, rhombohedral (BaTiO3 lt2, T = 40 K) crystal structure,

    P. Villars and K. Cenzual, “BaTiO 3, rhombohedral (BaTiO3 lt2, T = 40 K) crystal structure,” Springer- Materials, PAULING FILE Multinaries Edition – 2022 (2023), springer-Verlag Berlin Heidelberg; Material Phases Data System (MPDS), Switzerland; National In- stitute for Materials Science (NIMS), Japan

    work page 2022

  21. [21]

    Gajdoˇ s, K

    M. Gajdoˇ s, K. Hummer, G. Kresse, J. Furthm¨ uller, and F. Bechstedt, Phys. Rev. B73, 045112 (2006)

    work page 2006

  22. [23]

    Baroni, S

    S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Gian- nozzi, Rev. Mod. Phys.73, 515 (2001)

    work page 2001

  23. [24]

    Falletta, J

    S. Falletta, J. Wiktor, and A. Pasquarello, Phys. Rev. B 102, 041115 (2020)

    work page 2020

  24. [25]

    Erhart, A

    P. Erhart, A. Klein, D. ˚Aberg, and B. Sadigh, Phys. Rev. B90, 035204 (2014)

    work page 2014

  25. [26]

    ¨Osterbacka, P

    N. ¨Osterbacka, P. Erhart, S. Falletta, A. Pasquarello, and J. Wiktor, Chemistry of Materials32, 8393 (2020)

    work page 2020

  26. [27]

    J. P. Allen and G. W. Watson, Physical Chemistry Chem- ical Physics16, 21016 (2014)

    work page 2014

  27. [28]

    Wang and C.-C

    B. Wang and C.-C. Sun, Appl. Opt.40, 672 (2001)

    work page 2001

  28. [29]

    X. Wu, D. Vanderbilt, and D. R. Hamann, Phys. Rev. B72, 035105 (2005)

    work page 2005

  29. [30]

    G. C. Moore, M. K. Horton, E. Linscott, A. M. Ganose, M. Siron, D. D. O’Regan, and K. A. Persson, Phys. Rev. Mater.8, 014409 (2024)

    work page 2024

  30. [31]

    A. S. Alexandrov and J. T. Devreese,Advances in polaron physics, Vol. 159 (Springer, 2010)

    work page 2010

  31. [32]

    Emin,Polarons(Cambridge University Press, 2013)

    D. Emin,Polarons(Cambridge University Press, 2013)

    work page 2013

  32. [33]

    Erhart, A

    P. Erhart, A. Klein, D. ˚Aberg, and B. Sadigh, Physical Review B90, 035204 (2014)

    work page 2014

  33. [34]

    Falletta and A

    S. Falletta and A. Pasquarello, npj Computational Ma- terials8, 263 (2022)

    work page 2022

  34. [35]

    J. F. Janak, Phys. Rev. B18, 7165 (1978)

    work page 1978

  35. [36]

    Zorko, M

    A. Zorko, M. Pregelj, M. Gomilˇ sek, Z. Jagliˇ ci´ c, D. Paji´ c, M. Telling, I. Arˇ con, I. Mikulska, and M. Valant, Scien- tific Reports5, 7703 (2015)

    work page 2015

  36. [37]

    Ivanov, R

    S. Ivanov, R. Tellgren, F. Porcher, T. Ericsson, A. Mo- sunov, P. Beran, S. Korchagina, P. A. Kumar, R. Math- ieu, and P. Nordblad, Materials Research Bulletin47, 3253 (2012)

    work page 2012

  37. [38]

    I. W. Boateng, R. Tia, E. Adei, N. Y. Dzade, C. R. A. Catlow, and N. H. De Leeuw, Physical Chemistry Chem- ical Physics19, 7399 (2017)

    work page 2017

  38. [39]

    B¨ ottcher, H

    R. B¨ ottcher, H. Langhammer, T. M¨ uller, and H. Abicht, Journal of Physics: Condensed Matter17, 4925 (2005)

    work page 2005

  39. [40]

    K. A. M¨ uller, Physical review letters2, 341 (1959)

    work page 1959

  40. [41]

    Scharfschwerdt, O

    R. Scharfschwerdt, O. Schirmer, H. Hesse, and D. Rytz, Applied Physics B68, 807 (1999)

    work page 1999

  41. [42]

    Mukherjee, D

    S. Mukherjee, D. Phuyal, C. U. Segre, S. Das, O. Karis, T. Edvinsson, and H. Rensmo, The Journal of Physical Chemistry C125, 14910 (2021)

    work page 2021

  42. [43]

    Lu, Q.-L

    D.-Y. Lu, Q.-L. Liu, T. Ogata, X.-Y. Sun, and X.-F. Wang, Japanese Journal of Applied Physics50, 035806 (2011)

    work page 2011