Pith. sign in

REVIEW 3 major objections 5 minor 38 references

GW electronic structure calculations of Co doped ZnO

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper argues that the 1.74–1.88 eV luminescence in Co-doped ZnO comes from a cobalt atom on a zinc site with an oxygen interstitial nearby, and that the other common neutral point defects tested cannot produce that transition.

desk verdict Useful follow-up with solid neutral-defect formation energies, but the confirm/rule-out claim outruns the calculation because no optical transition energies are computed and the defect scan is neutral-only. read the letter →

arxiv 1908.04679 v1 pith:CH5RI4O2 submitted 2019-08-13 cond-mat.mtrl-sci

Dennis Franke , Michael Lorke , Thomas Frauenheim , A. L. da Rosa This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords cobalt-dopedzincoxideoxygeninterstitialdefectintra-3dluminescenceGWapproximationdensityfunctionaltheoryformationenergydilutedmagneticsemiconductorsorbital-resolvedbandstructure
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper sets out to identify the point defect responsible for the red intra-3d luminescence observed in cobalt-doped ZnO after ion implantation. By computing orbital-resolved band structures with density-functional theory and many-body GW corrections, it concludes that the only neutral defect among those tested that can support the observed d-d transition is a cobalt atom substituting a zinc atom with an oxygen interstitial nearby. It also concludes that isolated cobalt substitution, cobalt paired with a zinc vacancy, cobalt paired with an oxygen vacancy, and cobalt paired with a zinc interstitial can all be ruled out. The practical stake is that the optical center is now specified at the atomic level, which matters for using Co-doped ZnO as a tunable infrared emitter or a single-photon source.

What carries the argument

The central object is the neutral defect complex CoZn+Oint, a substitutional cobalt atom on a zinc site with an oxygen interstitial in its neighborhood. The argument is carried by orbital-resolved density of states from GW-corrected density-functional calculations, by the crystal-field splitting between the e and t2 groups of Co-3d states, and by the formation-energy formula of Eq. (1). The interstitial oxygen is the load-bearing feature: it reduces the e–t2 splitting, creates empty d states in the band gap, and produces the Co-d/Zn-4s overlap that the luminescence and the proposed ferromagnetism both require.

What would settle it

Grow or anneal Co-doped ZnO so that only the neutral CoZn+Oint complex should be present, measure photoluminescence, and look for the 1.74–1.88 eV line: if the line is absent, the assignment fails. Alternatively, compute the formation energies and optical transitions of charged versions of CoZn+Oint, CoZn+VZn, CoZn+VO, and CoZn+Znint; if any charged defect has a lower formation energy and a matching transition energy, the neutral-only conclusion is falsified.

Watch

Extended reading notes

Core claim

The authors find that adding the oxygen interstitial to CoZn changes the Co-3d electronic structure qualitatively: the e–t2 crystal-field splitting drops from 3.2 eV to 2.4 eV, the Co-d states hybridize more strongly with Zn-4s and O-p states, and empty d states appear in the gap so that the intra-3d transition at roughly 1.88 eV can occur. In contrast, CoZn alone has no empty d states in the gap, and the complexes CoZn+VZn, CoZn+VO, and CoZn+Znint either lack intra-gap Co states or show no clear optical transition. Formation energies computed for neutral defects show CoZn+Oint is among the stable complexes under oxygen-rich conditions, and because the isolated oxygen interstitial has a high formation energy but a low diffusion barrier, the authors argue it forms a stable complex once created. They also report a magnetic moment of 2.8 µB for the complex and enhanced Co-d/Zn-4s overlap, and suggest this could favor ferromagnetism.

Load-bearing premise

The calculations only consider neutral defects and complexes, so if the luminescence actually comes from a charged version of one of these defects, the confirmation of CoZn+Oint and the exclusion of the others would not follow.

Editorial extensions

