Kilohertz electron paramagnetic resonance spectroscopy of single nitrogen centers at zero magnetic field

See allHide authors and affiliations

Science Advances  27 May 2020:
Vol. 6, no. 22, eaaz8244
DOI: 10.1126/sciadv.aaz8244


Electron paramagnetic resonance (EPR) spectroscopy is among the most important analytical tools in physics, chemistry, and biology. The emergence of nitrogen-vacancy (NV) centers in diamond, serving as an atomic-sized magnetometer, has promoted this technique to single-spin level, even under ambient conditions. Despite the enormous progress in spatial resolution, the current megahertz spectral resolution is still insufficient to resolve key heterogeneous molecular information. A major challenge is the short coherence times of the sample electron spins. Here, we address this challenge by using a magnetic noise–insensitive transition between states of different symmetry. We demonstrate a 27-fold narrower spectrum of single substitutional nitrogen (P1) centers in diamond with a linewidth of several kilohertz, and then some weak couplings can be resolved. Those results show both spatial and spectral advances of NV center–based EPR and provide a route toward analytical (EPR) spectroscopy at the single-molecule level.


Electron paramagnetic resonance (EPR) spectroscopy, a technique for studying paramagnetic targets, is an indispensable component of magnetic resonance spectroscopy for investigations of molecular structures and fast dynamics (1). An important goal of this technique is to extract precise information from small-volume samples (2), which requires both high spatial and high spectral resolution. In the past decades, intensive efforts have been devoted to promote the spatial resolution, and ultimately, single-spin EPR has been realized by various approaches, such as magnetic resonance force microscopy (3), scanning tunneling microscopy (4), and nitrogen-vacancy (NV) center–based EPR spectroscopy (57). Among them, the NV center is more promising for biological applications because of the compatibility with ambient conditions (57). However, the current spectral resolution of those techniques is on the orders of megahertz (312). It is insufficient to resolve molecules of slightly different structures (13) or local polarity profiles (14). For example, the distance between paramagnetic centers studied by EPR techniques is in the range of 1.8 to 6 nm (15), corresponding to the dipolar coupling strength of megahertz to submegahertz, the variations of which induced by conformation changes will be even smaller.

The line broadening is generally attributed to the limited spin-state lifetime of sensors and the decoherence of target spins. The former one has recently been overcome by using quantum memories (16, 17). The latter, which is more fundamental (18), arises from the magnetic couplings to bath spins. To address spin decoherence, a simple but powerful strategy is the use of particular spin states that are naturally insensitive to external perturbations. This phenomenon exists in various physical systems, for example, so-called “clock transitions” in trapped ions (19, 20) and phosphorus donors in silicon (21, 22), or transitions between long-lived states in nuclear spin resonance spectroscopy (18, 23). The line-narrowing phenomenon has also been found long ago in conventional zero-field EPR spectroscopy (2426). Nonetheless, the detection via inductive pickup suffers from the very small thermal polarization at zero field and thus requires a large amount of samples (27), which prevents practical applications of zero-field EPR. Fortunately, the NV center is a good magnetic sensor at zero magnetic field (28). In a previous study (12), we have found that the NV center is a promising sensor for zero-field EPR without loss of sensitivity, because the statistical fluctuations of the spin polarization, which do not depend on the magnetic field, dominate in its nanoscale detection volumes rather than thermal polarization.

Here, we show that the magnetic noise–insensitive transition can be observed by a single NV center. To investigate the spectral figure, we develop a correlation method for zero-field EPR spectroscopy. We demonstrate the high-resolution nature on single substitutional nitrogen centers (P1 centers) in diamond, which are electron spins with coupling to 15N nuclear spins. The linewidth of the observed spectra is as narrow as 8.6 kHz, which is a 27-fold improvement compared with the ordinary spectrum. Our results show not only the spatial but also the spectral advances of NV centers in EPR detections.


Magnetic noise–insensitive transition at zero magnetic field

Our model consists of an optically probed NV center and a dark electron spin (Fig. 1A), which serve as the magnetic sensor and the target, respectively. In general, the target spin can be any spin-half electron spins (S = 1/2) with hyperfine coupling to a nearby spin half-integer nuclear spin (I = n/2). For simplicity, here we take the S = 1/2, I = 1/2 system as an example. At zero magnetic field, the spin Hamiltonian of this electron-nuclei system is determined solely by the hyperfine interaction and can be written as (24)H0=A(SxIx+SyIy)+ASzIz(1)where A and A are the hyperfine constants. The eigenstates consist of one antisymmetric singlet ∣S0⟩ with F = 0 and three symmetric triplet states ∣T0⟩ and ∣T±1⟩ with F = 1, where F = S + I is the total angular momentum. The corresponding eigenvalues areωS0=A2A4,ωT0=A2A4,ωT±1=A4(2)

