Suppression of spin-bath dynamics for improved coherence of multi-spin-qubit systems

Abstract

Multi-qubit systems are crucial for the advancement and application of quantum science. Such systems require maintaining long coherence times while increasing the number of qubits available for coherent manipulation. For solid-state spin systems, qubit coherence is closely related to fundamental questions of many-body spin dynamics. Here we apply a coherent spectroscopic technique to characterize the dynamics of the composite solid-state spin environment of nitrogen-vacancy colour centres in room temperature diamond. We identify a possible new mechanism in diamond for suppression of electronic spin-bath dynamics in the presence of a nuclear spin bath of sufficient concentration. This suppression enhances the efficacy of dynamical decoupling techniques, resulting in increased coherence times for multi-spin-qubit systems, thus paving the way for applications in quantum information, sensing and metrology.

ARTICLE

Received 9 Jan 2012 | Accepted 23 Apr 2012 | Published 22 May 2012 DOI: 10.1038/ncomms1856

N. Bar-Gill1,2, L.M. Pham3, C. Belthangady1, D. Le Sage1, P. Cappellaro4, J.R. Maze5, M.D. Lukin2, A. Yacoby2 & R. Walsworth1,2

Multi-qubit systems are crucial for the advancement and application of quantum science. Such systems require maintaining long coherence times while increasing the number of qubits available for coherent manipulation. For solid-state spin systems, qubit coherence is closely related to fundamental questions of many-body spin dynamics. Here we apply a coherent spectroscopic technique to characterize the dynamics of the composite solid-state spin environment of nitrogen-vacancy colour centres in room temperature diamond. We identify a possible new mechanism in diamond for suppression of electronic spin-bath dynamics in the presence of a nuclear spin bath of sufcient concentration. This suppression enhances the efcacy of dynamical decoupling techniques, resulting in increased coherence times for multi-spin-qubit systems, thus paving the way for applications in quantum information, sensing and metrology.

Suppression of spin-bath dynamics for improved coherence of multi-spin-qubit systems

1 Harvard-Smithsonian Center for Astrophysics, Cambridge, Massachusetts 02138, USA. 2 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA. 3 School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA. 4 Department of Nuclear Science and Engineering, MIT, Cambridge, Massachusetts 02139, USA. 5 Department of Physics, Ponticia Universidad Catolica de Chile, Santiago 7820436, Chile. Correspondence and requests for materials should be addressed to N.B.G. (email: [email protected]) or to

R.W. (email: [email protected]).

NATURE COMMUNICATIONS | 3:858 | DOI: 10.1038/ncomms1856 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1856

Understanding and controlling the coherence of multi- spin-qubit solid-state systems is crucial for quantum information science13, basic research on quantum many-body

dynamics4 and quantum sensing and metrology59. Examples of such systems include nitrogen-vacancy (NV) colour centres in diamond10, phosphorous donors in silicon11 and quantum dots12. In characterizing the potential usefulness of multi-qubit systems for quantum applications, the product of the qubit density (nqb) and the qubit coherence lifetime (T2) serves as a basic gure-of-merit (FOM), FOM = nqbT2. For example, in quantum measurements the phase-shi sensitivity scales as 1/ FOM 5,1315. Increasing this multi-qubit FOM requires an understanding of the sources of decoherence in the system and their interplay with qubit density. In particular, for solid-state spin qubits T2 is typically limited by interaction with an environment (that is, bath) of paramagnetic spin impurities, whereas nqb is limited by fabrication issues, the associated creation of additional spin impurities in the environment, and the ability to polarize (state prepare) the quantum spins.

The paradigm of a central spin coupled to a spin environment has been studied intensively for many years (see for example refs 16 and 17), and quantum control methods have been developed to extend the spin coherence lifetime by reducing the eective interaction with the environment. In particular, dynamical decoupling techniques pioneered in nuclear magnetic resonance have recently been applied successfully to extend the eective T2 of single NV-diamond

electronic spins by more than an order of magnitude1820.

Here we study experimentally the spin environment of NV colour centres in room temperature diamond and their multi- spin-qubit FOM. We apply a spectral decomposition technique2123 to characterize the dynamics of the composite solid-state spin bath, consisting of both electronic spin (N) and nuclear spin (13C) impurities. We study three dierent diamond samples with a wide range of NV densities and impurity spin concentrations (measuring both NV ensembles and single NV centres), and nd unexpectedly long correlation times for the electronic spin baths in two diamond samples with natural abundance (1.1%) of 13C nuclear spin impurities. We identify a possible new mechanism in diamond involving an interplay between the electronic and nuclear spin baths that can explain the observed suppression of electronic spin-bath dynamics. We show that this spin-bath suppression enhances the efficacy of dynamical decoupling techniques, enabling large values for the multi-qubit FOM ~21014 (ms cm 3; compare, for example, to ref. 24, and to the similar FOM obtained with a much higher NV density sample9).

