Chinese Physics Letters, 2021, Vol. 38, No. 4, Article code 044201 A Noise-Robust Pulse for Excitation Transfer in a Multi-Mode Quantum Memory Bo Gong (龚波), Tao Tu (涂涛)*, Ao-Lin Guo (郭奥林), Le-Tian Zhu (朱乐天), and Chuan-Feng Li (李传锋)* Affiliations Key Laboratory of Quantum Information, University of Science and Technology of China, Chinese Academy of Sciences, Hefei 230026, China Received 16 December 2020; accepted 3 February 2021; published online 6 April 2021 Supported by the National Natural Science Foundation of China (Grant No. 11974336), and the National Key Research and Development Program of China (Grant No. 2017YFA0304100).
*Corresponding authors. Email: tutao@ustc.edu.cn; licf@ustc.edu.cn
Citation Text: Gong B, Tu T, Guo A L, Zhu L T, and Li C F 2021 Chin. Phys. Lett. 38 044201    Abstract Multi-mode quantum memory is a basic element required for long-distance quantum communication, as well as scalable quantum computation. For on-demand readout and long storage times, control pulses are crucial in order to transfer atomic excitations back and forth into spin excitations. Here, we introduce noise-robust composite pulse sequences for high-fidelity excitation transfer in multi-mode quantum memory. These pulses are robust to the deviations in amplitude and the detuning parameters of realistic conditions. We show the efficiency of these composite pulses with a typical rare-earth ion-doped system. This approach could be applied to a variety of quantum memory schemes. DOI:10.1088/0256-307X/38/4/044201 © 2021 Chinese Physics Society Article Text Quantum memory comprises specially designed quantum processors which exploit the storage and manipulation of flying qubits and stationary qubits. They constitute the essential building blocks in various forms of quantum information processing: facilitating the synchronization of probabilistic photonic processes for quantum computing, functioning as quantum repeaters for long-distance quantum communication.[1–6] The implementation of quantum memory requires efficient interaction between flying qubits (photons) and stationary qubits (matter qubits). A variety of physical systems have been demonstrated as platforms for quantum memories, such as atomic ensembles,[7] single atoms in cavities,[8] spins in quantum dots[9,10] or NV centers,[11] and rare-earth ions doped in solids.[12] Rare-earth ion-doped systems are promising candidates for use in quantum memory, since on the one hand, they can achieve great efficiency by mapping photons onto collective atomic excitations, and on the other hand, solid state systems offer the potential for integration and scalability in conjunction with existing technology. Quantum memory offers a variety of advantages, in areas including storage time, fidelity, multi-mode capacity, and on-demand readout.[1,3] Rare-earth ions in solids exhibit excellent coherent properties in terms of spin states. Furthermore, they demonstrate a significant inhomogeneous broadening of optical transitions, which could be a useful resource for multiplexed quantum protocols, e.g., temporal or spectral multi-mode storage. Recent experimental advances in this system include a long coherence time, with time scales measured in hours,[13] and the mapping and retrieval of 64 multiplexed temporal modes,[14] or 26 multiplexed spectral modes,[15] as well as the storage of entanglement,[16,17] polarization[18–20] and the orbital angular momentum of photons.[21] The atomic frequency combs spin-wave (AFC-SW) protocol has been proposed as a useful resource, combining massive multiplexing, long storage time, and the capacity for on-demand readout of stored states.[22–24] On-demand readout can be achieved by transferring optical atomic excitations to collective spin excitations using control pulses. This also results in a much longer storage time, since the coherence time of spin excitations is far superior to that of atomic excitations. As such, the performance of quantum memory based on the AFC-SW protocol is significantly determined by the control pulses to transfer the excitations. In particular, the fidelity of the control pulses is crucial in applications where quantum memory is used as a quantum processing node in a network.[6] However, the achievement of fully functional quantum memory with high fidelity has so far remained elusive, due to limitations such as noise, and imperfections in the control pulses. In this work, we propose a scheme to design noise-robust pulses for excitation transfer with high fidelity. Firstly, we provide an overview of the AFC-SW protocol, focusing on the transfer process. We then develop a control model, taking into account the main sources of noise. We utilize a composite pulse approach to achieve the desired transfer operations, via analytic formulas and numerical simulations. This type of pulses can realize the desired population transfer between atomic excitations and spin excitations with high fidelity, even under realistic noise conditions. Many works have addressed the issue of composite pulses. The most common composite pulses are designed by employing various techniques, such as the Magnus expansion,[25,26] quaternion algebra, and geometric phases.[27] In addition, most composite pulses are designed to achieve particular classes of gate operations,[25–27] and to tackle a specific source of errors.[25,26] In contrast, our pulse sequences are constructed in a simple iterative manner, can achieve general quantum gate operations, and can systematically eliminate both amplitude and detuning errors. Therefore, our composite pulses are flexibly applicable to various quantum memory protocols.
cpl-38-4-044201-fig1.png
Fig. 1. Schematic of AFC-SW-based quantum memory. A signal photon field induces the optical transition between atomic levels $|g\rangle $ and $|e\rangle $. Next, a pair of control pulses transfer the excitations between levels $|e\rangle $ and $|s\rangle $, thereby facilitating on-demand readout and long-term storage. The capacity for multi-mode storage is based on the spectrally shaped combs on $|e\rangle $, which are composed of narrow peaks with width $\gamma $ and separation $\varDelta $.
We consider the typical system of rare-earth ions, doped in crystal, such as Pr$^{3+}$:Y$_{2}$SO$_{5}$,[23] or Eu$^{3+}$:Y$_{2}$SO$_{5}$.[24] Such a system can be described as an ensemble of atoms with three levels, as shown in Fig. 1. Since the levels are inhomogeneously broadened, the transition between the ground state $|g\rangle $ and the excited state $|e\rangle $ can be spectrally shaped. The atomic distribution could therefore consist of a frequency comb, where narrow peaks with a characteristic width $\gamma $ are separated by $\varDelta $, and span a large frequency range $\varGamma $. When a signal photon enters the crystal, it should be absorbed by the crystal. It would then be stored in an excitation mode, delocalized over all the atoms, which could be expressed by the Dicke-type state, as given in Ref. [22]: $\sum_{j=1}^{N}e^{-i\upsilon _{j}t}e^{-ik_{\rm s}z_{j}}|g_{1}\dots e_{j}\dots g_{N}\rangle $. Here, $N$ is the total number of atoms in the crystal, $\upsilon _{j}$ is the detuning parameter, $z_{j} $ is the position, $j$ corresponds to the $j$th atom, and $k_{\rm s}$ is the wave number of the absorbed signal photon field. If the comb peaks are sharp enough, the detuning parameter $\upsilon _{j}$ can be approximated as $\upsilon _{j}=m_{j}\varDelta $. After an evolution time $2\pi /\varDelta $, all the atomic components are in phase, leading to the collective emission of stored photons. The multi-mode capacity of this type of quantum memory can be implemented as a train of input signal modes. The duration of one mode is limited by the frequency range of the comb as $12\pi /\varGamma $, and the duration of the mode train is limited by the rephasing time as $2\pi /\varDelta $. Thus the total number of the input modes that can be stored could be given by the ratio between the train duration and the duration of one mode, as $n=\varGamma /6\varDelta $. In rare-earth ionic materials, it is easy to prepare hundreds of comb peaks (i.e., small $\varDelta$ with large $\varGamma $) making it possible to store many modes. As shown in Fig. 1, once the signal photon is absorbed, a “write” pulse is applied to transfer the atomic excitations in the $|g\rangle $ and $|e\rangle $ levels to the spin excitations in the $|g\rangle $ and $|s\rangle $ levels.[22,28] Since the spin states have a much longer coherence time, $T_{\rm s}$, than the optical coherence time, the protocol allows for long-term storage. Next, a “read” pulse is applied to transfer the spin excitations back to atomic excitations, and a re-emission of the stored photon occurs at time $2\pi /\varDelta +T_{\rm s}$. In this way, the whole process leads to on-demand storage and retrieval capability in the signal photon. The control pulses for excitation transfer play an important role in the AFC-SW protocol, which determines the remarkable features of quantum memory, i.e., the capacity for long-term storage and on-demand readout. In the AFC-SW protocol, the control pulses are applied to transfer collective excitations between the $|e\rangle $ and $|s\rangle $ levels. Under the rotating-wave approximation, the Hamiltonian of the system could be expressed as $$ H(t)=\frac{\omega (t)}{2}\sigma _{x}+\frac{\upsilon (t)}{2}\sigma _{z},~~ \tag {1} $$ where $\sigma _{x}$ and $\sigma _{z}$ are the Pauli operators, defined in the subspace of $|e\rangle $ and $|s\rangle $; $\upsilon (t)$ and $\omega (t)$ are the detuning and the amplitude of the control pulse, respectively. Thus, the dynamics of the system is determined by the propagator $$\begin{alignat}{1} R(\omega \hat{x}&+v\hat{z},\theta) =\exp \Big[-i\int_{0}^{t}dtH(t)\Big] \\ &=\exp \Big[-i\Big(\frac{\omega }{2}\sigma _{x}+\frac{\upsilon }{2}\sigma _{z}\Big)\frac{\theta }{\sqrt{\omega ^{2}+\upsilon ^{2}}}\Big].~~ \tag {2} \end{alignat} $$ Here $\theta $ is the rotating angle around the axis $\omega \hat{x}+v\hat{z}$ on the Bloch sphere. We begin with the case of the widely-used $\pi $ pulse. If the detuning maintains a value of $\upsilon (t)=0$, the propagator has an explicit form, as $$ R(\hat{x},\theta)= \begin{pmatrix} \cos (A/2) & -i\sin (A/2) \\ -i\sin (A/2) & \cos (A/2)\end{pmatrix} ,~~ \tag {3} $$ where $A=\int_{0}^{t}\omega (t)dt$ is the pulse area. When we set $A=\pi $, we achieve a $\pi $ pulse capable of realizing excitation transfer between the $|e\rangle $ and $|s\rangle $ levels. This pulse is a rotation around $x$ axis on the Bloch sphere. Since the $\pi $ pulse explicitly depends on the control parameters, it is not robust with respect to variations in these parameters. As such, the $\pi $ pulse struggles to achieve high fidelity under realistic conditions. In the following, we outline the following scheme to design noise-robust pulses for excitation transfer: Step 1: Choose the target operation. As described above, quantum memory based on the AFC-SW protocol utilizes a range of frequency detuning to produce the atomic frequency comb structure. As such, we require a control pulse at a fixed detuning, $v_{\rm c}$. The ordinary $\pi $ pulse around the $x$ axis is denoted by $R(\hat{x},\pi)$. We could achieve $R(\hat{x},\pi)$ via $$ R(\hat{x},\pi)=R(\hat{x}+\hat{z},\pi)R(\hat{z},\pi)R(\hat{x}+\hat{z},\pi).~~ \tag {4} $$ To implement a rotation $R(\hat{x}+\hat{z},\pi)$ around the axis $\hat{x}+\hat{z}$ with a $\pi $ angle, we could set $\omega (t)=\upsilon _{\rm c}$, and the time duration as $\pi /\sqrt{\omega (t)^{2}+\upsilon _{\rm c}^{2}}$. For a rotation $R(\hat{z},\pi)$ around the axis $\hat{z}$ with a $\pi $ angle, we could set $\omega (t)=0$ and time duration as $\pi /\sqrt{0^{2}+\upsilon _{\rm c}^{2}}$. Putting these three rotations together, we can obtain the target $\pi $ pulse around the $x$ axis. Step 2: Model the noise effects. There are two types of noise: fluctuation $\delta \omega $ in the pulse amplitude $\omega $, and fluctuation $\delta \upsilon $ in the detuning $\upsilon $. Inserting two kinds of noises into the Hamiltonian, we get $$ H^{1}(t)=\frac{\omega (t)+\delta \omega }{2}\sigma _{x}+\frac{\upsilon +\delta \upsilon }{2}\sigma _{z}.~~ \tag {5} $$ Since both amplitude and frequency noise are typically orders of magnitude slower than operation times, the resulting perturbations about $\omega $ and $v$ are treated as quasi-static noise in the following analysis. The evolution operator is now rewritten as $$\begin{alignat}{1} &U(\omega \hat{x}+v\hat{z},\theta)\\ ={}&\exp \Big[-i\int_{0}^{t}dtH^{1}(t)\Big] \\ ={}&\exp \Big[-i\Big(\frac{\omega +\delta \omega }{2}\sigma _{x}+\frac{\upsilon +\delta \upsilon }{2}\sigma _{z}\Big)\frac{\theta }{\sqrt{\omega ^{2}+\upsilon ^{2}}}\Big].~~ \tag {6} \end{alignat} $$
Table 1. Parameters for the designed composite pulse.
Rotation parameters $\omega_{0}$ $\omega_{1}$ $\omega_{2}$ $\omega_{3}$ $\omega_{4}$ $\omega_{5}$
Values (MHz) 8.1028 43.2337 0 30.9351 3.9232 0
In order to address the effects of noise, we need to separate the noise terms from the operation term. We therefore perform an expansion of the evolution operator to the first order of noise, given as $$\begin{alignat}{1} U(\omega \hat{x}+v\hat{z},\theta)={}&\exp \Big[-i\Big(\frac{\omega }{2}\sigma _{x}+\frac{\upsilon }{2}\sigma _{z}\Big)\frac{\theta }{\sqrt{\omega ^{2}+\upsilon ^{2}}}\Big] \\ &\times (I-i\sum_{k}\varDelta _{k}\sigma _{k}) \\ ={}&R(\omega \hat{x}+v\hat{z},\theta)(I-i\sum_{k}\varDelta _{k}\sigma _{k}).~~ \tag {7} \end{alignat} $$ Here, $R(\omega \hat{x}+v\hat{z},\theta)$ is the desired rotation, $I$ is the identity operator, $\sigma _{k}$ are the Pauli operators, and the error terms are given as $$\begin{alignat}{1} \varDelta _{x}& =\delta \omega \frac{\omega ^{2}\theta +\upsilon ^{2}\sin \theta }{2(\omega ^{2}+\upsilon ^{2})^{3/2}}+\delta \upsilon \frac{\omega \upsilon (\theta -\sin \theta)}{2(\omega ^{2}+\upsilon ^{2})^{3/2}}, \\ \varDelta _{y}& =\delta \omega \frac{\upsilon (\cos \theta -1)}{2(\omega ^{2}+\upsilon ^{2})}+\delta \upsilon \frac{\omega (1-\cos \theta)}{2(\omega ^{2}+\upsilon ^{2})}, \\ \varDelta _{z}& =\delta \omega \frac{\omega \upsilon (\theta -\sin \theta)}{2(\omega ^{2}+\upsilon ^{2})^{3/2}}+\delta \upsilon \frac{\upsilon ^{2}\theta +\omega ^{2}\sin \theta }{2(\omega ^{2}+\upsilon ^{2})^{3/2}}.~~ \tag {8} \end{alignat} $$ These formulas demonstrate how noise affects the evolution of the system, resulting in small errors in the target rotation, which are characterized as error terms $\varDelta _{k}$ in three degrees of freedom, $x,y,z$. Step 3: Construct a composite pulse to cancel the noise effects. We have a rotation around the $\hat{z}$ axis, together with two rotations around the $\hat{x}+\hat{z}$ axis to achieve the desired rotation, but these three rotations would induce various noise effects during the rotation process. How to eliminate noise is the key question in the field of quantum information processing. To resolve this problem, we exploit the idea of the composite pulse, as given in Refs. [29–32], to construct a composite $2\pi $ rotation. The composite $2\pi $ rotation can be viewed as a series of rotations with different axes and $2\pi $ angles. Specifically, we construct a five-level pulse sequence for the composite $2\pi $ rotation as: $$\begin{alignat}{1} \tilde{I}^{(5)} ={}&\,U\Big(\omega _{5}\hat{x}+v_{\rm c}\hat{z},\pi +\frac{3\pi }{2}\Big) \\ &\times U(\omega _{4}\hat{x}+v_{\rm c}\hat{z},\pi)U(\omega _{3}\hat{x}+v_{\rm c}\hat{z},\pi) \\ &\times U(\omega _{2}\hat{x}+v_{\rm c}\hat{z},\pi)U(\omega _{1}\hat{x}+v_{\rm c}\hat{z},\pi) \\ &\times U(\omega _{0}\hat{x}+v_{\rm c}\hat{z},2\pi) \\ &\times U(\omega _{1}\hat{x}+v_{\rm c}\hat{z},\pi)U(\omega _{2}\hat{x}+v_{\rm c}\hat{z},\pi) \\ &\times U(\omega _{3}\hat{x}+v_{\rm c}\hat{z},\pi)U(\omega _{4}\hat{x}+v_{\rm c}\hat{z},\pi) \\ &\times U\Big(\omega _{5}\hat{x}+v_{\rm c}\hat{z},\frac{\pi }{2}\Big).~~ \tag {9} \end{alignat} $$ Combined with this composite $2\pi $ rotation, we can then find the entire desired rotation, expressed as $$ U_{{\rm NRP}}=U(\hat{z}+\hat{x},\pi)\tilde{I}^{(5)}U(\hat{z},\pi)U(\hat{z}+\hat{x},\pi).~~ \tag {10} $$ This composite pulse sequence would also give rise to noise, as $$ U_{{\rm NRP}}=R(\hat{x},\pi)\Big(I-i\sum_{k}\eta _{k}\sigma _{k}\Big),~~ \tag {11} $$ where $\eta _{k}$ are error terms for the composite pulse sequence due to fluctuations in $\delta \omega $ and $\delta v$. We then select the five control parameters, $\omega _{i}$, to ensure that $\eta _{k}=0$. As shown in Eqs. (10) and (11), the composite $2\pi $ rotation plays an important role in this scheme: Firstly, it offers free parameters $\omega _{i}$ to allow for noise cancellation of the whole pulse sequence. Secondly, it is similar to an identity operator, in that it does not affect the overall target rotation. For concreteness, we consider a detailed example, based on the AFC-SW protocol in Pr$^{3+}$:Y$_{2}$SO$_{5}$,[23] or Eu$^{3+}$:Y$_{2}$SO$_{5}$.[24] On the one hand, the commonly used $\pi $ pulse for excitation transfer is shown in Fig. 2(a). The pulse amplitude is $\omega (t)= \pi $ MHz, and the pulse detuning is $v(t)=0$. On the other hand, we choose the frequency detuning parameter $v_{\rm c}=4\pi $ MHz, and numerically solve the coupled equations for the amplitude parameters, $\omega _{i}$. The results are shown in Table 1. With these parameter values, we are able to construct a noise-robust composite pulse for excitation transfer, as illustrated in Fig. 2(b). Figure 2(c) shows the population probability on the $|e\rangle $ and $|s\rangle $ levels as a function of composite pulse duration.
cpl-38-4-044201-fig2.png
Fig. 2. (a) Commonly used $\pi $ pulse with $\upsilon (t)=0$ and $\omega (t)=\pi $ MHz. (b) Noise robust pulse with detuning $v_{\rm c}=4\pi $ MHz. (c) Population transfer between levels $|e\rangle $ and $|s\rangle $, using the noise robust pulse.
Transfer fidelity is defined as the overlap between the ideal transferred state, and the realistic transferred state under conditions of noise. We use the formula $$ F={\rm Tr}\sqrt{\sqrt{\rho '}\rho \sqrt{\rho '}},~~ \tag {12} $$ where $\rho $ is the density matrix, as $\rho =|\psi \rangle \langle \psi |$ without noise, and $\rho '=|\psi '\rangle \langle \psi '|$ is the density matrix including noise. We calculate the infidelity $1-F$ as a function of amplitude noise $\delta \omega $ and detuning noise $\delta v$. We note that in realistic situations, the ratio of noise to the control field ($\frac{\delta \omega }{\omega}$ or $\frac{\delta \upsilon }{\upsilon}$) is between $0$ and $0.2$. Figure 3 shows the infidelity for the commonly used $\pi $ pulse (denoted as OP), and the noise-robust $\pi $ pulse (denoted as NRP), under different noise conditions. In Figs. 3(a) and 3(a) small variation $\delta \omega $ in the detuning induces an error in transfer fidelity on the order of $10^{-2}$ using the ordinary $\pi $ pulse, whereas the transfer infidelity with the noise-robust $\pi $ pulse remains one to two orders smaller than the ordinary $\pi $ pulse over a wide range of noise values. The noise-robust $\pi $ pulse permits fidelity beyond $99\%$ even with respect to a large noise value of $\delta \omega /\omega =0.2$. This high fidelity value is essential in many quantum network applications involving local operations requiring to be performed fault-tolerantly.[6] Figure 3(b) shows the infidelity as a function of driving detuning noise, where the noise-robust $\pi $ pulse remains more robust than the ordinary $\pi $ pulse.
cpl-38-4-044201-fig3.png
Fig. 3. Transfer infidelity versus various forms of noise for an ordinary $\pi $ pulse (denoted as OP), and a noise-robust pulse (denoted as NRP). (a) Infidelity as a function of fluctuations in amplitudes $\delta\omega /\omega $. The calculations are preformed based on a detuning noise of $\delta v/v=0.01$. (b) Infidelity as a function of fluctuations in detuning $\delta \upsilon /\upsilon $. The calculations are performed based on an amplitude noise of $\delta \omega /\omega =0.01$.
In conclusion, we have developed composite pulse sequences for high-fidelity excitation transfer in an AFC-SW based quantum memory. These composite pulses are designed to compensate for variations in various experimental parameters (e.g., intensity or amplitude, detuning or frequency shifts). We theoretically simulate the fidelity of the transfer process in realistic rare-earth ion-doped systems. In particular, our results show high fidelity beyond $99\%$ across a broad range of experimental parameters, and for different cases. The operational fidelity of these composite pulses are significantly larger compared to conventional $\pi $ pulses. The designed control pulses would be a highly accurate and robust tool for a wide variety of quantum memory applications (e.g., EIT,[33] DLCZ,[34] GEM[35] or microwave[36,37]), and would be particularly valuable under conditions involving significant experimental uncertainties.
References Quantum memoriesQuantum repeaters based on atomic ensembles and linear opticsProspective applications of optical quantum memoriesQuantum information transfer using photonsQuantum memories: emerging applications and recent advancesQuantum internet: A vision for the road aheadFunctional Quantum Nodes for Entanglement Distribution over Scalable Quantum NetworksAn elementary quantum network of single atoms in optical cavitiesObservation of entanglement between a quantum dot spin and a single photonQuantum-dot spin–photon entanglement via frequency downconversion to telecom wavelengthHeralded entanglement between solid-state qubits separated by three metresA solid-state light–matter interface at the single-photon levelOptically addressable nuclear spins in a solid with a six-hour coherence timeMapping multiple photonic qubits into and out of one solid-state atomic ensembleSpectral Multiplexing for Scalable Quantum Photonics using an Atomic Frequency Comb Quantum Memory and Feed-Forward ControlQuantum storage of photonic entanglement in a crystalBroadband waveguide quantum memory for entangled photonsQuantum Storage of Heralded Polarization Qubits in Birefringent and Anisotropically Absorbing MaterialsQuantum Storage of a Photonic Polarization Qubit in a SolidRealization of Reliable Solid-State Quantum Memory for Photonic Polarization QubitQuantum Storage of Three-Dimensional Orbital-Angular-Momentum Entanglement in a CrystalMultimode quantum memory based on atomic frequency combsSolid State Spin-Wave Quantum Memory for Time-Bin QubitsCoherent Spin Control at the Quantum Level in an Ensemble-Based Optical MemoryNMR population inversion using a composite pulseComposite pulses without phase distortionDesigning robust unitary gates: Application to concatenated composite pulsesUniversal quantum computer from a quantum magnetRobust Dynamical Decoupling for Quantum Computing and Quantum MemoryCorrection of Arbitrary Field Errors in Population Inversion of Quantum Systems by Universal Composite PulsesComposite pulses for robust universal control of singlet–triplet qubitsNoise-Resistant Control for a Spin Qubit ArrayStorage of Light in Atomic VaporLong-distance quantum communication with atomic ensembles and linear opticsEfficient quantum memory for lightQuantum Memory for Microwave Photons in an Inhomogeneously Broadened Spin EnsembleMeasurement of Spin Singlet-Triplet Qubit in Quantum Dots Using Superconducting Resonator *
[1] Simon C et al. 2010 Eur. Phys. J. D 58 1
[2] Sangouard N, Simon C, de Riedmatten H and Gisin N 2011 Rev. Mod. Phys. 83 33
[3] Bussieres F, Sangouard N, Afzelius M, de Riedmatten H, Simon C and Tittel W 2013 J. Mod. Opt. 60 1519
[4] Northup T E and Blatt R 2014 Nat. Photon. 8 356
[5] Heshami K, England D G, Humphreys P C, Bustard P J, Acosta V M, Nunn J and Sussman B J 2016 J. Mod. Opt. 63 2005
[6] Wehner1 S, Elkouss D and Hanson R 2018 Science 362 eaam9288
[7] Chow C W et al. 2007 Science 316 1316
[8] Ritter S et al. 2012 Nature 484 195
[9] Gao W B, Fallahi P, Togan E, Miguel-Sanchez J and Imamoglu A 2012 Nature 491 426
[10] De Greve K et al. 2012 Nature 491 421
[11] Bernien H et al. 2013 Nature 497 86
[12] de Riedmatten H, Afzelius M, Staudt M U, Simon C and Gisin N 2008 Nature 456 773
[13] Zhong M et al. 2015 Nature 517 177
[14] Usmani I, Afzelius M, de Riedmatten H and Gisin N 2010 Nat. Commun. 1 12
[15] Sinclair N et al. 2014 Phys. Rev. Lett. 113 053603
[16] Clausen C, Usmani I, Bussieres F, Sangouard N, Afzelius M, de Riedmatten H and Gisin N 2011 Nature 469 508
[17] Saglamyurek E, Sinclair N, Jin J, Slater J A, Oblak D, Bussieres F, George M, Ricken R, Sohler W and Tittel W 2011 Nature 469 512
[18] Clausen C, Bussieres F, Afzelius M and Gisin N 2012 Phys. Rev. Lett. 108 190503
[19] Gundogan M, Ledingham P M, Almasi A, Cristiani M and de Riedmatten H 2012 Phys. Rev. Lett. 108 190504
[20] Zhou Z Q, Lin W B, Yang M, Li C F and Guo G C 2012 Phys. Rev. Lett. 108 190505
[21] Zhou Z Q, Hua Y L, Liu X, Chen G, Xu J S, Han Y J, Li C F and Guo G C 2015 Phys. Rev. Lett. 115 070502
[22] Afzelius M, Simon C, de Riedmatten H and Gisin N 2009 Phys. Rev. A 79 052329
[23] Gundogan M, Ledingham P M, Kutluer K, Mazzera M and de Riedmatten H 2015 Phys. Rev. Lett. 114 230501
[24] Jobez P et al. 2015 Phys. Rev. Lett. 114 230502
[25] Levitt M H and Freeman R 1979 J. Magn. Reson. 33 473
[26] Tycko R, Cho H M, Schneider E and Pines A 1985 J. Magn. Reson. 61 90
[27] Ichikawa T, Bando M, Kondo Y and Nakahara M 2011 Phys. Rev. A 84 062311
[28] Minar J, Sangouard N, Afzelius M, de Riedmatten H and Gisin N 2010 Phys. Rev. A 82 052309
[29] Souza A M, Alvarez G A and Suter D 2011 Phys. Rev. Lett. 106 240501
[30] Genov G T, Schraft D, Halfmann T and Vitanov N V 2014 Phys. Rev. Lett. 113 043001
[31] Wang X et al. 2012 Nat. Commun. 3 997
[32] Kestner J P, Wang X, Bishop L S, Barnes E and Das S S 2013 Phys. Rev. Lett. 110 140502
[33] Phillips D F, Fleischhauer A, Mair A, Walsworth R L and Lukin M D 2001 Phys. Rev. Lett. 86 783
[34] Duan L M, Lukin M D, Cirac J I and Zoller P 2001 Nature 414 413
[35] Hedges M P, Longdell J J, Li Y and Sellars M J 2010 Nature 465 1052
[36] Julsgaard B, Grezes C, Bertet P and Molmer K 2013 Phys. Rev. Lett. 110 250503
[37] Zhu X Y, Tu T, Guo A L, Zhou Z Q and Guo G C 2020 Chin. Phys. Lett. 37 020302