---
title: Boron Vacancy (VB-) in hBN
url: https://www.emergentmind.com/topics/boron-vacancy-vb
type: topic
---

# Boron Vacancy (VB-) in hBN

The negatively charged boron vacancy, \(V_\text{B}^-\), is a point defect in hexagonal boron nitride (hBN) formed by a missing boron atom in the \(-1\) charge state. It is an optically addressable spin defect with a spin-triplet ground state, room-temperature optically detected magnetic resonance (ODMR), and a broad photoluminescence band centered near \(800\) nm, and it has emerged as a platform for quantum sensing at nanometer-scale or atomic-scale proximity in van der Waals materials [2111.10855] [2207.08357] [2305.05866].

## 1. Defect structure, symmetry, and spin Hamiltonians

In a single hBN layer, each boron atom is bonded to three \(sp^2\)-hybridized nitrogen atoms. The \(V_\text{B}^-\) defect is the absence of one boron atom, leaving three nitrogen dangling bonds equally spaced around the vacancy. In the symmetry-based description used for the ground state, the local point group is \(D_{3h}\), with the principal \(z\)-axis of the fine-structure tensor aligned with the crystallographic \(c\)-axis [2302.08831].

The ground-state electronic spin is \(S=1\). A standard effective Hamiltonian used across ODMR and ESR studies is
\[
H = hD\bigl(S_z^2 - S(S+1)/3\bigr) + hE\bigl(S_x^2 - S_y^2\bigr) + \mu_B g\,\mathbf{B}\cdot\mathbf{S},
\]
with \(D\) the axial zero-field splitting, \(E\) the transverse splitting, and \(g\) the electronic \(g\)-factor [2111.10855]. In related formulations the same physics is written as
\[
H = D\,S_z^2 + E(S_x^2-S_y^2)+g\mu_B\mathbf{B}\cdot\mathbf{S},
\]
or, when hyperfine structure is included explicitly, as a full spin Hamiltonian containing Zeeman, fine-structure, hyperfine, nuclear Zeeman, and nuclear quadrupole terms [2302.08831] [2307.06774].

Reported ground-state parameters are consistent at the level of a few percent but vary somewhat with method, sample, and temperature. Room-temperature ODMR work reported \(D_\mathrm{gs}=3.48\pm0.02\) GHz, \(E_\mathrm{gs}=48.0\pm7.1\) MHz, and \(g_\mathrm{gs}\approx2.00\) [2111.10855]. A cryogenic/room-temperature excited-state ODMR study reported \(D_\mathrm{GS}=3460\) MHz at \(293\) K and \(3684\) MHz at \(7\) K, with \(E_\mathrm{GS}=49\) MHz at room temperature and \(g_\mathrm{GS}\simeq(2.00\pm0.02)I\) [2112.06420]. Conventional ESR and high-frequency ENDOR reported \(D=3.600\pm0.010\) GHz and \(g_\mathrm{iso}=2.001(1)\) [2302.08831]. In ensemble studies the transverse splitting is typically much smaller than \(D\), often \(<100\) MHz, and is sensitive to rhombic distortion, strain, and local charge disorder [2507.16786] [2501.18481].

The charge state is a defining part of the defect’s functionality. The boron monovacancy can adopt several charge states, but only the \(-1\) state is optically active and spin-triplet, making \(V_\text{B}^-\) the relevant spin-photon interface for ODMR-based control and readout [2305.05866].

## 2. Hyperfine structure, quadrupole couplings, and spin-density localization

For the ground state coupled to the three nearest nitrogen nuclei, a symmetry-aware Hamiltonian is
\[
\hat H = \mu_B\,\mathbf{B}\cdot g\cdot \hat S + D\left[\hat S_z^2-\hat S(\hat S+1)/3\right]
+\sum_{k=1}^3 \hat S\cdot A^{(k)}\cdot \hat I^{(k)}
-g_n\mu_n\,\mathbf{B}\cdot\sum_k \hat I^{(k)}
+\sum_{k=1}^3 \hat I^{(k)}\cdot P^{(k)}\cdot \hat I^{(k)}.
\]
Here the three \(A^{(k)}\) tensors are related by \(120^\circ\) rotations about the \(c\)-axis, and the local principal axes of the hyperfine and quadrupole tensors follow the dangling-bond directions of the three nitrogen atoms [2302.08831].

