Scanning NV Center Spectroscopy

• Scanning NV spectroscopy is a technique that employs near-surface NV centers in diamond as local quantum sensors to detect magnetic fields, nuclear spins, and charge dynamics.
• It utilizes a range of protocols—including continuous-wave ODMR, correlation, and Fourier spectroscopy—to extract detailed spin Hamiltonians and local environmental parameters.
• Practical implementations feature engineered probe geometries such as nanopillars and shallow-implanted diamond chips, enabling nanometer-scale resolution and multiparametric analysis.

Scanning nitrogen-vacancy center spectroscopy denotes a family of spatially resolved spectroscopic methods in which a near-surface or tip-mounted nitrogen-vacancy (NV) center in diamond serves as a local quantum sensor while the probe is positioned with nanometric control over a target system. In the cited literature, this umbrella includes continuous-wave and pulsed optically detected magnetic resonance (ODMR), free-precession Fourier spectroscopy, correlation spectroscopy, double-electron–electron-resonance (DEER), amplitude-encoded high-field NMR protocols, spectral hole burning, level anti-crossing spectroscopy, and multidimensional coherent spectroscopy. Depending on geometry, the measured spectra encode local stray fields, nuclear Larmor frequencies, hyperfine tensors, scalar couplings, spin-bath composition, charge-state dynamics, or optical linewidths, with spatial resolution set primarily by the NV–sample separation rather than by diffraction (Welter et al., 2022, Xu et al., 6 Mar 2025, Boss et al., 2015).

1. Spin Hamiltonians and sensing principle

The spectroscopic role of the NV center is anchored in its spin-1 ground state and its anisotropic coupling to external and internal degrees of freedom. In the NV frame $(\bar x,\bar y,\bar z)$ with $\bar z \parallel \langle 111\rangle$, the ground-state Hamiltonian is written as

$$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$

with $D = 2\pi\cdot 2.87\ \mathrm{GHz}$, $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$, and $E$ a small transverse splitting. For $B_\perp \ll D$, the transition frequencies obey

$f_{\pm 1}=D\pm \gamma_e B_\parallel,$

while a transverse field produces the quadratic correction

$$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$

These relations determine both the usable bias-field geometry and the spectroscopic observable extracted from ODMR scans (Welter et al., 2022).

When the target is a nearby nuclear spin or spin cluster, the relevant Hamiltonian acquires conditional hyperfine terms. For the NV–$^{13}\mathrm{C}$–$\bar z \parallel \langle 111\rangle$0 complex, the cited formulation is

$\bar z \parallel \langle 111\rangle$1

with $\bar z \parallel \langle 111\rangle$2 and $\bar z \parallel \langle 111\rangle$3. For a single $\bar z \parallel \langle 111\rangle$4 nucleus under a static bias field, the free-precession Hamiltonian is

$\bar z \parallel \langle 111\rangle$5

where $\bar z \parallel \langle 111\rangle$6. These Hamiltonians underlie the central idea of scanning NV spectroscopy: the local environment shifts the NV or the coupled nuclear subsystem in a way that can be read out spectroscopically through fluorescence after tailored microwave and optical control (Laraoui et al., 2013, Boss et al., 2015).

2. Probe architectures, crystal orientation, and distance scales

Scanning implementations use either a single shallow NV embedded near the apex of a diamond nanopillar or an ensemble layer of NVs beneath a sample. A near-surface diamond chip implanted with $\bar z \parallel \langle 111\rangle$7 ions at $\bar z \parallel \langle 111\rangle$8 yields an approximately $\bar z \parallel \langle 111\rangle$9-deep NV layer, while scanning-probe variants mount a diamond nanopillar containing a near-surface NV on an AFM cantilever for raster scans (Boss et al., 2015). Correlation spectroscopy has been explicitly extended to a scanning geometry by proposing a shallow-implanted NV at less than $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$0 depth on the tip of a diamond nanopillar or AFM cantilever (Laraoui et al., 2013). At larger scale, high-field AERIS employs an ensemble of $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$1 NV centers a few microns below the surface under a picoliter sample in a $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$2 field (Munuera-Javaloy et al., 2022).

