An extended ab initio theory of the V$_{\text{B}}^-$ center in hBN: excited states, Jahn-Teller distortion, and pressure dependence

This paper employs high-level CASSCF-NEVPT2 calculations to comprehensively model the excited states, structural distortions, and stress-dependent properties of the V$_{\text{B}}^-$ center in hBN, thereby clarifying its complex optical cycle and establishing a theoretical foundation for advanced 2D quantum sensing applications.

The Gist

Imagine a tiny, invisible "quantum lightbulb" hidden inside a sheet of hexagonal boron nitride (hBN), which is essentially a super-thin, atomically flat layer of material. This lightbulb is a specific defect called the negatively charged boron vacancy ($V^-_B$). Scientists are excited about it because it can act as a sensor for magnetic fields and other tiny forces, working at room temperature and even fitting into ultra-thin 2D devices.

However, for a long time, scientists didn't fully understand how this lightbulb works. They knew it glowed and reacted to magnetic fields, but the internal mechanics were a mystery because the electrons involved are "strongly correlated"—a fancy way of saying they dance together in a complex, chaotic way that standard computer models can't easily predict.

This paper acts like a high-resolution manual, using advanced computer simulations to finally explain the inner workings of this quantum lightbulb. Here is the breakdown of their findings using simple analogies:

1. The Shape-Shifting Dance (Jahn-Teller Distortion)

When the lightbulb gets excited by a laser (like a green light), it doesn't just sit still. Imagine a perfectly round, equilateral triangle made of three nitrogen atoms. When the electron gets excited, this triangle suddenly gets "stretched" in one direction, turning into a lopsided shape.

• The Paper's Claim: This stretching is called a Jahn-Teller distortion. It's not a tiny wobble; it's a major structural change. The triangle becomes so distorted that it creates a "tricorn hat" shape in the energy landscape (imagine a hat with three distinct valleys).
• The Consequence: At low temperatures (below 200 K), the triangle gets "stuck" in one of these three valleys (a static state). But at room temperature, it has enough energy to hop between the valleys rapidly (a dynamic state). This hopping changes how the lightbulb behaves and how it splits its magnetic signals.

2. The "Ghost" of a Missing Atom

The defect is created because a boron atom is missing. This leaves behind six "dangling" electron orbitals on the neighboring nitrogen atoms.

• The Paper's Claim: The authors mapped out the energy levels of these electrons. They found that the lightbulb absorbs green light (about 2.3 eV) to get excited. However, when it relaxes back down, it doesn't just glow at a single sharp color. Instead, it emits a broad, fuzzy glow (a "phonon sideband") because the shape change is so drastic that it kicks out about five "sound waves" (phonons) for every photon of light it emits.
• The Result: The "pure" color of the light (the Zero-Phonon Line) is so faint (only 0.4% of the total light) that it's almost invisible, buried under the broad, fuzzy glow. This explains why experiments have struggled to see a sharp color peak.

3. The Secret Tunnel (Intersystem Crossing)

The magic of this lightbulb for sensing lies in its ability to switch between different "spin" states (think of these as different orientations of a tiny internal compass).

• The Paper's Claim: The authors discovered that the path the electron takes to switch spins depends heavily on its orientation ($m_S = 0$ vs. $m_S = \pm 1$).
  • One path is fast and direct.
  • The other path involves a "quasi-degenerate" state, where a singlet state (one type of spin) and a triplet state (another type) are so close in energy they almost touch.
• The Analogy: Imagine two parallel train tracks that are so close together the train can jump between them easily if the track shakes (vibrates). This "jumping" (Intersystem Crossing) is what allows the device to be read out optically. The paper suggests this jumping is highly sensitive to temperature and pressure.

4. Squeezing the Lightbulb (Pressure and Strain)

The researchers also tested what happens if you squeeze the material (apply pressure).

