pith. sign in

arxiv: 2607.01779 · v1 · pith:QUFE7F6Gnew · submitted 2026-07-02 · ❄️ cond-mat.mtrl-sci

Density functional study of native point defects in CaO

Yunhwa Jo , Minseok Choi This is my paper

Pith reviewed 2026-07-03 10:11 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords CaOnative point defectsoxygen vacancycalcium vacancydensity functional theoryformation energyoptical transitionvacancy complex
0
0 comments X

The pith

Oxygen vacancies dominate CaO under O-poor conditions while calcium vacancies dominate under O-rich conditions, and their stable complexes account for observed optical peaks.

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

The paper performs first-principles density-functional calculations to determine the formation, migration, and optical properties of native point defects in calcium oxide. It shows that the preferred vacancy type switches with the oxygen chemical potential during growth. The calculations further establish that complexes formed by these vacancies possess binding energies and migration barriers that allow them to remain stable through high-temperature annealing. Franck-Condon optical transition energies computed for the defects and complexes align with multiple experimentally reported absorption and emission features.

Core claim

First-principles density-functional calculations show that oxygen vacancies are favored under O-poor conditions, whereas calcium vacancies dominate under O-rich conditions. Calculated migration barriers and binding energies indicate that vacancy complexes are thermodynamically stable and can survive high-temperature annealing. Optical transition energies, evaluated using the Franck-Condon framework, suggest that several experimentally observed absorption and emission peaks can be attributed to negatively charged vacancy complexes as well as isolated oxygen vacancies.

What carries the argument

Density-functional calculations of defect formation energies under varying chemical potentials, together with migration barriers, binding energies, and Franck-Condon optical transition energies.

If this is right

  • Vacancy complexes remain thermodynamically stable after high-temperature annealing.
  • Negatively charged vacancy complexes contribute to observed optical absorption and emission spectra.
  • Isolated oxygen vacancies produce distinct optical transitions under O-poor growth conditions.
  • The dominant defect species can be selected by controlling the oxygen chemical potential during synthesis.

Where Pith is reading between the lines

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

  • Adjusting oxygen partial pressure during crystal growth offers a route to select between vacancy types for targeted electronic or optical behavior.
  • Accumulation of stable vacancy complexes may influence long-term material stability in high-temperature environments even after processing.
  • Extension of the same computational approach to related alkaline-earth oxides could identify common patterns in defect complex stability.

Load-bearing premise

The density-functional calculations correctly predict formation energies, migration barriers, binding energies, and optical transition energies without large errors from the chosen functional, supercell size, or other technical approximations.

What would settle it

Experimental measurement of defect concentrations or optical absorption/emission peaks under controlled oxygen-rich versus oxygen-poor conditions that contradict the calculated dominance switch or the assigned transition energies.

Figures

Figures reproduced from arXiv: 2607.01779 by Minseok Choi, Yunhwa Jo.

