Temperature dependence of the nitrogen-vacancy magnetic resonance in diamond V. M. Acosta,1, ∗ E. Bauch,1, 2 M. P. Ledbetter,1 A. Waxman,3 L.-S. Bouchard,4 and D. Budker1, 5, y 1 Department of Physics, University of California, Berkeley, CA 94720-7300 2 Technische Universität Berlin, Hardenbergstraße 28, 10623 Berlin, Germany 3Department of Physics, Ben-Gurion University, Be’er-Sheva, 84105, Israel 4Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095 5Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley CA 94720, USA (Dated: October 25, 2018) The temperature dependence of the magnetic resonance spectra of nitrogen-vacancy (NV-) ensembles in the range of 280-330 K was studied. Four samples prepared under different conditions were analyzed with NV- concentrations ranging from 10 ppb to 15 ppm. For all samples, the axial zero-field splitting (ZFS) parameter, D, was found to vary significantly with temperature, T , as dD=dT = −74:2(7) kHz=K. The transverse ZFS parameter, E, was non-zero (between 4 and 11 MHz) in all samples, and exhibited a temperature dependence of dE=(EdT ) = −1:4(3) × 10−4 K−1. The results might be accounted for by considering local thermal expansion. The temperature dependence of the ZFS parameters presents a significant challenge for diamond magnetometers and may ultimately limit their bandwidth and sensitivity. - 3 3 Magnetometers based on nitrogen-vacancy (NV) ing the NV centers’ A2 ! E optical transition ensembles in diamond [1–3] promise high-sensitivity, via a phonon sideband [19]. The same lens was rivaling those of superconducting quantum interfer- used to collect fluorescence from the diamond which ence devices (SQUIDs) [4] and alkali vapor magne- was then passed through a dichroic mirror and a tometers [5], in a scaleable solid state system that 650-800 nm bandpass filter and detected with a pho- can be operated over a wide range of temperatures. todiode. Noise due to laser power fluctuations was This remarkable combination of spatial resolution reduced by normalizing the fluorescence signal to [6, 7] and magnetic sensitivity [8] make diamond a reference photodiode which monitored the inci- magnetometers promising candidates for remote- dent laser power. The output of a microwave signal detection and low-field nuclear magnetic resonance generator was amplified, passed through a straight spectroscopy [9–12], nano-scale biological imaging ∼ 200-µm diameter copper wire of length ∼ 5 mm [6, 13–15], and studies of novel magnetic and super- placed within 500 µm of the focused light beam, conducting materials [3, 16]. Until now, the temper- and terminated with 50-Ω impedance. For temper- ature dependence of the magnetic resonance spectra ature control, the diamond was thermally connected has not been systematically studied and has only to a copper heatsink and placed inside an insulated briefly been mentioned in the literature [3, 17]. In aluminum housing. The temperatures of the heat this Letter, we report a striking temperature depen- sink and housing were controlled with separate ther- dence of the magnetic-resonance spectra of NV- en- moelectric (TE) elements. Unless otherwise stated, sembles in diamond over the temperature range of the results reported in this Letter were obtained < 280-330 K. These findings have important implica- with a magnetic field of ∼ 1 G, laser-light power of tions for the design of diamond magnetometers and ∼ 150 mW, and microwave power (after the wire) arXiv:0911.3938v4 [cond-mat.mes-hall] 17 Feb 2010 may ultimately limit their sensitivity. of ∼ 10 dBm. The resonance spectra were recorded using the For temperature scans, the temperature of the continuous-wave Fluorescence Detected Magnetic copper plate in direct thermal contact with the di- Resonance (FDMR) method [17, 18]. Light from amond was monitored with an AD590 sensor. The a 514-nm Argon-ion laser was focused with a 2:5 cm FDMR spectra were recorded with the temperature focal length lens onto the diamond samples, excit- stabilized so that temperature excursions were less than 0:05 K over 5 min. In order to avoid stray mag- ∗Electronic address: [email protected] netic fields when recording the spectrum, the cur- yElectronic address: [email protected] rents supplied to both TE elements were chopped 2 at a frequency of 2 Hz using photoMOS circuits, 1.000 ) and the spectra were recorded only when the TE . 0.999 m r currents were off. The process was repeated un- o n ( 0.998 til the temperature had been scanned through the e c n 0.997 280-330 K range several times in both directions. e c s The NV-ensemble magnetic-resonance spec- e 0.996 r o u troscopy has been described, for example, in Refs. l 0.995 [2, 18, 20–22]) and is only briefly summarized here. F 0.994 Optical pumping via a spin-selective decay path col- 2.855 2.860 2.865 2.870 2.875 2.880 lects NV centers (total spin S = 1) in the jms = 0i Microwave frequency [GHz] ground-state magnetic sublevel [19]. In the absence of external fields, the jms = 0i and jms = ±1i levels FIG. 1: Zero-field FDMR spectrum at 293 K for S3 are split by an energy equal to the axial zero-field and the corresponding fit based on Eq. 1 (solid green splitting (ZFS) parameter, D ≈ 2:87 GHz. For per- line). The six blue lines represent the fitted amplitudes fect C3v symmetry, the transverse ZFS parameter at each transition frequency, and the fitted linewidth was is E = 0 and the jms = ±1i levels remain degener- 3:3 MHz (full width at half maximum). The microwave ate. When the frequency of a microwave field that is power was reduced to ∼ −10 dBm to resolve the hyper- fine structure, resulting in the relatively small contrast transverse to the symmetry axis is tuned to the en- of ∼ 0:6%. The best-fit parameters for this scan are ergy splitting between the jms = 0i and jms = ±1i E = 4:1(2) MHz and D = 2866:8(2) MHz. levels, NV centers are transferred to the jms = ±1i sublevels, resulting in diminished fluorescence with a contrast as high as 30% [6]. In the presence of were fit to the spectra. Including residual magnetic an applied magnetic field, B, the jms = ±1i levels fields, measured by a commercial fluxgate magne- split, revealing resonances separated by 2gNV µBB, - tometer to be less than 1 G, into the model did not where gNV = 2:003 is the NV Landé factor [20, 23] significantly influence the fits. and µ is the Bohr magneton. For ensembles, there B Four single-crystal samples of mm-scale dimen- are four different NV orientations and, provided that sions were studied, which were labeled S2, S3, S5, g µ jBj D, only the projection of the magnetic NV B and S8 and characterized in Ref. [2]. Figure 1 shows field on the N-V axis affects the transition frequen- the FDMR spectrum at 293 K for S3, a sample syn- cies [22]. thesized by chemical vapor deposition (CVD) with The zero-field Hamiltonian for the ground state, [NV-]− ≈ 10 ppb [2]. As there was no applied mag- including hyperfine coupling to the 14N nucleus netic field, the splitting between resonance peaks is (spin I = 1), can be written as: due to non-zero E, induced by local strain [6, 22, 24]. 2 2 2 This feature is present in varying magnitudes for H0 ≈ DSz + E(Sx − Sy ) (1) all four samples. Even though all four NV orienta- +A S I + A (S I + S I ); k z z ? x x y y tions are present, the spectra are reasonably well- described by just six broad transitions, suggesting where Ak = −2:1 MHz and A? = −2:7 MHz are, respectively, the axial and transverse hyperfine con- that the strain splittings are spatially inhomogenous - stants [23]. Analysis of this Hamiltonian reveals six [17, 25]. As no correlation with NV concentration allowed microwave transitions for each N-V orienta- was observed (see Tab. I), further work is necessary tion. The relative intensities can be calculated by to determine the exact strain mechanism. treating the interaction with the microwave field, During each temperature scan, the spectrum was B~1 as a perturbation, H1 = gNV µBB~1 · S~, with ma- fit to an empirical function similar to the one de- trix elements that depend on the alignment of the scribed above, and the ZFS parameters were ex- microwave radiation with respect to the symmetry tracted. Figure 2(a) displays the spectra at two dif- and strain axes of each N-V center. However, since ferent temperatures for another sample, S8, a high- the exact geometry and the number of NV- centers pressure, high-temperature (HPHT) synthesized di- of each orientation were not known a priori, Gaus- amond with [NV-] ≈ 0:3 ppm, as well as the empir- sian functions with variable amplitudes and equal ical fits based on Eq. 1. Figures 2(b) and (c) show widths, centered about these transition frequencies, the ZFS parameters as a function of temperature for 3 ) . m - dD 1 dE −4 −1 r 1.00 5.3 # [NV ](ppm) (kHz=K) E(MHz) (10 K ) o dT E dT ] n ( z ç ç ç ççççç ç ç e 0.95 ç ç ç H ç ç çç ç ç ç çç ç S2 16 −71(1) 5:8(3) −1:7(5) c 5.2 ç ç ç ççç ç çç ççç çç çç ç ç ç ç ççççç n ç ç ççç ç ç çç ç ç M çç ç ç ç ç ç ç çç e ç ç ç [ ç ç c s 0.90 S3 0:01 −79(2) 4:3(2) 2(5) E 5.1 e r 326 K o u l 0.85 283 K a) 5.0 b) S5 12 −77(3) 11(1) −3:6(9) F 2.85 2.86 2.87 2.88 2.89 280 290 300 310 320 330 S8 0:3 −76(1) 5:2(1) −0:8(4) Microwave frequency [GHz] Temperature [K ] à 20 àà à àà ç àà à c) 2867 à à d) TABLE I: ZFS parameters and uncertainties for four ] àà 0 ç ç ààà à ç ç z à àà K dE/dT àà à à / H à à à à z -20 2866 ààà different samples.