Narrow-bandwidth sensing of high-frequency fields with continuous dynamical decoupling

TL;DR

This paper introduces a method to extend high sensitivity magnetic field sensing to high frequencies using continuous dynamical decoupling, demonstrated with nitrogen-vacancy centers in diamond, achieving detection at 1.6 GHz.

Contribution

The authors develop a general scheme integrating dynamical decoupling into high-frequency sensing, extending low-frequency techniques to the GHz range.

Findings

Achieved coherence times up to 1.43 ms in diamond NV centers.

Detected magnetic fields as weak as 4 nT at 1.6 GHz.

Observed increased coherence time due to the signal itself.

Abstract

State-of-the-art methods for sensing weak AC fields are only efficient in the low frequency domain (< 10 MHz). The inefficiency of sensing high frequency signals is due to the lack of ability to use dynamical decoupling. In this paper we show that dynamical decoupling can be incorporated into high frequency sensing schemes and by this we demonstrate that the high sensitivity achieved for low frequency can be extended to the whole spectrum. While our scheme is general and suitable to a variety of atomic and solid-state systems, we experimentally demonstrate it with the nitrogen-vacancy center in diamond. For a diamond with natural abundance of $^{13}$C we achieve coherence times up to 1.43 ms resulting in a smallest detectable magnetic field strength of 4 nT at 1.6 GHz. Attributed to the inherent nature of our scheme, we observe an additional increase in coherence time due to the signal itself.

Equations

$$H$$

$$H_I$$

$$H_{II}=\frac{{\Omega }_2}{4}{\sigma }_y+\frac{g}{4}\left(-i{\sigma }_{+}{e}^{-i\left(\frac{{\Omega }_2}{2}t+\phi \right)}+i{\sigma }_{-}{e}^{+i\left(\frac{{\Omega }_2}{2}t+\phi \right)}\right).$$

$$H_{III}=\frac{g}{8}\left({\sigma }_{+}{e}^{-i\phi }+{\sigma }_{-}{e}^{+i\phi }\right),$$

$$\delta B_{\text{min}}(t,\tau )=\frac{\delta S}{\text{max}\left|\frac{\partial S}{\partial B}\right|}=\frac{1}{{\gamma }_{\text{NV}}}\frac{\sigma (t)}{\alpha \tau C}.$$

Full text

dynamical decoupling, high frequency, quantum sensing, two-level system, nitrogen-vacancy, diamond

††preprint: AIP/123-QED

Improving the sensitivity of high frequency sensing schemes is of great significance, especially for classical fields sensing Chipaux et al. (2015); Kolkowitz et al. (2015); Shao et al. (2016), detection of electron spins in solids Sushkov et al. (2014); Hall et al. (2016) and nuclear magnetic resonance spectroscopy Kimmich and Anoardo (2004). The common method to detect high frequency field components is based on relaxation measurements, where the signal induces an observable effect on the lifetime, T1, of the probe system Sushkov et al. (2014); Schmid-Lorch et al. (2015); Hall et al. (2016). Nevertheless, the sensitivity of this method is limited by the pure dephasing time T${}_{2}^{*}$ of the probe system.

Pulsed dynamical decoupling Hahn (1950); Carr and Purcell (1954); Meiboom and Gill (1958) can substantially increase the coherence time Viola and Lloyd (1998); Biercuk et al. (2009); Du et al. (2009); Lange et al. (2010); Ryan et al. (2010); Naydenov et al. (2011); Wang et al. (2012); Bar-Gill et al. (2013). In order to carry out sensing with a decoupling scheme, the frequency of the decoupling pulses has to be matched with the frequency of the target field Taylor et al. (2008); Kotler et al. (2011). This largely restricts the approach to low frequencies, as the repetitive application of pulses is limited by the maximum available power per pulse Gordon et al. (2008). The same power restrictions are present for very rapid and composite pulse sequences aimed to decrease both external and controller noise Khodjasteh and Lidar (2005); Uhrig (2007); Souza et al. (2011); Yang et al. (2011); Farfurnik et al. (2015).

