From Mono- to Hexa-Interstitials: Computational Insights into Carbon Defects in Diamond
Pith reviewed 2026-05-17 00:20 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [Abstract] Abstract: convergence tests, functional benchmarks, and error estimates on formation energies are absent, which should be added for clarity.
- [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
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
-
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
-
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
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
free parameters (2)
- Exchange-correlation functional
- Interatomic potential parameters
axioms (2)
- standard math Born-Oppenheimer approximation for separating nuclear and electronic motion
- domain assumption Periodic supercell boundary conditions with finite-size corrections
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We calculated the formation energy per interstitial of the optimized defects (cf. Fig. 1) to show their (relative) stability in different charge states (see Fig. 3).
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The stability of various defect forms was evaluated based on the calculated formation energies. Specifically, we computed the formation energy E_q_f ...
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
-
Diamond crystal with Y-defects: spectroscopy and transmission electron microscopy
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
-
[1]
R. Kalish, Diamond as a unique high-tech electronic mate- rial: difficulties and prospects, Journal of Physics D: Applied Physics40, 6467 (2007)
work page 2007
-
[2]
R. J. Nemanich, J. A. Carlisle, A. Hirata, and K. Haenen, Cvd diamond—research, applications, and challenges, MRS Bul- letin39, 490–494 (2014)
work page 2014
- [3]
-
[4]
D. S. Misra, Nitrogen centers in single crystal diamond grown by chemical vapor deposition, Journal of Vacuum Science & Technology A41, 042703 (2023)
work page 2023
-
[5]
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)
work page 2020
-
[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)
work page 2009
- [7]
-
[8]
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)
work page 2021
- [9]
-
[10]
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)
work page 2008
-
[11]
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)
work page 2022
-
[12]
K. Iakoubovskii and G. Adriaenssens, Optical study of some interstitial-related centres in cvd diamond, physica status solidi (a)181, 59 (2000)
work page 2000
-
[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)
work page 1978
-
[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)
work page 1956
-
[15]
J. Walker, An optical study of the tr12 and 3h defects in irra- diated diamond, Journal of Physics C: Solid State Physics10, 3031 (1977)
work page 1977
- [16]
- [17]
-
[18]
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)
work page 2009
-
[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)
work page 2019
-
[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)
work page 2025
-
[21]
G. Davies and I. I. service),Properties and growth of diamond, EMIS datareviews series (Institution of Electrical Engineers, 1994)
work page 1994
-
[22]
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
work page 2002
-
[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)
work page 1995
-
[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)
work page 2003
-
[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)
work page 1907
-
[26]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett.77, 3865 (1996)
work page 1996
-
[27]
J. P. Perdew and A. Zunger, Self-interaction correction to density-functional approximations for many-electron systems, Phys. Rev. B23, 5048 (1981)
work page 1981
-
[28]
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)
work page 1988
-
[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)
work page 2000
-
[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
work page 2001
-
[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)
work page 2001
-
[32]
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)
work page 2020
-
[33]
A. Komarovskikh, Dft investigation of 100-split interstitial paramagnetic centers in diamond, Diamond and Related Ma- terials146, 111235 (2024)
work page 2024
-
[34]
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)
work page 2016
-
[35]
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
work page 1993
-
[36]
G. Kresse and J. Hafner, Ab initio molecular dynamics for liq- uid metals, Phys. Rev. B47, 558 (1993)
work page 1993
-
[37]
P. E. Bl ¨ochl, Projector augmented-wave method, Phys. Rev. B 50, 17953 (1994)
work page 1994
- [38]
-
[39]
N. W. Ashcroft and N. D. Mermin, Solid state physics brooks, Cole, Cengage Learning10(1976)
work page 1976
-
[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)
work page 2021
-
[41]
P. Erhart and K. Albe, Analytical potential for atomistic simu- lations of silicon, carbon, and silicon carbide, Phys. Rev. B71, 035211 (2005)
work page 2005
-
[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...
work page 2022
-
[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)
work page 2025
-
[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)
work page 1984
-
[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)
work page 1997
-
[46]
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)
work page 2018
-
[47]
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)
work page 1931
-
[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)
work page 2000
- [49]
-
[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)
work page 1996
- [51]
- [52]
-
[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)
work page 1995
-
[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)
work page 2019
-
[55]
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)
work page 2006
- [56]
-
[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)
work page 1996
-
[58]
G. P. Srivastava,The Physics of Phonons, 2nd ed. (CRC Press, Boca Raton, 2022)
work page 2022
-
[59]
T. Shimanouchiet al.,Tables of molecular vibrational frequen- cies(US Government Printing Office Washington, DC, 1973)
work page 1973
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.