Negatively Charged Boron Vacancy (VB-)

• V_B^- is a point defect in hBN created by a missing boron atom that traps an extra electron, resulting in unique spin-dependent photoluminescence.
• It exhibits distinct ground and excited spin Hamiltonians with measurable zero-field splitting and hyperfine interactions that enable ODMR-based quantum control.
• Engineering strategies like plasmonic enhancement and charge-state control improve its optical brightness and coherence, advancing applications in quantum sensing and photonics.

The negatively charged boron vacancy, $V_B^-$, is a point defect in hexagonal boron nitride (hBN) formed by a missing boron atom that traps an extra electron. In hBN it is the only boron-vacancy charge state known to show spin-dependent photoluminescence and therefore functions as the principal spin-photon interface among boron-vacancy configurations. Its combination of a room-temperature spin-triplet ground state, optically detected magnetic resonance (ODMR), broad red emission near $\sim 800\ \text{nm}$, and compatibility with atomically thin van der Waals devices has made it a central defect platform for quantum sensing, integrated photonics, and spin-based quantum control in hBN (Gale et al., 2023, Hennessey et al., 2023).

1. Defect identity, symmetry, and optical character

Structurally, $V_B^-$ is a boron vacancy surrounded by three neighboring nitrogen atoms. The defect is commonly described as having $D_{3h}$ point-group symmetry, with a ground triplet manifold $^3A_2'$, an excited triplet manifold $^3E''$, and a metastable singlet $^1E'$ in the three-manifold picture used for recent polarization and photodynamics analyses (Lee et al., 2024, Geng et al., 31 Aug 2025). The host material, hBN, is a wide-bandgap van der Waals crystal whose atomically thin geometry places defect spins in extreme proximity to external samples, which is why $V_B^-$ has been emphasized for nanoscale sensing of magnetic fields, temperature, pressure, strain, and nearby materials.

The optical response of $V_B^-$ is unusual among solid-state spin defects because the free-space photoluminescence is broad and largely featureless rather than being dominated by a narrow zero-phonon line. This broadening has been attributed to strong electron-phonon coupling and Jahn-Teller mixing of electronic states, and the emission is usually described as peaking around $800\ \text{nm}$ (Qian et al., 2022). Within cavity-enhanced measurements, the zero-phonon line was directly identified at $\sim 800\ \text{nm}$0 nm, in agreement with earlier theoretical estimates reported by Ivady et al., Reimers et al., and Libbi et al.; the dominant transition was assigned to $\sim 800\ \text{nm}$1 (Qian et al., 2022).

The defect’s basic spin signature is a triplet ground state with zero-field splitting of about $\sim 800\ \text{nm}$2. This enables optical initialization, microwave manipulation, and optical readout through ODMR. A recurring misconception in the broader boron-vacancy literature is to treat different charge states interchangeably. The available experimental record instead distinguishes them sharply: only the negatively charged state is known to exhibit the spin-dependent photoluminescence used for ODMR-based readout, whereas the neutral state $\sim 800\ \text{nm}$3 does not provide the same operational spin-photon interface (Gale et al., 2023).

2. Spin Hamiltonians and electronic-state spectroscopy

The ground-state spin Hamiltonian has been formulated in the standard triplet-defect form

$\sim 800\ \text{nm}$4

with additional nuclear Zeeman and quadrupole terms used in ENDOR treatments (Gracheva et al., 2023). In zero or weak magnetic field, ODMR resonances are commonly parameterized by

$\sim 800\ \text{nm}$5

where $\sim 800\ \text{nm}$6 is the axial zero-field splitting and $\sim 800\ \text{nm}$7 is the transverse splitting (Guo et al., 2021). For ion-implanted ensembles at room temperature, two Lorentzian ODMR dips were reported near $\sim 800\ \text{nm}$8 MHz and $\sim 800\ \text{nm}$9 MHz, giving $V_B^-$0 and $V_B^-$1 (Guo et al., 2021).

Direct experimental access to the excited-state spin structure was obtained through continuous-wave ODMR and then validated by pulsed ODMR. In continuous-wave ODMR, the known ground-state resonances near $V_B^-$2 were accompanied by a broad additional resonance near $V_B^-$3; in pulsed ODMR, those broad $V_B^-$4 GHz features disappeared when the microwave acted only while the laser was off, confirming their excited-state origin (Yu et al., 2021). The excited-state spin was modeled as

$V_B^-$5