Fig. 1 Schematic representation of the NV center–based zero-field EPR spectrometer.

(A) Geometry of the experimental setup. The sensor is a shallow NV center in diamond, which is observed by a confocal microscope with green laser excitation and red fluorescence collection. The Ω-shape waveguide radiates microwave and radio frequency (RF) to control the sensor and the target. The three-dimensional Helmholtz coils (only one of them is represented) are used to compensate for residual magnetic fields. The inset gives our model, where the target consists of a spin-1/2 electron spin and a spin-1/2 nuclear spin. (B) Energy levels of the target spin. The degeneracy of ∣T±1⟩ is lifted by a magnetic field B0 with a linearly dependent splitting, while ∣S0⟩ and ∣T0⟩ have zero first-order dependence on B0. The ST±1 and ST0 transitions can be driven by perpendicular RF1 and parallel RF2 pulses, respectively (section S1). Up/down arrows denote the high-field spin-up/down states.

In the presence of magnetic noise δb, the energy levels of the target spin will fluctuate, leading to line broadening. It can be described by a perturbation to the HamiltonianδH=j=x,y,zδjSj(3)where δj = γe · δbj, with γe being the gyromagnetic ratio of the electron spin. Here, we ignore the Zeeman terms of the nuclear spin because of the nearly three orders of magnitude smaller gyromagnetic ratio than electron spins. According to the perturbation theory, the energy level shifts can be simplified as (section S1)δωS0δx2+δy22(A+A)δz24A,δωT0δx2+δy22(AA)+δz24A,δωT±1±δz2(4)

It shows that the ∣S0⟩ and ∣T0⟩ states have zero first-order dependence on magnetic field, and the frequency fluctuation of the ∣S0⟩ → ∣ T0⟩ (denoted as ST0 hereinafter) transition is reduced to ∼δ2/A (Fig. 1B). Therefore, a line-narrowing phenomenon will appear.

Correlation detection protocol for zero-field EPR spectroscopy

Such a narrowed ST0 spectrum is challenging to observe by the previous microwave power-sweeping method (12), as it requires extreme power stability of the entire microwave circuits. The commonly used double electron-electron resonance (DEER) method (510) is also not suitable, because the interrogation time is limited by the coherence time T2 of the NV center. Inspired by the correlation spectroscopy of nuclear spins (16), we develop a modified correlation detection protocol for zero-field EPR spectroscopy. Then, the sensor’s lifetime can be increased to the spin-locking relaxation time T, which is usually much longer than T2 for shallow NV centers (12, 29).

As shown in Fig. 2A, our pulse protocol consists of two zero-field DEER sequences (see Materials and Methods) separated by a spin-locking sequence. After initialization of the NV center in ∣0⟩ state by a 532-nm laser excitation, a resonant microwave π pulse creates a superposition state (1+1)/2, which evolves during the first DEER period with an accumulated phase φ1, and thus becomes (eiφ11+eiφ11)/2. We rewrite it as cos φ1 ∣ ψ+⟩ + isin φ1 ∣ ψ⟩, where ψ±=(1±1)/2. During the spin-locking period, both the ∣ψ+⟩ and ∣ψ⟩ states can be locked by the continuous driving field, but the coherence between them vanishes, which can be regarded as an ensemble projection measurement (see Materials and Methods). In this period, we can perform any manipulations on the target spin, which determines the accumulated phase φ2 during the second DEER period. If the NV state is projected to ∣ψ+⟩ (or ∣ψ⟩) with probability of cos2φ1 (or sin2φ1) after the spin-locking period, the final state will be cos φ2 ∣ 0⟩ + i sin φ2 ∣ ψ⟩ (or sin φ2 ∣ 0⟩ − icos φ2 ∣ ψ⟩) with the population of ∣0⟩ read out by the photoluminescence (PL) rate of the NV center. The resulting correlation signal is given byScorr=12[1+cos 2φ1cos 2φ2](5)where the brackets denote statistical average.

Fig. 2 Correlation protocol for EPR measurements.

