arxiv: v1 [quant-ph] 8 Nov 2017

Transcription

1 Unambiguous nuclear spin detection using engineered quantum sensing sequence arxiv:7.94v [quant-ph 8 Nov 7 Ziun Shu,, Zhendong Zhang,, Qingun Cao,, Pengcheng Yang,, Martin B. Plenio, 3, Christoph Müller, 4 Johannes Lang, 4 Nikolas Tomek, 4 Boris Nadenov, 4 Liam P. McGuinness, 4 Fedor Jelezko, 4,,, and Jianming Cai School of Phsics, Huazhong Universit of Science and Technolog, Wuhan 4374, China International Joint Laborator on Quantum Sensing and Quantum Metrolog, Huazhong Universit of Science and Technolog, Wuhan, China, Institut für Theoretische Phsik & IQST, Albert-Einstein Allee, Universität Ulm, D-898 Ulm, German 4 Institut für Quantenoptik & IQST, Albert-Einstein Allee, Universität Ulm, D-898 Ulm, German (Dated: November 9, 7 Sensing, localising and identifing individual nuclear spins or frequenc components of a signal in the presence of a nois environments requires the development of robust and selective methods of dnamical decoupling. An important challenge that remains to be addressed in this context are spurious higher order resonances in current dnamical decoupling sequences as the can lead to the misidentification of nuclei or of different frequenc components of external signals. Here we overcome this challenge with engineered quantum sensing sequences that achieve both, enhanced robustness and the simultaneous suppression of higher order harmonic resonances. We demonstrate experimentall the principle using a single nitrogen-vacanc center spin sensor which we appl to the unambiguous detection of external protons. PACS numbers: 76.3.Mi, 76.7.Hb,7.55.Ge Introduction. The detection of single nuclear spins represents an important et challenging step towards single molecule magnetic resonance spectroscop and imaging [ which holds the promise for the observation of individual protein structures without the need for crstallization of large ensembles. The achievement of this goal ma have significant impact on structural biolog and medical imaging and ma also provide a new tool for the investigation of nuclear spin dnamics in non-trivial quantum biological processes [, 3. Recentl, research has made impressive progress in single nuclear spin detection using various phsical platforms [4 8, particularl in the highl challenging task of detecting single spins in external molecules [9 4. Among those phsical platforms, a sensor based on individual negativel charged nitrogen-vacanc (NV centers in diamond [5, 6 has demonstrated appealing prospects for applications in biolog and medicine, due to its biocompatibilit, nano-scale size and long coherence times under ambient conditions [7, 8. NV centers, shallowl implanted to within a few nanometers below the diamond surface, enable strong coupling between NV sensors and target nuclei, thereb promoting the detection sensitivit from a relative large ensemble of nuclei [9, to single nuclear spin sensitivit in small clusters of nuclear spins [ 4. Besides the interest in single molecule magnetic resonance spectroscop and imaging, single nuclear spin addressing also has an important role to pla in the precise coherent control of nuclear spin qubits, where it ma facilitate the realisation of quantum memories with long coherence times and nuclear spin based quantum information processors [9. One ultimate goal of single molecule magnetic resonance spectroscop using a single spin sensor is to detect a single nuclear spin and further infer the structure of a single molecule ianmingcai@hust.edu.cn [3. This would generall require a larger number of dnamical sensing pulses[4 such that even tin pulse errors ma accumulate and deteriorate significantl the magnetic resonance spectroscop signal [5 3. Furthermore, the effect of a finite pulse duration, during which the control pulses and the interaction with the nuclei/field simultaneousl act on the sensor, leads to spurious higher order resonances in the measurement data that ma lead to a misidentification of frequenc components in the signal. Indeed, as a prominent recent example these spurious resonances have recentl been identified as the source for the misidentification of nuclei in work towards nuclear magnetic resonance spectroscop of biological sstems [34. In particular, the ratio of gromagnetic moments µ H /µ 3C 4 of hdrogen ( H and carbon ( 3 C which leads to the overlap of different resonances preventing their unique attribution to specific nuclei [34. The presence of 3 C as part of the target molecule or the diamond itself is difficult to avoid making this an urgent problem to address [36. In this work, we both theoreticall and experimentall demonstrate quantum sensing pulse sequences with engineered phases that effectivel suppress spurious higher order resonances and are, at the same time, robust to pulse errors. In particular, we find that a tpe of MLEV- 8 pulse sequence [37 denoted as YY8 sequence ( π φ [ i N ( (π (π (π (π (π (π (π (π 3π as a representative example can eliminate unwanted higher order spuri- φ i ous harmonic response in nuclear magnetic resonance spectroscop and additionall can exhibit enhanced robustness. Using single nitrogen-vacanc spin sensor in diamond, we demonstrate the simultaneous suppression of both the nd and the 4 th order spurious resonance e.g. from the surrounding 3 C nuclei in diamond in proton sensing. We also observe the effect of the initial phase φ i in the XY8 famil pulse se- quences ( π [ N ( (π φ x (π (π x (π (π (π x (π (π x 3π i φ i [35. B choosing appropriate initial phase φ i, either the nd or

2 a b c d t x x x x ( FIG.. (Color online Robust quantum sensing sequences for single nuclear spin detection. (a XY8 pulse sequence consisting of repetitive [ (π x (π (π x (π (π (π x (π (π x pulses and (b YY8 pulse sequence consisting of [ (π (π (π (π (π (π (π (π pulses where each bar refers to a π rotation around ± ˆx, ±ŷ axis (with t π the pulse duration. (c-d A field with a frequenc of ω nd = ω L / and ω 4 th = ω c /4, where ω c = /(τ with τ the inter pulse time interval, ma lead to spurious high-order resonance signal. the 4 th order spurious resonance ma be suppressed individuall (though not simultaneousl in contrast with YY8 pulse sequence. Based on engineered quantum sensing pulse sequences that suppress spurious signals, we detect the hdrogen signal in an unambiguous wa using shallow implanted NV center spins. Spurious signal and pulse sequence engineering. Without loss of generalit, we take an NV center spin sensor in diamond to illustrate our idea. We remark that the present result is also applicable for quantum sensing with the other twolevel phsical sstems. B choosing e.g. {, }, one can construct a two-level quantum probe for sensing. B appling a sequence of N equidistant (π θk pulses (which represents a π rotation around the axis cos θ k ˆx + sin θ k ŷ on the NV center spin with the inter pulse time interval τ (see Fig., and assuming that the (π θk pulses are instantaneous, the transition probabilit from the initial state φ + = ( + e iφ / to the state φ = ( e iφ / is given b [, 4 P(τ = φ U(Nτ φ +, where φ is the relative phase in the initial state (i.e. the initial phase of the π pulse, U(Nτ is the sstem evolution operator under the oint action of sstem Hamiltonian and pulsed action. In the context of nuclear spin detection, nuclear spins undergo Larmor precession in an external magnetic field, and produce a contribution to the signal with a characteristic frequenc given b the Larmor frequenc. When the pulse interval matches the Larmor frequenc of nuclear spin, i.e. ω L = /(τ, a resonance signal corresponding to the specific nuclear species is observed. The NV center spin state is coherentl manipulated b appling microwave driving fields as H M (t = Ω(t cos (ωt [ cos φ(tσ x + sin φ(tσ, ( where σ = (σ x, σ, σ z is the Pauli operator in the subspace spanned b {, } that are the eigenstates of σ z =. A (θ φ rotation around the axis (cos φ, sin φ, realized in the time duration [t, t + t d is achieved b choosing Ω(t = Ω and φ(t = φ for t [t, t + t d, with t d = θ/ω. In the case that the microwave Rabi frequenc Ω is much larger than the hperfine interaction, the corresponding evolution operator of a (π φ rotation is Uφ π = cos φσ x + sin φσ. Under realistic conditions, due to the finite power of the microwave acting on the NV center, the π-pulses alwas have a finite width (as opposed to instantaneous pulses, which can lead to spurious harmonic response and thus the possible ambiguit in the identification of nuclei [34. Such an effect becomes particularl prominent in singlemolecule magnetic resonance spectroscop as a large number of pulses is required to get a resonance signal from weakl coupled nuclei. To overcome such a serious obstacle, we need to carefull engineer quantum sensing pulse sequence. In order to illustrate the basic idea, we first consider a classical oscillating field H ac = A ac sin(ω ndt σ z with a frequenc ω nd = ω L/, see Fig.(c. We prepare the NV spin sensor into the initial state as φ + = ( + e iφ /. As the pulse has a finite width, the corresponding evolution of the spin sensor state is not a perfect π-rotation around ˆx or ŷ. Instead, the real evolution can be described b a set of rotations in the form of Rˆnk = exp ( iθ k ˆn k σ/, where σ is the spin operator and θ k = π + π β k + O(β3 k, ˆn k = (cos β k cos φ k, cos β k sin φ k, sin β k and β k A ac sin(ωt k /Ω, t k is the moment when the pulse is applied. Such imperfection in pulses give rise to spurious resonance signal. For the XY8 pulse sequence, in the scenario that we are interested in, namel A ac /Ω, the total evolution for one ccle is U nd = 8 k= Rˆn k = ( 8β Î + i β (σ x σ, where β = A ac /Ω. This leads to a nd harmonic response signal P nd = 8 [ sin(φ β + O(β3 [34, 35, which spuriousl indicates a frequenc component that is twice of the real oscillating field frequenc. Similar analsis shows that the 4 th harmonic response signal is P 4 th = ( [ + sin(φ β + O(β3 [34, 35. It can be seen that b changing the initial phase φ, nd and 4 th harmonics can be suppressed respectivel, as P nd vanishes when φ = π/4 while P 4 th goes to zero when φ = 3π/4. This ma require extra experiment complexit. We find that a tpe of MLEV-8 pulse sequence [ N(for (π (π (π (π (π (π (π (π simplicit we call it as YY8 pulse sequence can suppress both nd and 4 th harmonic signals simultaneousl and thus reduce the experiment complexit. Using YY8 pulse sequence, the ccle evo-

3 FIG.. (Color online The spectrum signal of the 3 C nuclear spins observed b an NV center using XY8- and YY8- pulse sequence with different initial phases φ = (, π/ (, π (. The Rabi frequenc is Ω = (π5mhz, the magnetic field is 5G. lution of the NV center spin state for the nd and 4 th harmonic response are given b: P nd = 4 [ cos(φ β + O(β3 P 4 th = ( [ cos(φ β + O(β3 ( It can be seen that the nd and 4 th harmonic signal can thus be simultaneousl suppressed up to second-order when the initial phase is chosen as φ = [38. We remark that YY8 pulse sequence can suppress the other spurious resonances, including the (4/5 th, (/3 th and (4/3 th response, in contrast with XY8 pulse sequence [38. Unambiguous nuclear spin detection. Using an NV center spin sensor in diamond, we demonstrate experimentall the suppression of spurious high-order response signals in the detection of nuclear spins using engineered quantum sensing pulse sequences. The NV center is a point defect in diamond, which has a triplet ground state manifold with the spin quantum number m s =, ±. In NV-based nuclear magnetic resonance experiments, one applies a magnetic field B z along the NV axis that connects the nitrogen atom and the vacanc site, the NV center spin Hamiltonian is given as follows [5 H = DS z γ e B z S z, (3 where the zero-field splitting is D/π =.87GHz, γ e is the electron gromagnetic ratio, S is the spin- operator. As an example, we consider the detection of 3 C nucleus in bulk diamond. The hperfine interaction between the NV center spin and the nuclear spin is described b H I = S z A I = S z + a ˆ I x, (4 (a Îz where a and a are the longitudinal and transversal components of the hperfine interaction. The hperfine interaction will cause an effective fluctuating magnetic field with a characteristic frequenc that is equal to the Larmor frequenc of the nuclear spin. The transition probabilit of the NV[ spin state at time T is thus given b P(T = φ Tr I U( φ + φ + ρ N U φ, where U is the evolution operator, ρ N represents the initial state of the nuclear spins. The probabilit that the NV spin remains in the initial state P = P(Nτ is the signal observed in our NV-based nuclear magnetic resonance experiments, which provides us information on the nuclear spin, such as the tpe of nuclei and the hperfine interaction with the NV spin sensor. In our experiment, we polarize and read out the NV center spin in tpe IIa diamond (with natural abundance of 3 C using optical techniques. The microwave control pulses are generated using an arbitrar wave generator (Tektronix AWG7A and is amplified b an 8W microwave amplifier using a same channel. This allows us to avoid the amplitude imbalance between pulses, which otherwise would affect the robustness of YY8 pulse sequence [39. Therefore, YY8 pulses ma be more robust than XY8 sequence as is also evidenced b our numerical simulation [38. The amplitude and phase are tuned to implement the engineered quantum sensing pulse sequences. We appl a magnetic field (5G along the NV axis, under which condition the associated 4 N can be polarized to improve the sensing signal, and use m s = and m s = as quantum probe in our sensing experiments. We first measure the magnetic resonance signal of 3 C nuclear spin using XY8 and YY8 pulse sequences, see Fig.. It can be seen that the signal is almost independent on the initial phase, and both XY8 and YY8 pulse sequences demonstrate similar (I XY8 Sensing Sequence (II YY8 Sensing Sequence /3 /3 /3 /3 FIG. 3. (Color online Suppression of spurious high-order resonance signals using engineered quantum sensing pulse sequences. (a-b The nd and 4 th harmonic response using XY8-6 and XY8- pulse sequence with the initial phase φ = π/6 ( and 3π/4 ( respectivel. (c The amplitude of the nd ( and 4 th ( harmonic signal (i.e. the resonance contrast as a function of the initial phase using XY8-6 and XY8- pulse sequence. (d-e The nd and 4 th harmonic response using YY8-6 and YY8- pulse sequence with the initial phase φ = ( and π/ ( respectivel. (f The amplitude of the nd ( and 4 th ( harmonic signal as a function of the initial phase using YY8-6 and YY8- pulse sequence. The Rabi frequenc is Ω = (π5 MHz. Both spurious harmonics can be suppressed simultaneousl for YY8 pulse sequences which contrasts with XY8 pulse sequences.

4 FIG. 4. (Color online Detection of H magnetic resonance signal using YY8 pulse sequence. (a The spectrum as measured b a shallow NV center spin for various magnetic field. (b The frequenc of peak response as a function of the magnetic field B. The red line has a slope 4. ±.3 khz/g that gives an estimated gromagnetic ratio of H. The Rabi frequenc is Ω = (π5 MHz. performance. We further test the suppression of spurious high-order resonance signals using engineered pulse sequences. In Fig.3, we plot the dependence of the nd and 4 th harmonic signal on the initial phase φ using the YY8 pulse (d-f sequence as compared with the XY8 sequence (a-c. It can be seen that the spurious high-order resonance signals can be suppressed b choosing approximate initial phases of XY8 and YY8 pulse sequences, in which YY8 pulse sequence demonstrates better performance over XY8 pulse sequence. B choosing φ =, both nd and 4 th harmonics are suppressed simultaneousl using YY8 pulse sequence, see Fig.3 (d-f, while this cannot be achieved for an choice of φ for the XY8 sequence, see Fig.3(a-c. The amplitude of spurious high-order resonance signals are also relativel smaller using YY8 pulse sequence even when the initial phase is not perfect. These advantageous features of YY8 pulse sequence allows for the efficient unambiguous identification of nuclear spins, in particular it implies that YY8 pulse sequence can unambiguousl identif the signal from H, while avoiding the possible confusion with the 4 th harmonic signal of 3 C with no extra experiment requirement as compared with the alread widel used XY8 pulse sequence. We then use YY8 pulse sequence with the initial phase φ = to detect the magnetic resonance signal of H in oil on diamond surface. The shallowed implanted NV centers are created b a.5 kev N + ions into tpe IIa diamond. The Ramse and spin echo measurement indicates that the NV center shows T 3ns and T 5µs. In Fig.4(a, we measure the spectrum recorded b YY8 pulse sequence at different magnetic fields. The scaling of the peak frequenc with the magnetic field is shown in Fig.4(b. The linear slope of the frequenc 4. ±.3kHz/G fits well with the gromagnetic ratio of H (4.58 khz/g. As we show above, the YY8 pulse sequence can effectivel suppress the high-order spurious response, and therefore our measurement provides a wa to detect H unambiguousl. As H magnetic resonance spectroscop is ver important in various scenarios, the present engineered quantum sensing pulse sequence, particularl YY8 pulse sequence, ma find useful applications in single-molecule magnetic resonance spectroscop using quantum spin sensor. Conclusions We implement engineered quantum sensing pulse sequences that suppress spurious higher order resonances, a feature that will find application in unambiguous nuclear spin detection and more generall the identification of frequenc components of externall applied signals. Specificall, this eliminates the ambiguit in the nuclear magnetic resonance signal of single nuclear spin detection that is inherent to pulsed schemes due to the duration of individual pulses. The YY8 pulse sequences shows better performance in suppressing spurious response and are even more robust against the pulse errors than the widel used XY8 pulse sequences (e.g amplitude error. The simplicit and advantageous features of this tpe of pulse sequences make them an appealing choice for the applications in nuclear spin detection as well as the quantum control of nuclear spin registers. Acknowledgements We thank Jan Haase and Benedikt Tratzmiller for helpful comments and suggestions, and Michael Ferner and Manfred Bürzele for technical assistance. J.-M.C is supported b National Natural Science Foundation of China (5743, 693, 693. M. B. P. is supported b the DFG (FOR493 and SFB TR, the EU STREPs EQUAM, DIADEMS and HYPERDIAMOND, and the ERC Snerg grant BioQ. C. M., J. L., N. T., B. N., L. P. M., F. J. acknowledge support from DFG (FOR 493, SFB TR, SPP 93, VW Stiftung, BMBF, ERC, EU (DI- ADEMS, BW Stiftung, Ministr of Science and Arts, Center for Integrated quantum science and technolog (IQST. [ D. Rugar, R. Budakian, H. J. Mamin, and B. W. Chui, Nature 43, 39 (4 [ S. F. Huelga and M. B. Plenio. Contemp. Phs. 54, 8-7 (3. [3 N. Lambert, Y.-N. Chen, Y.-C. Cheng, C.-M. Li, G.-Y. Chen and F. Nori, Nature Phsics 9, (3. [4 R. Vincent, S. Klatskaa, M. Ruben, W. Wernsdorfer, and F. Balestro, Nature 488, 357 (. [5 J. J. Pla, K. Y. Tan, J. P. Dehollain, W. H. Lim, J. J. L. Morton, F. A. Zwanenburg, D. N. Jamieson, A. S. Dzurak, and A. Morello, Nature 496, 334 (3. [6 T. H. Taminiau, J. J. T. Wagenaar, T. van der Sar, F. Jelezko, V. V. Dobrovitski, and R. Hanson, Phs. Rev. Lett. 9, 376 (. [7 S. Kolkowitz, Q. P. Unterreithmeier, S. D. Bennett, and M. D. Lukin, Phs. Rev. Lett. 9, 376 (. [8 N. Zhao, J. Honert, B. Schmid, M. Klas, J. Isoa, M. Markham, D. Twitchen, F. Jelezko, R.-B. Liu, H. Fedder, and J. Wrachtrup, Nature Nanotechnolog 7, 657 (. [9 H. J. Mamin, M. Kim, M. H. Sherwood, C. T. Rettner, K. Ohno, D. D. Awschalom, and D. Rugar, Science 339, 557 (3. [ T. Staudacher, F. Shi, S. Pezzagna, J. Meier, J. Du, C. A. Meriles, F. Reinhard, and J. Wrachtrup, Science 339, 56 (3. [ C. Müller, X. Kong, J.-M. Cai, K. Melentievic, A. Stace, M.

5 5 Markham, J. Isoa, S. Pezzagna, J. Meier, J.-F. Du, M. B. Plenio, B. Nadenov, L. P. McGuinness and F. Jelezko. Nature Communications 5, 473 (4. [ A. O. Sushkov, I. Lovchinsk, N. Chisholm, R. L. Walsworth, H. Park, and M. D. Lukin, Phs. Rev. Lett. 3, 976 (4. [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, J.-F. Du, Science 347, 35 (5. [4 V. S. Perunicic, L. T. Hall, D. A. Simpson, C. D. Hill, L. C. L. Hollenberg, Phs. Rev. B 89, 5443 (4. [5 M. W. Dohert, N. B. Manson, P. Delane, F. Jelezko, J. Wrachtrup and L. C.L. Hollenberg, Phs. Reports 58, (3. [6 Y. Wu, F. Jelezko, M. B. Plenio, and T. Weil. Angewandte Chemie- International Edition Minireview 55, 6586 (6. [7 R. Schirhagl, K. Chang, M. Loretz, and C. L. Degen, Annu. Rev. Phs. Chem (4. [8 A. Ermakova, G. Pramanik, J.-M. Cai, G. Algara-Siller, U. Kaiser, T. Weil, H. C. Chang, L.P. McGuinness, M. B. Plenio, B. Nadenov and F. Jelezko, Nano Letters 3, 335 (3. [9 L. Jiang, J. S. Hodges, J. R. Maze, P. Maurer, J. M. Talor, D. G. Cor, P. R. Hemmer, R. L. Walsworth, A. Yacob, A. S. Zibrov, and M. D. Lukin, Science 36, 67 (9. [ P. C. Maurer, G. Kucsko, C. Latta, L. Jiang, N. Y. Yao, S. D. Bennett, F. Pastawski, D. Hunger, N. Chisholm, M. Markham, D. J. Twitchen, J. I. Cirac, and M. D. Lukin, Science 336, 83 (. [ G. Waldherr, Y. Wang, S. Zaiser, M. Jamali, T. Schulte- Herbrueggen, H. Abe, T. Ohshima, J. Isoa, P. Neumann, and J. Wrachtrup, Nature 567, 4 (4. [ J.-M. Cai, A. Retzker, F. Jelezko, and M. B. Plenio, Nature Phsics 9, 68 (3. [3 J.-M. Cai, F. Jelezko, M. B. Plenio, A. Retzker, New J. Phs. 5, 3 (3. [4 L. Cwiński, R.M. Lutchn, C.P. Nave, and S. Das Sarma, Phs. Rev. B 77, 7459 (8. [5 N. Zhao, J. Wrachtrup, R.-B. Liu, Phs. Rev. A 9, 339 (4. [6 W.-C. Ma, F.-Z. Shi, K.-B. Xu, P.-F. Wang, X.-K. Xu, X. Rong, C.-Y. Ju, C.-K. Duan, N. Zhao, and J.-F. Du, Phs. Rev. A 9, 3348 (5. [7 J. N. Greiner, D. D. Bhaktavatsala Rao, P. Neumann, J. Wrachtrup, arxiv: [8 A. Albrecht and M. B. Plenio, Phs. Rev. A 9, 34 (5 [9 J. Casanova, Z.-Y. Wang, J. F. Haase, and M. B. Plenio, Phs. Rev. A 9, 434 (5. [3 Z.-Y. Wang, J. F. Haase, J. Casanova, and M. B. Plenio, Phs. Rev. B 93, 744 (6. [3 T. Gullion, D. B. Baker, and M. S. Conradi, J. Magn. Reson. 89, 479 (99. [3 C. A. Ran, J. S. Hodges, and D. G. Cor, Phs. Rev. Lett. 5, 4 (. [33 J.-M. Cai, B. Nadenov, R. Pfeiffer, L. McGuinness, K. D. Jahnke, F. Jelezko, M. B. Plenio and A. Retzker, New J. Phs. 4, 33 (. [34 M. Loretz, J. M. Boss, T. Rosskopf, H. J. Mamin, D. Rugar, and C. L. Degen, Phs. Rev. X 5, 9 (5. [35 J. F. Haase, Z.-Y. Wang, J. Casanova, and M. B. Plenio, Phs. Rev. A 94, 33 (6 [36 T. Staudacher, N. Raatz, S. Pezzagna, J. Meier, F. Reinhard, C. A. Meriles, J. Wrachtrup, Nat. Commun. 6, 857 (5. [37 M. H. Levitt, R. Freeman, T. Frenkiel, in Advances in Magnetic Resonance (J. S. Waugh, Ed., Academic Press, New York, 983. [38 See supplementar material for the calculation details, which includes Refs. [4 46 [39 M. H. Levitt, R. Freeman, and T. Frenkiel, Advances in Magnetic Resonance (J. S. Waugh, Ed., Vol., Academic Press, New York (983. [4 G. de Lange, Z. H. Wang, D. Ristè, V. V. Dobrovitski, R. Hanson, Science, 6 (. [4 Z.-H. Xiao, L.-W. He, and W.-G. Wang.Phs.Rev.A83, 33 (. [4 I. Solomon, Phs. Rev. Lett., 3 (959. [43 C. D. Aiello, M. Hirose, P. Cappellaro, Nature Communications 4, 49 (3. [44 V. V. Mkhitaran and V. V. Dobrovitski, Phs. Rev. B 89, 44 (4. [45 V. V. Mkhitaran, F. Jelezko, V. V. Dobrovitski, Scientific Reports 5, 54 (5. [46 R. Cabrera, W.E. Balis, Phs. Lett. A. 368, 5 (7.