From field-dependent spectroscopy, the room-temperature excited-state parameters were extracted as $V_B^-$6, $V_B^-$7, and $V_B^-$8 (Yu et al., 2021). At cryogenic temperature, a closely related excited-state zero-field splitting of $V_B^-$9 was reported, together with excited-state ODMR contrast of $D_{3h}$0 at $D_{3h}$1 K; the excited-state $D_{3h}$2-factor was found to be similar to the ground-state value (Mu et al., 2021).

Manifold	Quantity	Reported value
Ground state	Zero-field splitting	about $D_{3h}$3
Excited state	Zero-field splitting	$D_{3h}$4 at room temperature
Excited state	Zero-field splitting	$D_{3h}$5 MHz at $D_{3h}$6 K
Excited state	$D_{3h}$7-factor	$D_{3h}$8
Excited state	Transverse anisotropy	$D_{3h}$9
Excited state	Hyperfine scales	$^3A_2'$0 MHz and $^3A_2'$1 MHz

Level anti-crossings are a central part of the $^3A_2'$2 spectroscopic phenomenology. Using the spin Hamiltonian, the excited-state LAC was predicted near $^3A_2'$3 G and the ground-state LAC near $^3A_2'$4 G at room temperature, while low-temperature measurements placed the excited-state LAC near $^3A_2'$5 G and the ground-state LAC near $^3A_2'$6 G (Yu et al., 2021, Mu et al., 2021). Near these crossings, strong spin mixing reduces the selectivity of optical pumping and lowers ODMR contrast; photoluminescence anomalies, angle dependence, and residual contrast suppression at nominal alignment were attributed to coherent spin precession and anisotropic relaxation, including residual mixing from hyperfine-induced effective tilted fields (Yu et al., 2021). A noteworthy contrast with the NV center in diamond is that in $^3A_2'$7 the excited-state LAC remains prominent at cryogenic temperature and the excited-state ODMR contrast becomes larger rather than disappearing (Mu et al., 2021).

3. Optical cycle, lifetimes, and spin polarization

Time-resolved photoluminescence established a spin-dependent optical cycle in which the excited-state decay rates differ strongly between the $^3A_2'$8 and $^3A_2'$9 manifolds. A seven-level description with ground and excited triplets plus an effective metastable singlet was used to analyze the dynamics. In this framework, optical pumping occurs at rate $^3E''$0, spin-conserving excited-state decay at rate $^3E''$1, intersystem crossing from the excited triplet to the singlet at rates $^3E''$2 and $^3E''$3, and singlet return to the ground triplet at rates $^3E''$4 and $^3E''$5 (Clua-Provost et al., 2024).

Using an all-optical protocol, room-temperature averages over 18 flakes yielded $^3E''$6, $^3E''$7, and spin polarization $^3E''$8. At $^3E''$9 K, the same analysis gave $^1E'$0, $^1E'$1, and $^1E'$2. From power-dependent transient fits, the metastable lifetime was inferred as $^1E'$3 at room temperature and $^1E'$4 at $^1E'$5 K (Clua-Provost et al., 2024). These data showed that the absolute lifetimes lengthen strongly on cooling, whereas the spin selectivity of the cycle remains substantial.

A complementary semiclassical analysis combined excited-state spectroscopy, power-dependent photoluminescence traces, and spin-resolved polarization dynamics to extract $^1E'$6, $^1E'$7, $^1E'$8, $^1E'$9, and $V_B^-$0 (Lee et al., 2024). On that basis, off-resonant optical pumping was predicted to produce electronic spin polarization $V_B^-$1 under ambient conditions. This suggests that the intrinsic polarization of $V_B^-$2 may be substantially higher than earlier rate models had implied, and it helps rationalize previously reported unusually large ODMR contrasts.

The metastable singlet lifetime was later measured directly in neutron-irradiated sub-micron flakes through pulse-pair recovery measurements. Averaging over 16 flakes gave $V_B^-$3 ns at room temperature (Escalante et al., 7 Apr 2025). In the same study, a conventional 7-level fit yielded $V_B^-$4 ns and failed to reproduce high-power quenching and the full thermal photoluminescence peak amplitude, whereas a 9-level model coupled to an additional 2-level subsystem produced $V_B^-$5 ns and matched the direct lifetime more closely. The additional subsystem was interpreted as another electronic state or possibly a charge-converted manifold, and larger flakes exhibited behavior consistent with optically induced conversion of $V_B^-$6 to another state, possibly $V_B^-$7 (Escalante et al., 7 Apr 2025). This is one of the clearest current model-refinement issues in the field: the seven-level picture is robust for many datasets, but it is not universally sufficient at higher excitation power or in all sample geometries.

