arxiv: v1 [quant-ph] 2 Apr 2018

Transcription

1 Longitudinal relaxation of a nitrogen-vacancy center in a spin bath by generalized cluster-correlation expansion method Zhi-Sheng Yang, Ming-Jie Tao, Mei Zhang, Qing Ai, and Fu-Guo Deng Department of Physics, Applied Optics Beijing Area Major Laboratory, Beijing Normal University, Beijing 875, China Wen Yang Beijing Computational Science Research Center, Beijing 94, China (Dated: April 4, 28) arxiv:84.8v [quant-ph] 2 Apr 28 We theoretically study the longitudinal relaxation of a nitrogen-vacancy (NV) center surrounded by a 3 C nuclear spin bath in diamond. By incorporating electron spin in the cluster, we generalize the cluster-correlation expansion (CCE) to theoretically simulate the population dynamics of electron spin of NV center. By means of the generalized CCE, we numerically demonstrate the decay process of electronic state induced by cross relaxation at the ambient temperature. It is shown that the CCE method is not only capable of describing pure-dephasing effect at large-detuning regime, but it can also simulate the quantum dynamics of populations in the nearly-resonant regime. I. INTRODUCTION Nitrogen-vacancy(NV) center is a pair of point defects in the diamond lattices,2. It is comprised of a substitutional nitrogen atom and a vacancy at the adjacent lattice. Because of possessing S = electron spin, the ground-state manifold of the negatively-charged NV center is triplet states and it can be optically initialized and measured. For a sample of high purity, the NV center possesses a long coherence time and can be further prolonged to approach one second by dynamical decoupling 3,4. It can be coherently controlled by microwave fields and is highly sensitive to external electric and magnetic fields. These superior characteristics make the NV center a promising candidate for quantum information processing 5 7 and quantum sensing of electric and magnetic fields 8, temperature, stress 2, biological structures 3 and chemical reactions 4. Such appealing applications require full knowledge of its decoherence behavior, including both transverse and longitudinal relaxation 5 7, over a broad range of physical parameters, such as temperature, magnetic field, type and concentration of impurities. Among quantum many-body theories for many-spin bath dynamics, e.g. the density matrix cluster expansion 8 2, the pair correlation approximation 2 23, and the linked-cluster expansion 24, the cluster-correlation expansion (CCE) 25,26 developed by one of the authors WY has been successfully applied to describe the pure-dephasing process of a central spin coupling to a many-spin bath in the largely-detuned regime 27. TheanomalousdecoherenceeffectofNVcenter inaquantumspinbathhasbeentheoreticallypredicted 28 and experimentally verified 29. Interestingly, the decoherence of NV center, which was previously considered as an obstacle to realize quantum computation 3,3, has been explored to detect a remote nuclear spin dimer 32. Meanwhile, the experimental investigation on the longitudinal relaxation of NV center has been delicately performed for a broad range of temperature, magnetic field, concentration of impurities, ever since The physical mechanism of longitudinal relaxation above K is attributed to Orbach and Raman phonon process 33,34. When the temperature reaches as low as a few tens of kelvins, T clearly demonstrates dependence on the concentration of magnetic impurities 33. In addition, when varying the magnetic field, T abruptly drops around some energy-level anti-crossing points 33. These observations motivate us to develop a theory to describe the longitudinal relaxation of NV center in a quantum many-spin bath. In this paper, by including the electron spin in the cluster expansion, we generalize the CCE to simulate the quantum dynamics of NV center exchanging population with a nuclear spin bath in the nearlyresonant regime. When astatic magneticfield isappliedtomakethe NV center close to the level anticrossing point, the hyperfine interaction of the NV center with the 3 C nuclear spins induces the NV electron spin relaxation. The numerical simulation with the CCE approach yields convergent result and thus is capable of fully describing the NV electron spin relaxation induced by the 3 C spin bath. This paper is organized as follows: In Sec. II, we present our system model. In Sec. III, we show the CCE approach after generalization. In Sec. IV, we simulate the longitudinal relaxation of NV center surrounded by 3 C nuclear bath with the generalized CCE. A detailed discussion and a summary are enclosed in Sec. V. II. MODEL As shown in Fig., a negatively-charged NV center is coupled to the 3 C nuclear spins randomly located at the diamond lattices, with the natural abundance of 3 C being.% 37. Here, since we focus on the longitudinal relaxation induced by 3 C nuclear spins far away from the NV center, in our numerical simulation we use the