Crystal cut fixes the accessible bias-field orientation. Standard $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$3-cut probes place the four $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$4 axes at $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$5 to the surface normal, which makes a purely in-plane bias field impossible without a large transverse component. By contrast, $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$6-cut diamond has two $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$7 axes exactly in the surface plane $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$8 and two perpendicular to it $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$9. Selecting an in-plane NV orientation allows application of a bias field purely in the sample plane without any transverse component, preserving ODMR contrast in large in-plane fields (Welter et al., 2022).

A practical description of spatial resolution requires three vertical distances (Xu et al., 6 Mar 2025):

Quantity	Definition	Representative values
$D = 2\pi\cdot 2.87\ \mathrm{GHz}$0	Mechanical stand-off between the lowest point of the diamond tip and the highest point of the sample surface	AM-AFM: $D = 2\pi\cdot 2.87\ \mathrm{GHz}$1; FM-AFM: $D = 2\pi\cdot 2.87\ \mathrm{GHz}$2
$D = 2\pi\cdot 2.87\ \mathrm{GHz}$3	Magnetic stand-off from the NV spin to the top of the magnetic sample	AM-AFM: $D = 2\pi\cdot 2.87\ \mathrm{GHz}$4, median $D = 2\pi\cdot 2.87\ \mathrm{GHz}$5; FM-AFM: $D = 2\pi\cdot 2.87\ \mathrm{GHz}$6, median $D = 2\pi\cdot 2.87\ \mathrm{GHz}$7
$D = 2\pi\cdot 2.87\ \mathrm{GHz}$8	Sub-surface NV depth below the diamond surface	$D = 2\pi\cdot 2.87\ \mathrm{GHz}$9, median $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$0; example $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$1

These distances are not interchangeable. The cited work shows that the stand-off distance is mainly limited by features on the surface of the diamond tip, not solely by implantation depth, and that frequency-modulated AFM feedback yields systematically lower and more consistent magnetic stand-offs than amplitude-modulated feedback. The same study reports a minimum NV-to-sample distance of $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$2 from $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$3 NMR on a polytetrafluoroethylene surface in soft contact (Xu et al., 6 Mar 2025).

3. Nuclear-spin spectroscopy protocols

Correlation spectroscopy extracts long-time nuclear dynamics by separating two sensing blocks with a variable hold interval. In the cited implementation, two consecutive Hahn echoes on the NV are interleaved by a storage or hold interval, optionally with phase storage in the host $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$4 nucleus. The first block accumulates a phase $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$5, the second a phase $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$6, and the measured signal is proportional to $\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$7. The correlation function is

$\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$8

and its Fourier transform yields resonances at

$\gamma_e = 2\pi\cdot 28\ \mathrm{GHz/T}$9

A central peak appears at $E$0 for distant $E$1, with side peaks shifted by $E$2 up to several hundred kHz for first-shell carbons. The spectral resolution is set by the NV $E$3, or by the $E$4 $E$5 when phase is stored in the $E$6, giving $E$7. In the scanning extension, the same two-block protocol is proposed to provide sub-$E$8 spatial resolution and sub-kHz spectral resolution (Laraoui et al., 2013).

Free-precession Fourier spectroscopy addresses the same target from a different angle. The sequence $E$9 prepares nuclear coherence, lets it evolve during a variable free-precession interval, and maps it back to the NV population. The signal has the form

$B_\perp \ll D$0

Because the protocol measures nuclear free precession directly, the Fourier spectrum displays a single peak at the true resonance, with no spurious harmonics. At $B_\perp \ll D$1, reported values include $B_\perp \ll D$2 at $B_\perp \ll D$3, $B_\perp \ll D$4 at $B_\perp \ll D$5, and $B_\perp \ll D$6 at $B_\perp \ll D$7. The same framework supports two-dimensional Fourier spectroscopy through two evolution intervals, correlating $B_\perp \ll D$8 peaks with $B_\perp \ll D$9 peaks to associate $f_{\pm 1}=D\pm \gamma_e B_\parallel,$0 and $f_{\pm 1}=D\pm \gamma_e B_\parallel,$1 of the same nucleus. Reported precisions include $f_{\pm 1}=D\pm \gamma_e B_\parallel,$2, $f_{\pm 1}=D\pm \gamma_e B_\parallel,$3, and $f_{\pm 1}=D\pm \gamma_e B_\parallel,$4, with scanning adaptation aimed at hyperspectral images of local NMR spectra at approximately $f_{\pm 1}=D\pm \gamma_e B_\parallel,$5 spatial resolution (Boss et al., 2015).