4. Hyperfine structure, quadrupole couplings, and the nuclear-spin environment

Electron-nuclear couplings in $V_B^-$8 are dominated by the three nearest nitrogen atoms around the vacancy. Conventional ESR and high-frequency ENDOR established that the nearest-neighbor hyperfine interaction is axially symmetric, with principal values

$V_B^-$9

and that the nuclear quadrupole interaction is characterized by

$V_B^-$0

The hyperfine principal axis aligns with the nitrogen dangling-bond direction rather than the crystallographic $V_B^-$1-axis, whereas the quadrupole tensor is nearly axially symmetric about the same local axis (Gracheva et al., 2023).

A major consequence of those measurements is the spatial picture of the defect wavefunction. Using a linear-combination-of-atomic-orbitals analysis together with the measured hyperfine values, the spin density on a single nearest nitrogen was estimated as about $V_B^-$2, implying that about $V_B^-$3 of the electronic spin density is localized on the three nearest nitrogen atoms (Gracheva et al., 2023). This established experimentally that $V_B^-$4 is not a broadly delocalized spin texture but a strongly localized defect state confined mainly to the first coordination shell within a single BN layer.

Excited-state spectroscopy revealed a related but distinct hyperfine pattern. Fine scans of the excited-state ODMR line showed a seven-peak structure analogous to the ground-state pattern of three equivalent nitrogen nuclei, but with a much larger splitting of about $V_B^-$5 MHz; an additional finer splitting of about $V_B^-$6 MHz was tentatively attributed to a boron nucleus farther from the vacancy (Yu et al., 2021). The stronger nitrogen-related excited-state splitting was interpreted as evidence of greater electron density near the nitrogen nuclei in the excited state.

ENDOR has also resolved more remote nuclei. High-field measurements detected $V_B^-$7 spins in the third nitrogen shell, N(3), approximately $V_B^-$8 nm from the vacancy, with $V_B^-$9 MHz, $800\ \text{nm}$0 MHz, $800\ \text{nm}$1 MHz, $800\ \text{nm}$2 MHz, and $800\ \text{nm}$3 (Mamin et al., 9 Apr 2025). Density-functional calculations reproduced these values closely and confirmed the shell assignment. This extended the role of $800\ \text{nm}$4 from a defect whose own hyperfine structure can be resolved to a local probe of remote nuclear magnetic moments in the hBN host.

5. Generation, charge-state control, and optical-brightness engineering

Ion irradiation is the dominant route to creating $800\ \text{nm}$5, but the relation between vacancy production and usable optically active centers is nontrivial. Large-area ion implantation of exfoliated hBN produced room-temperature ODMR with $800\ \text{nm}$6, $800\ \text{nm}$7, ODMR contrast up to $800\ \text{nm}$8, and a longest measured $800\ \text{nm}$9 of about $\sim 800\ \text{nm}$00. In that study, nitrogen implantation at $\sim 800\ \text{nm}$01 keV and $\sim 800\ \text{nm}$02 generated a strong photoluminescence band from about $\sim 800\ \text{nm}$03 to $\sim 800\ \text{nm}$04 nm centered near $\sim 800\ \text{nm}$05 nm (Guo et al., 2021). Fluence, energy, and ion species were all found to matter, but not in the same way: fluence controlled defect density and disorder, energy mainly affected brightness within the thin-flake regime studied, and heavier ions increased strain-related transverse splitting and reduced $\sim 800\ \text{nm}$06.

A later focused-ion-beam framework emphasized that robust fabrication in thin flakes is limited at least as much by impurity recoil implantation as by nominal vacancy creation. In flakes thinner than about $\sim 800\ \text{nm}$07, surface and interface effects dominate. In a capped-versus-uncapped comparison using $\sim 800\ \text{nm}$08 keV H$\sim 800\ \text{nm}$09 at $\sim 800\ \text{nm}$10, a region capped by $\sim 800\ \text{nm}$11 nm multilayer graphene before irradiation showed about a $\sim 800\ \text{nm}$12 reduction in $\sim 800\ \text{nm}$13 emission, whereas a control region covered with $\sim 800\ \text{nm}$14 nm multilayer graphene only after irradiation showed about a $\sim 800\ \text{nm}$15 reduction attributable to optical absorption; the difference was assigned to recoil implantation of carbon into hBN (Hennessey et al., 2023). The same work identified an optimal fluence in a thick-flake example irradiated with $\sim 800\ \text{nm}$16 keV H$\sim 800\ \text{nm}$17, where the photoluminescence intensity peaked at about $\sim 800\ \text{nm}$18.

