Inter-defect interactions, oxygen-vacancy distribution, and oxidation in acceptor-doped ABO3 perovskites
Pith reviewed 2026-05-16 20:12 UTC · model grok-4.3
The pith
In acceptor-doped ABO3 perovskites, oxygen vacancy-impurity interactions outweigh inter-vacancy correlations in shaping defect distribution, local coordination, and oxidation.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that, for realistic interaction energies, the interaction between oxygen vacancies and acceptor impurities exerts a stronger influence on defect thermodynamics, local ion coordination, and oxidation than correlations among the vacancies themselves. Inter-site vacancy repulsion grows noticeable only at high dopant content x, while on-site correlations appear in a narrow moderate-x window. Non-uniform impurity placement, which can arise during sample preparation, markedly changes the oxygen-vacancy distribution while affecting short-range order and oxidation more weakly. Overall, the inter-defect interactions reduce hole concentration, increase oxidation enthalpy, and can,
What carries the argument
Statistical theory of interacting defects combined with Monte Carlo lattice simulations that compute oxygen-vacancy distribution and oxidation thermodynamics in the presence of acceptor impurities.
If this is right
- At sufficiently high dopant fractions, inter-site repulsion between vacancies becomes detectable in the vacancy distribution.
- Non-uniform impurity placement produced by typical synthesis routes strongly perturbs oxygen-vacancy distribution while leaving short-range order and oxidation enthalpy less affected.
- Inter-defect interactions lower the equilibrium hole concentration and raise the oxidation enthalpy relative to the non-interacting case.
- The dependence of oxidation enthalpy and hole concentration on dopant level x can become non-monotonic once interactions are included.
- The calculated local ion coordination and hole conductivity trends reproduce experimental observations in acceptor-doped perovskites.
Where Pith is reading between the lines
- Synthesis protocols that deliberately homogenize or cluster the acceptor ions could be used to tune vacancy distribution and oxidation resistance without changing overall composition.
- The same interaction framework may be relevant to ionic conductivity in related perovskites employed as electrolytes or electrodes.
- Temperature-dependent simulations extending the present Monte Carlo runs would allow direct comparison with conductivity data taken at device operating temperatures.
- Assuming a perfectly random impurity distribution in device modeling could introduce systematic errors in predicted defect concentrations.
Load-bearing premise
The chosen numerical values for the defect interaction energies are realistic for the physical materials and that the lattice model captures all relevant effects without needing additional terms such as lattice relaxation or long-range strain.
What would settle it
A measured oxidation enthalpy versus dopant concentration x that deviates systematically from the non-monotonic or elevated values predicted by the model in a well-characterized ABO3 compound such as Sr-doped LaMnO3.
read the original abstract
The effects of inter-defect interaction and impurity disorder on defect thermodynamics, local ion coordination, and oxidation in acceptor-doped wide-gap ABO3 perovskites are explored using the developed statistical theory and Monte Carlo simulations. The results demonstrate that under realistic energy parameters the interaction between oxygen vacancies and impurities generally has a greater impact on the studied properties than inter-vacancy correlations. The influence of inter-vacancy interaction significantly depends on dopant content x: inter-site vacancy repulsion becomes noticeable at sufficiently high x, whereas on-site correlations can be pronounced within a narrow doping range at moderate x values. It is found that a non-uniform impurity allocation, which can result from a sample preparation procedure, considerably affects oxygen-vacancy distribution, and has a weaker effect on short-range order and oxidation. It is also shown that inter-defect interaction reduces the hole concentration, increases the oxidation enthalpy, and can result in their non-trivial dependence on x. The obtained results agree with experimental data on local ion coordination and help to explain the behavior of hole conductivity in the considered perovskites. The findings of this study contribute to understanding the fundamental properties of acceptor-doped oxides, facilitating the development of new materials for clean energy applications.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a statistical theory and Monte Carlo simulations to examine how inter-defect interactions and impurity disorder affect oxygen-vacancy thermodynamics, local coordination, and oxidation behavior in acceptor-doped ABO3 perovskites. It concludes that, under realistic energy parameters, vacancy-impurity interactions generally dominate inter-vacancy correlations, with doping-dependent effects on vacancy distribution, oxidation enthalpy, and hole concentration; non-uniform impurity allocation also influences results, and the findings align with experimental data on local ion coordination and hole conductivity.
Significance. If the interaction energies are independently validated, the work offers useful quantitative insight into defect correlations in perovskites and their impact on functional properties relevant to energy materials. The combination of analytical statistical theory with Monte Carlo sampling to capture both short- and long-range effects is a methodological strength, and the reported agreement with coordination experiments provides a concrete anchor. The doping-dependent crossover between interaction regimes is potentially falsifiable and could guide material optimization.
major comments (2)
- [Methods / Energy-parameter section] The central claim that vacancy-impurity interactions dominate inter-vacancy correlations rests on the specific numerical values chosen for the defect interaction energies. The manuscript must explicitly state whether these parameters were obtained from independent first-principles calculations, fitted to a subset of experiments with reported uncertainty, or selected by hand; without this, the dominance result risks being parameter-dependent rather than physically predictive. A sensitivity scan over plausible ranges of the vacancy-impurity binding versus vacancy-vacancy repulsion should be added.
- [Results / Oxidation and conductivity subsection] The reported increase in oxidation enthalpy and non-monotonic hole-concentration behavior with doping x are direct outputs of the chosen energies. These predictions should be compared quantitatively to measured oxidation enthalpies or conductivity data across a wider doping range, with explicit error bars, to test whether the model reproduces the experimental trends without post-hoc adjustment.
minor comments (3)
- Clarify the precise definition of 'on-site' versus 'inter-site' correlations in the statistical theory and how they are implemented in the Monte Carlo moves.
- Add a table or figure summarizing the numerical values of all interaction energies used, together with their literature sources or derivation method.
- Ensure that all simulation cell sizes, equilibration times, and averaging procedures are stated so that the Monte Carlo results can be reproduced.
Simulated Author's Rebuttal
We thank the referee for the constructive comments and positive assessment of our work. We address each major comment below and will revise the manuscript to incorporate the requested clarifications and additions.
read point-by-point responses
-
Referee: [Methods / Energy-parameter section] The central claim that vacancy-impurity interactions dominate inter-vacancy correlations rests on the specific numerical values chosen for the defect interaction energies. The manuscript must explicitly state whether these parameters were obtained from independent first-principles calculations, fitted to a subset of experiments with reported uncertainty, or selected by hand; without this, the dominance result risks being parameter-dependent rather than physically predictive. A sensitivity scan over plausible ranges of the vacancy-impurity binding versus vacancy-vacancy repulsion should be added.
Authors: We appreciate this clarification request. The interaction energies in our model are representative values drawn from independent first-principles DFT calculations reported in the literature for acceptor-doped ABO3 perovskites (e.g., SrTiO3 and related systems); they were not fitted to the experimental data presented here. We will explicitly state this origin and cite the relevant references in the revised Methods section. To address parameter dependence, we will add a sensitivity analysis that varies the vacancy-impurity binding energy and vacancy-vacancy repulsion over plausible ranges (approximately ±30% around the base values) and demonstrate that the dominance of vacancy-impurity interactions persists across this range, with only quantitative shifts in the reported trends. revision: yes
-
Referee: [Results / Oxidation and conductivity subsection] The reported increase in oxidation enthalpy and non-monotonic hole-concentration behavior with doping x are direct outputs of the chosen energies. These predictions should be compared quantitatively to measured oxidation enthalpies or conductivity data across a wider doping range, with explicit error bars, to test whether the model reproduces the experimental trends without post-hoc adjustment.
Authors: We agree that expanded quantitative validation is valuable. In the revised Results section we will include direct comparisons of our calculated oxidation enthalpies to experimental values from the literature for acceptor-doped ABO3 perovskites over a broader doping range (x = 0.01–0.20). Model predictions will be plotted against measured data, with error bars obtained from the standard deviation across independent Monte Carlo runs and with reported experimental uncertainties where available. We will also provide a quantitative discussion of the non-monotonic hole-concentration behavior in relation to published conductivity measurements, confirming that the trends align without any post-hoc parameter adjustment. revision: yes
Circularity Check
No significant circularity; central results are simulation outputs validated against external experiments
full rationale
The paper develops a statistical theory and Monte Carlo simulations whose outputs (vacancy-impurity dominance, doping dependence of interactions, oxidation enthalpy shifts) are generated from chosen but externally described 'realistic' energy parameters. These outputs are then compared to independent experimental data on local ion coordination and hole conductivity. No quoted equations or steps reduce by construction to fitted inputs, self-citations, or ansatzes; the parameter choice is presented as input rather than derived from the target claims. The derivation chain remains self-contained with external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- interaction energy parameters
axioms (1)
- domain assumption The developed statistical theory accurately describes defect thermodynamics and local coordination in ABO3 perovskites
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
under realistic energy parameters the interaction between oxygen vacancies and impurities generally has a greater impact... (abstract and §3.1–3.5)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Hamiltonian H0 = ∑ EVi nVi ... plus Hint = λ ∑<i,j> nVi nVj (Eqs. 1,3)
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]
C. Duan, J. Huang, N. Sullivan , R. O’Hayre, Proton-conducting oxides for energy conversion and storage, Appl. Phys. Rev. 7 (2020) 011314
work page 2020
-
[2]
A.V. Kasyanova, I.A. Zvonareva, N.A. Tarasova, L. Bi, D.A. Medvedev, Z. Shao, Electrolyte materials for protonic ceramic electrochemical cells: Main limitations and potential solutions, Mater. Rep.: Energy 2 (2022) 100158
work page 2022
-
[3]
L.A. Dunyushkina, S.A. Belyakov, N.M. Filatov, Proton -conducting alkaline earth hafnates: A review of manufacturing technologies, physicochemical properties and electrochemical performance, J. Eur. Ceram. Soc. 43 (2023) 6681-6698
work page 2023
- [4]
-
[5]
Kreuer, Proton-conducting oxides, Annu
K.D. Kreuer, Proton-conducting oxides, Annu. Rev. Mater. Res. 33 (2003) 333–359
work page 2003
-
[6]
Y. Yamazaki, F. Blanc, Y. Okuyama, L. B uannic, J.C. Lucio -Vega, C.P. Grey, S.M. Haile, Proton trapping in yttrium -doped barium zirconate , Nat. Mater. 12 (2013) 647–651
work page 2013
-
[7]
E. Gilardi, E. Fabbri, L. Bi, J.L.M. Rupp, T. Lippert, D. Pergolesi , E. Traversa, Effect of dopant–host ionic radii mismatch on acceptor-doped barium zirconate microstructure and proton conductivity, J. Phys. Chem. C 121 (2017) 9739−9747
work page 2017
-
[8]
D. Han, K. Shinoda, S. Sato, M. Majima , T. Uda, Correlation between electroconductive and structural properties of proton conductive a cceptor- doped barium zirconate, J. Mater. Chem. A 3 (2015) 1243–1250
work page 2015
- [9]
-
[10]
L.P. Putilov, N.A. Shevyrev, A.M. Mineev, A.S. Farlenkov, D.A. Medvedev , V.I. Tsidilkovski, Hydration of acceptor -doped BaSnO 3: Implications of the bound states of ionic defects, Acta Mater. 190 (2020) 70–80
work page 2020
-
[11]
I. Oika wa, H. Takamura, Correlation among oxygen vacancies, protonic defects, and the acceptor dopant in Sc -doped BaZrO3 studied by 45Sc nuclear magnetic resonance, Chem. Mater. 27 (2015) 6660–6667
work page 2015
-
[12]
M.E. Björketun, P.G. Sundell , G. Wahnström, Structure and thermodynamic stability of hydrogen interstitials in BaZrO 3 perovskite oxide from density functional calculations, Faraday Discuss. 134 (2007) 247–265
work page 2007
-
[13]
E. Bevillon, J. Hermet, G. Dezanneau , G. Geneste, How dopant size influences the protonic energy landscape in BaSn 1–xMxO3–x/2 (M = Ga, Sc, In, Y, Gd, La), J. Mater. Chem. A 2 (2014) 460–471
work page 2014
-
[14]
T.S. Bjørheim, E. Kotomin, J. Maier, Hydration entropy of BaZrO3 from first principles phonon calculations, J. Mater. Chem. A 3 (2015) 7639–7648
work page 2015
-
[15]
A. Løken, T.S. Bjørheim, R. Haugsrud, The pivotal role of the dopant choice on the thermodynamics of hydration and associations in proton conducting BaCe0.9X0.1O3– (X = Sc, Ga, Y, In, Gd and Er) , J. Mater. Chem. A 3 (2015) 23289–23298
work page 2015
-
[16]
Y. Yamazaki, A. Kuwabara , J. Hyodo, Y. Okuyama, C.A.J. Fisher , S.M. Haile, Oxygen affinity: The missing link enabling prediction of proton conductivities in doped barium zirconates , Chem. Mater. 32 (2020) 7292−7300
work page 2020
-
[17]
L.P. Putilov, V.I. Tsidilkovski, Impact of bound ionic defects on the hydration of acceptor-doped proton-conducting perovskites, Phys. Chem. Chem. Phys. 21 (2019) 6391–6406
work page 2019
-
[18]
L.P. Putilov , V.I. Tsidilkovski, Improving the performance of protonic ceramic fuel cells and electrolyzers: The role of acceptor impurities in oxide membranes, Energy Convers. Manage. 267 (2022) 115826
work page 2022
-
[19]
L.P. Putilov, V.I. Tsidilkovski, Proton conduction in acceptor-doped BaSnO3: The impact of the interaction between ionic defects and acceptor impurities , Materials 15 (2022) 4795
work page 2022
-
[20]
V.I. Tsidilkovski, M.Z. Uritsky, A.N. Varaksin , A.Y. Fishman, The influence of defects on proton diffusion in perovskites А IIВIV 1–xRIII xO3–δ: Monte Carlo study, Defect Diffus. Forum 258–260 (2006) 124–129
work page 2006
-
[21]
M.E. Björketun, P.G. Sundell , G. Wahnström, Effect of acceptor dopants on the proton mobility in BaZrO 3: A density functional investigation, Phys. Rev. B 76 (2007) 054307
work page 2007
-
[22]
C.Y.R. Vera, H. Ding, J. Urban -Klaehn, M. Li, Z. Zhao, F. Stewart, H. Tian, X. Liu, Y. Dong, J. Li, M. Zhou, H. Luo , D. Ding, Improving proton conductivity by navigating proton trapping in high scandium -doped barium zirconate electrolytes, Chem. Mater. 35 (2023) 5341–5352
work page 2023
-
[23]
L.P. Putilov , V.I. Tsidilkovski, The role of deep acceptor centers in the oxidation of acceptor -doped wide-band-gap perovskites ABO 3, J. Solid State Chem. 247 (2017) 147–155. 27
work page 2017
-
[24]
F.M. Draber, J.R. Denninger, P.C. Müller, I.K. Sommerfeld , M. Martin, The impact of nanoscale percolation in yttrium -doped BaZrO3 on the oxygen ion and proton conductivities: A density functional theory and kinetic Monte Carlo study, Adv. Energy Sustainability Res. 3 (2022) 2200007
work page 2022
-
[25]
A. Lindman, P. Erhart , G. Wahnström, Polaronic contributions to oxidation and hole conductivity in acceptor -doped BaZrO 3, Phys. Rev. B 94 (2016) 075204
work page 2016
-
[26]
P.G. Sundell, M.E. Bjorketun , G. Wahnstrom, Thermodynamics of doping and vacancy formation in BaZrO 3 perovskite oxide from density functional calculations, Phys. Rev. B 73 (2006) 104112
work page 2006
-
[27]
E. Bevillon, G. Dezanneau , G. Geneste, Oxygen incorporation in acceptor - doped perovskites, Phys. Rev. B 83 (2011) 174101
work page 2011
- [28]
-
[29]
L. Buannic, F. Blanc, I. Hung, Z. Gan, C.P. Grey, Probing the local structures and protonic conduction pathways in scandium substituted BaZrO 3 by multinuclear solid-state NMR spectroscopy, J. Mater. Chem. 20 (2010) 6322– 6332
work page 2010
-
[30]
S. Kasamatsu, O. Sugino, T. Ogawa , A. Kuwabara, Dopant arrangements in Y-doped BaZrO 3 under processing conditions and their impact on proton conduction: a large -scale first-principles thermodynamics study , J. Mater. Chem. A 8 (2020) 12674–12686
work page 2020
-
[31]
K. Hoshino, S. Kasamatsu, J. Hyodo, K. Yamamoto, H. Setoyama, T. Okajima, Y. Yamazaki, Probing local environments of oxygen vacancies responsible for hydration in Sc -doped barium zirconates at elevated temperatures: In situ X ‑ray absorption spectroscopy, thermogravimetry, and active learning ab initio replica exchange Monte Carlo simulations , Chem. Mat...
work page 2023
- [32]
- [33]
-
[34]
I.A. Zvonareva, A.M. Mineev, N.A. Tarasova, X. -Z. Fu, D.A. Medvedev, High-temperature transport properties of BaSn1−xScxO3−δ ceramic materials as promising electrolytes for protonic ceramic fuel cells , J. Adv. Ceram. 11 (2022) 1131–1143
work page 2022
-
[35]
C. De Dominicis, I. Giardina, Random fields and spin glasses. A field theory approach, Cambridge University Press, Cambridge, 2006
work page 2006
-
[36]
L.D. Landau , E.M. Lifshitz, Statistical Physics. Vol. 5, Butterworth - Heinemann, Oxford, 1980. 28
work page 1980
-
[37]
M.Z. Uritsky, V.I. Tsidilkovski, Role of an acceptor impurity in the proton transfer in proton-conducting oxides, Phys. Solid State 56 (2014) 2173–2180
work page 2014
-
[38]
M.Z. Uritskii, Distribution of oxygen vacancies in an external electric field and oxygen -ion conductivity in acceptor -doped fluorite -structured ionic conductors, Phys. Solid State 63 (2021) 937–941
work page 2021
-
[39]
Cowley, An approximate theory of order i n alloys, Phys
J.M. Cowley, An approximate theory of order i n alloys, Phys. Rev. 77 (1950) 669–675
work page 1950
- [40]
- [41]
- [42]
-
[43]
G. Geneste, B. Amadon, M. Torrent, G. Dezanneau , DFT+U study of self - trapping, trapping, and mobility of oxygen -type hole polarons in barium stannate, Phys. Rev. B 96 (2017) 134123
work page 2017
-
[44]
V.I. Tsidilkovski, L.P. Putilov, The role of deep acceptor levels in hydration and transport processes in BaZr 1−xYxO3–δ and related materials, J. Solid State Electrochem. 20 (2016) 629–643
work page 2016
-
[45]
L.P. Putilov , A.K. Demin, V.I. Tsidilkovski , P. Tsiakaras, Theoretical modeling of the gas humidification effect on the characteristics of proton ceramic fuel cells, Appl. Energy 242 (2019) 1448–1459
work page 2019
-
[46]
L.P. Putilov, V.I. Tsidilkovski, A.K. Demin, Revealing the e ffect of the cell voltage and external conditions on the characteristics of protonic ceramic fuel cells, J. Mater. Chem. A 8 (2020) 12641–12656
work page 2020
-
[47]
I.A. Zvonareva, G.N. Starostin, M.T. Akopian, G.K. Vdovin, X. -Z. Fu, D.A. Medvedev, Ionic and electronic transport of dense Y -doped barium stannate ceramics for high -temperature applications, J . Power Sources 565 (2023) 232883
work page 2023
-
[48]
G. Geneste, G. Dezanneau, Competition between elastic and chemical effects in the doping, defect association, and hydration of barium stannate , Solid State Ionics 308 (2017) 121–132
work page 2017
-
[49]
V.I. Tsidilkovski, I.A. Leonidov , A .A. Lakhtin, V.A. Mezrin, High- temperature equilibrium between high -Tc oxide and gas phase: The role of electrons and holes, Phys Status Solidi B 163 (1991) 231–240
work page 1991
-
[50]
V.I. Tsidilkovski, I.A. Leonidov, A.A. Lakhtin, V.A. Mezrin, The role of the electron-hole system in the thermodynamics of YBa2Cu3O7––gas equilibrium, Phys. Status Solidi B 168 (1991) 233–244
work page 1991
-
[51]
D. Han, K. Shinoda, T. Uda, Dopant Site Occupancy and Chemical Expansion in Rare Earth -Doped Barium Zirconate, J. Am. Ceram. Soc. 97 (2014) 643–650. 29
work page 2014
-
[52]
S. I mashuku, T. Uda, Y. Nose, G. Taniguchi, Y. Ito, Y. Awakura, Dependence of Dopant Cations on Microstructure and Proton Conductivity of Barium Zirconate, J. Electrochem. Soc. 156 (2009) B1-B8
work page 2009
-
[53]
L. Casillas-Trujillo, D.A. Andersson, B. Dorado, M. Nikl, K.E. Si ckafus, K.J. McClellan, C.R. Stanek, Intrinsic defects, nonstoichiometry, and aliovalent doping of A 2+B4+O3 perovskite scintillators, Phys. Status Solidi B 251 (2014) 2279–2286
work page 2014
-
[54]
L.P. Putilov, M.Z. Uritsky, V.I. Tsidilkovski, Oxygen-vacancy transport in acceptor-doped perovskites: effective medium approximation approach , Electrochem. Mater. Technol. 4 (2025) 20254055
work page 2025
-
[55]
L.P. Putilov, V.I. Tsidilkovski, Oxide-ion thermopower of acceptor -doped perovskites: The impact of inter -defect interactions , Electrochem. Mat er. Technol. 4 (2025) 20254056
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.