If this is right

  • If the assignment is correct, the 1.74–1.88 eV emission line is a fingerprint of the specific CoZn+Oint complex, so its intensity should track oxygen-rich preparation conditions.
  • The other neutral point defects tested can be excluded as sources of the observed intra-3d transition, simplifying the interpretation of Co-doped ZnO spectra.
  • Because CoZn+Oint is stable under oxygen-rich conditions and oxygen interstitials diffuse quickly, annealing in oxygen-rich atmospheres should favor the luminescence.
  • The enhanced overlap between Co-d and Zn-4s states in the complex suggests a route toward carrier-mediated ferromagnetism in Co-doped ZnO, triggered by oxygen excess.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, if the complex is the emitter, its exact transition energy should respond to strain or isotopic substitution of the interstitial oxygen, offering a way to experimentally fingerprint the defect.
  • The same logic could be tested in other transition-metal-doped oxides, where oxygen-interstitial complexes may serve as generic design centers for infrared emission.
  • Because only neutral charge states were tested, a charged version of CoZn+Oint or of a competing defect could still match the observed luminescence; a GW study of charged complexes would settle that.
  • The ferromagnetism suggestion can be tested by measuring the magnetization of Co-implanted ZnO annealed under oxygen-rich versus oxygen-poor conditions and looking for the predicted difference.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper uses PBE structural relaxations and PBE+GW0 electronic-structure calculations for 72-atom supercells of Co-doped ZnO containing CoZn, CoZn+Oint, CoZn+VZn, CoZn+VO, and CoZn+Znint defects. It compares orbital-resolved densities of states and band-projected charge densities with experimental intra-3d luminescence at 1.74-1.88 eV, reports PBE formation energies under O-rich and O-poor conditions, and concludes that the CoZn+Oint complex is the likely origin of the luminescence, that the other tested neutral defects can be ruled out, and that this complex may promote ferromagnetism. The central comparison is made at the level of one-particle eigenvalues; no optical transition energies, oscillator strengths, or charged-defect states are computed.

Significance. The manuscript addresses a concrete and timely materials problem: identifying the microscopic defect responsible for Co-related infrared emission in ZnO. Its systematic treatment of several defect complexes, tabulated formation energies, and orbital-decomposed DOS and charge-density plots are useful and provide a clear qualitative picture of defect-induced in-gap states. If the missing optical and charged-defect analyses were supplied, the identification would be significant for Co-implanted ZnO optoelectronics and for proposals of carrier-mediated ferromagnetism. At present, however, the quantitative 'confirm/rule-out' claims are stronger than what the reported calculations can support, so the significance is contingent on the additional evidence rather than established by the paper.

major comments (3)
  1. [Fig. 2 and the paragraph following it] The opening claim that the e-t2 splitting of about 2.4 eV for CoZn+Oint is 'in very good agreement with the experimentally observed luminescence signatures at 1.88 eV and 2.02 eV' is not derived from the reported data. A GW quasiparticle eigenvalue difference is not an optical excitation energy, especially for a localized d-d transition where self-energy and excitonic corrections do not cancel in a single-particle picture; the manuscript neither computes transition matrix elements nor performs a BSE or multiplet calculation. The text should either present an actual optical excitation calculation or explicitly downgrade the conclusion from 'confirm' to 'consistent with the earlier assignment'.
  2. [Eq. (1), Table I, and Fig. 5] The 'rule out' conclusion is overbroad because the study is restricted to neutral defect complexes. Equation (1) defines formation energies for neutral defects only, Table I lists only neutral complexes, and the DOS in Fig. 5 is computed only for neutral CoZn+VZn, CoZn+VO, and CoZn+Znint. Charged versions of these complexes and of isolated Oi, VO, and Zni are never considered. If the observed luminescence involves a charged defect, neither the confirmation of CoZn+Oint nor the exclusion of the alternatives follows; the conclusions should be restricted to neutral defects or complemented by charged-defect calculations.
  3. [Table I and the formation-energy discussion] The thermodynamic argument for CoZn+Oint being 'likely' is weakened by the PBE formation energies in Table I: under O-rich conditions CoZn+VZn (0.25 eV) is much more stable than CoZn+Oint (1.83 eV), and the distinction is made on electronic-structure grounds rather than on stability. Because PBE formation energies can shift by several tenths of an electronvolt with hybrid or GW total energies, the paper should explicitly state that the identification rests on the electronic-structure analysis and should discuss the implications of the competing CoZn+VZn complex for the proposed growth scenario.