2 2 (a) (b) ",* $ N B V! = " $! = " #$ % & '#( ",) $ " #,*$ " #,) $ FIG. : (color online) Schematic illustration of an NV center surrounded by a spin bath. (a) An NV center in the crystal structure of a diamond, where a static magnetic field is applied along the principal axis of NV center. (b) The hyperfine structure of an NV center and one 3 C nuclear spin in the bath. randomly-generated bath in which there is no 3 C nuclear spin in the vicinity of NV center 37. The NV center consists of a substitutional nitrogen atom adjacent to a vacancy. The ground state of NV center is a triplet state with m s = and m s = ± denoted by and ± respectively. When there is no magnetic field applied to the NV center, ± are degenerate and they are separated from by zero-field splitting D = 2.87GHz,2. In an external magnetic field B z along the NV axis, the degeneracy between m s = ± levels are lifted and thus the full Hamiltonian of the total system describing the NV center and spin bath is H = H NV +H bath +H int, () where the Hamiltonian of the NV center has the form H NV = DS 2 z γ eb z S z (2) with γ e =.76 rad s T being the electronic gyromagnetic ratio 37 and S z =. The Hamiltonian of the bath reads H bath = [ ] D ij I i I j 3(I i r ij )(r ij I j ) r 2 i<j ij γ c B z Ii z, (3) i where the gyromagnetic ratio of 3 C nuclear spin is 37 γ c = rad s T, ( ) D ij = µ γ 2 c 4πr 3 ij 3r ijr ij r 2 ij (4) is the magnetic dipole-dipole interaction between the ith andjth 3 Cnuclearspinswith r ij beingthe displacement vector from ith to jth spins and µ being the vacuum permeability. The electron spin and nuclear spins are coupled by hyperfine interaction. Because we are interested in a large numberof 3 CnuclearspinsfarawayfromtheNVcenter, the hyperfine interaction is mainly described by magnetic dipole-dipole interaction and thus the interaction Hamiltonian is H int = S A i I i, (5) where A i = µ γ c γ e 4πr 3 i i ( 3r ) ir i ri 2 (6) is the hyperfine coupling tensor, r i is the position vector of ith nuclear spin, and the position of NV center is chosen as the origin. III. GENERALIZED CLUSTER-CORRELATION EXPANSION In this section, by taking electron spin into consideration in the cluster, we generalize the CCE method 25,26 to describe the electron spin relaxation in a spin bath. Traditional CCE method has been successfully applied to calculating the pure dephasing of an NV center in the large-detuning regime. However, the central spin flip must be considered when the energy relaxation of the NV center is involved in the nearly-resonant regime. Assuming that the bath consists of spin i only, cf. Fig. 2(a), the probability P(t) of NV center in the state