(A) Schematics of the pulse sequence. The interrogations are two DEER sequences (n = 2) for detections of the target spin state. The detected phase signal is stored on the populations of NV dressed states and protected by spin locking. The correlation signal depends on the manipulations on the target spin, which is denoted by the black dashed box. MW, microwave. (B and C) Rabi oscillations for the ST±1 (B) and ST0 (C) transitions. Insets give the corresponding manipulations on P1 centers. The spin-locking time is fixed to 10 μs. The red lines are sine fittings. Error bars indicate SEM. arb., arbitrary.

Experimental demonstration of the high-resolution EPR spectroscopy

We perform the experiments on a coupled system of NV and P1 centers at zero magnetic field. We compensate for the residual magnetic field to ∼0.01 G by three-dimensional Helmholtz coils (see Materials and Methods and fig. S1). P1 center is another kind of defect in diamond, consisting of only a substitutional nitrogen atom (11,12). The 15N P1 center is an S = 1/2, I = 1/2 system, making it an ideal candidate for demonstrations of the high-resolution zero-field EPR spectroscopy. Both the NV and P1 centers are created by implantation of 15N2+ ions into a bulk diamond. With proper energy and dose (see Materials and Methods), the created NV centers can be well resolved by a home-built confocal microscope. For some NV centers, one can find a single adjacent P1 center by performing the zero-field DEER measurements (see section S2 and fig. S2). The hyperfine constants of 15N P1 center are A = 114 MHz and A = 160 MHz (30) and thus the expected transition frequencies of ST0 and ST±1 transitions are 114 and 137 MHz, respectively.

The P1 center can be fully controlled by resonant radio frequency (RF) pulses and reliably read out by the correlation detection protocol. Considering that the orientation of P1 center jumps between the four kinds of N─C bonds (fig. S1B) due to the Jahn-Teller effect (31), we adjust the RF direction to the vertical direction of the diamond surface by moving the P1 center to the central area of the Ω-shape waveguide, so that the effective control-field strengths for the different oriented P1 centers are the same to simplify the control. Figure 2B shows the Rabi oscillation between ∣S0⟩ and ∣T±1⟩, by varying the RF pulse length during the spin-locking period. The spin-locking power, denoted by the corresponding Rabi frequency of the NV center, is 30 MHz, resulting in T ∼ 150 μs. As a comparison, T2 is just 16 μs (fig. S3). This interrogation time can be further improved to T1 (∼milliseconds), although with compromise of signal contrast (see Materials and Methods). As described above, the ST0 transition is insensitive to the magnetic field, and thus also insensitive to the coupling with the NV center. To observe this mute transition, we use ∣T±1⟩ as auxiliary states. As shown in Fig. 2C, we can also observe the Rabi oscillation between ∣S0⟩ and ∣T0⟩. The detailed calculations of these correlation Rabi measurements are in section S3.

To obtain the EPR spectra, we perform the Ramsey experiments on the P1 center. Specifically, for the ST±1 transition, the manipulations during spin-locking period are two resonant π/2 RF pulses with varying separations t, as illustrated in Fig. 3A. Here, the RF pulses remain unchanged during the variation of t to avoid any spurious effect on the NV spin. Thus, the phase difference between them varies proportionally to t, resulting in oscillation signals (Fig. 3A). The Fourier transformation gives the resonance spectrum. The decay of the signal remarks the decoherence process, resulting in line broadening in the frequency domain. Similarly, the Ramsey experiment of the ST0 transition can also be performed with the assistance of the ST±1 transition, which shows a much slower decay. The detailed calculations of these correlation Ramsey measurements are in section S4. Figure 3B gives the Fourier transformation spectra and clearly shows the improvement of the spectral resolution. The Gaussian fittings give the linewidths [quoted as the full width at half maximum (FWHM)] of the spectra, the minimums of which are 230 ± 20 kHz and 8.6 ± 0.4 kHz for the ST±1 and ST0 transitions, respectively. The latter has (27 ± 3)–fold improvement. If the line broadening is solely induced by magnetic noises, the estimated improvement should be >130 (see Materials and Methods). This deviation suggests that other decoherence sources emerge when the magnetic noise is suppressed, such as electric or strain field noises (32).

Fig. 3 High-resolution EPR spectroscopy of single P1 centers.

