Experimental protection of quantum coherence by using a phase-tunable image drive.

The protection of quantum coherence is essential for building a practical quantum computer able to manipulate, store and read quantum information with a high degree of fidelity. Recently, it has been proposed to increase the operation time of a qubit by means of strong pulses to achieve a dynamical decoupling of the qubit from its environment. We propose and demonstrate a simple and highly efficient alternative route based on Floquet modes, which increases the Rabi decay time ([Formula: see text]) in a number of materials with different spin Hamiltonians and environments. We demonstrate the regime [Formula: see text] with [Formula: see text] the relaxation time, thus providing a route for spin qubits and spin ensembles to be used in quantum information processing and storage.

Introduction

In open quantum systems, coherence of spin qubits is limited by spin-spin interactions, spin diffusion, inhomogeneity of the static and microwave fields[1] as well as charge noise[2]. An increase in coherence time is achieved by dynamically decoupling (DD) qubits from their surroundings using distinct Electron Spin Resonance (ESR) pulses[3-7]. However, such pulses have inherent imperfections and fluctuations, thus requiring their own layer of DD, resulting in a doubly dressed qubit. The technique of concatenated DD[8,9] has been proposed for nitrogen vacancy (NV) centers up to the second order of dressing[8,10-12]. Here we demonstrate a pulse protocol based on Floquet modes which successfully increases the decoherence time, independently of qubit initial state, in a number of materials with different spin Hamiltonians and environments, such as low and high spin-orbit coupling for instance. Rather than focusing on decoupling from the bath by strong excitation, we use very weak pulses and alter the dynamics of the entire system. For short spin relaxation times accessible to our measurement setup (at around 40 K) one can do a direct comparison with the coherence time, and we demonstrate the regime . In magnetic diluted systems , e.g. of several ms in rare earth ions such as YSiO:Er[13] and YSiO:Yb[14] or Si:Bi with a tunable of thousand of seconds[15]. Our general method can thus lead to very long persistent Rabi oscillations, using a single circularly polarized image pulse.

The use of strong continuous microwave excitation has been proposed as a way to protect qubits[16,17] although the quantum gates would need proper redesigning. In related studies, complex pulse design using an arbitrary waveform generator, proved essential in studying Floquet Raman transitions[18,19] and quantum metric of a two-level system[20] in nitrogen vacancy (NV) centers. It is worth noting that in the case of concatenated DD, the frequency of the second order () excitation has to match the Rabi frequency of the first excitation (); also, the two excitations are linearly polarized and perpendicular to each other (the method extends to higher orders in n). Experimentally, the protocol quickly becomes complex and demanding in terms of pulse design and frequency stability, above the second order.

Our protocol uses two coherent microwave pulses: a main pulse drives the qubit Rabi precession while a low-power, circularly polarized (image) pulse continuously sustains the spin motion. The image drive has a frequency close to the main drive and its amplitude is 1-2 orders of magnitude smaller. In this way, a quantum gate could be driven by regular pulses, without the image pulse, while the time interval between gates could be filled with an integer number of Rabi nutations that use our protection protocol. Such scheme would protect the coherence of the qubit in-between usual quantum gates. We note that the initial phase difference between the two pulses allows to tune the spin dynamics by enhancing (or diminishing) the Floquet modes[21] of its second dressing.

The protocol described here can impact quantum sensing of magnetic fields seen as a perturbation of a given Rabi oscillation. Its Fast Fourier Transform (FFT) has a width and the condition significantly improves the field resolution.

The technical implementation is simple and can be generalized to any type of qubit, such as superconducting circuits or spin systems. In this paper we focus on the experimental implementation and simply observe that numerical simulations based on Bloch model with describe the final results very well.

Results and discussion