3 3 (a) 2-CCE (b) 3-CCE (c) 4-CCE NV NV NV FIG. 2: (color online) Schematic of generalized CCE. The bath of 3 C nuclear spins and the NV center forms an open quantum system. The 3 C atoms with /2 nuclear spin denoted by the small red arrows are randomly located on the diamond lattice. (a) 2-CCE: A cluster containing one bath spin and the NV center. (b) 3-CCE: A cluster containing two bath spins and the NV center. (c) 4-CCE: A cluster containing three bath spins and the NV center. The spins outside the cluster are assumed to be frozen when calculating the cluster contribution. is explicitly calculated as P {i} = P {i} Tr( ρ()e ih {i}t e ih {i}t ) Tr(ρ()e ihnvt e ihnvt ), (7) where ρ() is the initial state of the total system including the NV center and nuclear spin bath, H {i} = H NV γ c B z I z i +S A i I i (8) is obtained from Eq. () by droppingall terms other than spin i. Assuming that the bath consists of spin i and spin j, cf. Fig. 2(b), P(t) reads P {i,j} Tr( ρ()e ih {i,j}t e ih {i,j}t ) Tr(ρ()e ihnvt e ihnvt ), (9) where H {i,j} = H NV γ c B z +D ij [ α=i,j I z α +S α=i,j I i I j 3(I i r ij )(r ij I j ) r 2 ij A α I α ] () is calculated from Eq. () by dropping all terms other than spin i and spin j, and the spin-pair correlation is P {i,j} P {i,j} P {i} P{j}. () Assuming that the bath consists of three spins i and j and k, cf. Fig. 2(c), P(t) in this case becomes P {i,j,k} Tr( ρ()e ih {i,j,k}t e ) ih {i,j,k}t Tr(ρ()e ihnvt e ihnvt, (2) ) where H {i,j,k} = H NV γ c B z + α=i,j,k β>α [ D αβ α=i,j,k I z α +S α=i,j,k A α I α I α I β 3(I α r αβ )(r αβ I β ) r 2 αβ ] (3) is calculated from Eq. () by dropping all terms other than the three spins i and j and k, and the three-spin correlation is P {i,j,k} P {i,j,k} P {i} P{j} P{k} P{i,j} P{j,k} P{k,i}. (4) For the bath with arbitrary number of spins, P(t) is generalized as P c Tr( ρ()e ihct e ihct) Tr(ρ()e ihnvt e ihnvt ), (5) where H c is obtained from Eq. () by dropping all terms other than spins in the cluster c, and the spin-cluster correlation is P c P c c c P c. (6) Furthermore, in Eqs. (7,9,2,5), all of the denominators are equal to unity because the system is initially prepared at, which is also the eigenstate of H NV. Generally speaking, it is impossible to exactly calculate P(t) for a large number of spins as the dimension of Hilbert space scales exponentially with the number of spins. The M-CCE method approximates P(t) as P (M) = P c, (7) c M

