pith. sign in

arxiv: 2604.22393 · v2 · pith:3VBPNAK7new · submitted 2026-04-24 · ❄️ cond-mat.mtrl-sci · physics.comp-ph

Nature of point defects in bulk hexagonal diamond

Ling Zhu , Xuanxuan Zhang , Guliqinayi Alimu , Chen-Min Dai , Chunlan Ma , Zenghua Cai This is my paper

Pith reviewed 2026-05-19 18:05 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.comp-ph
keywords hexagonal diamondpoint defectscarbon vacancydopantsconductivitycolor centersspin statesfirst-principles calculations
0
0 comments X

The pith

First-principles calculations show carbon vacancies dominate intrinsic conductivity in hexagonal diamond while boron and nitrogen dopants enable p-type and n-type behavior.

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

The paper uses first-principles methods to map the formation energies, charge states, and electronic levels of intrinsic defects and extrinsic dopants in bulk hexagonal diamond. It establishes that the carbon vacancy is the dominant intrinsic defect controlling conductivity, while the carbon interstitial is unstable. Boron is identified as an effective acceptor that promotes p-type conductivity, and nitrogen or phosphorus act as donors for n-type conductivity. Defect complexes involving vacancies with magnesium or group-V elements introduce multiple spin and charge states inside the bandgap. These findings clarify how defects can be used to engineer conductivity and host quantum states in this recently synthesized carbon allotrope.

Core claim

In hexagonal diamond the carbon vacancy has the lowest formation energy among intrinsic point defects and introduces gap states that dominate intrinsic conductivity, whereas the carbon interstitial is unstable. Boron incorporates readily as a shallow acceptor that enhances p-type conductivity, while nitrogen and phosphorus function as effective donors for n-type conductivity; group-II and group-IV dopants exhibit high formation energies or remain neutral. Complexes such as VC, MgC, and XV display multiple charge and spin states within the hexagonal-diamond bandgap, positioning them as candidate color centers for qubit applications.

What carries the argument

First-principles calculations of defect formation energies and charge-transition levels for vacancies, interstitials, substitutional dopants, and their complexes in the hexagonal diamond lattice.

If this is right

  • Boron doping can be used to achieve controllable p-type conductivity in hexagonal diamond devices.
  • Nitrogen or phosphorus incorporation can produce n-type conductivity for electronic applications.
  • Vacancy-containing complexes such as VC, MgC, and XV can be engineered as color centers with multiple spin states for quantum information processing.
  • Group-II and group-IV dopants will have limited utility because of their high formation energies or neutral character.

Where Pith is reading between the lines

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

  • If the predicted defect levels are confirmed, hexagonal diamond could be integrated into hybrid quantum-classical devices that exploit both its mechanical hardness and tunable defect states.
  • Systematic comparison of the same defect calculations between hexagonal and cubic diamond would reveal whether the hexagonal structure offers distinct advantages for defect-based quantum sensing.
  • Controlled introduction of the identified complexes during synthesis could be tested as a route to stable single-photon sources or spin qubits in bulk material.

Load-bearing premise

The first-principles method used accurately computes formation energies, charge states, and stability of defects without large errors from functional choice or supercell size.

What would settle it

Experimental measurement of the dominant carrier type and activation energy in high-purity, undoped bulk hexagonal diamond samples that fails to match the predicted dominance of carbon-vacancy-related states.

read the original abstract

Hexagonal diamond (HD), an exotic carbon allotrope recently synthesized in bulk form, exhibits superior mechanical properties compared to cubic diamond (CD) and holds promise for advanced industrial and quantum applications. Using first-principles calcu-lations, we systematically investigate intrinsic defects, extrinsic dopants, and defect complexes in HD. Our study shows that VC dominates intrinsic conductivity, while Ci is unstable. Among extrinsic dopants, boron acts as a benign acceptor enhancing p-type conductivity, whereas nitrogen and phosphorus serve as effective donors for n-type conductivity. Group II and Group IV dopants, however, introduce high formation energies or neutral charge states with limited impact. Furthermore, VC, MgC and XV defect com-plexes display multiple spin and charge states within the HD band gap, highlighting their potential as color centers for hosting qubits. These results not only clarify the defect physics of HD but also demonstrate its broader implications for conductivity engineering and quantum technologies.

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 / 3 minor