The standard method to induce Rabi oscillations in a two-level system (TLS) is to apply an electro-magnetic pulse of frequency equal to the TLS level separation (resonance regime, , with the Bohr magneton, the static magnetic field and h the Planck constant ). The pulse will drive the spin population coherently between the two states. Experimentally, the drive is at a frequency (where is the Larmor frequency and is a small detuning away from the resonance condition) followed by read out pulses of frequency to record the state (see Fig. 1a). The method introduced here makes use of two coherent microwave pulses (see Fig. 1b,c): the drive pulse at of amplitude and length creates quantum Rabi oscillations while the second one sustain them using a very low power image of the drive (), operated at . In order to probe at the end of the Rabi sequence, we wait a time longer than , such that , followed by pulses to create a Hahn echo of intensity proportional to (Mn spins are readout without echo, as detailed in Sect. IIB of Supplementary Information). Figure 1b shows one way of coherently creating the drive and its image, by means of a mixer multiplying a pulse at frequency with an intermediate frequency (IF) cosine signal allowing to control the detuning , phase and the pulse length, shape and amplitude.

Figure 1

Comparison of Rabi oscillation with or without the image pulse. (a–c), Schematic of the microwave implementation. (a) The drive () and the readout () sources are independent. (b) Using the same source as the readout (), the drive pulse () is generated by a non-linear mixing with a low frequency signal , with , a process that creates a low-amplitude image drive () as well. The drive, image and readout pulses are coherent and the phase relationship is tunable. (c) Pulse sequence: the drive and the image pulses act at the same time on the qubit while the readout is sensitive to the projection by a spin-echo process. (d–f) Rabi oscillation of the CaWO:Gd, MgO:Mnand P1 defects respectively in the presence (blue) or absence (red) of the image pulse at the optimum condition . The green guideline in (f) shows the improvement of the coherence time when compared to a CPMG pulse sequence while the blue curve shows not only a long coherence but also beatings which are tunable via phase-tunable Floquet dynamics (see text).

Rabi oscillations of three different types of paramagnetic systems—a rare earth ion (Gd), a transition metal ion (Mn) and a defect in diamond (P1)—are shown in Fig. 1d–f, respectively. Their detuned Rabi oscillations induced by the drive pulse only (red curves) are of similar frequency ( MHz) and last for a small number of nutations (). The blue curve shows the Rabi oscillation when the image pulse is superimposed. The oscillations remain intense far beyond the decay time of the red curves and their number is dramatically increased. This effect has maximum impact when the frequency difference between the drive and the image pulses matches the Rabi frequency induced by the main drive, . In addition of the very long coherence time, we observe a slow amplitude modulation, depending on the phase and attributed to Floquet modes, as explained below.

We study the new decay time of the Rabi oscillation under image pumping, by tuning the relaxation time via temperature control and by applying the longest drive pulse available to us (Fig. 2). The length of the drive pulse is limited by the pulse power amplifier of the setup, with a maximum pulse length of 15 s. At 40 K, the relaxation time of the spin system MgO:Mn () is also s. The Rabi oscillation for the transition is shown in Fig. 2. Guidelines showing the exponential decays due to (in green) and (in orange) measured by Hahn echo and inversion recovery respectively, are added as well. While the amplitude of the Floquet mode (slow amplitude modulation) decreases with , the Rabi oscillations persists with a decay time .

Figure 2

Rabi decay under image pumping condition at 40 K. The decay of the Rabi oscillations (green) for a Mn spin is limited by the relaxation time (orange) when an image pulse is used to sustain the dynamics. The inset shows the long-lived oscillations, all the way up to the 15 s limit as imposed by the power amplifier of the setup.

Protection of quantum coherence for different initial states

We analyzed our protocol for different initial states. For initial and final states along , and we observe that the protection protocol is effective for times up to the maximum amplifier gate length of 15 s.

Pulse sequence

