Yuan Zhou^,1,∗, Dong-Yan Lü¹, Guang-Hui Wang², Yan-Hua Fu¹, Ming-Yao He¹, Hong-Tao Ren³¹School of Science, Hubei University of Automotive Technology, Shiyan 442002, China
²Department of of Automobile Engineering, Hubei University of Automotive Technology, Shiyan 442002, China
³School of Materials Science and Engineering, Liaocheng University, Liaocheng 252000, China

First author contact: ∗Author to whom any correspondence should be addressed.
Received:2021-01-17Revised:2021-03-5Accepted:2021-03-5Online:2021-04-09


Abstract
It remains a great challenge to realize direct manipulation of a nitrogen-vacancy (NV) spin at the single-quantum level with a microwave (MW) cavity. As an alternative, a hybrid system with the spin-phonon-photon triple interactions mediated by a squeezed cantilever-type harmonic resonator is proposed. According to the general mechanical parametric amplification of this in-between phonon mode, the direct spin-phonon and photon-phonon couplings are both exponentially enhanced, which can even further improve the coherent manipulation of a single NV spin and MW photon with a higher efficiency. In view of this triple system with enhanced couplings and the additional sideband adjustable designs, this scheme may provide a more efficient phonon-mediated platform to bridge or manipulate the MW quantum and a single electron spin coherently. It is also hoped to evoke wider applications in the areas of quantum state transfer and preparation, ultrasensitive detection and quantum nondestructive measurement, etc.
Keywords： NV center;mechanical parametric amplification;quantum manipulation