ResultsSpectral decomposition technique. Owing to coupling of the NV spins to their magnetic environment (Fig. 1a,b), coherence is lost over time with the general form C(t) = e (t), where the functional (t) describes the time dependence of the decoherence process (Supplementary Methods). In the presence of a modulation acting on the NV spins (for example, a resonant MW (microwave) pulse sequence), described by a lter function in the frequency domain Ft() (see below), the decoherence functional is given by25,26.

c p w w

w w

13C

V

13C

NV

N

N

Hahn -echo:

2 2

638800 nm

2

2

532 nm

MW

CPMG:

ms = +1 ms = 1

3E

2

2

3A2

Figure 1 | NV centre in diamond and applied spin-control pulse sequences. (a) Lattice structure of diamond with an NV colour centre. The NV electronic spin axis is dened by nitrogen and vacancy sites, in one of four crystallographic directions in the diamond lattice. NV orientation subsets can be spectrally selected by applying a static magnetic eld B0.

(b) Magnetic environment of NV centre electronic spin: 13C nuclear spin impurities and N electronic spin impurities. (c) Hahn-echo and multi-pulse (CPMG) spin-control sequences. (d) Negatively charged NV centre electronic energy level structure. Electronic spin polarization and readout is performed by optical excitation at 532 nm and red uorescence detection. Ground-state spin manipulation is achieved by resonant MW excitation. The ground-state triplet has a zero magnetic eld splitting 2.87 GHz.

technique2123. As seen from equation (1), if an appropriate modulation is applied to the NV spins such that F t

t ( )/( ) = ( )

2

,

that is, if a Dirac -function is localized at a desired frequency 0, then c w p

( ) = ( )/

0

t tS . Therefore, by measuring the time dependence of the qubit coherence C(t) when subjected to such a spectral -function modulation, we can extract the spin baths spectral component at frequency 0:

S C t t

( ) = ( ( ))/ .

0

w p

w w d w w

0

ln

This procedure can then be repeated for dierent values of 0 to provide complete spectral decomposition of the spin environment.

A close approximation to the ideal spectral lter function Ft() described above can be provided by a variation on the well-known

CarrPurcellMeiboomGill (CPMG) pulse sequence for dynamical decoupling of a qubit from its environment27 (Fig. 1c). The CPMG pulse sequence is an extension of the Hahn-echo sequence28 (Fig. 1c), with n equally spaced -pulses applied to the system aer initially rotating it into the x axis with a /2-pulse. We apply a deconvolution procedure to correct for deviations of this lter function from the ideal Dirac -function (Supplementary Fig. S1 and Supplementary Methods).

Spectral function of a spin bath. The composite solid-state spin environment in diamond is dominated by a bath of uctuating N electronic spin (S = 1/2) impurities, which causes decoherence of the probed NV electron-spin qubits through magnetic dipolar interactions. In the regime of low external magnetic elds and room temperature (relevant to the present experiments), the N bath spins

(2)(2)

( ) = 1 ( ) ( ),

0 2

t d S Ft

(1)(1)

where S() is the spectral function describing the coupling of the system to the environment. Equation (1) holds in the approximation of weak coupling of the NV spins to the environment, which is appropriate for systems with (dominantly) electronic spin baths17, as is the case with the diamond samples discussed here (see below).

S() can be determined from straightforward decoherence measurements of the NV spin qubits using a spectral decomposition

NATURE COMMUNICATIONS | 3:858 | DOI: 10.1038/ncomms1856 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1856

ARTICLE

are randomly oriented, and their ip-ops (spin-state exchanges) can be considered as random uncorrelated events17. Therefore, the resulting spectrum of the N baths coupling to the NV spins can be assumed to be Lorentzian25:

S( ) = 1

2 1 ( ) .

2

Table 1 | Comparison of key characteristics and extracted average-t Lorentzian spin-bath parameters for the NV-diamond samples studied in this work.

12C Apollo HPHT

Measured technique Wide eld Wide eld Confocal N Concentration, p.p.m. ~1 ~100 ~50NV density, cm 3 ~1014 ~1016 ~1012

13C Concentration, % 0.01 1.1 1.1 Echo (1-pulse) T2, s 240 (6) 2 (1) 5 (1)