The protocol proposed here was tested for different initial states and spin systems. In the following, the details of the experiment are given. The ground state of the spin is along axis and it is used as initial state. In Fig. 3, this is shown in orange and labeled “preparation”. The state is obtained by thermalization. States and are obtained from using a hard pulse around the y or x axes, respectively, able to excite the whole spin ensemble (or spectroscopic line). In MgO:Mnand P1, this is ensured by their narrow spectroscopic linewidths , while in CaWO:Gd, is of the same order of magnitude with the maximum excitation bandwidth. Therefore, some Gd spins of the spin packet might not be in a perfect or state. However, we didn’t notice any particular effect in the final results for Gd. We note that by combining rotations around any of the x, y and z axes, we can prepare the initial state in any position. Once prepared, the spin state is subjected to coherence protected Rabi protocol or to the usual Rabi drive for an integer number of Rabi flops (top and bottom “burst” panels in Fig. 3, respectively). Thus, the final state is along the same direction as the initial one[22]. The echo-based measurement for an initial state is shown in the green panel (with ). For initial states and (blue panel), one have to wait a time in order to let the spin packet defocus before applying a pulse and detect the subsequent echo signal.

Figure 3

Pulse sequence used for different initial states. First, the spin is prepared in the ground state or, after a , in the or states. A Rabi oscillation is performed during for an integer number of Rabi flops, with or without the image pulse (top and bottom panels respectively). The final state is along the same direction as the initial one, and the readout is done using an usual echo detection (blue and green areas).

Throughout this article, the and pulse lengths are 14 and 28 ns, respectively. The readout waiting times are approximately 5 s for P1 defects, 6 s for Mn and 10 s for Gd.

P1 defects in diamond

In the case of P1 defects in diamond, we present measurements of the spin echo signal when the initial and final states are along , with and without the protection protocol for different lengths of the Rabi pulse . Without image pulse (Fig. 4a) the signal is visible after a Rabi burst of 300 ns, but it is rapidly lost for times longer than 1 s. With the protection protocol in place (Fig. 4b), the signal is almost entirely conserved for times up to 15 s. This behavior shows a significant improvement over the CPMG method (see Supplementary Information 1,  s). For this experiment the temperature was set to 15 K and  ns.

Figure 4

P1 defects in diamond. Echo signal for initial and final state along measured after a Rabi burst of various length. Without image pulse (a) the coherence is quickly lost after 0.30 s as shown by an absence of signal above 1 s. With protection protocol in place (b), the signal is well conserved up to maximum burst length of 15 s.

CaWO

:Gd For Gd spins, we measured the spin-echo signal after an integer number of Rabi flops for initial/final states along the , and axes. The Rabi oscillations are shown in Fig. 5a as the real and imaginary part of the recorded signal (blue and orange lines, respectively). One notes that for the initial state along and the Rabi signal starts and ends at zero value, with the end time being indicated by the arrow in each inset. Figure 5b shows large echo signals after 10 s, much longer than  s (see Supplementary Information 1), for each of the initial state preparation. Using a combination of rotations around the x,y and z axes, we can create any initial state for the spin and thus use the image pulse protection for any arbitrary state. The experimental conditions for these measurements are: temperature  K and  ns.

Figure 5

CaWO:Gd. (a) Rabi oscillations shown as the real and imaginary part of the recorded signal (blue and orange lines, respectively). From left to right, the initial (and final) state of the spin are , and . After a time given by an integer number of Rabi flops, shown with an arrow in the inset, (b) the spin-echo signal is recorded for each spin component, using the measurement sequence described in Fig. 3.

Qubit dynamics

The qubit dynamics in the absence of a bath, is described by the spin Hamiltonian in the laboratory frame (see Supplementary Information 1):where is the Larmor frequency caused by the static field, and are the microwave drive and image field, respectively, , is a tunable phase (see Fig. 1) and is a small additional phase, possibly created by imperfections of the setup (as discussed in Supplementary Information 1). Variables and are expressed in units of MHz with . After using the rotating wave approximation (RWA) for a rotation with , the Hamiltonian (1) becomes:When the image field is absent, the Eq. (2) has no explicit time dependence and the Rabi frequency is simply . When is present, the dynamics of can be solved numerically, as it is shown in Fig. 6 for the case of CaWO:Gd.

Figure 6