minor comments (5)
  1. [Fig. 5 caption] The caption for Fig. 5 lists only panels (a) and (b), although the figure includes a third panel (c) for CoZn+Znint; the caption should describe all panels.
  2. [Paragraph following Fig. 1] The text refers to the 'ZnO octahedral crystal field' in wurtzite ZnO, but the Zn site is tetrahedrally coordinated and the e/t2 splitting described is the tetrahedral-field splitting; please correct the terminology.
  3. [Methods and figure captions] The method is called 'GW calculations' in the text while the figure captions specify 'PBE+GW0'; the authors should state whether a one-shot G0W0 calculation was used and provide convergence parameters for the GW step, including the number of unoccupied bands, k-point sampling, and frequency-grid parameters.
  4. [Reference list] The reference list contains duplicates: Refs. [3] and [8] are the same Sarsari et al. article, Refs. [30] and [33] are the same Lany and Zunger article, and Refs. [31] and [39] are the same Janotti and van de Walle article; these should be consolidated.
  5. [Throughout] There are numerous typographical errors, including 'inden-tified', 'interstital', 'fomation', 'cleary', and 'intestitial'; a careful proofreading pass is needed.

Circularity Check

1 steps flagged · score 4.0 of 10

Partially circular: the confirmation of CoZn+Oint leans on the authors' own prior assignment [1], while the rule-out and formation-energy calculations are independent.

  1. self citation load bearing [Introduction and Section 3 (DOS discussion, after Fig. 2)]
    "As reported previously, these results are in very good agreement with the experimentally observed luminescence signatures at 1.88 eV and 2.02 eV [1, 18, 19]."

    The paper's stated confirmation target is the defect assignment 'identified [1]' in the introduction. No optical transition energy, oscillator strength, or lifetime is computed; the only quantitative link offered for CoZn+Oint is a 2.4 eV e-t2 eigenvalue splitting that is then labeled as 'in very good agreement' with 1.88/2.02 eV luminescence solely by reference to the authors' prior paper [1]. Therefore the 'confirm' step is not an independent prediction: it takes the CoZn+Oint assignment as input from [1] and re-asserts it. The rule-out of CoZn+VZn, CoZn+VO, CoZn+Znint and the formation-energy table are new and do not reduce to [1].

full rationale

The paper is not wholly circular: it contains new, self-contained calculations for competing neutral defect complexes, a formation-energy analysis using the standard expression Ef = Etot_defect - Etot_bulk - sum_i n_i mu_i, and comparisons with independent experimental references [18,19]. However, the load-bearing 'confirmation' that CoZn+Oint is the source of the observed intra-3d luminescence is not derived from a computed optical transition; it rests on the authors' own prior assignment [1], reinforced by an eigenvalue-gap comparison that is asserted to agree with experiment 'as reported previously [1]'. That specific step reduces the central confirmation to a self-citation. The neutral-only charge-state restriction limits the scope of the rule-out but is a correctness concern rather than a separate circularity. Overall, the central claim has independent supporting content in the defect scan and formation energies, so a score of 4 is appropriate rather than a higher score indicating full reduction.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters are fitted; the only real inputs are the PBE and G0W0 approximations, standard chemical potentials, and the experimental PL energies used as benchmarks. The central claim rests on the assumption that neutral PBE+G0W0 defect levels can be directly compared with measured d-d luminescence, and that charged defects are irrelevant.

assumptions (4)
  • domain assumption PBE-relaxed geometries and PBE+G0W0 quasiparticle DOS provide a reliable description of defect levels in Co-doped ZnO.
    The conclusions about gap states and d-d splittings depend entirely on this approximation; no validation against higher-level methods or direct transition calculations is given.
  • domain assumption Only neutral charge states of the defects are relevant to the observed luminescence.
    Charged defects are not considered, yet many ZnO defects are stable in charged states; this assumption bounds the rule-out claim.
  • domain assumption The observed luminescence peaks at 1.88 and 2.02 eV can be identified with the computed e-t2 splitting of the Co d states rather than with a full many-body d-d excitation.
    The paper cites agreement with experiment without computing optical transition energies or matrix elements.
  • domain assumption The standard Van de Walle-Neugebauer formation-energy formalism applies to the defect complexes.
    The formation-energy results rely on this standard framework, including the chosen chemical potential limits.

how reviews work

0 comments
Cite this review

Pith. "Pith review of GW electronic structure calculations of Co doped ZnO." pith.science (2026). https://pith.science/paper/CH5RI4O2