T2 scaling n0.43(6) n0.65(5) n0.7 (1) Max. achieved T2, s 2.2 20 35 (expected) 60 kHz 6 MHz 3 MHz (measured) 30 (10) kHz 7 (3) MHz 1 (1) MHz c (expected), s 15 0.17 0.34c (measured), s 10 (5) 3 (2) 10 (5)

FOM (ms cm3) 21014 21014 1010

w t

p wt

c

(3)(3)

+

c

This spin-bath spectrum is characterized by two parameters: is the average coupling strength of the N bath to the probed NV spins, and c is the correlation time of the N bath spins with each other, which is related to their characteristic ip-op time. In general, the coupling strength is given by the average dipolar interaction energy between the bath spins and the NV spins, and the correlation time c is given by the inverse of the dipolar interaction energy between neighbouring bath-spins. As such spinspin interactions scale as 1/r3, where r is the distance between spins, it is expected that the coupling strength scales as the N bath spin density nspin (that is, nspin), and the correlation time scales as the inverse of this

density (that is, cnspin).

Note also that the multi-pulse CPMG sequence used in the spectral decomposition technique extends the NV spin coherence lifetime by suppressing the time-averaged coupling to the uctuating spin environment. In general, the coherence lifetime T2 increases with the number of pulses n used in the CPMG sequence. For a Lorentzian bath, in the limit of very short correlation times (cT2), the sequence is inefficient and T2n0 (no improvement with number of pulses). In the opposite limit of very long correlation times c T2, the scaling is T2n2/3 (refs 2931; see also recent work on quantum dots32). In the following we apply spectral decomposition to study the spin-bath dynamics and resulting scaling of T2 with n for NV centres in diamond.

Experimental application of spectral decomposition. The NV centre is composed of a substitutional nitrogen atom and a vacancy on adjacent lattice sites (Fig. 1a). The negatively charged state, which is the focus of this work, is believed to gain an electron from nearby nitrogen donor impurities. The energy-level diagram of an NV centre is depicted in Fig. 1d. The NV centre has an electronic spin-triplet ground state with a zero-magnetic-eld splitting 2.87 GHz between the ms = 0 and ms = 1 spin states, quantized along the NV axis. A small external magnetic eld applied along this axis lis the degeneracy of the ms = 1 energy levels with a Zeeman shi 2.8 MHz G 1. Optical transitions between the electronic ground and excited states have a characteristic zero-phonon line at 637 nm, although 532-nm light is typically used to drive excitation to a phonon-sideband, and NV centres uoresce at room temperature over a broad range of wavelengths that is roughly peaked around 700 nm. Optical cycling transitions between the ground and excited states are primarily spin conserving. However, there exists a decay path that preferentially transfers the ms = 1 excited state population to the ms = 0 ground-state through a metastable singlet state, without emitting a photon in the uorescence band. It is this decay channel that allows the NV centres spin-state to be determined from the uorescence signal, and also leads to optical pumping into the ms = 0 ground-state.

Experimentally, we manipulate the |0 |1 spin manifold of the NV triplet electronic ground state using a static magnetic eld and resonant MW pulses, and use a 532-nm laser to initialize and provide optical readout of the NV spin states (Fig. 1d). More specically, we optically initialize the NV spins to ms = 0, apply CPMG pulse sequences (Fig. 1c) with varying numbers of -pulses n and varying free precession times , and then measure the NV spin state using optical readout to determine the remaining NV multi-spin coherence (for details see the Methods and refs 13,14,10).

The measured coherence is then used to extract the corresponding spin-bath spectral component Sn() as described above.

We applied the spectral decomposition technique to extract the spin-bath parameters and c as well as the NV multi-qubit coherence T2 and FOM for three diamond samples with diering

NV densities and concentrations of electronic and nuclear spin impurities (Table 1). The rst sample, referred to as the 12C sample, is an isotopically pure 12C diamond sample grown by chemical vapour deposition (Element 6), which we studied using an NV wide-eld microscope13. This sample has a very low concentration of 13C nuclear spin impurities (0.01%), a moderate concentration of N electronic spin impurities (~1 p.p.m.), and a moderate NV density (~1014(cm 3)). We then measured the Apollo sample, a chemical vapour deposition diamond sample with a thin N-doped layer with natural 13C concentration (1.1%), high N concentration (~100 p.p.m.), and large NV density (~1016(cm 3)), which we also studied using the NV wide-eld microscope. The last sample, referred to as the high-pressure high-temperature (HPHT) sample, is a type 1b HPHT diamond sample (Element 6) with natural 13C concentration, high N concentration (~50 p.p.m.) and low NV density (~1012(cm 3)), which we observed under a confocal microscope able to measure single NV centres.

