# Patent

## Title
High sensitivity solid state magnetometer

## Abstract

A magnetometer for sensing a magnetic field may include a solid state electronic spin system, and a detector. The solid state electronic spin system may contain one or more electronic spins that are disposed within a solid state lattice, for example NV centers in diamond. The electronic spins may be configured to receive optical excitation radiation and to align with the magnetic field in response thereto. The electronic spins may be further induced to precess about the magnetic field to be sensed, in response to an external control such as an RF field, the frequency of the spin precession being linearly related to the magnetic field by the Zeeman shift of the electronic spin energy levels. The detector may be configured to detect output optical radiation from the electronic spin, so as to determine the Zeeman shift and thus the magnetic field.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application is based upon, and claims the benefit of priority under 35 U.S.C. §119(e) from U.S. Provisional Patent Application Ser. No. 60/991,872 (the “'872 provisional application”), filed Dec. 3, 2007, entitled “High Sensitivity Solid State Magnetometer”; and from U.S. Provisional Patent Application Ser. No. 61/058,366 (the “'366 provisional application”), filed Jun. 3, 2008, entitled “High Sensitivity Solid State Magnetometer.” The contents of the '872 provisional application and the '366 provisional application are incorporated herein by reference in their entireties as though fully set forth.
This application is related to PCT Application No. PCT/US2008/085428, filed concurrently herewith, entitled “Electronic Spin Based Enhancement of Magnetometer Sensitivity,” the content of which is hereby incorporated by reference in its entirety.


BACKGROUND
The detection of weak magnetic fields with high spatial resolution is an important problem in diverse areas ranging from fundamental physics and material science to data storage and biomedical science.
Over the past few decades, a wide variety of magnetic sensors have been developed using approaches that include, but are not limited to, SQUIDS (superconducting quantum interference devices), atomic vapor-based magnetometers, and magnetic resonance force microscopy. Even state-of-the-art systems have great difficulty, however, in detecting weak magnetic fields in small regions of space and under ambient environmental conditions, such as temperature.
There is a need for systems and methods for detecting magnetic fields with both higher precision and spatial resolution.
SUMMARY
A magnetometer for sensing a magnetic field may include a solid state electronic spin system, and a detector. The solid state electronic spin system may contain at least one electronic spin that is disposed within a solid state lattice, for example an NV (Nitrogen vacancy) center in diamond. The system may also contain many such electronic spins. The electronic spin or spins may be substantially free of interaction with the solid state lattice. The electronic spin or spins may be configured to align with the magnetic field in response to optical radiation applied thereto. The electronic spin or spins may be further configured to undergo a Zeeman shift in energy level that is proportional to the magnetic field to be sensed, in response to an external perturbation such as an RF (radio frequency) field. The detector may be configured to detect output optical radiation from the electronic spin or spins, e.g., by using confocal microscopy, and thereby determine the Zeeman shift and the magnetic field.



BRIEF DESCRIPTION OF DRAWINGS
The drawing figures depict one or more implementations in accord with the present concepts, by way of example only, not by way of limitations. The drawings disclose illustrative embodiments. They do not set forth all embodiments. Other embodiments may be used in addition or instead. In the figures, like reference numerals refer to the same or similar elements.
FIG. 1 is a schematic block diagram of a solid state magnetometer including a single NV center, in accordance with one embodiment of the present disclosure.
FIG. 2 illustrates the energy level structure of a single NV center.
FIG. 3 illustrates the crystal structure of diamond with a (111) NV center.
FIG. 4 illustrates the DC (direct current) and AC (alternating current) sensitivity to magnetic fields for a single NV center as a function of signal frequency.
FIG. 5 is a schematic flow chart of a method of detecting a magnetic field, in one embodiment of the present disclosure.
FIG. 6 is a schematic diagram of a magnetometer for probing and measuring localized sources of magnetic fields, including a nano-crystal of diamond attached to end of an optical waveguide, in accordance with one embodiment of the present disclosure.
FIG. 7 is a schematic block diagram of a solid state magnetometer including a relatively high density of NV centers within a macroscopic piece of diamond, in accordance with one embodiment of the present disclosure.
FIG. 8 is a graph of the simulated sensitivity as a function of NV center density, in an embodiment in which a solid state electronic spin system is used that includes a plurality of NV centers.
FIG. 9 illustrates a pulse sequence for dynamically decoupling unwanted interactions.



DETAILED DESCRIPTION
In the present disclosure, methods and systems are described for using solid-state electronic spin systems, including but not limited to defects in crystals and semiconductors, for precision sensing of weak magnetic fields and magnetic fields in small volumes. It is shown that solid state systems containing electronic spins can be used as high sensitivity magnetometers by manipulating the electronic spins and probing their behaviors with light and with RF (radiofrequency) electromagnetic fields. Nanometer-scale spatial resolution together with high magnetic field sensitivity can be obtained by probing a single electronic spin system. By probing a plurality of electronic spins (typically an ensemble of electronic spins), at a relatively high density, the sensitivity to magnetic fields may be substantially increased, even though spatial sensitivity may decrease somewhat. Augmentation with optical imaging techniques (e.g., using CCD arrays) can allow imaging of magnetic field spatial patterns. Integration with AFM (atomic force microscopy), magnetic field gradients, super-resolution (e.g., STED (stimulated emission depletion) microscopy) or other techniques can allow magnetic field patterns to be imaged with spatial resolution better than the conventional optical diffraction limit, at a scale near or better than ten nanometers
NV centers in the diamond lattice is one example of such solid state electronic spin systems, which are atom-like. As described below, these electronic spin systems can be used to detect magnetic fields with unprecedentedly high sensitivity, as well as to detect the variation in the magnetic field with very high spatial resolution.
FIG. 1 is a schematic block diagram of a high sensitivity solid state magnetometer 100 based on a NV center in diamond, in accordance with one embodiment of the present disclosure. In overview, the solid state magnetometer 100 includes an electronic spin system 110 that includes at least one Nitrogen vacancy 115 in a carbon lattice 117 of diamond; and a detector 120 configured to detect the Zeeman shift in the energy level of the electronic spin system 110 in response to light 135 and an RF excitation field 145 applied thereto. The magnetometer 100 may further include an optical source 130 and an RF field generator 140.
A NV center in diamond is a crystallographic defect in the structure of a diamond, i.e., an empty position in a diamond's lattice. The vacancies may interact with interstitial atoms, such as Nitrogen, and may act as color centers by absorbing visible light. The NV center is found as a defect in the lattice structure of a single crystal diamond. The NV impurity is based in the lattice of carbon atoms, where two adjacent sites are altered, because one carbon atom is replaced with a Nitrogen atom and the other space is left vacant. NV centers are visible as red spots when illuminated by laser.
The NV center is relatively insulated from magnetic interference from other spins. The quantum state of the spin of the NV center may be probed controlled at room temperature. NV centers in diamond, as well as systems involving other types of defects in regular solid state lattices, may provide electronic spins that have almost no interaction with the background lattice, i.e., nearly pure electronic spins that are practically frozen in space, with almost no corrupting interactions with the background lattice. Such electronic spins, which are practically motionless in free space, may be optically detectable with unique optical signatures that allow them to be used for magnetometry.
In the illustrated embodiment, the NV center 115 is an idealized electronic spin, as described above. The NV center 115 is configured to receive optical radiation 135 from the optical source 130, and to align itself with the magnetic field in response to the optical radiation 135.
The optical source 130 may be configured to generate light 135 which, when applied to the electronic spin system, i.e., the NV center 115, causes the electronic spin to align along one or more well-defined directions. For example, the optical source 130 may be a laser. In one embodiment, the laser may generate light 135 that has a wavelength of about 520 nm and has green color.
An external perturbation, for example an ESR (electron spin resonance) field 145 or pulses of RF radiation that are applied to the electronic spin or NV center 115 may cause the electronic spin 115 to precess. The external perturbation may cause a detectable splitting in the electronic spin energy levels, i.e. may cause a detectable Zeeman shift that is proportional to the magnetic field being measured. In the embodiment illustrated in FIG. 1, an RF source 140 may generate an RF field, or pulses of RF radiation.
In the present application, the terms “ESR field”, “RF field,” and “RF radiation” have the same meaning, and are used interchangeably.
The magnetometer 100 operates by detecting the accumulated phase resulting from the relative energy shift induced by a magnetic field b between two Zeeman sublevels,




  δω
  ∝
  
    
      
        g
        ⁢
        
          
        
        ⁢
        
          μ
          B
        
      
      ħ
    
    ⁢
    
      b
      .
    
  