With continuous dynamical decoupling (CDD) Fanchini et al. (2007); Bermudez et al. (2011, 2012); Cai et al. (2012a); Xu et al. (2012); Golter et al. (2014); Rabl et al. (2009); Gordon et al. (2008); Clausen et al. (2010); Laucht et al. (2017); Cai et al. (2012b); Cohen et al. (2015); Teissier et al. (2017) robustness to external and controller noise can be attained, especially for multi-level systems Baumgart et al. (2016); Timoney et al. (2011); Aharon et al. (2013, 2016). However, the significance of CDD for sensing high frequency fields remained elusive. Indeed, it was unclear whether it is possible to incorporate such a protection into the metrology task of sensing frequencies in the GHz domain. The first step towards this goal was done recently by integrating CDD in the sensing of high frequency fields with three level systems Aharon et al. (2016).

In this article, we propose, analyze and experimentally demonstrate for the first time a sensing scheme that is capable of probing high frequency signals with a coherence time, T2, limited sensitivity. Unlike relaxation measurements comprising a bandwidth $\propto 1/\text{T}_{2}^{*}$, determined by the pure dephasing time, T${}_{2}^{*}$, of the sensor (up to the MHz range), our protocol overcomes the imposed limitation by protecting the addressed two-level system (TLS) with an adapted concatenated CDD approach. We use and adjust it such that high frequency sensing becomes feasible even for not phase-matched signals. As a result, the proposed scheme is generic and works for many atomic or solid state TLS, in which the energy gap matches the frequency of the signal under interrogation. A remarkable feature of our scheme is the fact that the signal to be probed also works partially as a decoupling drive and thus improves further the sensitivity of the sensor.

We demonstrate the performance of continuous dynamical decoupling by applying it to a nitrogen-vacancy (NV) center in diamond with natural abundance of 13C (cf. Fig. 1a). Here, we utilize two of its ground sub-levels as the TLS. The states of the NV center can be read out and initialized by a 532 nm laser, which reveals spin dependent fluorescence between the two levels Jelezko et al. (2004a, b); Balasubramanian et al. (2008). The system can be manipulated by driving it with microwave fields. We show that by using a concatenation of two drives, an improvement in coherence time of the sensor by more than one order of magnitude is achieved. Taking into account the effect of an external signal, $g$, on the sensor during a concatenation of two drive fields, we obtain an improvement in bandwidth for high frequency sensing by three orders of magnitude in comparison to the relaxometry approach. Moreover, we report on the measurement of a weak high frequency signal with strength $g$, which relates to a smallest detectable magnetic field amplitude of $\delta B_{\text{min}}\approx 4$ nT.

The sensing scheme

The basic idea of utilizing concatenated continuous driving to create a robust qubit is illustrated in Fig. 1b. The concatenation of two phase-matched driving fields results in a robust qubit Cai et al. (2012b); Aharon et al. (2016). In what follows we show that such a robust qubit can be utilized as a sensor for frequencies in the range of the qubit’s energy separation and hence, dynamical decoupling can be integrated into the sensing task.

By the concatenated driving, the qubit is prepared in a state that allows for strong coherent coupling to the high-frequency signal to be probed (corresponding to the last TLS in Fig. 1c). In the total Hamiltonian, $H$, we consider the concatenation of two driving fields of strength (the Rabi frequency) $\Omega_{1}$ and $\Omega_{2}$, respectively. The Hamiltonians of the TLS, $H_{0}$, the protecting driving fields, $H_{\Omega_{1}},H_{\Omega_{2}}$ and the signal, $H_{s}$, are given by

[TABLE]

where $\omega_{0}$ is the energy gap of the bare states ($\hbar=1$), $\omega_{s}$ is the frequency of the signal, and $g$ is the signal strength which we want to determine. We tune the system, i.e., $\omega_{0}$, $\Omega_{1}$, and $\Omega_{2}$, such that $\omega_{s}=\omega_{0}+\Omega_{1}+\Omega_{2}/2$.