The results for the 12C sample are summarized in Fig. 2. In Fig. 2a we plot examples of the measured NV multi-spin coherence decay Cn(t) as a function of pulse sequence duration t for CPMG pulse sequences with dierent numbers of -pulses n. The measured Cn(t)

are well described by a stretched exponential e t T

p

( / )

2 , which is consistent with an electronic spin bath described by a Lorentzian spectrum (Supplementary Methods). We analysed such NV decoherence data using the spectral decomposition technique outlined above to extract the best-t Lorentzian spin-bath spectrum (t to the average of all data points), yielding a coupling strength of = 30 10 kHz and correlation time c = 10 5 s, which agrees well with the range of values we nd for the Lorentzian spin-bath spectra Sn() t to each CPMG pulse sequence individually, 30 to 50 kHz and c5 to 15 s (Fig. 2b, see Methods). These values are in reasonable agreement with the expected N dominated bath values for and c for this samples estimated concentrations of 13C and N spins (Table 1), indicating that N electronic spin impurities are the dominant source of NV decoherence. In Fig. 2c we plot the NV multi-spin coherence lifetime T2, determined from the measured coherence decay Cn(t), as a function of the number n of CPMG -pulses.

We repeated the spectral decomposition analysis for the two other diamond samples, with the results for all three samples summarized in Fig. 3 and Table 1. We nd reasonable agreement between the measured and expected values for the NV/spin-bath

NATURE COMMUNICATIONS | 3:858 | DOI: 10.1038/ncomms1856 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1856

1

107

Coherence n (t)

n = 1 Pulsen = 4 Pulsesn = 12 Pulses n = 64 Pulses n = 128 Pulses

0.5

S n ( ) (a.u.)

105

103

0

100 101 102 103

104 105 106 107

Time (s)

(rad s1)

10

b

103

S n ( ) (a.u.)

5

(3050) kHz, c (515) s

= 30(8) kHz, c = 10(5) s

0 0 20 40 60 80 100 120

104

T 2(s)

105

0 0 0.2 0.4 0.6 0.8 1

6

10 100 101 102

(106 rad s1)

n (number of CPMG pulses)

Figure 3 | Comparative application of the spectral decomposition technique to the 12C (blue), Apollo (green) and HPHT (red) NV-diamond samples. (a) Spin-bath spectral functions data for all CPMG pulse sequences (dots), and associated Lorentzian ts using the average-t method (solid lines) and the individual-t method (colour bands); a.u., arbitrary unit. (b) Scaling of T2 with the number of CPMG pulses n: derived from NV spin coherence decay data Cn(t) (dots); t of dots to a powerlaw (solid lines); and synthesized from the average-t Lorentzian spin-bath spectrum (open squares). The tted scalings are: T2n0.43(6) (12C),

T2n0.65(5) (Apollo), T2n0.7(1) (HPHT).

2

T 2(ms)

1

DataFit to data Calculated

n (number of CPMG pulses)

Figure 2 | Application of the spectral decomposition technique to an isotopically pure 12C NV-diamond sample. (a) Examples of measured NV multi-spin coherence as a function of time Cn(t) for CPMG pulse sequences with different numbers of pulses n. (b) Derived values for the spin-bath spectral functions Sn() for all CPMG pulse sequences (blue dots) and average values at each frequency (red crosses). Each blue dot represents a pulse sequence with a specic duration t and number of pulses n, such that the probed frequency is = n/t. Also shown is the best-t Lorentzian for the mean spectral function Sn()n (solid black line), and a range of best-t Lorentzians for the individual spectral functions Sn() for each CPMG pulse sequence (light-blue band); a.u., arbitrary unit. (c) Scaling of T2 with the number n of CPMG pulses: derived from NV spin coherence decay data Cn(t) (dots); t of data to a power law T2 = 250(40) sn0.43(6) (solid line); and synthesized from the average-t Lorentzian spin-bath spectrum with = 30 kHz, c = 10 s (open squares), which yields a consistent t

T2 = 270(100) sn0.4(1).

N

(4)(4)

coupling strength in all three NV-diamond samples, with scaling approximately linearly with the N concentration. As mentioned before, the measured and expected values for the spin-bath correlation time c agree well for the 12C sample. However, we nd a striking discrepancy between the measured and expected values of c for the two samples with 1.1% 13C concentration (Apollo and

HPHT): both samples have measured spin-bath correlation times that are more than an order of magnitude longer than that given by the simple electronic spin-bath model, although the relative values of c for these two samples scale inversely with N concentration, as expected.