• The Paper's Claim:
  • Squeezing from the top (Vertical pressure): This makes the layers of the material closer together. It speeds up the "spin-jumping" process significantly, making the lightbulb dimmer and its lifetime shorter.
  • Squeezing from the sides (Horizontal pressure): This changes the magnetic "splitting" (the D parameter) of the ground state.
• The Takeaway: The lightbulb is a very sensitive strain gauge. How it reacts to pressure depends on which way you squeeze it. The paper confirms that the changes in the magnetic signal under pressure are due to the physical compression of the atomic lattice.

5. What the Paper Doesn't Say

It is important to note what this paper does not claim:

• It does not claim to have built a working commercial sensor yet.
• It does not claim to have solved every mystery. The authors admit that the transition from the "zero spin" state to the singlet state is still too complex for their current models to calculate perfectly. They suggest that future work needs even more advanced simulation methods to fully understand that specific "jump."
• It does not discuss clinical uses or medical applications.

Summary

In short, this paper uses super-advanced computer modeling to draw a detailed map of the $V^-_B$ center. It explains that this quantum defect is a shape-shifter that distorts its own atomic structure when excited, creating a complex energy landscape. This distortion dictates how it glows, how it switches its magnetic spin, and how it reacts to being squeezed. This theoretical map provides the foundation needed to turn this defect into a reliable tool for nanoscale quantum sensing.

Technical Summary

Technical Summary: An Extended Ab Initio Theory of the V⁻B Center in hBN

Problem Statement
Negatively charged boron vacancy (V⁻B) centers in hexagonal boron nitride (hBN) have emerged as promising two-dimensional spin qubits for nanoscale quantum sensing, offering advantages over bulk diamond systems regarding surface proximity and spatial resolution. However, a comprehensive theoretical description of their optically detected magnetic resonance (ODMR) signal remains elusive. This challenge stems from the strongly correlated, multireference nature of the electronic states involved in the optical cycle. Previous theoretical efforts have struggled to fully capture the excited state fine structure, Jahn-Teller (JT) distortions, and intersystem crossing (ISC) pathways, limiting the ability to optimize the V⁻B center for integrated quantum sensors.

Methodology
To address these limitations, the authors employ a high-level, wavefunction-based ab initio approach, specifically Complete Active Space Self-Consistent Field (CASSCF) augmented by n-electron valence perturbation theory (NEVPT2). This method is chosen for its ability to accurately describe both static and dynamic electron correlation effects, which are critical for the V⁻B center's multireference character.

• Model System: Calculations are performed on a B₆₀N₆₀H₂₇ cluster model representing the V⁻B defect.
• Electronic Structure: The active space includes six defect orbitals localized on the first-neighbor nitrogen atoms (three in-plane $a'_1$ and $e'$ states, and three out-of-plane $a''_2$ and $e''$ states).
• Dynamics and Rates: Structural relaxations are computed at the CASSCF level. Transition rates for photoluminescence (PL) and intersystem crossing (ISC) are approximated using a 1D effective phonon model. Zero-field splitting (ZFS) parameters are derived via quasi-degenerate perturbation theory (QDPT).
• External Perturbations: The effects of hydrostatic pressure and unidirectional strain are modeled by deforming the molecular cluster geometry (adjusting interlayer distances and in-plane lattice constants) and re-evaluating electronic properties.

Key Contributions and Results

1. Electronic Structure and State Ordering:
   The study establishes the energy spectrum of the V⁻B center, identifying the ground state as $^3\bar{A}'_2$ ($D_{3h}$ symmetry). The first optically allowed transition occurs to the $^3\bar{E}'$ state. Upon excitation, the system undergoes rapid internal conversion to lower-lying triplet states. Crucially, the $^3\bar{E}''$ state undergoes a substantial pseudo-Jahn-Teller (pJT) distortion, lowering its symmetry to $C_{2v}$ and splitting into a lower-energy $^3\bar{A}_2$ state and a higher-energy $^3\bar{B}_2$ state. The $^3\bar{A}_2$ state is identified as the lowest-energy triplet excited state.