Summary. The manuscript reports first-principles calculations on intrinsic defects, extrinsic dopants, and defect complexes in bulk hexagonal diamond (HD). Key findings include the dominance of carbon vacancies (VC) in intrinsic conductivity, instability of carbon interstitials (Ci), boron as an effective acceptor for p-type doping, nitrogen and phosphorus as donors for n-type, and the potential of VC, MgC, and XV complexes as multi-state color centers for quantum applications.

Significance. This work is significant for advancing understanding of defect engineering in HD, a recently bulk-synthesized carbon allotrope with superior mechanical properties. The results on conductivity control and qubit candidates could impact materials design for industrial and quantum technologies. Strengths include systematic investigation of various defects; however, the impact is limited by potential methodological limitations in the DFT approach.

major comments (3)
  1. [Methods] Methods section: The exchange-correlation functional is not explicitly named or justified, and no convergence tests for supercell size, k-point sampling, or finite-size corrections for charged defects are reported. These details are load-bearing for the formation-energy diagrams and the identification of stable charge states for VC and the dopants.
  2. [Results on intrinsic defects] Results on intrinsic defects (formation-energy plots): The diagrams for VC and Ci do not mention image-charge corrections or potential alignment. Standard GGA functionals underestimate the HD gap by ~1.5 eV relative to experiment, which can shift thermodynamic transition levels and alter the conclusion that VC dominates intrinsic conductivity while Ci is unstable.
  3. [Dopant analysis] Dopant analysis section: Charge-transition levels for B, N, and P are presented without hybrid-functional benchmarks or comparison to cubic-diamond literature. This directly affects the claims that B is a benign acceptor and N/P are effective donors, as level positions within the gap are sensitive to the ~1 eV gap error.
minor comments (3)
  1. [Abstract] Abstract contains a line-break hyphenation error ('calcu-lations').
  2. [Figures] Figure captions should specify supercell sizes and whether spin polarization was included for the multi-state complexes.
  3. [Discussion] A brief comparison paragraph to known defects in cubic diamond would strengthen the discussion of HD-specific behavior.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment below and will make revisions to enhance the clarity and robustness of the presented results. Our core conclusions regarding the dominance of vacancies in conductivity and the roles of dopants remain supported by the calculations.

read point-by-point responses

  1. Referee: Methods section: The exchange-correlation functional is not explicitly named or justified, and no convergence tests for supercell size, k-point sampling, or finite-size corrections for charged defects are reported. These details are load-bearing for the formation-energy diagrams and the identification of stable charge states for VC and the dopants.

    Authors: We thank the referee for this important point. The manuscript employs the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) functional, which we will now explicitly name and justify in the revised Methods section by noting its widespread use in defect calculations for diamond structures. We performed convergence tests with respect to supercell size (up to 216-atom cells), k-point sampling (using Monkhorst-Pack grids), and energy cutoff, but these were not detailed in the original submission. We will add this information. Finite-size corrections were applied using the image charge correction method, and we will report the details of potential alignment in the revised version. revision: yes

  2. Referee: Results on intrinsic defects (formation-energy plots): The diagrams for VC and Ci do not mention image-charge corrections or potential alignment. Standard GGA functionals underestimate the HD gap by ~1.5 eV relative to experiment, which can shift thermodynamic transition levels and alter the conclusion that VC dominates intrinsic conductivity while Ci is unstable.

    Authors: We acknowledge the omission of explicit mention of corrections in the text and figures. Image charge corrections and potential alignment were indeed performed using established methods for charged defects in periodic boundary conditions. We will update the manuscript to include this information. On the band gap issue, although GGA underestimates the gap, the formation energies and relative stabilities are calculated consistently, and the conclusion that VC is dominant is based on its much lower formation energy compared to other defects like Ci. We will add a discussion noting this limitation and that absolute transition levels may require hybrid functionals for precise positioning, but the qualitative findings hold. revision: partial

  3. Referee: Dopant analysis section: Charge-transition levels for B, N, and P are presented without hybrid-functional benchmarks or comparison to cubic-diamond literature. This directly affects the claims that B is a benign acceptor and N/P are effective donors, as level positions within the gap are sensitive to the ~1 eV gap error.

    Authors: We agree that additional benchmarks would strengthen the claims. We will add comparisons to the cubic diamond literature, where similar GGA studies show B, N, and P behaving as acceptors and donors respectively. Furthermore, we conducted additional calculations with the HSE06 hybrid functional on selected dopant configurations in smaller supercells, which support the p-type and n-type doping effectiveness. These will be included in the revision to address concerns about the gap error. The identification of B as benign and N/P as effective is based on their low formation energies and suitable charge transition levels within our computational gap. revision: yes

