We use magnetic-field-dependent features in the photoluminescence of negatively charged nitrogen-vacancy centers to measure magnetic fields without the use of microwaves. In particular, we present a magnetometer based on the level anti-crossing in the triplet ground state at 102.4 mT with a demonstrated noise floor of 6 nT/Hz, limited by the intensity noise of the laser and the performance of the background-field power supply. The technique presented here can be useful in applications where the sensor is placed close to conductive materials, e.g., magnetic induction tomography or magnetic field mapping, and in remote-sensing applications since principally no electrical access is needed.

The negatively charged nitrogen-vacancy (NV) center in diamond has emerged as a unique nanoscale sensor full of interesting applications and has been extensively researched in the past years resulting in numerous technological breakthroughs. It forms the basis for sensors to detect magnetic fields1 (MFs), temperature,2,3 strain,4 rotation,5,6 electric fields,7 and quantum geometrical phases.8 In particular, the use of the optically detected magnetic resonance (ODMR) technique9 to probe the magnetically sensitive ground state of the NV center has proven to be a successful tool for sensitive MF measurements, both with single and ensembles of NV centers.1,10,11 Realization of ODMR sensing protocols typically involves green pump light for optical polarization of the NV centers, the application of microwave (MW) fields for the manipulation of their spin state, and an optical readout step involving either detection of NV-photoluminescence (PL)1 or absorption of 1042 nm radiation resonant with the singlet transition.12–15 The relevant energy levels can be seen in Fig. 1.1 When the applied MW fields are resonant with the splitting of the Zeeman sub-levels of the NV center, the transfer of spin populations results in an observable change in PL or absorption.

FIG. 1.

(a) Schematic of the experimental setup. (b) NV-center energy level schematic.

FIG. 1.

(a) Schematic of the experimental setup. (b) NV-center energy level schematic.

Close modal

There are cases where the use of strong MW fields proves to be detrimental to the sensing protocol and therefore can prohibit the employment of an NV-based sensor. An example is the detection of MFs generated by eddy currents in conductive materials in the context of magnetic induction tomography16,17 (MIT), a research application currently undertaken in our laboratory,18 where the presence of a conductive structure under examination will heavily affect the application of MW to the diamond. Another example is MF mapping of conductive, magnetic structures.19 

There have been several demonstrations of MW-free, and all-optical, diamond-based magnetic sensors, initially implemented with single NV centers attached to scanning atomic force microscopes,20–22 and more recently with ensembles of NV centers.19 These MW-free magnetometric protocols have been realized by exploiting either the properties of the NV centers' PL or their decoherence properties under the influence of external MFs. So far, these protocols remain either qualitative, requiring complicated setups to achieve high spatial resolution, or lack high magnetic-field sensitivities and bandwidth.

In this letter, we demonstrate the principles of a sensitive MW-free magnetometer by exploiting the properties of the ground-state level anti-crossing (GSLAC) of the NV center in diamond. We note that the presented technique can be extended to other magnetically sensitive features in the NV-PL or absorption, as discussed later, as well as features associated with other spin defects in solid-state systems. In particular, for the NV center, a ∼102.4 mT background MF causes ground-state Zeeman-sublevel degeneracy and mixing (the anti-crossing), which is visible as a drop in NV-PL under optical pumping. Any additional external MF will perturb the anti-crossing condition and, thus, result in a PL change that can be monitored and used for sensitive detection of the perturbing MF.

Using this technique, we demonstrate a MW-free magnetometer with a 6 nT/Hz MF sensitivity, a bandwidth exceeding 300 kHz, and a projected 0.43 nT/Hz sensitivity limited by the photon shot noise of the PL detection.23 

A schematic of the experimental setup is shown in Fig. 1. We used a single-crystal (111)-cut (2.1×2.3×0.6)mm3 diamond, synthesized using a high-temperature high-pressure (HPHT) method (Element 6). The diamond with an initial nitrogen concentration of <200 ppm was electron-irradiated at 14 MeV (dose: 1018cm2) and then annealed at 700 °C for 3 h. The resulting NV centers are randomly oriented along all four <111> crystallographic axes of the diamond.

The diamond is placed within a custom-made electromagnet (EM) [Fig. 1(a)] and can be rotated around the z-axis. The EM produces 11.2 mT/A; the current is provided by a computer-controlled power supply (Statron Typ 3257.1). The EM can be moved with a computer-controlled 3D translation stage (Thorlabs PT3-Z8) to align the MF with the respect to the diamond. Additionally, it can be rotated around the y-axis. Therefore, all degrees of freedom to place the diamond in the center of the magnet and to align the [111] NV-axis parallel to the MF can be addressed and optimized.