Discussion

We explain this suppression of spin-bath dynamics as a result of random, relative detuning of electronic spin energy levels due to interactions between proximal electronic (N) and nuclear (13C)

spin impurities (similar to processes identied by Bloembergen33 and Portis34 for other solid-state spin systems). The ensemble average eect of such random electronic-nuclear spin interactions is to induce an inhomogeneous broadening E of the resonant electronic spin transitions in the bath, which reduces the electronic spin ip-op rate R (~1/c) given by35,36

R E

p9 ,

2

where N is the dipolar interaction between N electronic spins. In this physical picture, E is proportional to the concentration of 13C impurities and to the N-13C hyperne interaction energy, whereas N is proportional to the N concentration. Given the magnetic moments and concentrations of the N and 13C spin impurities, and the large N13C hyperne interaction in diamond37,38, we estimate 10 MHz and N~1 MHz for the Apollo and HPHT samples.

These values imply an order of magnitude suppression of R compared with the bare electronic spin ip-op rate ignoring N13C interactions (Rbare~N), which is consistent with our experimental results for c~1/R (Table 1). The estimate of E is based on a Monte-

Carlo simulation of the relative hyperne interactions of 1,000 N spin pairs with a random distribution of 13C spin impurities at the natural abundance concentration of 1.1% (Supplementary Methods). For the

12C sample, the low concentration of 13C spins greatly reduces the magnitude of E, and hence the spin-bath suppression is negligible.

Importantly, we found large values for the multi-qubit FOM ~21014(ms cm 3) for both the 12C and Apollo samples. (The HPHT sample has a relatively low FOM, as expected given its low NV density.) The long coherence time obtained for the Apollo sample, with high NV density, natural 1.1% 13C abundance and large N concentration,

NATURE COMMUNICATIONS | 3:858 | DOI: 10.1038/ncomms1856 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1856

ARTICLE

is a result of the combined eect of suppression of electronic spin-bath dynamics and dynamical decoupling: the extended spin-bath correlation time enhances the eectiveness of dynamical decoupling, such that T2 is increased by a factor of ~10 despite the large concentration of N impurities (which would naively suggest almost no dynamic decoupling improvement of T2 for this sample30).

In summary, we applied the spectral decomposition technique to three NV-diamond samples with dierent composite-spin environments to characterize the spin-bath dynamics and determine the NV multi-spin-qubit FOM. For samples with a nite concentration of 13C nuclear spin impurities, this technique revealed an order-of-magnitude suppression of the N electronic spin-bath dynamics, which can be explained by random interactions between proximal electronic and 13C nuclear spin impurities. This spin-bath suppression enhances the efficacy of dynamical decoupling for samples with high N impurity concentration, enabling increased NV spin coherence times and thus realization of a FOM ~21014 (ms cm 3), which is within an order of magnitude of the state-of-the-art spin coherence FOM achieved in atomic systems15. Further optimization of the multi-spin-qubit FOM may be possible by engineering this spin-bath suppression, for example, with 13C concentration higher than the natural value. Such optimization, along with spectral decomposition studies of samples with varying 13C concentrations, at low temperatures and at high magnetic elds, will be pursued in future work. The present results, together with the possibility of single qubit addressability through AFM3941 or super-resolution optical techniques42 and intrinsic qubitqubit interactions1, pave the way for quantum information, sensing and metrology applications in a robust, multi-qubit solid-state architecture. Finally, the spectral decomposition technique presented here, based on well-known pulse sequences and a simple reconstruction algorithm, can be applied to other composite solid-state spin systems, such as quantum dots and phosphorous donors in silicon. Such measurements could provide a powerful approach for the study of many-body dynamics of complex spin environments, potentially exhibiting non-trivial eects related to the interplay between nuclear and electronic spin baths.

Methods

Single NV confocal microscope. Single NV measurements (HPHT sample) were performed using a custom-built confocal microscope. Optical excitation was provided by a 300 mW 532 nm diode pumped solid-state laser (Changchun New Industries Optoelectronics Tech, MLLIII532-300-1), focused onto the sample using a 100, NA = 1.3 oil-immersion objective (Nikon CFI Plan Fluor 100 oil). NV uorescence was collected through the same optics, and separated from the excitation beam using a dichroic lter (Semrock LM01-552-25). The light was additionally ltered (Semrock LP02-633RS-25) and focused onto a single-photon counting module (Perkin-Elmer, SPCM-ARQH-12). The excitation laser was pulsed by focusing it through an acousto-optic modulator (Isomet 1205C-2).