Rabi dynamics as a function of detuning at fixed drive level . (a) Contour plot showing as a function of pulse length and microwave detuning . At (dashed), the Rabi oscillations are optimally sustained by the image pulse. (b) The FFT of (a) shows an intense peak where the optimum protection is achieved (dashed lines). (c) A simulated FFT using the rotating wave approximation and the parameters of (b), is in agreement with the experimental observations.

For a fixed power of the drive pulse, Rabi oscillations are measured as a function of the detuning . As shown in the contour plot of Fig. 6a, vanishes after few oscillations except when the condition is met. At this Floquet resonance, keeps oscillating for a very long time (s). Its Fast Fourier Transform (FFT) is presented in Fig. 6b. The free (unprotected) Rabi oscillations mode is rather weak and broad showing the large damping caused by the environment. However, when the mode crosses the frequency of the image pulse, indicated by the vertical white dashed line, the peak becomes intense and narrower, as the qubit protection from environment is activated. The condition (or ) gives the most efficient protection of the Rabi oscillation (see Supplementary Information 1).

The general condition is showing a comensurate motion of the qubit and on the Bloch sphere. In other words, the qubit and the image field share a synchronous dynamics generated by the torques of the two drives (see also Sec.III D of the Supplementary Information 1).

We can compare the experimental result to the model described by the Hamiltonian (2) which can be rewritten as in Eq. (S14) in Supplementary Information : . When the “image” pulse is not applied, the Hamiltonian is time independent and the propagator is simply the matrix exponential of the Hamiltonian: . When the image pulse is present (), the Hamiltonian becomes explicitly time-dependent. Although a second canonical transformation RWA could remove the time dependence if , it is importat to leave as a free parameter since the methods works at resonance as well (). Thus, for the sake of generality, we solved numerically the explicit time-dependent differential equations using QuTIP[23]. The parameters used in the simulation have been measured independently: the microwave drive field has been calibrated using the frequency of Rabi oscillations at no detuning (), the image drive was measured by a spectrum analyzer directly connected to the output of the AWG (), relaxation () and decoherence () times were measured by inversion recovery and Carr-Purcell-Meiboom-Gill (CPMG) protocol, respectively (see Supplementary Information 1). We used QuTiP implementation of Lindblad’s master equation with as collapse operator which is equivalent to the phenomenological Bloch model for the case . Figure 6c shows the FFT of computed using the time evolution of . The Hamiltonian (2) describes very well the protection of the coherence by means of the image pulse. Note the existence of a Floquet mode at =7.5 MHz of frequency  MHz, visible in both the experimental and theoretical contour plots of Fig. 6.

The Floquet mode appears as beatings of the Rabi frequency and is -tunable. Similarly to the case of Gd, the qubit protection and Floquet mode dynamics is obtained for the spin of MgO:Mn, here measured under the experimental conditions of Fig. 2. Rabi oscillations and corresponding FFT spectra are shown in Fig. 7 for two values of the initial phase, (green) and (gold), while simulations are shown in black. The decay times are much larger than for both values (here  s, see Fig. 2); however, the dynamics is strikingly different. When the drive and image pulses have the same initial phase (one can consider the initial time in Eq. (2) as without loss of generality), the Rabi oscillations have maximum visibility, with almost no beatings (see Supplementary Information 1). At , the spin torques generated by the and fields induce strong beatings or a Floquet mode creating two additional modes of the Rabi frequency. The left panel shows Rabi splittings equal to the Floquet frequency for and a single Rabi oscillation for .

Figure 7

Phase-tunable Floquet dynamics. Rabi oscillations of Mn at 40 K for two values of . For (green symbols), and pulses have the same initial phases resulting in optimally sustained Rabi oscillation. For (orange symbols), strong beatings are observed ( terms of are ). The solid lines are simulations using our model with no fitting parameters. The right panel shows the experimental FFT traces. The green one has a well-defined peak at while in the other one two peaks (top inset) are separated by twice the Floquet mode frequency (bottom inset).