Experimentally, the nearest-neighbor hyperfine interaction is axially symmetric to good approximation. ESR and \(94\) GHz ENDOR gave \(A_\perp=45.5\pm0.9\) MHz and \(A_\parallel=87.0\pm0.5\) MHz for each nearest \(^{14}\)N, while the nuclear quadrupole interaction is characterized by \(C_q=1.96\pm0.05\) MHz and \(\eta=0.070\pm0.005\) [2302.08831]. A linear-combination-of-atomic-orbitals analysis based on these values yielded the result that approximately \(84\%\) of the total \(V_\text{B}^-\) electron-spin density resides on the three nearest nitrogen atoms, directly demonstrating strong in-plane localization within a single BN layer [2302.08831].

Isotopic control modifies both line multiplicity and linewidth. In natural \(^{14}\)N material, the splitting per nitrogen is \(\mathcal A_{zz}(^{14}\mathrm N)=44.3(2)\) MHz, giving the familiar seven-line ESR pattern. In \(^{15}\)N-purified \(h^{10}\mathrm B^{15}\mathrm N\), the longitudinal coupling becomes \(|\mathcal A_{zz}(^{15}\mathrm N)|=64.1(2)\) MHz and only four ESR lines appear. The unresolved broadening from second-neighbor boron nuclei also depends on isotope: the ESR full width at half maximum changes from \(\Delta(^{11}\mathrm B)=52.9(3)\) MHz to \(\Delta(^{10}\mathrm B)=44.3(3)\) MHz, so \(^{10}\)B purification narrows the linewidth by about \(8.6\) MHz [2307.06774].

The nuclear-spin degrees of freedom are not restricted to nearest neighbors. ENDOR measurements identified \(^{14}\)N spins in the third nearest-neighbor shell, denoted N(3), at a distance of approximately \(0.4\) nm from the vacancy. For these sites the measured principal hyperfine values are \(A_{xx}=4.41\) MHz, \(A_{yy}=4.43\) MHz, and \(A_{zz}=6.86\) MHz; the corresponding isotropic and dipolar parts are \(A_\mathrm{iso}\simeq5.23\) MHz and \(A_\mathrm{dd}\simeq0.82\) MHz [2504.06798]. This establishes \(V_\text{B}^-\) as a probe of more remote nuclear moments in the host lattice.

Optical pumping can also polarize the host nuclei. In \(h^{10}\mathrm B^{15}\mathrm N\), optically induced \(^{15}\)N polarization was demonstrated near the ground-state level anticrossing. A four-level rate model gave
\[
\mathcal P_{\max}=\frac{\gamma_+-\gamma_-}{\gamma_++\gamma_-}
=\frac{1-|A_-/A_+|}{1+|A_-/A_+|},
\]
and the experimental power-dependent fit at \(B_z=92\) mT yielded \(\mathcal P_{\max}=22.6(7)\%\) and \(P_\mathrm{sat}=1.1(2)\) mW [2307.06774].

## 3. Excited-state spectroscopy and the symmetry question

A defining development in the \(V_\text{B}^-\) literature has been the direct observation of spin resonance in the orbital excited state. Continuous-wave ODMR revealed an additional set of resonances distinct from the ground-state transitions, and pulsed ODMR then verified their excited-state origin by eliminating any microwave-driven spin rotation after the system had fully relaxed into the ground state [2111.10855] [2112.02912] [2112.06420].