It is an important feature that phase matching between the signal and the control is not required, which means that the signal phase $\varphi$ can be unknown and moreover, it may vary between experimental runs. In addition, we make the assumption that $\omega_{0}\gg\Omega_{1}\gg\Omega_{2}\gg g$. Moving to the interaction picture (IP) with respect to $H_{0}=\frac{\omega_{0}}{2}\sigma_{z}$ and making the rotating-wave-approximation, we obtain

[TABLE]

This picture incorporates the effect of $\Omega_{1}$ onto a TLS and express the new system in eigenstates of $\sigma_{x}$, the $|\pm\rangle$ (dressed) states, which separates the contributions from $\Omega_{2}$ and $g$. For a large enough drive $\Omega_{1}$, the $|\pm\rangle$ eigenstates are decoupled (in first order) from magnetic noise, $\delta B\sigma_{z}$, because $\langle\pm|\sigma_{z}|\pm\rangle=0$. However, power fluctuations $\delta\Omega_{1}$ of $\Omega_{1}$ limit the coherence time of the dressed states. The resulting IP is illustrated in the second TLS in Fig. 1c.

We continue by moving to a second IP with respect to $H_{01}=\frac{\Omega_{1}}{2}\sigma_{x}$, which leads to

[TABLE]

Once again, we incorporate $\Omega_{2}$ into the dressed states, so that solely the contribution of the signal $g$ becomes obvious, which is depicted in the last TLS of Fig. 1c. The second drive, $\Omega_{2}$, which is larger than $\delta\Omega_{1}$, creates effectively doubly-dressed states (the $\sigma_{y}$ eigenstates). These doubly-dressed states are immune to power fluctuations of $\Omega_{1}$ and hence prolong the coherence time (see Supplementary Sec. II for more details). Moving to the third IP with respect to $H_{02}=\frac{\Omega_{2}}{4}\sigma_{y}$ results in

[TABLE]

where we can clearly see that the signal $g$ induces rotations in the robust qubit subspace (either with $\sigma_{+}$ or $\sigma_{-}$). These rotations are obtained for any value of an arbitrary phase $\varphi$ and the bandwidth ($\propto 1/\text{T}_{2}$) is now limited by the coherence time, T2, of the sensor. Hence, if a given TLS exhibits the possibility of manipulating it via drive fields $H_{\Omega_{1}}$ and $H_{\Omega_{2}}$, we can achieve a high frequency sensor in the range of $\omega_{0}$.

By this, we overcome the low frequency limit that is common to state-of-the-art pulsed dynamical decoupling sensing methods. In addition, we present an analog pulsed version of our scheme, where the pulsing rate is much lower than the frequency of the signal (see Supplementary Sec. IX). However it is not a direct measurement of the signal, but based on a signal demodulation approach. Compared to the pulsed schemes, continuous dynamical decoupling does not suffer from being susceptible to higher harmonics of the decoupling window appearing naturally from the periodic character of the pulsed sequence Loretz et al. (2015). Eventually, less power per unit time is used in the continuous scheme leading to a smaller overall noise contribution from the drive.

Results

After determining the optimal drive parameters, $\Omega_{1}$ and $\Omega_{2}$, for the concatenated sensing sequence, and thereby maximize the coherence times, T${}_{2}^{\Omega_{1}}$ and T${}_{2}^{\Omega_{1},\Omega_{2}}$, respectively, of the sensor (see Supplementary Sec. VI), we apply an external high frequency signal (according to $H_{s}$ in equation (1)) tuned to one of the four appearing energy gaps $\omega_{s}$ of the doubly-dressed states. In these energy gaps an effective population transfer can occur between the states of the robust TLS, evidenced by signal induced Rabi oscillations at a rate $g^{\prime\prime}=g^{\prime}/2=g/4$ in the double drive case (see Supplementary Sec. II).

