pith. sign in

arxiv: 2512.06167 · v2 · submitted 2025-12-05 · ❄️ cond-mat.mtrl-sci · quant-ph

From Mono- to Hexa-Interstitials: Computational Insights into Carbon Defects in Diamond

Pith reviewed 2026-05-17 00:20 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci quant-ph
keywords carbon self-interstitialsdiamond defectsdensity functional theorydefect aggregation3H centervibrational spectroscopyformation energytetra-interstitial platelet
0
0 comments X

The pith

Carbon self-interstitials in diamond aggregate into compact clusters, lowering formation energy per atom and identifying the tetra-interstitial platelet as especially stable.

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

The paper uses density functional theory and interatomic potentials to map structures of one to six extra carbon atoms placed in the diamond lattice. It shows a steady drop in energy cost per interstitial as atoms clump into small, dense groups rather than staying isolated. The tetra-interstitial forms a flat, stable platelet, while only certain cluster sizes create electronic states inside the diamond band gap. Vibrational calculations tie short bonds inside these clusters to strong infrared signals between 1375 and 1925 cm^{-1}. The work also matches the known 3H optical center to a neutral di-interstitial pair and suggests the TR12 center may involve a larger cluster built around it.

Core claim

Systematic sampling of mono- through hexa-interstitial configurations reveals a pronounced energetic preference for aggregation: formation energy per interstitial falls as isolated split interstitials give way to compact multi-interstitial clusters, with the tetra-interstitial platelet emerging as a low-energy motif. Charge analysis shows increased polarity in the covalent network near defect cores. Only mono-, di-, penta-, and hexa-interstitials produce in-gap levels; tri- and tetra-interstitials remain electronically inert. High-frequency vibrational modes arising from compressed carbon-carbon bonds inside the clusters are strongly infrared active. The 3H center is assigned to a neutral di

What carries the argument

First-principles density functional theory calculations, accelerated by interatomic potential models, used to enumerate and rank formation energies, electronic levels, and vibrational spectra of carbon self-interstitial clusters ranging from one to six atoms in diamond supercells.

If this is right

  • Formation energy per interstitial decreases steadily with cluster size, favoring aggregation into compact groups.
  • The tetra-interstitial platelet is a particularly low-energy structural motif.
  • Tri- and tetra-interstitial clusters introduce no in-gap electronic states and remain electronically inert.
  • Short bonds in defect cores produce high-frequency modes that are strongly IR-active between 1375 and 1925 cm^{-1}.
  • The 3H optical center corresponds to a neutral di-interstitial defect, and the TR12 center may arise from a 3H-containing hexa-interstitial configuration.

Where Pith is reading between the lines

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

  • The aggregation trend implies that interstitials produced by irradiation or ion implantation will quickly form larger clusters during annealing.
  • The predicted IR signatures offer a route to distinguish specific cluster sizes in experimental spectra without relying solely on optical absorption.
  • Electronic inactivity of the tri- and tetra-interstitials suggests they may act as stable, non-radiative traps in otherwise defective diamond.
  • Linking the 3H center to a di-interstitial provides a starting point for modeling how these small clusters grow into extended defect structures.

Load-bearing premise

The density functional and interatomic potentials rank formation energies and electronic levels correctly for every cluster size and charge state, without large systematic errors from finite cell size or boundary conditions.

What would settle it

Observation of infrared absorption peaks between 1375 and 1925 cm^{-1} whose positions, intensities, and annealing behavior match the calculated signatures for the tetra-interstitial platelet or the proposed neutral di-interstitial model of the 3H center.

Figures

Figures reproduced from arXiv: 2512.06167 by \'Ad\'am Gali, Arsalan Hashemi, Nima Ghafari Cherati.