At large magnetic fields, the AERIS protocol shifts otherwise inaccessible high-frequency NMR information into the NV detection band. An RF $f_{\pm 1}=D\pm \gamma_e B_\parallel,$6 pulse triggers nuclear precession, the nuclei then evolve during a free-precession interval $f_{\pm 1}=D\pm \gamma_e B_\parallel,$7 while the NV is idle, and during an induced-rotation interval $f_{\pm 1}=D\pm \gamma_e B_\parallel,$8 a continuous RF drive and an NV XY4 sequence map the amplitude of the nuclear field onto the sensor. The stroboscopic record

$f_{\pm 1}=D\pm \gamma_e B_\parallel,$9

reveals chemical shifts and $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$0-splittings after discrete Fourier transform. Since the NV does not participate during $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$1, the ultimate linewidth is set by the nuclear $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$2, not the NV $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$3: $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$4 For ethanol at $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$5, $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$6, $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$7, and $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$8, simulations give Lorentzian peaks with full-width at half-maximum $$f_{\pm 1}\approx D\pm \gamma_e B_\parallel + \frac{3\gamma_e^2 B_\perp^2}{2D}.$$9, resolving $^{13}\mathrm{C}$0 chemical shifts and the $^{13}\mathrm{C}$1 $^{13}\mathrm{C}$2-coupling. The same paper states that the protocol can be grafted onto scanning confocal or scanning-probe NV microscopes, where a single NV at approximately $^{13}\mathrm{C}$3 stand-off enables sub-$^{13}\mathrm{C}$4 chemical NMR imaging (Munuera-Javaloy et al., 2022).

4. Electron-spin, defect, and spin-bath spectroscopy

DEER turns the NV into a local electron-spin spectrometer by embedding a bath-spin $^{13}\mathrm{C}$5-pulse into an NV Hahn echo. In the single-NV implementation, the normalized signal is

$^{13}\mathrm{C}$6

with

$^{13}\mathrm{C}$7

Sweeping the second microwave frequency $^{13}\mathrm{C}$8 yields a local ESR spectrum of the bath. In a type-Ib crystal at $^{13}\mathrm{C}$9, five dips at approximately $\bar z \parallel \langle 111\rangle$00, $\bar z \parallel \langle 111\rangle$01, $\bar z \parallel \langle 111\rangle$02, $\bar z \parallel \langle 111\rangle$03, and $\bar z \parallel \langle 111\rangle$04 match the hyperfine-split transitions of substitutional nitrogen (P1). By comparing experiment with Monte Carlo simulation, the number of detected spins was estimated to be $\bar z \parallel \langle 111\rangle$05, and site-to-site variability among four NVs showed that each NV probes its own microscopic bath volume (Abeywardana et al., 2015).

For sub-micron, $\bar z \parallel \langle 111\rangle$06-oriented diamond layers, DEER has been extended to quantify both P1 and NVH$\bar z \parallel \langle 111\rangle$07 defects. The cited contrast model for a single bath species is

$\bar z \parallel \langle 111\rangle$08

and for two species the exponent generalizes to $\bar z \parallel \langle 111\rangle$09. High-resolution spectra acquired with a bath $\bar z \parallel \langle 111\rangle$10-pulse of approximately $\bar z \parallel \langle 111\rangle$11 reveal the four P1($\bar z \parallel \langle 111\rangle$12) lines, two forbidden transitions, and four narrow NVH$\bar z \parallel \langle 111\rangle$13 lines between $\bar z \parallel \langle 111\rangle$14. Simultaneous spectral fits gave, for sample A1, $\bar z \parallel \langle 111\rangle$15, $\bar z \parallel \langle 111\rangle$16, and $\bar z \parallel \langle 111\rangle$17, $\bar z \parallel \langle 111\rangle$18; for sample B1, $\bar z \parallel \langle 111\rangle$19 and $\bar z \parallel \langle 111\rangle$20. The same work reports ppm-level sensitivity to dark defects in sub-micron layers and frames the method as a fast feedback loop for CVD recipe optimization (Findler et al., 2023).

A complementary route is magnetic-field scanning through level anti-crossings. In that method, the applied field is modulated as

$\bar z \parallel \langle 111\rangle$21