Microwaves were delivered to the sample through a 20-m thick wire soldered across it. The wire was driven by an amplied (Mini-circuits ZHL-16W-43-S + ) MW synthesizer (Agilent E4428C). Phase modulation of the MW pulses was achieved using an in-phase/quadrature (IQ) mixer (Marki IQ-1545). MW and optical pulses were controlled using a computer-based digital delay generator (SpinCore PulseBlaster ESR400). A static magnetic eld of 50 G was applied to the sample to li the degeneracy of the ms = 1 levels.

NV wide-eld uorescence microscope. NV multi-spin measurements (for samples 12C and Apollo) were performed using a custom-built wide-eld uorescence microscope (Fig. 4a), which has previously been used for two-dimensional magnetic eld imaging with thin-layer NV-diamond samples13. Optical excitation was provided by a 3-W 532-nm laser (Laser Quantum Opus), which illuminateda large region of interest on the diamond surface through a 20 NA = 0.75 objective (Nikon). An acousto-optic modulator (Isomet M1133-aQ80 L-H) acted asan optical switch, to pulse the laser with precise timing in order to prepare and detect the NV spin states. A loop antenna was positioned near the diamond surface and connected to the amplied (Mini-circuits ZHL-16W-43-S + ) output of a microwave signal generator (Agilent E8257D) to generate a homogeneous B1

eld over the region of interest. Fast-switching of the MW eld (ZASWA-2-50DR) enabled coherent manipulation of the NV spin states necessary for coherence decay measurements using CPMG pulse sequences, in order to perform the spectral

decomposition technique. Optical and MW pulse timings were controlled through a computer-based digital delay generator (SpinCore PulseBlaster PRO ESR500).

NV uorescence was collected by the objective, ltered, and imaged onto a cooled charge-coupled device camera (Starlight Xpress SXV-H9). As the duration of a single measurement is shorter than the minimum exposure time of the camera, the measurement was repeated for several thousand averages within a single exposure and synchronized to an optical chopper placed before the camera in orderto block uorescence from the optical preparation pulse (Fig. 4b). Repeating the measurement without the MW control pulses provided a reference for long-term dris in the uorescence intensity.

In the measurements described here, the loop of wire delivered 3.07 GHz MW pulses to the sample, resonant with the NV |0 |1 spin transition for the applied static magnetic eld 70 G, to manipulate the NV spin coherence and implement CPMG spin-control pulse sequences.

Spectral decomposition deconvolution procedure. The coherence of a two-level quantum system can be related to the magnitude of the o-diagonal elementsof the systems density matrix. Specically, we deal here with NV electronic spin qubits in a nite external magnetic eld, which can be treated as eective two-level spin systems with quantization (z) axis aligned with the NV axis. When the NV spins are placed into a coherent superposition of spin eigenstates, for example, aligned with the x axis of the Bloch sphere, the measurable spin coherence is given by C t Tr t Sx

( ) = [ ( ) ]

r .

The lter function for the n-pulse CPMG control sequence FCPMG() covers a narrow frequency region (given by /t, where t is total length of the sequence), which is centred at 0 = n/t, and is given by43:

F t t

t n t n

CCD

Optical chopper

Red filters

Dichroic mirror

AOM

Objective

Thin NV layer on High-purity diamond

Laser

Optical chopper

Optical excitation

Fluorescence blocked

Pump Readout

/2 /2

Microwave

CPMG

Figure 4 | Details of ensemble measurement setup and pulse sequences. (a) Diagram of the NV wide-eld uorescence microscope. The laser produces 532-nm light, which is switched by the acousto-optic modulator (AOM) and directed through a dichroic mirror and an objective ontothe diamond sample. The uorescence from the sample passes through the dichroic and, following an optical chopper and appropriate lters, is collected by a charge-coupled device (CCD). (b) Pulse sequence used to perform spectral decomposition measurements with the NV wide-eld microscope. Optical pulses are used to rst initialize the NV and thento readout its spin state. The optical chopper is synced such that the initialization pulse is blocked from the CCD, while the readout pulse is recorded. The microwave pulses are applied between the initialization and readout optical pulses.

nCPMG( ) = 8 2

4

2

w w

w

w

2

sin

cos

4

2

sin

.

(5)(5)

NATURE COMMUNICATIONS | 3:858 | DOI: 10.1038/ncomms1856 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1856

21. Almog, I. et al. Direct measurement of the system environment coupling as a tool for understanding decoherence and dynamical decoupling. J. Phys. B: Atomic, Mol Opt. Phys. 44, 154006 (2011).