At room temperature, one study determined the excited-state parameters to be \(D_\mathrm{es}=2.11\pm0.03\) GHz, \(E_\mathrm{es}=74\pm42\) MHz, and \(g_\mathrm{es}\approx2.00\) [2111.10855]. A related room-temperature study reported \(D_\mathrm{es}=2.06\pm0.04\) GHz, \(E_\mathrm{es}=93.1\pm23.3\) MHz, and \(g_\mathrm{es}\simeq2.04\), and also resolved two types of excited-state hyperfine splitting, approximately \(90\) MHz from the three equivalent \(^{14}\)N nuclei and \(18.8\) MHz from a more distant \(^{11}\)B nucleus [2112.02912]. At cryogenic temperature, a separate ODMR study reported \(D_\mathrm{ES}=(2160\pm50)\) MHz at \(7\) K, \(D_\mathrm{ES}=2117\) MHz at \(293\) K, and an excited-state ODMR contrast of \(12\%\) at cryogenic temperature [2112.06420].

Level anticrossings provide a second route to excited-state characterization. In photoluminescence-versus-field measurements without microwaves, dips appear near the excited-state level anticrossing at about \(750\) G or \(769\) G, and near the ground-state anticrossing at about \(1240\) G, \(1307\) G, or \(1330\) G, depending on the dataset and fitting procedure [2111.10855] [2112.06420] [2112.02912]. These anticrossings are associated with field-induced mixing of \(\lvert m_S=0\rangle\) and \(\lvert m_S=-1\rangle\) under the spin Hamiltonian.

The symmetry of the excited state has been a central issue. One theoretical possibility is \(D_{3h}\) symmetry with a \({}^3E''\) orbital doublet; another is lowered \(C_{2v}\) symmetry with a \({}^3B_1\) orbital singlet. The experimental observation of a single excited-state spin-triplet resonance from \(10\) to \(300\) K, with no loss of contrast at low temperature, was interpreted as consistent with an orbital-singlet excited state. In that analysis, the robust low-temperature ODMR and the strong temperature dependence of the transverse splitting \(E_\mathrm{es}\) pointed to lowered \(C_{2v}\) symmetry and a \({}^3B_1\) orbital-singlet excited state [2111.10855]. By contrast, other excited-state treatments use a triplet \({}^3E\) description as an effective spectroscopic model [2112.06420]. The literature therefore contains both notations.

Temperature dependence strengthens the distinction between axial and transverse excited-state terms. From \(10\) to \(300\) K, \(D_\mathrm{es}\) remains essentially constant at \(2.1\) GHz within experimental error, whereas \(E_\mathrm{es}\) changes strongly, from \(258\pm2\) MHz at \(10\) K to approximately \(74\) MHz at \(300\) K [2111.10855]. This behavior underlies proposals for intrinsic thermometry and strain sensing using the excited-state manifold.

## 4. Creation pathways, charge-state control, and defect engineering

Several routes have been used to create \(V_\text{B}^-\) in hBN, and the outcome depends strongly on ion species, dose, flake thickness, cleaning, and substrate environment. In exfoliated hBN, focused-ion-beam implantation with \(12\) keV \( \mathrm N^+\) or \( \mathrm O^+\) ions reproducibly generated the \(800\) nm \(V_\text{B}^-\) photoluminescence band and ODMR near \(D\approx3.47\) GHz, whereas commercial \(40\) keV implantation did not yield the \(800\) nm spectral signature under the tested conditions, and electron irradiation up to \(300\) keV did not efficiently create \(V_\text{B}^-\) in the tested flakes [2305.05241].

Thickness and cleaning are important in exfoliated samples. For similar nitrogen fluence \(\Phi_N\simeq5\)–\(6\times10^{14}\,\mathrm{cm^{-2}}\), a \(150\) nm region gave signal-to-noise ratio \(\approx20\) whereas a \(100\) nm region gave \(\approx9\). Forming-gas annealing at \(375^\circ\)C for \(4\) h, or hot-plate annealing in air at \(400^\circ\)C for \(15\) min, restored \(V_\text{B}^-\) photoluminescence in as-exfoliated flakes, and the authors recommended keeping fluence below approximately \(10^{16}\,\mathrm{cm^{-2}}\) to avoid amorphization or sputtering [2305.05241].