Experimentally, we can continuously vary the value of and analyze the frequency and intensity of the Floquet mode. As an example, a comparison between theory and experiment is shown in Fig. 8 for the case of CaWO:Gdfor  MHz. For even and odd multiples of , a single and a splitted Rabi mode is observed, respectively. The Rabi splitting is the Floquet mode and is constant as a function of but its intensity oscillates with a period of . The effect is evident in simulations as well, since the terms in of are along the same direction or orthogonal, for and 45 respectively.

Figure 8

Phase-dependence of the Floquet mode. Experimental (top) and simulated (bottom) FFT of Rabi oscillations for  MHz as a function of (amplitude is normalized to the highest peak). For , Rabi oscillations have a single mode while for a splitting of twice the Floquet mode frequency is observed. While its frequency is fixed, the intensity of the Floquet mode changes gradually, with a period of .

In regards to decoherence sources, it is safe to assume that the main contribution comes from the spin bath made by nuclear spins surrounding the central spin, as well as other electronic spins located in its closed vicinity. Such scenario is the typical situation in spin systems operated at low enough temperatures to reduce the role of the phonon bath on . The details of the entangled qubit-bath dynamics are outside the scope of the current study. We do observe that the final results are well described when only dissipation is the source of decoherence, leading to . This may indicate that the image pulse is able to control the dynamics and thus the decoherence of the spin bath.

The qubit rotation is thus tunable by using a pre-selected value of , allowing to create complex rotations. With a decoherence time approaching spin lifetime , the value of can be changed while qubit control is still ongoing. Our study demonstrates a sustained quantum coherence using a general protocol that can be readily implemented to any type of qubit. Our approach can be used in other detection schemes, such as sensitive spin detection using on-chip resonance techniques[24-26].

While preparing our manuscript for resubmission following initial refereeing, a  related protocol applied to NV centers has been published[27].

Methods

Spectrometer setup

The measurements have been performed on a conventional pulse ESR spectrometer Bruker E680 equipped with an incoherent electron double resonance (ELDOR) bridge and a coherent arbitrary waveform generator (AWG) bridge. In the ELDOR bridge (Fig. 1a), the drive and the read out pulses come from two independent sources while with the AWG bridge (Fig.1b) all the pulses are generated using the same microwave source and thus they are all phase coherent. The drive frequency is generated by mixing the source (used as a local oscillator) with a low frequency and phase controllable signal through an in-phase quadrature (IQ) mixer. Ideally, the output of the mixer is monochromatic with the frequency . In reality, the output consists of a principal frequency (the drive) and of lower amplitude images (see Supplementary Information 1 for more information). Since the effect of the image is the central part of this paper, we have characterized the AWG bridge using a spectrum analyzer, right before the power amplification stage. An example of spectrum is presented in Fig. S6 of the Supplementary Information. The power of the image is lower by  dB than . Consequently, an amplitude ratio of the MW magnetic fields around 0.12 is used in simulations.

Pulse sequence

First, the system is set to be in resonance condition . The drive pulse of amplitude , frequency and length induces Rabi oscillation in detuning regime. At the same moment the image pulse (generated through the IQ mixer) of amplitude , frequency and the same length also irradiate the spins. In order to probe at the end of the Rabi sequence, we wait a time longer than , such that , followed by pulses to create a Hahn echo of intensity proportional to .

Spin systems

The methodology presented here is demonstrated on different spin systems: the nitrogen substitution in diamond P1 defect () (concentration :100 ppm), Mn impurities in MgO ()[28-31] with a concentration of 10 ppm and Gd impurities in CaWO ()[24,32] with a concentration of 50 ppm. Despite the large Hilbert space of the Mn and Gd spin Hamiltonians, the orientation of the magnetic field and the frequency and power of the microwave excitation are chosen to avoid multiple level transitions and thus select only one resonance[28]. Therefore, the spin systems can be considered as effective two-level systems undergoing coherent Rabi rotations. The spin Hamiltonians, operating parameters (fields and frequencies) as well as characteristic times for these materials are given in the Supplementary Information 1.

Supplementary Information.