A secondary coil with four turns is wound around the diamond mount to apply small modulation of the MF. The additional oscillating component Bm is produced with the power amplified output (amplifier: AE Techron 7224-P) of a function generator (Tektronix AFG2021) that is also used as the local oscillator (LO) for the lock-in amplifier (LIA, SRS 865).

The NV centers in diamond are optically spin-polarized with 220 mW of 532 nm light taken from a 12 W laser (Coherent Verdi). The power is adjusted with a half-wave plate (HWP) and a polarizing beam-splitter (PBS) cube. Before the diamond, the light is sent trough an acousto-optical modulator (AOM, AA Electronics MT350-A0) to enable power modulation. Part of the laser light is split-off and measured on a photodiode (PD, Thorlabs PDA36A). The signal is input into a proportional–integral–derivative controller (PID, SRS SIM960) to stabilize the beam power. After the AOM, the beam is focused with a 40 mm focal-length lens into the diamond. The red/near-infrared NV-PL is collected with a 30 mm focal-length lens (numerical aperture 0.5). The collimated PL is separated from the green transmission with a dichroic mirror and a band-stop filter for 532 nm light before being focused onto a PD (Thorlabs PDA36A). The PD signal is sent to the LIA and demodulated at the LO frequency or measured at dc. After initial alignment and calibration of the magnet, the field was scanned from 0 mT to 120 mT in 5 s and the PL was monitored. Figure 2(a) shows the PL measured with the PD as a function of the MF. Figure 2(b) shows the corresponding LIA signal. The modulation frequency of the field was 100 kHz, the modulation depth ∼0.1 mT, and the LIA time constant 30 μs; 64 scans were averaged.

FIG. 2.

(a) NV-PL as a function of the applied MF normalized to the PL at 80 mT. (b) Derivative signal of (a) as given by the in-phase output (X) of the LIA. (c) Detail of the PL for fields around 51.4 mT showing additional features. (d) Detail of the PL for fields around the GSLAC at 102.4 mT. The data (red) weighted to ignore the side peaks are fitted with a Lorentzian function (green). The residuals of the fit are shown in black.

FIG. 2.

(a) NV-PL as a function of the applied MF normalized to the PL at 80 mT. (b) Derivative signal of (a) as given by the in-phase output (X) of the LIA. (c) Detail of the PL for fields around 51.4 mT showing additional features. (d) Detail of the PL for fields around the GSLAC at 102.4 mT. The data (red) weighted to ignore the side peaks are fitted with a Lorentzian function (green). The residuals of the fit are shown in black.

Close modal

Both plots contain several features previously discussed in the literature. The initial gradual decrease in PL is associated with a reduction in emission of the non-aligned NV centers due to spin-mixing.24 Around 51.4 mT, the observed features [1–7 in Fig. 2(c)] correspond to cross-relaxation events between the NV center and single substitutional nitrogen (P1) centers.24–26 The feature at 60 mT [9 in Fig. 2(c)] is attributed to cross relaxation with NV centers that are not aligned along the MF.24 At ∼102.4 mT [6 in Fig. 2(d)] is the GSLAC. Several additional features are visible. They can be associated with cross-relaxation events with either the nuclear spins of nearby P1 centers [2–5 and 7–10 in Fig. 2(d)]24–26 or nuclear spins of 13C atoms [1 in Fig. 2(d)].

The angles α and β between the NV-axis and the applied MF need to be precisely controlled [Fig. 1(a)] within ∼1 mrad.27 Misalignment causes a transverse field component which couples the mS=1 and the mS = 0 magnetic sublevels, broadens the observed GSLAC feature, and therefore leads to a reduction in magnetometric sensitivity. To optimize the GSLAC feature parameters, the angles and the position of the magnet were aligned until a minimum full width at half maximum of 1.2 mT and an optimal contrast of 4.5% was observed [Figs. 2(b) and 2(c)].

An important characteristic of any magnetometer is the sensitivity to ac MFs. For example, in eddy current sensing experiments, oscillating MFs need to be detected. To document the capacity of the MW-free magnetometer to detect these fields, the background MF is scanned around the GSLAC feature while a small oscillating MF (Bm0.09 mT) is applied at a given frequency. The frequency is then stepped from 300 Hz to 300 kHz (limited by the bandwidth of the power amplifier). The PL is measured with a PD and its oscillating component is read-out with the LIA. On the slopes of the GSLAC feature, the magnetometer is most sensitive to oscillating MFs. An example of such scan is shown in Fig. 3(a). For each frequency, the peak-to-peak response signal of the LIA amplitude output (R) is recorded and normalized with the oscillating current through the driving coil. This current is measured via its voltage drop over a 8 Ω resistor. Initially, the transfer function of the amplifier was measured, and during the experiment coarsely compensated for by adjusting the drive amplitude voltage. This way the current amplitude through the driving coil was effectively the same for all frequencies.