with $\bar z \parallel \langle 111\rangle$22, and the lock-in output approximates $\bar z \parallel \langle 111\rangle$23. The main NV ground-state LAC appears near $\bar z \parallel \langle 111\rangle$24, with additional lines produced by coupling to P1, NV$\bar z \parallel \langle 111\rangle$25, NV$\bar z \parallel \langle 111\rangle$26, or other defects. Fits of the multispin Hamiltonian reproduced the seven lines around $\bar z \parallel \langle 111\rangle$27 using the known P1 hyperfine tensor and identified an additional center $\bar z \parallel \langle 111\rangle$28 with $\bar z \parallel \langle 111\rangle$29, $\bar z \parallel \langle 111\rangle$30, $\bar z \parallel \langle 111\rangle$31, $\bar z \parallel \langle 111\rangle$32, $\bar z \parallel \langle 111\rangle$33, and $\bar z \parallel \langle 111\rangle$34. This suggests that scanning-field spectroscopy can identify paramagnetic impurities through avoided-crossing structure even when ODMR alone does not isolate the defect Hamiltonian (Anishchik et al., 2016).

5. Resolution limits, line narrowing, and operating-field constraints

The literature distinguishes several different resolution ceilings. In free-precession Fourier spectroscopy, the ultimate spectral resolution is $\bar z \parallel \langle 111\rangle$35, where $\bar z \parallel \langle 111\rangle$36 is the longest usable free-precession interval, and the cited value $\bar z \parallel \langle 111\rangle$37 implies $\bar z \parallel \langle 111\rangle$38 (Boss et al., 2015). In correlation spectroscopy, the resolution is instead tied to the NV $\bar z \parallel \langle 111\rangle$39 or the $\bar z \parallel \langle 111\rangle$40 $\bar z \parallel \langle 111\rangle$41 when the host nucleus is used as a memory (Laraoui et al., 2013). In AERIS, the linewidth is ultimately set by the coherence of the nuclear spin signal itself, $\bar z \parallel \langle 111\rangle$42, so Hz-scale linewidths are compatible with NV-based detection at large field (Munuera-Javaloy et al., 2022). These distinctions matter because they show that high spectral resolution in scanning NV spectroscopy is not generically limited by NV $\bar z \parallel \langle 111\rangle$43.

Hole-burning spectroscopy narrows ODMR lines by saturating a selected subensemble. For the ground-state transition frequencies

$\bar z \parallel \langle 111\rangle$44

pump–probe correlations distinguish magnetic-field broadening from strain broadening. Continuous-wave hole burning modulates the pump at approximately $\bar z \parallel \langle 111\rangle$45 and isolates a narrow spectral hole in the probe transition; pulsed hole burning uses a narrow-band $\bar z \parallel \langle 111\rangle$46-pulse to transfer a selected spectral slice. Reported ordinary ODMR linewidths range from $\bar z \parallel \langle 111\rangle$47, while hole-burning reduces them to $\bar z \parallel \langle 111\rangle$48, consistent with removal of the $\bar z \parallel \langle 111\rangle$49 contribution of approximately $\bar z \parallel \langle 111\rangle$50. A $\bar z \parallel \langle 111\rangle$51 Fourier-limited hole was observed in pulsed mode, and approximately $\bar z \parallel \langle 111\rangle$52 is stated to be possible with longer pulses. The same pump–probe relation yields a magnetic-field-insensitive thermometry channel, with $\bar z \parallel \langle 111\rangle$53 and a temperature sensitivity on the order of $\bar z \parallel \langle 111\rangle$54 (Kehayias et al., 2014).

At optical frequencies, multidimensional coherent spectroscopy isolates homogeneous linewidth, phonon dephasing, and spectral diffusion of the NV zero-phonon line. Three $\bar z \parallel \langle 111\rangle$55 resonant pulses in box geometry generate a heterodyne-detected photon echo, and the cross-diagonal width of the resulting two-dimensional spectrum gives the homogeneous dephasing rate. The cited fit over $\bar z \parallel \langle 111\rangle$56 yields $\bar z \parallel \langle 111\rangle$57, $\bar z \parallel \langle 111\rangle$58, and $\bar z \parallel \langle 111\rangle$59, implying $\bar z \parallel \langle 111\rangle$60 at $\bar z \parallel \langle 111\rangle$61. A spectral-diffusion-induced slope $\bar z \parallel \langle 111\rangle$62 at one location and $\bar z \parallel \langle 111\rangle$63 at another corresponds to a characteristic diffusion time of order $\bar z \parallel \langle 111\rangle$64. The same work reports a temperature-dependent Stark splitting consistent with internal transverse electric fields of approximately $\bar z \parallel \langle 111\rangle$65 at $\bar z \parallel \langle 111\rangle$66 and $\bar z \parallel \langle 111\rangle$67 at $\bar z \parallel \langle 111\rangle$68 (Liu et al., 2020).