In suspended, large-area CVD-grown hBN, the dominant created defect type depends on projectile mass. Light ions and neutrons favor boron vacancies, while \(30\) keV electrons predominantly generate a \(650\) nm emitter assigned to \(N_BV_N\). For \(30\) keV bombardment of suspended \(120\) nm hBN, the defect-yield proxy \(Y\equiv \mathrm{PL}_{800}/E_{2g}\) peaked at about \(0.7\) for He\(^+\) at \(\Phi\approx100\) ions nm\(^{-2}\), about \(0.8\) for Ne\(^+\) at \(\Phi\approx10\) ions nm\(^{-2}\), and about \(0.3\) for Ar\(^+\) at \(\Phi\approx1\) ion nm\(^{-2}\) [2503.18894]. Supported hBN behaves differently because ions penetrating through thin hBN can generate secondary particles in the substrate, which then promote other defect species; in that work, boron-vacancy purity depended strongly on hBN thickness relative to the end-of-range of the projectile [2503.18894].

Focused He-ion creation has also been quantified at the level of charge-state yield. Using \(30\) keV \(^4\)He\(^+\) irradiation of circular \(400\) nm patches in \(\sim80\) nm hBN on graphite, followed by ODMR-based extraction of local electric-field-induced splitting, a lower bound of \(0.2\%\) was obtained for the fraction of all vacancies that end up in the optically active, negatively charged state \(V_\text{B}^-\) [2501.18481]. This separates total vacancy creation from creation of the specific charge state relevant for spin optics.

The charge state can also be actively manipulated. In layered heterostructures under simultaneous laser and electron-beam excitation, deterministic and reversible switching between \(V_\text{B}^-\) and \(V_\text{B}^0\) was demonstrated through
\[
V_B^- \rightleftharpoons V_B^0 + e^-.
\]
The ionization rate \(k_\mathrm{ion}\) scales linearly with local electron-beam flux, the recovery rate \(k_\mathrm{rec}\) scales approximately linearly with \(532\) nm laser power, and the steady-state negative-state fraction follows
\[
\frac{N(V_B^-)}{N_\mathrm{tot}}=\frac{k_\mathrm{rec}(P)}{k_\mathrm{rec}(P)+k_\mathrm{ion}(f)}.
\]
Under \(-6\) V bottom-gate bias, quenching under \(f=7\times10^{15}\,\mathrm{e^-\,cm^{-2}\,s^{-1}}\) was \(30\%\), compared with \(57\%\) at \(+6\) V, showing that excess electrons stabilize the \(-1\) charge state [2305.05866].

## 5. Coherence, decoherence, and relaxation dynamics

Broadband microwave control has enabled direct measurements of \(T_2^*\) and \(T_2\) in isotopically enriched hBN. In a \( \lesssim30\) nm thick \(h^{10}\mathrm B^{15}\mathrm N\) flake stamped onto a \(1.5\,\mu\)m-wide gold stripline and irradiated with \(30\) keV He\(^+\), the \(m_s=0\leftrightarrow -1\) transition near \(2.49\) GHz at \(34.5\) mT displayed four overlapping hyperfine lines from the three nearest \(^{15}\)N nuclei. Sub-GHz Rabi driving gave \(f_R=421\pm1\) MHz at moderate power, corresponding to a \(\pi\)-pulse of about \(1.2\) ns, and coherent oscillations at \(1.26\) GHz were observed at three times that power [2604.27843].