(A) Correlation signals of Ramsey experiments for the ST±1 (top) and ST0 (bottom) transitions. Insets give the corresponding manipulations on P1 centers. Each experiment consists of two measurements, π/2 → − π/2 and π/2 → π/2, and the differential signal is presented. The top data have a 10× magnification in horizontal axis and still show faster decay. All the signals are undersampled, and the actual frequencies can be recovered with the prior knowledge of the rough resonance frequencies. (B) Fourier transformations of the time-domain data. The frequencies are relative to 137 MHz (ST±1) and 114 MHz (ST0). The points are experimental results, while the solid lines are two-Gaussian fittings. The fitting FWHM is 230±20 kHz (left) and 260±20 kHz (right) for the ST±1 spectra and 11.6±0.6 kHz (left) and 8.6±0.4 kHz (right) for the ST0 spectra.

Both of the spectra show clear line splitting, but the reasons are different. For the ST±1 transition, the line splitting of 370 kHz is induced by the coupling with a nearby 13C nuclear spin, which is widespread in natural-isotope diamonds. As a comparison, we repeat the measurements on a 12C isotopically purified diamond, where the line splitting disappears but remains in the ST0 spectrum (fig. S4). Actually, the 13C coupling will not induce obvious splitting for the ST0 transition because of the quadratic dependence on the magnetic field according to Eq. 4. A probable reason for the ST0 line splitting is the existence of local electric or strain fields (33). Similar with donor electron spins in silicon (21, 34), the hyperfine coupling of P1 centers should also depend on electric or strain fields. Given that the hyperfine coupling here is anisotropic with jumping orientations, up to four ST0 lines can emerge in a local static electric or strain field. The observation of only two lines suggests that the field is along a symmetric direction, probably perpendicular to the surface, as reported in a recent paper (35). Different from the magnetic field, the electric or strain field induces the same-order frequency shifts of the ST0 and ST±1 transitions because both of them are linearly dependent on the hyperfine constants. Therefore, such weak couplings cannot be resolved by the ST±1 spectrum or previous non–zero-field measurements.

To further investigate the spectral figure, we apply a small magnetic field B and observe the spectral variance. As shown in Fig. 4A, further splittings appear with increasing B, which suggests that the ST±1 line splitting at zero field comes from the coupling with a two-level system (i.e., 13C nuclear spin) rather than a local static magnetic field. On the other hand, the line shape of the ST0 transition remains unchanged, despite the overall small shifts and line broadening (Fig. 4B). The line splitting is almost constant versus B but varies for different P1 centers (fig. S4), which suggests different local electric or strain environments (33). The resonance frequencies clearly show the linear and quadratic dependence on B for the ST±1 and ST0 transitions, respectively, which is well predicted by Eq. 4. The almost changeless peak positions in Fig. 4B also reveal the line-narrowing nature, which is degrading with increasing B because of the growing fluctuation of (γeB + δ)2.

Fig. 4 Magnetic-field dependence of the line shape.

The magnetic field B is applied along the vertical direction of the diamond surface, which is a symmetric direction of the four P1 axes. All the points are experimental data, fitted with m-Gaussian (m = 2 − 4) functions (solid lines). Insets give the fitting peak positions (fres) versus magnetic fields with error bars defined by the fitting errors. (A) ST±1 spectrum. Each peak at zero field splits independently with increasing B, and the expected splitting (solid lines in the inset) is calculated according to Eq. 4. (B) ST0 spectrum. The two-peak pattern remains unchanged with increasing B, despite overall peak shifts and line broadening. Each group of peak positions (solid lines in the inset) is calculated according to Eq. 4.


In conclusion, we have presented a method for measuring EPR spectra with high spectral resolution by using a magnetic noise–insensitive singlet-triplet transition at zero magnetic field. The experiment demonstrates that the EPR linewidth of P1 centers could be reduced by more than one order of amplitude to several kilohertz, and then some weak nonmagnetic effects can be found. Our results show that, as a sensor for EPR spectroscopy, the NV center has the potential to simultaneously achieve high spatial and high spectral resolution.