@misc{pith2026190804679,
  author       = {Pith},
  title        = {Pith review of: GW electronic structure calculations of Co doped ZnO},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/CH5RI4O2}},
  note         = {Machine review of arXiv:1908.04679}
}
abstract

Recently the point defect responsible for the intra-3$d$ luminescence of cobalt in doped ZnO samples has been indentified\,\cite{pssb2019}. In this work we further extend our investigation to other point defects in Co-doped ZnO. We use density-functional theory and GW calculations to determined the orbital-resolved band structure of cobalt doped zinc oxide (ZnO). We show that mainly O-p and Co-d orbitals take part in the process and confirm that an oxygen interstitial nearby a cobalt atom is a likely defect to occur in Co-implanted ZnO samples. We also rule out that other common point defects in ZnO can be responsible for the observed intra-3$d$ transition. Finally, we suggest that defect complexes involving oxygen interstitials could be used to promote ferromagnetism in cobalt doped ZnO samples.

Figures

Figures reproduced from arXiv: 1908.04679 by the authors.

Figure 1
Figure 1. FIG. 1. Atomic structure around the a) Co [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Orbital projected density of states calculated with [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Band projected charge density for the Co [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Atomic structure around the a) Co [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Total and atom projected density of states for a) Co [PITH_FULL_IMAGE:figures/full_fig_p013_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

38 extracted references · 32 canonical work pages

  1. [1]

    Röder, M

    R. Röder, M. Zapf, S. Geburt, D. Franke, M. Lorke, A. L. da R osa, T. Frauenheim, and C. Ronning, phys. stat. sol. (b) 256, 1800604 (2019)

    work page 2019

  2. [2]

    Djerdj, Z

    I. Djerdj, Z. Jagli ˇci´c, D. Ar ˇcon, and M. Niederberger, Nanoscale 2, 1096 (2010)

    work page 2010

  3. [4]

    Janisch, P

    R. Janisch, P . Gopal, and N. A. Spaldin, Journal of Physic s: Condensed Matter 17, R657 (2005). 8

    work page 2005

  4. [5]

    H. A. Weakliem, The Journal of Chemical Physics 36, 2117 (1962)

    work page 1962

  5. [6]

    K. C. V erma and R. K. Kotnala, Phys. Chem. Chem. Phys. 18, 17565 (2016)

    work page 2016

  6. [7]

    W. Li, G. Wang, C. Chen, J. Liao, and Z. Li, Nanomaterials 7, 20 (2017)

    work page 2017

  7. [8]

    I. A. Sarsari, C. D. Pemmaraju, H. Salamati, and S. Sanvit o, Phys. Rev. B 87, 245118 (2013)

    work page 2013

  8. [9]

    A. L. Schoenhalz and G. M. Dalpian, Phys. Chem. Chem. Phys 15, 15863 (2013)

    work page 2013

Show all 38 references
  1. [10]

    C. H. Patterson, Phys. Rev. B 74, 144432 (2006)

    2006

  2. [11]

    K. Yim, J. Lee, D. Lee, M. Lee, E. Cho, H. S. Lee, H.-H. Nahm , and S. Han, Sci. Rep. 7, 40907 (2017)

    2017

  3. [12]

    Das Sarma, American Scientist 89, 516 (2001)

    S. Das Sarma, American Scientist 89, 516 (2001)

    2001

  4. [13]

    Feynman, Foundation of Physics 16, 507 (1986)

    R. Feynman, Foundation of Physics 16, 507 (1986)

    1986

  5. [14]

    Ronning, C

    C. Ronning, C. Borschel, S. Geburt, R. Niepelt, S. Mülle r, D. Stichtenoth, J. P . Richters, A. Dev, T. V oss, L. Chen, et al., physica status solidi (b) 247, 2329 (2010), ISSN 1521-3951, URL http://dx.doi.org/10.1002/pssb.201046192

    2010 doi

  6. [15]

    Geburt, M

    S. Geburt, M. Lorke, A. L. da Rosa, T. Frauenheim, R. Röde r, T. V oss, U. Kaiser, W. Heimbrodt, and C. Ronning, Nano Letters 14, 4523 (2014)

    2014

  7. [16]

    Geburt, R

    S. Geburt, R. Röder, U. Kaiser, L. Chen, M.-H. Chu, J. Seg ura-Ruiz, G. Martínez-Criado, W. Heim- brodt, and C. Ronning, physica status solidi (RRL) – Rapid Re search Letters 7, 886 (2013)

    2013

  8. [17]

    Segura-Ruiz, G

    J. Segura-Ruiz, G. Martínez-Criado, M. H. Chu, S. Gebur t, and C. Ronning, Nano Letters 11, 5322 (2011)

    2011

  9. [18]

    Z. Jin, T. Fukumura, M. Kawasaki, K. Ando, H. Saito, T. Se kiguchi, Y . Z. Y oo, M. Murakami, Y . Mat- sumoto, T. Hasegawa, et al., Applied Physics Letters 78, 3824 (2001)

    2001

  10. [19]

    C. A. Johnson, T. C. Kaspar, S. A. Chambers, G. M. Salley, and D. R. Gamelin, Phys. Rev. B 81, 125206 (2010)

    2010

  11. [20]

    Hohenberg and W

    P . Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964)

    1964

  12. [21]

    Kohn and L

    W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965)

    1965

  13. [22]

    Hedin, Phys

    L. Hedin, Phys. Rev. 139, A796 (1965)

    1965

  14. [23]

    Kresse and D

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

    1999

  15. [24]

    P . E. Blöchl, Phys. Rev. B 50, 17953 (1994)

    1994

  16. [25]

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

    1996

  17. [26]

    Lorke, T

    M. Lorke, T. Frauenheim, and A. L. da Rosa, Phys. Rev. B 93, 115132 (2016)

    2016

  18. [27]

    Franke, A

    D. Franke, A. L. da Rosa, M. Lorke, and T. Frauenheim, phy s. stat. sol. (b) 256, 1800455 (2019). 9

    2019

  19. [29]

    D. Wang, G. Xing, M. Gao, L. Y ang, J. Y ang, and T. Wu, The Jo urnal of Phys- ical Chemistry C 115, 22729 (2011), http://pubs.acs.org/doi/pdf/10.1021/jp 204572v, URL http://pubs.acs.org/doi/abs/10.1021/jp204572v

    2011 doi

  20. [31]

    Janotti and C

    A. Janotti and C. G. van de Walle, Phys. Rev. B 76, 165202 (2007)

    2007

  21. [32]

    S. J. Clark, J. Robertson, S. Lany, and A. Zunger, Phys. R ev. B 81, 115311 (2010), URL http://link.aps.org/doi/10.1103/PhysRevB.81.115311

    2010 doi

  22. [33]

    Lany and A

    S. Lany and A. Zunger, Phys. Rev. B 81, 113201 (2010)

    2010

  23. [34]

    C. G. V an de Walle and J. Neugebauer, Journal of Applied P hysics 95, 3851 (2004)

    2004

  24. [35]

    F. Oba, A. Togo, I. Tanaka, J. Paier, and G. Kresse, Phys. Rev. B 77, 245202 (2008)

    2008

  25. [36]

    S. B. Zhang, S.-H. Wei, and A. Zunger, Phys. Rev. B 63, 075205 (2001)

    2001

  26. [37]

    W. M. Haynes, ed., CRC handbook of chemistry and physics: A ready-reference bo ok of chemical and physical data (CRC Press, Boca Raton, 2016), 97th ed

    2016

  27. [38]

    Kaxiras, Atomic and Electronic Structure of Solids(Cambridge University Press, Cambridge, 2003)

    E. Kaxiras, Atomic and Electronic Structure of Solids(Cambridge University Press, Cambridge, 2003)

    2003

  28. [39]

    Janotti and C

    A. Janotti and C. G. V an de Walle, Phys. Rev. B 76, 165202 (2007)

    2007

  29. [40]

    Huang, C.-Y

    G.-Y . Huang, C.-Y . Wang, and J.-T. Wang, J. Phys.: Condens. Matter. 21, 195403 (2009)

    2009

  30. [41]

    G. S. Chang, E. Z. Kurmaev, D. W. Boukhvalov, L. D. Finkel stein, S. Colis, T. M. Pedersen, A. Moewes, and A. Dinia, Phys. Rev. B 75, 195215 (2007). 10 -8 -6 -4 -2 0 2 4 6 8 Energy (eV) -8 -6 -4 -2 0 2 4 6 8DOS (states/eV) total Zn-s O-px O-py O-pz Co-px Co-py Co-pz Co-dxy Co-...

    2007

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.