FIG. 3.

(a) Example of the LIA amplitude output (R) for a modulation frequency of 1 kHz. For noise reduction, a moving average was applied to the data. (b) Measurement of the frequency response of the magnetometer from 300 Hz to 300 kHz.

FIG. 3.

(a) Example of the LIA amplitude output (R) for a modulation frequency of 1 kHz. For noise reduction, a moving average was applied to the data. (b) Measurement of the frequency response of the magnetometer from 300 Hz to 300 kHz.

Close modal

Figure 3(b) shows the peak lock-in amplitude as a function of the modulation frequency. After an initial drop in the response to oscillating fields, the spectrum appears basically flat between 60 kHz and 300 kHz. The initial drop in the response and the slight increase for higher frequencies can be attributed to mutual inductance of the driving coil and the surrounding background field magnet as well as induction in the aluminum magnet mount. Observation of the induced current in the EM is consistent with the frequency characteristic of this feature.

Around the GSLAC feature, the derivative fluorescence signal as detected in the properly phased LIA X output depends linearly on the MF and can therefore be used for precise MF measurements. The calibration signal is shown in Fig. 4(a); the modulation frequency (100 kHz), modulation depth (∼0.1 mT), alignment, and laser power were optimized to maximize the slope, and therefore, the sensitivity of the magnetometer. The data near the zero-crossing are fitted with a straight line to translate the LIA output signal amplitude into MF. Then, the background MF is set to the center of the GSLAC feature (102.4 mT) where the magnetometer is maximally sensitive to external MFs. Figure 4(b) shows a time trace of the magnetometer response to a square-wave-modulated MF of ∼45 μT peak-to-peak amplitude applied with an additional external coil. The standard deviation of the data for a single step level is 1.8 μT (∼80 ms, 4800 samples), so that the steps in the MF can be observed with a signal-to-noise ratio of 25. For noise measurements, the LIA voltage output is recorded for 1 s and translated into MF variations. A fast-Fourier transform is performed to extract information of the MF noise. The data are displayed in Fig. 4(c) with the green pump-power stabilization (red) and without it (green). For comparison, similar data are collected at a MF of around 80 mT (blue). At this field, the setup is insensitive to MF variations and the data can be used to understand the technical noise level of the magnetometer. The noise floor is flat and around 6 nT/Hz. The electronic noise floor without green pump light and therefore without PL (black) is about 0.25 nT/Hz.

FIG. 4.

Magnetometer noise characterization. (a) Detail of the GSLAC feature around 102.4 mT, fitted linearly. (b) Response of the magnetometer to a square-wave MF modulation. (c) Noise of the magnetometer: magnetically sensitive and pump intensity stabilized (red), magnetically sensitive and pump not stabilized (green), insensitive to MFs with pump stabilized (blue) and electronic noise without the pump light (black).

FIG. 4.

Magnetometer noise characterization. (a) Detail of the GSLAC feature around 102.4 mT, fitted linearly. (b) Response of the magnetometer to a square-wave MF modulation. (c) Noise of the magnetometer: magnetically sensitive and pump intensity stabilized (red), magnetically sensitive and pump not stabilized (green), insensitive to MFs with pump stabilized (blue) and electronic noise without the pump light (black).

Close modal

Fundamentally, the magnetometer is limited by the shot-noise of the collected PL. For the given setup, the photon shot-noise limit is calculated to be 0.43 nT/Hz (Ref. 23) (resulting from 300 μW collected PL). However, this limit could be reduced by orders of magnitude by maximizing the amount of emitted and collected PL, possible by saturating the NV-PL and increasing the numerical aperture of the collection optics. The 1/f MF noise in Fig. 4(b) is attributed to the main field power supply. In an actual device, however, scanning of the MF would not be necessary, so that a small permanent magnet with less MF noise could be used. The frequency spikes at the line voltage frequency and its higher harmonics are also attributed to the power supply. They are the dominating noise component in Fig. 4(b). The roll-off for frequencies above 3.5 kHz is a result of filtering by the LIA. The time-constant for the measurements in Fig. 4(c) was 300 μs and the filter slope 24 dB/octave.