Figure 1
Figure 1. Figure 1: FIG. 1: Formation energies of native point defects as a function of the Fermi level at (a) O-poor limit, (b) Ca-poor limit, and [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Local atomic structure for O vacancies in CaO in [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: (a) Local atomic structure and (b) side veiw for Ca [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Local atomic structures for O [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: (a) Local structure for Ca [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: (a) Optical absorption and (b) emission associated with native defects in CaO. The depicted levels are not zero-phonon [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗
read the original abstract

We investigate the structural, electronic, and optical properties of native point defects in CaO using first-principles density-functional calculations. Oxygen vacancies are favored under O-poor conditions, whereas calcium vacancies dominate under O-rich conditions. Calculated migration barriers and binding energies indicate that vacancy complexes are thermodynamically stable and can survive high-temperature annealing. Optical transition energies, evaluated using the Franck-Condon framework, suggest that several experimentally observed absorption and emission peaks can be attributed to negatively charged vacancy complexes as well as isolated oxygen vacancies.

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

3 major / 1 minor

Summary. The manuscript reports first-principles density-functional calculations of native point defects in CaO. It concludes that oxygen vacancies are favored under O-poor conditions while calcium vacancies dominate under O-rich conditions. Migration barriers and binding energies indicate that vacancy complexes are thermodynamically stable and survive high-temperature annealing. Optical transition energies computed in the Franck-Condon framework are used to attribute several experimental absorption and emission peaks to negatively charged vacancy complexes as well as isolated oxygen vacancies.

Significance. If the computed formation energies, barriers, binding energies, and optical levels prove accurate, the work would provide a useful link between defect thermodynamics, stability after annealing, and observed optical spectra in CaO, an ionic oxide of technological interest. The explicit connection drawn to experimental peaks via the Franck-Condon model is a positive feature when the underlying energies are reliable.

major comments (3)
  1. [Computational Methods] The central claims on thermodynamic preference (O-poor vs. O-rich) and on survival of complexes after annealing rest on the absolute accuracy of formation energies and migration barriers, yet the manuscript supplies no hybrid-functional benchmarks, no explicit supercell-size extrapolation, and no finite-size correction details for charged defects. These omissions are load-bearing because semi-local functionals are known to underestimate the CaO gap by ~2 eV and to misplace defect levels.
  2. [Results and Discussion] The assignment of specific experimental absorption/emission peaks to vacancy complexes and isolated oxygen vacancies depends on the Franck-Condon optical transition energies being within ~0.3–0.5 eV of experiment. No comparison table or quantitative error analysis against measured peak positions is presented, leaving the attribution unsupported if the DFT levels carry the typical GGA error.
  3. [Results] No convergence data, dielectric-constant values used for image-charge corrections, or tests of the chosen supercell size appear for the formation-energy calculations. In an ionic material such as CaO these technical choices directly affect the reported stability ordering between oxygen and calcium vacancies.
minor comments (1)
  1. [Abstract] The abstract states the methods and conclusions but contains no numerical values, error bars, or key energies, which reduces its utility as a standalone summary.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough review and valuable suggestions. We address the major comments point by point below, indicating the revisions we plan to make to strengthen the manuscript.

read point-by-point responses

  1. Referee: [Computational Methods] The central claims on thermodynamic preference (O-poor vs. O-rich) and on survival of complexes after annealing rest on the absolute accuracy of formation energies and migration barriers, yet the manuscript supplies no hybrid-functional benchmarks, no explicit supercell-size extrapolation, and no finite-size correction details for charged defects. These omissions are load-bearing because semi-local functionals are known to underestimate the CaO gap by ~2 eV and to misplace defect levels.

    Authors: We agree that additional technical details would improve the clarity and robustness of the results. In the revised manuscript, we will provide the finite-size correction details used for charged defects and include explicit tests of supercell-size convergence for the formation energies. We will also report the dielectric constant values employed. Regarding hybrid-functional benchmarks, performing such calculations is beyond the current scope of this work; however, we will add a discussion section addressing the known limitations of semi-local functionals for defect levels in wide-gap oxides like CaO and justify the use of our approach based on consistency with experimental trends. revision: partial

  2. Referee: [Results and Discussion] The assignment of specific experimental absorption/emission peaks to vacancy complexes and isolated oxygen vacancies depends on the Franck-Condon optical transition energies being within ~0.3–0.5 eV of experiment. No comparison table or quantitative error analysis against measured peak positions is presented, leaving the attribution unsupported if the DFT levels carry the typical GGA error.

    Authors: We will include a new table in the revised manuscript that directly compares the computed optical transition energies with the experimental peak positions. This table will be accompanied by a quantitative discussion of the agreement and an assessment of possible errors arising from the functional choice. revision: yes

  3. Referee: [Results] No convergence data, dielectric-constant values used for image-charge corrections, or tests of the chosen supercell size appear for the formation-energy calculations. In an ionic material such as CaO these technical choices directly affect the reported stability ordering between oxygen and calcium vacancies.

    Authors: As noted in response to the first comment, we will add the requested convergence data, dielectric constants, and supercell size tests to the revised manuscript to demonstrate that the stability ordering is robust. revision: yes

Circularity Check

0 steps flagged

No circularity: forward DFT computations independent of target observables

full rationale

The paper reports standard first-principles DFT calculations of formation energies, migration barriers, binding energies, and Franck-Condon optical transition energies for native defects in CaO. These quantities are obtained directly from total-energy differences and electronic-structure evaluations within the chosen functional and supercell setup; none are defined by or fitted to the experimental absorption/emission peaks that the results are later compared against. No self-citation chain, ansatz smuggling, or renaming of known results is invoked as a load-bearing step in the derivation. The reported thermodynamic preferences and peak attributions therefore remain independent forward predictions rather than tautological restatements of inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the unstated technical validity of standard DFT for defect formation and optical energies in this ionic oxide; no free parameters or new entities are mentioned in the abstract.

axioms (1)
  • domain assumption Density functional theory approximations are sufficiently accurate for formation energies, migration barriers, and Franck-Condon optical transitions of native point defects in CaO.
    The abstract states that first-principles density-functional calculations were used to obtain all reported quantities.

pith-pipeline@v0.9.1-grok · 5603 in / 1231 out tokens · 46096 ms · 2026-07-03T10:11:42.646539+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

37 extracted references · 4 canonical work pages · 1 internal anchor

  1. [1]

    VO thus acts as a deep donor, which agrees with previous GGA [3] and HSE study [8]

    Vacancies As shown in Fig 1, O vacancies introduce two charge- state transition levels in the band -gap: a (+/0) level at 0.95 eV and (2+/+) level at 2.65 eV below the CBM. VO thus acts as a deep donor, which agrees with previous GGA [3] and HSE study [8]. The charge density of the defect states in the band-gap and the local atomic relaxations are shown i...

  2. [2]

    Interstitials The oxygen interstitial, O i, has modest formation en- ergies at O-rich limit (Fig. 1(b)). However, the defects may form only when the Fermi level is close to the con- duction band. The defect has a (0/2–) level at 2.29 eV above the VBM, which indicates that the interstitials are deep acceptors. Oi can have two distinct configurations. In th...

  3. [3]

    The defect has a (0/2–) level at 0.05 eV, (2–/3–) level at 3.42 eV, and (3–/4–) level at 3.71 eV above the VBM

    Antisite-like defects Although O Ca has low formation energies under Ca- poor condition (Fig 1 (b)), the defect may be not likely to form in experiment since its formation energies are quite high as seen in Fig 1 (c). The defect has a (0/2–) level at 0.05 eV, (2–/3–) level at 3.42 eV, and (3–/4–) level at 3.71 eV above the VBM. OCa shows significant offce...

  4. [4]

    A recent HSE study em- ploying a tuned mixing parameter reported some unusual aspects of defect charges in CaO

    Comparison with hybrid functional To further examine the reliability of the predicted SCAN results, we compared the SCAN results with hybrid-functional calculations. A recent HSE study em- ploying a tuned mixing parameter reported some unusual aspects of defect charges in CaO. For example, although Ca vacancies are well known acceptors, the HSE predicts t...

    2024

  5. [5]

    Whited, C

    R. Whited, C. J. Flaten, and W. Walker, Solid State Communications13, 1903 (1973)

    1903

  6. [6]

    Albuquerque and M

    E. Albuquerque and M. Vasconcelos, inJournal of Physics: Conference Series(2008), vol. 100, p. 042006

    2008

  7. [7]

    K. Wu, T. Liu, R. Sun, J. Song, and C. Shi, The Euro- pean Physical Journal D74, 209 (2020)

    2020

  8. [8]

    M. L. Franckel, S. Falletta, W. Chen, and A. Pasquarello, Physical Review B112, 104101 (2025)

    2025

  9. [9]

    Kanai, F

    S. Kanai, F. J. Heremans, H. Seo, G. Wolfowicz, C. P. Anderson, S. E. Sullivan, M. Onizhuk, G. Galli, D. D. Awschalom, and H. Ohno, Proceedings of the National Academy of Sciences119, e2121808119 (2022)

    2022

  10. [10]

    Davidsson, M

    J. Davidsson, M. Onizhuk, C. Vorwerk, and G. Galli, Nature Communications15, 4812 (2024)

    2024

  11. [11]

    Y. Chen, V. Orera, R. Gonzalez, R. Williams, G. Williams, G. Rosenblatt, and G. Pogatshnik, Phys- ical Review B42, 1410 (1990)

    1990

  12. [12]

    Yuan and G

    Z. Yuan and G. Hautier, Applied Physics Letters124 (2024)

    2024

  13. [13]

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, The Journal of chemical physics118, 8207 (2003)

    2003

  14. [14]

    A. V. Krukau, O. A. Vydrov, A. F. Izmaylov, and G. E. Scuseria, The Journal of chemical physics125(2006)

    2006

  15. [15]

    Elfimov, S

    I. Elfimov, S. Yunoki, and G. Sawatzky, Physical review letters89, 216403 (2002)

    2002

  16. [16]

    Osorio-Guill´ en, S

    J. Osorio-Guill´ en, S. Lany, S. Barabash, and A. Zunger, Physical review letters96, 107203 (2006)

    2006

  17. [17]

    P. E. Bl¨ ochl, Physical review B50, 17953 (1994)

    1994

  18. [18]

    Kresse and J

    G. Kresse and J. Hafner, Physical Review B48, 13115 (1993)

    1993

  19. [19]

    J. Sun, A. Ruzsinszky, and J. P. Perdew, Physical review letters115, 036402 (2015)

    2015

  20. [20]

    Freysoldt, B

    C. Freysoldt, B. Grabowski, T. Hickel, J. Neugebauer, G. Kresse, A. Janotti, and C. G. Van de Walle, Reviews of modern physics86, 253 (2014)

    2014

  21. [21]

    Kumagai, Physical Review B107, L220101 (2023)

    Y. Kumagai, Physical Review B107, L220101 (2023)

    2023

  22. [22]

    Henkelman, B

    G. Henkelman, B. P. Uberuaga, and H. J´ onsson, The Journal of chemical physics113, 9901 (2000)

    2000

  23. [23]

    G. H. Vineyard, Journal of Physics and Chemistry of Solids3, 121 (1957)

    1957

  24. [24]

    Ohsato, T

    H. Ohsato, T. Sugimura, and K. Kageyama, Journal of Crystal Growth48, 459 (1980). 9

    1980

  25. [25]

    D. K. Smith and H. Leider, Applied Crystallography1, 246 (1968)

    1968

  26. [26]

    Wriedt, Bulletin of Alloy phase diagrams6, 337 (1985)

    H. Wriedt, Bulletin of Alloy phase diagrams6, 337 (1985)

    1985

  27. [27]

    M. W. Chase, J. of Physical and Chemical Refernce Data p. 1529 (1998)

    1998

  28. [28]

    Momma and F

    K. Momma and F. Izumi, Applied Crystallography44, 1272 (2011)

    2011

  29. [29]

    Y. Cui, Y. Pan, L. Pascua, H. Qiu, C. Stiehler, H. Kuh- lenbeck, N. Nilius, and H.-J. Freund, Physical Review B 91, 035418 (2015)

    2015

  30. [30]

    Janotti and C

    A. Janotti and C. G. Van de Walle, Phys. Rev. B 76, 165202 (2007), URLhttps://link.aps.org/doi/ 10.1103/PhysRevB.76.165202

    work page doi:10.1103/physrevb.76.165202 2007

  31. [31]

    M. Choi, M. E. Turiansky, B. Zhao, J. D. Thompson, and C. G. V. de Walle,Native defects and erbium impurities in cawo4(2026), 2605.24809, URLhttps://arxiv.org/ abs/2605.24809

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  32. [32]

    Alkauskas, P

    A. Alkauskas, P. Broqvist, and A. Pasquarello, Phys. Rev. Lett.101, 046405 (2008)

    2008

  33. [33]

    Alkauskas and A

    A. Alkauskas and A. Pasquarello, Phys. Rev. B84, 125206 (2011)

    2011

  34. [34]

    Henderson, Y

    B. Henderson, Y. Chen, and W. A. Sibley, Phys. Rev. B6, 4060 (1972), URLhttps://link.aps.org/doi/10. 1103/PhysRevB.6.4060

    1972

  35. [35]

    L. S. Welch, A. E. Hughes, and G. P. Summers, Journal of Physics C: Solid State Physics13, 1791 (1980), URL https://doi.org/10.1088/0022-3719/13/9/024

    work page doi:10.1088/0022-3719/13/9/024 1980

  36. [36]

    Subramanian, R

    M. Subramanian, R. Shannon, B. Chai, M. Abraham, and M. Wintersgill, Physics and chemistry of minerals 16, 741 (1989)

    1989

  37. [37]

    J. E. Medvedeva, E. N. Teasley, and M. D. Hoffman, Phys. Rev. B76, 155107 (2007), URLhttps://link. aps.org/doi/10.1103/PhysRevB.76.155107

    work page doi:10.1103/physrevb.76.155107 2007