Field geometry imposes an additional constraint. Transverse fields mix NV levels and can quench ODMR contrast entirely, whereas properly aligned $\bar z \parallel \langle 111\rangle$69-cut probes retain more than $\bar z \parallel \langle 111\rangle$70 of their zero-field contrast up to $\bar z \parallel \langle 111\rangle$71; standard $\bar z \parallel \langle 111\rangle$72-cut probes lose more than $\bar z \parallel \langle 111\rangle$73 contrast above $\bar z \parallel \langle 111\rangle$74 and become unusable by $\bar z \parallel \langle 111\rangle$75 (Welter et al., 2022). Correlation spectroscopy also exhibits a field threshold: the amplitude of $\bar z \parallel \langle 111\rangle$76 is large and field-independent for $\bar z \parallel \langle 111\rangle$77, then collapses abruptly for $\bar z \parallel \langle 111\rangle$78, which the cited discussion attributes to a quantum-to-classical crossover in the spin-cluster dynamics (Laraoui et al., 2013).

6. Charge-state and singlet-manifold spectroscopy of the NV center

Scanning NV spectroscopy also includes spectroscopy of the NV center’s own optical and charge dynamics. Field-quenching photoluminescence spectroscopy addresses the $\bar z \parallel \langle 111\rangle$79 singlet manifold by comparing field-on and field-off spectra under selected excitation wavelengths, powers, and temperatures. In a heavily nitrogen-doped CVD sample, increased NV$\bar z \parallel \langle 111\rangle$80 photoluminescence and decreased NV$\bar z \parallel \langle 111\rangle$81 zero-phonon-line width were observed in the presence of an applied magnetic field, indicating ionization from the long-lived $\bar z \parallel \langle 111\rangle$82 singlet state. The measured single-photon ionization threshold from $\bar z \parallel \langle 111\rangle$83 brackets

$\bar z \parallel \langle 111\rangle$84

corresponding to $\bar z \parallel \langle 111\rangle$85, and

$\bar z \parallel \langle 111\rangle$86

The same work gives a temperature-dependent parametrization $\bar z \parallel \langle 111\rangle$87 with $\bar z \parallel \langle 111\rangle$88, and argues that direct knowledge of $\bar z \parallel \langle 111\rangle$89 underpins spin-to-charge-conversion readout schemes in scanning probes (Blakley et al., 2023).

A later temperature-resolved study refined this threshold with magnetically mediated spin-selective photoluminescence quenching. Measurements were performed for excitation wavelengths between $\bar z \parallel \langle 111\rangle$90 and $\bar z \parallel \langle 111\rangle$91 in $\bar z \parallel \langle 111\rangle$92 increments, and for temperatures from about $\bar z \parallel \langle 111\rangle$93 to $\bar z \parallel \langle 111\rangle$94 in $\bar z \parallel \langle 111\rangle$95 increments. The reported ionization energy lies between $\bar z \parallel \langle 111\rangle$96 and $\bar z \parallel \langle 111\rangle$97, with about a two-fold reduction in uncertainty, and no statistically significant shift within $\bar z \parallel \langle 111\rangle$98 over $\bar z \parallel \langle 111\rangle$99. The same protocol associates negative contrast with $$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$00 shelving below threshold and positive contrast with net NV$$\frac{\hat H}{\hbar} = D S_z^2 + \gamma_e(B_x S_x + B_y S_y + B_z S_z) + E(S_x^2-S_y^2),$$01 production above threshold (Ung et al., 14 Jul 2025).

Taken together, these optical studies show that scanning nitrogen-vacancy center spectroscopy is not restricted to detecting external magnetic signals. It also encompasses spectroscopic interrogation of the NV defect’s internal manifolds, charge conversion, and spin-selective ionization pathways. A plausible implication is that future scanning implementations will combine nanoscale magnetic sensing with charge-state-resolved optical spectroscopy in the same probe, especially where spin-to-charge conversion, photoelectric detection, or single-defect chemical contrast are required.