PDF (0KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Yuan Zhou, Dong-Yan Lü, Guang-Hui Wang, Yan-Hua Fu, Ming-Yao He, Hong-Tao Ren. Improvement on the manipulation of a single nitrogen-vacancy spin and microwave photon at single-quantum level. Communications in Theoretical Physics, 2021, 73(6): 065101- doi:10.1088/1572-9494/abec3a

1. Introduction

Implementing quantum manipulations and quantum simulations in the microwave (MW) domain has stimulated huge progress in recent times [1-3]. Designing a hybrid quantum system in this special electromagnetic frequency region and implementing efficient quantum information processing (QIP) is gradually becoming recognized as one of the most reliable technological routes [4-7], whereupon the most central task lies in the achievement of perfect manipulation of atoms or artificial atoms utilizing these kinds of MW quantum devices [8], i.e. the MW cavity [9, 10], the MW waveguide [11, 12], the superconductive circuit [13, 14], etc. In addition, it is also regarded as a promising attempt to perform well-controlled quantum measurement or quantum sensing with ultrasensitivity [15, 16] through fabricating MW-based hybrid setups or devices [17]. With the rapid development of high-quality MW devices integrating technology, the prospect of using MW to dominate QIP is also becoming clearer [18, 19].

Regarding today’s mainstream schemes on QIP, and comparing all kinds of quantum units applicable for manipulation, the nitrogen-vacancy (NV) center in diamond is particularly attractive due to its excellent spin properties, i.e. ultrasensitivity to magnetic fields, fast MW manipulation, optical preparation and detection, and long coherence time even at ambient conditions [20-23]. The long spin coherence of NVs can also evoke many intriguing applications at the single-photon as well as macroscopic level. For example, it is a high-quality MW emitter in a diamond maser scheme proposed by Jin et al [24], which was also experimentally demonstrated by Breeze group in 2018 [25], etc. However, due to the inherent limitation of magnetic coupling mechanism and magnetic field localization degree, it remains a considerable huge challenge to realize a direct strong magnetic coupling to single NV spin in a given MW device [26, 27], especially in the MW cavity [28].

In a different way, here we discuss a different approach and propose a hybrid system with phonon-mediated spin-phonon-photon triple interactions, in which a squeezed quantum mechanical mode can not only interact with an NV spin through the first-order gradient magnetic field, but also capacitively couple to the MW cavity [29-32]. Utilizing the general mechanical parametric amplification (MPA), we can significantly enhance two pairs of direct spin-phonon and photon-phonon couplings both with the exponential rate (∼e^r) in a new squeezed frame [29, 33-38]. For a given squeezing parameter r, we can further adjust both coupled pairs to obey the Jaynes-Cummings (JC)-type or anti-Jaynes-Cummings-type interactions, only through modulating the parameters Δ_c/q. We state both Δ_c/q are conveniently tunable parameters with respect to the processes of MW cavity linearization and dressed-state NVs. In view of this hybrid quantum system with adjustable sidebands and enhanced couplings, our scheme may provide a more efficient phonon-mediated quantum interface to bridge the MW quanta and a single electron spin coherently, and even provide a fruitful route towards the manipulation of a single NV spin and MW photon at the single-quantum level with an improved efficiency. We also hope this result may evoke wider applications in the areas of quantum state transfer [39, 40], ultrasensitive detection [41], quantum nondestructive measurement (QND) [42], etc.

2. Setup

The whole schematic design and relevant physical image are all illustrated in figures 1(a)-(c). Working as a mediator, a cantilever-type mechanical resonator (MR) interacts with the MW cavity and a single NV spin, and its spring constant is also periodically tunable, so we can obtain the second-order nonlinear interaction on this mechanical mode [29, 38, 43-45]. In a new squeezed frame, the quantum fluctuation of this MR amplitude is amplified exponentially with ∼e^r, known as the MPA [29]. As specified in figure 1(a), there are two essential direct strong-coupling parts to form this triple photon-phonon-spin interaction system. For the first part, a sharp nanomagnetic tip is attached at the end of the cantilever with a designed tunable spring constant, just above which a single NV center is set near to this nanomagnet; we can then obtain the strong magnetic spin-phonon coupling with the dressed-state expression $\sim \lambda {{\rm{e}}}^{r}(\hat{b}+{\hat{b}}^{\dagger })({\hat{\sigma }}^{+}+{\hat{\sigma }}^{-})$ in a squeezed frame [29]. We also state another available strain coupling scheme is also equivalent to this magnetic spin-phonon coupling setup, i.e. the investigations [46-50]. For the second part, we can also perform the strong coupling through the capacitive interaction between the MW cavity and the MR [32, 51]. As shown in figure 1(a), a superconducting coplanar cavity is connected to the MR through the capacitor circuit. Then, we can obtain the so-called cavity optomechanical coupling system, and acquire an exponentially enhanced interaction as $\sim g{{\rm{e}}}^{r}(\hat{b}+{\hat{b}}^{\dagger })(\hat{a}+{\hat{a}}^{\dagger })$, after the cavity mode is linearized [32, 51].

Figure 1.

New window|Download| PPT slide
Figure 1.(Color online) (a) Potential photon-phonon-spin hybrid setup with respect to the triple coupling system mediated by a squeezed mechanical mode. A sharp nano magnetic tip is attached at the end of the silicon cantilever with its tunable spring constant, above which a single NV center is set near to this magnet. A superconducting coplanar MW cavity also couples to this mechanical cantilever through capacitive interactions. As a result, it is available for this MR to achieve capacitive coupling to the MW cavity and magnetic coupling to the single NV spin, respectively. (b) Schematic of enhancing the indirect coupling between the MW cavity and single NV spin through MPA. (c) Ground-state energy diagram for a single NV center in the basis of bare state and dressed state.

3. Mechanism of this triple coupling system

As illustrated in figure 1(a), we assume that the spring constant of this mechanical cantilever is a periodically tunable parameter with the frequency ω_χ. According to the statement in the detail of appendix A, we can obtain the Hamiltonian of this quantized cantilever with the second-order nonlinear interaction
(1)$\begin{eqnarray}{\hat{H}}_{\mathrm{mr}}\simeq {\omega }_{m}{\hat{b}}^{\dagger }\hat{b}-\displaystyle \frac{{{\rm{\Omega }}}_{\chi }}{2}({\hat{b}}^{\dagger 2}{{\rm{e}}}^{2{\rm{i}}{\omega }_{\chi }t}+{\hat{b}}^{2}{{\rm{e}}}^{2{\rm{i}}{\omega }_{\chi }t}+{\rm{h}}.{\rm{c}}.).\end{eqnarray}$
In equation (1), we have discarded the energy-shift item caused by this nonlinear drive amplitude Ω_χ.

According to figure 1(c), for the first NV-MR coupling pair, we apply the dichromatic MW fields to dress the ground state of the NV center, and acquire the efficient spin-phonon transition channel according to appendix B. In a selected dressed-state subspace {∣d⟩, ∣g⟩}, the effective spin-phonon interaction can be rewritten as the JC-type Hamiltonian
(2)$\begin{eqnarray}{\hat{H}}_{\mathrm{sp}}^{d}\simeq \displaystyle \frac{{\omega }_{0}}{2}{\hat{\sigma }}_{z}+\lambda ({\hat{b}}^{\dagger }{\hat{\sigma }}^{-}+\hat{b}{\hat{\sigma }}^{+}).\end{eqnarray}$
Here, we define ω₀ ≡ ω_dg, ${\hat{\sigma }}_{z}\equiv {\hat{\sigma }}_{{dd}}-{\hat{\sigma }}_{{gg}}$, and ${\hat{\sigma }}^{\pm }\equiv {\hat{\sigma }}_{{dg}}^{\pm }$ for simplicity.

For another pair of MR-cavity coupling, we assume the MW cavity interacts with the MR capacitively, which forms the so-called cavity optomechanical coupling pair. This type of interaction is described as
(3)$\begin{eqnarray}\begin{array}{rcl}{\hat{H}}_{\mathrm{op}} & \simeq & [{\omega }_{c}+{g}_{0}({\hat{b}}^{\dagger }+\hat{b})]{\hat{a}}^{\dagger }\hat{a}\\ & & +\,\beta ({\hat{a}}^{\dagger }{{\rm{e}}}^{{\rm{i}}{\omega }_{L}t}+{\rm{h}}.{\rm{c}}.).\end{array}\end{eqnarray}$
In equation (3), the first two items describe the general optomechanical coupling with the fundamental cavity frequency ω_c and the optomechanical coupling strength g₀, and the last item means the classical MW laser drive to this cavity mode with driving amplitude $\beta$ and frequency ω_L. Utilizing the standard linearization method with $\hat{a}=\alpha +\delta \hat{a}$ and ${\hat{a}}^{\dagger }={\alpha }^{* }+\delta {\hat{a}}^{\dagger }$ [32, 51], we can linearize this optomechanical interaction (3) as
(4)$\begin{eqnarray}{\hat{H}}_{\mathrm{op}}^{L}\approx {\delta }_{c}{\hat{a}}^{\dagger }\hat{a}+g({\hat{a}}^{\dagger }\hat{b}+\hat{a}{\hat{b}}^{\dagger }),\end{eqnarray}$
in which the parameters are assumed as δ_c = ω_c − ω_L, $\alpha$^* = $\alpha$, and g = g₀$\alpha$, and we have set $\delta {\hat{a}}^{(\dagger )}\equiv {\hat{a}}^{(\dagger )}$ for simplicity.

As a result, together with equations (1), (2), and (4), we can obtain the whole Hamiltonian for describing this phonon-mediated three-body hybrid system:
(5)$\begin{eqnarray}{\hat{H}}_{w}={\hat{H}}_{\mathrm{mr}}+{\hat{H}}_{\mathrm{op}}^{L}+{\hat{H}}_{\mathrm{sp}}^{d}.\end{eqnarray}$
Then, we apply an unitary transformation ${\hat{U}}_{0}=\exp [{\rm{i}}{\hat{H}}_{0}t]$ to equation (5) with ${\hat{H}}_{0}={\omega }_{\chi }({\hat{b}}^{\dagger }\hat{b}+{\hat{a}}^{\dagger }\hat{a}+{\hat{\sigma }}_{z}/2)$. In this new rotating frame, we can obtain
(6)$\begin{eqnarray}\begin{array}{rcl}{\hat{H}}_{w}^{{\prime} } & \simeq & {\delta }_{m}{\hat{b}}^{\dagger }\hat{b}-\displaystyle \frac{{{\rm{\Omega }}}_{\chi }}{2}({\hat{b}}^{\dagger 2}+{\hat{b}}^{2})+{{\rm{\Delta }}}_{c}{\hat{a}}^{\dagger }\hat{a}+\displaystyle \frac{{{\rm{\Delta }}}_{q}}{2}{\hat{\sigma }}_{z}\\ & & +\,g({\hat{a}}^{\dagger }\hat{b}+\hat{a}{\hat{b}}^{\dagger })+\lambda ({\hat{b}}^{\dagger }{\hat{\sigma }}^{-}+\hat{b}{\hat{\sigma }}^{+}),\end{array}\end{eqnarray}$
with δ_m = ω_m − ω_χ, Δ_c = δ_c − ω_χ, and Δ_q = ω₀ − ω_χ.

Considering this Hamiltonian (6) accompanied by the second-order nonlinear interaction of this mechanical freedom, we can utilize the standard unitary transformation ${\hat{U}}_{s}(r)=\exp [r({\hat{b}}^{2}-{\hat{b}}^{\dagger 2})/2]$ to diagonalize this mechanical mode belonging to ${\hat{H}}_{w}^{{\prime} }$, where the squeezing parameter r is defined by $\tanh 2r={{\rm{\Omega }}}_{\chi }/{\delta }_{m}$. In this squeezed frame, we can obtain the Hamiltonian with the new expression
(7)$\begin{eqnarray}\begin{array}{rcl}{\hat{H}}_{w}^{S} & \approx & {{\rm{\Delta }}}_{m}{\hat{b}}^{\dagger }\hat{b}+{{\rm{\Delta }}}_{c}{\hat{a}}^{\dagger }\hat{a}+\displaystyle \frac{{{\rm{\Delta }}}_{q}}{2}{\hat{\sigma }}_{z}\\ & & +\,g{{\rm{e}}}^{r}(\hat{b}+{\hat{b}}^{\dagger })(\hat{a}+{\hat{a}}^{\dagger })+\lambda {{\rm{e}}}^{r}(\hat{b}+{\hat{b}}^{\dagger }){\hat{\sigma }}_{x},\end{array}\end{eqnarray}$
with ${{\rm{\Delta }}}_{m}={\delta }_{m}/\cosh 2r$. Here, it is noted that we have omitted the weak disturbance items during this diagonalization process, because of their diminishing coefficients (∝ e^−r) as we increase the squeezing parameter r. In the interaction picture (IP), we obtain
(8)$\begin{eqnarray}\begin{array}{rcl}{\hat{H}}_{\mathrm{IP}}^{S}(t) & \approx & {G}_{1}(\hat{b}\hat{a}{{\rm{e}}}^{-{\rm{i}}{\delta }_{1}t}+{\hat{b}}^{\dagger }\hat{a}{{\rm{e}}}^{-{\rm{i}}{\delta }_{2}t}+{\rm{h}}.{\rm{c}}.)\\ & & +\,{G}_{2}(\hat{b}{\hat{\sigma }}^{-}{{\rm{e}}}^{-{\rm{i}}{\delta }_{3}t}+{\hat{b}}^{\dagger }{\hat{\sigma }}^{-}{{\rm{e}}}^{-{\rm{i}}{\delta }_{4}t}+{\rm{h}}.{\rm{c}}.),\end{array}\end{eqnarray}$
with G₁ = ge^r, G₂ = λe^r, δ₁ = Δ_c + Δ_m, δ₂ = Δ_c − Δ_m, δ₃ = Δ_q + Δ_m, and δ₄ = Δ_q − Δ_m. For a realistic condition, the dynamical evolution of this hybrid system is dominated by the master equation effectively as follows:
(9)$\begin{eqnarray}\begin{array}{rcl}\displaystyle \frac{{\rm{d}}\hat{\rho }}{{\rm{d}}t} & = & {\rm{i}}[\hat{\rho },{\hat{H}}_{w}^{S}({\hat{H}}_{\mathrm{IP}}^{S})]+{\kappa }_{1}{ \mathcal D }[\hat{a}]\hat{\rho }\\ & & +\,{\kappa }_{\mathrm{eff}}^{S}{ \mathcal D }[\hat{b}]\hat{\rho }+\gamma { \mathcal D }[{\hat{\sigma }}_{z}]\hat{\rho },\end{array}\end{eqnarray}$
where ${ \mathcal D }[\hat{o}]\hat{\rho }=\hat{o}\hat{\rho }{\hat{o}}^{\dagger }-\hat{\rho }{\hat{o}}^{\dagger }\hat{o}/2-{\hat{o}}^{\dagger }\hat{o}\hat{\rho }/2$, $\hat{\rho }$ is the density matrix of this three-body quantum system, and κ₁, ${\kappa }_{\mathrm{eff}}^{S}$, and γ correspond to the MW cavity dissipation rate, the effective squeezed MR damping rate, and the NV dephasing rate, respectively.

In this squeezed frame, the mechanical mode is fast and strongly damped with a large effective dissipation rate ${\kappa }_{\mathrm{eff}}^{S}$ [38, 52]. In this case, we can inevitably maintain this phonon mode in its ground state ∣0⟩ and neglect the populations of any excited modes. This squeezed mechanical mode is then considered as an artificial or engineered reservoir of the heat phonon. In contrast, we can also apply an auxiliary cooling system or design to suppress its effective dissipation rate at the same magnitude as ∼ κ₁, according to investigations [52-56].

Considering Hamiltonian (8), we also note, for a given cantilever with a certain squeezing parameter r, that Δ_m is very small and untunable. However, Δ_c and Δ_q are both tunable parameters, which can be implemented via the cavity-resonator linearization process and NV-resonator MW dressed process. Through modulating both relevant parameters, we can obtain an effective triple Hamiltonian with distinct expressions. Utilizing this, we can also carry out different simulations for multiple applications. For example, in case (i), when Δ_c − Δ_m ≃ Δ_q − Δ_m ≃ 0, the Hamiltonian reduces to the beam-split-type interaction effectively, with $\sim {\hat{b}}^{\dagger }[{G}_{1}\hat{a}+{G}_{2}{\hat{\sigma }}^{-}]+{\rm{h}}.{\rm{c}}.$, which is applicable for establishing a quantum interface to bridge the MW photons and NV spin efficiently. In case (ii), when Δ_c + Δ_m ≃ Δ_q − Δ_m ≃ 0 (or Δ_c − Δ_m ≃ Δ_q + Δ_m ≃ 0), the Hamiltonian reduces to the types of $\sim {\hat{b}}^{\dagger }[{G}_{1}{\hat{a}}^{\dagger }\,+{G}_{2}{\hat{\sigma }}^{-}]+{\rm{h}}.{\rm{c}}.$ (or $\sim {\hat{b}}^{\dagger }[{G}_{1}\hat{a}+{G}_{2}{\hat{\sigma }}^{+}]+{\rm{h}}.{\rm{c}}.$), which is suitable for engineering the MW photons and NV spin into the nonclassical steady states, such as certain entanglement or quantum correlation. In case (iii), once Δ_c + Δ_m ≃Δ_q + Δ_m ≃ 0, the Hamiltonian is also reduced to the form $\sim {\hat{b}}^{\dagger }[{G}_{1}{\hat{a}}^{\dagger }+{G}_{2}{\hat{\sigma }}^{+}]+{\rm{h}}.{\rm{c}}.$ effectively, which can also be applied to the preparation of special nonclassical states; i.e. governed by this type of interaction, an investigation on generating the Bell-type maximum entanglement through a fast dissipation process was recently reported [57]. Taking cases (i) and (ii) as examples, we present specific discussions and numerical simulations of relevant applications in the following section.

4. Applications

4.1. Efficiency improvement on quantum state transfer between the NV spin and MW photon

For case (i), the effective Hamiltonian reduces to
(10)$\begin{eqnarray}{\hat{H}}_{\left(({\rm{i}})\right)}^{S}\simeq ({G}_{1}{\hat{b}}^{\dagger }\hat{a}+{G}_{2}\hat{b}{\hat{\sigma }}^{+}+{\rm{h}}.{\rm{c}}.).\end{eqnarray}$
Utilizing this interaction, we can implement the target of quantum state transfer from the MW photons to an NV spin mediated by this squeezed phonon mode, and vice versa [39, 40]. Therefore, the first application of this scheme is focused on establishing a quantum spin-photon interface with further improved efficiency.

For an ideal condition without considering any decoherence factors, this triple-interaction system is governed by the Schrödinger equation ${\rm{i}}\partial | {\rm{\Psi }}(t)\rangle /\partial t={\hat{H}}_{\left(({\rm{i}})\right)}^{S}| {\rm{\Psi }}(t)\rangle $. Taking the quantum state transfer from the MW photon to NV spin at the single-quantum level as an example, we first assume this photon mode is reduced to the lowest subspace {∣0⟩, ∣1⟩} for simplicity. Without loss of generality, the initial state of the MW cavity is assumed as ∣ψ⟩_MW = $\alpha$∣0⟩ + $\beta$∣1⟩ with the normalizing condition ∣$\alpha$∣² + ∣$\beta$∣² = 1, and the mechanical mode and NV spin are both initialized in their ground states ∣0⟩_m and ∣g⟩_s. Therefore, the initial state of this whole system is initialized in the assumed state ∣Ψ(0)⟩ = ∣ψ⟩_MW ⨂ ∣0⟩_m ⨂ ∣g⟩_s, and the relevant quantum information to be transferred to NV spin is initially coded in the freedom of this MW mode with the given coefficients $\alpha$ and $\beta$.

Here we assume g ≃ λ and G₁ = − G₂ = G = ge^r for simplicity, and this whole system is dominated by Hamiltonian (10). Through solving this Schrödinger equation, we can obtain its evolved state ∣Ψ(t)⟩ analytically at any arbitrary time t. At the special moment with τ_s ∼ π/G = π/(ge^r), we can obtain the target state with ∣Ψ(τ_s)⟩ = ∣0⟩_MW ⨂ ∣0⟩_m ⨂ ∣ψ⟩_s and ∣ψ⟩_s = $\alpha$∣g⟩ + $\beta$∣d⟩. Then, we can reach the target of transferring the arbitrary quantum information of MW mode to NV spin and vice versa. Furthermore, its transfer efficiency will be further improved significantly by this squeezed mechanical mode owing to the amplification of both coupling strength and rate e^r [29].

To indicate this physical process more clearly and intuitively, we need to carry out the relevant numerical simulation on it. Considering all adverse decoherence factors in our numerical simulation, its dynamical evolution is governed by the master equation ${\rm{d}}\hat{\rho }/{\rm{d}}t\,={\rm{i}}[\hat{\rho },{\hat{H}}_{\left(({\rm{i}})\right)}^{S}]+{\kappa }_{1}{ \mathcal D }[\hat{a}]\hat{\rho }+{\kappa }_{\mathrm{eff}}^{S}{ \mathcal D }[\hat{b}]\hat{\rho }+\gamma { \mathcal D }[{\hat{\sigma }}_{z}]\hat{\rho }$. In general, here we set $\alpha$ = 0.6, $\beta$ = 0.8, κ₁ ∼ γ ∼ 0.01g, and ${\kappa }_{\mathrm{eff}}^{S}\sim 0.1g$ for simplicity; then, we can solve this master equation numerically and achieve dynamical fidelity according to the definition $F=\langle {\rm{\Psi }}({\tau }_{s})| \hat{\rho }(t)| {\rm{\Psi }}({\tau }_{s})\rangle $, namely the overlap between the target state ∣Ψ(τ_s)⟩ and the dynamically evolved states $\hat{\rho }(t)=| {\rm{\Psi }}(t)\rangle \langle {\rm{\Psi }}(t)| $ at arbitrary moment t.

The corresponding numerical results are plotted in figure 2, and these numerical results are highly consistent with our analysis above. We note that the fidelity can be near to 1 infinitely in the first few periods of oscillation in time, regardless of whether the squeezing parameter is set as r = 1, r = 3, or r = 5, and even with all of the decoherence factors being considered in this master equation. In addition, the relevant state-transfer efficiency and velocity can also be further improved significantly as we increase the squeezing parameter r from r = 0 to r = 5, and its shortest time interval for transfer satisfies τ_s ∝ e^−r.

Figure 2.

New window|Download| PPT slide
Figure 2.(Color online) Dynamical fidelity for this photon-spin state transfer process. The parameters are set as g = λ, r = 0 the black solid line with solid spheres, r = 1 the olive solid line with solid squares, r = 3 the blue solid line with open stars, and r = 5 the orange solid line, κ₁ ∼ γ ∼ 0.01g, and ${\kappa }_{\mathrm{eff}}^{S}\sim 0.1g$.

4.2. Engineering NV spin and MW photon into the quantum correlated state indirectly by the squeezed phonon mode

In contrast, for case (ii), the effective Hamiltonian obeys the forms of
(11)$\begin{eqnarray}{\hat{H}}_{(\mathrm{ii})}^{S(1)}\simeq ({G}_{1}{\hat{b}}^{\dagger }{\hat{a}}^{\dagger }+{G}_{2}\hat{b}{\hat{\sigma }}^{+}+{\rm{h}}.{\rm{c}}.),\end{eqnarray}$
or
(12)$\begin{eqnarray}{\hat{H}}_{(\mathrm{ii})}^{S\left(2\right)}\simeq ({G}_{1}{\hat{b}}^{\dagger }\hat{a}+{G}_{2}{\hat{b}}^{\dagger }{\hat{\sigma }}^{+}+{\rm{h}}.{\rm{c}}.).\end{eqnarray}$
Dominated by both types of triple interactions, we can also carry out the task of engineering the MW photon and NV spin into the steady nonclassical state, utilizing different routes such as the adiabatic evolution or fast dissipation, etc [38, 52, 57-60]. Taking Hamiltonian (11) as an example, we first assume G₁ is opposite to G₂ in phase and G₁ ≠ G₂; then, we can rewrite this Hamiltonian as ${\hat{H}}_{(\mathrm{ii})}^{S(1)}\simeq \hat{b}({G}_{2}{\hat{\sigma }}^{+}-{G}_{1}\hat{a})+{\rm{h}}.{\rm{c}}.$. We note this type of interaction is quite similar to the general Hamiltonian for engineering the two-mode quantum correlated state assisted by fast dissipation [60-63]. In this scheme, we stress that the photon-phonon and spin-phonon coupling can be exponentially enhanced in a squeezed frame; however, the realistic dissipation rate of this squeezed phonon mode will be inevitably amplified without any additional cooling design. To engineer this pair of NV and MW photon into the relevant nonclassical state, we are confident that the fast dissipation route is a more suitable choice in this scheme.

Starting from the initial separated state of this whole system ∣Ψ(0)⟩_W = ∣0⟩_MW ⨂ ∣0⟩_m ⨂ ∣g⟩_s, since this mechanical mode possesses a large damping rate caused by MPA, its steady state will be obviously damped to the vacuum phonon state. The rest freedom with respect to the NV spin and MW photon will also be engineered into a steady state with the quantum correlation [60]. To illustrate this physical process more intuitively and clearly, we need to carry out the relevant numerical simulation on this dynamical evolution, which is dominated by the master equation ${\rm{d}}\hat{\rho }/{\rm{d}}t\,=\,{\rm{i}}[\hat{\rho },{\hat{H}}_{(\mathrm{ii})}^{S(1)}]\,+\,{\kappa }_{1}{ \mathcal D }[\hat{a}]\hat{\rho }\,+\,{\kappa }_{\mathrm{eff}}^{S}{ \mathcal D }[\hat{b}]\hat{\rho }\,+\,\gamma { \mathcal D }[{\hat{\sigma }}_{z}]\hat{\rho }$. By setting the parameters as λ = − 0.7g, r = 2, κ₁ ∼ γ ∼ 0.01g, and ${\kappa }_{\mathrm{eff}}^{S}\sim 200g$, we plot the corresponding simulation results in figure 3 through solving this master equation numerically. Therein, we use three different variables to characterize this nonclassical correlation between the NV spin and MW photon: the joint variance $V({\hat{X}}_{a},{\hat{\sigma }}_{x})$ with $V{({\hat{X}}_{a},{\hat{\sigma }}_{x})\equiv [\langle {\hat{X}}_{a}-{\hat{\sigma }}_{x}{\rangle }^{2}-\langle ({\hat{X}}_{a}-{\hat{\sigma }}_{x})}^{2}\rangle ]/2$ and ${\hat{X}}_{a}\,=\hat{a}+{\hat{a}}^{\dagger }$ [62, 63], the Von Neumann entropy (VNE) with $S(\hat{\varrho })\equiv -\mathrm{Tr}(\hat{\varrho }{\mathrm{log}}_{1}\hat{\varrho })$ and the reduced density matrices for a subsystem of the NV spin and MW photon $\hat{\varrho }={\mathrm{Tr}}_{\mathrm{phonon}}[\hat{\rho }(t)]$ [64, 65], and the traditional quantum correlation function defined as $C({\hat{X}}_{a},{\hat{\sigma }}_{x})\equiv [\langle {\hat{X}}_{a}\cdot {\hat{\sigma }}_{x}\rangle -\langle {\hat{X}}_{a}\rangle \langle {\hat{\sigma }}_{x}\rangle ]/2$ [61, 66].

Figure 3.

New window|Download| PPT slide
Figure 3.(Color online) Time-dependent quantum correlations of the MW photon and NV spin. The parameters are set as λ = − 0.7g, r = 2, κ₁ ∼ γ ∼ 0.01g, and ${\kappa }_{\mathrm{eff}}^{S}\sim 200g$, where the black solid line means the joint variance $V({\hat{X}}_{a},{\hat{\sigma }}_{x})$, the red dashed line means the VNE $S(\hat{\varrho })$, and the blue short dashed line means the quantum correlation $C({\hat{X}}_{a},{\hat{\sigma }}_{x})$.


The results in figure 3 clearly show that, for a given squeezing parameter r = 2, we can achieve an entangled steady state of this pair of NV spin and MW photon as the whole system is dynamically evolved to the moment of ∼ 10/g. We also note that this steady state is not the maximum entanglement, and all three different characterizations are highly consistent with our analysis above. Therefore, although there is no direct strong interaction between this pair of NV spin and MW cavity, they will be fast-engineered into a steady quantum correlated state with the assistance of this strongly damped MR. This type of application may also provide another efficient route toward the manipulation of a single NV spin with an MW photon at the single-quantum level.

5. Experimental estimations

To examine the feasibility of our scheme in realistic conditions, we now discuss the relevant experimental parameters. We consider a silicon cantilever with dimensions (l = 11, w = 0.05, t = 0.04) $\mu$m; the fundamental frequency and the zero-field fluctuation can be expressed as ${\omega }_{m}\sim 3.516\times (t/{l}^{2})\sqrt{E/12\varrho }\sim 2\pi \times 2.5$ MHz and ${a}_{0}\,=\sqrt{{\hslash }/2m{\omega }_{m}}\sim 5.0\times {10}^{-13}$ m, respectively [67-73], the quality factor of this mechanical mode is estimated as about Q_m ∈ [10³, 10⁶] with Young’s modulus E ∼ 1.3 × 10¹¹ Pa, the mass density ϱ ∼ 2.33 × 10³kg/m³, and the effective mass of this magnet-attached cantilever resonator is estimated as m ∼ 1.4 × 10⁻¹⁷ kg [72, 74, 75]. With the rapid development of advanced nanotechnology, the first-order gradient magnetic field caused by the nanomagnet is expected to reach near G_m ∼ 10⁶ T/m, and we can also obtain that the magnetic coupling strength between the MR and NV spin satisfies g ∈ 2π × [1, 100] kHz [74, 76, 77]. On the other side, the fundamental frequency of this MW cavity with quality more than 10⁷ is also estimated as ω_c ∈ 2π × [1, 20] GHz [3, 4, 18, 19].

6. Conclusion

Addressing a key scientific challenge with respect to realizing direct manipulation of an NV spin at the single-quantum level with an MW cavity, we propose a hybrid system with spin-phonon-photon triple interactions, mediated by a squeezed quantum harmonic resonator. Applying the general MPA to this in-between phonon mode, we can exponentially enhance both direct spin-phonon and photon-phonon couplings, which can even further improve the coherent manipulation of an NV spin and MW photons with higher efficiency. In view of this three-body system with enhanced couplings and adjustable sidebands, this scheme may provide a more efficient phonon-mediated platform to bridge or manipulate the MW quanta and a single electron spin coherently. We also hope it can evoke wider applications in the areas of quantum state transfer and preparation, ultrasensitive detection, QND, etc. We stress that this proposed scheme is general and can be applied to others electron spins [78], such as silicon- [79, 80], germanium- [81], and tin-vacancy centers [82], and its original magnetic spin-phonon coupling mechanism can also be applied in other devices with strain coupling or electric coupling [46-48, 50, 83, 84].

Appendix A. Realizing the second-order nonlinear interaction of this MR

The main mechanism and physical image of this proposal is illustrated in figures 1(a) and (b). Working as an in-middle mediator, the Hamiltonian of this MR with a time-dependent spring constant is described as
(A1)$\begin{eqnarray}{H}_{\mathrm{mr}}=\frac{{\hat{p}}_{z}^{2}}{2m}+\frac{1}{2}{k}_{0}{\hat{z}}^{2}+\frac{1}{2}{k}_{\chi }(t){\hat{z}}^{2},\end{eqnarray}$
where $\hat{z}$ and ${\hat{p}}_{z}$ are the displacement and momentum operators in the z direction, m means the effective mass of this MR, and ${k}_{0}={\omega }_{m}^{2}m$ and ${k}_{\chi }(t)={\rm{\Delta }}k\cos 2{\omega }_{\chi }t$ correspond to the unperturbed fundamental spring constant (the fundamental frequency ω_m) and the time-varying spring constant, respectively. We also state that Δk is the variation of the spring constant. Applying the definition of the displacement operator $\hat{z}={z}_{{zf}}({\hat{b}}^{\dagger }+\hat{b})$ and the zero-field quantum fluctuation ${z}_{{zf}}=\sqrt{{\hslash }/2m{\omega }_{m}}$, we can quantize the MR’s Hamiltonian ${\hat{H}}_{\mathrm{mr}}$ as (ℏ = 1),
(A2)$\begin{eqnarray}{\hat{H}}_{\mathrm{mr}}={\omega }_{m}{\hat{b}}^{\dagger }\hat{b}-{{\rm{\Omega }}}_{\chi }\cos {(2{\omega }_{\chi }t)({\hat{b}}^{\dagger }+\hat{b})}^{2},\end{eqnarray}$
where ${{\rm{\Omega }}}_{\chi }=-{\rm{\Delta }}{{kz}}_{\mathrm{zf}}^{2}/2$ is the second-order nonlinear drive amplitude.

Appendix B. Magnetic coupling of the MR to an NV center in the dressed-state basis

For the second spin-phonon coupling part, this vibrating magnet tip produces the time-dependent gradient magnetic field ${{\boldsymbol{B}}}_{g}(t)={{\boldsymbol{B}}}_{g}\cos {\omega }_{m}t$ and then acts on the NV spin magnetically, with the gradient magnetic field vector components B_g = (B_x, B_y, B_z), and the MR fundamental frequency ω_m. Here, we note ω_m (∼MHz) is much smaller than the NV ground-state transition frequency (D ± δ/2∼ GHz); these far-off resonant interactions between the NV spin and the gradient magnetic fields along the x and y directions are negligible. In the rotating frame at the frequency ω_m, the Hamiltonian for describing the magnetic interaction between the mechanical mode and the single NV center is
(B1)$\begin{eqnarray}{\hat{H}}_{\mathrm{sp}}={\mu }_{B}{g}_{e}{G}_{m}\hat{z}{\hat{S}}_{z}={\lambda }_{0}({\hat{b}}^{\dagger }+\hat{b}){\hat{S}}_{z},\end{eqnarray}$
where λ₀ = $\mu$_Bg_eG_mz_zf is the magnetic coupling strength with the Bohr magneton $\mu$_B = 14 GHz/T, the NV Landé factor g_e ≃ 2, and the first-order magnetic field gradient G_m = ∂B_z/∂z. As illustrated in figure 1(c), in the bare-state basis, we can not achieve the effective spin-phonon transition interaction because of ω_m ≪ D ± δ/2, but we can utilize the dressed-state method to resolve this problem.

As shown in figure 1, to achieve the effective quantum transition interaction between the MR and the NV spin, we apply the dichromatic MW classical fields ${B}_{x}^{\pm }(t)$ (with frequencies ω₊ and ω₋ and polarized in the x direction) to dress the state transitions ∣0⟩ ⇔ ∣ ± 1⟩. The Hamiltonian for describing the single NV center dressed by the dichromatic MW fields is expressed as ${\hat{H}}_{\mathrm{NV}}=D{\hat{S}}_{z}^{2}+\delta {\hat{S}}_{z}/2\,+{\mu }_{B}{g}_{e}[{B}_{x}^{+}(t)+{B}_{x}^{-}(t)]{\hat{S}}_{x}$, with the classical periodic driving fields ${B}_{x}^{\pm }(t)={B}_{0x}^{\pm }\cos ({\omega }_{\pm }t+{\phi }_{\pm })$. Then, we can rewrite its Hamiltonian under the rotation frame with the MW frequencies of ω_±,
(B2)$\begin{eqnarray}{\hat{H}}_{\mathrm{NV}}=\displaystyle \sum _{j=\pm }-{{\rm{\Delta }}}_{j}| j\rangle \langle j| +\displaystyle \frac{{{\rm{\Omega }}}_{j}}{2}(| 0\rangle \langle j| +| j\rangle \langle 0| ),\end{eqnarray}$
where Δ_± ≡ ∣D − ω_± ± δ/2∣ and ${{\rm{\Omega }}}_{\pm }\equiv {g}_{e}{\mu }_{B}{B}_{0x}^{\pm }/\sqrt{2}$. For simplicity, we set Δ_± = Δ and Ω_± = Ω in the following discussions. The Hamiltonian (A1) couples the state ∣0⟩ to a 'bright'state $| b\rangle =(| +1\rangle +| -1\rangle )/\sqrt{2}$, while the 'dark'state $| d\rangle =(| +\rangle -| -\rangle )/\sqrt{2}$ is decoupled. The resulting eigen basis of ${\hat{H}}_{\mathrm{NV}}$ is therefore given by ∣d⟩ and the two dressed states $| g\rangle =\cos \theta | 0\rangle -\sin \theta | b\rangle $ and $| e\rangle =\cos \theta | b\rangle +\sin \theta | 0\rangle $, with the definition $\tan (2\theta )\,=-\sqrt{2}{\rm{\Omega }}/{\rm{\Delta }}$. In this dressed basis, we acquire the eigenfrequencies: ω_d = − Δ, and ${\omega }_{e/g}=(-{\rm{\Delta }}\pm \sqrt{{{\rm{\Delta }}}^{2}+2{{\rm{\Omega }}}^{2}})/\sqrt{2}$. The energy-level diagram of the bare states and the dressed states are both illustrated in figure 1(c). With the purpose of building a efficient spin-phonon transition passage, we stress that its parameters Ω and Δ are tunable here, and then we can obtain the available energy levels to match the frequency ω_m of this MR. Therefore, under the rotating wave approximation, we can rewrite the spin-phonon Hamiltonian from equations (B1)-(B2):
(B3)$\begin{eqnarray}{\hat{H}}_{\mathrm{sp}}^{d}=\displaystyle \sum _{r=e,d}{\omega }_{r}{\hat{\sigma }}_{\mathrm{rr}}+{\lambda }_{d}{\hat{b}}^{\dagger }{\hat{\sigma }}_{\mathrm{dg}}^{-}+{\lambda }_{e}{\hat{b}}^{\dagger }{\hat{\sigma }}_{\mathrm{ed}}^{-}+{\rm{h}}.{\rm{c}}.,\end{eqnarray}$
where the parameters are expressed as ω_e ≡ ω_e − ω_g, ω_d ≡ ω_d − ω_g, ${\hat{\sigma }}_{\mathrm{ee}}\equiv | e\rangle \langle e| $, ${\hat{\sigma }}_{\mathrm{dd}}\equiv | d\rangle \langle d| $, ${\hat{\sigma }}_{\mathrm{gg}}\equiv | g\rangle \langle g| $, ${\hat{\sigma }}_{\mathrm{dg}}^{-}\equiv | g\rangle \langle d| $, ${\hat{\sigma }}_{\mathrm{ed}}^{-}\equiv | d\rangle \langle e| $, ${\lambda }_{d}=-{\lambda }_{0}\sin \theta $, and ${\lambda }_{e}\,={\lambda }_{0}\cos \theta $.

Acknowledgments

Y. Z. is supported by the National Natural Science Foundation of China under Grants No. 11 774 285, No. 11 774 282, and No. 11 504 102, the Natural Science Foundation of Hubei Province under Grants No. 2020CFB748, and No. 2019CFB788, the Research Project of Hubei Education Department under Grants No. D20201803, No. B2020079, and No. B2020078, the Doctoral Scientific Research Foundation of Hubei University of Automotive Technology (HUAT) under Grants No. BK201906, and No. BK202008, and the Foundation of Discipline Innovation Team of HUAT. H. T. R. is supported by the the Doctoral Scientific Research Foundation of Liaocheng University under Grant No. 318 052 054. Parts of the simulations are coded in python using the qutip library [64, 65].

Reference View Option By original order
By published year
By cited within times
By Impact factor

[1]Georgescu I M, Ashhab S, Nori Franco 2014Quantum simulation
Rev. Mod. Phys. 86 153 185
DOI:10.1103/RevModPhys.86.153 [Cited within: 1]

[2]Xiang Z-L, Ashhab S, You J Q, Nori F 2013Hybrid quantum circuits: Superconducting circuits interacting with other quantum systems
Rev. Mod. Phys. 85 623 653
DOI:10.1103/RevModPhys.85.623

[3]Carusotto I, Houck A A, Kollár A J, Roushan P, Schuster D I, Simon J 1919Photonic materials in circuit quantum electrodynamics
Nat. Phys. 16 268 279
DOI:10.1038/s41567-020-0815-y [Cited within: 2]

[4]Blais A, Girvin S M, Oliver W D 1919Quantum information processing and quantum optics with circuit quantum electrodynamics
Nat. Phys. 16 247 256
DOI:10.1038/s41567-020-0806-z [Cited within: 2]

[5]Ladd T D, Jelezko F, Laflamme R, Nakamura Y, Monroe C, O’Brien J L 2010Quantum computers
Nature 464 45
DOI:10.1038/nature08812

[6]Ren W et al. 1919Simultaneous excitation of two noninteracting atoms with time-frequency correlated photon pairs in a superconducting circuit
Phys. Rev. Lett. 125 133601
DOI:10.1103/PhysRevLett.125.133601

[7]Liu Y-X, You J Q, Wei L F, Sun C P, Nori F 2005Optical selection rules and phase-dependent adiabatic state control in a superconducting quantum circuit
Phys. Rev. Lett. 95 087001
DOI:10.1103/PhysRevLett.95.087001 [Cited within: 1]

[8]Buluta I, Ashhab S, Nori F 2011Natural and artificial atoms for quantum computation
Rep. Prog. Phys. 74 104401
DOI:10.1088/0034-4885/74/10/104401 [Cited within: 1]

[9]Lecocq F, Ranzani L, Peterson G A, Cicak K, Jin X Y, Simmonds R W, Teufel J D, Aumentado J 2021Efficient qubit measurement with a nonreciprocal microwave amplifier
Phys. Rev. Lett. 126 020502
DOI:10.1103/PhysRevLett.126.020502 [Cited within: 1]

[10]Li P-B, Liu Y-C, Gao S-Y, Xiang Z-L, Rabl P, Xiao Y-F, Li F-L 2015Hybrid quantum device based on NV centers in diamond nanomechanical resonators plus superconducting waveguide cavities
Phys. Rev. Appl. 4 044003
DOI:10.1103/PhysRevApplied.4.044003 [Cited within: 1]

[11]Romero G, García-Ripoll J J, Solano E 2009Microwave photon detector in circuit QED
Phys. Rev. Lett. 102 173602
DOI:10.1103/PhysRevLett.102.173602 [Cited within: 1]

[12]Song G-Z, Kwek L-C, Deng F-G, Long G-L 2019Microwave transmission through an artificial atomic chain coupled to a superconducting photonic crystal
Phys. Rev. A 99 043830
DOI:10.1103/PhysRevA.99.043830 [Cited within: 1]

[13]Song C et al. 201710-qubit entanglement and parallel logic operations with a superconducting circuit
Phys. Rev. Lett. 119 180511
DOI:10.1103/PhysRevLett.119.180511 [Cited within: 1]

[14]Li X et al. 2018Perfect quantum state transfer in a superconducting qubit chain with parametrically tunable couplings
Phys. Rev. Appl. 10 054009
DOI:10.1103/PhysRevApplied.10.054009 [Cited within: 1]

[15]Clerk A A, Devoret M H, Girvin S M, Florian Marquardt, Schoelkopf R J 2010Introduction to quantum noise, measurement, and amplification
Rev. Mod. Phys. 82 1155 1208
DOI:10.1103/RevModPhys.82.1155 [Cited within: 1]

[16]Degen C L, Reinhard F, Cappellaro P 2017Quantum sensing
Rev. Mod. Phys. 89 035002
DOI:10.1103/RevModPhys.89.035002 [Cited within: 1]

[17]Forn-Díaz P, Lamata L, Rico E, Kono J, Solano E 2019Ultrastrong coupling regimes of light-matter interaction
Rev. Mod. Phys. 91 025005
DOI:10.1103/RevModPhys.91.025005 [Cited within: 1]

[18]Clerk A A, Lehnert K W, Bertet P, Petta J R, Nakamura Y 2020Hybrid quantum systems with circuit quantum electrodynamics
Nat. Phys. 16 257 267
DOI:10.1038/s41567-020-0797-9 [Cited within: 2]

[19]Haroche S, Brune M, Raimond J M 2020From cavity to circuit quantum electrodynamics
Nat. Phys. 16 243 246
DOI:10.1038/s41567-020-0812-1 [Cited within: 2]

[20]Doherty M W, Manson N B, Delaney P, Jelezko F, Wrachtrup J, Hollenberg L C L 2013The nitrogen-vacancy colour centre in diamond
Phys. Rep. 528 1 45
DOI:10.1016/j.physrep.2013.02.001 [Cited within: 1]

[21]Cai J, Retzker A, Jelezko F, Plenio M B 2013A large-scale quantum simulator on a diamond surface at room temperature
Nat. Phys. 9 168
DOI:10.1038/nphys2519

[22]Casola F, van der Sar T, Yacoby A 2018Probing condensed matter physics with magnetometry based on nitrogen-vacancy centres in diamond
Nat. Rev. Mater. 3 17088
DOI:10.1038/natrevmats.2017.88

[23]Guo H et al. 1919NV center pumped and enhanced by nanowire ring resonator laser to integrate a 10 $\mu$m-scale spin-based sensor structure
Nanotechnology 32 055502
DOI:10.1088/1361-6528/abc20b [Cited within: 1]

[24]Jin L, Pfender M, Aslam N, Neumann P, Yang S, Wrachtrup J, Liu R-B 2015Proposal for a room-temperature diamond maser
Nat. Commun. 6 8251
DOI:10.1038/ncomms9251 [Cited within: 1]

[25]Breeze J D, Salvadori E, Sathian J, Alford N M, Kay C W M 2018Continuous-wave room-temperature diamond maser
Nature 555 493 496
DOI:10.1038/nature25970 [Cited within: 1]

[26]Zhu X et al. 2011Coherent coupling of a superconducting flux qubit to an electron spin ensemble in diamond
Nature 478 221
DOI:10.1038/nature10462 [Cited within: 1]

[27]Kubo Y et al. 2011Hybrid quantum circuit with a superconducting qubit coupled to a spin ensemble
Phys. Rev. Lett. 107 220501
DOI:10.1103/PhysRevLett.107.220501 [Cited within: 1]

[28]Kubo Y et al. 2010Strong coupling of a spin ensemble to a superconducting resonator
Phys. Rev. Lett. 105 140502
DOI:10.1103/PhysRevLett.105.140502 [Cited within: 1]

[29]Li P-B, Zhou Y, Gao W-B, Nori F 1919Enhancing spin-phonon and spin-spin interactions using linear resources in a hybrid quantum system
Phys. Rev. Lett. 125 153602
DOI:10.1103/PhysRevLett.125.153602 [Cited within: 6]

[30]Rabl P, Cappellaro P, Gurudev Dutt M V, Jiang L, Maze J R, Lukin M D 2009Strong magnetic coupling between an electronic spin qubit and a mechanical resonator
Phys. Rev. B 79 041302
DOI:10.1103/PhysRevB.79.041302

[31]Li P-B, Xiang Z-L, Rabl P, Nori F 2016aHybrid quantum device with nitrogen-vacancy centers in diamond coupled to carbon nanotubes
Phys. Rev. Lett. 117 015502
DOI:10.1103/PhysRevLett.117.015502

[32]Li P-B, Li H-R, Li F-L 2016bEnhanced electromechanical coupling of a nanomechanical resonator to coupled superconducting cavities
Sci. Rep. 6 19065
DOI:10.1038/srep19065 [Cited within: 4]

[33]Lü X-Y, Wu Y, Johansson J R, Jing H, Zhang J, Nori F 2015Squeezed optomechanics with phase-matched amplification and dissipation
Phys. Rev. Lett. 114 093602
DOI:10.1103/PhysRevLett.114.093602 [Cited within: 1]

[34]Nation P D, Johansson J R, Blencowe M P, Nori F 2012Colloquium: Stimulating uncertainty: Amplifying the quantum vacuum with superconducting circuits
Rev. Mod. Phys. 84 1 24
DOI:10.1103/RevModPhys.84.1

[35]Muñoz C S, Lara A, Puebla J, Nori F 2018Hybrid systems for the generation of nonclassical mechanical states via quadratic interactions
Phys. Rev. Lett. 121 123604
DOI:10.1103/PhysRevLett.121.123604

[36]Qin W, Miranowicz A, Li P-B, Lü X-Y, You J Q, Nori F 2018Exponentially enhanced light-matter interaction, cooperativities, and steady-state entanglement using parametric amplification
Phys. Rev. Lett. 120 093601
DOI:10.1103/PhysRevLett.120.093601

[37]Leroux C, Govia L C G, Clerk A A 2018Enhancing cavity quantum electrodynamics via antisqueezing: synthetic ultrastrong coupling
Phys. Rev. Lett. 120 093602
DOI:10.1103/PhysRevLett.120.093602

[38]Lemonde M-A, Didier N, Clerk A A 2016Enhanced nonlinear interactions in quantum optomechanics via mechanical amplification
Nat. Commun. 7 11338
DOI:10.1038/ncomms11338 [Cited within: 4]

[39]Blum S, O’Brien C, Lauk N, Bushev P, Fleischhauer M, Morigi G 2015Interfacing microwave qubits and optical photons via spin ensembles
Phys. Rev. A 91 033834
DOI:10.1103/PhysRevA.91.033834 [Cited within: 2]

[40]Li B, Li P-B, Zhou Y, Ma S-L, Li F-L 2017Quantum microwave-optical interface with nitrogen-vacancy centers in diamond
Phys. Rev. A 96 032342
DOI:10.1103/PhysRevA.96.032342 [Cited within: 2]

[41]Taylor J M, Cappellaro P, Childress L, Jiang L, Budker D, Hemmer P R, Yacoby A, Walsworth R, Lukin M D 2008High-sensitivity diamond magnetometer with nanoscale resolution
Nat. Phys. 4 810 816
DOI:10.1038/nphys1075 [Cited within: 1]

[42]Tang Y-C, Li Y-S, Hao L, Hou S-Y, Long G L 2012Quantum-nondemolition determination of an unknown Werner state
Phys. Rev. A 85 022329
DOI:10.1103/PhysRevA.85.022329 [Cited within: 1]

[43]Szorkovszky A, Clerk A A, Doherty A C, Bowen W P 2014Mechanical entanglement via detuned parametric amplification
New J. Phys. 16 063043
DOI:10.1088/1367-2630/16/6/063043 [Cited within: 1]

[44]Szorkovszky A, Doherty A C, Harris G I, Bowen W P 2011Mechanical squeezing via parametric amplification and weak measurement
Phys. Rev. Lett. 107 213603
DOI:10.1103/PhysRevLett.107.213603

[45]Rugar D, Grütter P 1991Mechanical parametric amplification and thermomechanical noise squeezing
Phys. Rev. Lett. 67 699
DOI:10.1103/PhysRevLett.67.699 [Cited within: 1]

[46]Teissier J, Barfuss A, Appel P, Neu E, Maletinsky P 2014Strain coupling of a nitrogen-vacancy center spin to a diamond mechanical oscillator
Phys. Rev. Lett. 113 020503
DOI:10.1103/PhysRevLett.113.020503 [Cited within: 2]

[47]Meesala S, Sohn Y-I, Atikian H A, Kim S, Burek M J, Choy J T, Lončar M 2016Enhanced strain coupling of nitrogen-vacancy spins to nanoscale diamond cantilevers
Phys. Rev. Appl. 5 034010
DOI:10.1103/PhysRevApplied.5.034010

[48]MacQuarrie E R, Otten M, Gray S K, Fuchs G D 2017Cooling a mechanical resonator with nitrogen-vacancy centres using a room temperature excited state spin-strain interaction
Nat. Commun. 8 14358
DOI:10.1038/ncomms14358 [Cited within: 1]

[49]Cai J, Jelezko F, Plenio M B 2014Hybrid sensors based on colour centres in diamond and piezoactive layers
Nat. Commun. 5 4065
DOI:10.1038/ncomms5065

[50]Ovartchaiyapong P, Lee K W, Myers B A, Jayich A C B 2014Dynamic strain-mediated coupling of a single diamond spin to a mechanical resonator
Nat. Commun. 5 4429
DOI:10.1038/ncomms5429 [Cited within: 2]

[51]Aspelmeyer M, Kippenberg T J, Marquardt F 2014Cavity optomechanics
Rev. Mod. Phys. 86 1391
DOI:10.1103/RevModPhys.86.1391 [Cited within: 3]

[52]Kronwald A, Marquardt F, Clerk A A 2013Arbitrarily large steady-state bosonic squeezing via dissipation
Phys. Rev. A 88 063833
DOI:10.1103/PhysRevA.88.063833 [Cited within: 3]

[53]Kienzler D, Lo H-Y, Keitch B, de Clercq L, Leupold F, Lindenfelser F, Marinelli M, Negnevitsky V, Home J P 2015Quantum harmonic oscillator state synthesis by reservoir engineering
Science 347 53
DOI:10.1126/science.1261033

[54]Wollman E E, Lei C U, Weinstein A J, Suh J, Kronwald A, Marquardt F, Clerk A A, Schwab K C 2015Quantum squeezing of motion in a mechanical resonator
Science 349 952
DOI:10.1126/science.aac5138

[55]Pirkkalainen J-M, Damskägg E, Brandt M, Massel F, Sillanpää M A 2015Squeezing of quantum noise of motion in a micromechanical resonator
Phys. Rev. Lett. 115 243601
DOI:10.1103/PhysRevLett.115.243601

[56]Cirac J I, Parkins A S, Blatt R, Zoller P 1993Dark squeezed states of the motion of a trapped ion
Phys. Rev. Lett. 70 556
DOI:10.1103/PhysRevLett.70.556 [Cited within: 1]

[57]Ma S-L, Li X-K, Liu X-Y, Xie J-K, Li F-L 2019Stabilizing Bell states of two separated superconducting qubits via quantum reservoir engineering
Phys. Rev. A 99 042336
DOI:10.1103/PhysRevA.99.042336 [Cited within: 2]

[58]Lü D-Y, Wang G H, Zhou Y, Xu L, Wang Q L 2021Collective decay induce quantum phase transition in a well-controlled hybrid quantum system
Results Phys. 21 103832
DOI:10.1016/j.rinp.2021.103832

[59]Zhou Y, Li B, Li X-X, Li F-L, Li P-B 2018Preparing multiparticle entangled states of nitrogen-vacancy centers via adiabatic ground-state transitions
Phys. Rev. A 98 052346
DOI:10.1103/PhysRevA.98.052346

[60]Li P 2008Generation of two-mode field squeezing through selective dynamics in cavity QED
Phys. Rev. A 77 015809
DOI:10.1103/PhysRevA.77.015809 [Cited within: 3]

[61]Zhou T, Cui J, Long G L 2011Measure of nonclassical correlation in coherence-vector representation
Phys. Rev. A 84 062105
DOI:10.1103/PhysRevA.84.062105 [Cited within: 1]

[62]Ma S-L, Li Z, Fang A-P, Li P-B, Gao S-Y, Li F-L 2014Controllable generation of two-mode-entangled states in two-resonator circuit QED with a single gap-tunable superconducting qubit
Phys. Rev. A 90 062342
DOI:10.1103/PhysRevA.90.062342 [Cited within: 1]

[63]Li P-B, Gao S-Y, Li F-L 2012Engineering two-mode entangled states between two superconducting resonators by dissipation
Phys. Rev. A 86 012318
DOI:10.1103/PhysRevA.86.012318 [Cited within: 2]

[64]Johansson J R, Nation P D, Nori F 2012Qutip: An open-source Python framework for the dynamics of open quantum systems
Comput. Phys. Commun. 183 1760 1772
DOI:10.1016/j.cpc.2012.02.021 [Cited within: 2]

[65]Johansson J R, Nation P D, Nori F 2013Qutip 2: A Python framework for the dynamics of open quantum systems
Comput. Phys. Commun. 184 1234 1240
DOI:10.1016/j.cpc.2012.11.019 [Cited within: 2]

[66]Zeng B, Chen X, Zhou D-L, Wen X-G 2015Quantum information meets quantum matter—from quantum entanglement to topological phase in many-body systems
arXiv:1508.02595
[Cited within: 1]

[67]Sidles J A, Garbini J L, Bruland K J, Rugar D, Züger O, Hoen S, Yannoni C S 1995Magnetic resonance force microscopy
Rev. Mod. Phys. 67 249 265
DOI:10.1103/RevModPhys.67.249 [Cited within: 1]

[68]Imboden M, Mohanty P 2014Dissipation in nanoelectromechanical systems
Phys. Rep. 534 89 146
DOI:10.1016/j.physrep.2013.09.003

[69]Atatüre M, Englund D, Vamivakas N, Lee S-Y, Wrachtrup J 2018Scalable architecture for a room temperature solid-state quantum information processor
Nat. Rev. Mater. 3 38
DOI:10.1038/s41578-018-0008-9

[70]Li M, Tang H X, Roukes M L 2007Ultra-sensitive NEMS-based cantilevers for sensing, scanned probe and very high-frequency applications
Nat. Nanotechnol. 2 114
DOI:10.1038/nnano.2006.208

[71]Ekinci K L, Roukes M L 2005Nanoelectromechanical systems
Rev. Sci. Instrum. 76 061101
DOI:10.1063/1.1927327

[72]Yang J, Ono T, Esashi M 2000Surface effects and high quality factors in ultrathin single-crystal silicon cantilevers
Appl. Phys. Lett. 77 3860 3862
DOI:10.1063/1.1330225 [Cited within: 1]

[73]Brantley W A 1973Calculated elastic constants for stress problems associated with semiconductor devices
J. Appl. Phys. 44 534 535
DOI:10.1063/1.1661935 [Cited within: 1]

[74]Tao Y, Eichler A, Holzherr T, Degen C L 2016Ultrasensitive mechanical detection of magnetic moment using a commercial disk drive write head
Nat. Commun. 7 12714
DOI:10.1038/ncomms12714 [Cited within: 2]

[75]Tao Y, Boss J M, Moores B A, Degen C L 2014Single-crystal diamond nanomechanical resonators with quality factors exceeding one million
Nat. Commun. 5 3638
DOI:10.1038/ncomms4638 [Cited within: 1]

[76]Meyer E, Rast S 2007Magnetic tips probe the nanoworld
Nat. Nanotechnol. 2 267 268
DOI:10.1038/nnano.2007.115 [Cited within: 1]

[77]Mamin H J, Poggio M, Degen C L, Rugar D 2007Nuclear magnetic resonance imaging with 90-nm resolution
Nat. Nanotechnol. 2 301 306
DOI:10.1038/nnano.2007.105 [Cited within: 1]

[78]Thiering G, Gali A 2018Ab initio magneto-optical spectrum of group-IV vacancy color centers in diamond
Phys. Rev. X 8 021063
DOI:10.1103/PhysRevX.8.021063 [Cited within: 1]

[79]Hepp C et al. 2014Electronic structure of the silicon vacancy color center in diamond
Phys. Rev. Lett. 112 036405
DOI:10.1103/PhysRevLett.112.036405 [Cited within: 1]

[80]Sipahigil A, Jahnke K D, Rogers L J, Teraji T, Isoya J, Zibrov A S, Jelezko F, Lukin M D 2014Indistinguishable photons from separated silicon-vacancy centers in diamond
Phys. Rev. Lett. 113 113602
DOI:10.1103/PhysRevLett.113.113602 [Cited within: 1]

[81]Bhaskar M K et al. 2017Quantum nonlinear optics with a germanium-vacancy color center in a nanoscale diamond waveguide
Phys. Rev. Lett. 118 223603
DOI:10.1103/PhysRevLett.118.223603 [Cited within: 1]

[82]Iwasaki T, Miyamoto Y, Taniguchi T, Siyushev P, Metsch M H, Jelezko F, Hatano M 2017Tin-vacancy quantum emitters in diamond
Phys. Rev. Lett. 119 253601
DOI:10.1103/PhysRevLett.119.253601 [Cited within: 1]

[83]Rabl P, Kolkowitz S J, Koppens F H L, Harris J G E, Zoller P, Lukin M D 2010A quantum spin transducer based on nanoelectromechanical resonator arrays
Nat. Phys. 6 602 608
DOI:10.1038/nphys1679 [Cited within: 1]

[84]Dolde F, Fedder H, Doherty M W, Nöbauer T, Rempp F, Balasubramanian G, Wolf T, Reinhard F, Hollenberg L C L, Jelezko F 2011Electric-field sensing using single diamond spins
Nat. Phys. 7 459 463
DOI:10.1038/nphys1969 [Cited within: 1]