4 4 where c is the number of spins in the cluster c. For example, the first-order truncation of P(t) yields P () = i P {i}. (8) And the second-order truncation of P(t) reads P (2) = P {i} P {i,j}. (9) i i,j IV. DECAY OF NV CENTER INDUCED BY NUCLEAR SPIN BATH In this section, we discuss the longitudinal relaxation of electron spin of NV center induced by coupling to the nuclear spin bath. When the energy gap between and approaches the energy gap of 3 C nuclear spins, e.g. by tuning the magnetic field, the electron spin exchanges polarization with nuclear spins, cf. Fig. (c), and the longitudinal relaxation of electronic spin occurs. For simplicity, we assume that the decay of electron spin is mainly due to hyperfine interaction with nuclear spins, and the effect of diamond lattice phonons on the relaxation is beyond the scope of the present investigation. The state of electronic spin can be optically initialized into, corresponding to the density matrix ρ NV =. (2) In principle, the 3 C nuclear spin bath should be in the thermal equilibrium state since the typical experimental temperature is much higher than the nuclear spin Zeeman splitting, even in a strong magnetic field, e.g. several T. Therefore, without loss of generality, we take the density matrix of the 3 C bath as ρ B = 2 N N i= ( i + i ), (2) where i ( i ) is the spin-up (spin-down) state of the ith 3 C nuclearspin, with the quantizationaxis alongthe N-V symmetry axis. As a result, the density matrix of total system at the initial time is ρ() = ρ NV ρ B. (22) The evolution of the coupled system is given by ρ(t) = e iht ρ()e iht. (23) By partially tracing over the degrees of freedom of the bath, we could obtain the survival probability of the initial state of the NV electron spin as P(t) = Tr B e iht ρ()e iht. (24) For a small spin bath, we can exactly calculate the longitudinal relaxation of electron spin via Eq. (24). However, with increasing size of the nuclear spin bath, this approach quickly becomes unfeasible, because the dimension of the Hilbert space grows exponentially with the numberofnuclei 38,39. Inthiscase,anapproximatemanybody theory with high performance has to be considered. In previous works, the CCE theory has been successfully appliedtodescribethepuredephasingofacentralspinin aspinbath 25,26,28,32. Here, wegeneralizetheccetheory to deal with the longitudinal relaxation of the electron spin in a spin bath. In Fig. 3(a), we numerically simulate the survival probability of the initial state of electron spin surrounded by N = 5 nuclear spins under the magneticfieldintensityb z = GalongtheNVaxis. When the order of CCE method is increased, the results quickly converge, e.g., the results from 4th-order CCE and 5th-order CCE are indistinguishable. This suggests that 4th-order CCE already offers a reliable solution to the quantum dynamics of central spin under the influence of a many-spin bath. Figure 3(a) shows the population of electron spin experiences a damped oscillation with a small amplitude because there is still a large energy gap between the electron spin and nuclear spin as compared to their hyperfine interaction. When we reduce the energy gap of electron spin by tuning the magnetic field, cf. Fig. 3(b), both the amplitude and period of damped oscillation becomes larger as the electron and nuclear spins are getting close to resonance. If we further tune the electronic level to the nearly-resonant case, as shown in Fig. 3(c), the population of electronic initial state decays to zero around µs without collapse and revival. When the electron spin resonantly interacts with a nuclear spin, its population will demonstrate Rabi oscillation with full amplitude and frequency A i. Since all nuclear spins interact with the electron spin with different hyperfine interactions, the electronic population will irreversibly decay and no revival will occur. Furthermore, we explore the effect of the bath s size on the relaxation. In Fig. 3(d), the population dynamics of the electron spin is investigated for 3 C baths of different sizes with a large energy gap. When the size of the bath is increased, the difference before and after the change becomes more and more negligible. As further enlarging the size will not significantly decrease the difference but result in much more computation time, the 4-CCE theory with a bath of N = 5 nuclear spins already yields a reliable result, as confirmed in Fig. 3(e) and (f). V. DISCUSSION AND SUMMARY In this paper, by embodying the electron spin in the cluster, we generalize the CCE approach to deal with the longitudinal relaxation of a central electron spin due to its couplings to a nuclear spin bath in diamond. In the largely-detuned case, because the electron spin can not effectively exchange energy with the nuclear spins, the couplings to the nuclear spin bath manifest as the puredephasing phenomenon of the NV center 37. Since the pure-dephasing effect is sensitive to the local magnetic

5 5 Qubit Population P (a).4.2 (b).4.2 (c).4.2 2CCE 3CCE 4CCE 5CCE Qubit Population P (d).4.2 (e).4.2 (f).4.2 N=2 N=5 N= FIG. 3: (color online) The longitudinal relaxation process of the NV center under different magnetic fields. Survival probability of initial state of electron spin by different orders of CCE with N = 5 and (a) B z = 25.G, (b) B z = 24.99G, (c) B z = 24.97G. Survival probability with different bath sizes and (d) B z = 25.G, (e) B z = 24.99G, (f) B z = 24.97G. field at the NV center, it can be fully utilized to detect magnetic source with internal structure 32. When getting close to the degenerate point, we show diverse decay phenomena of initial state s population by tuning the magnetic field. When the energy gap of electronic levels are small but not too small, only a few nuclear spins can exchange population with the electron spin. It results in a damped oscillation with a very small decay rate. If we further reduce the energy gap, a more remarkable damped oscillation be observed as more nuclear spins can be probed by the electron spin. When the electron spin is finally tuned to be in resonance with all nuclear spins, the population of initial state will monotonically decay and no revival will take place. When we compare the steady-state populations of the above three cases, an interesting phenomenon is present. The equilibrium-state population of the excited state gradually decreases to zero along with the shrink of energy gap. This is in sharp contrast with the Boltzmann distribution, which predicts a rising in excited state population when decreasing the energy gap. We consider the feasibility with the current accessible parameters in the experiments. The state of NV center can be initialized by applying 532 nm laser. In order to bring the electronic spin close to the degenerate point, a magnetic field with strength B z 25 G should be applied to NV center, which is available at current laboratory conditions 4. After the optical pumping, the NV center will be left to evolve under the influence of nuclear spin bath for an interval and then readout by detecting the fluorescence after spin-dependent optical excitation 4. In this paper, because of the generalization, the CCE approach is capable of describing the quantum evolution of the whole density matrix of the system including both the off-diagonal and diagonal terms. Recently, quantum coherent energy transfer has attracted broad interest from different disciplines. Due to couplings to environment, optimal and efficient excited-state energy transfer is assisted by the noise. Therein, the environment is made up of protein scaffold with charge, which manifests forced vibration when the neighbouring pigment is in the excited state. Since the spin degree of freedom of electron is of different nature as compared to the harmonic oscillator, the excited-state energy transfer in a spin bath would be substantially different from that in a bath of harmonic oscillators. In this case, a different approach other than the traditional master equation 42 is