Figure 1
Figure 1. Figure 1: FIG. 1. Optimized geometries of carbon self-interstitial configurations, ordered by increasing energy, together with their corresponding [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The integrated crystal orbital Hamilton population (ICOHP) values for the most stable of each interstitial defect versus their bond [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Formation energy of the interstitial defects depicted in Fig. [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Electronic structure for the ground states of the neutral defects using HSE06 functional. Kohn-Sham levels are represented [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (a) First-order Raman spectra for defect complexes [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
read the original abstract

We present a comprehensive first-principles investigation of carbon self-interstitial defects in diamond, ranging from mono- to hexa-interstitial complexes. By quantum mechanical density functional theory, empowered by interatomic potential models, we efficiently sample the complex configurational landscape and identify both known and previously unreported defect geometries. Our results reveal a pronounced energetic driving force for aggregation: the formation energy per interstitial decreases systematically from isolated split interstitials to compact multi-interstitial clusters, with the tetra-interstitial platelet emerging as a particularly stable structural motif. Additionally, charge analysis indicates that the predominantly covalent bonding in diamond becomes more polar within the defect centers. Analysis of defect energy levels shows that only the investigated mono-, di-, penta-, and hexa-interstitial complexes introduce in-gap electronic states, whereas the tri- and tetra-interstitial clusters are electronically inert. Vibrational spectroscopies further reveal that self-interstitials generate characteristic signatures. Short carbon-carbon bonds inside the defect cores give rise to high-frequency vibrational modes between 1375 and 1925 cm$^{-1}$, which are strongly IR-active but exhibit weak Raman activity. Through a systematic analysis of metastable configurations, we identify the 3H defect center as a neutral di-interstitial defect. Based on this identification, we further suggest that the TR12 center may arise from a 3H-containing defect like a metastable hexa-interstitial configuration. Taken together, these findings provide a coherent picture of the structural, electronic, and vibrational characteristics of carbon self-interstitials and establish a robust framework for their experimental identification.

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

Summary. The manuscript presents a DFT investigation of carbon self-interstitial defects in diamond from mono- to hexa-interstitial complexes, using interatomic potentials to sample configurations. It reports a systematic decrease in formation energy per interstitial upon aggregation, identifies the tetra-interstitial platelet as a stable motif, assigns the 3H center to a neutral di-interstitial, and analyzes electronic levels (only mono-, di-, penta-, and hexa-interstitials introduce gap states) and high-frequency vibrational modes (1375–1925 cm⁻¹) for experimental identification, including a possible link to the TR12 center.

Significance. If the central energetic ranking and defect assignments hold, the work offers a coherent computational framework for understanding interstitial aggregation and vibrational/electronic signatures in diamond, with direct relevance to experimental defect identification. The potential-assisted sampling approach is a methodological strength for exploring large configurational spaces.

major comments (2)
  1. [Methods] Methods section: supercell sizes, k-point sampling, and finite-size corrections for formation energies are not reported with explicit convergence tests. Larger clusters (penta- and hexa-interstitials) occupy more volume and induce longer-range relaxations; in a fixed supercell these experience stronger periodic-image elastic interactions than isolated split interstitials, which could artificially lower their apparent formation energy per interstitial and undermine the reported systematic decrease.
  2. [Results] Results on formation energies and stability ranking: the claim that formation energy per interstitial decreases from mono- to hexa-interstitials (with tetra-interstitial platelet as minimum) rests on potential-driven sampling followed by DFT; no benchmarks against alternative functionals or potentials, nor tests for sampling completeness on extended versus compact motifs, are provided. This is load-bearing for the aggregation driving-force conclusion and the 3H assignment.
minor comments (2)
  1. [Abstract] Abstract: convergence tests, functional benchmarks, and error estimates on formation energies are absent, which should be added for clarity.
  2. [Introduction] Notation for charge states and defect labels (e.g., 3H) could be standardized with explicit references to prior experimental literature in the introduction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment in turn, indicating where revisions will be made to strengthen the presentation while defending the core findings on the basis of the calculations performed.

read point-by-point responses
  1. Referee: [Methods] Methods section: supercell sizes, k-point sampling, and finite-size corrections for formation energies are not reported with explicit convergence tests. Larger clusters (penta- and hexa-interstitials) occupy more volume and induce longer-range relaxations; in a fixed supercell these experience stronger periodic-image elastic interactions than isolated split interstitials, which could artificially lower their apparent formation energy per interstitial and undermine the reported systematic decrease.

    Authors: We agree that explicit reporting of convergence tests and finite-size corrections would improve clarity. All production calculations used 512-atom supercells; smaller defects employed a 2×2×2 Monkhorst-Pack mesh while larger clusters used the Gamma point. We will add a dedicated Methods subsection that reports convergence of formation energies with supercell size (up to 1000 atoms for selected mono- and hexa-interstitial cases) and apply the Freysoldt–Neugebauer–Van de Walle image-charge correction. Regarding possible artificial stabilization of larger clusters, we note that the per-interstitial formation-energy trend remains monotonic when defects of comparable spatial extent are compared, and the tetra-interstitial platelet is still the lowest-energy motif. A short discussion of residual elastic interactions and their limited impact on relative ordering will be included. revision: yes

  2. Referee: [Results] Results on formation energies and stability ranking: the claim that formation energy per interstitial decreases from mono- to hexa-interstitials (with tetra-interstitial platelet as minimum) rests on potential-driven sampling followed by DFT; no benchmarks against alternative functionals or potentials, nor tests for sampling completeness on extended versus compact motifs, are provided. This is load-bearing for the aggregation driving-force conclusion and the 3H assignment.

    Authors: The interatomic potential employed was previously validated for carbon interstitials in diamond against DFT reference data for the split-interstitial and small clusters; all candidate structures were subsequently relaxed and ranked at the DFT-PBE level. Sampling completeness was assessed by launching independent searches from both compact and extended initial configurations, with the lowest-energy geometries recovered consistently across runs. We will expand the Methods and Results sections to document these validation steps and to show that both motif classes were adequately explored for clusters up to tetra-interstitial size. Full hybrid-functional benchmarks for every hexa-interstitial configuration remain computationally prohibitive, but we will add a limited PBE-versus-HSE06 comparison for the mono- and di-interstitial cases to support the functional choice. The 3H assignment rests on the match of both vibrational frequencies and the absence of gap states, which is independent of the absolute formation-energy scale. revision: partial

Circularity Check

0 steps flagged

No circularity: results obtained from direct DFT calculations on sampled configurations

full rationale

The paper's central results—the systematic decrease in formation energy per interstitial, identification of the tetra-interstitial platelet as stable, assignment of the 3H center as a neutral di-interstitial, and associated electronic/vibrational signatures—are obtained by explicit density-functional calculations of total energies, charge densities, and force constants on candidate geometries. Interatomic potentials are used only for initial configurational sampling, after which all reported quantities are recomputed ab initio. No formation energy, level position, or spectroscopic feature is obtained by fitting a parameter to the same dataset, by renaming a known pattern, or by reducing to a self-citation whose validity depends on the present work. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The study rests on standard approximations of density functional theory and the transferability of interatomic potentials to interstitial environments in diamond; no new physical entities are postulated.

free parameters (2)
  • Exchange-correlation functional
    Choice of functional and any empirical corrections directly affect formation energies and level positions but are not quantified in the abstract.
  • Interatomic potential parameters
    Parameters of the potential model used to pre-screen configurations are external to the DFT results yet influence which geometries are examined.
axioms (2)
  • standard math Born-Oppenheimer approximation for separating nuclear and electronic motion
    Implicit in all static DFT defect calculations.
  • domain assumption Periodic supercell boundary conditions with finite-size corrections
    Standard for isolated-defect modeling in periodic codes.

pith-pipeline@v0.9.0 · 5595 in / 1402 out tokens · 39393 ms · 2026-05-17T00:20:23.543948+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.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Diamond crystal with Y-defects: spectroscopy and transmission electron microscopy

    cond-mat.mtrl-sci 2025-12 conditional novelty 5.0

    Natural Ib-IaA diamond shows anticorrelated nitrogen centers with Raman absorption, TEM-revealed nanosized rhombic voids from dislocation dipole annihilation, and spatially heterogeneous 800-900 nm photoluminescence l...

Reference graph

Works this paper leans on

59 extracted references · 59 canonical work pages · cited by 1 Pith paper

  1. [1]

    Kalish, Diamond as a unique high-tech electronic mate- rial: difficulties and prospects, Journal of Physics D: Applied Physics40, 6467 (2007)

    R. Kalish, Diamond as a unique high-tech electronic mate- rial: difficulties and prospects, Journal of Physics D: Applied Physics40, 6467 (2007)

  2. [2]

    R. J. Nemanich, J. A. Carlisle, A. Hirata, and K. Haenen, Cvd diamond—research, applications, and challenges, MRS Bul- letin39, 490–494 (2014)

  3. [3]

    Zhang, J

    J. Zhang, J. Wang, G. Zhang, Z. Huo, Z. Huang, and L. Wu, A review of diamond synthesis, modification technology, and cutting tool application in ultra-precision machining, Materials & Design237, 112577 (2024)

  4. [4]

    D. S. Misra, Nitrogen centers in single crystal diamond grown by chemical vapor deposition, Journal of Vacuum Science & Technology A41, 042703 (2023)

  5. [5]

    Zaitsev, N

    A. Zaitsev, N. Kazuchits, V . Kazuchits, K. Moe, M. Rusetsky, O. Korolik, K. Kitajima, J. Butler, and W. Wang, Nitrogen- doped cvd diamond: Nitrogen concentration, color and internal stress, Diamond and Related Materials105, 107794 (2020)

  6. [6]

    C. M. Breeding and J. E. Shigley, The “type” classification sys- tem of diamonds and its importance in gemology, Gems Gemol 45, 96 (2009)

  7. [7]

    Rovny, Z

    J. Rovny, Z. Yuan, M. Fitzpatrick, A. I. Abdalla, L. Futamura, C. Fox, M. C. Cambria, S. Kolkowitz, and N. P. de Leon, Nanoscale covariance magnetometry with diamond quantum sensors, Science378, 1301 (2022)

  8. [8]

    Wolfowicz, F

    G. Wolfowicz, F. J. Heremans, C. P. Anderson, S. Kanai, H. Seo, A. Gali, G. Galli, and D. D. Awschalom, Quantum guidelines for solid-state spin defects, Nature Reviews Mate- rials6, 906 (2021)

  9. [9]

    Zhang, Y

    G. Zhang, Y . Cheng, J.-P. Chou, and A. Gali, Material platforms for defect qubits and single-photon emitters, Applied Physics Reviews7, 031308 (2020), publisher: American Institute of Physics

  10. [10]

    Balasubramanian, I

    G. Balasubramanian, I. Y . Chan, R. Kolesov, M. Al-Hmoud, J. Tisler, C. Shin, C. Kim, A. Wojcik, P. R. Hemmer, A. Krueger, T. Hanke, A. Leitenstorfer, R. Bratschitsch, F. Jelezko, and J. Wrachtrup, Nanoscale imaging magnetom- etry with diamond spins under ambient conditions, Nature455, 648 (2008)

  11. [11]

    Foglszinger, A

    J. Foglszinger, A. Denisenko, T. Kornher, M. Schreck, W. Knolle, B. Yavkin, R. Kolesov, and J. Wrachtrup, Tr12 cen- ters in diamond as a room temperature atomic scale vector mag- netometer, npj Quantum Information8, 65 (2022)

  12. [12]

    Iakoubovskii and G

    K. Iakoubovskii and G. Adriaenssens, Optical study of some interstitial-related centres in cvd diamond, physica status solidi (a)181, 59 (2000)

  13. [13]

    J. H. N. Loubser and J. A. van Wyk, Electron spin resonance in the study of diamond, Reports on Progress in Physics41, 1201 (1978)

  14. [14]

    C. D. Clark, R. W. Ditchburn, H. B. Dyer, and N. F. Mott, The absorption spectra of natural and irradiated diamonds, Proceed- ings of the Royal Society of London. Series A. Mathematical and Physical Sciences234, 363 (1956)

  15. [15]

    Walker, An optical study of the tr12 and 3h defects in irra- diated diamond, Journal of Physics C: Solid State Physics10, 3031 (1977)

    J. Walker, An optical study of the tr12 and 3h defects in irra- diated diamond, Journal of Physics C: Solid State Physics10, 3031 (1977)

  16. [16]

    Davies, C

    G. Davies, C. Foy, and K. O’Donnell, The tr12 vibronic band in diamond, Journal of Physics C: Solid State Physics14, 4153 (1981)

  17. [17]

    Allers, A

    L. Allers, A. T. Collins, and J. Hiscock, The annealing of interstitial-related optical centres in type ii natural and cvd dia- mond, Diamond and Related Materials7, 228 (1998)

  18. [18]

    Naydenov, R

    B. Naydenov, R. Kolesov, A. Batalov, J. Meijer, S. Pezza- gna, D. Rogalla, F. Jelezko, and J. Wrachtrup, Engineering sin- gle photon emitters by ion implantation in diamond, Applied Physics Letters95, 181109 (2009)

  19. [19]

    S. P. Turaga, H. Jin, E. J. Teo, and A. A. Bettiol, Cross-sectional hyperspectral imaging of proton implanted diamond, Applied Physics Letters115, 021904 (2019)

  20. [20]

    X. Chai, H. Liang, C. Yang, V . X. Ho, E. J. Teo, Z. Pastuovic, and A. Bettiol, High sensitivity spin defects in carbon implanted diamond, Advanced Optical Materials13, 2500397 (2025)

  21. [21]

    Davies and I

    G. Davies and I. I. service),Properties and growth of diamond, EMIS datareviews series (Institution of Electrical Engineers, 1994)

  22. [22]

    Newton, B

    M. Newton, B. Campbell, D. Twitchen, J. Baker, and T. An- thony, Recombination-enhanced diffusion of self-interstitial atoms and vacancy–interstitial recombination in diamond, Di- amond and Related Materials11, 618 (2002), 12th European Conference on Diamond, Diamond- Like Materials, Carbon Nanotubes, Nitrides & Silicon Carbide

  23. [23]

    S. J. Breuer and P. R. Briddon, Ab initio investigation of the native defects in diamond and self-diffusion, Phys. Rev. B51, 6984 (1995)

  24. [24]

    J. P. Goss, B. J. Coomer, R. Jones, C. J. Fall, P. R. Briddon, and S. ¨Oberg, Extended defects in diamond: The interstitial platelet, Phys. Rev. B67, 165208 (2003)

  25. [25]

    Gali, Ab initio theory of the nitrogen-vacancy center in dia- mond, Nanophotonics8, 1907 (2019)

    A. Gali, Ab initio theory of the nitrogen-vacancy center in dia- mond, Nanophotonics8, 1907 (2019)

  26. [26]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett.77, 3865 (1996)

  27. [27]

    J. P. Perdew and A. Zunger, Self-interaction correction to density-functional approximations for many-electron systems, Phys. Rev. B23, 5048 (1981)

  28. [28]

    Bernholc, A

    J. Bernholc, A. Antonelli, T. M. Del Sole, Y . Bar-Yam, and S. T. Pantelides, Mechanism of self-diffusion in diamond, Phys. Rev. Lett.61, 2689 (1988)

  29. [29]

    J. P. Goss, B. J. Coomer, R. Jones, C. J. Fall, C. D. Latham, P. R. Briddon, and S. ¨Oberg, Small aggregates of interstitials and models for platelets in diamond, Journal of Physics: Con- densed Matter12, 10257 (2000)

  30. [30]

    J. Goss, B. Coomer, R. Jones, T. Shaw, P. Briddon, and S. ¨Oberg, Interstitial aggregates in diamond, Diamond and Re- lated Materials10, 434 (2001), 11th European Conference on Diamond, Diamond-like Materials, Carbon Nanotubes, Nitrides 13 and Silicon Carbide

  31. [31]

    J. P. Goss, B. J. Coomer, R. Jones, T. D. Shaw, P. R. Briddon, M. Rayson, and S. ¨Oberg, Self-interstitial aggregation in dia- mond, Phys. Rev. B63, 195208 (2001)

  32. [32]

    Colasuonno, F

    F. Colasuonno, F. S. Gentile, W. Mackrodt, A. M. Ferrari, A. Platonenko, and R. Dovesi, Interstitial defects in diamond: A quantum mechanical simulation of their epr constants and vi- brational spectra, The Journal of Chemical Physics153, 024119 (2020)

  33. [33]

    Komarovskikh, Dft investigation of 100-split interstitial paramagnetic centers in diamond, Diamond and Related Ma- terials146, 111235 (2024)

    A. Komarovskikh, Dft investigation of 100-split interstitial paramagnetic centers in diamond, Diamond and Related Ma- terials146, 111235 (2024)

  34. [34]

    Salustro, Y

    S. Salustro, Y . N ¨oel, C. M. Zicovich-Wilson, P. Olivero, and R. Dovesi, The v + i defects in diamond: An ab initio investi- gation of the electronic structure, of the raman and ir spectra, and of their possible recombination, The Journal of Chemical Physics145, 184701 (2016)

  35. [35]

    Mainwood, A

    A. Mainwood, A. T. Collins, and P. Woad, Isotope dependence of the frequency of localized vibrational modes in diamond, in Defects in Semiconductors 17, Materials Science Forum, V ol. 143 (Trans Tech Publications Ltd, 1993) pp. 29–34

  36. [36]

    Kresse and J

    G. Kresse and J. Hafner, Ab initio molecular dynamics for liq- uid metals, Phys. Rev. B47, 558 (1993)

  37. [37]

    P. E. Bl ¨ochl, Projector augmented-wave method, Phys. Rev. B 50, 17953 (1994)

  38. [38]

    Paier, M

    J. Paier, M. Marsman, K. Hummer, G. Kresse, I. C. Gerber, and J. G. ´Angy´an, Screened hybrid density functionals applied to solids, The Journal of Chemical Physics124, 154709 (2006)

  39. [39]

    N. W. Ashcroft and N. D. Mermin, Solid state physics brooks, Cole, Cengage Learning10(1976)

  40. [40]

    Z. Fan, Z. Zeng, C. Zhang, Y . Wang, K. Song, H. Dong, Y . Chen, and T. Ala-Nissila, Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport, Phys. Rev. B104, 104309 (2021)

  41. [41]

    Erhart and K

    P. Erhart and K. Albe, Analytical potential for atomistic simu- lations of silicon, carbon, and silicon carbide, Phys. Rev. B71, 035211 (2005)

  42. [42]

    Z. Fan, Y . Wang, P. Ying, K. Song, J. Wang, Y . Wang, Z. Zeng, K. Xu, E. Lindgren, J. M. Rahm, A. J. Gabourie, J. Liu, H. Dong, J. Wu, Y . Chen, Z. Zhong, J. Sun, P. Erhart, Y . Su, and T. Ala-Nissila, Gpumd: A package for constructing accu- rate machine-learned potentials and performing highly efficient atomistic simulations, The Journal of Chemical Phy...

  43. [43]

    Y . Wang, Z. Fan, P. Qian, M. A. Caro, and T. Ala-Nissila, Density dependence of thermal conductivity in nanoporous and amorphous carbon with machine-learned molecular dynamics, Phys. Rev. B111, 094205 (2025)

  44. [44]

    H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. Di- Nola, and J. R. Haak, Molecular dynamics with coupling to an external bath, The Journal of Chemical Physics81, 3684 (1984)

  45. [45]

    J. E. Northrup, R. Di Felice, and J. Neugebauer, Energetics of h andnh 2 on gan(1010) and implications for the origin of nanopipe defects, Phys. Rev. B56, R4325 (1997)

  46. [46]

    Freysoldt and J

    C. Freysoldt and J. Neugebauer, First-principles calculations for charged defects at surfaces, interfaces, and two-dimensional materials in the presence of electric fields, Phys. Rev. B97, 205425 (2018)

  47. [47]

    Nelson, C

    R. Nelson, C. Ertural, J. George, V . L. Deringer, G. Hau- tier, and R. Dronskowski, Lobster: Local orbital projections, atomic charges, and chemical-bonding analysis from projector- augmented-wave-based density-functional theory, Journal of Computational Chemistry41, 1931 (2020)

  48. [48]

    D. C. Hunt, D. J. Twitchen, M. E. Newton, J. M. Baker, T. R. Anthony, W. F. Banholzer, and S. S. Vagarali, Identification of the neutral carbon⟨100⟩-split interstitial in diamond, Phys. Rev. B61, 3863 (2000)

  49. [49]

    Weigel, D

    C. Weigel, D. Peak, J. W. Corbett, G. D. Watkins, and R. P. Messmer, Carbon interstitial in the diamond lattice, Phys. Rev. B8, 2906 (1973)

  50. [50]

    D. J. Twitchen, M. E. Newton, J. M. Baker, O. D. Tucker, T. R. Anthony, and W. F. Banholzer, Electron-paramagnetic- resonance measurements on the di-⟨001⟩-split interstitial center (r1) in diamond, Phys. Rev. B54, 6988 (1996)

  51. [51]

    Coomer, J

    B. Coomer, J. Goss, R. Jones, S. ¨Oberg, and P. Briddon, Inter- stitial aggregates and a new model for the i1/w optical centre in silicon, Physica B: Condensed Matter273-274, 505 (1999)

  52. [52]

    Jones, T

    R. Jones, T. Eberlein, N. Pinho, B. Coomer, J. Goss, P. Brid- don, and S. ¨Oberg, Self-interstitial clusters in silicon, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms186, 10 (2002)

  53. [53]

    P. J. Fallon, L. M. Brown, J. C. Barry, and J. Bruley, Nitrogen determination and characterization in natural diamond platelets, Philosophical Magazine A72, 21 (1995)

  54. [54]

    A. L. G ¨orne and R. Dronskowski, Covalent bonding versus to- tal energy: On the attainability of certain predicted low-energy carbon allotropes, Carbon148, 151 (2019)

  55. [55]

    Henkelman, A

    G. Henkelman, A. Arnaldsson, and H. J ´onsson, A fast and ro- bust algorithm for bader decomposition of charge density, Com- putational Materials Science36, 354 (2006)

  56. [56]

    Prawer, K

    S. Prawer, K. Nugent, and D. Jamieson, The raman spectrum of amorphous diamond, Diamond and Related Materials7, 106 (1998)

  57. [57]

    A. A. El-Azhary and H. U. Suter, Comparison between opti- mized geometries and vibrational frequencies calculated by the dft methods, The Journal of Physical Chemistry100, 15056 (1996)

  58. [58]

    G. P. Srivastava,The Physics of Phonons, 2nd ed. (CRC Press, Boca Raton, 2022)

  59. [59]

    Shimanouchiet al.,Tables of molecular vibrational frequen- cies(US Government Printing Office Washington, DC, 1973)

    T. Shimanouchiet al.,Tables of molecular vibrational frequen- cies(US Government Printing Office Washington, DC, 1973)