The P1 centers in diamond are ideal systems for demonstrations of this technique, while practical applications require target spins external to the diamond surface. Although this line-narrowing phenomenon also exists in many other paramagnetic materials, such as 15N nitroxide spin labels (7) and some transition ions (25, 26), the observations with NV centers are more challenging. First, fast diffusions of the target spins will cancel their dipolar couplings with the NV sensor in the absence of external magnetic field, making the detections impossible. Thus, the sample should be solid or quasi-solid. Fortunately, the zero-field powder spectra do not have the inhomogeneous line-broadening issue (12, 2426). Second, the spectral improvement will degrade with increasing magnetic noise strength as ∼δ/A and will vanish when the noise strength δ is comparable to or even larger than the hyperfine constant A. It means that the coherence properties of the target spins should not be too worse relative to 1/A. This may be a problem for the transition ions, whose coherence times are usually much worse than free-radical spin labels. Nevertheless, it is recently reported that, with proper synthetic design, the coherence times of some transition metal complexes can even approach milliseconds (36). Third, the current measurement is still time-consuming. The widely used nitroxide spin labels may be bleached during the strong laser illumination (7, 8). A possible solution is improving the detection efficiency by implanting the spin-to-charge readout technique (37) or extending the lifetimes of spin labels by lowering the temperature and moving to a vacuum environment (8).

If successfully generalizing to external paramagnetic targets, this technique will benefit many EPR studies with substantially improved precision. For example, zero-field EPR spectroscopy is a powerful tool to unambiguously extract the hyperfine constants (12), which can reflect the local polarity profiles (14). By improving spectral resolution, more detailed microenvironment information can be explored. In addition, such magnetic field–insensitive singlet and triplet states can serve as a narrowband filter, which filters out magnetic fields of all frequencies except the one resonant with the transition. If considering dual spin targets in the presence of spin bath, all the target-environment spin dipolar couplings are killed, but the target-target coupling survives. Thus, the latter can be precisely detected by deploying our zero-field EPR technique, which offers interesting avenues for enhanced accuracy and detection range in distance measurement of spin-modified biomolecules (15).


Experimental setup and diamond samples

The optical part of our setup is a home-built confocal microscopy, with a diode laser (CNI Laser MGL-III-532) used for illumination and an avalanche photodiode (Perkin Elmer SPCM-AQRH-14) used for photon collection. All the microwave and RF pulses are generated by an arbitrary waveform generator (Keysight M8190a) and amplified by two amplifiers accordingly (Mini-Circuits ZHL-16W-43+ and LZY-22+). The Ω-shape waveguide has a diameter of 50 μm with a radiation field perpendicular to the plane in the central area.

Two diamond samples are involved, denoted as diamonds A and B. Both of them are obtained commercially, 100-oriented, electronic grade, and implanted with 15N2+ ions. Diamond A is implanted with an energy of 30 keV and a dose of 3 × 1010 cm−2. We fabricated nanopillars on diamond A to enhance the photon-collection efficiency (38). Diamond B has an extra growth layer with 99.9% 12C isotopic purity. It is implanted with an energy of 5 keV and a dose of 1.7 × 1010 cm−2. The simulated ion straggling (∼7 nm for diamond A and ∼4 nm for diamond B) is much smaller than the mean spacing of ions (∼60 nm for diamond A and ∼80 nm for diamond B), so the detected signal is dominated by a single P1 center arising from the adjacent nitrogen ion. The data in fig. S4 are measured on diamond B, and all the other data are measured on diamond A.

Compensation of residual magnetic field

The NV center can detect magnetic fields via optically detected magnetic resonance (ODMR) spectrum. In a small magnetic field (γeBD), the line splitting is proportional to the component parallel to the N-V axis. However, the smallest resolvable magnetic field is limited by the ODMR linewidth. To reliably compensate for the residual magnetic field, we take advantage of the symmetric nature of the ODMR spectrum versus magnetic fields. As shown in fig. S1A, with increasing currents applied to the coils, the ODMR linewidth first decreases and then increases. The symmetric center point marks the zero axial magnetic field. By using three differently oriented NV centers, the residual vector magnetic field B can be totally compensated. Figure S1B gives the schematics of the three NV centers, among which NV A/B and C are oriented perpendicular to the X and Y axis, respectively. We perform the compensation as follows: (i) Compensation of Bz. The ODMR linewidths of NV A and B are measured with sweeping Z-coil currents, and two symmetric center points IZ, A and IZ, B are recorded (fig. S1C). Then, the Z-coil compensation current is IZ = (IZ, A + IZ, B)/2 owing to the symmetry of NV A and B with respect to the Z axis. (ii) Compensation of By. Since the NV A and B are insensitive to Bx, the compensation of By is straightforward. As shown in fig. S1D, the symmetric center point is just the Y-coil compensation current. (iii) Compensation of Bx. Similar with (ii), the symmetric center point of NV C is just the X-coil compensation current (fig. S1E). The maximum fit error of the symmetric center point is 1 mA, while the stability of the power supplies is ∼4 mA. The magnetic field linearly depends on the current with a coefficient of 2.8 G/A. Therefore, the magnetic field after compensation is estimated to be ∼0.01 G.