In conclusion, we demonstrated a MW-free, NV-center based magnetometer with a 6 nT/Hz noise floor and a bandwidth exceeding 300 kHz. This device can be useful in applications where MW spectroscopy cannot be performed on the NV centers. This is the case where the diamond-based sensor is placed in proximity to conductive objects, and as such, is of particular relevance for spatially resolved conductivity measurements in the context of MIT. The ability of the present technique to detect nuclear spins (seen as side-features near the GSLAC peak) with high signal-to-noise ratio indicates a possibility of applications in sensing spins external to the NV centers. If a layer of shallow-implanted NV-centers is used, spins external to diamond can be probed.

Future investigations will involve a thorough study of the lineshape and width of the signal near the GSLAC, as well as of the additional features around it, with the aim of understanding the fundamental sensitivity and bandwidth limitations28 of our sensing protocol. In addition, combination of the presented MW-free magnetometer with an absorption-based protocol will allow for MF sensing with a sensitivity exceeding the PL shot noise limit.14,15

We acknowledge support by the DFG through the DIP Program (FO 703/2-1). H.Z. is a recipient of a fellowship through GRK Symmetry Breaking (DFG/GRK 1581). N.L. acknowledges support from a Marie Curie International Incoming Fellowship within the 7th European Community Framework Programme. L.B. is supported by a Marie Curie Individual Fellowship within the second Horizon 2020 Work Programme. D.B. and A.J. acknowledge support from the AFOSR/DARPA QuASAR Program. V.M.A. acknowledges support from NSF Grant No. IIP-1549836. We thank P. R. Nelson, J. W. Blanchard, and D. Twitchen for useful discussions.