Using those pulses, Ramsey interference and Hahn echo gave
\[
C_R(\tau)=\exp[-(\tau/T_2^*)^2],\qquad
C_E(2\tau)=\exp[-(2\tau/T_2)^\alpha].
\]
The measured values were \(T_2^*=13.8\pm0.5\) ns with Gaussian-like decay, and \(T_2=108.7\pm5.5\) ns with stretch factor \(\alpha=1.25\pm0.11\) [2604.27843]. The same work attributed the short \(T_2^*\) primarily to strong inhomogeneous hyperfine fields from nearby \(^{10}\)B and \(^{15}\)N spins, and emphasized that sub-GHz driving reduces pulse errors that can otherwise distort echo-based coherence estimates.

Longitudinal relaxation has been mapped over a much wider field and temperature range. Measurements of \(T_1\) from \(15\) to \(250\) K and fields up to \(7\) T revealed three regimes. At low temperature and low field, \(T\lesssim50\) K and \(H_\mathrm{ext}\lesssim0.05\) T, the decay is stretched exponential with \(\beta<1\), consistent with spin-spin interactions and disorder; at \(T=30\) K and \(H=0.03\) T, \(\beta\approx0.7\). At fields above about \(2\) T and \(T\lesssim100\) K, the relaxation rate follows an approximately direct-process scaling
\[
\Gamma_\mathrm{direct}\propto T\,\omega_0^{\,n_1}\approx T(H_\mathrm{ext})^{1.6},
\]
with \(n_1\approx1.6\). At \(T\gtrsim125\) K, the relaxation becomes nearly field-independent and follows
\[
\Gamma_\mathrm{Raman}\propto A_2\,T^{\,n_2},\qquad n_2\approx2.
\]
A combined model used
\[
\Gamma(H_\mathrm{ext},T)=\Gamma^\mathrm{spin-spin}+\Gamma^\mathrm{spin-ph},
\]
with a spin-bath correlation time \(\tau_c\approx100\) ps [2507.16786].

Defect density also affects \(T_1\). In CVD-grown hBN irradiated with He\(^+\), all-optical measurements at room temperature gave \(T_1\approx15\,\mu\)s at low \(V_\text{B}^-\) density and less than \(1\,\mu\)s at high density. The transverse splitting \(E/h\) increased from about \(51\) MHz at low dose to about \(95\) MHz at the highest dose, indicating that dense irradiation also strengthens static disorder or local symmetry breaking [2503.18894].

First-principles many-body simulations further predict a field-driven crossover in Hahn-echo decoherence. Below a transition boundary \(B_T\), decoherence is dominated by independent nuclear-spin dynamics and occurs within submicrosecond timescales; above \(B_T\), pairwise flip-flop processes dominate and the coherence time reaches tens of microseconds. The predicted transition boundaries are \(5020\) G for \(h\)-\(^{10}\)B\(^{14}\)N and \(2050\) G for \(h\)-\(^{11}\)B\(^{14}\)N, in excellent agreement with the numerical results in that study [2505.03306]. This suggests that magnetic field and isotopic composition jointly determine whether the dense nuclear environment acts primarily through single-spin modulation or pairwise dynamics.

## 6. Photophysics, emission engineering, heterostructures, and sensing modalities

The intrinsic optical signature of \(V_\text{B}^-\) is a broad photoluminescence band centered near \(800\) nm, with strong electron-phonon coupling and no easily resolved room-temperature zero-phonon line in uncoupled ensemble measurements [2305.05866] [2202.10980]. Coupling to a high-\(Q\) Si\(_3\)N\(_4\) nanobeam cavity revealed the room-temperature zero-phonon line at \(773\pm2\) nm and showed a maximum cavity-enhanced ZPL intensity increase of \(21\times\) at the field antinode, consistent with Purcell-enhanced redistribution of the emission [2202.10980].

The defect’s weak brightness has motivated several enhancement strategies. In a low-loss plasmonic nano-patch antenna, the raw intensity enhancement averaged over the laser spot was about \(250\times\), while the enhancement localized to the nano-patch region was calculated to be about \(1685\times\). The ODMR contrast was preserved, increasing from about \(3\%\) on bare SiO\(_2\)/Si to about \(6\%\) in the antenna, with linewidth \(\Delta\nu\sim100\) MHz unchanged, and a continuous-wave ODMR sensitivity of about \(138\) nT/\(\sqrt{\mathrm{Hz}}\) was estimated at \(30\,\mu\)W laser power [2207.08357].