This allows for the precise estimation of an applied DC or AC magnetic field b. The ultimate sensitivity of such sensors is determined by the linewidth of the spin transition and by the signal-to-noise ratio.

A detector 120 detects output optical radiation from the NV center 115 after the light 135 from the laser 130 and the RF or ESR field 145 have passed therethrough. In some embodiments, the detector may be configured to create an image of the optical radiation emitted from electronic spins in different regions of the solid-state system, thereby allowing the magnetic field to be imaged. The detector may also employ super-resolution (for example, STED), atomic force microscopy, magnetic field gradients or other techniques to provide location of the electronic spin or spins (and hence the magnetic field) better than the conventional optical diffraction limit. The detector 120 thus may include, but is not limited to, combinations of CCD (charge coupled device) arrays, confocal microscopy, super-resolution optical techniques (e.g., STED), optical fibers, plasmonic waveguides, atomic force microscopy, and magnetic field gradients. A control and data processing system, which may be integrated with the detector 120, may determine the magnitude of the Zeeman shift and thus calculate the magnetic field, from the detected output radiation from the NV center 115.
FIGS. 2-4 illustrate further details of magnetometry in accordance with the embodiment illustrated in FIG. 1, namely an embodiment in which magnetometers are used that contain at least one Nitrogen-vacancy center(s) in diamond.
FIG. 2 illustrates the level structure of a single NV center. As shown in FIG. 2, the fine structure of the electronic ground state of a NV center is a spin triplet (having |1>, |0>, and |−1> states), with a Δ=2.88 GHz crystal field splitting and a Zeeman shift δω. As indicated in FIG. 2, under the application of excitation light 135 having a wavelength of about  ̃530 nm, the NV center exhibits SDPL (spin-dependent photo-luminescence) near the zero phonon line at about 638 nm, even at room temperature. This allows for optical detection of ESR and optical pumping of the spin state of the NV center, as well as allowing phase shifts resulting from the magnetic field to be sensed: the excitation light helps put the spins into the |1>, |0>, or |−1> states, in varying proportions.
The crystal field splits the ms=±1 sub-levels from the mst=0 sub-level by Δ=2π×2.88 GHz, allowing resonance techniques at vanishing external field. At low field, Δ is the highest energy scale and sets the quantization axis to be parallel to the Nitrogen-to-vacancy direction.
The secular Hamiltonian, including a small external field {right arrow over (B)}(t)≦(Bx, By, Bz), is

H=ħΔSz2+gμBBzSz 

where Bz, is the component of the magnetic field along the NV center's axis and Sz, takes the values ms=0, ±1. Terms proportional to the perpendicular field are suppressed to order Bx,y2/Δ, and may be neglected.
By using the ms=±1 manifold in a low external magnetic field, the NV center acts as a vector magnetometer, sensitive only to fields aligned with its quantization axis. This manifold is twice as sensitive to the signal field than the 0-1 manifold, can operate with high-power microwave ESR fields even at vanishing external magnetic fields, and is insensitive to phase errors in the microwave source. Furthermore, the ±1 manifold is less affected by nuclear spin-induced decoherence at low fields, since inter-nuclei interactions are effectively quenched by the large hyperfine field and dephasing can be refocused by spin echo.
FIG. 3 illustrates the crystal structure of a diamond lattice 117 having a (111) NV center 115. As seen in FIG. 3, there are four crystallographic classes of centers, each oriented along a different (111) axis. For this reason, a transverse (DC) magnetic field B⊥≧0.3 mT must be applied to detune the other three classes. Thus, the ESR field used for quantum control only affects the class whose crystalline axis is aligned with the external magnetic field.
In one embodiment of the present disclosure, coherent control of the NV center's states may be obtained via an ESR magnetic field, oscillating at the angular frequency Δ. In some embodiments, a variety of ESR control techniques may be implemented. It is known that the NV center's spin triplet is an inverted lambda system. An external linearly polarized microwave magnetic field tuned to the crystal field splitting drives transitions between |0 and the superposition |+=(|1+|−1)/√{square root over (2)}, while the state |−=(|1−|−1)/√{square root over (2)} is dark—it is decoupled from the field due to destructive interference.
Application of a small magnetic field aligned with the NV center z axis breaks this interference effect, and allows for complete quantum control of the spin triplet manifold. An echo sequence would have the traditional π/2−π−π/2 structure replaced by π−2π−π: the first pulse creates |+, the second induces a relative π−phase between |+ and |−, and the third converts |+ to |0 while leaving |− population trapped in the ms=±1 manifold.
Linearly polarized ESR pulses may rotate between the two dimensional sub-space of |0and |+=(|1+|−1)/√{square root over (2)}. Additional phase control may be provided by a background oscillating reference field └Bref sin(2π/τ)┘ along the z-axis. In one embodiment, choosing Bref=(ħ/gμB)π2/8τ may yield an optimal phase offset to achieve a signal linear in the field strength.
FIG. 4 illustrates the sensitivity as a function of the signal frequency for both AC and DC detection, for the case of a single NV center. For diamond with natural abundance of C-13 nuclei, T2*≈1 μs and T2≈0.3 ms, sensitivities of about ηDC≈200 nT Hz−1/2 and ηAC≈10 nT Hz−1/2 may be obtained for a single NV center with this approach.
In one embodiment of the present disclosure, a typical approach to the measurement of magnetic fields may be adopted, in which a Ramsey-type decoupling sequence is used, as illustrated in the right hand side of FIG. 4. As shown in FIG. 4, a π/2-pulse creates a superposition of two Zeeman levels, which acquire a relative phase




  φ
  =
  
    δωτ
    ∝
    
      
        
          g
          ⁢
          
            
          
          ⁢
          
            μ
            B
          
        
        ħ
      
      ⁢
      b
    
  