6 6 desired. In our forthcoming paper, we will subtly make the comparison by using two numerically exact approach, i.e. the hierarchical equation of motion 43 and CCE. ACKNOWLEDGMENTS We thank L. P. Yang for helpful discussions. F. G. Deng was supported by the National Natural Science Foundation of China under Grant No and the Fundamental Research Funds for the Central Universities under Grant No. 25KJJCA. W. Yang was supported by the National Natural Science Foundation of China under Grant No and No , and the MOST under Grant No. 24CB8487. M. Zhang was supported by the National Natural Science Foundation of China under Grant No Q. Ai was supported by the National Natural Science Foundation of China under Grant No. 557, and the Open Research Fund Program of the State Key Laboratory of Low-Dimensional Quantum Physics, Tsinghua University under Grant No. KF252. Electronic address: M. W. Doherty, N. B. Manson, P. Delaney, F. Jelezko, J. Wrachtrup, and L. C. L. Hollenberg, The nitrogen-vacancy colour centre in diamond, Phys. Rep. 528, (23). 2 V. V. Dobrovitski, G. D. Fuchs, A. L. Falk, C. Santori, and D. D. Awschalom, Quantum control over single spins in diamond, Annu. Rev. Condens. Matter Phys. 4, 23 (23). 3 G. de Lange, Z. H. Wang, D. Ristè, V. V. Dobrovitski, and R. Hanson, Universal dynamical decoupling of a single solid-state spin from a spin bath, Science 33, 6 (2). 4 N. Bar-Gill, L. M. Pham, A. Jarmola, D. Budker, and R. L. Walsworth, Solid-state electronic spin coherence time approaching one second, Nat. Commun. 4, 743 (23). 5 L. Childress, J. M. Taylor, A. S. Sørensen, and M. D. Lukin, Fault-tolerant quantum communication based on solid-state photon emitters, Phys. Rev. Lett. 96, 754 (26). 6 M. V. G. Dutt, L. Childress, L. Jiang, E. Togan, J. Maze, F. Jelezko, A. S. Zibrov, P. R. Hemmer, and M. D. Lukin, Quantum register based on individual electronic and nuclear spin qubits in diamond, Science 36, 32 (27). 7 X. K. Song, H, Zhang, Q. Ai, J. Qiu, and F.-G. Deng, Shortcuts to adiabatic holonomic quantum computation in decoherence-free subspace with transitionless quantum driving algorithm, New J. Phys. 8, 23 (26). 8 F. Dolde, H. Fedder, M. W. Doherty, T. Nobauer, F. Rempp, G. Balasubramanian, T. Wolf, F. Reinhard, L. C. L. Hollenberg, F. Jelezko, and J. Wrachtrup, Electricfield sensing using single diamond spins, Nat. Phys. 7, 459 (2). 9 J. R. Maze, P. L. Stanwix, J. S. Hodges, S. Hong, J. M. Taylor, P. Cappellaro, L. Jiang, M. V. G. Dutt, E. Togan, A. S. Zibrov, A. Yacoby, R. L. Walsworth, and M. D. Lukin,, Nanoscale magnetic sensing with an individual electronic spin in diamond, Nature 45, 644 (28). J. P. Tetienne, T. Hingant, L. Rondin, A. Cavaillès, L. Mayer, G. Dantelle, T. Gacoin, J. Wrachtrup, J. F. Roch, and V. Jacques, Spin relaxometry of single nitrogenvacancy defects in diamond nanocrystals for magnetic noise sensing, Phys. Rev. B 87, (23). D. M. Toyli, C. F. de las Casas, D. J. Christle, V. V. Dobrovitski, and D. D. Awschalom, Fluorescence thermometry enhanced by the quantum coherence of single spins in diamond, Proc. Natl. Acad. Sci. U.S.A, 847 (23). 2 F. Grazioso, B. R. Patton, P. Delaney, M. L. Markham, D. J. Twitchen, and J. M. Smith, Measurement of the full stress tensor in a crystal using photoluminescence from point defects: The example of nitrogen vacancy centers in diamond, App. Phys. Lett. 3, 95 (23). 3 F. Z. Shi, Q. Zhang, P. F. Wang, H. B. Sun, J. R. Wang, X. Rong, M. Chen, C. Y. Ju, F. Reinhard, H. W. Chen, J. Wrachtrup, J. F. Wang, and J. F. Du, Single-protein spin resonance spectroscopy under ambient conditions, Science 347, 35 (25). 4 H. B. Liu, M. B. Plenio, and J. M. Cai, Scheme for detection of single-molecule radical pair reaction using spin in diamond, Phys. Rev. Lett. 8, 242 (27). 5 L. T. Hall, P. Kehayias, D. A. Simpson, A. Jarmola, A. Stacey, D. Budker, and L. C. L. Hollenberg, Detection of nanoscale electron spin resonance spectra demonstrated using nitrogen-vacancy centre probes in diamond, Nat. Commun. 7, 2 (26). 6 H. J. Wang, C. S. Shin, S. J. Seltzer, C. E. Avalos, A. Pines, and V. S. Bajaj, Optically detected cross-relaxation spectroscopy of electron spins in diamond, Nat. Commun. 5, 435 (25). 7 P. Wang and W. Yang, Theory of nuclear spin dephasing and relaxation by optically illuminated nitrogen-vacancy center, New J. Phys. 7, 34 (25). 8 W. M. Witzel and S. Das Sarma, Quantum theory for electron spin decoherence induced by nuclear spin dynamics in semiconductor quantum computer architectures: Spectral diffusion of localized electron spins in the nuclear solidstate environment, Phys. Rev. B 74, (26). 9 W. M. Witzel and S. Das Sarma, Multiple-pulse coherence enhancement of solid state spin qubits, Phys. Rev. Lett. 98, 776 (27). 2 W. M. Witzel and S. Das Sarma, Concatenated dynamical decoupling in a solid-state spin bath, Phys. Rev. B 76, 2433(R) (27). 2 W. Yao, R. B. Liu, and L. J. Sham, Theory of electron spin decoherence by interacting nuclear spins in a quantum dot, Phys. Rev. B 74, 953 (26). 22 W. Yao, R. B. Liu, and L. J. Sham, Restoring coherence lost to a slow interacting mesoscopic spin bath, Phys. Rev. Lett. 98, 7762 (27). 23 R.B. Liu, W. Yao, andl. J. Sham, Control of electron spin decoherence caused by electroncnuclear spin dynamics in a quantum dot, New J. Phys. 9, 226 (27). 24 S. K. Saikin, W. Yao, and L. J. Sham, Single-electron spin decoherence by nuclear spin bath: Linked-cluster expansion approach, Phys. Rev. B 75, 2534 (27). 25 W. Yang and R. B. Liu, Quantum many-body theory of qubit decoherence in a finite-size spin bath, Phys. Rev. B

