Simultaneous wide-field imaging of phase and magnitude of AC magnetic signal using diamond quantum magnetometry.

Spectroscopic analysis of AC magnetic signal using diamond quantum magnetometry is a promising technique for inductive imaging. Conventional dynamic decoupling like XY8 provides a high sensitivity of an oscillating magnetic signal with intricate dependence on magnitude and phase, complicating high throughput detection of each parameter. In this study, a simple measurement scheme for independent and simultaneous detection of magnitude and phase is demonstrated by a sequential measurement protocol. Wide-field imaging experiment was performed for an oscillating magnetic field with approximately [Formula: see text]-squared observation area. Single pixel phase precision was [Formula: see text] for [Formula: see text] AC magnetic signal. Our method enables potential applications including inductive inspection and impedance imaging.

Introduction

A negatively charged nitrogen-vacancy (NV) center in diamond offers a promising material platform for quantum sensing[1,2]. Spin-state manipulation with a state-selective microwave (MW) pulse, combined with spin-dependent fluorescence has been utilized for spectroscopic measurement of magnetic field. Fabrication of a dense ensemble of NV centers [3-6] allows these measurements to be applied in wide-field imaging modality[7-10], and local current characterizations using NV center have been achieved[11-16]. Wide-field imaging using NV center paves the way for inductive inspection[17,18], which is still challenging with micron-scale resolution by established methods. Previous studies indicate that dynamical decoupling (DD) protocols like XY8 achieve magnetic field spectroscopy with high sensitivity for magnitude and phase of such a signa[19-22]. Measuring magnitude and phase of oscillating magnetic signal is at the heart of Eddy-current inspection based upon inductive and impedance sensing, enabling conductivity measurement[23,24]. Despite its high sensitivity, DD can only measure the output of fluorescence intensity with an intricate dependence on a magnitude and a phase of a signal. Therefore, a magnitude with a known phase or a phase with a known magnitude can only be measured with DD protocol[25]. Independent and simultaneous measurements of magnitude and phase for magnetic field spectroscopy in wide-field open up new diamond applications such an accurate inductive sensing.

In this study, we propose a stroboscopic measurement termed iQdyne[26], a wide-field modality of Qdyne[27,28] enabled with lock-in detection, as a simple sensing scheme for magnitude and phase of an oscillating magnetic field. The iQdyne provides an orthogonal measurement for magnitude and phase; it involves two input parameters and time-series outputs . Fourier analysis easily extracts two resulting parameters, magnitude and phase , which are separated from each other and are readily interpretable.

We implemented this stroboscopic protocol on a wide-field microscope and demonstrated an imaging experiment of an oscillating magnetic field generated from a current pattern fabricated on a diamond substrate. An oscillating current generates an oscillating magnetic field due to Ampère’s law. The wide-field observation area was approximately , and we estimated that single pixel precision for phase sensing was with AC magnetic signal. This demonstration is a fundamental part of a local current investigation technique like local impedance microscopy and inductive testing method. Wide-field modality provides us an optically high spatial resolution and a wide observation area, leading to high throughputs for such measurements.

Principle

In this section, we describe our proposed scheme based on iQdyne and compare it with the conventional XY8 measurement. The XY8 sequence is a well-studied quantum sensing scheme, a sort of DD technique, providing excellent sensitivity for oscillating magnetic signals (Fig. 1a). The DD sequence generally consists of spin-flip pulses enabling to decouple the sensor spin from the environment and coupling to the environment selectively. The XY8 sequence consists of equally-spaced MW pulses, which detects an oscillating signal if the half of the oscillation period matches the pulse cycle. The iQdyne protocol is a quantum-classical hybrid measurement protocol comprising sequential measurements using DD with a precisely disciplined interval (Fig. 1a). Splicing the Fourier spectrum at the signal frequency and moving to the IQ diagram, the magnitude of the spectrum depends only on the signal magnitude, and the phase of spectrum directly accords with a signal phase (Fig. 1b).

Figure 1

