Bright oxygen- and vacancy-derived spin-singlet diamond color centers with metastable spin triplets: OV$^{2+}$ and VOV$^{2+}$

Using embedded multiconfigurational quantum mechanical theory, this study identifies the long-sought ST1 diamond color center as the metastable spin-triplet, spin-singlet ground-state V$_C$O$_C$V$_{C}^{2+}$ defect, whose calculated optical transitions closely match experimental observations.

The Gist

The Big Picture: Finding the "Ghost" in the Diamond

Imagine a diamond not just as a shiny rock, but as a giant, perfectly organized city made of carbon atoms. Sometimes, the city has "construction sites" (vacancies) where a building is missing, or "foreigners" (impurities) like oxygen atoms living in the wrong houses.

Scientists have been hunting for a specific type of "ghost" in this city for over ten years. They call it the ST1 defect. They know this ghost exists because it glows with a specific color of light and has a special "magnetic heartbeat" (spin) that makes it perfect for building quantum computers and super-sensitive sensors.

However, despite years of searching, no one could figure out exactly what this ghost looked like. Was it one oxygen atom next to one empty spot? Two empty spots? A specific arrangement? The paper you shared is the moment the detective finally solves the case.

The Detective Work: Trying on Different Outfits

The author, John Mark P. Martirez, acted like a digital tailor. He used powerful computer simulations to "try on" different outfits (structures) for the ST1 ghost to see which one fit the experimental clues.

He focused on two main suspects:

1. Suspect A (OCVC): One oxygen atom sitting right next to one empty spot (vacancy). Think of this as a lonely guest sitting next to an empty chair.
2. Suspect B (VCOCVC): One oxygen atom sitting in the middle of two empty spots. Think of this as a guest sitting between two empty chairs, holding hands with both sides.

The Clues: Why Suspect B is the Winner

The author compared these two suspects against five major clues gathered from real-world experiments:

1. The "Ground State" (The Default Mood)

• The Clue: The ST1 ghost is naturally calm and quiet (a "spin singlet"), but it can get excited into a magnetic state (a "triplet").
• The Verdict: Suspect A (OCVC) is naturally magnetic (a triplet) when neutral, which doesn't match. Suspect B (VCOCVC), however, is naturally calm (a singlet) and can easily switch to a magnetic triplet. Suspect B wins.

2. The "Glow" (Optical Color)

• The Clue: When you shine a light on ST1, it absorbs energy at a very specific color (around 2.2 to 2.3 electron-volts) and glows back.
• The Verdict: Suspect A absorbs light that is too blue (too much energy). Suspect B absorbs light at the exact right color that matches the experiment. Suspect B wins.

3. The "Magnetic Signature" (Zero-Field Splitting)