The measurements take place in the laboratory frame, i.e. all three contributions $\Omega_{1},\Omega_{2}$ and $g$ to the population dynamics of the TLS will be visible. In order to see solely the effect of $g$ on the TLS, we alter the modulation of the second drive in equation (1) to $\cos(\Omega_{1}t+\pi/2)$. This does not change the performance of the scheme, but only changes the axis of rotation to $\sigma_{z}$ for the second drive $\Omega_{2}$. Since the readout laser is effectively projecting the population in the $\sigma_{z}$ eigenbasis, we can make the Rabi rotations of $\Omega_{2}$ invisible to the readout. To remove the effect of $\Omega_{1}$ in the data, we can simply sample the measurement at multiple times of $\tau_{\Omega_{1}}=2\pi/\Omega_{1}$, i.e. we measure at times $t=N\tau_{\Omega_{1}}$ ($N\in\mathbb{N}$). This procedure reveals directly $g^{\prime\prime}$($g^{\prime}$) as the signal induces Rabi oscillations of the robust qubit under double (single) drive (Fig. 2). Alternatively, we could have applied at the end of the drive a correction pulse in order to complete the full $\Omega_{1}$ and $\Omega_{2}$ rotations so that just the effect of $g^{\prime\prime}$ remains.

Without a signal, $g$, we achieve coherence times of T${}_{2}^{\Omega_{1}}\approx 60$ µs with a single drive ($\Omega_{2}=0$, compare Fig. S3 in Sec. VI.A in Supplementary) and T${}_{2}^{\Omega_{1},\Omega_{2}}\approx 393$  µs with a double drive (compare Fig. S4 in Sec. VI.B in Supplementary). The results for long and slow Rabi oscillations induced by an external signal, $g$, under single and double drive ($\Omega_{2}/\Omega_{1}\approx 0.15$) are shown in Fig. 2. These illustrate a significant increase of the coherence time of the sensor by two orders of magnitude, from T${}_{2}^{\Omega_{1}}\approx 60$  µs to a lifetime limited coherence time of (T${}_{1}/2\approx$) T${}_{2}^{\Omega_{1},\Omega_{2},g}\approx 1.43$ ms. It should be noted that the signal itself can be considered as an additional drive (cf. equation 3), correcting external errors $\delta\Omega$ of the previous drive and thereby prolonging the coherence time even further. Consequently, we can improve the bandwidth for high frequency sensing by almost three orders of magnitude from $\sim 900\,$kHz (for T${}_{2}^{*}\approx 1.1$  µs) to $\sim 700$ Hz (for a T${}_{2}^{\Omega_{1},\Omega_{2},g}\approx 1.43$ ms). Moreover, in the Supplementary Sec. III we discuss an improved version of our scheme which has the potential to push the coherence time of the sensor further towards the lifetime limit.

To benchmark the double drive scheme against a standard single drive approach, we determine the smallest magnetic field which can be sensed after an accumulation time $t$. The smallest measurable signal $S$ is eventually bounded by the smallest measurable magnetic field change $\delta B_{\text{min}}$, which is found to be

[TABLE]

Here, $\gamma_{\text{NV}}/2\pi=28.8\,$GHz T*-1* is the gyromagnetic ratio of the NV defect, $\sigma(t)$ is the standard deviation of the measured normalized fluorescence counts after time $t$, $\alpha$ accounts for a different phase accumulation rate depending on the decoupling scheme and $C$ is the contrast of the signal (see Supplementary Sec. V for detailed derivation). Since the photon counting is shot noise limited, we have $\sigma(t)=1/\sqrt{N_{ph}\cdot N}$, with $N_{\text{ph}}$ being the number of photons measured in $\tau$ and $N=t/\tau$ is the number of sequence repetitions. With this, equation  (5) will transform in the commonly known form Itano et al. (1993); Budker and Romalis (2007) with some measurement dependent constants.