1.
L.
Rondin
,
J.-P.
Tetienne
,
T.
Hingant
,
J.-F.
Roch
,
P.
Maletinsky
, and
V.
Jacques
,
Rep. Prog. Phys.
77
,
056503
(
2014
).
2.
G.
Kucsko
,
P. C.
Maurer
,
N. Y.
Yao
,
M.
Kubo
,
H. J.
Noh
,
P. K.
Lo
,
H.
Park
, and
M. D.
Lukin
,
Nat. Phys. Lett.
500
,
54
(
2013
).
3.
V. M.
Acosta
,
E.
Bauch
,
M. P.
Ledbetter
,
A.
Waxman
,
L.-S.
Bouchard
, and
D.
Budker
,
Phys. Rev. Lett.
104
,
070801
(
2010
).
4.
P.
Ovartchaiyapong
,
K. W.
Lee
,
B. A.
Myers
,
A. C. B.
Jayich
,
A.
Stacey
,
D.
Budker
, and
L. C. L.
Hollenberg
,
Nat. Commun.
5
,
4429
(
2014
).
5.
M. P.
Ledbetter
,
K.
Jensen
,
R.
Fischer
,
A.
Jarmola
, and
D.
Budker
,
Phys. Rev. A
86
,
052116
(
2012
).
6.
A.
Ajoy
and
P.
Cappellaro
,
Phys. Rev. A
86
,
062104
(
2012
).
7.
F.
Dolde
,
H.
Fedder
,
M. W.
Doherty
,
T.
Nobauer
,
F.
Rempp
,
G.
Balasubramanian
,
T.
Wolf
,
F.
Reinhard
,
L. C. L.
Hollenberg
,
F.
Jelezko
, and
J.
Wrachtrup
,
Nat. Phys.
7
,
459
463
(
2011
).
8.
D.
Maclaurin
,
M. W.
Doherty
,
L. C. L.
Hollenberg
, and
A. M.
Martin
,
Phys. Rev. Lett.
108
,
240403
(
2012
).
9.
E.
van Oort
,
P.
Stroomer
, and
M.
Glasbeek
,
Phys. Rev. B
42
,
8605
(
1990
).
10.
T.
Wolf
,
P.
Neumann
,
K.
Nakamura
,
H.
Sumiya
,
T.
Ohshima
,
J.
Isoya
, and
J.
Wrachtrup
,
Phys. Rev. X
5
,
041001
(
2015
).
11.
J. M.
Taylor
,
P.
Cappellaro
,
L.
Childress
,
L.
Jiang
,
P.
Neumann
,
D.
Budker
,
P. R.
Hemmer
,
A.
Yacoby
,
R.
Walsworth
, and
M. D.
Lukin
,
Nat. Phys.
4
,
810
(
2008
).
12.
V. M.
Acosta
,
E.
Bauch
,
A.
Jarmola
,
L. J.
Zipp
,
M. P.
Ledbetter
, and
D.
Budker
,
Appl. Phys. Lett.
97
,
174104
(
2010
).
13.
V. M.
Acosta
,
A.
Jarmola
,
E.
Bauch
, and
D.
Budker
,
Phys. Rev. B
82
,
201202
(
2010
).
14.
Y.
Dumeige
,
M.
Chipaux
,
V.
Jacques
,
F.
Treussart
,
J.-F.
Roch
,
T.
Debuisschert
,
V. M.
Acosta
,
A.
Jarmola
,
K.
Jensen
,
P.
Kehayias
, and
D.
Budker
,
Phys. Rev. B
87
,
155202
(
2013
).
15.
K.
Jensen
,
N.
Leefer
,
A.
Jarmola
,
Y.
Dumeige
,
V. M.
Acosta
,
P.
Kehayias
,
B.
Patton
, and
D.
Budker
,
Phys. Rev. Lett.
112
,
160802
(
2014
).
16.
A.
Wickenbrock
,
S.
Jurgilas
,
A.
Dow
,
L.
Marmugi
, and
F.
Renzoni
,
Opt. Lett.
39
,
6367
(
2014
).
17.
C.
Deans
,
L.
Marmugi
,
S.
Hussain
, and
F.
Renzoni
,
Appl. Phys. Lett.
108
,
103503
(
2016
).
18.
A.
Wickenbrock
,
N.
Leefer
,
J. W.
Blanchard
, and
D.
Budker
,
Appl. Phys. Lett.
108
,
183507
(
2016
).
19.
D. A.
Simpson
,
J.-P.
Tetienne
,
J. M.
McCoey
,
K.
Ganesan
,
L. T.
Hall
,
S.
Petrou
,
R. E.
Scholten
, and
L. C. L.
Hollenberg
,
Sci. Rep.
6
,
22797 EP
(
2016
).
20.
P.
Maletinsky
,
S.
Hong
,
M. S.
Grinolds
,
B.
Hausmann
,
M. D.
Lukin
,
R. L.
Walsworth
,
M.
Loncar
, and
A.
Yacoby
,
Nat. Nano
7
,
320
(
2012
).
21.
L.
Rondin
,
J.-P.
Tetienne
,
P.
Spinicelli
,
C.
Dal Savio
,
K.
Karrai
,
G.
Dantelle
,
A.
Thiaville
,
S.
Rohart
,
J.-F.
Roch
, and
V.
Jacques
,
Appl. Phys. Lett.
100
,
153118
(
2012
).
22.
J.-P.
Tetienne
,
L.
Rondin
,
P.
Spinicelli
,
M.
Chipaux
,
T.
Debuisschert
,
J.-F.
Roch
, and
V.
Jacques
,
New J. Phys.
14
,
103033
(
2012
).
23.
V. M.
Acosta
,
K.
Jensen
,
C.
Santori
,
D.
Budker
, and
R. G.
Beausoleil
,
Phys. Rev. Lett.
110
,
213605
(
2013
).
24.
S.
Armstrong
,
L. J.
Rogers
,
R. L.
McMurtrie
, and
N. B.
Manson
,
Phys. Procedia
3
,
1569
(
2010
); in
Proceedings of the Tenth International Meeting on Hole Burning, Single Molecule and Related Spectroscopies: Science and Applications-HBSM
(2009).
25.
L. T.
Hall
,
P.
Kehayias
,
D. A.
Simpson
,
A.
Jarmola
,
A.
Stacey
,
D.
Budker
, and
L. C. L.
Hollenberg
,
Nat. Commun.
7
,
10211
(
2016
).
26.
J. D. A.
Wood
,
D. A.
Broadway
,
H. L. T.
,
A.
Stacey
,
D. A.
Simpson
,
J.-P.
Tetienne
, and
L. C. L.
Hollenberg
, e-print arXiv:1604.00160.
27.
R. J.
Epstein
,
F. M.
Mendoza
,
Y. K.
Kato
, and
D. D.
Awschalom
,
Nat. Phys.
1
,
94
(
2005
).
28.
C. S.
Shin
,
C. E.
Avalos
,
M. C.
Butler
,
D. R.
Trease
,
S. J.
Seltzer
,
J.
Peter Mustonen
,
D. J.
Kennedy
,
V. M.
Acosta
,
D.
Budker
,
A.
Pines
, and
V. S.
Bajaj
,
J. Appl. Phys.
112
,
124519
(
2012
).