• The Clue: When the ghost is in its magnetic triplet state, it has a specific "fingerprint" in magnetic resonance tests. Crucially, this fingerprint shows it is not perfectly symmetrical (it's lopsided).
• The Verdict: Suspect A is perfectly symmetrical (like a tripod), so its fingerprint would be boring and symmetrical. Suspect B is lopsided (like a seesaw), creating the exact messy, non-symmetrical fingerprint seen in the lab. Suspect B wins.

4. The "Social Life" (Thermodynamics)

• The Clue: Nature prefers things that are stable and easy to form.
• The Verdict: The paper shows that oxygen atoms love to sit between two empty spots. It's like a magnet; the oxygen has "lonely hands" (lone pairs of electrons) that desperately want to grab onto two empty spots at once. This makes Suspect B form naturally and easily, while Suspect A is less stable. Suspect B wins.

5. The Charge State

• The Clue: The ghost needs to be positively charged (+2) to show its best behavior.
• The Verdict: The simulations show that Suspect B (VCOCVC) can happily lose two electrons to become a +2 charged defect, whereas the other options struggle to do so.

The Final Solution: The "Double-Vacancy" Diamond

The paper concludes that the mysterious ST1 defect is actually two empty spots with an oxygen atom sandwiched right in the middle, carrying a +2 charge.

The Analogy:
Imagine a row of three seats in a theater.

• Suspect A is a person sitting in the middle seat with an empty seat next to them.
• Suspect B is a person sitting in the middle seat, but both seats next to them are empty.

The paper proves that the "ghost" is Suspect B. The oxygen atom acts like a bridge connecting two empty spaces. This specific arrangement creates the perfect conditions for the diamond to act as a quantum bit (qubit), capable of storing information without losing it to noise.

Why Does This Matter?

Now that we know exactly what this defect looks like, scientists can stop guessing and start engineering. They can:

• Build better Quantum Computers: By creating more of these specific "ghosts," we can build faster, more stable quantum processors.
• Create Super-Sensors: These defects are so sensitive that they can detect tiny magnetic fields, potentially allowing us to image single molecules or even brain activity with incredible precision.

In short, this paper solved a 10-year-old mystery, turning a blurry photo of a quantum defect into a crystal-clear blueprint for the future of technology.

Technical Summary

1. Problem Statement

The ST1 diamond color center is a promising candidate for quantum information applications (qubits and quantum sensors) due to its spin singlet ground state, metastable spin triplet excited state, and lack of nuclear spin decoherence (composed of $^{16}$O and $^{12}$C). Discovered in 2013, ST1 exhibits long spin resonance lifetimes and high-contrast optical readout. However, despite over a decade of experimental and theoretical efforts, its exact atomic structure and charge state remain unsolved.

Previous hypotheses suggested various oxygen-vacancy complexes (e.g., $O_C$, $O_CV_C$, $O_CV_CV_C$, or hydrogenated variants), but none could simultaneously explain ST1's specific optical properties (Zero Phonon Line at 2.2–2.3 eV), magnetic properties (Zero-Field Splitting parameters), and symmetry constraints. The primary challenge is identifying a defect model that possesses a spin-singlet ground state, a metastable triplet, and optical transitions matching the experimental Zero Phonon Line (ZPL).

2. Methodology

The authors employed a rigorous, multi-scale quantum mechanical modeling approach to investigate two primary candidate configurations:

1. $O_CV_C$: A substitutional oxygen atom adjacent to a single carbon vacancy (1:1 ratio).
2. $V_CO_CV_C$: A substitutional oxygen atom flanked by two carbon vacancies along the [110] axis (1:2 ratio).

Computational Framework:

• Embedding Theory: To accurately model isolated defects within the bulk diamond crystal, the authors used Capped Density Functional Embedding Theory (capped-DFET). This involves carving finite clusters ($C_{15}O$ for $O_CV_C$ and $C_{26}O$ for $V_CO_CV_C$) from periodic supercells and capping the edges with F and O atoms.
• Environment Potential: The effect of the periodic diamond environment on the clusters was modeled using an optimized effective potential derived from DFT-HSE06 calculations on the full periodic system.
• Excited State Theory: To capture the multiconfigurational nature of the defects (essential for accurate spin states and excitation energies), the authors performed embedded Complete Active Space Self-Consistent Field (emb-CASSCF) calculations followed by N-Electron Valence State Perturbation Theory (emb-NEVPT2).
• Properties Calculated:
  • Vertical Excitation Energies (VEEs) and Zero Phonon Lines (ZPL).
  • Transition dipole moments ($\mu_t$), oscillator strengths ($f$), and spontaneous emission lifetimes ($\tau$).
  • Zero-Field Splitting (ZFS) parameters ($D_{ZFS}$ and $E_{ZFS}$) derived from spin-spin dipolar interactions.
  • Thermodynamic coupling energies and ionization potentials.

3. Key Contributions

• Structural Identification: The study definitively identifies the ST1 defect as the doubly positively charged $V_CO_CV_C^{2+}$ complex (denoted as $VOV^{2+}$ in the title, though chemically $V_CO_CV_C^{2+}$).
• Mechanism of Formation: The authors explain why the 1:2 vacancy-to-oxygen ratio is favored over the 1:1 ratio, attributing it to the electronic stabilization provided by the two oxygen lone pairs, analogous to how a nitrogen lone pair stabilizes the $NV^-$ center.
• Charge State Stabilization: The paper elucidates the thermodynamic conditions required to stabilize the $+2$ charge state, specifically the role of hole-doped environments or electron traps in facilitating ionization.

4. Key Results

A. Structural and Electronic Properties

• Symmetry: The $V_CO_CV_C^{2+}$ defect adopts $C_{2v}$ symmetry, consistent with experimental observations of the ST1 axis. In contrast, the $O_CV_C^{2+}$ candidate has $C_{3v}$ symmetry.
• Ground State: Both $O_CV_C^{2+}$ and $V_CO_CV_C^{2+}$ possess a spin-singlet ground state ($^1A_1$), satisfying the requirement for ST1. However, the neutral $O_CV_C^0$ is a triplet, and $V_CO_CV_C^0$ has nearly degenerate singlet/triplet states, making them less likely candidates for the specific ST1 behavior observed.

B. Optical Properties (The "Smoking Gun")

• Excitation Energies:
  • $V_CO_CV_C^{2+}$: The calculated vertical excitation energies for the first two bright singlet excited states ($^1B_2$ and $^1A_1$) are 2.38 eV and 2.52 eV (or 2.29/2.47 eV with different state averaging). These values closely match the experimental ST1 ZPL of 2.2–2.3 eV.
  • $O_CV_C^{2+}$: The first excitation energy is 2.83 eV, which is significantly higher than the experimental ST1 ZPL, effectively ruling out this structure.
• Brightness and Lifetime: The $V_CO_CV_C^{2+}$ transitions are bright (oscillator strengths $f \approx 0.20–0.36$) with calculated lifetimes of 4–9 ns, consistent with experimental readout capabilities.

C. Magnetic Properties (ZFS)

• Zero-Field Splitting: Experimental ST1 exhibits non-zero $E_{ZFS}$ ($0.139$ GHz), indicating a lack of axial symmetry.
  • $O_CV_C^{2+}$: Due to its $C_{3v}$ symmetry, it predicts $E_{ZFS} = 0$, which contradicts experimental data.
  • $V_CO_CV_C^{2+}$: The calculated ZFS parameters for its metastable triplet states ($^3B_2$ and $^3A_1$) yield $D_{ZFS} \approx 2.6$ GHz and $E_{ZFS} \approx 0.9$ GHz (depending on the specific triplet state and averaging). While the absolute values vary, the presence of a significant non-zero $E_{ZFS}$ and the alignment of the principal axes with the [110] direction strongly support the $V_CO_CV_C^{2+}$ model.

D. Thermodynamics and Ionization

• Coupling Energy: The formation of $V_CO_CV_C$ from $O_CV_C$ and a vacancy is highly exothermic ($-4.28$ eV), suggesting that under conditions of high vacancy concentration, the 1:2 complex is thermodynamically favored over the 1:1 complex.
• Ionization: The study shows that while $V_CO_CV_C$ is harder to ionize to $+1$ than $O_CV_C$, the second ionization to $+2$ is energetically accessible in hole-doped diamond or the presence of electron traps (like Boron). This explains why the $+2$ state is observed in specific experimental preparations (e.g., reactive ion etching with oxygen plasma).

5. Significance

• Resolution of a Decade-Old Mystery: This work provides the first theoretically consistent and experimentally validated identification of the ST1 defect structure as $V_CO_CV_C^{2+}$.
• Quantum Technology Impact: By confirming the structure, the study enables the rational engineering of ST1 centers for quantum computing and sensing. Understanding the charge state stability ($+2$) allows researchers to optimize synthesis protocols (e.g., controlling vacancy concentration and doping) to maximize ST1 yield.
• Methodological Benchmark: The successful application of capped-DFET with embedded multiconfigurational wavefunction theory demonstrates a robust pathway for solving complex defect problems in wide-bandgap semiconductors where standard DFT fails to capture excited-state physics and spin multiplicities accurately.
• Design Principles: The finding that oxygen's lone pairs favor the formation of double-vacancy complexes ($V_CO_CV_C$) over single-vacancy complexes ($O_CV_C$) offers new design principles for engineering other oxygen-based quantum defects in diamond.