We recorded $\sigma(t)$ as a function of time and use this to determine $\delta B_{\text{min}}$. The results of this measurement for both the single and double drive are summarized in Fig. 3. The sensitivity can be obtained by $\eta(\tau)=\delta B_{\text{min}}(t,\tau)\sqrt{t}$, which is optimal in the vicinity of the coherence time of the sensor, $\tau\approx\text{T}_{2}$. With our system, we achieve a sensitivity of $\eta_{\Omega_{1},\Omega_{2},g}\lesssim 1$  µT Hz*-0.5* in the double drive case at $\sim 1.6$ GHz, which should be compared to $\eta_{\Omega_{1},g}\lesssim 20$  µT Hz*-0.5* for a single drive approach.

Both traces in Fig. 3 were recorded while a signal $g$ was applied. Apart from the mere fact, that the number of driving fields are different, the specific choice for $\tau$ will also determine the magnitude of the smallest measurable magnetic field change $\delta B_{\text{min}}$. Obtaining the coherence time without a signal, $g$, (which are T${}_{2}^{\Omega_{1}}$ and T${}_{2}^{\Omega_{1},\Omega_{2}}$), is a common practice in the field, but will not result in a correct choice of $\tau$ for the sensitivity measurement and also for equation (5), since the signal has an impact on the sensor’s sensitivity. However, if $\delta B_{\text{min}}$ shall be evaluated correctly, then the non-linearity, i.e. the coherence time prolonging effect of the signal, has to be taken into account. Otherwise an even worse $\delta B_{min}$ will be measured as it is exemplarily shown for the single drive case in Fig. 3, where $\delta B_{min}$ was evaluated and measured under the naive assumption that the signal has no effect on the coherence time of the sensor (i.e. we measure at $\tau=\text{T}_{2}^{\Omega_{1}}$ and not at $\tau=\text{T}_{2}^{\Omega_{1},g}$). This effect was included in the double drive case.

To examine the signal protection effect more in detail, the coherence time of the sensor is measured as a function of signal strength, $g$, in a single drive configuration (see also Supplementary Sec. VII). From these measurements we project the sensitivity associated with a specific signal strength (Fig. 4), assuming $\sigma(t)$ is unchanged for the same repetition $N$. This is a reasonable assumption given that the only difference between measurements is the signal strength, $g$, and sequence length, $\tau$.

Eventually, this phenomenon, which seems to be an inherent part of this continuous scheme, can be used to further increase the performance of the sensor by fine tuning the controlled parameters (static bias field $B_{\text{bias}}$ and thereby changing $\omega_{0}$, $\Omega_{1}$ and $\Omega_{2}$) to match the signal frequency, $\omega_{s}$, and strength, $g$ (see Supplementary Sec. III for further discussions).

Conclusion

We have demonstrated for the first time that dynamical decoupling can be used in the context of sensing high frequency fields. In contrast to state-of-the-art pulsed dynamical decoupling protocols, we can show that continuous dynamical decoupling can be simultaneously integrated into the sensing task. By utilizing a NV center in diamond we have demonstrated by pure concatenation of two drives a coherence time of $\sim 393$ µs which constitutes an improvement of more than two orders of magnitude over T${}_{2}^{*}$, and an increase of resolution from the MHz to a few kHz. The application of this method for wireless communication Baylis et al. (2014) could have a transformative effect due to the high resolution of the protocol. Since the protocol is applicable to a variety of solid-state, molecular, and atomic systems, we believe that it has a great potential to have a significant impact on many fields and tasks that involve high frequency sensing (up to frequencies in the THz range). Eventually, this method could also be used to improve the coupling to quantum systems Cai et al. (2012a).

Acknowledgements

Acknowledgements.

The experiments presented here were supported by the Qudi Software Suite Binder et al. (2017). A. S., A. H. and U. L. A. acknowledge funding from the Innovation Foundation Denmark through the project EXMAD and the Qubiz center, and the Danish Research Council via the Sapere Aude project (DIMS). T. U. and F.J. acknowledge the Volkswagenstiftung. A. R. acknowledges the support of the Israel Science Foundation (grant no. 1500/13), the support of the European commission, EU Project DIADEMS. This project has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 667192 Hyperdiamond and Research Cooperation Program and DIP program (FO 703/2-1).