Zero-field DEER

For NV center–based EPR, the DEER sequence is widely used for the detection of electron spins (57, 9, 10), where the NV center itself is also an electron spin (S = 1). Because of the degeneracy of ∣ ± ⟩1 at zero magnetic field, the zero-field DEER sequence is slightly different (fig. S2). Specifically, a π pulse can flip the NV state from ∣0⟩ to a superposition state (1+1)/2, which serves as an interferometer, and a 2π pulse can switch ∣1⟩ and ∣− 1⟩, which is a decoupling operation. During the evolution, the superposition state accumulates a phase ϕ = a · τ, where a is the dipolar coupling strength, and becomes (eiϕ+1+eiϕ1)/2. Last, a second π pulse reverses this state back to ∣0⟩ with the population cos2ϕ read out by the PL rate.

The dipole-dipole coupling between the sensor and the target can be approached byHddC(nNV,ntar,r)SzNVSzzT(6)where the coupling strength C depends on the orientations of the NV center nNV, the target spin ntar, and the spacing vector r between them. SzzT is a reduced spin operator of the target spin in the singlet-triplet basis {∣T+1⟩, ∣S0⟩, ∣T0⟩, ∣T−1⟩}SzzT=12(1001)(7)

The detailed derivation process is in section S1. Now, one can directly see that the coupling will disappear if the target spin is in the ∣S0⟩ or ∣T0⟩ states, with a magnetic field–insensitive nature. Hence, the DEER measurement cannot directly capture the signal of the ST0 transition. Otherwise, if the target spin is in the ∣T±1⟩ states, the coupling strength a = ± C/2. Similar to the NV center, a π RF pulse will flip the target spin from ∣S0⟩ to a superposition state of ∣T±1⟩, while a 2π RF pulse will switch ∣T+1⟩ and ∣T−1⟩. The latter leads to a stronger equivalent coupling strength (fig. S2), which is more favorable in the presence of decoherence. As calculated in section S2, the zero-field DEER signal isS(τ)=34+14e(τ/T2,NV)pcos Cτ(8)where the decoherence of the NV center is described by a stretched exponential decay, and p is in the range of 1 to 3 determined by the dynamic of bath (39).

Zero-field spin locking

The Hamiltonian of the NV center at zero magnetic field isH0=DSz2(9)where S(S = 1) is the NV electron spin operators and D is the zero-field splitting. During the spin-locking period, we apply a phase-modulated microwave of the form (40)H1=Ω1cos [Dt+2Ω2Ω1sin Ω1t]Sx(10)where Ω1 is the corresponding Rabi frequency and Ω2 is the phase modulation strength. In this experiment, Ω2 = 0.3Ω1. By moving to the interaction picture, the Hamiltonian becomesHI=eif(t)Sz2(H0+H1)eif(t)Sz2f(t)Sz2=Ω12Sx2Ω2cos Ω1tSz2(11)where f(t) = Dt + (2Ω21) sin Ω1t, and we ignore the high-frequency items. This is just the common Rabi model. Moving again to the second interaction picture and ignoring the high-frequency items, we can writeHII=eiΩ1t2SxHIeiΩ1t2SxΩ12Sx=Ω22(Sz2Sy2)(12)

For an initial state ∣ψ(0)⟩, the evolution of the NV state can be written asψ(t)=eif(t)Sz2eiΩ12SxteiHIItψ(0)(13)

By choosing an integral-period evolution time, i.e., t=k·4πΩ1,k=1,2,, the evolution is simplified toψ(t)=eiHIItψ(0)(14)where the phase item exp [iDtSz2] is absorbed in the normal rotating reference frame. After the first DEER period, the NV state becomes cos ϕ1 ∣ψ+⟩ + i sin ϕ1 ∣ ψ⟩. Note that both ∣ψ+⟩ and ∣ψ⟩ are eigenstates of HII, so the spin-locking process is just a free induction decay process in the interaction picture. The populations of ∣ψ±⟩ are protected while the coherence between them is destroyed by the environmental noise. Therefore, the NV state becomes a mixed state of ∣ψ+⟩ and ∣ψ⟩ after the spin-locking period.