due to the external field b during the free evolution interval τ, where μB is the Bohr magneton and g is about ≈2 for NV centers. Another π/2-pulse transforms the phase into a population difference, from which the Zeeman shift is inferred by SDPL. For small φ, the signal depends linearly on the field:





  s
  ≈
  
    
      
        g
        ⁢
        
          
        
        ⁢
        
          μ
          B
        
      
      ħ
    
    ⁢
    b
    ⁢
    
      
    
    ⁢
    
      τ
      .
    
  




In the embodiment illustrated in FIG. 4, the following relaxation and measurement times, appropriate for the case of NV centers where the bath is constituted by the 1.1% natural abundant Carbon-13 spins, were used: T2*=1 μs6, T2=250 μs7, and tM=1 μs. In addition, a measurement efficiency parameter C was set to C=0.5.
Higher frequency signals may be measured using composite pulse sequences, which include but are not limited to CP (Carr-Purcell) sequences, with an improvement in sensitivity. Ideally, the sensitivity could improve indefinitely at higher frequencies by increasing the number of pulses, as seen by the dotted line. In practice, pulse errors will eventually reduce the achievable sensitivity, as shown by the solid line. A pulse error probability of 0.002 was used in the illustrated embodiment.
The control sequences illustrated in the top right hand side of FIG. 4 is the DC pulse sequence. Starting with a polarized spin, an initial π/2 pulse about the x axis is shown as rotating the spin to the y axis, where it can precess for a free evolution interval τ. This free evolution leads to a phase rotation about the z axis by an amount φ=δωτ. A final π/2 pulse about the y axis is shown as rotating the x component of the spin to the z axis, where measurement of the spin reveals oscillations in the z component ∂ sin(φ).
In the bottom right hand side of FIG. 4, the AC pulse sequence is illustrated, consisting of π/2|x−π|x−π/2|y. For small accumulated phases, a signal linear in the field can also be obtained with all the pulses along the x direction if a reference field Bref sin(2πt/τ) is added. The measurement of the population in the |0 state at the end of the sequence yields P=1⁄2[1+cos(δφ)], where




  
    δφ
    22
  
  =
  
    
      
        
          g
          ⁢
          
            
          
          ⁢
          
            μ
            B
          
        
        
          ħ
          /
          τ
        
      
      ⁡
      
        [
        
          
            4
            ⁢
            
              
                B
                ref
              
              /
              π
            
          
          +
          
            b
            ⁢
            
              
            
            ⁢
            
              f
              ⁡
              
                (
                
                  ντ
                  ,
                  φ
                
                )
              
            
          
        
        ]
      
    
    .
  




The reference field provides an offset of the cosine function to its point of maximal slope at





  
    B
    ref
  
  =
  
    
      
        π
        2
      
      ⁢
      ħ
    
    
      8
      ⁢
      g
      ⁢
      
        
      
      ⁢
      
        μ
        B
      
      ⁢
      τ
    
  




resulting in





  
    
      
        P
        -
        
          1
          /
          2
        
      
      ≈
      
        
          
            g
            ⁢
            
              
            
            ⁢
            
              μ
              B
            
            ⁢
            τ
          
          ħ
        
        ⁢
        b
        ⁢
        
          
        
        ⁢
        
          f
          ⁡
          
            (
            
              ντ
              ,
              φ
            
            )
          
        
      
    
    ]
  
  .




Increasing τ improves the sensitivity until random perturbations, such as environmental perturbations, lead to decoherence and the accompanying decay of the free-precession signal. A principle difference between isolated atomic systems and solid-state electronic spins is that the latter couple strongly to their environment, resulting in relatively short coherence times and a correspondingly short τ. In the case of electronic spin systems, the coherence is limited by interactions with lattice nuclei, resulting in an ensemble dephasing time T2*. Furthermore, the measurement efficiency, denoted by C, may be less than unity, as further described below. Thus the sensitivity of a magnetometer based a single electronic spin is given by




  
    η
    
      D
      ⁢
      
        
      
      ⁢
      C
    
  
  ≈
  
    
      ħ
      
        g
        ⁢
        
          
        
        ⁢
        
          μ
          B
        
        ⁢
        C
        ⁢
        
          
            T
            2
            *
          
        
      
    
    .
  




The sensitivity of the magnetometer may be improved for AC fields with coherent control techniques. Due to the long correlation times of the lattice nuclei, a spin echo may dramatically extend the relevant coherence times. This approach, illustrated in FIG. 4, is a Hahn-echo sequence, in which the Ramsey free evolution is interrupted by an additional microwave π pulse τ/2. The resulting averaged evolution removes slowly varying fields. At the same time, a signal field b(t) oscillating in-phase with the pulse sequence produces an additive phase shift. The phase accumulated in the two time intervals is




  δφ
  =
  
    
      
        
          g
          ⁢
          
            
          
          ⁢
          
            μ
            B
          
        
        ħ
      
      ⁡
      
        [
        
          
            
              ∫
              0
              
                τ
                /
                2
              
            
            ⁢
            
              
                b
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                d
                t
              
            
          
          -
          
            
              ∫
              
                τ
                /
                2
              
              τ
            
            ⁢
            
              
                b
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                d
                t
              
            
          
        
        ]
      
    
    .
  




For a signal field b(t)=bsin(νt+l0), with





  
    δφ
    =
    
      
        
          g
          ⁢
          
            
          
          ⁢
          
            μ
            B
          
        
        ħ
      
      ⁢
      b
      ⁢
      
        
      
      ⁢
      τ
      ⁢
      
        
      
      ⁢
      
        f
        ⁡
        
          (
          
            
              ν
              ⁢
              
                
              
              ⁢
              t
            
            +
            
              l
              0
            
          
          )
        
      
    
  
  ,