Quantifying the yield of optically active $\sim 800\ \text{nm}$19 remains a separate issue from counting all vacancies. By comparing fluence-dependent ODMR splittings with a microscopic charge model and molecular-dynamics vacancy counts for $\sim 800\ \text{nm}$20 keV He irradiation, a lower bound of about $\sim 800\ \text{nm}$21 was obtained for the fraction of all vacancies that enter the optically active negatively charged state (Carbone et al., 30 Jan 2025). This result strongly indicates that charge-state stabilization and environmental control are as important as vacancy production itself.

Charge-state control is therefore fundamental. Reversible switching between $\sim 800\ \text{nm}$22 and $\sim 800\ \text{nm}$23,

$\sim 800\ \text{nm}$24

was demonstrated under $\sim 800\ \text{nm}$25 keV electron-beam irradiation and $\sim 800\ \text{nm}$26 nm laser excitation (Gale et al., 2023). Electron-beam exposure quenched the $\sim 800\ \text{nm}$27 photoluminescence by about $\sim 800\ \text{nm}$28 at low flux and by more than $\sim 800\ \text{nm}$29 around $\sim 800\ \text{nm}$30, whereas the laser drove recovery toward the negative state. In an FLG/hBN/FLG heterostructure, the balance could be tuned electrically: at $\sim 800\ \text{nm}$31 V the quenching was about $\sim 800\ \text{nm}$32, at $\sim 800\ \text{nm}$33 V it increased to about $\sim 800\ \text{nm}$34, and without the electron beam the bias had no significant effect on photoluminescence (Gale et al., 2023). The physical interpretation was that holes promote $\sim 800\ \text{nm}$35 formation and electrons promote recovery of $\sim 800\ \text{nm}$36.

Because the intrinsic photoluminescence of $\sim 800\ \text{nm}$37 is weak, substantial effort has gone into brightness enhancement without destroying spin functionality. Several strategies have now been demonstrated. Low-loss plasmonic nano-patch antennas produced overall intensity enhancement up to $\sim 800\ \text{nm}$38, corresponding to an estimated actual enhancement of $\sim 800\ \text{nm}$39 after correcting for laser spot size, while preserving ODMR contrast (Xu et al., 2022). In suspended hBN, implantation-induced local deformation correlated with photoluminescence enhancement of up to $\sim 800\ \text{nm}$40 relative to supported regions, while the bright spot retained linewidth $\sim 800\ \text{nm}$41 MHz, ODMR contrast $\sim 800\ \text{nm}$42, and $\sim 800\ \text{nm}$43, comparable to darker supported regions (Geng et al., 31 Aug 2025). In PbI$\sim 800\ \text{nm}$44/hBN heterostructures, donor-assisted energy transfer increased $\sim 800\ \text{nm}$45 photoluminescence by $\sim 800\ \text{nm}$46–$\sim 800\ \text{nm}$47 and improved continuous-wave ODMR sensitivity by $\sim 800\ \text{nm}$48 at $\sim 800\ \text{nm}$49 nm and $\sim 800\ \text{nm}$50 at $\sim 800\ \text{nm}$51 nm, with a best sensitivity of about $\sim 800\ \text{nm}$52 at $\sim 800\ \text{nm}$53 (Mayner et al., 2 Feb 2026). Taken together, these results show that weak emission is not a fixed property of the defect alone; it can be strongly modified by photonic environment, strain symmetry breaking, and heterostructure-mediated energy transfer.

6. Relaxation, coherence, and quantum-technology roles

The longitudinal and transverse spin dynamics of $\sim 800\ \text{nm}$54 are strongly environment-dependent. At room temperature, single- and double-quantum relaxation measurements gave $\sim 800\ \text{nm}$55 and $\sim 800\ \text{nm}$56 at $\sim 800\ \text{nm}$57 K, with both rates increasing from $\sim 800\ \text{nm}$58 to $\sim 800\ \text{nm}$59 K and the double-quantum channel growing more rapidly (Xie et al., 17 Jun 2025). Above $\sim 800\ \text{nm}$60 K, the double-quantum rate was reported to become much greater than the single-quantum rate and may dominate the decoherence channel. A second-order spin-phonon model using effective phonon modes at $\sim 800\ \text{nm}$61, $\sim 800\ \text{nm}$62, and $\sim 800\ \text{nm}$63 meV reproduced the temperature dependence and implicated higher-energy phonons as especially important (Xie et al., 17 Jun 2025).