A van der Waals sensitizer strategy has also been demonstrated. In hBN–PbI\(_2\) heterostructures, \(V_\text{B}^-\) photoluminescence was enhanced by \(5\)–\(45\times\), depending on temperature and excitation, while ODMR contrast remained roughly preserved. Under \(488\) nm excitation the PL count rate increased by about \(4\times\), under \(555\) nm excitation by about \(3\times\), and the reported magnetic sensitivity improved by about \(190\%\) and \(325\%\), respectively, with a best-case \(\eta\sim76\,\mu\mathrm T\,\mathrm{Hz}^{-1/2}\) at \(160\,\mu\)W [2602.02256]. That work attributed the enhancement to type-I band alignment, efficient fluorescence resonance energy transfer, and a DFT-predicted photon-ratcheting mechanism with minimal hybridization of the spin defect states.

Mechanical deformation can produce a different kind of enhancement. Suspended regions of hBN exhibited photoluminescence up to \(30\times\) brighter than substrate-supported regions, while ODMR contrast, linewidth, and spin lifetime were well preserved. In the bright suspended regions, the measured \(T_1\) was \(10.3(3)\,\mu\)s, compared with \(9(1)\,\mu\)s in darker supported regions. The increase in brightness correlated with larger local deformation, increased zero-field splitting \(2E/h\), appearance of the deformation-activated \(E'\) Raman peak near \(1290\) cm\(^{-1}\), and shifts in Kelvin-probe surface potential [2509.00912].

At the same time, proximity to other 2D materials does not necessarily quench the optical readout. A quantitative FRET study found that for hBN sensing layers thicker than approximately \(3\) nm, photoluminescence quenching by monolayer graphene is negligible, with \(Q_\mathrm{tot}<1.1\), meaning \(\Gamma_\mathrm{FRET}\ll\Gamma_\mathrm{nr}\). Below about \(3\) nm, the FRET rate increases rapidly, and at \(d=2\) nm a quenching factor \(Q_\mathrm{tot}\approx3\) was observed. By simultaneously fitting thickness-dependent photoluminescence and lifetime data, that study extracted an intrinsic radiative rate
\[
\Gamma_0=(1.35\pm0.68)\times10^5\,\mathrm{s^{-1}},
\]
corresponding to an intrinsic radiative lifetime \(\tau_0\simeq7.4\,\mu\)s, and concluded that \(V_\text{B}^-\) is unusually tolerant of integration into ultrathin van der Waals heterostructures because of its intrinsically low quantum yield [2512.03970]. The same work directly verified that \(2\) nm hBN on MoS\(_2\) showed no lifetime shortening, whereas MoTe\(_2\) did quench the defect photoluminescence, consistent with band-gap-controlled acceptor behavior [2512.03970].

These optical and spin properties motivate multiple sensing modes. High ODMR contrast in both ground and excited states, strong temperature dependence of \(E_\mathrm{es}\), and sensitivity of \(E\) to in-plane strain support magnetic, thermal, and strain sensing [2111.10855]. At high magnetic fields up to \(7\) T, the ground-state splitting reaches approximately \(200\) GHz, which has been proposed as a route to sub-terahertz sensing through field-dependent \(T_1\) relaxometry [2507.16786]. ENDOR measurements of remote nitrogen moments further indicate that \(V_\text{B}^-\) can serve as a probe of host-lattice nuclear spins in van der Waals materials [2504.06798]. The two-dimensional host geometry, together with preserved spin functionality near surfaces and in heterostructures, is the recurring basis for its use as a surface-native quantum sensor [2207.08357] [2512.03970].

Source: https://www.emergentmind.com/topics/boron-vacancy-vb