We provide an overview of the experimental techniques, measurement modalities, and diverse applications of the quantum diamond microscope (QDM). The QDM employs a dense layer of fluorescent nitrogen-vacancy (NV) color centers near the surface of a transparent diamond chip on which a sample of interest is placed. NV electronic spins are coherently probed with microwaves and optically initialized and read out to provide spatially resolved maps of local magnetic fields. NV fluorescence is measured simultaneously across the diamond surface, resulting in a wide-field, two-dimensional magnetic field image with adjustable spatial pixel size set by the parameters of the imaging system. NV measurement protocols are tailored for imaging of broadband and narrowband fields, from DC to GHz frequencies. Here we summarize the physical principles common to diverse implementations of the QDM and review example applications of the technology in geoscience, biology, and materials science.
Nitrogen-vacancy (NV) centers in diamond are a leading modality for sensitive, high-spatial-resolution, wide-field-of-view imaging of microscopic magnetic fields. NV-diamond sensors operate in a wide range of conditions, from cryogenic to well above room temperature, and can serve as broadband detectors of slowly varying magnetic fields or as narrowband detectors of magnetic fields over a wide range of frequencies from near DC to GHz. Full vector magnetic field sensing is possible using the distribution of NV orientations along the four crystallographic directions in diamond.
NV centers function at ambient conditions, and have magnetically, electrically, and thermally sensitive electronic spin ground states with long coherence lifetimes. The NV spin state can be initialized, and the evolution of the spin states can be detected optically, thus allowing precision sensing of magnetic fields and other effects. Magnetic field sensitivity and spatial resolution are determined by the number of NVs in the sensing volume, the resonance linewidth, the resonance spin-state fluorescence contrast, the collected NV fluorescence intensity, and the NV-to-sample separation.
Variation of the experimental setup and measurement protocol allows NV-diamond magnetic imaging to be adapted for a wide range of applications in different fields of research. Although the desired capabilities for each magnetic imaging application vary widely, common requirements include good field sensitivity within a defined frequency range, fine spatial resolution, large field of view, quantitative vector magnetometry, wide-field and frequency dynamic range, and flexibility in the bias field and temperature during measurement. For example, imaging for geoscience  and cell biology  applications generally require high sensitivity to DC magnetic fields, spatial resolution at the optical diffraction limit, and room-temperature operation. In contrast, microelectronics magnetic field imaging  can require magnetic field frequency sensitivity up to the GHz range. Applications that do not require simultaneous imaging over a wide-field of view can also leverage scanning magnetometry using single NV centers at the tips of monolithic diamond nanopillars, or in nanodiamonds at the ends of atomic force microscopy cantilevers , , .
With proper optimization, NV-diamond magnetometry can offer combinations of the above capabilities unattainable using alternative magnetic imaging techniques. The magnetic force microscope (MFM) , while offering higher spatial resolution, is limited by small (<100 μm) fields of view, worse DC field resolution (>10 μT), and potential complications due to sensor-sample interactions. Meanwhile, the superconducting quantum interference device (SQUID) microscope, when measuring a sample at room temperature, can only achieve spatial resolution of >150 μm, although with excellent DC sensitivity . Finally, other techniques such as magneto-optical Kerr effect (MOKE) ,  and other Faraday effect-based magneto-optical imaging cannot produce reliable, quantitative maps of the vector magnetic field.
This review article provides an overview of the quantum diamond microscope (QDM), a common approach to ensemble NV wide-field magnetic imaging, and describes specific optimization of the QDM for several applications , . Schematics of typical QDM setups are shown in Figure 1. The QDM uses an optical microscope and a camera to measure the fluorescence from a thin ensemble NV layer at the surface of the diamond sensor chip, with the sample placed near to or in contact with the diamond. The local magnetic field of the sample is extracted from each camera pixel measurement, and a wide-field map of the magnetic field is constructed from the pixel array. We present the methods needed to image static and dynamic magnetic fields with the QDM, and briefly discuss imaging of temperature and electric fields. For each type of sample field – narrowband, broadband, etc. – we describe the quantum control procedures and hardware choices that are necessary for ideal imaging, and emphasize the design tradeoffs in optimal sensitivity and resolution limits that can be achieved.
2 NV physics relevant to QDMs
QDM implementation, including assembly and method of operation, depends on the intended application and the characteristics of the sample fields. However, there are principles of NV physics relevant to all QDM experiments. These principles rely on single-NV spin properties and their ensemble behavior. The QDM has optical, static magnetic, and microwave (MW) fields that are applied to manipulate the NV electronic- and spin-state populations in a controlled manner. An unknown sample field modifies the NV spin states and is detectable by changes in NV fluorescence. The three QDM driving fields are chosen to optimize coupling between the sample field and the NV spin state.
2.1 NV ground electronic state in the absence of external fields
Quantum control of NV centers with the QDM driving fields is possible because of the NV electronic- and spin-level structure , . An NV center consists of a substitutional nitrogen and an adjacent lattice vacancy defect in a diamond crystal. A negatively charged NV has six electrons, with two electrons from nitrogen, one electron from each of the three carbon atoms, and an additional electron from the lattice. These electrons occupy four sp3 atomic orbitals with electronic spin quantum number S=1. These sp3 orbitals linearly combine to form four molecular orbitals , comprising the ground electronic configuration. The lowest energy state of the ground configuration is the orbital singlet, spin triplet state 3A2, which has fine, Zeeman, and hyperfine structures, as shown in Figure 2. The four molecular orbitals also give rise to electronic excited states: orbital-doublet spin-triplet 3E, and spin-singlet orbital-singlet 1E, and 1A1 shown in Figure 3.
NV magnetometry uses fluorescence from electronic state transitions to detect changes to the 3A2 ground state configuration that result from coupling to a sample field. Therefore, focus is placed on the physics of the 3A2 Hamiltonian. NV centers have C3v point-group symmetry, and are spatially invariant under the C3v symmetry transformations (the identity, two 120° rotations about the NV axis, and three vertical reflection planes). NV centers also have a built-in quantization axis along the NV axis (called the NV z-axis, or the crystallographic  direction). The 3A2 electronic ground state is an orbital singlet and spin triplet manifold, with ground state Hamiltonian 
where h is Planck’s constant, and and are the dimensionless electron and nitrogen nuclear spin operators, respectively. The first term is the fine structure splitting due to the electronic spin-spin interaction, with the fine structure tensor D . The second term is the hyperfine interaction between NV electrons and the nitrogen nucleus, with the hyperfine tensor A. The third term is the nuclear electric quadrupole interaction, with the electric quadrupole tensor Q. Under the C3v symmetry of the NV center, D, A, and Q are diagonal in the NV coordinate system , , and can be written as 
D(T) is the fine-structure term called the zero-field splitting (ZFS), A‖ and A⊥ are the axial and transverse hyperfine terms, and P is the nuclear electric quadrupole component. Two important features of the ground state are evident from the Hamiltonian. First, the 3A2ms=±1 magnetic sublevels and the ms=0 level have D(T) difference in energy. D(T) is temperature dependent because the spin-spin interaction changes with the lattice constant , , with D≈2.87 GHz and dD/dT=−74.2 kHz/K at room temperature. Second, the 3A2 electronic states have an additional hyperfine energy splitting due to the nitrogen nucleus. I=1 for a 14N nucleus, while I=1/2 for a 15N nucleus. The energy level diagrams for 14N and 15N are shown in Figure 2. The hyperfine parameters are and .
Here, Mx, My, Mz, Nx, and Ny are stress-dependent amplitudes. The Mz term contributes to the ZFS, while the other terms may be negligible or suppressed depending on the experimental conditions (such as an applied magnetic field ). The NV spin sensitivity to this spin-stress-induced interaction can be used to image internal or external diamond stress , which is important for diamond material characterization. However, for imaging external magnetic fields, we consider NV stress sensitivity as a potential limitation.
2.2 NV electronic transitions
The first excited electronic configuration of the NV has an orbital-doublet, spin-triplet state, 3E, shown in Figure 3. The two orbital states and three spin states of 3E combine to form six fine-structure states that reduce to three states at room temperature , resembling the 3A2 state. 3E is coupled to the 3A2 ground state by an optical 637-nm zero-phonon line (ZPL). The 3E↔3A2 is a radiative transition that generally conserves the electron spin state ms as a result of weak spin-orbit interaction . The 3E→3A2 (3A2→3E) transition works for longer (shorter) wavelengths in fluorescence (absorption) as a result of the phonon sideband (PSB). This behavior is similar to Stokes and anti-Stokes shifted transitions . Figure 3 also shows the radiative spin-conserving 1E↔1A1 transition, which has an infrared 1042-nm ZPL and its own sideband structure.
Nonradiative transitions between states of different spin multiplicity exist between 3E and 1A1 and between 1E and 3A2. These nonradiative transitions are caused by an electron-phonon-mediated intersystem crossing (ISC) mechanism and do not conserve spin. The probability of the ISC transition occurring for 3E to 1A1 is only non-negligible for the ms=±1 states of 3E and is characterized by the ISC rate of transition . Similarly, the ISC transition probability from 1E to the ms=0 state of 3A2 is approximately 1.1–2 times that of the ISC transition from 1E to the ms=±1 states of 3A2 , . These state-selective differences in the ISC transition rate allow for spin polarization of the NV under optical excitation from 532-nm laser illumination.
2.3 Optical pumping and spin polarization
An optical driving field from a pump laser is applied in order to spin polarize the NV electronic state at the start of a QDM measurement. This pump laser is also used at the end of a measurement to read out the final NV spin state through the fluorescence intensity. NV optical pumping takes advantage of the ms-selective nonradiative ISC decay pathway , . An NV that is optically excited from 3A2 to 3E state by a 532-nm photon, decays along either the optically radiative 3E→3A2 pathway or the non-optical, ISC-mediated 3E→1A1→1E→3A2 pathway. The ms-selectivity of the ISC will preferentially depopulate the ms=±1 spin projection states. NVs starting in the 3A2ms=±1 sublevel are eventually pumped (on average, after a few pump photon absorption cycles) into the 3A2ms=0 sublevel. Typically only ∼80% of NVs in an ensemble can be initilallized into the ms=0 state , where they remain in a cycling transition. The 1E state is metastable with a ~200 ns lifetime at room temperature , . The 3E upper state has a t3E≈13 ns lifetime , , and the 3A2→3E absorption cross section at λ=532 nm  is σ=3.1×10−17 cm2 (although there is disparity in the reported 532-nm absorption cross-section value and saturation intensity , ). These correspond to (hc)/(λσt3E)≈0.9 MW/cm2 saturation intensity, where c is the speed of light.
The ISC is also responsible for the reduced fluorescence intensities of NVs in the ms=±1 sublevels, since they emit fewer optical photons when returning to the 3A2 state through the ISC-mediated pathway. The fractional fluorescence difference between NVs in the ms=±1 sublevels and NVs in the ms=0 sublevel is called the fluorescence contrast, which can be as large as ~20% for a single NV . The fluorescence intensity from an optically pumped NV diamond chip therefore indicates the percentage of the NVs in the ms=0 state or in the ms=± 1 states. A transition of NVs from the ms=0 to the ms=±1 state, e.g. induced by a resonant MW field, drops the fluorescence as more NVs follow the ISC-mediated decay transition. This is the mechanism underlying optical readout for QDM imaging.
2.4 Microwave driving field
A MW driving field resonant with the ms=0 to +1 or −1 transitions induces Rabi oscillations, transferring the NV population from one sublevel to the other and creating superpositions of ms states. Either a continuous-wave (CW) or a pulsed MW field can be used. The length of the MW pulse determines its impact on the NV population: π pulses are of sufficient duration to transfer the NV population from ms=0 to ms=1 when the NVs are initialized in the ms=0 state; π/2 pulses are of duration to create an equal superposition of ms states. The utility of pulsed MW fields for QDM detection of different types of sample fields will be discussed below.
Applying resonant CW MWs simultaneously while optically pumping of the NVs to the ms=0 sublevel results in MW-induced transfer of the NV population out of the ms=0 sublevel, spoiling the optical spin polarization and decreasing the emitted fluorescence intensity. Measuring the NV fluorescence intensity as a function of the probing MW frequency is called optically detected magnetic resonance (ODMR) spectroscopy . Simulations of ensemble NV ODMR spectra for NVs with 14N and 15N isotopes are shown in Figure 2. The known dependence of the 3A2 sublevel energy on external fields allows conversion of these ODMR spectra into magnetic field, electric field, temperature, and crystal stress information.
2.5 Static magnetic bias field and Zeeman splitting
A static magnetic field B0 causes a Zeeman interaction in the 3A2 state, written as
Here, μB=9.27×10−24 J/T is the Bohr magneton, g is the electronic g-factor tensor (which is nearly diagonal), ge≈2.003 is equal to the NV center’s electronic g-factor , and γ=geμB/h is the NV gyromagnetic ratio. The Zeeman interaction lifts the degeneracy between the ms=±1 sublevels, and for |B0| along the NV axis, the ms=±1 sublevel energies split linearly with |B0| while the ms=0 sublevel is unaffected. The nuclear Zeeman terms are considered negligible and have been excluded.
A sufficiently large bias magnetic field makes the Zeeman term dominant in the Hamiltonian. Otherwise, terms including stress and electric field would dominate with the Zeeman term acting as a perturbation, reducing magnetic field sensitivity and complicating the data analysis. For magnetic imaging, both the static bias fields, which are part of the QDM, and the sample magnetic fields contribute to the Zeeman interaction.
2.6 Sample fields
QDM experiments create a two-dimensional image of the magnetic fields from a sample containing a distribution of magnetic field sources. It is also possible to image a sample’s temperature distribution and electric fields.
The sample magnetic field is generated by field sources, such as current densities or magnetic dipoles, with either known or unknown distributions. Measurement of the sample magnetic field can be used for the inverse problem of estimating an unknown source distribution under certain conditions , , . The form of the sample magnetic field in terms of its sources is
where J is the current density of the sample, R=r–r′ is the distance from a magnetic source at r′ to an observation point at r, R=|r–r′|, and Equation (5) is the Biot-Savart law for static fields and applies to fields in the quasi-static regime for which the characteristic system size is small compared to the electromagnetic wavelength. A sample consisting of small magnetic particles will have a magnetic field composed of single dipole fields of the form
Here, m is the magnetic moment, n=r/r, and the delta function only contributes to the field at the site of the dipole r=0. Other typical sample fields, such as the narrowband magnetic field from the Larmor precession of protons, can also be derived . Figure 4 shows examples of the magnetic fields for a current distribution and a distribution of magnetic dipoles.
The time dependence of the sample magnetic field will determine the QDM measurement protocol. Static and quasi-static sample fields will contribute to the NV Hamiltonian by an additional term in Eq. (4):
where Bs is the magnetic field of the sample, which can take the forms given in Eqs. (5) and (6). The magnitude of the sample field along the NV axis is therefore determined by changes in separation in the ODMR resonance features that result from Bs in addition to the effect of the bias field. ODMR measurements, with and without the sample, then allow the determination of the unknown Bs field. Reasonable assumptions can be made to determine Bs without having to take multiple measurements . A time-varying sample magnetic field with frequency components near the 2.87-GHz ZFS will in turn induce NV spin transitions if the B0 bias field has been tuned to the appropriate Zeeman splitting. Sweeping the B0 field will then locate the frequency of the sample fields, with the magnitude determined by the ODMR contrast depth and linewidth.
Electric field and temperature distributions from the sample will also change the NV spin states. The external sample electric field Es=(Esx, Esy, Esz) adds to the internal local electric fields  Eloc=(Eloc,x, Eloc,y, Eloc,z) in the diamond, e.g. induced by a high density of P1 (nitrogen) centers, such that Etot=Es+Eloc contributes to the Hamiltonian in Eq. (2) :
Here, d‖ and d⊥ are the coupling constants related to the NV electric dipole moment, and d‖=d⊥ with d‖ =3.5×10−3 Hz/(V/m) and d⊥ =0.17 Hz/(V/m) , . For the typical scale of sample electric fields, coupling to the NVs is small compared to sample magnetic fields of interest. Hence electric fields do not cause noticeable shifts in ODMR resonances for most QDM magnetic imaging experiments. External temperature variations, e.g. from the sample, couple to the NV by the temperature dependence D(T) of the ZFS . Changes in temperature of the diamond due to the sample temperature field will therefore result in a common mode shift of the ODMR resonance, which is distinct from the effect of magnetic fields.
2.7 NV ground-state Hamiltonian
Detecting the resultant spatial distribution of changes in the ODMR spectra across an NV ensemble due to spatially varying sample fields is the principle underlying QDM high-resolution imaging. The ground-state Hamiltonian necessary to capture the relevant dynamics of single NVs for QDM imaging can be summarized by combining Eqs. (2), (3), (7), and (8):
Equation (9) summarizes the interaction between an NV center and temperature, magnetic field, and electric field. This equation demonstrates that NVs can in principle be used to image all of these quantities. For simplicity, Eq. (9) does not include the nuclear electric quadrupole interaction and the comparatively negligible terms of the crystal strain interaction. The hyperfine splitting terms from Eq. (2), also excluded for simplicity, are important and visibly evident in the hyperfine splitting of the ODMR resonances in Figure 14. Further simplifications can be made to Eq. (9) depending on the magnitude of the parameter of interest and the bias fields.
2.8 Behavior of NV ensembles
QDMs use ensembles of NVs to obtain simultaneous measurements over a wide-field of view. Ensembles yield stronger signal than single NVs because of the larger number of NVs contributing photons to overall fluorescence but introduce ensemble behavior that can worsen contrast compared to single NV performance. Intuition about single NVs also does not simply extend to the ensemble case: the NV spin-ensemble behavior can be substantially different from single-NV behavior. The complex spin bath environment of diamond contributes several mechanisms to NV spin ensemble dephasing and decoherence, which ultimately limit the magnetic field sensitivity of NV ensembles. These ensemble dynamics are characterized by the longitudinal (spin-lattice) relaxation time T1, the transverse (decoherence) relaxation time T2, and inhomogeneously broadened transverse (dephasing) relaxation time . In particular, understanding and minimizing ensemble NV dephasing is critical for high-performance broadband and static field magnetic imaging with QDMs, informing both diamond material design and quantum control techniques. This topic has been treated extensively in the literature , , and will be discussed later in this article, after introduction of NV measurement protocols.
3 NV measurement protocols relevant to QDMs
The large toolbox of QDM sensing protocols allows for imaging magnetic fields over a wide range of characteristic timescales. Figure 5 shows a schematic of the most commonly used sensing protocols. The interplay and timing of the laser pulses and MW pulses determines the basic properties of the techniques introduced in this section.
3.1 DC magnetometry: static and low-frequency fields
Three established QDM sensing protocols exist for measuring static (DC) and slowly varying magnetic fields: CW ODMR, pulsed ODMR, and Ramsey magnetometry. These protocols have been used to sample time-varying magnetic fields up to 1 MHz in a single-pixel experiment .
3.1.1 CW ODMR
CW ODMR is a robust and simple method that can image the vector components of a magnetic field in the QDM modality. Because of its easy implementation, CW ODMR is the most common technique used for QDM applications. Continuous laser pumping, MW driving, and fluorescence readout are performed simultaneously, as shown in Figure 5. The laser is used to both pump the NVs into the ms=0 spin state and probe the spin states of the population via NV fluorescence. The frequency of the MW drive is swept in time and synchronized with the readout. A decrease in fluorescence occurs when the MW frequency matches the NV resonance due to the spin-state dependence of NV photon emission described in Section 2.2.
Figure 6A shows an example where a change in B0 shifts the line center of the resonance feature. For an NV ensemble, the resonance lineshape – often modeled as a Lorentzian or Gaussian – is parameterized by the center frequency, linewidth, and fluorescence contrast. The center frequencies of every NV resonance feature are fit to the appropriate Hamiltonian to extract the desired magnetic field, strain, temperature, and electric field. In a magnetic imaging experiment, this analysis yields B0+Bs, from which the magnetic field of the sample can be determined , , .
Measuring the entire resonance spectrum in CW ODMR limits the sensitivity and the temporal resolution of the measurements because of the significant fraction of experiment time spent interrogating with probe frequencies that yield no signal contrast. Sparse sampling of the resonance spectrum can improve the sensitivity of the measurement by minimizing the dead time. An extreme version of sparse sampling can be achieved using a lock-in modality where the probe frequency is modulated between the points of maximum slope of an ODMR resonance feature . This technique has been extended to monitor multiple ODMR features simultaneously to extract the vector magnetic field by modulating at different frequencies . Frequency-modulated ODMR has been performed with bandwidths up to 2 MHz, but was demonstrated on a small volume and required high laser and MW intensity beyond that typically employed with QDMs .
3.1.2 Pulsed ODMR
CW ODMR suffers from laser repumping of the NV spins through the entire measurement. This simultaneous laser pumping and MW drive spoils the measurement sensitivity as a result of the competing processes of initializing the spin state (laser) and driving transitions (MW drive) . In order to mitigate this power-broadening, a pulsed ODMR protocol uses a temporally separated laser initialization, an MW control π pulse, and a laser readout pulse, as demonstrated in Figure 5. This leads to the decreased linewidths shown in Figure 6B as compared to CW ODMR. Alteration of the MW power changes the necessary duration of a π pulse, and must be optimized to balance the linewidth and contrast of ODMR resonance features .
3.1.3 Ramsey spectroscopy
Ramsey spectroscopy  determines the magnitude of a DC magnetic field by measuring the relative phase accumulation between the different electronic spin states prepared in a superposition with a π/2 pulse , . A green laser pulse initializes the spin state into the ms=0 state to begin the sequence. Next, a resonant MW π/2 pulse prepares the spin into a superposition of the ms=0 and ms=±1 spin states (or ms=−1 depending on the drive frequency). The system is allowed to evolve under the relevant Hamiltonian for a free precession time τ. In the simplified scenario where the dynamics are dominated by the magnetic field, the NV superposition state accumulates a phase ϕ=2πγ(|B0+Bs|)τ. A second MW π/2 pulse is applied to project the accumulated phase information onto the relative population of ms=0 and ms=+1 spin states. A second laser pulse is applied to measure the spin state population through the spin-dependent fluorescence of the NV.
To obtain information about the magnetic field |B0+Bs|, a Ramsey pulse sequence is repeated several times, sweeping the free precession interval such that each measurement is taken for different τ values. The resultant fluorescence contrast signal as a function of τ is known as the Ramsey fringes, illustrated in Figure 6C. By taking the Fourier transform of these fringes, one can locate the position of the dominant frequencies and determine the deviations from those set by the bias magnetic field that result from the sample field.
Mapping out the Ramsey fringes is inefficient with respect to the speed of measurement, similar to the inefficiency of the frequency sweep for ODMR. Instead of mapping out the full fringe and taking the Fourier transform, the free precession time τsense is fixed to sample the Ramsey fringe at the point of maximum sensitivity, which is the point of maximum slope closest to This process maps out a magnetometry curve, as illustrated in Figure 6D. The steeper the slope of the magnetometry curve, the more sensitive the protocol.
A key feature of Ramsey magnetometry is having both the laser and MWs switched off when the NV electronic spin is accumulating phase via interaction with the magnetic field. The Ramsey protocol is consequently not vulnerable to the power-broadening that impacts CW ODMR and allows the use of high MW and laser power to increase the sensitivity .
Other benefits of Ramsey magnetometry over CW and pulsed ODMR is that it more efficiently leverages protocols that mitigate dephasing such as spin-bath driving and is compatible with sensing in the double quantum basis .
3.2 AC magnetometry: narrowband fields
A QDM can measure narrowband oscillating magnetic fields using AC magnetometry sequences, including Hahn echo and dynamical decoupling. These pulse sequence protocols act as frequency filters and allow the QDM to operate as a sensitive lock-in detector  of these AC fields. The frequency range of narrowband signals that are detectable with NV AC magnetometry is ~1 kHz to ~10 MHz, limited at the low end by NV decoherence and at the high end by the amplitude of fast MW pulses that can be realistically applied to an NV ensemble.
3.2.1 Hahn echo
The addition of a π pulse into the middle of a Ramsey sequence mitigates environmental perturbations that are slow compared to the free precession interval between pulses . This pulse sequence is known as the Hahn Echo sequence , , and results in the refocusing of NV ensemble dephasing such that the limiting measurement timescale becomes the decoherence time T2 rather than the dephasing time The consequence is improved magnetic field sensitivity (discussed in Section 4.1), especially for lower frequency signals, because T2 typically exceeds by over an order of magnitude . Figure 7A demonstrates a decoherence curve when using a Hahn echo pulse protocol. The spacing between MW pulses acts as a narrowband filter in frequency space. The width of this filter is given by the filter response function . Hahn echo uses only one π pulse and therefore has a fairly broad filter, allowing for the sensing of a wide bandwidth of magnetic field frequencies.
To optimally sense external oscillating fields, the spin evolution time is set to be ~T2; however, the frequency of the sensed magnetic field can lead to operation with a non-optimal τ . For a fixed spin evolution time, a change in magnetic field will lead to a difference in phase accumulation that maps onto the total fluorescence, Figure 7B.
3.2.2 Dynamical decoupling
Building upon the Hahn echo sequence, dynamical decoupling techniques commonly apply multiple refocusing pulses with spacing determined by the period of the sample field , . These additional refocusing pulses result in an advantageous extension of the decoherence time by narrowing the width of the filter response function and reducing sensitivity to magnetic noise outside the bandwidth. In particular, decoupling of the NV from spin-bath-induced magnetic noise improves with additional pulses, though at the cost of making the technique sensitive to a narrower range of frequencies , . The extension in the decoherence time, Figure 7A, can lead to a dramatic improvement in magnetic field sensitivity, Figure 7B. Dynamical decoupling also increases the time during which NVs can interrogate the sample field Bs optimally towards the extended decoherence time.
The Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence is a dynamical decoupling sequence employing π pulses, which rotate the NV about the same axis as is polarized by the initial π/2 pulse. Another common sequence, XY8, extends this by choosing the rotation axis for each π pulse in order to suppress the effects of pulse errors. A large family of similar sequences exist, many well known in NMR, to improve NV sensing through more efficient, robust control of the NV electronic spin state .
3.3 Resonant coupling to external GHz fields
Applications that require measurement of GHz-scale oscillations can leverage interactions between the NV and magnetic signals near the NV resonance as a probe , , . CW ODMR constitutes a simple protocol that can be used in this manner. Measurements of the contrast and linewidth enable the determination of the optical power and MW power that broaden the lines, Figure 8A, in addition to other mechanisms that contribute to the inhomogeneous dephasing of the ensembles . However, this method is not very sensitive and difficult to quantify because of the various ways the contrast and linewidth can vary over a field of view . Alternative methods to CW ODMR include Rabi driving and T1 relaxometry.
3.3.1 Rabi driving
Use of a Rabi sequence provides a more direct way to determine local magnetic fields oscillating at or near GHz frequencies as compared to CW ODMR . Similar to the previously discussed protocols, the NV spin state is initiated to the ms=0 state with a green laser light. A MW drive is left on for a varying amount of time. If the MW drive is in resonance with the NV spin state transitions, for example, between ms=0 and ms=±1 states, the population will be driven back and forth between the spin states. The strength of this GHz drive determines the rate at which the transition is driven. This rate is called the Rabi frequency, and scales with the square root of the input microwave power. Figure 8B illustrates the increase in Rabi frequency as a function of increasing amplitude of the MW driving field.
3.3.2 T1 relaxometry
When the NV is initially polarized into the ms=0 state with green illumination, there is a characteristic timescale over which the spin population decays back to a thermally mixed state. This timescale is the longitudinal (spin-lattice) relaxation time T1 and can be up to 6 ms when dominated by phonon interactions at room temperature . However, T1 can be spoiled by the presence of magnetic frequency noise or other paramagnetic spins at the NV resonance frequency, as shown in Figure 8C . The local bias field can be swept to change the noise frequency to which the measurement is sensitive.
4 QDM performance
QDM performance characteristics include magnetic field sensitivity, temporal resolution, frequency bandwidth, spatial resolution, and field of view of the sample field. These characteristics depend on the sensing protocols of the QDM discussed in Section 3, which in turn are determined by the spectral and spatial qualities of the sample fields to be imaged. This section focuses on the physical limits to performance; the performance impact resulting from the use of different experimental components for the QDM is treated in Section 5.
4.1 Magnetic field sensitivity
The minimum detectable field difference is defined as the change in magnetic field magnitude δB for which the resulting change in a given measurement of the field equals the standard deviation of a series of identical measurements. However, characterizing the minimum detectable field difference must consider the total measurement duration as well as the total number of NVs that contribute to the measurement for meaningful determination of sensor performance. The magnetic field sensitivity scales as the square root of the number of detected photons. The number of photons collected over unit time from a unit volume of NVs increases proportionally with time and volume. To account for measurement time, sensitivity is represented as with units of T Hz−1/2, where tmeas is the total measurement time. To account for the number of NV spins required to reach a given sensitivity, a volume-normalized sensitivity is defined as with units T μm3/2 Hz−1/2, where V is the volume for a fixed density of NVs , .
CW ODMR magnetometry is the most widely used QDM measurement technique because of its simplicity. The sensitivity of a CW ODMR magnetometry sequence is characterized by the slope of the resonance line ∂I/∂ν0, with fluorescence intensity I and frequency ν0, and the rate R of photon detection from a cubic micrometer of NVs. The CW ODMR shot-noise-limited sensitivity is
where C is the contrast and Δν is the linewidth of the ODMR resonance. The resonance line shape is typically fit by a Lorentzian, giving the factor for the maximum slope. The relationship between the ODMR linewidth and the previously defined dephasing time is approximated by , , , . When performing ensemble measurements, many mechanisms can contribute to the linewidth, as demonstrated in Figure 9.
The sensitivity of CW ODMR magnetometry is limited by laser and MW ODMR lineshape power broadening. Solving the Bloch equation for a simplified two-level model yields the contrast, linewidth, and volume-nomralized magnetic sensitivity, shown in Figure 10. The calculations are based on  for CW ODMR using parameters from Table 1, Figure 10 displays a broad range of laser and MW powers to indicate how these affect the sensitivity. The tradeoff between laser and MW power limits the achievable volume-normalized sensitivity of CW ODMR, precluding simultaneous optimal contrast and narrow linewidth. Applications that require higher temporal and spatial resolution must use techniques more sensitive than CW ODMR.
|Diamond [N]||nN||NV/N||nNV,SA||Single axis (counts/s)||Photon rate (counts/s)||Counts per readout||Contrast||(μs)||T2 (μs)|
Properties of two notional diamonds used for performance simulations for a 1-μm3 QDM voxel. [N] is the concentration of nitrogen in the diamond. nN is the number of nitrogen atoms per 1-μm3 voxel, NV/N=10% of the N atoms are NV centers. A single NV axis is considered, giving nNV,SA with 105 fluorescence counts/s for each NV in a 1-μm3 QDM voxel. More NV centers per voxel increases the magnetic field sensitivity because the rate of photons emitted scales with nNV−. Counts per readout are for an assumed readout time of 300 ns. Assumed scaling of and T2 are from Ref. .
Ramsey magnetometry achieves the best DC magnetic field sensitivity of the QDM protocols because of its pulse scheme: the NV spins interrogate the sample fields during an interval without simultaneous interaction of the optical and MW driving fields. The shot-noise-limited sensitivity for DC magnetic fields using a Ramsey pulse is 
where N is the number of noninteracting NVs contributing to the measurement S=1/2 spins, and Δms is the generalization to >1 spin state difference used for measurement (e.g. Δms=2 for the NV ms=−1 to ms=1 transition when operating with a double-quantum coherence ), C is the resonance contrast, n is the average number of photons collected per NV per measurement, τ is the spin interrogation time, and t1 and tR are the optical spin-state initialization and readout times, respectively (tmeas=tI+τ+tR). The spin-projection-noise-limited sensitivity is given by the first two terms of Eq. (11). It is evident that longer interrogation time τ and larger number of spins N allow for better sensitivity to small magnetic fields. However, several factors cause Ramsey magnetometry to fall short of this limit: a decrease in sensitivity due to spin dephasing with characteristic time is accounted for in the exponential term with parameter p depending on the origin of dephasing; imperfect readout contributes the first square root term; and the reduced fraction of total measurement time allocated for spin interrogation due to the overhead time from t1 and tR is accounted for in the last term. Optimal DC sensitivity is achieved for . Figure 11 compares the sensitivity of Ramsey magnetometry as a function of the frequency of the field being measured for the two diamonds in Table 1.
The sensitivity for measurement of AC magnetic fields using the Hahn echo protocol is
Hahn echo magnetometry builds on the Ramsey protocol as discussed in Section 3.2.1, resulting in similar physics underlying the AC magnetic field sensitivity to that of DC fields. The additional MW π pulse in the Hahn echo sequence refocuses the dephasing NV ensemble such that the sensing duration τ approaches the spin decoherence time T2. Because T2 is at least an order of magnitude longer than the spin dephasing time , the sensing duration increases translating to an improvement in sensitivity. AC sensing protocols are thus limited by T2, whereas DC sensing protocols are limited by because the AC protocols can generally achieve better sensitivity than DC protocols. However, the benefit of being T2-limited can be degraded by coherent interactions between the NV spin ensemble and other spin impurities, which decrease the T2 coherence time. The optimal spin interrogation time τ for Hahn echo magnetometry is τ~T2; additionally, τ should match the period of the AC magnetic field TAC. Consequently, maximum sensitivity is achieved for AC magnetic fields with TAC~T2 of the diamond.
CPMG pulse sequences improve the sensitivity by extending T2 even further :
where k is the number of pulses, τ=k/(2f0) is the full spin evolution time, and f0 is the frequency of the sample field. The optimal number of pulses for a given frequency is given by with s~2/3 and p~3/2, and is set mostly by the spin bath dynamics . The measurement time increases linearly with increased number of pulses, whereas the decoherence time T2 increases sublinearly; the optimal number of pulses balances these effects . Extensions of Eq. (13) exist to take into account multipulse dynamical decoupling protocols.
4.2 Temporal resolution and frequency bandwidth
QDM’s temporal resolution is defined as the time required between subsequent measurements of the sample field. The physical limitation determining the fastest temporal resolution is set by the time it takes for the NVs to react to a change in the sample field. The temporal resolution can never be faster than ~5 MHz (the maximum optical pumping rate), which is limited by the 1E metastable state lifetime of 200 ns, discussed in Section 2.3. The same is true for pulsed measurements, because NVs are optically reinitialized to the ms=0 state before each measurement. For a measurement with continuous laser illumination and MW field, the NV temporal resolution is set by the optical pumping rate and the MW Rabi frequency. There is also a practical limit to the temporal resolution, set by signal-to-noise ratio (SNR) tolerance: faster temporal resolution gives worse SNR per measurement.
The NV sensor frequency bandwidth is the range of sample frequencies that can be interrogated with the same experimental protocol. A DC magnetometry experiment has a frequency bandwidth spanning DC to the NV temporal resolution cutoff. A dynamical decoupling AC magnetometry experiment has an approximate frequency bandwidth which is roughly equal to 1/Ttot (the Fourier limit), where Ttot is the duration of the dynamical decoupling pulse sequence. An AC magnetometry measurement based on driving the spin population between ms sublevels (Rabi) or spoiling of the initialized spin state (T1 relaxometry) has a frequency bandwidth corresponding to the NV resonance linewidth, i.e. the frequency span over which the NVs are in resonance with the MW field, which is
AC magnetometry sequences that are based on pulse control of the NV spin state have a frequency bandwidth dictated by the filter function S(f) for the specific pulse sequence being used. The center frequency and bandwidth are defined by the number of pulses, k, and the spacing between the π pulses, τ . The center frequency of the filter is given by f0=1/2τ. For a sequence of k pulses, with total measurement time T=kτ, the width of the filter function is given by Δf=1/T=1/kτ. The filter function S(f) depends upon a protocol-specific response F(fT):
where an example response function for the CPMG protocol is 
Figure 12A demonstrates the need to change the number of pulses in order to operate at the sensitivity limit shown in the CPMG curves in Figure 11. Figure 12B gives the filter functions for the most sensitive points along the curves presented in Figure 11.
It is tempting to conflate temporal bandwidth and frequency bandwidth, but they in fact represent different characteristics. For example, an NV T1 measurement can detect magnetic noise across a few MHz frequency bandwidth around a central frequency ranging from near zero to many GHz (depending on the applied B0), but measurements may only be repeatable at <1 kHz (temporal resolution). Only in the case of DC magnetometry protocols do the temporal and frequency bandwidth correspond to the same sensor property.
4.3 Spatial resolution and field of view
QDM magnetic imaging seeks to resolve magnetic fields with high spatial variation over a wide-field of view and to successfully invert the magnetic field measurements to generate a map of closely separated magnetic sources. Both goals have fundamental and sensor-specific limitations. It is ideal to operate at the limit of magnetic field inversion and not to be limited by the sensor properties such as resolution and field of view.
The magnetic inversion problem does not generally have a unique solution. Only if the current distribution is limited to two dimensions (2D) can the inverse problem be solved uniquely from a planar measurement of the magnetic field. A magnetometer must sample the field at discrete points in a 2D plane with a sufficient sampling density to recover the continuous magnetic field created by the sample sources. The spatial resolution that can be obtained from this 2D map of the field is then limited by the offset distance between the measurement plane and the sources, as well as by the noise in the data , . In general, the offset distance should be as small as, or smaller than, the characteristic length scale of the magnetic field sources, as shown in Figure 4, for reliable inversion of the magnetic image to the source distribution. In analogy to the Nyquist sampling theorem, the pixel size sets the maximum spatial (k-space) frequency. The size of the field of view sets the spatial frequency resolution (again by a Fourier transform argument). Both of these effects impact the ability to perform magnetic field inversions and map the underlying sources .
The in plane pixel size is made too small, the noise level could preclude detection of the magnetic fields of interest. This is similar to the negative impact to δB that can result from pushing the temporal resolution, discussed in the previous section. On the other hand, if the pixel size is too large, then small-length-scale signals of interest will be blurred out and the fidelity of the magnetic field amplitude will be degraded. Figure 13 illustrates an example of this tradeoff for magnetic fields simulated in Figure 4.
The QDM’s spatial resolution is set by the following:
NV-sample standoff distance. As the standoff distance Δz increases, the 2D magnetic map is convolved with a Lorentzian of width Δz, reducing the ability to resolve closely separated magnetic sources . Reducing the standoff distance improves the field strength and sometimes the spatial resolution.
NV layer thickness. A thick NV layer has a layer-sample separation Δz corresponding to somewhere between the NVs in the layer nearest to and farthest from the sample. An NV layer that is thin compared to the sample could have better sensitivity than with a thicker NV layer; the far-standoff NVs will measure a Bs comparable to that of the near-standoff NVs, and the photon shot noise improves with a thicker layer. Conversely, an NV layer that is thick compared to the sample will have far-standoff NVs that measure almost no field but add background fluorescence and can cause deleterious artifacts .
Optical diffraction limit, set by the numerical aperture (NA) of the microscope objective λ/(2NA) for a typical fluorescence wavelength of λ≈700 nm. This assumes that the camera pixel size is small compared to the diffraction-limited spot size in the image plane. The spherical aberration from the diamond chip or other optics can also degrade the resolution.
QDM magnetic field imaging is best used for applications that need high spatial resolution over a wide-field of view and can afford small NV-sample separation. The intuitive rule of thumb is to have NV layer thickness, standoff distance, and sample thickness of comparable sizes.
5 QDM components and design choices
The goal of an experimental designer is to make sure that the application space of an experiment is limited by the fundamental physics of the system and not the equipment used. However, this is not always possible because of the availability and current state of technology. This section considers equipment choice and its impact on reaching the field sensitivity, temporal resolution, and spatial resolution presented in Section 4 for different protocols.
Optimal performance for a given QDM target application can only be achieved with informed equipment choices. These choices include proper selection of the diamond, bias magnetic field, MW field, optical illumination, optics, camera, and sample mounting. The QDM components and their impact on the QDM performance is presented here with focus on informed hardware choices that enable operating the QDM at the physics-limited sensitivity and performance.
Properties of the diamond chip that impact QDM performance, include NV layer thickness, NV concentration, isotope and impurity concentration, and diamond cut. These properties are controlled during the diamond fabrication process. Single-crystal diamond substrates used as the platform for QDM imaging are grown in one of two ways. One technique, high-pressure high-temperature (HPHT) growth, resembles natural diamond formation, requires an anvil press at ~1700 K and 5 GPa, and produces diamond samples with ~100 ppm nitrogen density. The second technique, chemical vapor deposition (CVD), grows diamond substrates layer by layer from a plasma, and yields diamond samples with low ppb nitrogen concentration.
Imaging a thin two-dimensional magnetic sample is optimal when the NV layer thickness is comparable to the magnetic sample thickness as discussed in Section 4.3.
The typical NV layer thickness for QDMs ranges from ~10 nm to ~100 μm. There are several methods available to make NV layers of varying thickness.
N+ or is implanted in a type IIa diamond with ppb impurity density to create a ~10–100 nm shallow layer. Annealing the diamond improves the NV yield and NV density .
A ppm-density nitrogen-rich layer is grown on top of a type IIa diamond substrate using CVD. After growth, electron irradiation of the diamonds introduces vacancies, and annealing improves the NV yield by converting substitutional nitrogen atoms (P1 centers) into NVs with a ~10% conversion rate . The nitrogen-rich layer can be from several micrometers down to several nanometers in thickness .
Instead of irradiating in #2, naturally formed NVs are used and can be preferentially oriented along one of the crystallographic directions (instead of equal NV fractions along all four orientations). Removing three of the NV orientations can improve the magnetic field sensitivity by ~2×, but can come at the expense of reduced NV density and fluorescence .
Similar to #2, nitrogen is temporarily introduced during CVD diamond growth to create a nitrogen-rich layer of a few nanometers, followed by a nitrogen-free diamond cap layer. NV centers are then created by electron irradiation and annealing. This technique is called delta doping . The cap layer adds to the standoff distance, so the surface layer version in #2 is often preferred, or the cap layer is etched away .
The NV density in the NV layer is optimized to achieve a desired magentic field sensitivity. High NV density yields more fluorescence intensity and good photon shot noise. However, the greater density of P1 paramagnetic impurities – required for high NV yield – contributes to magnetic inhomogeneity, thereby broadening ODMR resonances and spoiling magnetic field sensitivity. Optimal sensitivity therefore requires balancing the ODMR linewidth and contrast with the NV density in Eq. (10). Conditions for a favorable ratio of the two NV charge states, NV−/ NV0, are also needed to ameliorate the NV0 contribution to background fluorescence, which spoils the NV− contrast used for imaging .
The performance of diamonds with different C and N isotopes is an important consideration. The 15NV (spin-1/2 nucleus) is more favorable for QDM imaging because it gives greater ODMR contrast and requires a narrower range of MW probe frequencies than the more common 14NV (spin-1 nucleus). However, because 15N is the less abundant isotope, diamonds fabricated without special procedures for isotopic control will typically be dominated by 14N.
Magnetic inhomogeneity from 13C (spin-1/2) and paramagentic P1 defect centers limits the NV thus, isotopically purified 12C (I=0) diamonds are ideal , . For diamonds with a 1.1% natural abundance of 13C present, it is advantageous to increase the P1 density, resulting in larger NV density without contributing too much to the P1-limited . An NV layer fabricated in an isotopically enriched 12C layer can reduce the ODMR linewidth. However, this may be irrelevant for NVs shallower than ~10 nm because of magnetic inhomogeneity introduced by electrons on the diamond surface .
Synthetic diamond chips used in QDMs are available in several cuts. The most common are diamonds with the top face along the  plane and the sides along the  or  planes (Figure 14A). The NVs in these diamonds point roughly 35 ° out of the plane. Less common diamond cuts include  and  top faces. The former has two NV orientations in the plane, while the latter has one NV orientation pointing normal to the face. Other more exotic diamond cuts exist: for instance, Ref.  used a diamond with a  NV layer. The choice of diamond cut does not usually impact the QDM performance. However, different cuts of diamond have different availability and pricing because of the challenge of producing crystals that are not grown along diamond’s preferential growth axis. Surface termination effects can be of impact .
The impact of diamond characteristics on specific QDM techniques is summarized as follows:
For CW ODMR imaging, the laser and MW linewidth broadening should match the diamond (Eq. (10)).
For Ramsey imaging, the diamond limits the phase accumulation time.
For dynamical decoupling imaging, the diamond T2 limits the phase accumulation time (depending on the magnetic noise spectrum and pulse sequence).
For Rabi and T1 imaging, the diamond sets the spectral filter function. The intrinsic NV T1 depends on the NV density and depth.
A QDM typically uses a 532-nm solid-state laser for optical pumping because of its availability and performance. The green pump laser intensity is weak, typically ~10–1000 W/cm2, when illuminating a field of view of a few millimeters, which can be a limitation for pulsed readout techniques. The NV 3A2→3E optical transition spans hundreds of nanometers because of the phonon sideband as discussed in Section 2.2, which allows for laser excitation wavelengths ranging from 637 to ~470 nm . Past experiments have pumped NVs with 532-nm frequency-doubled Nd:YAG and Nd:YVO4 lasers, 637-nm and 520-nm diode lasers and LEDs, 594-nm HeNe lasers, argon-ion laser lines (457, 476, 488, 496, and 514 nm), and supercontinuum lasers with an acousto-optic tunable filter , , . There have been attempts to find the illumination wavelength with the most favorable cross-section and NV−/NV0 charge-state ratio . Since the NV readout measures a fluorescence intensity, fluorescence intensity instability from the laser or the optics must be minimized for the QDM magnetic sensitivity to reach the photon shot noise limit.
Increasing the illumination intensity improves the NV fluorescence intensity, the photon shot noise, and sometimes the ODMR lineshape. The 3A2→3E optical transition is dipole-allowed when illuminating with light polarized in the xy-plane of the NV coordinate system defined in Section 2.1 . Thus, in a projection magnetic microscopy experiment (Figure 14C), a laser polarization is chosen that favors the optical absorption selection rules for the selected NV orientation. If all NV orientations are interrogated, a laser polarization is selected that addresses all NV orientations with comparable strength. Increasing the laser illumination power increases diamond and sample heating on approximately linear scaling, while the photon shot noise limit only increases as the square root of the laser power. Furthermore, as the fractional photon shot noise improves, the analog-to-digital conversion bit depth must also improve to avoid being limited by quantization noise.
The available laser intensity affects the various QDM techniques in the following ways:
For CW ODMR imaging, varying the laser intensity affects the ODMR linewidth (Figure 10).
For pulsed imaging experiments, ideally the laser intensity should be close to optical saturation. Weaker laser intensity, longer tI, and longer tR will worsen the experiment time resolution.
5.3 Microwave source
The simplest way to apply a MW field to the NVs is with a piece of wire connected to a coaxial cable. The QDM MW field is ideally uniform across the NV layer’s field of view, and there are a variety of alternative engineered MW antennas that aim to optimize the MW field homogeneity, efficiency, or bandwidth , , , , , , . By the transition selection rules, the transitions between the 3A2 sublevels require left-circularly or right-circularly polarized MW . One QDM MW antenna option is a MW loop as shown in Figure 1; another option is a pair of crossed MW stripline resonators . The striplines are excited in phase (or 90° out of phase) to produce a linearly (or circularly) polarized MW field as needed for a given sensing protocol.
Increasing the MW power improves the contrast in a CW ODMR measurement but also increases the linewidth between 3A2 resonances. Figure 10 demonstrates that, for simulated CW ODMR measurements, optimizing QDM magnetic sensitivity implies tradeoffs of laser and MW power . Choosing an optical pumping rate much greater than the MW transition rate results in weak contrast, since the laser quickly repumps any NV population fraction removed by resonant MWs. Increasing the MW field amplitude improves the fluorescence contrast but also increases the ODMR linewidth. The MW intensity noise also affects the QDM sensitivity by influencing the ODMR contrast and linewidth in a manner similar to fluorescence intensity noise.
Selecting an appropriate MW frequency sweep rate is critical in CW ODMR measurements, where the probe MW frequency is swept across the NV resonance. The NV reaction time depends on the NV optical pumping rate and MW transition rate, and a sufficient response time is needed for the NVs to re-equilibrate to the updated conditions after updating in the MW probe frequency. This also applies to experiments using lock-in detection to combat fluorescence intensity noise: the MW modulation rate must be slower than the NV reaction time, typically set by the optically pumping rate .
When deciding how to apply the microwave field, some of the options affect the various QDM modalities differently:
For CW ODMR imaging, increasing the MW power broadens the ODMR line but also improves the contrast.
For Ramsey, Hahn echo, and dynamical decoupling imaging, spatial MW inhomogeneity and pulse errors can reduce the NV contrast and worsen the sensitivity.
5.4 Static magnetic field
The QDM bias magnetic field B0 can be provided by electromagnets (Helmholtz coil sets, solenoids, and C-frame/H-frame electromagnets) or permanent magnets , ,  as shown in Figure 1. Electromagnets allow the selection of any arbitrary |B0| up to a few tesla. However, they require a stable current supply, may need water cooling for the magnet, and can add to sample and system heating. Permanent magnets allow higher B0 in a more compact instrument, though the applied B0 can drift with temperature.
The choice of the bias field amplitude |B0| depends on the samples being measured. Soft magnetic samples that might have their magnetization changed by an applied magnetic field require |B0| to be minimized. This has the added benefit that small |B0| typically implies a small |B0| gradient across the field of view. A large |B0| can be beneficial when imaging paramagnetic minerals, since the magnetization from paramagnetic particles scales with |B0| until saturation . For Rabi imaging or T1 magnetometry, |B0| is chosen such that the NV spin transition frequency matches the AC sample frequency being interrogated , . Due to nitrogen nuclear polarization, |B0| ~30–50 mT improves the NV fluorescence contrast , . Finally, |B0|=0 is an intuitive choice for NV thermometry or electrometry experiments (Figure 14B).
The direction of B0 also factors into the specific QDM application . Alignment of B0 along the NV axis ( crystallographic direction), allow for interrogating the NVs along this direction, ignoring the other three NV directions, Figure 14C. This approach allows the optimization of the other measurement parameters, e.g. the optical polarization, to maximize the fluorescence and contrast from the selected NV orientation. Alternatively, the B0 magnitude and direction are chosen such that each NV orientation has different resonance frequencies and nonoverlapping spectra, Figure 14D. This approach allows the reconstruction of vector magnetic field information from the eight NV resonance frequencies. B0 can also be aligned with the crystallographic  or  directions, such that the resonance frequencies for different NV orientations are degenerate, leading to improved contrast. If B0 is aligned along the  direction with a diamond cut along , the magnetic field projection direction is normal to the chip, though the Zeeman shift is times weaker than for B0 along the  direction (Figure 14E). Finally, there may be some experiments where the choice of B0 is forced by the sample being tested. This could cause the NV ODMR lines to overlap, making it difficult to resolve the resonance frequencies and extract vector magnetic field information. This difficulty can be ameliorated using the 3A2→3E optical polarization selection rules to distinguish the light contributions from each NV orientation .
The B0 field should ideally be as uniform as possible. B0 inhomogeneity can cause the following problems:
In pulsed NV experiments, B0 inhomogeneity will cause spatially dependent pulse errors, which limit the NV contrast and sensitivity.
For all experiments, a B0 gradient on top of a the desired Bs is something that should be subtracted out. A uniform B0, allows for subtraction of a constant offset.
In an extreme case, B0 inhomogeneity can contribute to NV linewidth broadening within each pixel.
QDMs employ various ways to illuminate the NV layer with the pump laser light, depending on other experimental constraints. Side illumination of the diamond chip  is a good method for QDMs using a low-magnification (long working distance) objective with a large field of view since the beam will have enough space to avoid clipping the objective and also illuminate a large area. Another approach is to illuminate through the objective by focusing the pump laser at the back aperture to get parallel rays out of the objective . This method works better for QDMs operating with high-magnification microscopes. The laser polarization is easier to control, but focusing the laser at the objective back aperture can lead to burns. Techniques to avoid illuminating the sample as well as the NVs include illumination via total internal reflection in the diamond, shaping the pump laser beam into a light sheet using cylindrical lenses, or coating the NV surface with a reflective layer to reduce the optical intensity through the diamond chip , , .
Optimal photon collection efficiency requires the largest achievable NA for the microscope objective. In practice, the NA for a given magnification is limited, and high-NA objectives are often also high-magnification objectives with a short working distance (sometimes shorter than the diamond thickness). Imaging NV fluorescence through the diamond chip may cause optical aberrations that can spoil the image quality, though the authors are unaware of any QDM experiment that corrects for this. As with other optical microscopes, a QDM images a broadband NV fluorescence (~637–800 nm), so chromatic aberration in the microscope optics is also important to mitigate. Pulsed NV experiments commonly use an acousto-optic modulator (AOM) as an optical switch. For the AOM, the rise time, extinction ratio, and efficiency are the parameters to consider for a given application.
QDM camera selection for a targeted application requires consideration of the expected photon collection rate from the NV layer, camera read noise and dark current noise, well depth, global/rolling shutter capability, software and external triggering, frame rate, data transfer rate, pixel size, and quantum efficiency . For experiments with a high photon count rate, the camera must handle enough photoelectrons per second without saturating. Here, the pixel well depth, number of pixels, quantum efficiency, and frame rate are the important quantities to consider, because they determine the maximum photon count rate for fluorescence detection.
The camera frame rate can limit the experimentally realizable temporal resolution. Increasing the camera frame rate is possible by using only a fraction of the sensor. However, the resulting product between the frame rate and the number of pixels usually decreases, indicating that use of the full camera sensor is better for maximizing the number of photoelectrons per second. Alternatively, if the photon count rate is low, parameters such as the read noise and dark current noise should be minimized, while the quantum efficiency is maximized. For pulsed experiments, a low camera frame rate can throttle the experiment repetition rate and sensitivity.
The camera sensor size determines the microscope magnification for a desired field of view. The microscope spatial resolution can be set by the camera pixel size (rather than the optical diffraction limit) if the camera pixels are too widely spaced for the microscope magnification. The diffraction-limited spatial resolution should be oversampled by at least 2× to avoid having the pixel size spoil the diffraction-limited spatial resolution. A given choice of microscope optics has a finite effective image area, and the camera image can have darkened corners (vignetting) if the camera sensor area is too large.
As previously mentioned, the optical readout time tR must be balanced with the minimum camera exposure time and the maximum camera frame rate for pulsed QDM experiments. Specifically:
Sensitivity is lost for experiments with a measurement time tmeas faster than the camera frame rate because the camera is too slow to acquire a new frame at the rate it takes to do each experiment.
Experiments for which the minimum camera exposure time is longer than tR require the readout laser to be off for the duration of the time difference.
5.7 Diamond mounting and configuration
There are two primary ways to prepare the diamond sensor chip and the sample in the QDM. The first method is to fix the diamond chip in the optical microscope setup and move the sample independently with kinematic stages. This way, the diamond chip position (and all other optics) are permanent, keeping the relative positions of the optics, diamond location, and orientation, MW field, and magnetic field constant for all measurements to improve reproducibility and enable faster setup time for new samples. The second method is to mount the diamond chip directly on the sample, and then move the diamond and the sample together within the microscope field of view. This integrated diamond/sample approach offers more certainty that the NV-sample separation is minimized. Generally, sample mounting and manipulation in a QDM is easier with an upright microscope setup rather than with an inverted microscope.
5.8 General design considerations
Table 2 summarizes equipment parameters that optimize QDM operation. While some of the above specifications are technique- or application-specific, this table describes general design choices that affect all QDM instruments.
|Diamond||–||NV density affects the sensitivity|
|–||Inhomogeneity in strain, density, and magnetic environment spoils the sensitivity|
|–||Match the NV layer thickness and sample thickness|
|Laser||–||Field of view sets laser intensity|
|–||Laser intensity noise can limit sensitivity|
|–||Laser polarization addresses different NV orientations|
|–||Homogeneous illumination is desirable|
|Microwaves||–||Amplitude and phase instability affect sensitivity|
|–||Amplitude homogeneity is desirable|
|B0 field||–||A B0 gradient can cause varying sensitivity uniformity|
|Optics||–||Microscope objective sets the collection efficiency and optical diffraction limit|
|–||Microscope components set the magnification and field of view size|
|Camera||–||Pixel size should oversample other resolution limitations (e.g. optical diffraction)|
|–||Frame rate×well depth×number of pixels set the best possible sensitivity|
|–||Transfer rate and buffer size can throttle the maximum experiment rate|
|–||Camera efficiency is worse than photodiode efficiency|
|Diamond mount||–||Aim for high thermal and mechanical stability during an experiment|
|Magnetic environment||–||Mitigate background field (e.g. Earth’s field, electronics, etc.)|
6 QDM applications
QDM magnetic field imaging has been applied to a diverse range of applications across numerous fields of research. For every given application, the appropriate experimental protocol must be chosen for optimal performance, including the desired temporal resolution and magnetic frequency range. This subsequently dictates the QDM component implementation. Table 3 lists the application target areas for the various QDM techniques and the respective frequencies. To more easily motivate future unrealized QDM applications, the following sections includes examples of successful QDM applications for each frequency range and application area.
|Broadband 0–1 MHz||Narrowband ∼1 kHz–20 MHz||Narrowband 10 MHz–100 GHz|
|Applications||Paleomagnetism and rock magnetism|
Overview of techniques and potential applications for widefield magnetic imaging.
6.1 Broadband imaging of 0–1 MHz magnetic fields
CW ODMR imaging experiments of static magnetic fields is among the most successful QDM imaging applications to date. The relatively low MW and laser power requirements and simplicity of the experimental control allow for imaging of static magnetic fields up to a 4×4 mm region (limited by the size of a diamond substrate.) Most experiments up to this point have chosen to focus on large quasi-static magnetic fields because of the relatively loose requirements for performance of the QDM. Figure 15 shows several examples.
6.1.1 Current distributions
Imaging magnetic fields from 2D current distributions was among the first demonstrations of a QDM system because of the flexibility in choice of the magnetic field amplitude, temporal profile, and spatial structure, making it well suited for verifying the fidelity of magnetic field imaging experiments , . For sufficiently simple wire patterns, one can simulate the expected current distribution and field map to compare with and validate the QDM measurement. In these early experiments, the large sample current amplitudes were needed because of limited diamond sensitivity values at the time.
Magnetic field imaging for determining the current flow along circuit traces was one of the first demonstrations  of imaging the vector component from a nontrivial current distribution. However, the sensitivity of the system was not optimized, and nontrivial temporal dynamics of current flow were not investigated.
Magnetic field imaging can allow the interrogation of nontrivial current flow in 2D materials. Magnetic imaging has been utilized in probing the spatial dynamics of current flow in graphene. In an initial demonstration, current was passed through graphene, and defects in the 2D graphene were apparent as a result of the current flowing around them . In a separate experiment, magnetic imaging experiments were performed to probe the viscous Dirac fluid nature of current near the Dirac point. High-resolution magnetic field imaging allowed the direct measurement of the parabolic current profile associated with the hydrodynamic behavior of this Dirac fluid .
Looking to the future, the application space can dramatically improve if Ramsey imaging is implemented and optimized. With an optimized version of a Ramsey imaging system, there is projected to be sufficient volume-normalized sensitivity to image the propagation of activity-associated current in a mammalian neuron in real time .
6.1.2 Magnetic particles and domains
Measuring the DC component of magnetic particles and domains has yielded some of the most transformative applications of wide-field magnetic imaging to date. Examples in the literature span from magnetotactic bacteria  and magnetically labeled cells ,  to remanent magnetization in geological samples  and thin magnetic films . Success in these applications has been due in part to the generally static (enabling signal averaging) and large magnetic fields produced by these sources, which together reduce the need to push the state of the art on sensitivity.
In the earliest biological QDM experiment, the intrinsic magnetite inside magnetotactic bacteria was measured , as shown in Figure 15D. Other works have been performed to look at the intrinsic magnetite in chiton teeth to study iron mineralization  (see Figure 15G) and malarial hemozoin nanocrystals .
Magnetically labeling cells is a promising technique for tracking and identifying rare cell types , . Other groups have followed up on this work with extrinsic magnetic particles in applications relating to probing the origin of contrast agents in magnetic resonance imaging (MRI) , as shown in Figure 15H, and furthering the imaging resolution and sensitivity on magnetic particle imaging , .
QDMs have proven to be a valuable component in the toolbox of remanent magnetization studies in geological samples (see Figure 15F) . Initial demonstrations  were performed on the Semarkona meteorite to assist in determining the magnetic field present during planetary formation. Follow-up work with QDMs have demonstrated their utility in imaging magnetization carriers at the grain scale. Recent example applications have included the imaging of large magnetite grains to visualize multidomain structure  and of zircons , ,  to understand and constrain the history of Earth’s dynamo. The full potential of QDM as a rock magnetic instrument is only beginning to be explored, with experiments on terrestrial and extraterrestrial rock types being pursued.
QDMs have extended their range to condensed matter to probe thin magnetic films such as magnetic memory (see Figure 15E)  and explore questions related to the origins and properties of vortices in superconductors .
6.2 Narrowband imaging of ~1 kHz–20 MHz magnetic fields
Narrowband magnetic imaging in an intermediate frequency range is mostly applicable for imaging magnetic fields originating from current distributions and the magnetic field from precessing spins in nuclear magnetic resonance (NMR) applications (see Figure 16). Much of the early work in NV magnetic imaging pushed the state of the art in these regimes, but more development is needed to explore the full range of applications.
6.2.1 Current distributions
Similarly to the broadband case, current distributions were initially used to validate the fidelity and effectiveness of AC magnetometry pulse sequences in an imaging modality . In this demonstration, current with frequencies ranging from 4 kHz to 100 kHz were sent through wires fabricated on the diamond.
One promising application of this technique is imaging magnetic fields that oscillate near the clock frequency of circuits for side channel attack analysis . NV diamond can allow for the ability to include spatial information. Extending the sensing frequency beyond ~20 MHz is challenging for several reasons. For sensing high frequencies, the MW π pulse duration should be short compared to the period of the sensing signal. Short π pulses require strong MW fields to achieve high Rabi frequencies. Strong, uniform MW pulses over a large area are a difficult engineering challenge requiring more sophisticated MW antenna design. Furthermore, even if these requirements are fulfilled, the strong MW fields can interfere with or damage the sample being sensed.
6.2.2 NMR signals
NMR spectroscopy allows label-free detection and quantification of molecules with excellent chemical specifity. The use of narrowband AC magnetic imaging techniques to record local NMR signals in individual QDM pixels opens the possibility of highly multiplexed 2D spatial density mapping of arbitrary molecular species. Potential applications include imaging small-molecule concentrations in neuronal slice preparations or bacterial biofilms , spatially resolved battery electrochemistry , detection and determination of the chemical composition of proteins , or possibly a readout for molecular data storage .
Even without the high spectral resolution required to distinguish molecular species (typically 1 ppm of the nuclear Larmor frequency or better, which places stringent technical constraints on the magnitude, stability, and homogeneity of the bias magnetic field B0), the combination of QDM imaging with correlation spectroscopy techniques  and/or strong pulsed magnetic gradients  can provide spatial maps of physical properties of samples, such as fluid density, net flow velocity fields, and/or local diffusion rates . This could have applications in the study of porous media in petrochemistry, filtration, or catalysis.
As with broadband Ramsey spectroscopy, the pulse sequences used for narrowband AC magnetic imaging necessitate efficient temporal segmentation of NV fluorescence data at fast timescales, which is challenging for standard scientific imaging cameras. For this reason, there have been few reported demonstrations of NMR signal imaging using QDMs reported in the literature to date , , and none with the spectral resolution required to distinguish molecular species. Nevertheless, we anticipate that ongoing work to integrate broadband Ramsey spectroscopy into the QDM platform can be directly extended to narrowband AC signal detection and ultimately to high spectral-resolution NMR readout techniques , , , .
6.3 Narrowband imaging of 10 MHz–100 GHz magnetic fields
6.3.1 Microwave imaging
QDM imaging of the microwave field from wires, resonators, and structures is possible by measuring the Rabi frequency in a pulsed experiment , or by using the fluorescence contrast in a CW ODMR experiment  (see Figure 17A and B). An initial step is to compare the NV measurement to a predicted magnetic field map from a finite element method (FEM) calculation. One goal is to use NV microwave imaging as a tool to validate that the FEM or the fabrication is what is expected for more nontrivial devices such as atom chips.
6.3.2 T1 imaging for paramagnetic spins
Just as coherent resonant microwaves drive transitions between the NV 3A2 sublevels, external paramagnetic spins can have the same effect, spoiling the NV T1. This can be used for imaging, as shown in Figure 17C and D. Paramagnetic spins with short T1 can produce broadband magnetic noise that spoils the NV T1, while long-lived paramagnetic spins can spoil the NV T1 for particular |B0| where there is level-crossing between the NVs and the external spins. Previous experiments have examined NV T1 relaxation due to external paramagnetic spins often used as MRI contrast agents (e.g. Gd3+, Mn2+), Cu2+ and iron ions in ferritin. The motivation was to monitor the concentration in a microfluidic device over time , , . Analyzing NV T1 data as a function of |B0| can generate the magnetic noise spectrum, identify specific paramagnetic species, and yield the paramagnetic concentrations. Further work may investigate imaging paramagnetic spins using double electron-electron resonance (DEER) .
7 Conclusion and outlook
In recent years, QDM has addressed important scientific questions in diverse fields, which further motivates interest in this technology. Fortunately, a QDM is relatively straightforward to build, and the technology is sufficiently mature that running a QDM experiment from start to finish is straightforward. As diamond characteristics and NV sensing techniques improve, a growing range of QDM capabilities and applications can be expected, including in extreme environments, e.g. high-pressure, high-temperature, and cryogenic , , .
QDM imaging of magnetic fields is well established with a rapidly expanding application space. The sensitivity of NVs to temperature distributions and electric fields indicate that QDMs should also be applicable to imaging temperature and electric field. However, imaging temperature inhomogeneity is challenging since temperature gradients dissipate quickly at micrometer length scales in most materials, and an in-contact diamond chip will accentuate the heat dissipation from the sample being tested due to the excellent thermal conductivity of diamond, thereby modifying the temperature profile being measured. Compared to magnetic sensing, electric sensing has the drawback that generally the electric susceptibility of a material is larger than the magnetic susceptibility, meaning that materials are often effectively transparent to magnetic fields while screening or significantly modifying electric fields. Nonetheless, QDM electric field imaging is an exciting direction that is largely unexplored. Finally, NV imaging of stress within a diamond chip (which was previously done with a single NV in an atomic force microscopy setup) is now being pursued in widefield experiments , . NV stress measurements can provide information about internal and external tensile and shear stress felt by the NVs, and could eventually be used to image pressure or to measure nuclear recoil tracks for particle physics experiments .
We acknowledge support for this work from the MITRE Corporation through the MITRE Innovation Program (MIP) and the MITRE Corporation Research Agreement No. 124787 with Harvard University. P.K. acknowledges support from the Sandia National Laboratories (SNL) Truman Fellowship Program, funded by the Laboratory Directed Research and Development (LDRD) Program at SNL, operated by National Technology and Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International, Inc., for the U.S. NNSA under contract DE-NA0003525.
 Glenn DR, Fu RR, Kehayias P, et al. Micrometer-scale magnetic imaging of geological samples using a quantum diamond microscope. Geochem Geophys Geosyst 2017;18:3254–67.10.1002/2017GC006946 Search in Google Scholar
 Schirhagl R, Chang K, Loretz M, Degen CL. Nitrogen-vacancy centers in diamond: nanoscale sensors for physics and biology. Ann Rev Phys Chem 2014;65:83–105.10.1146/annurev-physchem-040513-103659 Search in Google Scholar
 Rondin L, Tetienne J-P, Spinicelli P, et al. Nanoscale magnetic field mapping with a single spin scanning probe magnetometer. Appl Phys Lett 2012;100:153118.10.1063/1.3703128 Search in Google Scholar
 Maletinsky P, Hong S, Grinolds MS, et al. A robust scanning diamond sensor for nanoscale imaging with single nitrogen-vacancy centres. Nat Nanotechnol 2012;7:320–4.2250470810.1038/nnano.2012.50 Search in Google Scholar
 Fong LE, Holzer JR, McBride KK, Lima EA, Baudenbacher F, Radparvar M. High-resolution room-temperature sample scanning superconducting quantum interference device microscope configurable for geological and biomagnetic applications. Rev Sci Instrum 2005;76:053703.10.1063/1.1884025 Search in Google Scholar
 Allwood DA, Xiong G, Cooke MD, Cowburn RP. Magneto-optical kerr effect analysis of magnetic nanostructures. J Phys D: Appl Phys 2003;36:2175–82.10.1088/0022-3727/36/18/001 Search in Google Scholar
 Doherty MW, Manson NB, Delaney P, Jelezko F, Wrachtrup J, Hollenberg LCL. The nitrogen-vacancy colour centre in diamond. Phys Rep 2013;528:1–45.10.1016/j.physrep.2013.02.001 Search in Google Scholar
 Rondin L, Tetienne J-P, Hingant T, Roch J-F, Maletinsky P, Jacques V. Magnetometry with nitrogen-vacancy defects in diamond. Rep Prog Phys 2014;77:056503.2480149410.1088/0034-4885/77/5/056503 Search in Google Scholar
 Bauch E, Hart CA, Schloss JM, et al. Ultralong dephasing times in solid-state spin ensembles via quantum control. Phys Rev X 2018;8:031025. Search in Google Scholar
 Shin CS, Butler MC, Wang H-J, et al. Optically detected nuclear quadrupolar interaction of 14N in nitrogen-vacancy centers in diamond. Phys Rev B 2014;89:205202.10.1103/PhysRevB.89.205202 Search in Google Scholar
 Acosta VM, Bauch E, Ledbetter MP, Waxman A, Bouchard L-S, Budker D. Temperature dependence of the nitrogen-vacancy magnetic resonance in diamond. Phys Rev Lett 2010;104:070801.2036686810.1103/PhysRevLett.104.070801 Search in Google Scholar
 Felton S, Edmonds AM, Newton ME, et al. Hyperfine interaction in the ground state of the negatively charged nitrogen vacancy center in diamond. Phys Rev B 2009;79:075203.10.1103/PhysRevB.79.075203 Search in Google Scholar
 Udvarhelyi P, Shkolnikov VO, Gali A, Burkard G, Pályi A. Spin-strain interaction in nitrogen-vacancy centers in diamond. Phys Rev B 2018;98:075201.10.1103/PhysRevB.98.075201 Search in Google Scholar
 Barfuss A, Kasperczyk M, Kölbl J, Maletinsky P. Spin-stress and spin-strain coupling in diamond-based hybrid spin oscillator systems. Phys Rev B 2019;88:174102. Search in Google Scholar
 Jamonneau P, Lesik M, Tetienne JP, et al. Competition between electric field and magnetic field noise in the decoherence of a single spin in diamond. Phys Rev B 2016;93:024305.10.1103/PhysRevB.93.024305 Search in Google Scholar
 Kehayias P, Turner MJ, Trubko R, et al. Imaging crystal stress in diamond using ensembles of nitrogen-vacancy centers. 2019;arXiv:1908.09904. Search in Google Scholar
 Rogers LJ, McMurtrie RL, Sellars MJ, Manson NB. Time-averaging within the excited state of the nitrogen-vacancy centre in diamond. New J Phys 2009;11:063007.10.1088/1367-2630/11/6/063007 Search in Google Scholar
 Batalov A, Jacques V, Kaiser F, et al. Low temperature studies of the excited-state structure of negatively charged nitrogen-vacancy color centers in diamond. Phys Rev Lett 2009;102:195506.1951897310.1103/PhysRevLett.102.195506 Search in Google Scholar
 Henderson B, Imbusch GF. Optical spectroscopy of inorganic solids. New York, Oxford University Press, 2006. Search in Google Scholar
 Robledo L, Bernien H, van der Sar T, Hanson R. Spin dynamics in the optical cycle of single nitrogen-vacancy centres in diamond. New J Phys 2011;13:025013.10.1088/1367-2630/13/2/025013 Search in Google Scholar
 Tetienne J-P, Rondin L, Spinicelli P, et al. Magnetic-field-dependent photodynamics of single NV defects in diamond: an application to qualitative all-optical magnetic imaging. New J Phys 2012;14:103033.10.1088/1367-2630/14/10/103033 Search in Google Scholar
 Harrison J, Sellars MJ, Manson NB. Measurement of the optically induced spin polarisation of NV centres in diamond. Diam Relat Mater 2006;15:586–8.10.1016/j.diamond.2005.12.027 Search in Google Scholar
 Rogers LJ, Armstrong S, Sellars MJ, Manson NB. Infrared emission of the NV centre in diamond: zeeman and uniaxial stress studies. New J Phys 2008;10:103024.10.1088/1367-2630/10/10/103024 Search in Google Scholar
 Wee T-L, Tzeng Y-K, Han C-C, et al. Two-photon excited fluorescence of nitrogen-vacancy centers in proton-irradiated type ib diamond. J Phys Chem A 2007;111:9379–86.10.1021/jp073938o17705460 Search in Google Scholar
 Chapman R, Plakhotnik T. Quantitative luminescence microscopy on nitrogen-vacancy centres in diamond: saturation effects under pulsed excitation. Chem Phys Lett 2011;507:190–4.10.1016/j.cplett.2011.03.057 Search in Google Scholar
 Peppers J, Martyshkin DV, Mirov SB, et al. Laser spectroscopy of highly doped NV− centers in diamond. Proc. SPIE 10511, Solid State Lasers XXVII: Technology and Devices, 2018. Search in Google Scholar
 Dréau A, Lesik M, Rondin L, et al. Avoiding power broadening in optically detected magnetic resonance of single NV defects for enhanced dc magnetic field sensitivity. Phys Rev B 2011;84:195204.10.1103/PhysRevB.84.195204 Search in Google Scholar
 Gruber A, Dräbenstedt A, Tietz C, Fleury L, Wrachtrup J, von Borczyskowski C. Scanning confocal optical microscopy and magnetic resonance on single defect centers. Science 1997;276:2012–4.10.1126/science.276.5321.2012 Search in Google Scholar
 Lima EA, Weiss BP. Ultra-high sensitivity moment magnetometry of geological samples using magnetic microscopy. Geochem Geophys Geosyst 2016;17:3754–74.10.1002/2016GC006487 Search in Google Scholar
 Baratchart L, Chevillard S, Hardin DP, Leblond J, Lima EA, Marmorat J-P. Magnetic moment estimation and bounded extremal problems. Inverse Probl Imaging 2019;13:29. Search in Google Scholar
 Glenn DR, Bucher DB, Lee J, Lukin MD, Park H, Walsworth RL. High-resolution magnetic resonance spectroscopy using a solid-state spin sensor. Nature 2018;555:351–4.10.1038/nature2578129542693 Search in Google Scholar
 Mittiga T, Hsieh S, Zu C, et al. Imaging the local charge environment of nitrogen-vacancy centers in diamond. Phys Rev Lett 2018;121:246402.3060873210.1103/PhysRevLett.121.246402 Search in Google Scholar
 Slichter CP. Principles of magnetic resonance. In: Springer series in solid-state sciences. Berlin, Heidelberg: Springer, 1996. Search in Google Scholar
 Barry JF, Schloss JM, Bauch E, et al. Sensitivity optimization for NV-diamond magnetometry. 2019;arXiv:1903.08176. Search in Google Scholar
 Shin CS, Avalos CE, Butler MC, et al. Room-temperature operation of a radiofrequency diamond magnetometer near the shot-noise limit. J Appl Phys 2012;112:124519.10.1063/1.4771924 Search in Google Scholar
 Barry JF, Turner MJ, Schloss JM, et al. Optical magnetic detection of single-neuron action potentials using quantum defects in diamond. Proc Natl Acad Sci 2016;113:14133–8.10.1073/pnas.1601513113 Search in Google Scholar
 Schloss JM, Barry JF, Turner MJ, Walsworth RL. Simultaneous broadband vector magnetometry using solid-state spins. Phys Rev Appl 2018; 10:034044.10.1103/PhysRevApplied.10.034044 Search in Google Scholar
 de Lange G, Wang ZH, Ristè D, Dobrovitski VV, Hanson R. Universal dynamical decoupling of a single solid-state spin from a spin bath. Science 2010;330:60–3.2082945210.1126/science.1192739 Search in Google Scholar
 Bauch E, Singh S, Lee J, et al. Decoherence of dipolar spin ensembles in diamond. 2019;arXiv:1904.08763. Search in Google Scholar
 de Sousa R. Electron spin as a spectrometer of nuclear-spin noise and other fluctuations. In: Fanciulli M, editor. Electron spin resonance and related phenomena in low-dimensional structures, volume 115 of topics in applied physics. Berlin, Heidelberg: Springer, 2009:183–220. Search in Google Scholar
 Szankowski P, Ramon G, Krzywda J, Kwiatkowski D, Cywinski L. Environmental noise spectroscopy with qubits subjected to dynamical decoupling. J Phys Condens Matter 2017;29:333001.2856923910.1088/1361-648X/aa7648 Search in Google Scholar
 Shao L, Liu R, Zhang M, et al. Wide-field optical microscopy of microwave fields using nitrogen-vacancy centers in diamonds. Adv Opt Mater 2016;19:E2347. Search in Google Scholar
 Horsley A, Appel P, Wolters J, et al. Microwave device characterization using a widefield diamond microscope. Phys Rev Appl 2018;10:044039.10.1103/PhysRevApplied.10.044039 Search in Google Scholar
 Allred JC, Lyman RN, Kornack TW, Romalis MV. High-sensitivity atomic magnetometer unaffected by spin-exchange relaxation. Phys Rev Lett 2002;89:130801.1222501310.1103/PhysRevLett.89.130801 Search in Google Scholar
 Dobrovitski VV, Feiguin AE, Hanson R, Awschalom DD. Decay of rabi oscillations by dipolar-coupled dynamical spin environments. Phys Rev Lett 2009;102:237601.1965897310.1103/PhysRevLett.102.237601 Search in Google Scholar
 Hämäläinen M, Hari R, Ilmoniemi RJ, Knuutila Ju, Lounasmaa OV. Magnetoencephalography – theory, instrumentation, and applications to noninvasive studies of the working human brain. Rev Mod Phys 1993;65:413–97.10.1103/RevModPhys.65.413 Search in Google Scholar
 Lima EA, Weiss BP. Obtaining vector magnetic field maps from single-component measurements of geological samples. J Geophys Res: Solid Earth 2009;114. Search in Google Scholar
 Tetienne J-P, Broadway DA, Lillie SE, et al. Proximity-induced artefacts in magnetic imaging with nitrogen-vacancy ensembles in diamond. Sensors 2018;18.29690603 Search in Google Scholar
 Naydenov B, Reinhard F, Lämmle A, et al. Increasing the coherence time of single electron spins in diamond by high temperature annealing. Appl Phys Lett 2010;97:242511.10.1063/1.3527975 Search in Google Scholar
 Pezzagna S, Naydenov B, Jelezko F, Wrachtrup J, Meijer J. Creation efficiency of nitrogen-vacancy centres in diamond. New J Phys 2010;12:065017.10.1088/1367-2630/12/6/065017 Search in Google Scholar
 Ozawa H, Tahara K, Ishiwata H, Hatano M, Iwasaki T. Formation of perfectly aligned nitrogen-vacancy-center ensembles in chemical-vapor-deposition-grown diamond (111). Appl Phys Express 2017;10:045501.10.7567/APEX.10.045501 Search in Google Scholar
 Loretz M, Pezzagna S, Meijer J, Degen CL. Nanoscale nuclear magnetic resonance with a 1.9-nm-deep nitrogen-vacancy sensor. Appl Phys Lett 2014;104:033102.10.1063/1.4862749 Search in Google Scholar
 Acosta VM, Bauch E, Ledbetter MP, et al. Diamonds with a high density of nitrogen-vacancy centers for magnetometry applications. Phys Rev B 2009;80:115202.10.1103/PhysRevB.80.115202 Search in Google Scholar
 Smits J, Damron J, Kehayias P, et al. Two-dimensional nuclear magnetic resonance spectroscopy with a microfluidic diamond quantum sensor. Sci Adv 2019;5:eaaw7895.3136076910.1126/sciadv.aaw7895 Search in Google Scholar
 Kleinsasser EE, Stanfield MM, Banks JKQ, et al. High density nitrogen-vacancy sensing surface created via He+ ion implantation of 12C diamond. Appl Phys Lett 2016;108:202401.10.1063/1.4949357 Search in Google Scholar
 Aude Craik DPL, Kehayias P, Greenspon AS, et al. A microwave-assisted spectroscopy technique for determining charge state in nitrogen-vacancy ensembles in diamond. 2018;arXiv:1811.01972. Search in Google Scholar
 Grinolds MS, Warner M, De Greve K, et al. Subnanometre resolution in three-dimensional magnetic resonance imaging of individual dark spins. Nat Nanotechnol 2014;9: 279–84.10.1038/nnano.2014.3024658168 Search in Google Scholar
 Rosskopf T, Dussaux A, Ohashi K, et al. Investigation of surface magnetic noise by shallow spins in diamond. Phys Rev Lett 2014;112:147602.2476601510.1103/PhysRevLett.112.147602 Search in Google Scholar
 Kilin SYa, Nizovtsev AP, Maevskaya TM, Dräbenstedt A, Wrachtrup J. Spectroscopy on single n-v defect centers in diamond: tunneling of nitrogen atoms into vacancies and fluorescence spectra. J Lumin 2000;86:201–6.10.1016/S0022-2313(00)00164-2 Search in Google Scholar
 Aslam N, Waldherr G, Neumann P, Jelezko F, Wrachtrup J. Photo-induced ionization dynamics of the nitrogen vacancy defect in diamond investigated by single-shot charge state detection. New J Phys 2013;15:013064.10.1088/1367-2630/15/1/013064 Search in Google Scholar
 Bauch E. Nitrogen-vacancy defects in diamond for sub-millimeter magnetometry. Master’s thesis. Berlin: TU Berlin, 2010. Search in Google Scholar
 Stürner FM, Brenneis A, Kassel J, et al. Compact integrated magnetometer based on nitrogen-vacancy centres in diamond. Diam Relat Mater 2019;93:59–65.10.1016/j.diamond.2019.01.008 Search in Google Scholar
 Beha K, Batalov A, Manson NB, Bratschitsch R, Leitenstorfer A. Optimum photoluminescence excitation and recharging cycle of single nitrogen-vacancy centers in ultrapure diamond. Phys Rev Lett 2012;109:097404.2300288310.1103/PhysRevLett.109.097404 Search in Google Scholar
 Bayat K, Choy J, Farrokh Baroughi M, Meesala S, Loncar M. Efficient, uniform, and large area microwave magnetic coupling to NV centers in diamond using double split-ring resonators. Nano Lett 2014;14:1208–13.10.1021/nl404072s24571706 Search in Google Scholar
 Zhang N, Zhang C, Xu L, et al. Microwave magnetic field coupling with nitrogen-vacancy center ensembles in diamond with high homogeneity. Appl Magn Reson 2016;47:589–99.10.1007/s00723-016-0777-5 Search in Google Scholar
 Labanowski D, Bhallamudi VP, Guo Q, et al. Voltage-driven, local, and efficient excitation of nitrogen-vacancy centers in diamond. Sci Adv 2018;4.30202783 Search in Google Scholar
 Jia W, Shi Z, Qin X, Rong X, Du J. Ultra-broadband coplanar waveguide for optically detected magnetic resonance of nitrogen-vacancy centers in diamond. Rev Sci Instrum 2018;89:064705.2996052710.1063/1.5028335 Search in Google Scholar
 Sasaki K, Monnai Y, Saijo S, et al. Broadband, large-area microwave antenna for optically detected magnetic resonance of nitrogen-vacancy centers in diamond. Rev Sci Instrum 2016;87:053904.2725043910.1063/1.4952418 Search in Google Scholar
 Eisenach ER, Barry JF, Pham LM, Rojas RG, Englund DR, Braje DA. Broadband loop gap resonator for nitrogen vacancy centers in diamond. Rev Sci Instrum 2018;89:094705.3027872410.1063/1.5037465 Search in Google Scholar
 Mrozek M, Mlynarczyk J, Rudnicki DS, Gawlik W. Circularly polarized microwaves for magnetic resonance study in the ghz range: application to nitrogen-vacancy in diamonds. Appl Phys Lett 2015;107:013505.10.1063/1.4923252 Search in Google Scholar
 Alegre TPM, Santori C, Medeiros-Ribeiro G, Beausoleil RG. Polarization-selective excitation of nitrogen vacancy centers in diamond. Phys Rev B 2007;76:165205.10.1103/PhysRevB.76.165205 Search in Google Scholar
 Hall LT, Kehayias P, Simpson DA, et al. Detection of nanoscale electron spin resonance spectra demonstrated using nitrogen-vacancy centre probes in diamond. Nat Commun 2016;7:10211.10.1038/ncomms1021126728001 Search in Google Scholar
 Fescenko I, Laraoui A, Smits J, et al. Diamond magnetic microscopy of malarial hemozoin nanocrystals. Phys Rev Appl 2019;11:034029.10.1103/PhysRevApplied.11.03402931245433 Search in Google Scholar
 Steiner M, Neumann P, Beck J, Jelezko F, Wrachtrup J. Universal enhancement of the optical readout fidelity of single electron spins at nitrogen-vacancy centers in diamond. Phys Rev B 2010;81:035205.10.1103/PhysRevB.81.035205 Search in Google Scholar
 Farchi E, Ebert Y, Farfurnik D, Haim G, Shaar R, Bar-Gill N. Quantitative vectorial magnetic imaging of multi-domain rock forming minerals using nitrogen-vacancy centers in diamond. SPIN 2017;07:1740015.10.1142/S201032471740015X Search in Google Scholar
 Wojciechowski AM, Karadas M, Huck A, et al. Contributed review: camera-limits for wide-field magnetic resonance imaging with a nitrogen-vacancy spin sensor. Rev Sci Instrum 2018;89:031501.10.1063/1.501028229604724 Search in Google Scholar
 Tetienne J-P, Dontschuk N, Broadway DA, Stacey A, Simpson DA, Hollenberg LCL. Quantum imaging of current flow in graphene. Sci Adv 2017;3:e1602429.10.1126/sciadv.160242928508040 Search in Google Scholar
 Ku MJH, Zhou TX, Li Q, et al. Imaging viscous flow of the dirac fluid in graphene using a quantum spin magnetometer. 2019;arXiv:1905.10791. Search in Google Scholar
 McCoey JM, Matsuoka M, de Gille RW, et al. Quantum magnetic imaging of iron biomineralisation in teeth of the chiton Acanthopleura hirtosa. 2019;arXiv:1902.09637. Search in Google Scholar
 Davis HC, Ramesh P, Bhatnagar A, et al. Mapping the microscale origins of mri contrast with subcellular NV diamond magnetometry. Nat Commun 2018;9:131.10.1038/s41467-017-02471-7 Search in Google Scholar
 Nowodzinski A, Chipaux M, Toraille L, Jacques V, Roch JF, Debuisschert T. Nitrogen-vacancy centers in diamond for current imaging at the redistributive layer level of integrated circuits. Microelectron Reliab 2015;55:1549–53.10.1016/j.microrel.2015.06.069 Search in Google Scholar
 Gould M, Barbour RJ, Thomas N, Arami H, Krishnan KM, Fu K-MC. Room-temperature detection of a single 19nm super-paramagnetic nanoparticle with an imaging magnetometer. Appl Phys Lett 2014;105:072406.10.1063/1.4893602 Search in Google Scholar
 Toraille L, Aïzel K, Balloul É, et al. Optical magnetometry of single biocompatible micromagnets for quantitative magnetogenetic and magnetomechanical assays. Nano Lett 2018;18:7635–41.10.1021/acs.nanolett.8b0322230380877 Search in Google Scholar
 Fu RR, Weiss BP, Lima EA, et al. Evaluating the paleomagnetic potential of single zircon crystals using the bishop tuff. Earth Planet Sci Lett 2017;458:1–13.10.1016/j.epsl.2016.09.038 Search in Google Scholar
 Tang F, Taylor RJM, Einsle JF, et al. Secondary magnetite in ancient zircon precludes analysis of a hadean geodynamo. Proc Natl Acad Sci 2019;116:407–12.10.1073/pnas.1811074116 Search in Google Scholar
 Weiss BP, Fu R, Einsle J, et al. Secondary magnetic inclusions in detrital zircons from the jack hills, western australia and implications for the origin of the geodynamo. Geology 2018;46:427–30.10.1130/G39938.1 Search in Google Scholar
 Schlussel Y, Lenz T, Rohner D, et al. Wide-field imaging of superconductor vortices with electron spins in diamond. Phys Rev Appl 2018;10:034032.10.1103/PhysRevApplied.10.034032 Search in Google Scholar
 DeVience SJ, Pham LM, Lovchinsky I, et al. Nanoscale NMR spectroscopy and imaging of multiple nuclear species. Nat Nanotechnol 2015;10:129–34.10.1038/nnano.2014.31325559712 Search in Google Scholar
 Ziem F, Garsi M, Fedder H, Wrachtrup J. Quantitative nanoscale MRI with a wide field of view. 2018;arXiv:1807.08343. Search in Google Scholar
 Agrawal D, Archambeault B, Rao JR, Rohatgi P. The EM side-channel(s). In: Kaliski BS, Koç CK, Paar C, editors. Cryptographic hardware and embedded systems – CHES 2002. Berlin, Heidelberg: Springer, 2003:29–45. Search in Google Scholar
 Teschler JK, Zamorano-Sánchez D, Utada AS, et al. Living in the matrix: assembly and control of vibrio cholerae biofilms. Nat Rev Microbiol 2015;13:255–68.10.1038/nrmicro343325895940 Search in Google Scholar
 Lovchinsky I, Sushkov AO, Urbach E, et al. Nuclear magnetic resonance detection and spectroscopy of single proteins using quantum logic. Science 2016;351:836–41.2684754410.1126/science.aad8022 Search in Google Scholar
 Cafferty BJ, Ten AS, Fink MJ, et al. Storage of information using small organic molecules. ACS Cent Sci 2019;5:911–6.31139727 Search in Google Scholar
 Laraoui A, Dolde F, Burk C, Reinhard F, Wrachtrup J, Meriles CA. High-resolution correlation spectroscopy of 13c spins near a nitrogen-vacancy centre in diamond. Nat Commun 2013;4:1651.10.1038/ncomms2685 Search in Google Scholar
 Arai K, Belthangady C, Zhang H, et al. Fourier magnetic imaging with nanoscale resolution and compressed sensing speed-up using electronic spins in diamond. Nat Nanotechnol 2015;10:859–64.10.1038/nnano.2015.17126258549 Search in Google Scholar
 Staudacher T, Raatz N, Pezzagna S, et al. Probing molecular dynamics at the nanoscale via an individual paramagnetic centre. Nat Commun 2015;6:8527.2645601710.1038/ncomms9527 Search in Google Scholar
 Schmitt S, Gefen T, Stürner FM, et al. Submillihertz magnetic spectroscopy performed with a nanoscale quantum sensor. Science 2017;356:832–7.10.1126/science.aam553228546208 Search in Google Scholar
 Mizuno K, Nakajima M, Ishiwata H, Masuyama Y, Iwasaki T, Hatano M. Wide-field diamond magnetometry with millihertz frequency resolution and nanotesla sensitivity. AIP Advances 2018;8:125316.10.1063/1.5048265 Search in Google Scholar
 Simpson DA, Ryan RG, Hall LT, et al. Electron paramagnetic resonance microscopy using spins in diamond under ambient conditions. Nat Commun 2017;8:458.10.1038/s41467-017-00466-y28878240 Search in Google Scholar
 Ziem FC, Götz NS, Zappe A, Steinert S, Wrachtrup J. Highly sensitive detection of physiological spins in a microfluidic device. Nano Lett 2013;13:4093–8.10.1021/nl401522a Search in Google Scholar
 Lesik M, Plisson T, Toraille L, et al. Magnetic measurements on micron-size samples under high pressure using designed NV centers. 2018;arXiv:1812.09894. Search in Google Scholar
 Hsieh S, Bhattacharyya P, Zu C, et al. Imaging stress and magnetism at high pressures using a nanoscale quantum sensor. 2018;arXiv:1812.08796. Search in Google Scholar
 Broadway DA, Johnson BC, Barson MSJ, et al. Microscopic imaging of elastic deformation in diamond via in-situ stress tensor sensors. 2018;arXiv:1812.01152. Search in Google Scholar
 Rajendran S, Zobrist N, Sushkov AO, Walsworth RL, Lukin M. A method for directional detection of dark matter using spectroscopy of crystal defects. Phys Rev D 2017;96:035009.10.1103/PhysRevD.96.035009 Search in Google Scholar
©2019 Edlyn V. Levine et al., published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 Public License.