Field- and temperature-dependent relaxometry over a broader regime revealed three distinct relaxation regimes. Between $\sim 800\ \text{nm}$64 and $\sim 800\ \text{nm}$65 K and up to $\sim 800\ \text{nm}$66 T, low-temperature low-field behavior was dominated by spin-spin interactions and disorder-induced stretched exponential relaxation, high fields activated a first-order direct single-phonon process, and at $\sim 800\ \text{nm}$67–$\sim 800\ \text{nm}$68 a Raman-like two-phonon process with approximately $\sim 800\ \text{nm}$69 scaling dominated (Solanki et al., 22 Jul 2025). In the high-field regime the relaxation rate scaled approximately as

$\sim 800\ \text{nm}$70

and the spin transition frequency could be pushed from a few GHz to roughly $\sim 800\ \text{nm}$71 by magnetic field alone. This suggests a route toward broadband and sub-terahertz relaxometry using atomically thin hBN sensors.

Direct coherence measurements with broadband microwave control showed that short coherence remains a major limitation in many current samples. Using an isotopically enriched $\sim 800\ \text{nm}$72 thin film and sub-GHz Rabi driving, Ramsey interference gave $\sim 800\ \text{nm}$73 ns with Gaussian-like decay, while Hahn echo gave $\sim 800\ \text{nm}$74 ns and stretch factor $\sim 800\ \text{nm}$75 (Nakamura et al., 30 Apr 2026). The same work argued that strong, broadband pulses are necessary for reliable coherence extraction because the hyperfine-broadened ODMR spectrum otherwise produces substantial pulse errors.

First-principles many-body simulations indicate that the dominant decoherence mechanism changes with magnetic field in the dense nuclear-spin bath of hBN. A transition boundary was predicted at about $\sim 800\ \text{nm}$76 G for h-$\sim 800\ \text{nm}$77B$\sim 800\ \text{nm}$78N and $\sim 800\ \text{nm}$79 G for h-$\sim 800\ \text{nm}$80B$\sim 800\ \text{nm}$81N. Below that boundary, decoherence is governed by independent nuclear-spin dynamics and occurs within submicrosecond timescales; above it, pair-wise flip-flop transitions dominate and Hahn-echo $\sim 800\ \text{nm}$82 reaches tens of microseconds, with reported values at $\sim 800\ \text{nm}$83 T of $\sim 800\ \text{nm}$84 for h-$\sim 800\ \text{nm}$85B$\sim 800\ \text{nm}$86N and $\sim 800\ \text{nm}$87 for h-$\sim 800\ \text{nm}$88B$\sim 800\ \text{nm}$89N (Lee et al., 6 May 2025). This is a theoretical result, but it provides a concrete operating principle: isotope engineering and sufficiently large magnetic fields should move $\sim 800\ \text{nm}$90 into a more favorable coherence regime.

These relaxation and coherence properties frame the defect’s present and prospective applications. High optical polarization, including the model-based estimate $\sim 800\ \text{nm}$91, directly benefits sensing and initialization protocols (Lee et al., 2024). Excited-state level anti-crossings and hyperfine couplings are central to dynamic nuclear polarization, and the experimentally resolved remote $\sim 800\ \text{nm}$92N shells show that the defect can serve as a local NMR-like probe inside hBN itself (Yu et al., 2021, Mamin et al., 9 Apr 2025). On the quantum-information side, proposals for synchronous nuclear-spin control around the defect use the electron spin as a mediator for collective $\sim 800\ \text{nm}$93, $\sim 800\ \text{nm}$94, and Hadamard gates on the three nearest nuclei, with GHZ-state preparation reported at fidelity $\sim 800\ \text{nm}$95 in numerical analyses that include decoherence effects (Sakuldee et al., 2024). A plausible implication is that the long-term significance of $\sim 800\ \text{nm}$96 will depend not on any single metric—brightness, polarization, or coherence—but on the degree to which fabrication, charge-state stabilization, nuclear-spin control, and photonic engineering can be made mutually compatible in realistic hBN devices.