22. Bylander, J. et al. Noise spectroscopy through dynamical decoupling with a superconducting ux qubit. Nat. Phys. 7, 565570 (2011).

23. lvarez, G. A. & Suter, D. Measuring the spectrum of colored noise by dynamical decoupling. Phys. Rev. Lett. 107, 230501 (2011).

24. Du, J. et al. Preserving electron spin coherence in solids by optimal dynamical decoupling. Nature 461, 12651268 (2009).

25. Klauder, J. R. & Anderson, P. W. Spectral diusion decay in spin resonance experiments. Phys. Rev. 125, 912932 (1962).

26. Gordon, G., Erez, N. & Kurizki, G. Universal dynamical decoherence control of noisy single- and multi-qubit systems. J. Phys. B 40, S75 (2007).

27. Meiboom, S. & Gill, D. Modied spin-echo method for measuring nuclear relaxation times. Rev. Sci. Instr. 29, 688691 (1958).28. Hahn, E. L. Spin echoes. Phys. Rev. 80, 580594 (1950).29. Slichter, C. P. Principles of Magnetic Resonance (Springer-Verlag, 1990).30. de Sousa, R. Electron spin as a spectrometer of nuclear spin noise and other uctuations. In: Electron Spin Resonance and Related Phenomena in Low-Dimensional Structures, of Topics in Applied Physics (ed. Fanciulli, M.) 115, 183220 (2009).

31. de Lange, G., Rist, D., Dobrovitski, V. V. & Hanson, R. Single-spin magneto-metry with multipulse sensing sequences. Phys. Rev. Lett. 106, 080802 (2011).

32. Medford, J. et al. Scaling of dynamical decoupling for spin qubits. Phys. Rev. Lett. 108, 086802 (2012).

33. Bloembergen, N. On the interaction of nuclear spins in a crystalline lattice. Physica 15, 386426 (1949).

34. Portis, A. M. Spectral diusion in magnetic resonance. Phys. Rev. 104, 584588 (1956).

35. de Sousa, R. & Das Sarma, S. Electron spin coherence in semiconductors: considerations for a spin-based solid-state quantum computer architecture. Phys. Rev. B 67, 033301 (2003).

36. Witzel, W. M., Carroll, M. S., Morello, A., Cywinski, L. L. & Das Sarma, S. Electron spin decoherence in isotope-enriched silicon. Phys. Rev. Lett. 105, 187602 (2010).

37. Cox, A., Newton, M. E. & Baker, J. M. 13C, 14N and 15N endor measurements on the single substitutional nitrogen centre (p1) in diamond. J. Phys.: Condens. Matter 6, 551 (1994).

38. Barklie, R. C. & Guven, J. 13C hyperne structure and relaxation times of the p1 center in diamond. J. Phys. C 14, 3621 (1981).

39. Grinolds, M. S. et al. Quantum control of proximal spins using nanoscale magnetic resonance imaging. Nat. Phys. 7, 687692 (2011).

40. Arcizet, O. et al. A single nitrogen-vacancy defect coupled to a nanomechanical oscillator. Nat. Phys. 7, 879883 (2011).

41. Maletinsky, P. et al. A robust scanning diamond sensor for nanoscale imaging with single nitrogen-vacancy centres. Nat. Nano. doi:10.1038/nnano.2012.50 (2012).

42. Maurer, P. C. et al. Far-eld optical imaging and manipulation of individual spins with nanoscale resolution. Nat. Phys. 6, 912918 (2010).

43. Cywiski, L. L., Lutchyn, R. M., Nave, C. P. & Das Sarma, S. How to enhance dephasing time in superconducting qubits. Phys. Rev. B 77, 174509 (2008).

Acknowledgements

We gratefully acknowledge fruitful discussions with Patrick Maletinsky and Shimon Kolkowitz; and assistance with samples by Daniel Twitchen and Matthew Markham (Element 6) and Patrick Doering and Robert Linares (Apollo). This work was supported by NIST, NSF, US Army Research Office and DARPA (QuEST and QuASAR programs). PC was partially funded by NSF under grant no. DMG-1005926.

Author contributions

All authors contributed to all aspects of this work.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/ naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/ reprintsandpermissions/

How to cite this article: Bar-Gill, N. et al. Suppression of spin-bath dynamics for improved coherence of multi-spin-qubit systems. Nat. Commun. 3:858 doi: 10.1038/ncomms1856 (2012).

We note that the narrow-band feature of the CPMG lter essentially denes the eectiveness of this sequence for dynamical decoupling.

The lter function F t

nCPMG( )