(a) Measurement protocol with sequential measurements recording the signal in undersampling. Each green represents quantum sensing iterated by times and a CCD readout. (b) Principle of an oscillating field analysis involving the three-dimensional Fourier spectrum: frequency axis, real, and imaginary part of the Fourier coefficient. Splicing out at the signal frequency and moving to the IQ diagram, the Fourier spectrum represents the magnitude and the phase of the oscillating signal. (c) Overview schematic of experimental setup.

At first, we formulated the detection of an oscillating magnetic field using XY8 and pointed out inherent challenges. Let the magnetic field oscillation obey a cosine, and assuming the time at the origin at the end of the former half pulse,where , , and denote the magnitude, initial phase, and frequency of the oscillating field, respectively. Assuming that the interpulse spacing is half the reciprocal of the signal frequency and ignoring the finite length effect of the MW pulses, the electron spin phase acquired is expressed asafter the XY8 sequence[25,29], where , , and represent the gyromagnetic ratio of the electron spin of the NV center, the number of pulses in an XY8, and the proportional coefficient, respectively. The spin state is initialized along the x-axis of the Bloch sphere and read out by rotating around the x-axis, with the resulting fluorescence of XY8 given aswhere and C represent the average count and maximal fluorescence amplitude of the NV center, respectively. The resulting XY8 sequence signals intricately depend on the magnitude and phase of the oscillating magnetic field. In a working condition of , it is completely insensitive for a phase shift but most sensitive for magnitude, and vice versa in . Therefore, XY8 requires a careful pre-adjustment of the working condition depending on the parameter of interest. Moreover, if the two parameters change simultaneously, distinguishing their effect becomes challenging.

Next, we show that the dependence of the iQdyne signal on the parameters is simple. Assuming as an interval between sequential measurements, the oscillating signal phase advances with from one measurement to the next. Since the sampling interval is longer than the oscillation period , the resulting signal includes a low-frequency component due to Nyquist folding. Let be the apparent frequency of the oscillation, another notation of the advancing phase isThen, an initial phase of the oscillating signal of the mth measurement is . The iQdyne result of the mth measurement is expressed as Using iQdyne protocol, the oscillating signal at appears as several oscillating signals with the apparent frequencies given by , the magnitudes by , and the initial phases by , where is the first-kind Bessel function. These are directly observable using Fourier transform, and the derivation details are provided in the method.

Unlike XY8, the iQdyne protocol involves two output parameters, which the signal magnitude depends only on the oscillation magnitude , and the phase depends only on the initial phase , i.e., two parameters are orthogonal. The resulting phase corresponds linearly to the oscillating phase , and its susceptibility is independent of the oscillation magnitude. The iQdyne protocol provides a simple and simultaneous measurement of magnitude and phase. Furthermore, this means that the iQdyne protocol does not need a pre-adjustment of the working condition, since the iQdyne measurement provides the magnitude and phase concurrently, ensuring interpretation is possible under any working condition.

We used an iQdyne protocol on a charge-coupled device (CCD) camera setup (Fig. 1c). Each DD measurement of iQdyne was times iterated DD sequences and the CCD readout (Fig. 1a), adjusting the iteration period to lock a multiple of the signal period. This modification enhances sensitivity like the lock-in amplifier and compensates the overhead due to the long readout time of the CCD, but the principle remains unchanged. We have reported details of this protocol and its frequency characteristics elsewhere[26].

Result

Oscillating field measurement