Circularity Check

0 steps flagged

No circularity in first-principles defect calculations for hexagonal diamond

full rationale

The paper derives defect formation energies, charge states, and stability conclusions directly from standard DFT-based first-principles calculations applied to the hexagonal diamond structure. These computations follow established quantum mechanical protocols without reducing to fitted parameters from prior HD data, self-definitional loops, or load-bearing self-citations that presuppose the target results. Claims about VC dominance, Ci instability, and dopant roles emerge from the computed energies and transition levels rather than being constructed into the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The study relies on standard DFT assumptions for defect formation energies and electronic structure. No free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Standard DFT approximations (e.g., PBE or hybrid functional) sufficiently describe defect energetics in carbon allotropes.
    Invoked implicitly for all formation-energy and charge-state results.

pith-pipeline@v0.9.0 · 5705 in / 1113 out tokens · 33994 ms · 2026-05-19T18:05:18.326094+00:00 · methodology

Review history (2 revisions) →

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

29 extracted references · 29 canonical work pages

  1. [1]

    S. Lai, X. Yang, J. Shi, S. Liu, Y . Guo, L. Yan, J. Zang, Z. Zhang, Q. Jia, J. Sun, S. Cheng, and C. Shan, Bulk hexagonal diamond, Nature 651, 621-625 (2026)

    work page 2026

  2. [2]

    L. Yang, K. C. Lau, Z. Zeng, D. Zhang, H. Tang, B. Yan, G. Niu, H. Gou, Y . Yang, W. Yang, D. Luo, and H. K. Mao, Synthesis of bulk hexagonal diamond, Nature 644, 370-375 (2025)

    work page 2025

  3. [3]

    Qingkun, S

    L. Qingkun, S. Yi, L. Zhiyuan, and Z. Yu, Lonsdaleite – A material stronger and stiffer than diamond, Scr. Mater. 65, 229-232 (2011)

    work page 2011

  4. [4]

    K. A. Abdelghafar, D. S. Choi, and K. A. Askar, Luminescence lineshapes of nitrogen vacancy center in lonsdaleite and dual structure of diamond/lonsdaleite: a DFT study, Sci. Rep. 15, 15334 (2025)

    work page 2025

  5. [5]

    R. J. Narayan, R. D. Boehm, and A. V . Sumant, Medical applications of diamond particles & surfaces, Mater. Today 14, 154-163 (2011)

    work page 2011

  6. [6]

    Chatterjee, P

    A. Chatterjee, P. Stevenson, S. De Franceschi, A. Morello, N. P. de Leon, and F. Kuemmeth, Semiconductor qubits in practice, Nat. Rev. Phys. 3, 157-177 (2021)

    work page 2021

  7. [7]

    Gali, Ab initio theory of the nitrogen‐vacancy center in diamond, Nanophotonics 8, 1907- 1943 (2019)

    Á. Gali, Ab initio theory of the nitrogen‐vacancy center in diamond, Nanophotonics 8, 1907- 1943 (2019)

    work page 1907

  8. [8]

    Manian, M

    A. Manian, M. O. d. Vries, D. Stavrevski, Q. Sun, S. P. Russo, and A. D. Greentree, Nitrogen-vacancy centres in lonsdaleite: a novel nanoscale sensor?, Phys. Chem. Chem. Phys. 27, 21837-21849 (2025)

    work page 2025

  9. [9]

    Z. Pan, H. Sun, Y . Zhang, and C. Chen, Harder than diamond: superior indentation strength of wurtzite BN and lonsdaleite, Phys. Rev. Lett. 102, 055503 (2009)

    work page 2009

  10. [10]

    Freysoldt, B

    C. Freysoldt, B. Grabowski, T. Hickel, J. Neugebauer, G. Kresse, A. Janotti, and C. G. Van de Walle, First-principles calculations for point defects in solids, Rev. Mod. Phys. 86, 253-305 (2014)

    work page 2014

  11. [11]

    C. Qiu, H. P. Liang, S. Geng, Y . Chen, X. L. Yan, H. X. Deng, and S. H. Wei, Chemical Control of Zero -Phonon Line Wavelength in Diamond Color Centers for Telecom -Band Emission, J. Am. Chem. Soc. 148, 13454-13462 (2026)

    work page 2026

  12. [12]

    Muruganathan, and H

    M. Muruganathan, and H. Mizuta, Boron vacancy color center in diamond: Ab initio study, Diam. Relat. Mater. 114, 108341 (2021)

    work page 2021

  13. [13]

    Thiering, and A

    G. Thiering, and A. Gali, Ab Initio Magneto-Optical Spectrum of Group-IV Vacancy Color Centers in Diamond, Phys. Rev. X 8, (2018)

    work page 2018

  14. [14]

    C. E. Dreyer, A. Alkauskas, J. L. Lyons, A. Janotti, and C. G. Van de Walle, First-Principles Calculations of Point Defects for Quantum Technologies, Annu. Rev. Mater. Res. 48, 1 -26 (2018)

    work page 2018

  15. [15]

    K. A. Abdelghafar, D. S. Choi, and K. Askar, Ab initio study of N2V color center in lonsdaleite/diamond dual structure: A two -site Hubbard model, Comput. Mater. Sci. 263, (2026)

    work page 2026

  16. [16]

    K. A. Abdelghafar, D. S. Choi, and K. Askar, A density functional theory study of nitrogen vacancy center in lonsdaleite, J. Appl. Phys. 58, (2024)

    work page 2024

  17. [17]

    Kresse, and J

    G. Kresse, and J. Furthmüller, Efficiency of ab -initio total energy calculations for metals and semiconductors using a plane-wave basis set, Comput. Mater. Sci. 6, 15-50 (1996)

    work page 1996

  18. [18]

    H. J. Monkhorst, and J. D. Pack, Special points for Brillouin-zone integrations, Phys. Rev. B 13, 5188-5192 (1976)

    work page 1976

  19. [19]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 77, 3865 (1996)

    work page 1996

  20. [20]

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid functionals based on a screened Coulomb potential, J. Chem. Phys. 118, 8207-8215 (2003)

    work page 2003

  21. [21]

    Gao, Band gaps and dielectric functions of cubic and hexagonal diamond polytypes calculated by many-body perturbation theory, Phys

    S.-P. Gao, Band gaps and dielectric functions of cubic and hexagonal diamond polytypes calculated by many-body perturbation theory, Phys. Status Solidi B 252, 235-242 (2015)

    work page 2015

  22. [22]

    See Supplemental Material at for the structural configurations and defect states of CaC, SrC, Ni1, and for the Benchmark calculation, formation energy calculations of all the key competing secondary phases and all the allowed chemical potential results of C and extrinsic impurity elements.,

    work page

  23. [23]

    R. Bai, D. Liu, S. Chen, and Y .-N. Wu, Origin of the unidentified color centers in 4H-SiC from first principles, Phys. Rev. B 111, 014106 (2025)

    work page 2025

  24. [24]

    Huang, Z

    M. Huang, Z. Zheng, Z. Dai, X. Guo, S. Wang, L. Jiang, J. Wei, and S. Chen, DASP: Defect and Dopant ab-initio Simulation Package, J. Semicond. 43, 042101 (2022)

    work page 2022

  25. [25]

    Freysoldt, J

    C. Freysoldt, J. Neugebauer, and C. G. Van de Walle, Fully Ab Initio Finite -Size Corrections for Charged-Defect Supercell Calculations, Phys. Rev. Lett. 102, 016402 (2009)

    work page 2009

  26. [26]

    De, and C

    A. De, and C. E. Pryor, Electronic structure and optical properties of Si, Ge and diamond in the lonsdaleite phase, J. Phys. Condens. Matter 26, 045801 (2014)

    work page 2014

  27. [27]

    C. A. Coulson, and F. P. Larkins, Electronic structure of the neutral isolated divacancy in diamond, J. Phys. Chem. Solids 30, 1963-1972 (1969)

    work page 1963

  28. [28]

    Amato, T

    M. Amato, T. Kaewmaraya, and A. Zobelli, Extrinsic Doping in Group IV Hexagonal - Diamond-Type Crystals, J. Phys. Chem. C 124, 17290-17298 (2020)

    work page 2020

  29. [29]

    Koizumi, and M

    S. Koizumi, and M. Suzuki, n -Type doping of diamond, Phys. Status Solidi A 203, 3358- 3366 (2006)

    work page 2006