w associated with the CPMG pulse sequence deviates from a -function by the nite width of the main spectral peak and by the presence of higher harmonics (Supplementary Fig. 1). In addition, the central frequency of the lter function changes with the duration of the experiment and with the number of pulses. Therefore, reconstructing the spin-bath spectral function S() from the decoherence data requires a solution of the Fredholm type equation to extract S() equation (1):

c p w w

w w

( ) = 1 ( ) ( ) .

0 2

t d S Ft

(6)(6)

This extraction is accomplished by assuming that the spectral function decays to zero at high frequency and by noting that the lter function includes high harmonics but negligible low frequency artefacts. Thus, the procedure starts at the high frequency values of C(), and works back to lower frequencies by subtracting the eect of higher harmonics using the reconstructed high frequency values of S() and the analytical form of the lter function.

Owing to the high frequency components of the lter function and the monotonically decreasing nature of the spectral function, a naive reconstruction of the spectral function assuming a -function form of the lter function produces results that are biased to lower values by ~15%. This bias is removed using the reconstruction algorithm presented above.

References

1. Jiang, L. et al. Repetitive readout of a single electronic spin via quantum logic with nuclear spin ancillae. Science 326, 267272 (2009).

2. Neumann, P. et al. Single-shot readout of a single nuclear spin. Science 329, 542544 (2010).

3. McCamey, D. R., Van Tol, J., Morley, G. W. & Boehme, C. Electronic spin storage in an electrically readable nuclear spin memory with a lifetime >100 s. Science 330, 16521656 (2010).

4. Zurek, W. H. Decoherence, einselection, and the quantum origins of the classical. Rev. Mod. Phys. 75, 715775 (2003).

5. Taylor, J. M. et al. High-sensitivity diamond magnetometer with nanoscale resolution. Nat. Phys. 4, 810816 (2008).

6. Maze, J. R. et al. Nanoscale magnetic sensing with an individual electronic spin in diamond. Nature 455, 644647 (2008).

7. Balasubramanian, G. et al. Nanoscale imaging magnetometry with diamond spins under ambient conditions. Nature 455, 648651 (2008).

8. Goldstein, G. et al. Environment-assisted precision measurement. Phys. Rev. Lett. 106, 140502 (2011).

9. Acosta, V. M. et al. Diamonds with a high density of nitrogen-vacancy centers for magnetometry applications. Phys. Rev. B 80, 115202 (2009).

10. Childress, L. et al. Coherent dynamics of coupled electron and nuclear spin qubits in diamond. Science 314, 281285 (2006).

11. Tyryshkin, A. M., Lyon, S. A., Astashkin, A. V. & Raitsimring, A. M. Electron spin relaxation times of phosphorus donors in silicon. Phys. Rev. B 68, 193207 (2003).

12. Hanson, R., Kouwenhoven, L. P., Petta, J. R., Tarucha, S. & Vandersypen, L. M. K. Spins in few-electron quantum dots. Rev. Mod. Phys. 79, 12171265 (2007).13. Pham, L. M. et al. Magnetic eld imaging with nitrogen-vacancy ensembles. New J. Phys. 13, 045021 (2011).

14. Stanwix, P. L. et al. Coherence of nitrogen-vacancy electronic spin ensembles in diamond. Phys. Rev. B 82, 201201 (2010).

15. Allred, J. C., Lyman, R. N., Kornack, T. W. & Romalis, M. V. High-sensitivity atomic magnetometer unaected by spin-exchange relaxation. Phys. Rev. Lett. 89, 130801 (2002).

16. Dobrovitski, V. V., Feiguin, A. E., Awschalom, D. D. & Hanson, R. Decoherence dynamics of a single spin versus spin ensemble. Phys. Rev. B 77, 245212 (2008).

17. Barnes, E., Cywinski, L. L. & Das Sarma, S. Master equation approach to the central spin decoherence problem: Uniform coupling model and role of projection operators. Phys. Rev. B 84, 155315 (2011).

18. de Lange, G., Wang, Z. H., Rist, D., Dobrovitski, V. V. & Hanson, R. Universal dynamical decoupling of a single solid-state spin from a spin bath. Science 330, 6063 (2010).

19. Naydenov, B. et al. Dynamical decoupling of a single-electron spin at room temperature. Phys. Rev. B 83, 081201 (2011).

20. Ryan, C. A., Hodges, J. S. & Cory, D. G. Robust decoupling techniques to extend quantum coherence in diamond. Phys. Rev. Lett. 105, 200402 (2010).

NATURE COMMUNICATIONS | 3:858 | DOI: 10.1038/ncomms1856 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.