李蓬勃的个人主页空间 - 李 蓬勃 - 教师个人主页.pdf
PHYSICAL REVIEW APPLIED 4, 044003 (2015) Hybrid Quantum Device Based on NV Centers in Diamond Nanomechanical Resonators Plus Superconducting Waveguide Cavities Peng-Bo Li,1 Yong-Chun Liu,2 S.-Y. Gao,1 Ze-Liang Xiang,3 Peter Rabl,3 Yun-Feng Xiao,2 and Fu-Li Li1 1 Department of Applied Physics, Institute of Quantum Optics and Quantum Information, Xi’an Jiaotong University, Xi’an 710049, China 2 State Key Laboratory for Mesoscopic Physics, School of Physics, Peking University; Collaborative Innovation Center of Quantum Matter, Beijing 100871, China 3 Institute of Atomic and Subatomic Physics, TU Wien, Stadionallee 2, 1020 Wien, Austria (Received 25 August 2014; revised manuscript received 4 August 2015; published 8 October 2015) We propose and analyze a hybrid device by integrating a microscale diamond beam with a single built-in nitrogen-vacancy (NV) center spin to a superconducting coplanar waveguide (CPW) cavity. We find that under an ac electric field the quantized motion of the diamond beam can strongly couple to the single cavity photons via a dielectric interaction. Together with the strong spin-motion interaction via a large magneticfield gradient, it provides a hybrid quantum device where the diamond resonator can strongly couple both to the single microwave-cavity photons and to the single NV center spin. This enables coherent information transfer and effective coupling between the NV spin and the CPW cavity via mechanically dark polaritons. This hybrid spin-electromechanical device, with tunable couplings by external fields, offers a realistic platform for implementing quantum information with single NV spins, diamond mechanical resonators, and single microwave photons. DOI: 10.1103/PhysRevApplied.4.044003 I. INTRODUCTION Hybrid quantum architectures (for a review, see Ref. [1]) take the advantages and strengths of different components, involving degrees of freedom of completely different nature, which are promising for developing quantum technologies and discovering rich physics. Furthermore, the construction of hybrid quantum devices can benefit greatly from the progress achieved so far in the fields of atomic physics, quantum optics, condensed-matter physics, and nanoscience. A growing interest is emerging for exploring hybrid quantum architectures that could find applications in implementing quantum technologies. Prominent examples include superconducting waveguide cavities or mechanical resonators coupled to cold atoms [2–7], polar molecules [8–10], and quantum dots [11–15], as well as other solid-state spin systems [16–36]. A key challenge in the field of hybrid quantum systems is the realization of a controlled interface between a superconducting circuit and a single solid-state spin qubit. It would allow the realization of long-lived quantum memories for superconducting qubits, without the big problem of inhomogeneous broadening encountered in spin ensembles [37]. However, the direct magnetic coupling between a single spin and a single microwave photon is typically only a few hertz [2], much smaller than the relevant decoherence rates. To overcome this problem, it has been proposed to enhance this coupling via an intermediate persistent current loop [21,22], but the decoherence rates of such small superconducting loops are still unknown and will to a large extent compensate the achievable increase in coupling strength. 2331-7019=15=4(4)=044003(13) In this work, we propose and analyze a practical design for a coherent quantum interface between a single spin qubit and a microwave resonator, which makes use of the recent advances in the fabrication of single-crystal diamond nanoresonators. Over the past years, mechanical resonators made out of diamond have received considerable attention for the study of fundamental physics, as well as for applications in quantum science and technology [25,38–44]. Diamond, in addition to possessing desirable mechanical and optical properties, can host defect color centers [45], whose highly coherent electronic spins are particularly useful for quantum-information processing at room temperature. Single-crystal diamond cantilevers with exceptional quality factors exceeding one million have been demonstrated in a recent experiment [38]. Moreover, hybrid quantum systems, consisting of single-crystal diamond cantilevers with embedded NV center spins, have been realized in recent experiments [39,40]. Such on-chip devices, enabling direct coupling between mechanical and spin degrees of freedom, can be used as key components for hybrid quantum systems. In the setup investigated in this work, a doubly clamped diamond microbeam with a single built-in NV center is placed in the near field of a coplanar waveguide (CPW) cavity (see Fig. 1). Compared to previously considered capacitive coupling schemes [46] usually employed in cavity electromechanics [47–49], we here consider the case where the coupling between the quantized motion of the diamond beam and the microwave-cavity photons results from the dielectric interaction—a fundamental 044003-1 © 2015 American Physical Society LI et al. PHYS. REV. APPLIED 4, 044003 (2015) mechanism that any polarizable body placed in an inhomogeneous electric field will experience a dielectric force. This dielectric coupling is employed in the classical regime for the on-chip actuation of thin mechanical beams [50] and is scalable to large arrays of diamond mechanical resonators. For the present purpose of a quantum interface, it is important that this dielectric coupling is fully controlled by a tunable driving field but opposed to the capacitive coupling [47–49]; it results in large photon-phonon interactions without driving the cavity mode itself. When combined with a magnetic spin-phonon interaction in the presence of a strong magnetic-field gradient [51], a tunable coupling for exchanging a quantum state between the single spin memory and the superconducting microwave cavity can be implemented. We show that under realistic conditions the resulting effective interactions between a single spin qubit and a single microwave photon can exceed 10 kHz, roughly 3 orders of magnitude stronger than the direct coupling. The state transfer can further be optimized by employing a resonant state transfer scheme via mechanically dark polaritons, and we discuss the expected transfer fidelities for different parameters. II. DESCRIPTION OF THE DEVICE A. The setup B. Dielectric interaction We first describe the coupling between the diamond beam and the localized electric field of the CPW cavity. 6 x (µm) As shown in Figs. 1(a) and 1(b), a doubly clamped microscale diamond beam embedding a single spin-1 NV center is positioned along the z axis at a distance x0 from the gap of the CPW cavity surface. The configuration shown in Figs. 1(c) and 1(d) is somewhat equivalent to that of Figs. 1(a) and 1(b). Through small tip electrodes [52,53], ~ p ðtÞ (with frequency ωp and a strong ac electric field E amplitude Ep ) is applied to the diamond beam, which induces a large macroscopic electric-dipole moment. The diamond beam of length l has a circular cross section of radius r (r ≪ l). As the beam vibrates, x0 changes by the beam’s effective transverse displacement, and restricted to the lowest vibrational mode it can be modeled as a harmonic oscillator with frequency ωm and bosonic mode operator b̂. For a thin beam, the resonance frequency ωm can be calculated by the Euler-Bernoulli theory, i.e., ωm ¼ pffiffiffiffiffiffiffiffiffiffiffiffiffiffi k20 EI=ρA [54], where k0 ¼ 4.73=l is the wave number of the fundamental mode, E is the Young’s modulus, I is the moment of inertia, ρ is the density of diamond, and A is the cross-section area. The CPW cavity consisting of a central conductor stripe plus two ground planes is fabricated on a dielectric substrate which supports quasi-TEM microwave fields strongly confined near the gaps between the conductor and the ground planes (Fig. 2). For a CPW cavity of stripline length L and electrode distance d (with effective cavity volume V c ∼ πd2 L), the single-mode electric-field operator of the CPW cavity (with frequency ωc ) can be written as ~ˆ c ð~r; tÞ ¼ E 0 e~tr ðx; yÞðâe−iωc t þ ↠eiωc t Þ cosðπz=LÞ, [54] E pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi where E 0 ¼ ℏωc =ϵ0 V c is the field amplitude, e~tr ðx; yÞ the dimensionless transverse mode function (see the simulations in Fig. 2), and â the destruction operator for the microwave-cavity photons. Then, the free Hamiltonian for the CPW cavity field reads Ĥ c ¼ ℏωc ð↠â þ 12Þ. (a) 4 2 0 −20 1 −10 0 y (µm) x = 0.5 µm x = 5 µm x = 15 µm 0 −20 FIG. 1. Schematic of a hybrid device for coupling a doubly clamped diamond microbeam with a single, well-controlled embedded NV center spin to the CPW cavity field. (a) Top view of the configuration of the diamond beam positioned near z0 ∼ L above the gap of the CPW cavity. (b) Side view of the device shown in (a). (c) Top view of a different device where the beam couples to the electric field between the central conductor stripe of the CPW cavity and an external electrode. (d) Side view of the device shown in (c). −10 (c) 0 y (µm) 10 20 y = 0.2 µm y = 5 µm 0.5 0 0 20 (b) 0.5 1 10 5 10 x (µm) 15 20 FIG. 2. (a) Distribution of the transverse electric field (per photon) above the CPW cavity. (b),(c) The electric field (per photon) at various positions above the surface of the CPW cavity. 044003-2 HYBRID QUANTUM DEVICE BASED ON NV CENTERS … Generally, the electrostatic interaction between a dielectric object and electric fields can be described by Ĥpol ¼ R ~ rÞ · Eð~ ~ rÞd~r, where Pð~ ~ rÞ is the polarization − 12 V Pð~ ~ rÞ. In the linear response induced by the electric field Eð~ regime, the polarization field responds linearly to the ~ ¼ αE, ~ where α is the polarelectric field, such that P izability tensor, which depends on the symmetry of the diamond beam and its relative orientation with the cavity. The diamond beam is assumed to be placed in parallel with the gap of the cavity, i.e., along the z direction as shown in Fig. 1(a). The total electric field affected by the diamond beam is a transverse field, which can be written ~ ⊥ ð~r; tÞ ¼ E ~ p ðtÞ þ E ~ˆ c ð~r; tÞ. as E Considering the case where the dimension of the diamond beam is much smaller than the wavelength of the field, its dielectric response is well approximated by a point dipole, and the components of the polarizability tensor can be approximated by those induced by a uniform electric field. The polarization responds linearly to the ~ rÞ ¼ α⊥ E ~ ⊥ þ αz E ~ z, electric field, approximating to Pð~ α⊥ ¼ ϵ0 ðϵr − 1Þ=½1 þ N ⊥ ðϵr − 1Þ, αz ¼ ϵ0 ðϵr − 1Þ=½1 þ N z ðϵr − 1Þ for a dielectric microbeam [55], where ϵ0 is the free space permittivity, ϵr ¼ ϵ=ϵ0 , N z ¼ 1−e2 ðlnð1 þ eÞ=ð1 − eÞ − 2eÞ, N ⊥ ¼ 12 ð1 − N z Þ, and 2e3 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi e ¼ 1 − r2 =l2 . Here the transversal and longitudinal directions are defined by the coordinate axes as shown in Fig. 1(a). As the beam vibrates, the cavity electric field affected by the beam is modulated by the vibration. We first derive the interaction between the CPW cavity field and the diamond beam depicted in Figs. 1(a) and 1(b). Expanding the cavityfield operator around the position of the beam up to first order in the transverse displacement operator q̂x, we obtain ~ p · ∂ xE ~ c ð~r; tÞq̂x , with V the volume of [54] Ĥpol ¼ −Vα⊥ E the diamond beam. By assuming Δ ¼ ωp − ωc ≪ ωp ; ωc and neglecting all rapidly oscillating terms, we obtain the linear photon-phonon coupling [54] results in a dramatic enhancement of the field per photon E 0 and a large field gradient. Now we consider the coupling between the CPW cavity and the diamond beam as depicted in Figs. 1(c) and 1(d). In this device, the beam couples to the electric field between the central conductor stripe of the CPW cavity and an external electrode. Therefore, we can use the voltage distributionpofffiffiffiffiffiffiffiffiffiffiffiffiffiffi the cavity uðzÞ ¼ u0 cosðπz=LÞð↠þ âÞ, with u0 ¼ ℏωc =C and C the total capacitance of the cavity. The electric field from the electrodes is approximated as jEc j ∼ u0 ζ=h, and j∂Ec =∂xj ∼ u0 ζ=h2 , where h is the height of the beam above the substrate and ζ is a dimensionless constant of order unity set by the electrode geometry [6]. Then the coupling between the beam and the cavity in the interaction picture can be approximated as ∂Ec q̂ ð↠e−iΔt þ âeiΔt Þ ∂x x ¼ ℏgðâeiΔt þ ↠e−iΔt Þðb̂e−iωm t þ b̂† eiωm t Þ Ĥpol ¼ −Vα⊥ Ep ð3Þ with uζ g ¼ −Vα⊥ Ep 02 h sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 : 2mℏωm ð4Þ The linear coupling between the CPW cavity field and the vibrational mode of the beam is analogous to the effect of radiation pressure on a moving mirror of an optical cavity [56,57] or the capacitive coupling between the motion of a mechanical resonator and an electrical circuit [47–49]. However, different from the previous studies in optomechanics with linear optomechanical coupling, here the coupling is at the single-photon and -phonon level. If the detuning is chosen as Δ ¼ ωp − ωc ∼ ωm ≫ g, then under the rotating-wave approximation we can obtain the beam-splitter Hamiltonian Ĥpol ¼ ℏgðâeiΔt þ ↠e−iΔt Þðb̂e−iωm t þ b̂† eiωm t Þ: ð1Þ H1 ¼ ℏΔ↠â þ ℏωm b̂† b̂ þ ℏg↠b̂ þ ℏgâb̂† : Here sffiffiffiffiffiffiffiffiffiffiffiffiffi 1 ℏ ~ p · ½∂ x e~tr ðx; yÞðx ;y Þ g ¼ − Vα⊥ E 0 E ; 0 0 ℏ 2mωm PHYS. REV. APPLIED 4, 044003 (2015) ð5Þ C. Cooling of the vibration mode ð2Þ where m is the effective mass of the mechanical resonator. This coupling strength is proportional to the classic electric-field amplitude and the cavity-field gradient along the transverse direction. It can be greatly enhanced, since the classical electric-field amplitude Ep can be very large for strong enough fields, and the high concentration of cavity-field energy near the surface of the CPW cavity We now discuss ground-state cooling of the vibration mode of the diamond beam through the cavity-assisted sideband cooling approach in the resolved sideband regime [58–60], where ωm > κ, with κ being the cavity decay rate. For the mechanical mode of frequency ωm =2π ∼ 320 kHz, assuming an environmental temperature T ∼ 20 mK, the thermal phonon number is about nth ∼ 103 . Thus, additional cooling of the vibration mode is needed. Taking the dissipations of mechanical motion and cavity photons into consideration, the electromechanical system is governed by the quantum master equation 044003-3 LI et al. PHYS. REV. APPLIED 4, 044003 (2015) dρ̂ i ¼ − ½H1 ; ρ̂ þ κD½âρ̂ dt ℏ þ nth γ m D½b̂† þ ðnth þ 1Þγ m D½b̂; ð6Þ with κ the cavity photon loss rate, γ m the mechanical damping rate of the beam due to clamping, nth ¼ ðeℏωm =kB T − 1Þ−1 the thermal phonon number at the environment temperature T, and D½ôρ̂ ¼ ô ρ̂ ô† − 12 ô† ô ρ̂ − 1 † 2 ρ̂ô ô for a given operator ô. We focus on the resolved sideband regime ωm ≫ κ and set Δ ¼ ωm , in which the beam-splitter interaction is on resonance. To realize cooling, the cooperativity 4g2 =γκ ≫ 1 and the dynamical stability condition from the Routh-Hurwitz criterion 2g < ωm should be satisfied. To calculate the mean phonon number nm ¼ hb̂† b̂i, we need to determine the mean values of all the secondorder moments h↠âi;hb̂† b̂i;h↠b̂i;hâ b̂i;hâ2 i;hb̂2 i. Starting from the master equation, we obtain a set of differential equations for the mean values of the second-order moments as d † hâ âi ¼ −ighð↠− âÞðb̂† þ b̂Þi − κh↠âi; dt d † hb̂ b̂i ¼ −ighð↠þ âÞðb̂† − b̂Þi − γ m hb̂† b̂i þ γ m nm ; dt d † γ þκ † hâ b̂i ¼ − m hâ b̂i − ig½hð↠Þ2 i − hðb̂† Þ2 i dt 2 − hb̂† b̂i þ h↠âi; d γm þ κ hâ b̂i ¼ −iðωm þ ΔÞ − hâ b̂i − ig½1 þ hb̂2 i dt 2 obtain the final phonon number nf ∼ 0.3 for the vibration mode, which is well in the quantum ground state. This result is also verified based on numerically solving the quantum master equation (6) (Fig. 3), from which we find that the mechanical mode will enter the quantum ground state at the time t ∼ 100 μs [Fig. 3(a)]. In the experiment, the diamond beam will need to be sufficiently isolated from other mechanical modes that it can be cooled to its ground state. We now consider the influence of other mechanical modes on the cooling process of the fundamental mode. From the EulerBernoulli theory, we know that the first two vibrating modes have the resonance frequencies [54] sffiffiffiffiffiffi EI ; ρA sffiffiffiffiffiffi 7.852 EI ω1 ¼ 2 : ρA l 4.732 ωm ¼ 2 l ð10Þ Therefore, if the detuning between these modes is sufficiently large, i.e., δ ¼ ω1 − ωm ≫ g, then we can safely ignore the coupling between the cavity mode and other mechanical modes. One can find that ω1 ∼ 3ωm , and δ ∼ 2ωm ∼ 4 MHz, which is much larger than the vibration-photon coupling strength. Thus, we need only to consider the coupling between the fundamental mechanical mode and the cavity photons. Provided that ω1 − ωm ≫ g, the diamond resonator will be sufficiently isolated from other mechanical modes. Cooling of the mechanical resonators to the quantum ground state has been realized in a þ hâ2 i þ h↠âi þ hb̂† b̂i; d 2 hâ i ¼ −2ighðb̂† þ b̂Þâi − ½κ þ 2iΔhâ2 i; dt d 2 hb̂ i ¼ −2ighð↠þ âÞb̂i − ½γ m þ 2iωm hb̂2 i: dt ð7Þ The steady-state solution of these equations then yields an analytical formula for the final occupancy of the mechanical mode. In the weak coupling regime, i.e., g ≪ κ, the final phonon number is nf ≃ γ m nth κ2 þ ; Γ þ γ m 16ω2m Γ ¼ 4g2 =κ; ð8Þ while in the strong coupling regime, i.e., κ ≪ g ≪ ωm , the final phonon number is nf ≃ γ m nth g2 þ : κ þ γ m 2½ω2m − 4g2 ð9Þ In this work, we focus on the strong coupling regime. With the given experimental parameters in the main text, we FIG. 3. (a) Time evolution of the mean phonon number from numerically solving the master equation (6) under realistic parameters. (b) The final phonon number in a logarithmic scale as a function of g=κ. The parameters are chosen as those in (a). 044003-4 HYBRID QUANTUM DEVICE BASED ON NV CENTERS … variety of experiments exploiting the cavity-assisted sideband cooling approach [61,62]. We believe that cooling the diamond resonators via the same approach, as described in this work, could be realized in the near future with the stateof-the-art technology. D. Spin-motion couplings We now turn to considering the coupling between the NV center spin and the mechanical motion. Quantum interface between the spin and the CPW cavity can be achieved through a Jaynes-Cummings (JC) spin-motion interaction under a large magnetic-field gradient [51,63–65]. NV centers have a S ¼ 1 ground state, with zero-field splitting D ¼ 2π × 2.87 GHz between the jms ¼ 1i and jms ¼ 0i states. We consider a single NV center embedded in the thin diamond beam, with its NV axis along the z axis. An external magnetic field composed of a homogeneous bias field and a local linear gradient is applied to ~ NV ¼ B0 e~z þ ð∂B=∂xÞx~ex . The homothe system, i.e., B geneous bias field is used to lift the degeneracy of states jms ¼ þ1i and jms ¼ −1i, while the local magnetic-field gradient is used to couple the mechanical motion of the beam to the S ¼ 1 spin of the NV center. The magnetic-field gradient can be achieved by a single-domain ferromagnet such as Co nanobars near the NV center [51,63–65]. In a realistic setup, an objective lens and a laser are needed for optically pumping the spin state of the NV center. The Hamiltonian for the spin-motion coupling system ~ NV , with then is ĤNV ¼ ℏDS2z þ ℏωm b̂† b̂ þ gNV μB S~ · B gNV ¼ 2 the Landé factor of the NV center, μB the Bohr magneton, and S~ the spin operator for the NV center. When D − gNV μB B0 =ℏ − ωm ¼ δ ≪ gNV μB B0 =ℏ, we can make the rotating-wave approximation to describe the nearresonance interaction between the NV spin and the vibration mode and neglect the far-out-of-resonance state jms ¼ þ1i. In this case, the JC Hamiltonian describing the spin-motion dynamics reads [54] 1 H2 ¼ ℏωþ σ̂ z þ ℏωm b̂† b̂ þ ℏλb̂σ̂ þ þ ℏλb̂† σ̂ − ; 2 ð11Þ PHYS. REV. APPLIED 4, 044003 (2015) 1 H ¼ ℏωþ σ̂ z þ ℏωm b̂† b̂ þ ℏΔ↠â 2 þ ℏg↠b̂ þ ℏgâb̂† þ ℏλb̂σ̂ þ þ ℏλb̂† σ̂ − : ð12Þ The first three terms describe the free Hamiltonian for the spin, the phonons, and the photons, while the last four terms describe the coherent coupling of the phonons to both the spin and the photons. In a realistic experimental situation, we need to consider cavity photon loss, mechanical dissipation, and spin dephasing. The full dynamics of our system that takes these incoherent processes into account is described by the master equation dρ̂ðtÞ i ¼ − ½H; ρ̂ þ κD½âρ̂ þ γ s D½σ̂ z ρ̂ dt ℏ þ nth γ m D½b̂† þ ðnth þ 1Þγ m D½b̂ ð13Þ with γ s the single spin dephasing rate of the NV center. To enter the strong coupling regime requires that fg; λg > fnth γ m ; κ; γ s g. Let us consider the experimental feasibility in the configuration as shown in Fig. 1 and the appropriate parameters to achieve strong coupling. (i) For a CPW cavity with strip-line length L ∼ 1 cm, electrode distance d ∼ 5 μm, and effective dielectric constant ϵeff ∼ 6, the mode frequency for the CPW cavity is ωc ¼ pffiffiffiffiffiffiffi πc=L ϵeff ∼ 2π × 6 GHz. With the above parameters for the CPW cavity, the electric-field amplitude of a single photon is E 0 ∼ 0.76 V=m. If the beam is positioned at a distance of about 1 μm above the CPW surface gap, then we can estimate the mode function as j~etr ðx0 ; y0 Þj ∼e−0.2 and ½∂ x e~tr ðx; yÞðx0 ;y0 Þ ∼ ð2 μmÞ−1 , respectively [see Fig. 2(c)]. (ii) We consider a diamond microbeam of crosssectional radius r ∼ 100 nm under a strong ac electric field with Ep ∼ 10 V=μm and a gradient magnetic field with ∂ x B ∼ 107 T=m [54]. Figure 4 shows the calculated coupling strengths as a function of the length of the diamond beam with the given parameters. We find that the optimal length of the beam under the given parameters is about where ωþ ¼ D − gNV μB B0 =ℏ, σ̂ z ¼ j − 1ih−1j −pj0ih0j, ffiffiffi σ̂ þ ¼ j −p 1ih0j, σ̂ ¼ j0ih−1j, and λ ¼ ðg μ = 2ℏÞ× − NV B ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ð∂B=∂xÞ ℏ=2mωm . Using dressed state qubits [51] would be an equivalent alternative approach for implementing this model, which however would need microwave driving of the NV spin states. III. REALISTIC CONSIDERATIONS AND EXPERIMENTAL PARAMETERS Putting everything together, the total Hamiltonian describing the spin-mechanics-cavity hybrid tripartite system is FIG. 4. Photon-motion coupling strength g and spin-motion coupling strength λ as a function of the length of the diamond beam. The relevant parameters are r ∼ 100 nm, E 0 ∼ 0.76 V=m, ½∂ x e~tr ðx;yÞðx0 ;y0 Þ ∼ð2 μmÞ−1 , Ep ∼10 V=μm, and ∂ x B∼107 T=m. 044003-5 LI et al. PHYS. REV. APPLIED 4, 044003 (2015) (a) (b) (c) (d) T 1 can be several minutes at low temperatures. The effect of other centers coupled to the mechanical motion of the beam can be neglected [54]. Figure 5 shows the numerical simulations of quantum dynamics of the spin-mechanics-cavity system through solving the master equation (6). We find that, with the given parameters, the coherent interactions can dominate the decoherence processes in the hybrid setup, which enables the strong coupling regime to be entered. In this regime, the mechanical motion becomes strongly coupled to the spin and the cavity photons in direct analogy to strong coupling of cavity QED. FIG. 5. (a) Vacuum Rabi oscillations of the hybrid system where the mechanical resonator couples to the spin and the cavity without dissipations. The initial state of the beam is a thermal state with nm ∼ 0.3, while the spin is initially in the spin-up state and the cavity in the ground state. (b) Photon distribution in the cavity after a half period of Rabi oscillations. (c),(d) The same as (a),(b) but with dissipations for the spin, the mechanical resonator, and the cavity. The relevant parameters are g=2π ∼ 16 kHz, λ=2π ∼ 16 kHz, κ=2π ∼ 6 kHz, nth ∼ 1000, γ m =2π ∼ 3.2 Hz, and γ s =2π ∼ 2 kHz. l ∼ 80 μm. The vibration frequency for the fundamental mode is ωm =2π ∼ 320 kHz. Then, the photon-motion and spin-motion coupling strengths can reach g=2π ∼ 16 kHz and λ=2π ∼ 16 kHz, respectively. We next turn to damping of the mechanical motion. For a diamond beam with frequency ωm and quality factor Q, the mechanical damping rate is γ m ¼ ωm =Q. The recent demonstration of high-quality single-crystal diamond beams or cantilevers with embedded NV centers leads to a mechanical quality factor exceeding 105 [38–44]. For our doubly clamped diamond beam with vibration frequency ωm =2π ∼ 320 kHz, the damping rate is about γ m =2π ∼ 3.2 Hz and to ensure the strong coupling condition thus requires that nth < 5000. For the mechanical mode of frequency ωm =2π ∼ 320 kHz, assuming an environmental temperature T ∼ 20 mK in a dilution refrigerator, the thermal phonon number is about nth ∼ 103 . To keep the mechanical motion in the quantum ground state thus needs additional cooling. Finally, we consider the cavity photon loss and the dephasing of the NV spin. For a realistic value of the quality factor for the CPW cavity Q ∼ 106, the photon decay rate is about κ=2π ∼ 6 kHz. In the experiment, superconducting cavities are able to maintain high Q even at applied in-plane magnetic fields > 200 mT [18–20]. Therefore, we can safely ignore the effect of ultrastrong nanomagnets in close proximity to superconducting CPW cavities. When it comes to the NV center, the dephasing time T 2 induced by the fluctuations in the nuclear spin bath can be increased to several milliseconds in ultrapure diamond [66]. We can ignore single spin relaxation as IV. APPLICATIONS A. Quantum state transfer via dark polaritons In general, the spin-electromechanical hybrid tripartite system modeled by the Hamiltonian (12) can find use in many aspects of quantum science and technology. In the following, we propose to realize coherent information transfer between the single spin and the microwave cavity, using mechanically dark polaritons. This method is very efficient and robust against mechanical noise compared to the direct-transfer process, since, during the state conversion process, the mechanical mode is decoupled from the dark polaritons composed of the spin and the cavity modes. We proceed by assuming that ωþ ∼ Δ ∼ ωm ¼ Ω. Then, it can readily be verified that the Hamiltonian (12) can take a compact form in terms of spin-photon polariton and spin-photon-phonon polaron operators H ¼ ℏΩP †d P d þ ℏΩþ P †þ P þ þ ℏΩ− P †− P − ; ð14Þ with P d ¼ cos θσ̂ − − sin θâ; tan θ ¼ λ=g being the polariton operator, describing quasiparticles formed by combinations pffiffiffi of spin and photon excitations, and P ¼ ð1= 2ÞðP b b̂Þ, P b ¼ cos θâ þ sin θσ̂ − , describing polarons formed by combinations of polariton and phonon excitations. pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiThe frequencies of the polarons are Ω ¼ Ω g2 þ λ2 . We refer to P d as the mechanically dark polariton operator, due to the fact that it is decoupled from the mechanical mode, independently of the couplings. The mechanically dark polaritons can find use in quantum state conversion between the spin and cavity photons. This task can be accomplished through using an adiabatic passage approach, similar to the well-known STIRAP scheme [67–70]. By simply making gðtÞ initially bigger than λ, but keeping λ finite, one can modulate the coupling strengths gðtÞ slowly to ensure that the system adiabatically follows the dark polaritons. This will adiabatically rotate the mixing angle θ from 0 to π=2, leading to a complete and reversible transfer of the spin state to the photonic state; i.e., the dark polariton operator adiabatically evolves from being σ̂ − at t ¼ 0 to −â at the end of the protocol at a time t ¼ tf . 044003-6 HYBRID QUANTUM DEVICE BASED ON NV CENTERS … In Fig. 6, we display the numerical results for the state conversion process using dark polaritons through solving the p master equation (6). The initial spin state is chosen as ffiffiffi ð1= 2Þðj0i þ j − 1iÞ, while the cavity state is initially in the ground state and the mechanical mode is initially in the thermal state with nm ¼ 0.1. At the end p offfiffiffithe transfer process, the cavity state is steered into ð1= 2Þðj0i þ j1iÞ with a fidelity above 90%. This fidelity could further be improved with optimized pulses for gðtÞ and detunings. In the simulations, we choose the time dependence for the coupling g as exponential, but this pulse form is not necessarily required. The polarization coupling strength ~ p ðtÞ, which g is proportional to the classical electric field E can be tailored to give the desired time-dependent form. The spin state thus can be transferred from the NV center to the photonic state in the CPW cavity using the mechanical motion but without actually populating it. Therefore, this approach offers a distinct feature that the state transfer process exploiting dark polaritons is highly immune to mechanical noises. B. Effective strong coupling between the spin and the CPW cavity Alternatively, one can detune the mechanical mode and couple the spin to the microwave field via virtual motional excitations. In general, the direct coupling between a single NV spin and the microwave-cavity field is inherently rather weak. Here we propose to reach the effective strong coupling regime with this hybrid spin-electromechanical system. The induced effective strong coupling offers great potential for single-photon manipulation in the microwave frequency domain with this hybrid architecture. FIG. 6. (a) Adiabatic population transfer without dissipations. (b) The same as (a) but with dissipations for the spin, the mechanical resonator, and the cavity. Relevant parameters are 2 g ¼ g0 e−t =4 , g0 ¼ 1.8λ, κ ¼ 0.1λ; γ m ¼ 0.0001λ, nth ¼ 1000, and γ s ¼ 0.1λ. PHYS. REV. APPLIED 4, 044003 (2015) In the frame rotating at the spin’s resonance frequency ωþ , the Hamiltonian of the hybrid system is given by H ¼ ℏΔ1 b̂† b̂ þ ℏΔ2 ↠â þ ℏg↠b̂ þ ℏgâb̂† þ ℏλb̂σ̂ þ þ ℏλb̂† σ̂ − ; ð15Þ where Δ1 ¼ ωm − ωþ and Δ2 ¼ Δ − ωþ . Taking the dissipations into consideration, the system is described by the quantum master equation dρ̂ðtÞ i ¼ − ½H; ρ̂ þ κD½âρ̂ þ γ s D½σ̂ z ρ̂ þ nth γ m D½b̂† dt ℏ þ ðnth þ 1Þγ m D½b̂: ð16Þ By adiabatically eliminating the mechanical mode b̂ for large detunings, Δ1 ; Δ2 ≫ g; λ, virtual excitations of the mechanical mode result in an effective interaction between the NV spin and the CPW cavity mode â, with the effective Hamiltonian 1 Heff ¼ ℏðΔ2 − β2 Δ1 Þ↠â − α2 Δ1 σ̂ z þ ℏgeff ↠σ̂ − 2 þ ℏgeff âσ̂ þ ; ð17Þ where the parameters α ¼ λ=Δ1 , β ¼ g=Δ1 , and the effective spin-photon coupling strength geff ¼ αg. Then the reduced density matrix for the spin-cavity system will satisfy the effective master equation dϱ̂ðtÞ i ¼ − ½Heff ; ϱ̂ þ κ 1eff D½âϱ̂ þ γ s D½σ̂ z ρ̂ þ γ 1eff D½σ̂ − ϱ̂ dt ℏ þ κ 2eff D½â† ϱ̂ þ γ 2eff D½σ̂ þ ϱ̂: ð18Þ The effective decay rates of the cavity mode and the NV spin are described by κ 1eff ¼ κ þ β2 ðnth þ 1Þγ m, κ2eff ¼ β2 nth γ m and γ 1eff ¼ α2 ðnth þ 1Þγ m , γ 2eff ¼ α2 nth γ m , respectively. It can be easily found that the effective coupling strength geff depends linearly on α, while the effective decay rates κieff and γ ieff are quadratic functions of the parameters α and β. Therefore, the spin-cavity coupled system can be steered into the strong coupling regime if fα; βg ≪ 1 and fκ; γ s g < geff ; i.e., the effective coupling strength can exceed the decay rates geff > γ s ; κieff ; γ ieff , i ¼ 1; 2. We now consider the relevant parameters to reach the effective strong coupling. Taking Δ1 ≃ Δ2 ¼ Δ; g ≃ λ; Δ ¼ 10g, one can get geff ¼ 0.1g, κ 1eff ≃ κ, κ2eff ≃ 0, γ 1eff ≃ γ 2eff ≃ 0. Then the effective master equation describing the spin-cavity system reads 044003-7 dϱ̂ðtÞ i ¼ − ½Heff ; ϱ̂ þ κD½âϱ̂ þ γ s D½σ̂ z ρ̂: dt ℏ ð19Þ LI et al. PHYS. REV. APPLIED 4, 044003 (2015) So the effective strong coupling regime requires geff > κ; γ s , which means g; λ ≥ 10κ; 10γ s . We now estimate the coupling strengths in the configuration shown in Figs. 1(c) and 1(d). For the CPW cavity considered in this work, u0 ∼ 5 μV. We consider a diamond microbeam of cross-sectional radius r ∼ 50 nm and length approximately 1 μm. Taking ζ ∼ 0.4 and h ∼ 100 nm, then we have g=2π ∼ 60 kHz and λ=2π ∼ 40 kHz. Then if we take Δ=2π ∼ 300 kHz, we can obtain geff =2π ∼ 10 kHz. This effective coupling strength can exceed the decay rates of a CPW cavity with a quality factor Q > 106. Magnetic coupling between a spin ensemble and a superconducting cavity has been reported recently [18–20], but the coupling between a single NV spin and the microwave-cavity field is inherently rather weak. Here we have proposed an efficient method to reach the effective strong coupling regime with this hybrid spin-electromechanical system. The resulting effective coupling strength can be significantly enhanced by approximately 3 orders of magnitude. Related schemes have been investigated before for interfacing a single spin to a transmission line cavity [22,46]. Different fundamentally from these proposals, here we specifically exploit the dielectric interaction through an ac electric field and state transfer schemes via mechanically dark polaritons. Both techniques are, in particular, useful to realize such interactions with NV centers in diamond beams, a system which is currently very actively explored in this context [38–44]. Furthermore, our proposed device just requires placing the diamond beam above the CPW surface, with no need to integrate the diamond resonator into the tiny coupling capacitor. This design is thus much easier to implement in practice and possesses the advantage of scalability, particularly for much bigger diamond microbeams. ACKNOWLEDGMENTS V. CONCLUSIONS We present a spin-mechanics-cavity hybrid device where a vibrating diamond beam with implanted single NV spins is coupled to a superconducting CPW cavity. We show that, under an ac electric field, the diamond beam can strongly couple to the CPW cavity through dielectric interaction. Together with the strong spin-motion interaction via a large magnetic-field gradient, it provides a hybrid quantum device where the diamond resonator can strongly couple both to the single microwave-cavity photons and to the single NV center spin. The distinct feature of this device is that it is on chip and scalable to large arrays of mechanical resonators coupled to the same CPW cavity. As for applications, we propose to use this hybrid setup to implement quantum-information transfer between the NV spin and the CPW cavity via mechanically dark polaritons. This hybrid spin-electromechanical device can offer a realistic platform for implementing quantum information with single NV spins, mechanical resonators, and single microwave photons. This work is supported by the NSFC under Grants No. 11474227 and No. 11534008 and the Fundamental Research Funds for the Central Universities. Part of the simulations are coded in PYTHON using the QUTIP library [71]. Work at Vienna is supported by the WWTF, the Austrian Science Fund (FWF) through SFB FOQUS, and the START Grant No. Y 591-N16 and the European Commission through Marie Sklodowska-Curie Grant No. IF 657788. APPENDIX A: DIELECTRIC COUPLINGS 1. Fundamental vibration mode of the diamond beam For a thin beam, the Euler-Bernoulli elastic theory is valid [72,73]. We consider a doubly clamped diamond beam with dimensions l ≫ r. The equation for the lateral vibration of a thin beam is ρA ∂2 ∂4 ϕðz; tÞ þ EI ϕðz; tÞ ¼ 0; ∂t2 ∂z4 ðA1Þ where ϕðz; tÞ is the lateral displacement in the x direction, A is the beam cross section, and I is the moment of inertia, I ¼ πr4 =8 for a cylindrical beam. The solutions to this equation are ϕðz; tÞ ¼ uðzÞe−iωt , with the mode function uðzÞ ¼ C1 ðcoskz − coshkzÞ þ C2 ðsinkz − sinhkzÞ; ðA2Þ which satisfies the boundary conditions uð0Þ ¼ uðlÞ ¼ 0; u0 ð0Þ ¼ u0 ðlÞ ¼ 0 for a doubly clamped beam. The frequency equation is given by cos kl cosh kl ¼ 0: ðA3Þ The first five nontrivial consecutive roots of this equation are given below: k0 l k1 l k2 l k3 l k4 l 4.730 7.853 10.996 14.137 17.279 and the corresponding eigenfrequencies are sffiffiffiffiffiffi EI : ωn ¼ k2n ρA ðA4Þ Therefore, the fundamental mode pffiffiffiffiffiffiffiffiffiffiffiffiffi ffi has the vibration fre2 2 quency ωm ¼ ð4.73 =l Þ EI=ρA. 2. CPW cavity-field operator For a CPW cavity as shown in Fig. 1 of the main text, the coplanar waveguide problem can be reduced to a rectangular waveguide problem by inserting magnetic walls at 044003-8 HYBRID QUANTUM DEVICE BASED ON NV CENTERS … y ¼ 0 and y ¼ b and electric walls at z ¼ 0 and z ¼ L [74,75]. If the central conductor is interrupted by two capacitors or gaps with a strip-line distance L, the cavity modes will be standing waves in the axial direction. Then, the classical electric-field components are given by [74,75] X 1 sin nπδ nπ δ̄ nπy −γn x 2 Ex ¼ − E0 cos e sin Fn nπ 2 b 2 δ n mπz × cos ; L X sin nπδ nπ δ̄ nπy 2 −γ n x Ey ¼ E 0 nπ sin sin e 2 b 2 δ n mπz × cos ; L Ez ¼ 0; ðA5Þ ðA6Þ ~ˆ r; tÞ ¼ E 0 e~tr ðx; yÞðâe−iωc t þ ↠eiωc t Þ cos πz ; Eð~ L 3. Cavity-resonator couplings For a thin beam with a circular cross section, it lacks a closed-form expression for the polarizability tenser. However, it has been shown that the analytical expression for the polarizability of a spheroid can be very close to that of a cylinder of the same permittivity ϵ and aspect ratio [55]. In the following, we employ the polarizability tenser of a thin prolate spheroid instead. Considering the case where the dimension of the diamond beam is much smaller than the wavelength of the electric field, its dielectric response is well approximated by a point dipole: ~ ⊥ ð~rÞδð~r − ~r0 Þ: ~ ð~r0 Þ ¼ Vα⊥ E p nπ δ̄ nπy −γn x cos e e~x sin e~tr ðx;yÞ ¼ − Fn nπ 2 b 2 δ n Xsin nπδ nπ δ̄ nπy 2 −γ n x ~ sin e ey þ ðA9Þ sin nπ 2 b 2 δ n X 1 sin nπδ 2 and e~x and e~y are the unit vectors for the x and y axes, respectively. Then the electric-field operator in the position of the diamond beam is ðA10Þ ðA11Þ In this case, the Hamiltonian describing the electrostatic interaction between the microbeam and the electric field is 1 ~ ⊥ ð~rdm ; tÞj2 : Ĥpol ¼ − Vα⊥ jE 2 ðA12Þ As the beam vibrates, the cavity electric field affected by the beam will be modulated by the vibration. Expanding the cavity-field operator around the position of the beam up to first order in the transverse displacement operator q̂x and neglecting rapidly oscillating and other higher-order terms, the Hamiltonian describing the coupled system reads ~ c ð~r0 ; tÞ ~p · E Ĥ pol ¼ −Vα⊥ E ~ p · ∂ xE ~ c ð~r; tÞq̂x : − Vα⊥ E ðA8Þ where the transverse mode function is ~ˆ rdm ; tÞ ¼ E 0 e~tr ðx0 ; y0 Þðâe−iωc t þ ↠eiωc t Þ; Eð~ from numerical simulations we find that the mode function can be approximated as j~etr ðx0 ; y0 Þj ∼ e−1 , and ½∂ x e~tr ðx; yÞðx0 ;y0 Þ ∼ γ ∼ ð5 μmÞ−1 . ðA7Þ pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi where δ ¼ d=b, δ̄∼δ, Fn ¼ðbγ n =nπÞ¼ 1þð2bv=nλ0 Þ2 , pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi v ¼ ðλ0 =λc Þ2 − 1, and λ0 is the free space wavelength for the mode frequency ωc. The cavity wavelength λc is related to the free space wavelength λ0 with the expression pffiffiffiffiffiffiffi λc ¼ λ0 = ϵeff , where ϵeff is the effective relative dielectric constant. The diamond beam positioned a few micrometers above the gap will experience the very strong localized electric field of the CPW cavity. In this work, we take the halfwavelength mode with m ¼ 1, in which case the cavity wavelength is λc ¼ 2L. Following the standard procedure for quantizing the electromagnetic fields, we obtain the quantized form of the single-mode electric-field operator for the CPW cavity: PHYS. REV. APPLIED 4, 044003 (2015) ðA13Þ The first term corresponds to the driving of the cavity mode by the electric dipole, while the second term describes the optomechanical coupling between the vibration and the cavity mode. In order to eliminate the first term, we can simply drive the cavity with a second field with an appropriately chosen amplitude and which is π out of ~ p . Then, the coupling between this phase with respect to E dipole and the cavity can cancel the first term in Eq. (A13) as a result of destructive interference between these two coupling terms. In this case, the pure electromechanical coupling of the vibration to the cavity mode can be derived as ~ p · ½∂ x e~tr ðx; yÞðx ;y Þ ðeiωp t þ e−iωp t Þ Ĥpol ¼ −Vα⊥ E 0 E 0 0 where we assume that the beam is positioned at the maximum slope of the standing wave mode, i.e., cosðπz0 =LÞ ¼ 1. For the case that the beam is positioned at a distance of several micrometers above the cavity gap, 044003-9 × ðâe−iωc t þ ↠eiωc t Þq̂x ~ p · ½∂ x e~tr ðx; yÞðx ;y Þ ðâeiΔt þ ↠e−iΔt Þq̂x : ¼ −Vα⊥ E 0 E 0 0 ðA14Þ LI et al. PHYS. REV. APPLIED 4, 044003 (2015) After quantizing the vibration mode of the diamond beam, ffi † pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi i.e., q̂x ¼ ℏ=2mωm ðb̂ þ b̂Þ, we have sffiffiffiffiffiffiffiffiffiffiffiffiffi ℏ ~ p · ½∂ x e~tr ðx; yÞðx ;y Þ Ĥpol ¼ −Vα⊥ E 0 E 0 0 2mωm 2. Large magnetic-field gradient induced by a micromagnet × ðâeiΔt þ ↠e−iΔt Þðb̂e−iωm t þ b̂† eiωm t Þ ¼ ℏgðâe iΔt † −iΔt þ â e −iωm t Þðb̂e † iωm t þ b̂ e Þ: ðA15Þ APPENDIX B: SPIN-MOTION COUPLINGS 1. Detailed derivation of the spin-motion interaction Hamiltonian We consider a single NV center embedded in the microscale diamond beam, with its NV axis along the z direction. NV centers have a S ¼ 1 ground state, with zerofield splitting D ¼ 2π × 2.87 GHz between the jms ¼ 1i and jms ¼ 0i states. The Hamiltonian describing the NV ~ NV has the form center in an external magnetic field B ~ NV : ĤNV ¼ ℏDS2z þ gNV μB S~ · B Ĥ NV ¼ ℏDS2z þ gNV μB B0 Sz þ ΛSx ðb̂ þ b̂† Þ ðB2Þ with the spin-motion pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi coupling strength Λ ¼ gNV μB ð∂B=∂xÞ ℏ=2mωm . In the basis defined by the eigenstates of Sz , i.e., fjms i; ms ¼ 0; 1g, with Sz jms i ¼ ms jms i, we have ĤNV ¼ ðℏD þ gNV μB B0 Þj þ 1ihþ1j þ ðℏD − gNV μB B0 Þj − 1ih−1j þ ℏωm b̂† b̂ þ Λðb̂ þ b̂† Þ½hþ1jSx j0ij þ 1ih0j þ h−1jSx j0ij − 1ih0j þ H:c:: ðB3Þ When D − gNV μB B0 =ℏ − ωm ¼ δ ≪ gNV μB B0 =ℏ, we can make the rotating-wave approximation to describe the near-resonance interaction between the NV spin and the mechanical motion and neglect the far-out-of-resonance state jms ¼ þ1i. Then we can obtain the Hamiltonian describing the spin-motion dynamics 1 H2 ¼ ℏωþ σ̂ z þ ℏωm b̂† b̂ þ ℏλb̂σ̂ þ þ ℏλb̂† σ̂ − ; 2 In the main text, the spin-motion coupling between the NV center and the vibration mode needs a large magneticfield gradient along the x direction. This can be generated ~ oriented along the x by a micromagnet of magnetization M axis near the NV center. The micromagnet can be created via lithographic processes, which is a single magnetic domain whose magnetic moment is spontaneously oriented along its long axis due to the shape anisotropy. We consider Co nanobars with dimensions ðl; w; tÞ ¼ ð200; 50; 50Þ nm. Approximating the bar by a magnetic dipole, we have 4 ~ and d0 is the ~ ~ ¼ lwtM ð∂B=∂xÞ ¼ 3μ0 jmj=4πd 0 , where m distance between the tip and the NV spin. Taking M ¼ 1.5 × 106 A=m, d0 ¼ 60 nm, we obtain ð∂B=∂xÞ ∼ 107 T=m. A large field gradient of such a value is reported in magnetic resonance force microscopy experiments [76]. ðB1Þ We assume that the external magnetic field is composed of a homogeneous bias field and a linear gradient, ~ NV ¼ B0 e~z þ ∂ x Bx~ex . Then, the Hamiltonian (B1) i.e., B becomes þ ℏωm b̂† b̂; where ωþ ¼ D − gNV μB B0 =ℏ, σ̂ z ¼ j − 1ih−1j −pj0ih0j, ffiffiffi σ̂ þ ¼ j −p 1ih0j, σ̂ ¼ j0ih−1j, and λ ¼ ðg μ = 2ℏÞ× − NV B ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ð∂B=∂xÞ ℏ=2mωm . ðB4Þ APPENDIX C: EFFECTS OF OTHER DECOHERENCE PROCESSES In this Appendix, we consider some other decoherence processes that are less important and have been ignored in the main text. These decoherence processes include scattered photon-recoil heating by dipole radiation and spin relaxation due to strain-induced coupling to the nearby NV centers. 1. Scattered photon-recoil heating by dipole radiation We now consider the scattered photon-recoil heating of the vibration mode due to dipole radiation. Since the electric dipole moment induced by the strong ac electric field in the diamond beam is oscillating in time, it will radiate photons from the beam in the form of dipole radiation, which in turn causes decoherence to the mechanical motion due to momentum recoil kicks. We assume that each scattered photon contributes the maximum possible momentum kick of ℏk along the x axis, giving rise to a momentum diffusion process dhp2x i=dt ¼ Γsc ðℏkÞ2 , where Γsc is the photon scattering rate. The decay rate γ sc associated with scattered photon-recoil heating can be approximated as γ sc ≃ ðωr =ωm ÞΓsc , with ωr ¼ ℏk2 =ð2mÞ the recoil frequency. For dipole radiation, the photon scattering rate Γsc can be calculated by taking the power ~ and dividing by the energy radiated by the dipole strength p per photon, i.e., Γsc ¼ ðc2 Z0 k4 =12πℏωÞj~ pj2 . For the given parameters, the photon-recoil heating rate is about γ sc ∼ 10−5 Hz, which thus can be neglected. 044003-10 HYBRID QUANTUM DEVICE BASED ON NV CENTERS … 2. Spin decoherence due to strain-induced coupling between nearby NV centers Let us consider the effect of nearby NV centers on this NV center spin via strain-induced couplings. When the beam vibrates, it strains the diamond lattice, which in turn couples directly to the spin-triplet states in the NV electronic ground state. It has been shown that [25] this strain-induced spin-phonon coupling can lead to effective spin-spin interactions mediated by virtual phonons. This phonon-mediated spin-spin coupling strength is about 2g2 =Δ, where g is the coupling strength between a single NV spin and a single phonon via strains and Δ is frequency detuning. For the diamond beam considered in this work, the single phonon coupling is very weak, g ∼ 1 Hz. Thus, the phonon-mediated spin-spin couplings can be completely ignored. [1] Ze-Liang Xiang, Sahel Ashhab, J. Q. You, and Franco Nori, Hybrid quantum circuits: Superconducting circuits interacting with other quantum systems, Rev. Mod. Phys. 85, 623 (2013). [2] J. Verdú, H. Zoubi, Ch. Koller, J. Majer, H. Ritsch, and J. Schmiedmayer, Strong Magnetic Coupling of an Ultracold Gas to a Superconducting Waveguide Cavity, Phys. Rev. Lett. 103, 043603 (2009). [3] Peng-Bo Li and Fu-Li Li, Controlled generation of field squeezing with cold atomic clouds coupled to a superconducting transmission line resonator, Phys. Rev. A 81, 035802 (2010). [4] Anders S. Sørensen, Caspar H. van der Wal, Lilian I. Childress, and Mikhail D. Lukin, Capacitive Coupling of Atomic Systems to Mesoscopic Conductors, Phys. Rev. Lett. 92, 063601 (2004). [5] D. Petrosyan and M. Fleischhauer, Quantum Information Processing with Single Photons and Atomic Ensembles in Microwave Coplanar Waveguide Resonators, Phys. Rev. Lett. 100, 170501 (2008). [6] D. Kielpinski, D. Kafri, M. J. Woolley, G. J. Milburn, and J. M. Taylor, Quantum Interface between an Electrical Circuit and a Single Atom, Phys. Rev. Lett. 108, 130504 (2012). [7] Z. Darázs, Z. Kurucz, O. Kálmán, T. Kiss, J. Fortágh, and P. Domokos, Parametric Amplification of the Mechanical Vibrations of a Suspended Nanowire by Magnetic Coupling to a Bose-Einstein Condensate, Phys. Rev. Lett. 112, 133603 (2014). [8] A. André, D. DeMille, J. M. Doyle, M. D. Lukin, S. E. Maxwell, P. Rabl, R. J. Schoelkopf, and P. Zoller, A coherent all-electrical interface between polar molecules and mesoscopic superconducting resonators, Nat. Phys. 2, 636 (2006). [9] P. Rabl, D. DeMille, J. M. Doyle, M. D. Lukin, R. J. Schoelkopf, and P. Zoller, Hybrid Quantum Processors: Molecular Ensembles as Quantum Memory for Solid State Circuits, Phys. Rev. Lett. 97, 033003 (2006). PHYS. REV. APPLIED 4, 044003 (2015) [10] Karl Tordrup, Antonio Negretti, and Klaus Mølmer, Holographic Quantum Computing, Phys. Rev. Lett. 101, 040501 (2008). [11] L. Childress, A. S. Sørensen, and M. D. Lukin, Mesoscopic cavity quantum electrodynamics with quantum dots, Phys. Rev. A 69, 042302 (2004). [12] Audrey Cottet and Takis Kontos, Spin Quantum Bit with Ferromagnetic Contacts for Circuit QED, Phys. Rev. Lett. 105, 160502 (2010). [13] Pei-Qing Jin, Michael Marthaler, Alexander Shnirman, and Gerd Schon, Strong Coupling of Spin Qubits to a Transmission Line Resonator, Phys. Rev. Lett. 108, 190506 (2012). [14] P. Zhang, Ze-Liang Xiang, and Franco Nori, Spin-orbit qubit on a multiferroic insulator in a superconducting resonator, Phys. Rev. B 89, 115417 (2014). [15] J. Stehlik, Y.-Y. Liu, C. M. Quintana, C. Eichler, T. R. Hartke, and J. R. Petta, Fast Charge Sensing of a Cavity-Coupled Double Quantum Dot Using a Josephson Parametric Amplifier, Phys. Rev. Applied 4, 014018 (2015). [16] Atac Imamoǧlu, Cavity QED Based on Collective Magnetic Dipole Coupling: Spin Ensembles as Hybrid Two-Level Systems, Phys. Rev. Lett. 102, 083602 (2009). [17] J. H. Wesenberg, A. Ardavan, G. A. D. Briggs, J. J. L. Morton, R. J. Schoelkopf, D. I. Schuster, and K. Mølmer, Quantum Computing with an Electron Spin Ensemble, Phys. Rev. Lett. 103, 070502 (2009). [18] D. I. Schuster, A. P. Sears, E. Ginossar, L. DiCarlo, L. Frunzio, J. J. L. Morton, H. Wu, G. A. D. Briggs, B. B. Buckley, D. D. Awschalom, and R. J. Schoelkopf, HighCooperativity Coupling of Electron-Spin Ensembles to Superconducting Cavities, Phys. Rev. Lett. 105, 140501 (2010). [19] Y. Kubo, F. R. Ong, P. Bertet, D. Vion, V. Jacques, D. Zheng, A. Dreau, J.-F. Roch, A. Auffeves, F. Jelezko, J. Wrachtrup, M. F. Barthe, P. Bergonzo, and D. Esteve, Strong Coupling of a Spin Ensemble to a Superconducting Resonator, Phys. Rev. Lett. 105, 140502 (2010). [20] R. Amsüss, Ch. Koller, T. Nöbauer, S. Putz, S. Rotter, K. Sandner, S. Schneider, M. Schramböck, G. Steinhauser, H. Ritsch, J. Schmiedmayer, and J. Majer, Cavity QED with Magnetically Coupled Collective Spin States, Phys. Rev. Lett. 107, 060502 (2011). [21] D. Marcos, M. Wubs, J. M. Taylor, R. Aguado, M. D. Lukin, and A. S. Sørensen, Coupling Nitrogen-Vacancy Centers in Diamond to Superconducting Flux Qubits, Phys. Rev. Lett. 105, 210501 (2010). [22] J. Twamley and S. D. Barrett, Superconducting cavity bus for single nitrogen-vacancy defect centers in diamond, Phys. Rev. B 81, 241202(R) (2010). [23] P. Rabl, S. J. Kolkowitz, F. H. L. Koppens, J. G. E. Harris, P. Zoller, and M. D. Lukin, A quantum spin transducer based on nanoelectromechanical resonator arrays, Nat. Phys. 6, 602 (2010). [24] Brian Julsgaard, Cécile Grezes, Patrice Bertet, and Klaus Mølmer, Quantum Memory for Microwave Photons in an Inhomogeneously Broadened Spin Ensemble, Phys. Rev. Lett. 110, 250503 (2013). 044003-11 LI et al. PHYS. REV. APPLIED 4, 044003 (2015) [25] S. D. Bennett, N. Y. Yao, J. Otterbach, P. Zoller, P. Rabl, and M. D. Lukin, Phonon-Induced Spin-Spin Interactions in Diamond Nanostructures: Application to Spin Squeezing, Phys. Rev. Lett. 110, 156402 (2013). [26] S. Carretta, A. Chiesa, F. Troiani, D. Gerace, G. Amoretti, and P. Santini, Quantum Information Processing with Hybrid Spin-Photon Qubit Encoding, Phys. Rev. Lett. 111, 110501 (2013). [27] Martin Leijnse and Karsten Flensberg, Coupling Spin Qubits via Superconductors, Phys. Rev. Lett. 111, 060501 (2013). [28] Alexey A. Kovalev, Amrit De, and Kirill Shtengel, Spin Transfer of Quantum Information between Majorana Modes and a Resonator, Phys. Rev. Lett. 112, 106402 (2014). [29] L. J. Zou, D. Marcos, S. Diehl, S. Putz, J. Schmiedmayer, J. Majer, and P. Rabl, Implementation of the Dicke Lattice Model in Hybrid Quantum System Arrays, Phys. Rev. Lett. 113, 023603 (2014). [30] Christopher O’Brien, Nikolai Lauk, Susanne Blum, Giovanna Morigi, and Michael Fleischhauer, Interfacing Superconducting Qubits and Telecom Photons via a Rare-Earth-Doped Crystal, Phys. Rev. Lett. 113, 063603 (2014). [31] E. R. MacQuarrie, T. A. Gosavi, N. R. Jungwirth, S. A. Bhave, and G. D. Fuchs, Mechanical Spin Control of Nitrogen-Vacancy Centers in Diamond, Phys. Rev. Lett. 111, 227602 (2013). [32] S. D. Bennett, S. Kolkowitz, Q. P. Unterreithmeier, P. Rabl, A. C. Bleszynski Jayich, J. G. E. Harris, and M. D. Lukin, Measuring mechanical motion with a single spin, New J. Phys. 14, 125004 (2012). [33] Z. Y. Xu, Y. M. Hu, W. L. Yang, M. Feng, and J. F. Du, Deterministically entangling distant nitrogen-vacancy centers by a nanomechanical cantilever, Phys. Rev. A 80, 022335 (2009). [34] Mohamed Almokhtar, Masazumi Fujiwara, Hideaki Takashima, and Shigeki Takeuchi, Numerical simulations of nanodiamond nitrogen-vacancy centers coupled with tapered optical fibers as hybrid quantum nanophotonic devices, Opt. Express 22, 20045 (2014). [35] Bo-Bo Wei, Christian Burk, Jörg Wrachtrup, and Ren-Bao Liu, Magnetic ordering of nitrogen-vacancy centers in diamond via resonator-mediated coupling, EPJ Quantum Technology 2, 18 (2015). [36] Burm Baek, William H. Rippard, Matthew R. Pufall, Samuel P. Benz, Stephen E. Russek, Horst Rogalla, and Paul D. Dresselhaus, Spin-Transfer Torque Switching in Nanopillar Superconducting-Magnetic Hybrid Josephson Junctions, Phys. Rev. Applied 3, 011001 (2015). [37] Y. Kubo, C. Grezes, A. Dewes, T. Umeda, J. Isoya, H. Sumiya, N. Morishita, H. Abe, S. Onoda, T. Ohshima, V. Jacques, A. Dréau, J.-F. Roch, I. Diniz, A. Auffeves, D. Vion, D. Esteve, and P. Bertet, Hybrid Quantum Circuit with a Superconducting Qubit Coupled to a Spin Ensemble, Phys. Rev. Lett. 107, 220501 (2011). [38] Y. Tao, J. M. Boss, B. A. Moores, and C. L. Degen, Singlecrystal diamond nanomechanical resonators with quality factors exceeding one million, Nat. Commun. 5, 3638 (2014). [39] Preeti Ovartchaiyapong, Kenneth W. Lee, Bryan A. Myers, and Ania C. Bleszynski Jayich, Dynamic strain-mediated coupling of a single diamond spin to a mechanical resonator, Nat. Commun. 5, 4429 (2014). [40] J. Teissier, A. Barfuss, P. Appel, E. Neu, and P. Maletinsky, Strain Coupling of a Nitrogen-Vacancy Center Spin to a Diamond Mechanical Oscillator, Phys. Rev. Lett. 113, 020503 (2014). [41] Patrik Rath, Svetlana Khasminskaya, Christoph Nebel, Christoph Wild, and Wolfram H. P. Pernice, Diamondintegrated optomechanical circuits, Nat. Commun. 4, 1690 (2013). [42] Michael J. Burek, Daniel Ramos, Parth Patel, Ian W. Frank, and Marko Loncar, Nanomechanical resonant structures in single-crystal diamond, Appl. Phys. Lett. 103, 131904 (2013). [43] P. Ovartchaiyapong, L. M. A. Pascal, B. A. Myers, P. Lauria, and A. C. Bleszynski Jayich, High quality factor singlecrystal diamond mechanical resonators, Appl. Phys. Lett. 101, 163505 (2012). [44] Behzad Khanaliloo, Harishankar Jayakumar, Aaron C. Hryciw, David P. Lake, Hamidreza Kaviani, and Paul E. Barclay, Diamond nanobeam waveguide optomechanics, arXiv:1502.01788v1. [45] Marcus W. Doherty, Neil B. Manson, Paul Delaney, Fedor Jelezko, Jörg Wrachtrup, and Lloyd C. L. Hollenberg, The nitrogen-vacancy colour centre in diamond, Phys. Rep. 528, 1 (2013). [46] M. Gao, C. W. Wu, Z. J. Deng, W. J. Zou, L. Zhou, C. Z. Li, and X. B. Wang, Controllable strong coupling between individual spin qubits and a transmission line resonator via nanomechanical resonators, Phys. Lett. A 376, 595 (2012). [47] C. A. Regal and K. W. Lehnert, From cavity electromechanics to cavity optomechanics, J. Phys. Conf. Ser. 264, 012025 (2011). [48] N. Didier and R. Fazio, Putting mechanics into circuit quantum electrodynamics, C.R. Phys. 13, 470 (2012). [49] J. D. Teufel, Dale Li, M. S. Allman, K. Cicak, A. J. Sirois, J. D. Whittaker, and R. W. Simmonds, Circuit cavity electromechanics in the strong-coupling regime, Nature (London) 471, 204 (2011). [50] Quirin P. Unterreithmeier, Eva M. Weig, and Jörg P. Kotthaus, Universal transduction scheme for nanomechanical systems based on dielectric forces, Nature (London) 458, 1001 (2009). [51] P. Rabl, P. Cappellaro, M. V. Gurudev Dutt, L. Jiang, J. R. Maze, and M. D. Lukin, Strong magnetic coupling between an electronic spin qubit and a mechanical resonator, Phys. Rev. B 79, 041302(R) (2009). [52] T. Faust, J. Rieger, M. J. Seitner, J. P. Kotthaus, and E. M. Weig, Coherent control of a classical nanomechanical twolevel system, Nat. Phys. 9, 485 (2013). [53] Simon Rips and Michael J. Hartmann, Quantum Information Processing with Nanomechanical Qubits, Phys. Rev. Lett. 110, 120503 (2013). [54] See Appendixes A, B, and C for more details. [55] J. Venermo and A. Sihvola, Dielectric polarizability of circular cylinder, J. Electrost. 63, 101 (2005). 044003-12 HYBRID QUANTUM DEVICE BASED ON NV CENTERS … [56] T. J. Kippenberg and K. J. Vahala, Cavity optomechanics: Back-action at the mesoscale, Science 321, 1172 (2008). [57] F. Marquardt and S. M. Girvin, Optomechanics, Physics 2, 40 (2009). [58] I. Wilson-Rae, N. Nooshi, W. Zwerger, and T. J. Kippenberg, Theory of Ground State Cooling of a Mechanical Oscillator Using Dynamical Backaction, Phys. Rev. Lett. 99, 093901 (2007). [59] F. Marquardt, J. P. Chen, A. A. Clerk, and S. M. Girvin, Quantum Theory of Cavity-Assisted Sideband Cooling of Mechanical Motion, Phys. Rev. Lett. 99, 093902 (2007). [60] Yong-Chun Liu, Yun-Feng Xiao, Xingsheng Luan, and Chee Wei Wong, Dynamic Dissipative Cooling of a Mechanical Resonator in Strong Coupling Optomechanics, Phys. Rev. Lett. 110, 153606 (2013). [61] J. D. Teufel, T. Donner, Dale Li, J. W. Harlow, M. S. Allman, K. Cicak, A. J. Sirois, J. D. Whittaker, K. W. Lehnert, and R. W. Simmonds, Sideband cooling of micromechanical motion to the quantum ground state, Nature (London) 475, 359 (2011). [62] Jasper Chan, T. P. Mayer Alegre, Amir H. Safavi-Naeini, Jeff T. Hill, Alex Krause, Simon Gröblacher, Markus Aspelmeyer, and Oskar Painter, Laser cooling of a nanomechanical oscillator into its quantum ground state, Nature (London) 478, 89 (2011). [63] S. Hong, M. S. Grinolds, P. Maletinsky, R. L. Walsworth, M. D. Lukin, and A. Yacoby, Coherent, mechanical control of a single electronic spin, Nano Lett. 12, 3920 (2012). [64] O. Arcizet, V. Jacques, A. Siria, P. Poncharal, P. Vincent, and S. Seidelin, A single nitrogen-vacancy defect coupled to a nanomechanical oscillator, Nat. Phys. 7, 879 (2011). [65] Zhang Qi Yin, Tong Cang Li, Xiang Zhang, and L. M. Duan, Large quantum superpositions of a levitated PHYS. REV. APPLIED 4, 044003 (2015) nanodiamond through spin-optomechanical coupling, Phys. Rev. A 88, 033614 (2013). [66] G. Balasubramanian, P. Neumann, D. Twitchen, M. Markham, R. Kolesov, N. Mizuochi, J. Isoya, J. Achard, J. Beck, J. Tissler, V. Jacques, P. R. Hemmer, F. Jelezko, and J. Wrachtrup, Ultralong spin coherence time in isotopically engineered diamond, Nat. Mater. 8, 383 (2009). [67] K. Bergmann, H. Theuer, and B. W. Shore, Coherent population transfer among quantum states of atoms and molecules, Rev. Mod. Phys. 70, 1003 (1998). [68] K. Eckert, M. Lewenstein, R. Corbalán, G. Birkl, W. Ertmer, and J. Mompart, Three-level atom optics via the tunneling interaction, Phys. Rev. A 70, 023606 (2004). [69] K. Eckert, J. Mompart, R. Corbalán, M. Lewenstein, and G. Birkl, Three level atom optics in dipole traps and waveguides, Opt. Commun. 264, 264 (2006). [70] F. Dreisow, A. Szameit, M. Heinrich, R. Keil, S. Nolte, A. Tünnermann, and S. Longhi, Adiabatic transfer of light via a continuum in optical waveguides, Opt. Lett. 34, 2405 (2009). [71] J. R. Johansson, P. N. Nation, and F. Nori, Qutip 2: A python framework for the dynamics of open quantum systems, Comput. Phys. Commun. 184, 1234 (2013). [72] W. Weaver, S. P. Timoshenko, and D. H. Young, Vibration Problems in Engineering, 5th ed. (Wiley, New York, 1990). [73] L. D. Landau and E. M. Lifshitz, Theory of Elasticity (Butterworth-Heinemann, Oxford, 1986). [74] R. N. Simons, Coplanar Waveguide Circuits, Components, and Systems (Wiley, New York, 2004). [75] R. N. Simons and R. K. Arora, Coupled slot line field components, IEEE Trans. Microwave Theory Tech. 30, 1094 (1982). [76] C. L. Degen, H. J. Mamin, M. Poggio, and D. Rugar, Nuclear magnetic resonance imaging with 90-nm resolution, Nat. Nanotechnol. 2, 301 (2007). 044003-13