The present disclosure generally relates to the field of magnetometers, such as methods and systems for resolving the natural ambiguity of diamond nitrogen vacancy magnetic sensors.
Some embodiments relate to a system. The system may comprise a nitrogen vacancy (NV) diamond material comprising a plurality of NV centers; a magnetic field source; a radio frequency (RF) excitation source configured to provide RF excitation to the NV diamond material; an optical excitation source configured to provide optical excitation to the NV diamond material; an optical detector configured to receive an optical signal emitted by the NV diamond material; and a controller. The controller may be configured to control the RF excitation source to provide pulsed RF excitation to the NV diamond material, and determine a sign of the magnetic field vector at the NV diamond material based on a received light detection signal from the optical detector. The controller may be configured to control the optical excitation source to provide continuous wave optical excitation to the NV diamond. The controller may be further configured to identify Lorentzian peaks in a received light detection signal from the optical detector as a function of RF excitation frequency. The controller may be configured to determine a sign of the magnetic field vector based on an equilibration time for a pair of the identified Lorentzian peaks.
Other embodiments relate to a system. The system may comprise a magneto-optical defect center material; a magnetic field source; a radio frequency (RF) excitation source configured to provide RF excitation to the magneto-optical defect center material; an optical excitation source configured to provide optical excitation to the magneto-optical defect center material; an optical detector configured to receive an optical signal emitted by the magneto-optical defect center material; and a controller. The controller may be configured to control the RF excitation source to provide pulsed RF excitation to the magneto-optical defect center material, control the optical excitation source to provide optical excitation to the magneto-optical defect center material, and determine a sign of the magnetic field vector at the magneto-optical defect center material based on a received light detection signal from the optical detector. The controller may be configured to control the optical excitation source to provide continuous wave optical excitation to the magneto-optical defect center material. The controller may be further configured to identify Lorentzian peaks in a received light detection signal from the optical detector as a function of RF excitation frequency. The controller may be configured to determine a sign of the magnetic field vector based on an equilibration time for a pair of the identified Lorentzian peaks.
Other embodiments relate to a system. The system may comprise a nitrogen vacancy (NV) diamond material comprising a plurality of NV centers; a magnetic field source; a radio frequency (RF) excitation source configured to provide RF excitation to the NV diamond material; an optical excitation source configured to provide optical excitation to the NV diamond material; an optical detector configured to receive an optical signal emitted by the NV diamond material; and a controller. The controller may be configured to determine a first equilibration time for a first peak of a Lorentzian pair based on a received light detection signal from the optical detector, determine a second equilibration time for a second peak of the Lorentzian pair based on a received light detection signal from the optical detector, and determine a sign of the magnetic field vector at the NV diamond material based on the first equilibration time and the second equilibration time. The controller may be configured to assign a positive spin state to the peak of the Lorentzian pair with the longer equilibration time. The first equilibration time and the second equilibration time may be determined by measuring the time to reach 60% of a normalized equilibrium intensity after the beginning of an RF pulse, wherein the normalized equilibrium intensity is determined based on the intensity in the absence of the RF pulse and the equilibrium intensity in the presence of the RF pulse.
It is possible to resolve a magnetic field vector from a diamond nitrogen vacancy magnetic field sensor. The method of determining the sign of the magnetic field vector resolved by the DNV magnetic field sensor described herein may resolve a natural ambiguity of the magnetic field sensor with regard to the sign of the vector. The ability to resolve the sign of the resolved magnetic field vector expands the applications in which the DNV magnetic field sensor may be employed.
NV Center, its Electronic Structure, and Optical and RF Interaction
The nitrogen vacancy (NV) center in diamond comprises a substitutional nitrogen atom in a lattice site adjacent a carbon vacancy as shown in
The NV center may exist in a neutral charge state or a negative charge state. Conventionally, the neutral charge state uses the nomenclature NV0, while the negative charge state uses the nomenclature NV, which is adopted in this description.
The NV center has a number of electrons including three unpaired electrons, each one from the vacancy to a respective of the three carbon atoms adjacent to the vacancy, and a pair of electrons between the nitrogen and the vacancy. The NV center, which is in the negatively charged state, also includes an extra electron.
The NV center has rotational symmetry, and as shown in
Introducing an external magnetic field with a component along the NV axis lifts the degeneracy of the ms=±1 energy levels, splitting the energy levels ms=±1 by an amount 2 gμBBz, where g is the g-factor, μB is the Bohr magneton, and Bz is the component of the external magnetic field along the NV axis. This relationship is correct to a first order and inclusion of higher order corrections is a straightforward matter and should not affect the computational and logic steps in the systems and methods described below.
The NV center electronic structure further includes an excited triplet state 3E with corresponding ms=0 and ms=±1 spin states. The optical transitions between the ground state 3A2 and the excited triplet 3E are predominantly spin conserving, meaning that the optical transitions are between initial and final states which have the same spin. For a direct transition between the excited triplet 3E and the ground state 3A2, a photon of red light is emitted with a photon energy corresponding to the energy difference between the energy levels of the transitions.
There is, however, an alternate non-radiative decay route from the triplet 3E to the ground state 3A2 via intermediate electron states, which are thought to be intermediate singlet states A, E with intermediate energy levels. Significantly, the transition rate from the ms=±1 spin states of the excited triplet 3E to the intermediate energy levels is significantly greater than that from the ms=0 spin state of the excited triplet 3E to the intermediate energy levels. The transition from the singlet states A, E to the ground state triplet 3A2 predominantly decays to the ms=0 spin state over the ms=±1 spin states. These features of the decay from the excited triplet 3E state via the intermediate singlet states A, E to the ground state triplet 3A2 allows that if optical excitation is provided to the system, the optical excitation will eventually pump the NV center into the ms=0 spin state of the ground state 3A2. In this way, the population of the ms=0 spin state of the ground state 3A2 may be “reset” to a maximum polarization determined by the decay rates from the triplet 3E to the intermediate singlet states.
Another feature of the decay is that the fluorescence intensity due to optically stimulating the excited triplet 3E state is less for the ms=±1 states than for the ms=0 spin state. This is so because the decay via the intermediate states does not result in a photon emitted in the fluorescence band, and because of the greater probability that the ms=±1 states of the excited triplet 3E state will decay via the non-radiative decay path. The lower fluorescence intensity for the ms=±1 states than for the ms=0 spin state allows the fluorescence intensity to be used to determine the spin state. As the population of the ms=±1 states increases relative to the ms=0 spin, the overall fluorescence intensity will be reduced.
NV Center, or Magneto-Optical Defect Center, Magnetic Sensor System
The RF excitation source 330 may be a microwave coil, for example. The RF excitation source 330 when emitting RF radiation with a photon energy resonant with the transition energy between ground ms=0 spin state and the ms=+1 spin state excites a transition between those spin states. For such a resonance, the spin state cycles between ground ms=0 spin state and the ms=+1 spin state, reducing the population in the ms=0 spin state and reducing the overall fluorescence at resonance. Similarly resonance occurs between the ms=0 spin state and the ms=−1 spin state of the ground state when the photon energy of the RF radiation emitted by the RF excitation source is the difference in energies of the ms=0 spin state and the ms=−1 spin state. At resonance between the ms=0 spin state and the ms=−1 spin state, or between the ms=0 spin state and the ms=+1 spin state, there is a decrease in the fluorescence intensity.
The optical excitation source 310 may be a laser or a light emitting diode, for example, which emits light in the green, for example. The optical excitation source 310 induces fluorescence in the red, which corresponds to an electronic transition from the excited state to the ground state. Light from the NV diamond material 320 is directed through the optical filter 350 to filter out light in the excitation band (in the green for example), and to pass light in the red fluorescence band, which in turn is detected by the detector 340. The optical excitation light source 310, in addition to exciting fluorescence in the diamond material 320, also serves to reset the population of the ms=0 spin state of the ground state 3A2 to a maximum polarization, or other desired polarization.
For continuous wave excitation, the optical excitation source 310 continuously pumps the NV centers, and the RF excitation source 330 sweeps across a frequency range which includes the zero splitting (when the ms=±1 spin states have the same energy) photon energy of 2.87 GHz. The fluorescence for an RF sweep corresponding to a diamond material 320 with NV centers aligned along a single direction is shown in
In general, the diamond material 320 will have NV centers aligned along directions of four different orientation classes.
While
The RF excitation source 630 may be a microwave coil, for example. The RF excitation source 630 is controlled to emit RF radiation with a photon energy resonant with the transition energy between the ground ms=0 spin state and the ms=±1 spin states as discussed above with respect to
The optical excitation source 610 may be a laser or a light emitting diode, for example, which emits light in the green, for example. The optical excitation source 610 induces fluorescence in the red, which corresponds to an electronic transition from the excited state to the ground state. Light from the NV diamond material 620 is directed through the optical filter 650 to filter out light in the excitation band (in the green for example), and to pass light in the red fluorescence band, which in turn is detected by the optical detector 640. The EMI filter 660 is arranged between the optical filter 650 and the optical detector 640 and suppresses conducted interference. The optical excitation light source 610, in addition to exciting fluorescence in the NV diamond material 620, also serves to reset the population of the ms=0 spin state of the ground state 3A2 to a maximum polarization, or other desired polarization.
The controller 680 is arranged to receive a light detection signal from the optical detector 640 and to control the optical excitation source 610 and the RF excitation source 630. The controller may include a processor 682 and a memory 684, in order to control the operation of the optical excitation source 610 and the RF excitation source 630. The memory 684, which may include a nontransitory computer readable medium, may store instructions to allow the operation of the optical excitation source 610 and the RF excitation source 630 to be controlled.
According to some embodiments of operation, the controller 680 controls the operation such that the optical excitation source 610 continuously pumps the NV centers of the NV diamond material 620. The RF excitation source 630 is controlled to continuously sweep across a frequency range which includes the zero splitting (when the ms=±1 spin states have the same energy) photon energy of 2.87 GHz. When the photon energy of the RF radiation emitted by the RF excitation source 630 is the difference in energies of the ms=0 spin state and the ms=−1 or ms=+1 spin state, the overall fluorescence intensity is reduced at resonance, as discussed above with respect to
As noted above, the diamond material 620 will have NV centers aligned along directions of four different orientation classes, and the component Bz along each of the different orientations may be determined based on the difference in energy between the ms=−1 and the ms=+1 spin states for the respective orientation classes. In certain cases, however, it may be difficult to determine which energy splitting corresponds to which orientation class, due to overlap of the energies, etc. The bias magnet 670 provides a magnetic field, which is preferably uniform on the NV diamond material 620, to separate the energies for the different orientation classes, so that they may be more easily identified.
Natural Ambiguity of NV Center Magnetic Sensor System
The NV center magnetic sensor that operates as described above is capable of resolving a magnetic field to an unsigned vector. As shown in
In some circumstances, real world conditions allow the intelligent assignment of a sign to the unsigned magnetic field vector determined from the fluorescence spectra described above. If a known bias field is used that is much larger than the signal of interest, the sign of the magnetic field vector may be determine by whether the total magnetic field, cumulative of the bias field and the signal of interest, increases or decreases. If the magnetic sensor is employed to detect submarines from a surface ship, assigning the calculated magnetic field vector a sign that would place a detected submarine above the surface ship would be nonsensical. Alternatively, where the sign of the vector is not important a sign can be arbitrarily assigned to the unsigned vector.
It is possible to unambiguously determine a magnetic field vector with a DNV magnetic field sensor. The method of determining the signed magnetic field vector may be performed with a DNV magnetic field sensor of the type shown in
As shown in
As described above, the application of RF excitation to the DNV material produces a decrease in fluorescence intensity at the resonant RF frequencies for the ms=−1 and the ms=+1 spin states. For this reason, at RF frequencies that excite electrons to the ms=−1 and the ms=+1 spin states, an equilibrium fluorescence intensity will be lower than the equilibrium fluorescence intensity in the absence of the applied RF excitation. The time it takes to transition from the equilibrium fluorescence intensity in the absence of RF excitation to the equilibrium fluorescence intensity with the application of RF excitation may be employed to calculate an “equilibration time.”
An “equilibration time” as utilized herein refers to the time between the start of an RF excitation pulse and when a predetermined percentage of the equilibrium fluorescence intensity is achieved. The predetermined amount of the equilibrium fluorescence at which the equilibration time is calculated may be about 20% to about 80% of the equilibrium fluorescence, such as about 30%, 40%, 50%, 60%, or 70% of the equilibrium fluorescence. The equilibration time as shown in
A shown in
The difference in the relaxation time between the electrons of the ms=−1 and the ms=+1 spin states may be measured due to the different RF excitation resonant frequencies for each spin state. As shown in
The Lorentzian pair of the fluorescence spectra which are located furthest from the zero splitting energy may be selected to calculate the equilibration time. These peaks include the least signal interference and noise, allowing a more reliable measurement. The preferred Lorentzian pair is boxed in
A plot of the fluorescence intensity for a single RF pulse as a function of time is shown in
The equilibrium fluorescence intensity under the application of the RF excitation may be set by any appropriate method. According to some embodiments, the RF excitation may be maintained until the intensity becomes constant, and the constant intensity may be considered the equilibrium intensity value utilized to calculate the equilibration time. Alternatively, the equilibrium intensity may be set to the intensity at the end of an RF excitation pulse. According to other embodiments, a decay constant may be calculated based on the measured fluorescence intensity and a theoretical data fit employed to determine the equilibrium intensity value.
The peak in the Lorentzian pair that exhibits the higher measured equilibration time is associated with the higher energy level electron spin state. For this reason, the peak of the Lorentzian pair with the longer equilibration time is assigned the ms=+1 spin state, and the other peak in the Lorentzian pair is assigned the ms=−1 spin state. The signs of the peaks in the other Lorentzian pairs in the fluorescence spectra of the DNV material as a function of RF frequency may then be assigned, and the signed magnetic field vector calculated.
To demonstrate that the equilibration time of each peak in a Lorentzian pair does indeed vary with magnetic field direction, the equilibration time for a single peak in a Lorentzian pair was measured under both a positive and a negative magnetic bias field which were otherwise equivalent. As shown in
The method of determining a sign of a magnetic field vector with a DNV magnetic sensor described herein may be performed with the DNV magnetic field sensor shown in
The controller of the magnetic field sensor may be programmed to determine the location of peaks in a fluorescence spectra of a DNV material as a function of RF frequency. The equilibration time for the peaks of a Lorentzian pair located the furthest from the zero field energy may then be calculated. The controller may be programmed to provide a pulsed RF excitation energy by controlling a RF excitation source and also control an optical excitation source to excite the DNV material with continuous wave optical excitation. The resulting optical signal received at the optical detector may be analyzed by the controller to determine the equilibration time associated with each peak in the manner described above. The controller may be programmed to assign a sign to each peak based on the measured equilibration time. The peak with the greater measured equilibration time may be assigned the ms=+1 spin state.
The method of assigning a sign to a magnetic field vector described above may also be applied to magnetic field sensors based on magneto-optical defect center materials other than DNV.
The DNV magnetic field sensor described herein that produces a signed magnetic field vector may be especially useful in applications in which the direction of a measured magnetic field is important. For example, the DNV magnetic field sensor may be employed in magnetic field based navigation or positioning systems.
The embodiments of the concepts disclosed herein have been described in detail with particular reference to preferred embodiments thereof, but it will be understood by those skilled in the art that variations and modifications can be effected within the spirit and scope of the described concepts.
The present application is related to co-pending U.S. application Ser. No. ______, Attorney Docket No. 111423-1046, filed Jan. 21, 2016, titled “APPARATUS AND METHOD FOR RECOVERY OF THREE DIMENSIONAL MAGNETIC FIELD FROM A MAGNETIC DETECTION SYSTEM”, which is incorporated by reference herein in its entirety.