7 7 78, 8535 (28). 26 W. Yang and R. B. Liu, Quantum many-body theory of qubit decoherence in a finite-size spin bath. II. Ensemble dynamics, Phys. Rev. B 79, 532 (29). 27 L. P. Yang, C. Burk, M. Widmann, S. Y. Lee, J. Wrachtrup, and N. Zhao, Electron spin decoherence in silicon carbide nuclear spin bath, Phys. Rev. B 9, 2423(R) (24). 28 N. Zhao, Z. Wang, and R. B. Liu, Anomalous decoherence effect in a quantum bath, Phys. Rev. Lett. 6, 2725 (2). 29 P. Huang, X. Kong, N. Zhao, F. Z. Shi, P. F. Wang, X. Rong, R. B. Liu, and J. F. Du, Observation of an anomalous decoherence effect in a quantum bath at room temperature, Nature Commun. 2, 57 (2). 3 V. Jacques, P. Neumann, J. Beck, M. Markham, D. Twitchen, J. Meijer, F. Kaiser, G. Balasubramanian, F. Jelezko, and J. Wrachtrup, Dynamic polarization of single nuclear spins by optical pumping of nitrogen-vacancy color centers in diamond at room temperature, Phys. Rev. Lett. 2, 5743 (29). 3 X. K. Song, Q. Ai, J. Qiu, and F.-G. Deng, Physically feasible three-level transitionless quantum driving with multiple Schrödinger dynamics, Phys. Rev. A 93, (26). 32 N. Zhao, J. L. Hu, S. W. Ho, J. T. K. Wan, and R. B. Liu, Atomic-scale magnetometry of distant nuclear spin clusters via nitrogen-vacancy spin in diamond, Nature Nanotech. 6, 242 (2). 33 A. Jarmola, V. M. Acosta, K. Jensen, S. Chemerisov, and D. Budker, Temperature- and magnetic-field-dependent longitudinal spin relaxation in nitrogen-vacancy ensembles in diamond, Phys. Rev. Lett. 8, 976 (22). 34 D. A. Redman, S. Brown, R. H. Sands, and S. C. Rand, Spin dynamics and electronic states of N-V centers in diamond by EPR and four-wave-mixing spectroscopy, Phys. Rev. Lett. 67, 342 (99). 35 A. Jarmola, A. Berzins, J. Smits, K. Smits, J. Prikulis, F. Gahbauer, R. Ferber, D. Erts, M. Auzinsh, and D. Budker, Longitudinal spin-relaxation in nitrogen-vacancy centers in electron irradiated diamond, App. Phys. Lett. 7, (25). 36 M. Mrózek, D. Rudnicki, P. Kehayias, A. Jarmola, D. Budker, and W. Gawlik, Longitudinal spin relaxation in nitrogen-vacancy ensembles in diamond, EPJ Quantum Technology 2, 22 (25). 37 N. Zhao, S. W. Ho, and R. B. Liu, Decoherence and dynamical decoupling control of nitrogen-vacancy centre electron spins in nuclear spin baths, Phys. Rev. B 85, 533 (22). 38 V. V. Dobrovitski, A. E. Feiguin, R. Hanson, and D. D. Awschalom, Decay of Rabi oscillations by dipolarcoupled dynamical spin environments, Phys. Rev. Lett. 2, 2376 (29). 39 W. Zhang, N. Konstantinidis, K. A. A. Hassanieh, and V. V. Dobrovitski, Modelling decoherence in quantum spin systems, J. Phys. Condens. Matter 9, 8322 (27). 4 S. E. Lillie, D. A. Broadway, J. D. A. Wood, D. A. Simpson, A. Stacey, J.-P. Tetienne, and L. C. L. Hollenberg, Environmentally mediated coherent control of a spin qubit in diamond, Phys. Rev. Lett. 8, 6724 (27). 4 A. Jarmola, V. M. Acosta, K. Jensen, S. Chemerisov, and D. Budker, Observation of coherent oscillations in a single electron spin, Phys. Rev. Lett. 92, 764 (24). 42 Q. Ai, T. C. Yen, B. Y. Jin, and Y.-C. Cheng, Clustered geometries exploiting quantum coherence effects for efficient energy transfer in light harvesting, J. Phys. Chem. Lett. 4, 2577 (23). 43 A. Ishizaki and G. R. Fleming, Theoretical examination of quantum coherence in a photosynthetic system at physiological temperature. Proc. Natl. Acad. Sci. U.S.A. 6, 7255 (29).