with





  
    f
    ⁡
    
      (
      
        x
        ,
        
          l
          0
        
      
      )
    
  
  =
  
    
      
        
          sin
          2
        
        (
        
          
            x
            /
            4
          
          ⁢
          
            
          
          ⁢
          
            cos
            (
            
              
                x
                /
                2
              
              +
              
                l
                0
              
            
          
        
      
      
        x
        /
        4
      
    
    .
  




For maximal response to CW signals (assuming small b), τ=2λ/ν and l0=0 were found to be optimal. In essence, the spin echo allows the limit associated with T*2 to be extended to an intrinsic dephasing time T*2, at the cost of a reduced bandwidth and insensitivity to frequencies ≦1/T2. For a bath of dipolar-coupled spins, the signal is expected to decay as exp[−((τ/T2)3], and the sensitivity is improved by √{square root over (T2*/T2:)}




  
    η
    AC
  
  ≈
  
    
      πħ
      
        2
        ⁢
        
          
        
        ⁢
        g
        ⁢
        
          
        
        ⁢
        
          μ
          B
        
        ⁢
        C
        ⁢
        
          
            T
            2
          
        
      
    
    .
  




It is possible to further improve the sensitivity for higher-frequency signals, since the interrogation time need not be limited by the period of the oscillation but can be multiples of it. By using a sequence of 2ηcπ-pulses, the interrogation time can be extended to 2πnc/ν≈TCP, where TCP is the effective decay time under the pulse sequence. This decay time is usually longer than T2: for example, TCP ̃nc2/3T2 for a Lorentzian noise spectrum. This results in further enhancement of sensitivity that is, in principle, limited only by errors associated with refocusing pulses.
FIG. 5 is a schematic flow chart of a method 500 of detecting a magnetic field, in one embodiment of the present disclosure. In operation, the method 500 may include an act 510 of applying optical radiation to a solid state electronic spin system that contains at least one electronic spin within a solid state lattice, thereby aligning the electronic spin with the magnetic field. The method 500 may further include an act 520 of applying one or more pulses of RF radiation to the solid state electronic spin system so as to rotate the electronic spins and cause a Zeeman shift in the electronic spin energy level, proportional to the magnetic field.
The method 500 may further include an act 530 of detecting output radiation from the solid state electronic spin system after the optical excitation radiation and the RF radiation have passed through the electronic spins in the solid state electronic spin system, so as to determine the Zeeman shift and the magnetic field from the output radiation.
Several approaches may be envisioned for implementing magnetic sensors using the above-described methods and systems. For example, a diamond nano-sized crystal containing one or several NV centers can be attached to a tip of a tapered fiber, an atomic force microscope, or a plasmonic waveguide. This may be used as a scanning magnetic probe, capable of probing localized sources of magnetic fields. Magnetic field gradients may be applied to increase the spatial resolution of the sensor by localizing the region of electron spin resonance with the applied RF radiation. Confocal microscopy and super-resolution optical techniques (such as STED) may be used to image NV centers in diamond nano-crystals with spatial resolution better than the conventional diffraction limit.
FIG. 6 is a schematic diagram of a magnetometer 600 for probing and measuring localized sources of magnetic fields, in accordance with one embodiment of the present disclosure. In the embodiment illustrated in FIG. 6, the magnetometer 600 may include a nano-crystal of diamond 607 that is attached at the and of an optical waveguide 610 for photon transmission and collection. The crystal 607 may include a plurality of NV centers 615, each of which acts as a substantially free electronic spin, and each of which is configured to receive optical excitation radiation 135 (optical source not shown for clarity) and RF radiation (not shown for clarity). The optical waveguide 610 may be configured to collect and transmit therethrough photons of the optical radiation and the RF radiation.
Upon receiving the photons of optical radiation that have transmitted through the optical waveguide 610, the electronic spins or NV centers 615 may be configured to align themselves with the magnetic field being sensed. Upon receiving the pulses of RF radiation that have passed through the optical waveguide, the electronic spins 615 may be configured to accumulate a phase resulting from the Zeeman shift induced by the magnetic signal field, which is detected by a detector (not shown) which detects output radiation 670 from the electronic spins 615.
In the embodiment illustrated in FIG. 6, the resolution is limited by the size of the crystal 607. By way of example, a  ̃25 nm diameter nano-crystal of diamond with about 10 NV centers, attached to the end of an optical fiber or plasmonic waveguide 610 shown and discussed above in conjunction with FIG. 6 would provide sub-100 nm magnetic field resolution while achieving orders of magnitude higher sensitivity than magnetic force microscopy. Provided the waveguide 610 yield high collection efficiency (≈50%), the sensitivity may approach 5 nT Hz−1/2, an order of magnitude better than Hall bar- or SQUID-based probes with 10 times better spatial resolution.
This may be particularly important as the field from a magnetic dipole decays as the cube of the distance. Thus, smaller sensors may be brought closer to the source, yielding an improved sensitivity. For example, the magnetic field from a single proton is  ̃10 nT at 20 nm separation, which an NV nanocrystal magnetometer would be able to detect.
In an augmentation to the embodiment illustrated in FIG. 6, a magnetic field gradient could be applied across the diamond nano-crystal, in order to improve spatial resolution by localizing the spatial region in which the NV center's ESR transition is on resonance with the applied RF radiation. With magnetic field gradients provided by existing technologies, such as nano-scale magnetic tips, this technique would allow spatial localization to <1 nm.
In an alternative embodiment, a bulk diamond sample with a high density of NV centers can be used to sense fields created by extended objects or remote objects, with ultra-high sensitivity.
One principal advantage of using solid-state electronic spins is the potential to achieve extremely high density of sensing spins. This directly translates into an improvement of the sensitivity to the average field over the magnetometer volume, since the shot-noise decreases with increasing electronic spin density n as 1/√{square root over (n)}. Although atomic magnetometers cannot achieve high densities or small volumes without collision-induced reductions in coherence properties, spins in a lattice have fixed locations which allow for substantially higher densities before the deleterious effects of spin-spin interactions impair the magnetometer's operation.
FIG. 7 is a schematic block diagram of a solid state magnetometer 700 including a relatively high density of NV centers 715 within a macroscopic piece of diamond 707, in accordance with one embodiment of the present disclosure. In the embodiment illustrated in FIG. 7, the resolution is limited by the wavelength of light. This embodiment allows for very high spatial resolution and signal-to-noise.
In the illustrated embodiment, green laser light 135 produces SDPL, which is detected by measuring red light 770 from the NV centers imaged onto a CCD 720. A dichroic mirror 750 is configured to direct the optical radiation 135 from the laser toward the solid state electronic spin system, i.e. the bulk diamond sample 707, and to direct the output radiation 770 from the solid state electronic spin system 707 toward the CCD 720. The microwave source for spin-resonance is not shown.
In natural diamonds, the density of the desired (negatively charged) NV centers is low (<1 ppm), but it may be artificially increased. To create NV centers in a controlled way, Nitrogen ions may be implanted into diamonds that contain low concentrations of native Nitrogen, followed by annealing to recombine the Nitrogens with vacancies. Assuming a Nitrogen concentration of about 250 ppm with about a 10% conversion, a diamond crystal may be created with a NV center density exceeding about  ̃1018 cm−3, with an average distance of less than about 10 nm.
At higher densities, the couplings among the electronic spins may no longer be neglected, as they may start to be comparable to 1/T2 (4 kHz at 25 nm separation, corresponding to densities  ̃1017 cm−3). The coupling of a NV center electronic spin to other NV centers ({right arrow over (S)}k) and paramagnetic impurities ({right arrow over (I)}k, such as Nitrogen) is given by the magnetic dipolar interaction.
To first order in 1/Δ only the terms in the dipolar Hamiltonian commuting with the large crystal-field splitting are effective: Hspin-spin=Hzz+Hepr, with Hzz=ΣjkSz,j{right arrow over (D)}jk·zŜz,k and Hepr=ΣjkSz,j{right arrow over (D)}jk·zŜz,k. The dipole interaction vector is




  
    
      
        D
        →
      
      jk
    
    =
    
      
        μ
        ⁢
        
          
        
        ⁢
        
          og
          2
        
        ⁢
        
          
            μ
            B
            2
          
          ⁡
          
            [
            
              3
              ⁢
              
                (
                
                  
                    
                      τ
                      ^
                    
                    ⁢
                    
                      jk
                      ·
                      
                        
                          r
                          ^
                        
                        jk
                      
                    
                  
                  -
                  
                    z
                    ^
                  
                
                )
              
            
            ]
          
        
      
      
        4
        ⁢
        πħ
        ⁢
        
          
        
        ⁢
        
          r
          jk
          3
        
      
    
  
  ,




with the {circumflex over (z)} axis set by the N-V crystal axis of the sensing spin centers. Spin echo removes, to first order, the effects of the paramagnetic impurities and spectator NV centers (aligned along different crystalline axes), which are out of resonance with the ESR field due to the biasing field B⊥. Additional effects, due to the rapid precession caused by B⊥, lead to a small reduction of the echo envelope, which only causes a slight time- and signal-independent decrease of the spin-echo signal.

However, the interaction Hzz between NV centers with the same crystalline axis is unaffected by the spin echo. This interaction results in a decay of the AC magnetometer signal given by Πk′ cos2(Djkτ)≈exp(−D2τ2/2) with a dephasing rate averaged over NV center positions




  
    ∼
    
      
        〈
        
          D
          2
        
        〉
      
    
  
  =
  
    
      
        
          1
          N
        
        ⁢
        
          
            ∑
            jk
          
          ⁢
          
            
          
          ⁢
          
            D
            jk
            2
          
        
      
    
    .
  




Since the dipolar interaction is proportional to 1/r3, D2 scales as the density squared and the sensitivity per root volume ηACV=ηAC√{square root over (V)} is





  
    
      η
      AC
      V
    
    =
    
      
        ħ
        
          g
          ⁢
          
            
          
          ⁢
          
            μ
            B
          
        
      
      ⁢
      
        
          π
          ⁢
          
            
          
          ⁢
          
            e
            
              
                (
                
                  τ
                  /
                  
                    T
                    2
                  
                
                )
              
              3
            
          
        
        
          C
          ⁢
          
            
              n
              ⁢
              
                
              
              ⁢
              τ
            
          
        
      
      ⁢
      
        e
        
          
            an
            2
          
          ⁢
          
            τ
            2
          
        
      
    
  
  ,




where α is on the order of




  
    
      
        (
        
          
            
              μ
              0
            
            
              4
              ⁢
              π
            
          
          ⁢
          
            
              
                (
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  
                    μ
                    B
                  
                
                )
              
              2
            
            ħ
          
        
        )
      
      2
    
    ≈
    
      1.1
      ×
      
        10
        
          -
          25
        
      
      ⁢
      
        s
        
          -
          2
        
      
      ⁢
      
        
      
      ⁢
      
        cm
        6
      
    
  
  ,




and where the fact that the sensing centers are only one fourth of the total NV centers has been taken into account.

FIG. 8 illustrates the achievable sensitivity as a function of density, i.e. shows a graph of the simulated sensitivity as a function of NV center density, in an embodiment in which a solid state electronic spin system is used that includes a plurality of NV centers.
In FIG. 8, the simulated sensitivity in aT Hz−1/2 cm3/2 is shown as a function of NV center density, with dipolar interactions refocused by either a standard echo technique (plotted using stars) or a modified MREV8 sequence (plotted using circles). The simulation was performed for 10 spins. The decoupling sequence was repeated 3 times (102 pulses) with a τ spacing of ≈2 μs and the π-pulse duration 50 ns. The magnetic field frequency was 10 kHz (close to the optimal frequency for T2=250 μs and a 3-echo measurement) yielding a sensitivity as low as 4 aT Hz−1/2 cm3/2. The dotted lines indicate the analytical results from the moment expansion calculations.
To push the sensitivity limits beyond the cutoff imposed by spin-spin interactions and take advantage of the high NV center densities, ideas may be borrowed from dynamical decoupling, of which spin echo is an archetypal sequence. By applying external, time-dependent controls, this technique causes the system to evolve under an effective time-averaged Hamiltonian that is an appropriate symmetrization of the undesired interactions. For example, by using collective rotations, the interaction Hzz can be successively rotated through the x, y and z axes for the same time duration, so that on average the spins will experience the Hamiltonian:




  
    H
    iso
  
  =
  
    
      ∑
      jk
    
    ⁢
    
      
    
    ⁢
    
      
        
          D
          _
        
        jk
      
      ⁡
      
        (
        
          
            
              S
              
                j
                ,
                x
              
            
            ⁢
            
              S
              
                k
                ,
                x
              
            
          
          +
          
            
              S
              
                y
                ,
                l
              
            
            ⁢
            
              S
              
                k
                ,
                y
              
            
          
          +
          
            
              S
              
                j
                ,
                z
              
            
            ⁢
            
              S
              
                k
                ,
                z
              
            
          
        
        )
      
    
  




Since this isotropic Heisenberg Hamiltonian commutes with the signal perturbation and the external magnetic field, assuming that the magnetic field is spatially homogeneous over the magnetometer volume, the Hamiltonian causes no harm to the spin evolution necessary for magnetometry. The control sequence also affects the spins interaction with the signal field, whose average is rotated and rescaled.
Embedding this control sequence within a spin echo sequence simultaneously transforms the dipolar interaction into the harmless isotropic interaction, while averaging to zero all time in-variant linear terms ∂Sz, as required to avoid the dephasing caused by nuclear spins and other paramagnetic impurities. Simple sequences perform the averaging only to first order in time but more complex concatenated sequences can achieve averaging of higher order terms. The convergence to the averaged interaction holds provided the pulse switching frequency is fast on the scale set by the internal couplings and the bath spectral density, a condition that is expected to be accessible experimentally.
In one embodiment, an MREV8-type sequence may achieve these goals, as shown in FIG. 9. FIG. 9 illustrates a Ramsey-type pulse sequence for dynamically decoupling unwanted interactions. Dynamical decoupling of unwanted interactions may be obtained by coherent averaging. By means of an external control the internal Hamiltonian may be made time-dependent. Using cyclic MPS (multiple pulse sequences) and considering only stroboscopic measurements, the evolution may be described by an effective Hamiltonian that, to leading order in time, is given by the time average of the modulated internal Hamiltonian. Higher order terms in the effective Hamiltonian expansion, Heff= H(1)+ H(2)+ . . . , may also be calculated using average Hamiltonian theory and the MPS tailored to produce the desired evolution.
In the exemplary decoupling sequence shown in FIG. 9, the narrow bars are π/2 pulses around different axes in the ms=±1 manifold, while the wide bars are π-pulses. The overall pulse sequence illustrated in FIG. 9 is comprised of 34 pulses, with a cycle time of 48τ. The variable b in FIG. 9 stands for the AC field to be measured. During the total averaging interval T, T/τ measurements can be made, yielding a shot-noise-limited sensitivity η given by the minimum detectable field,




  
    
      b
      min
    
    ≡
    
      η
      /
      
        T
      
    
  
  =
  
    
      ħ
      
        g
        ⁢
        
          
        
        ⁢
        
          μ
          B
        
      
    
    ⁢
    
      
        1
        
          T
        
      
      /
      
        
          T
        
        .
      
    
  




Increasing the number of sensing spins results in further improvement of sensitivity.

In the inset in FIG. 9, the MREV8 sequence is shown, which comprises 8 pulses with cycle time 12τ. The z-component of the spin vector is cyclically rotated through the y, x and z direction so that the average spin vector is along the [1, 1, 1] axis, as seen in the sphere illustrated in FIG. 2. The Hzz Hamiltonian becomes Hzz(1)=Hiso/3 and the linear Hamiltonian Sz→(Sz+Sx)/3. {tilde over (H)}, shown in the inset, is the direction of the internal Hamiltonian in the interaction frame.
By preparing the spins in a direction ≈π/4 in the z-x plane the accumulated phase is √{square root over (2/3)}≈0.47 of that acquired during an echo sequence of equivalent duration and the most important error contribution from the residual dipolar Hamiltonian is a third order correction term. Thus the sensitivity improves to




  
    η
    AC
    V
  
  =
  
    
      ħ
      
        g
        ⁢
        
          
        
        ⁢
        
          μ
          B
        
      
    
    ⁢
    
      
        3
        ⁢
        π
        ⁢
        
          
        
        ⁢
        
          e
          
            
              (
              
                τ
                /
                
                  T
                  2
                
              
              )
            
            3
          
        
      
      
        C
        ⁢
        
          
            2
            ⁢
            n
            ⁢
            
              
            
            ⁢
            τ
          
        
      
    
    ⁢
    
      
        e
        
          
            α
            ~
          
          ⁢
          
            n
            6
          
          ⁢
          
            τ
            2
          
        
      
      .
    
  




The coefficient {tilde over (α)} can be estimated to be of order




  
    
      
        (
        
          
            
              μ
              0
            
            
              4
              ⁢
              π
            
          
          ⁢
          
            
              
                (
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  
                    μ
                    B
                  
                
                )
              
              2
            
            ħ
          
        
        )
      
      6
    
    ⁢
    
      δτ
      4
    
  
  ,




where δτ is the interval between pulses and τ the total measurement time (∂48 δτ). In the embodiment shown in FIG. 9, the dynamics of a small number of spins (N=10) was simulated, under the influence of the modified MREV8 sequence as a function of spin density, taking into account the finite width of the pulses.

The improvement in sensitivity with respect to a simple echo measurement is clearly seen in FIG. 9: the sensitivity can be as low as 4 aT Hz−1/2 cm−3/2 at 10 kHz. By varying the number of cycles over which the signal is recorded, using concatenated decoupling sequences, it is possible to obtain similar sensitivities over a wide range of frequencies. The maximum frequency is constrained by the ability to perform fast ESR pulses while the minimum frequency is set by the T2 time.
In the illustrated embodiment, the MREV8 sequence is chosen for its robustness against pulse errors. The linear term Sz (and therefore the signal) is scaled by √{square root over (2/3)} and rotated to the z+x axis. A π-pulse in between two MREV8 sequences may cause any quasi-static term proportional to z to go to zero, while retaining the AC signal. Time symmetrization of the sequence may bring to zero the second order terms.
The dephasing function may then given by the leading order error Hzz(3) (neglecting the second order correction given by the cross-terms of the external field and the dipolar Hamiltonian). Its effects may be captured by a moment expansion I° to second order of the effective Hamiltonian, D6=Tr[[ Hzz(3),J⊥]/Tr[J⊥2], where J⊥=ΣkS⊥,k is the collective spin in a direction perpendicular to z. D3 may have the character of a sixth moment and its value is a function of the sixth power of the local field generated by the dipolar interaction. The sensitivity decay rate is thus proportional to the sixth order of the density and the square of the total time, with a coefficient




  
    a
    ~
  
  =
  
    
      
        D
        6
      
      /
      
        n
        6
      
    
    ∼
    
      
        
          (
          
            
              
                μ
                0
              
              
                4
                ⁢
                π
              
            
            ⁢
            
              
                
                  (
                  
                    g
                    ⁢
                    
                      
                    
                    ⁢
                    
                      μ
                      B
                    
                  
                  )
                
                2
              
              ħ
            
          
          )
        
        6
      
      ⁢
      
        
          τ
          4
        
        .
      
    
  




The MREV8 sequence entails a large number of control pulses. For many typical errors, for example phase-lag and overshoot/undershoot, the refocusing may only be affected at higher order. However, depolarizing pulse errors occurring with probability p lead to a reduction of contrast: C C′=C(1−p)k for k pulses. Using the MREV8 sequence with echo gives k=34 and a requirement p≦0.0002 for contrasts near unity.
More complex detection techniques, such as multiple-pulse echo sequences embedded in the dipolar coupling refocusing sequence, may allow good sensitivity to be maintained at higher frequency signal fields, ranging from tens of kHz up to MHz. Furthermore, integration of a thin crystal of diamond with a CCD detector may allow for individual “pixel” elements, corresponding to individual pixels of the CCD. For a 10 μm thick crystal with (2 μm)2 “pixels”, each pixel could achieve 1 pT Hz−1/2 sensitivity with a few micron spatial resolution of the signal field.
When the signal field is macroscopic, such as in low field MRI, larger diamond crystals may provide extremely sensitive detection. For example, in the detection of proton NMR (426 kHz in a 10 mT field), a 0.01 cm3 diamond with 1018 centers per cm−3 would provide a sensitivity η≈60 aT Hz−1/2 for measurement in ambient conditions. This compares with η≈800 aT Hz−1/2 for SQUIDS and η≈1−0.5 fT Hz−1/2 for atomic magnetometers, while having a substantially smaller (by a factor of 10) detector volume.
Increasing even further the diamond volume could potentially lead to the detection of atto Tesla fields, opening the prospect of improved tests of fundamental symmetries and physical laws. The sensitivity could be also improved by working with synthesized, isotopically purified diamond containing lower fraction of C-13. Other paramagnetic centers in various materials may also be considered. Many potential applications of such nano-sized centers may be found in bio-science.
These considerations indicate that coherent control of electronic spins in diamond can be used to create a magnetic field sensor of unprecedented sensitivity and spatial resolution. Further refinements could include using non-classical spin states, such as squeezed states induced by the spin-spin coupling. Potential applications in bio-science, materials science, fundamental physics and single electron and nuclear spin detection can be foreseen.
High-field magnetometry is another regime of operation for a NV center-based magnetometer. In particular, for fields such that the Larmor precession of Carbon-13 is much faster than the hyperfine interaction with the electronic spins of the NV center, ESEEM (Electron Spin Echo Envelope Modulation) is also highly suppressed. Thus, for fields above  ̃100 mT, the NV system may be still be used for vector magnetometer with the 0-1 manifold.
In some embodiments, the effects of different NV center orientations may be considered. In order to use the NV centers as a vector magnetometer, the signal should originate only from one of the four different crystallographic axes (for example the (111) axis). Under application of a DC transverse magnetic field B⊥{circumflex over (x)}, the spectator centers (with crystalline axis {circumflex over (n)} have the |±1) levels split by gμBB⊥{circumflex over (x)}·{circumflex over (n)}. This detunes the other three classes of centers from the microwave field used for preparing and manipulating the ms=±1 subspace, with ideal choice of {circumflex over (x)} aligned with the




  (
  
    
      1
      2
    
    ,
    
      1
      2
    
    ,
    
      1
      _
    
  
  )




axis. For pulse errors to be smaller than the assumed measurement errors for the desired (111) axis, it is required that Ω≧3π/T2*. This translates to a requirement that gμBB⊥>3ħΩ√{square root over (3/2)} for the other three axes. For T2*=1 μs, it is required that B⊥≧0.3 mT.

In some embodiments of a NV center based magnetometry in accordance with the present disclosure, the four crystallographic classes of NV centers may be exploited to provide a full 3D vector and scalar magnetometer. This may be achieved by changing the direction of the biasing transverse field B⊥ in between measurements. By alternating between three crystalline axes, a full vector measurement of the signal is, in principle, possible. This may reduce the overall integration time for a given axis by a factor of three, as well as requiring additional calibration of the ESR field strengths and measurement efficiencies.
In the high density regime, the spectator centers may also cause unwanted couplings to the active NV centers. In general, the interaction between NV centers and any paramagnetic (epr) impurity reduces the signal by a multiplicative factor




  
    
      1
      -
      
        ɛ
        jk
      
    
    =
    
      
        1
        -
        
          2
          ⁢
          
            
              
                D
                
                  jk
                  ,
                  z
                
                2
              
              +
              
                D
                
                  jk
                  ,
                  y
                
                2
              
            
            
              
                
                
                  
                    g
                    ⁢
                    
                      
                    
                    ⁢
                    
                      μ
                      B
                    
                    ⁢
                    B
                  
                  ⊥
                
                
              
              2
            
          
        
      
      ≈
      
        exp
        ⁡
        
          (
          
            -
            
              ɛ
              jk
            
          
          )
        
      
    
  
  ,




provided a biasing field gμBB⊥>>|{right arrow over (Djk)}|| is applied. For many randomly located impurities, the factor becomes exp(−εj) with





  
    
      ɛ
      jk
    
    =
    
      
        
          ∑
          k
        
        ⁢
        
          
        
        ⁢
        
          ɛ
          jk
        
      
      ≈
      
        
          
            
              (
              
                
                  μ
                  0
                
                ⁢
                g
                ⁢
                
                  
                
                ⁢
                
                  μ
                  B
                
              
              )
            
            2
          
          ⁢
          
            n
            epr
          
        
        
          12
          ⁢
          π
          ⁢
          
            
          
          ⁢
          
            r
            0
            3
          
          ⁢
          
            B
            ⊥
            2
          
        
      
    
  
  ,




where r0 is the closest epr impurity spacing that does not detune the NV center off resonance (r0≈3 nm) and nepr is the density of epr centers. As the bias field is increased in strength, the paramagnetic impurities get averaged by their fast rotation, leading to an improved signal for the NV center system. For B⊥=0.3 mT, ε≈nepr×5×10−33 cm3, which remains negligible even for nN>1019 cm−3.

In sum, the present disclosure provides methods and systems for enhancing the spatial resolution and magnetic field sensitivity of solid-state magnetometers based on electronic spin systems, such as defects in crystals and semiconductors. A particular example has been described of such a solid-state magnetometer system, namely NV-diamond, i.e. the electronic spins associated with Nitrogen vacancy centers in diamond.
Techniques described in the present disclosure include, but are not limited to: (i) coherent control schemes (e.g., using RF pulses) to decouple electronic spins from the local environment and from each other, leading to a substantial improvement in sensitivity to external time-varying magnetic fields and magnetometer operation at high spin densities; (ii) attaching a nano-sized crystal containing one or several spins (e.g., NV centers) to the tip of an atomic force microscope, optical fiber or plasmonic waveguide, to be used as a scanning magnetic probe capable of detecting localized sources of magnetic fields with nanoscale resolution; (iii) integration with super-resolution optical techniques such as STED, for combined far-field imaging and magnetometry with resolution better than the conventional diffraction limit; (iv) use of magnetic field gradients to improve spatial resolution to <1 nm by localizing the region of electronic spin resonance with the applied RF radiation; (v) operation as a vector magnetometer by use of spins in specific crystallographic classes and/or quantization axis; and (vi) use of the measured magnetometer signal variance as the magnetic field probe, e.g., when the magnetic field frequency and/or phase vary in time.
The extremely high magnetic field sensitivity in a small volume offered by solid state spin-cubits such as NV centers in diamond can find a wide range of applications, from fundamental physics tests or quantum computing applications to detection of NMR signals, surface physics and material science, structural biology, intracellular bioimaging, and medical imaging and biomagnetism. This robust technology could be invaluable both in nanoscale magnetic field imaging and in macroscopic field detection scenarios, such as low-field MRI and mapping dynamics of neural circuits.
The methods and systems described in the present disclosure would enable superior spatial resolution and magnetic field sensitivity to existing state-of-the-art techniques, including SQUIDs, atomic vapor-based magnetometers, and magnetic resonance force microscopy. Intended or possible future uses of such methods and systems include, but are not limited to the detection and mapping of AC and DC magnetic fields with high-spatial-resolution, as small as about 0.1 nm, and high-sensitivity, as low as about 1 atto tesla. Such detection and mapping may be performed in a robust, solid-state package that operates at ambient conditions, such as room temperature. The methods and systems described in the present disclosure may find wide applications in improved sensing, diagnostics and guidance of nano-scale manufacturing, for example in biology and medicine. Applications in biology and medicine may include, but are not limited to, structural biology, intracellular bioimaging, mapping neural activity and low-field NMR/MRI, semiconductors and electronics, materials analysis, and detection of rare isotopes, for example for homeland security.
One or more of the systems and methods described above may be implemented using a processing system, including but not limited to a control and data processing system. The methods in the present disclosure are not described with reference to any particular programming language. It will be appreciated that a variety of programming languages may be used to implement the teachings of the present disclosure. The system may be selectively configured and/or activated by a computer program stored in the computer.
Such a computer program may be stored in any computer readable storage medium, including but not limited to, any type of disk including floppy disks, optical disks, CD-ROMs, and magnetic-optical disks, ROMs (read-only memories), RAMs (random access memories), EPROMs (erasable programmable read-only memories), EEPROMs (electrically erasable programmable read-only memories), magnetic or optical cards, or any type of media suitable for storing electronic instructions.
The methods and systems presented herein are not inherently related to any particular computer, processor or other apparatus. Various general purpose systems may be used with different computer programs in accordance with the teachings herein. Any of the methods and systems described in the present disclosure may be implemented in hard-wired circuitry, by programming a general purpose processor, a graphics processor, or by any combination of hardware and software.
It is contemplated that the subject matter described herein may be embodied in many forms. Accordingly, the embodiments described in detail below are illustrative embodiments, and are not to be considered limitations. Other embodiments may be used in addition or instead.
The components, steps, features, objects, benefits and advantages that have been discussed are merely illustrative. None of them, nor the discussions relating to them, are intended to limit the scope of protection in any way. Numerous other embodiments are also contemplated, including embodiments that have fewer, additional, and/or different components, steps, features, objects, benefits and advantages. The components and steps may also be arranged and ordered differently.
The phrase “means for” when used in a claim embraces the corresponding structures and materials that have been described and their equivalents. Similarly, the phrase “step for” when used in a claim embraces the corresponding acts that have been described and their equivalents. The absence of these phrases means that the claim is not limited to any of the corresponding structures, materials, or acts or to their equivalents.
Nothing that has been stated or illustrated is intended to cause a dedication of any component, step, feature, object, benefit, advantage, or equivalent to the public, regardless of whether it is recited in the claims.
In short, the scope of protection is limited solely by the claims that now follow. That scope is intended to be as broad as is reasonably consistent with the language that is used in the claims and to encompass all structural and functional equivalents.

## Claims

1. A magnetometer for sensing an unknown external magnetic field, comprising:
a solid state electronic spin system containing one or more a plurality of electronic spins that are disposed within a solid state lattice, and that are substantially free of interaction with the solid state lattice;
wherein the plurality of electronic spins are configured to align with the unknown external magnetic field in response to optical excitation radiation applied thereto by the optical excitation system; and
wherein the plurality of electronic spins are further responsive to the external control system so as to induce a precession of the electronic spins about the unknown external magnetic field to be sensed, the frequency of the precession being linearly related to the unknown external magnetic field by the Zeeman shift of the electronic spin energy levels; and
a detector configured to detect output optical radiation correlated with the plurality of electronic spins, after the plurality of electronic spins have been subject to the optical excitation radiation and the external control, and wherein the detector is further configured to determine the Zeeman shift and use the determined Zeeman shift to calculate the unknown external magnetic field,
wherein the plurality of electronic spins include different sets of spin centers having different crystallographic orientations within the solid state lattice, and
wherein the magnetometer is configured to select one or more of the sets of spin centers in one or more crystallographic orientations and use the selected spin defects to determine the unknown external magnetic field.

2. The magnetometer of claim 1, further comprising an optical source configured to generate the optical excitation radiation that causes the electronic spins to align with the unknown external magnetic field.

3. The magnetometer of claim 2, wherein the optical source comprises a laser.

4. The magnetometer of claim 1, wherein the solid state electronic spin system comprises a diamond, and wherein the plurality of electronic spins are NV (Nitrogen vacancy) centers.

5. The magnetometer of claim 3, wherein the optical excitation radiation is laser light having a wavelength of about 530 nanometers.

6. The magnetometer of claim 1, wherein the external control comprises RF radiation, and wherein each electronic spin is configured to precess around an axis of the magnetic field in response to the RF radiation.

7. The magnetometer of claim 6, further comprising an RF source configured to generate the RF radiation that induces the precession of the electronic spins about the unknown external magnetic field when applied to the aligned electronic spins.

8. The magnetometer of claim 7, wherein the RF source comprises an RF pulse generator that generates an excitation RF field that can be applied in any desired direction relative to the unknown external magnetic field.

9. The magnetometer of claim 4, wherein the sensitivity of the magnetometer for a single NV center is about 200 nT Hz−1/2 when the unknown external magnetic field is a DC magnetic field, and about 10 nT Hz−1/2 when the unknown external magnetic field is an AC magnetic field.

10. The magnetometer of claim 1, wherein the detector comprises a CCD (charge-coupled device) and one or more associated optical components.

11. The magnetometer of claim 1, wherein the detector comprises one of:
a confocal microscope; and
a super-resolution imaging system configured to implement one or more techniques for imaging below the conventional diffraction limit.

12. The magnetometer of claim 11, wherein the one or more techniques comprise STED (stimulated emission depletion).

13. The magnetometer of claim 6, wherein the detector is configured to employ magnetic field gradients to improve spatial resolution by localizing the spatial region in which the electronic spins are on resonance with the applied RF radiation.

14. The magnetometer of claim 1, wherein the detector comprises a control and data processing system configured to determine the Zeeman shift and the unknown external magnetic field from the detected output optical radiation, and wherein the Zeeman frequency shift δω is proportional to the unknown external magnetic field according to a formula given by:




  δω
  ∝
  
    
      
        g
        ⁢
        
          
        
        ⁢
        
          μ
          B
        
      
      ħ
    
    ⁢
    b
  



where δω represents the angular frequency shift between two Zeeman sublevels of the electronic spins;
b represents the unknown external magnetic field;
g represents the spin g-factor;
μB represents the Bohr magneton; and
ħ represents Planck's constant divided by 2π.

15. The magnetometer of claim 2, further comprising one or more optical components configured to direct the optical excitation radiation from the optical source toward the electronic spin system, and further configured to direct the output optical radiation from the electronic spin system toward the detector.

16. The magnetometer of claim 1, wherein the density of the electronic spins is about 1018 cm−3, and wherein an average distance between the electronic spins is about 10 nm.

17. The magnetometer of claim 15, wherein the one or more optical components comprises a dichroic mirror.

18. The magnetometer of claim 1, wherein the magnetometer is configured to select the one or more sets of spin centers by using a magnetic field to detune the other sets of spin centers.

19. The magnetometer of claim 1, wherein the magnetometer is configured to select each set of spin centers in turn to take magnetic field measurements thereby providing a full 3D vector and scalar magnetometer using more than one set of spin centers within the solid state lattice.

20. A magnetometer for sensing an unknown external magnetic field, comprising:
a solid state electronic spin system containing one or more electronic spins that are disposed within a solid state lattice, and that are substantially free of interaction with the solid state lattice;
wherein the one or more electronic spins are configured to align with the unknown external magnetic field in response to optical excitation radiation applied thereto by the optical excitation system; and
wherein the one or more electronic spins are further responsive to the external control system so as to induce a precession of the one or more electronic spins about the unknown external magnetic field to be sensed, the frequency of the precession being linearly related to the unknown external magnetic field by the Zeeman shift of the electronic spin energy levels; and
a detector configured to detect output optical radiation correlated with the plurality of electronic spins, after the one or more electronic spins have been subject to the optical excitation radiation and the external control, and wherein the detector is further configured to determine the Zeeman shift and use the determined Zeeman shift to calculate the unknown external magnetic field,
wherein the solid state electronic spin system is configured such that the one or more electronic spins are disposed close to a surface of the solid state lattice to enable the one or more electronic spins to be positioned sufficiently close to a source of the unknown external magnetic field to provide nanometer or micron scale spatial resolution.

21. The magnetometer of claim 20, wherein the solid state electronic spin system is in the form of a nanocrystal.

22. The magnetometer of claim 20, wherein the solid state electronic spin system is in the form of a thin crystal plate.

23. The magnetometer of claim 20, further comprising an AFM (atomic force microscope), and wherein the solid state spin system is mounted at the tip of the AFM.

24. The magnetometer of claim 20, wherein the magnetometer is configured to probe a single electronic spin close to the surface of the solid state lattice.

25. The magnetometer of claim 20, wherein the magnetometer is configured to probe a plurality of electronic spins close to the surface of the solid state lattice, the plurality of electronic spins including different sets of spin centers having different crystallographic orientations within the solid state lattice, and wherein the magnetometer is configured to select one or more of the sets of spin centers in one or more crystallographic orientations and use the selected spin defects to determine the unknown external magnetic field.

26. The magnetometer of claim 20, wherein the solid state electronic spin system comprises a diamond, and wherein the one or more electronic spins are NV (Nitrogen vacancy) centers.