2. Jahn-Teller Distortion and Potential Energy Surface:
   The $^3\bar{E}'' \to ^3\bar{A}_2$ relaxation is driven by a strong pJT effect, resulting in a stabilization energy of 355 meV and a structural distortion where the inner nitrogen triangle stretches (N-N distances: 2.49 Å and 2.65 Å). The potential energy surface is characterized as a "tricorn hat" with three minima separated by a barrier of $\delta_{JT} \approx 308$ meV.

   • Temperature Dependence: Below ~200 K, the system behaves as a static JT system in a $C_{2v}$ configuration. Above 200 K, transitions between the three JT minima occur within the excited state lifetime, leading to a dynamic JT regime. This dynamic behavior explains the temperature dependence of the transverse zero-field splitting (ZFS) parameter ($E$).

3. Optical Properties and Zero-Phonon Line (ZPL):
   The dominant radiative pathway is the $^3\bar{A}_2 \to ^3\bar{A}'_2$ transition. The calculated ZPL energy is 1.88 eV. However, the Huang-Rhys factor is calculated to be $S = 5.4$, indicating significant electron-phonon coupling and a large phonon sideband (PSB). The Debye-Waller factor is extremely low ($f_{DW} \approx 0.4\%$), suggesting the ZPL is obscured by the PSB, consistent with experimental observations of a broad emission spectrum. The calculated PSB maximum (1.60 eV) aligns well with experimental data.

4. Intersystem Crossing (ISC) and Spin Polarization:
   The study details distinct ISC pathways for the $m_S = 0$ and $m_S = \pm 1$ spin sublevels, contrasting with the NV center in diamond.

   • The $m_S = \pm 1$ states decay rapidly ($\tau \approx 3.2$ ns) to the singlet sector via the $^3\bar{A}_2 \to ^1\bar{E}'$ channel.
   • The $m_S = 0$ channel involves quasi-degenerate singlet-triplet states ($^3\bar{A}_2$ and $^1\bar{A}_2$), making the transition rate highly sensitive to temperature and strain.
   • The system populates a metastable singlet shelving state ($^1\bar{E}'$) with a lifetime of ~340 ns before decaying back to the ground state, facilitating spin initialization and readout.

5. Strain and Pressure Dependence:
   The authors model the response of ODMR parameters to hydrostatic pressure and unidirectional strain.

   • Ground State: The $D$ tensor of the ground state ($^3\bar{A}'_2$) shows a positive shift with pressure, primarily driven by in-plane (horizontal) compression rather than interlayer (vertical) compression. The calculated slope (~38 MHz/GPa) agrees qualitatively with experimental values.
   • Excited State Dynamics: The ISC rate for the $m_S = \pm 1$ channel is highly pressure-dependent, increasing by approximately a factor of two under 2.7 GPa of hydrostatic pressure. This is attributed to the sensitivity of the spin-orbit coupling matrix element and the reduction of the energy gap under vertical compression. This finding theoretically explains the experimentally observed decrease in photoluminescence intensity and excited-state lifetime under pressure.

Significance and Claims
The paper claims to provide the first extended, high-level theoretical framework for the V⁻B center that successfully reconciles complex multireference electronic effects with experimental ODMR observations. By clarifying the role of the pseudo-Jahn-Teller effect, the static-to-dynamic transition of the excited state, and the distinct spin-dependent decay pathways, the work establishes a theoretical foundation for understanding the fundamental behavior of the V⁻B center.

The authors modestly conclude that while their results agree well with experimental ZFS and photoluminescence parameters, a complete quantitative description of non-radiative decay (particularly the $m_S=0$ channel) requires further development of multi-surface dynamic simulation methods beyond the Born-Oppenheimer approximation. Nevertheless, the study provides crucial theoretical support for the deployment of V⁻B centers as integrated 2D quantum sensors, particularly in interpreting strain and pressure sensing capabilities.