In this section, demonstration of stroboscopic measurement using iQdyne on a wide-field microscope setup is described. The sensor involved a shallow and dense NV center ensemble on a IIa (100) diamond substrate (See Fig. 1c and sample #1 in the method). A laser excitation via epi- and Köhler-illumination initializes and reads the spin states of the electrons, while a uniform MW radiation via a micro loop coil controls the electron spins. An oscillation signal as a sensing target at (oscillation period ) was generated from a current path created on the surface by photolithography. Each DD protocol involved the pulsed XY8 sequence. The sensing conditions were as follows: measurement interval , number of iterations , sampling interval , and measurement length . Under these conditions, the target signal appears at .

The Qdyne response of the first peak magnitude against the applied field magnitude is displayed in Fig. 2a. The inset shows a typical spectrum measured with . The black dots represents the magnitude at the first peak. As derived earlier, this behavior is consistent with the Bessel function , assuming the oscillating field is proportional to an applied voltage (red curve). We swept the oscillation phase (Fig. 2b) at a working condition of (). Since the proportional coefficient between the assigned initial and measured phases is exactly one, this also agrees with the theoretical description.

Figure 2

The response of iQdyne-based measurement for an oscillating magnetic field. (a) The magnitude of the NV center fluorescence of the first peak (black dots) obeys the first-kind Bessel function (red) with the oscillation field strength. Inset: a measured spectrum with . (b) The measured phase of the first peak corresponds one-to-one to the oscillation phase, with an experimental condition . Inset: plotted on an IQ diagram.

Wide-field imaging of an oscillating field along a current path

For stroboscopic imaging of the magnitude and phase of an oscillating magnetic field, we designed another wide-field microscope setup (Fig. 3a). The sensor was a shallow and dense NV center ensemble on a IIa (111) diamond substrate. The NV center was fabricated by ion implantation with acceleration and dose, creating an approximately -thick layer (see sample #2 in the method). We created a situation with different magnitude and phase of the oscillating field () in an approximately observation area. The oscillating field was generated from a meandering path deposited on the substrate by photolithography (center of Fig. 3a). The fluorescence of the NV center (right of Fig. 3a) corresponds to an intensity distribution of a green laser via a lens array. The thin sensor layer involved an NV ensemble oriented perpendicular to the substrate surface chosen by selective MW, so that measured signals represent a perpendicular component of the magnetic field vector to the surface. According to Ampère’s law, the phase measured at the left side of the path should be inverted compared with the right side. Each DD protocol involved the pulsed XY8 sequence, with sensing conditions as follows: measurement interval , number of iterations , sampling interval , and measurement length . A typical spectrum displaying the frequency versus magnitude is shown in Fig. 3b, and under the stated conditions, the target signal appears at .

Figure 3

Simultaneous wide-field imaging experiment involving (a) the sample setup with the MW resonator to the left, the optical image of the current path as the signal source in the center, and fluorescence image to the right, (b) typical measured spectrum showing frequency versus magnitude, (c) the magnitude map, and (d) phase map of the first peak at each pixel. Arrows in (c) and (d) indicate zero-magnitude lines and phase inversion boundaries corresponding to the zero of Bessel function.

Maps of magnitude and phase of the first peak, respectively, are shown in Fig. 3c,d. Although magnitude and phase have simultaneously and spatially distributed, our iQdyne protocol can distinguish their changes independently. Higher magnitudes appear near the current path because the induced magnetic field obeys Ampère’s law. Comparing the right and left sides of the current path, the measured phases are inverted due to Ampère’s law. The inside of the meandering corner exhibits constructive interference of the magnetic field, with the zero-magnitude points (black arrows in Fig. 3c) corresponding to zero of the Bessel function. This zero line coincides with the boundary of phase inversion (white arrows in Fig. 3d) . Beyond the zero line, the phase is also inverted because the strength of the magnitude suffices for the Bessel function to be minus.

Discussion

Considering its application for spectroscopy and testing, small phase shift detection by protocol is crucial. In this section, we discuss and define the noise floor for phase sensing, namely, precision of the phase shift. Figure 4a shows the typical iQdyne spectrum in Fig. 3b plotted on an IQ diagram. The red point represents the measured Fourier coefficient of the first peak while the semitransparent circles around the origin are Fourier coefficients except for this peak, corresponding to a noise floor. Assuming the noises contained in the measurements is Gaussian, its distribution on the IQ diagram is a two-dimensional (2D) Gaussian distribution with the origin as the mean. The probability distribution of the noise’s magnitude is a Rayleigh distribution[30], and one of the phases is the uniform distribution on the interval . This forms a circle with radius as the standard deviation of the 2D Gaussian representing the precision of the data point. Then, we define the precision of the phase as a projected angle of this noise circle viewed from the data point (Fig. 4b). Through a geometric consideration, the precision is formulated aswhere and r are the distance between the origin and the data point and the radius of the noise circle, respectively. This geometric definition is also justified through statistical estimation theory (see “Methods” section).

Figure 4

(a) Measured Fourier coefficient of the first peak plotted in an IQ diagram, with the inset representing a histogram of noise magnitudes and (b) definition of the phase precision as a projection angle of noise around the origin from a measured point.

For the data in Fig. 4a, the phase precision is , producing a measured phase of , that was acquired at , in within a active area corresponding to a single pixel. This phase precision means a noise floor of the out-of-phase magnetic field (imaginary part or I-axis) with .

Moreover, we noted the robustness of phase sensing. The magnitude map (Fig. 3c) was affected by the intensity distribution of laser excitation (Fig. 3a, right), because the magnitude corresponded to the fluorescence intensity of NV centers. However, the phase map (Fig. 3d) indicated no effect from the excitation inhomogeneity. The reason is that the Fourier transform extracts the phase from the time-development of the fluorescence. This suggests an interesting approach; the phase sensing by our protocol composes a robust measurement with a signal-to-phase encoding translating some signal of interest as in the oscillation phase. Developing a concrete method for such encoding is a topic of future work.

Conclusions and outlook

We propose a new measurement protocol for an oscillating magnetic field based on iQdyne. The approach enables simultaneous detection of magnitude and phase in a wide-field setup without pre-adjusting the measurement condition. We verified that our protocol provides an orthogonal and readily interpretable measurement of magnitude and phase of the oscillation. Furthermore, we demonstrated an oscillating magnetic field imaging around a current path involving concurrent distribution of magnitude and phase with an approximately wide-field observation area. We indicated that the single-pixel phase precision was , corresponding to a -squared region.

The simultaneous measurement of the magnitude and phase of an oscillating magnetic field with high spatial resolution is useful, paving the way to localized current distribution microscopy. In particular, we note that the high spatial resolution of the NV center is suitable for a phenomenon that is zero in sum, but locally non-zero like a vortex[15]. Such a technique has significant applications for Eddy-current inspection, investigating materials or devices by high-frequency impedance imaging, current distribution on topological or 2D materials, and tiny scale inductive testing, which remain elusive with the existing methods.

Methods

Sample preparation

We fabricated two high-purity diamond samples named as sample #1 and #2. Both samples are single crystalline diamond substrates including natural abundance of atoms. Ensemble NV layers were produced by nitrogen ion implantation and subsequent annealing in vacuum, with uniform spatial distribution attributed to the larger ion beam relative to the substrate. The metal pattern on the surface is a -thick Ti/Cu/Au fabricated by photolithography and electron beam deposition. The patterns were connected via gold wire bonding.

Sample #1 has a (100) top facet, with ions implanted at energy and dose at . Annealing was implemented at for . The metal pattern includes -shaped MW loop antennas and linear patterns for MW and RF radiation. Sample #1 and MW antennas are the same as our previous work[26]. The bias magnetic field was around .

Sample #2 has (111) top facet, with ions implanted at energy and dose at . Annealing was implemented at for 2 h. The metal pattern includes meander and linear patterns for RF radiation. The MW control for the sample #2 experiment was applied via a planar resonator with the diamond substrate mounted. This antenna used is a modified version of that reported in our previous work[21]. The MW distribution is almost uniform, but the MW strength is slightly stronger at points close to the current path due to coupling between the current path and MW resonator. The bias magnetic field was around .

Detail of experimental setup

The excitation light is generated from an optically pumped semiconductor laser (Coherent, Verdi G5) and chopped by an acousto-optic-modulator (Gooch & Housego, 35250-.2-.53-XQ). The laser illumination involves epi- and Köhler-illumination optics for initialization and readout of electron spins. The laser is focused onto the back focal plane of an objective lens to enlarge the excitation area. The lens array reduces the spatial coherence and mitigates the interference. The laser passes through an objective lens and fluorescence from the NV center is collected by the same objective and detected by the EMCCD camera (Andor, iXon3 860). Most of the sampling interval accounts for the readout process (approximately ) of CCD pixels, which is a main factor of the iQdyne bandwidth. The objectives are a air objective (Olympus, PlanApoN 60XO) for #1 and a oil objective (Olympus, MPLFLN20X) for #2. The MW pulses are generated by a signal generator with a quadrature (IQ) modulation (Keysight, N5182B), switched by Mini-Circuits, ZASWA-2-50DR+, and amplified by R&K, CGA701M602-4444R. The target signals are generated by a function generator (NF Corporation, WF1967). All timings are controlled by a data timing generator (Tektronix, DTG5274). For the Qdyne-type protocol, stabilizing the sampling interval is crucial. A frequency standard (Stanford Research Systems, FS725) disciplines clocks of all instruments.

Derivation of dependencies of XY8 and Qdyne on parameters

In this section, we formulate the dependencies of each measurement protocol on the magnitude and initial phase . Let C and be the average fluorescence intensity of an NV center and the amplitude of intensity, respectively, these are then expressed as where is the fluorescence intensity when the spin state of the NV center is . We selected the readout axis to be the same direction as the initial superposition state, with the resulting intensity of the XY8 protocol given aswhere is the accumulated phase on the spin state[25]. Equally space dynamical decouplings including XY8 consist of interrogation time , flip pulse, interrogation time, and their repetitions: . Assuming a target field is an oscillating signal , XY8 provides high sensitivity when its interpulse spacing equals a half of the oscillation cycle (). Every pulse makes a virtual flip of a signal which the electron is detecting. For example, the electron spin senses a magnetic field in and in with a MW pulse at . The accumulation by -pulsed XY8 sequence is expressed as with as the proportional coefficient. Hence, the resulting intensity of XY8 is given as:Their susceptibilities are as follows Next, we derive the dependency of Qdyne by assuming M points measurement with a sampling interval , with a series of measurements recorded at timings . During a sampling interval from to , the oscillating signal phase advances byThen, the initial phase of the oscillating signal of the mth measurement is . The resulting intensity of the mth Qdyne measurement is obtained from Here, we used the following relationThe Qdyne signal is expressed as the sum of odd order components of the signal frequency. We can write the th term using the magnitude , phase , and frequency asSince equivalent parameter sets like exist, we impose constraints and . Therefore, the sign of the phase is determined by the magnetic field and the signal frequency . Alternatively, the sign of is determined by two factors: the number of times that crosses zero of the Bessel function and the number of times that the Nyquist folding occurs. Eventually, the Qdyne signal is given as followsSo, we obtain and . The resulting magnitude depends only on the magnetic field , whereas the resulting phase relies only on the initial phase , i.e., the Qdyne protocol achieves an orthogonalized measurement regarding the magnitude and phase of an oscillating signal.

Finaly, each frequency component appears as Nyquist-folded signal since the sampling rate is slower than the frequencies (See the next section).

Nyquist folding and apparent frequency

In the undersampling condition when the sampling interval is longer than half of the target frequency, the oscillating signal with the frequency and initial phase resemble another oscillating signal with an apparent frequency and an apparent phase due to Nyquist folding. At first, the phase advancing is congruent with modulo with the relationship expressed as:In the Fourier spectrum, the frequency axis spans , but the apparent frequency is in . Considering the negative frequency, we obtain a modified phase advancing given by where is a floor function and represents the fractional part of x. This produced the apparent frequency asConsidering a case where the frequency increases from 0, Nyquist folding occurs under the condition that the product is an integer or a half-integer. The apparent phase flips after each folding, and therefore, an apparent phase is obtained from the following

Definition of phase precision

In the main text, we define phase precision by a geometric and intuitive illustration. In this section, we justify this definition through the statistical estimation theory.

Our measurement scheme involves estimation of two parameters, radius r and angle from a complex observed value Z with noise. Let and be the real and imaginary parts, respectively. Assuming obey a normal distribution with a variance independently, the probability distribution functions are as follows: Then, a likelihood function and Fisher information matrix are obtained as: According to the Cramèr–Rao inequality, the standard deviation of any unbiased estimator of is bounded by the reciprocal of the square root of the Fisher information, expressed as:Meanwhile, in this situation, our definition of the phase precision is:Therefore, our definition is above the Cramèr–Rao bound and equals the bound asymptotically with a small noise .