Alternative protocol for correlation detection

Alternatively, we can also avoid the spin-locking microwave in some further practical applications, such as sensing in a living cell, where long microwave pulses may be potentially harmful to biological tissues. This correlation protocol consists of two DEER sequence separated by a free evolution (fig. S5). Similar to the analysis in Correlation detection protocol for zero-field EPR spectroscopy, the NV state starts from ∣0⟩ and becomes cos φ1 ∣0⟩ + i sin φ1 ∣ψ⟩ after the first DEER period. During the free evolution, the coherence vanishes with NV state projecting to either ∣0⟩ with a probability of cos2φ11/2sin2φ1, or a totally mixed state (∣0〉〈0∣ + ∣1〉〈1∣ + ∣ − 1〉〈−1∣)/3 with a probability of 3/2sin2φ1. The former has a similar evolution to the first DEER period and becomes cos φ2∣0⟩ + isin φ2∣ψ⟩ after the second DEER period, while the latter has an idle evolution. Therefore, the correlation signal isScorr(φ1,φ2)=(cos2φ112sin2φ1)cos2φ2+12sin2φ1=18[3+cos 2φ1+cos 2φ2+3cos 2φ1cos 2φ2](15)

Here, one can see that the signal contrast is reduced by 25%, which is the cost of removing the spin-locking microwave. Nevertheless, the limitation of interrogation duration can be further released to T1, NV.

Spectral linewidth analysis

According to Eq. 4, the fluctuations of the transition frequencies can be approached byδωST±1±12δz,δωST0AA2A2(δx2+δy2)+12Aδz2(16)

For Ramsey measurements, the magnetic noise is dominated by the low-frequency component, which can be modeled as a quasi-static random variable characterized by a normal distributionP(δ=x)=12πσex22σ2(17)with SD σ used to characterize the noise amplitude. Accordingly, the energy levels will fluctuate and can also be modeled as some kind of distributions, of which the SDs can be used to estimate the linewidth. According to Eq. 16, the SD of ωST±1 can be directly written asσST±1=12σ(18)

Considering isotropic magnetic noise, i.e., σx = σy = σz = σ, the SD of ωST0 can be calculated asσST0=(δωST0)2δωST02=σ24A2(A2A2)2+12A2(19)where 〈•〉 = ∫ ∫ ∫ • Px)Py)Pz)dδxdδydδz. Thus, the spectral-resolution improvement is estimated to beχσST±1σST0=[64A2(A2A2)2+8A2]12·(σST±1)1(20)

For the data in Fig. 3, the linewidth of the ST±1 spectrum is 230 kHz, which is defined by the FWHM. Thus, σST±1=FWHM/8ln 298 kHz, and χ ∼ 130. Actually, this improvement is underestimated because the distribution of the ST0 transition is not Gaussian but a sharper type (fig. S6). The quadratic dependence on magnetic field is responsible to the asymmetric line shape. Both the disappearance of the asymmetric pattern and the degrading spectral-resolution improvement of the experimentally measured ST0 spectrum suggest the existence of other decoherence resources. The magnetic noise can be increased by applying noise currents to the Helmholtz coils, and then the asymmetric pattern appears (fig. S6).


Supplementary material for this article is available at

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank R.-B. Liu for helpful discussions. Funding: This work was supported by the National Key Research and Development Program of China (grant nos. 2018YFA0306600 and 2016YFA0502400), the National Natural Science Foundation of China (grant nos. 81788101, 91636217, 11722544, 11761131011, and 31971156), the CAS (grants nos. GJJSTD20170001, QYZDY-SSW-SLH004, and YIPA2015370), the Anhui Initiative in Quantum Information Technologies (grant no. AHY050000), the CEBioM, the National Youth Talent Support Program, the Fundamental Research Funds for the Central Universities, the China Postdoctoral Science Foundation (grant no. BX20180294), and USTC Research Funds of the Double First-Class Initiative (grant no. YD2340002004). Author contributions: J.D. and F.S. supervised the entire project. F.K., F.S., and J.D. conceived the experiments. P.Z. and F.K. prepared the setup and performed the experiments, together with Z.Q. and Z.W. P.Y. and M.W. performed the ion implantation and nanofabrication. F.K. and Z.H. carried out the calculations. F.K., P.Z., and F.S. wrote the manuscript. All authors analyzed the data, discussed the results, and commented on the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article