Next Article in Journal
Mesoscopic States of Light for the Detection of Weakly Absorbing Objects
Next Article in Special Issue
Method for the Quantum Metric Tensor Measurement in a Continuous Variable System
Previous Article in Journal
Manipulating Transverse Spin Angular Momentum with Symmetrically Modulated Hybridly Polarized Vector Optical Field
Previous Article in Special Issue
Integer and Fractional Floquet Resonances in a Driven Three-Well System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phase-Controlled Entanglement in a Four-Mode Optomechanical System

School of Physics and Electronic Electrical Engineering, Huaiyin Normal University, Huai’an 223300, China
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(11), 818; https://doi.org/10.3390/photonics9110818
Submission received: 10 October 2022 / Revised: 24 October 2022 / Accepted: 26 October 2022 / Published: 29 October 2022

Abstract

:
We present a scheme for realizing phase-controlled entanglement in a microwave optomechanical system comprising two microwave cavities and two mechanical oscillators. Under specific driving conditions, we show that this optomechanical interface can be exploited to generate simultaneously the stationary cavity–cavity entanglement, mechanical–mechanical entanglement, and cavity–mechanical entanglement. Due to the closed loop interaction, we find that the entanglement can be controlled flexibly by tuning the phase difference between the optomechanical coupling strengths. The dependence of the entanglement on the amplitudes of the optomechanical coupling strengths is also explored in detail. Moreover, the bipartite entanglements are robust against temperature, and it is shown that the mechanical oscillators are cooled to the ground state in the parameter regimes for observing entanglement.

1. Introduction

Entanglement, one essential feature of quantum mechanics, is an indispensable resource for quantum information processing. This phenomenon has been demonstrated in various systems, such as superconducting qubits [1], atomic ensembles [2,3], individual atoms [4] and ions [5], and electron spins [6]. In the past few decades, mechanical oscillators with high resonance frequencies and quality factors have been extensively explored. They can be coupled to a variety of quantum systems, including superconducting qubits [7], quantum dots [8], a single-spin qubit [9], Bose–Einstein condensate [10], and to photons at nearly arbitrary wavelength. Recently, entanglement between mechanical oscillator and atomic spin ensemble [11] as well as the phonon-mediated entanglement between two superconducting qubits [12] have been demonstrated. Therefore, hybrid mechanical systems provide a promising platform for building quantum networks [13].
Cavity optomechanical systems, where the mechanical oscillator interacts with the electromagnetic cavity via radiation pressure, can also serve as an excellent candidate for the interface that connects different components in quantum networks [14,15,16,17,18,19,20,21,22]. Recent milestones in this field include cooling of the mechanical oscillator to the quantum ground state [23,24,25] and observation of quantum effects such as optomechanical squeezing of light and mechanical motion [26,27,28,29]. Moreover, a myriad of schemes have been proposed to generate the entanglement based on optomechanical systems [30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48]. Recent optomechanical experiments have successfully demonstrated the entanglement between photons and a mechanical oscillator [49], between two mechanical oscillators [50,51,52,53], and between photons at both optical [54] and microwave [55] frequencies. Very recently, significant attention has been paid to the nonreciprocal transmission in multimode optomechanical systems with closed loop interaction [56,57,58,59,60,61,62,63,64,65,66,67,68,69,70], where the phase difference between the optomechanical coupling strengths can be exploited to control the direction of light propagation. In addition, it has been shown that the light-vibration entanglement can be controlled by the phase difference in three-mode optomechanical systems [71,72].
Here, we study the phase-controlled entanglement in a four-mode optomechanical system, which consists of two cavity modes coupled via two mechanical oscillators [64,65,66,67,68]. To generate the entanglement quantified by logarithmic negativity, the cavities are driven by four tones, including two blue tones and two red tones. By tuning the phase difference and the optomechanical coupling strengths, we show that this hybrid system provides a platform to generate simultaneously the cavity–cavity entanglement, mechanical–mechanical entanglement, and cavity–mechanical entanglement. Furthermore, we find that the stationary entanglement is robust against the environmental temperature and both the mechanical oscillators can be cooled to the ground state. Our scheme could be realized based on the current experiments [64,65,66,67].

2. Model and Theory

We consider the optomechanical system schematically shown in Figure 1a, where two cavity modes are, respectively, coupled to two mechanical oscillators via radiation pressure. The Hamiltonian of the system is given by
H / = i = 1 2 ω i a i a i + i = 1 2 Ω i b i b i + i = 1 2 j = 1 2 g 0 , i j a i a i ( b j + b j ) + i = 1 2 j = 1 2 ε i j [ a i e i ( ω d , i j t + ϕ i j ) + a i e i ( ω d , i j t + ϕ i j ) ] ,
where the first term describes the energy of the cavity modes with creation (annihilation) operators a ( a ) and resonance frequency ω i ; the second term represents the mechanical modes with creation (annihilation) operators b ( b ) and resonance frequency Ω i ; the third term denotes the optomechanical coupling between the cavity modes a i and the mechanical modes b j with g 0 , i j being the vacuum optomechanical coupling strength; and the last term corresponds to the interaction between the cavity modes and the driving fields with amplitude ε i j , frequency ω d , i j , and phase ϕ i j . We linearize the Hamiltonian by expanding cavity modes as the sum of their steady-state field amplitudes and the small fluctuations, i.e., a i a i + j = 1 2 α i j e i ω d , i j t , where α i j is the coherent state amplitude produced in cavity i due to the drive with detuning Δ i j = ω i ω d , i j . Here, we choose Δ 11 = Ω 1 + δ 1 , Δ 12 = ( Ω 2 + δ 2 ) , Δ 21 = ( Ω 1 + δ 1 ) , and Δ 22 = Ω 2 + δ 2 with δ i being the extra detuning from the exact red (blue) sidebands to study the entanglement in this system, as shown in Figure 1b. Moving into the rotating frame with respect to H 0 / = i = 1 2 ω i a i a i + ( Ω i + δ i ) b i b i and neglecting the counter-rotating and off-resonant terms under the condition of Ω i , | Ω 1 Ω 2 | > { κ i , G 11 , G 22 , J 12 , J 21 } , the linearized Hamiltonian can be derived as (see details in Appendix A)
H / = i = 1 2 δ i b i b i + ( G 11 e i ϕ 11 a 1 b 1 + J 12 e i ϕ 12 a 1 b 2 + J 21 e i ϕ 21 a 2 b 1 + G 22 e i ϕ 22 a 2 b 2 + H . c ) ,
where G 11 = g 0 , 11 | α 11 | , G 22 = g 0 , 22 | α 22 | , J 12 = g 0 , 12 | α 12 | , and J 21 = g 0 , 21 | α 21 | are the field-enhanced optomechanical coupling strengths. Here, we denote G i j ( J i j ) as the red (blue) sidebands of cavity a i on the mechanical mode b j . Due to the loop interaction in this optomechanical system, only the phase difference between the optomechanical coupling strengths has physical effects. By redefining the operators a 1 a 1 e i ϕ 11 , b 2 b 2 e i ( ϕ 12 ϕ 11 ) , a 2 a 2 e i ϕ 21 , the Hamiltonian Equation (2) becomes
H / = δ 1 b 1 b 1 δ 2 b 2 b 2 + G 11 ( a 1 b 1 + a 1 b 1 ) + J 12 ( a 1 b 2 + a 1 b 2 ) + J 21 ( a 2 b 1 + a 2 b 1 ) + G 22 ( a 2 b 2 e i ϕ + a 2 b 2 e i ϕ ) ,
with ϕ = ϕ 11 + ϕ 21 ϕ 12 ϕ 22 being the phase difference.
According to the linearized Hamiltonian (3) and Heisenberg equation of motion, the quantum Langevin equations (QLEs) can be derived by adding the damping and input noise terms phenomenologically [64,65,67,68], which yields
a 1 ˙ = κ 1 2 a 1 i G 11 b 1 i J 12 b 2 + κ 1 a 1 , in , a 2 ˙ = κ 2 2 a 2 i J 21 b 1 i G 22 b 2 e i ϕ + κ 2 a 2 , in , b 1 ˙ = γ 1 2 i δ 1 b 1 i G 11 a 1 i J 21 a 2 + γ 1 b 1 , in , b 2 ˙ = γ 2 2 i δ 2 b 2 i J 12 a 1 i G 22 a 2 e i ϕ + γ 2 b 2 , in ,
where κ i ( i = 1 , 2 ) is the decay rate of cavity mode a i with zero-mean noise operator a i , in , and γ i is the damping rate of mechanical mode b i with zero-mean operator b i , in . Defining the quadratures U i = ( a i + a i ) / 2 , V i = ( a i a i ) / i 2 , X i = ( b i + b i ) / 2 , P i = ( b i b i ) / i 2 , and the corresponding input noise operators U i , in = ( a i , in + a i , in ) / 2 , V i , in = ( a i , in a i , in ) / i 2 , X i , in = ( b i , in + b i , in ) / 2 , P i , in = ( b i , in b i , in ) / i 2 , the linearized QLEs can be written in the following compact matrix form:
μ ˙ ( t ) = A μ ( t ) + n ( t ) ,
where μ ( t ) = ( U 1 ( t ) , V 1 ( t ) , U 2 ( t ) , V 2 ( t ) , X 1 ( t ) , P 1 ( t ) , X 2 ( t ) , P 2 ( t ) ) T is the vector of the fluctuation operators, n ( t ) = ( κ 1 U 1 , in ( t ) , κ 1 V 1 , in ( t ) , κ 2 U 2 , in ( t ) , κ 2 V 2 , in ( t ) , γ 1 X 1 , in ( t ) , γ 1 P 1 , in ( t ) , γ 2 X 2 , in ( t ) , γ 2 P 2 , in ( t ) ) T is the vector of the corresponding noises, and the coefficient matrix
A = κ 1 / 2 0 0 0 0 G 11 0 J 12 0 κ 1 / 2 0 0 G 11 0 J 12 0 0 0 κ 2 / 2 0 0 J 21 G 22 sin ϕ G 22 cos ϕ 0 0 0 κ 2 / 2 J 21 0 G 22 cos ϕ G 22 sin ϕ 0 G 11 0 J 21 γ 1 / 2 δ 1 0 0 G 11 0 J 21 0 δ 1 γ 1 / 2 0 0 0 J 12 G 22 sin ϕ G 22 cos ϕ 0 0 γ 2 / 2 δ 2 J 12 0 G 22 cos ϕ G 22 sin ϕ 0 0 δ 2 γ 2 / 2 .
The formal solution of Equation (5) is
μ ( t ) = M ( t ) μ ( 0 ) + 0 t d s M ( s ) n ( t s ) ,
where
M ( t ) = exp { A t } .
The system is stable when all of the eigenvalues of A have negative real parts, which can be derived from the Routh–Hurwitz criterion [73]. When the stability conditions are fulfilled, one can obtain M ( ) = 0 in the steady state.
The steady state of the system is a zero-mean Gaussian state due to the linearized dynamics and the Gaussian nature of the input quantum noises, which is fully characterized by an 8 × 8 correlation matrix (CM) V, with its components
V i j = [ μ i ( ) μ j ( ) + μ j ( ) μ i ( ) ] / 2 .
When the system is stable, by substituting Equation (7) into Equation (9), we can obtain
V i j = k , l 0 d s 0 d s M i k ( s ) M j l ( s ) Φ k l ( s s ) ,
where Φ k l ( s s ) = ( n k ( s ) n k ( s ) + n l ( s ) n k ( s ) ) / 2 = D k l δ ( s s ) . According to the correlation functions that the noise operators satisfy,
a i , in ( s ) a i , in ( s ) = δ ( s s ) , a i , in ( s ) a i , in ( s ) = 0 , b i , in ( s ) b i , in ( s ) = ( n m , i + 1 ) δ ( s s ) , b i , in ( s ) b i , in ( s ) = n m , i δ ( s s ) ,
where n m , i = { exp [ Ω i / ( k B T ) ] 1 } 1 denotes the thermal phonon number of the mechanical mode b i with k B being the Boltzmann constant and T being the environment temperature, we can obtain D = Diag [ κ 1 / 2 , κ 1 / 2 , κ 2 / 2 , κ 2 / 2 , ( 2 n m , 1 + 1 ) γ 1 / 2 , ( 2 n m , 1 + 1 ) γ 1 / 2 , ( 2 n m , 2 + 1 ) γ 2 / 2 , ( 2 n m , 2 + 1 ) γ 2 / 2 ] , and Equation (10) becomes V = 0 d s M ( s ) D M ( s ) T . The time evolution of the correlation matrix V satisfies V ˙ ( t ) = A V ( t ) + V ( t ) A T + D [36]. In the steady state, the correlation matrix V fulfills the Lyapunov equation [30]:
A V + V A T = D ,
which is a linear equation for V and can be straightforwardly solved for given parameters in matrix A and matrix D; but the general analytical expression is too cumbersome to be reported here. The bipartite entanglement between different modes can be quantified by the logarithmic negativity [30,31,74]
E N = max [ 0 , ln 2 η ] ,
where η 2 1 / 2 { Σ ( V bp ) [ Σ ( V bp ) 2 4 det V bp ] 1 / 2 } 1 / 2 , with Σ ( V bp ) det B + det B 2 det C . The reduced CM V bp contains the entries of V, and it can be obtained by selecting the rows and columns of the interesting mode in V. By writing the reduced CM V bp in a 2 × 2 block form, we have
V bp B C C T B .
Note that matrix A in Equation (6) contains the phase difference ϕ ; hence, the correlation matrix V and the logarithmic negativity E N depend on ϕ . We will study how to control the entanglement quantified by Equation (13) in this optomechanical system by tuning the phase difference ϕ .
Furthermore, based on the correlation matrix V, the final average phonon numbers in the mechanical modes b 1 and b 2 can be calculated by
b 1 b 1 = n 1 f = V 5 , 5 + V 6 , 6 1 2 , b 2 b 2 = n 2 f = V 7 , 7 + V 8 , 8 1 2 ,
which affect the quantum entanglement in the hybrid system.

3. Results and Discussion

In this section, we choose the parameters from a recent experiment to demonstrate the numerical results about the phase-controlled entanglement in this system. The parameters are [65] Ω 1 / 2 π = 4.34 MHz, Ω 2 / 2 π = 5.64 MHz, κ 1 / 2 π = 2.42 MHz, κ 2 / 2 π = 1.98 MHz, γ 1 / 2 π = 40 Hz, γ 2 / 2 π = 80 Hz, and T = 20 mK. The system works in the resolved sideband regime and the rotating-wave approximation is valid.
In Figure 2, we plot the logarithmic negativity E N between (a) the two cavity modes as well as (b) the two mechanical modes versus the optomechanical coupling strength G 22 / 2 π and the phase difference ϕ / π . At G 22 = 0 , the optomechanical interactions cannot form a closed loop and the phase-dependent effect is absent. With the increase in the coupling strength G 22 , it is evident that the quantum entanglement depends on the phase difference ϕ . Note that cavity mode a 1 is driven on the blue (red) sideband of the mechanical mode b 2 ( b 1 ), while cavity mode a 2 is driven on the blue (red) sideband of the mechanical mode b 1 ( b 2 ). Therefore, the parametric-amplifier interaction J 12 ( a 1 b 2 + a 1 b 2 ) entangles a 1 and b 2 , and the interaction J 21 ( a 2 b 1 + a 2 b 1 ) entangles a 2 and b 1 . Moreover, the beam-splitter interaction G 11 ( a 1 b 1 + a 1 b 1 ) swaps a 1 and b 1 , and the interaction G 22 ( a 2 b 2 e i ϕ + a 2 b 2 e i ϕ ) swaps a 2 and b 2 . Consequently, the two cavity modes can be entangled via the routes a 1 b 1 a 2 and a 1 b 2 a 2 , and the two mechanical modes can also be entangled via the routes b 1 a 1 b 2 and b 1 a 2 b 2 , where the phase difference ϕ plays a vital role. At ϕ = π , we can see that both the cavity–cavity and mechanical–mechanical entanglement firstly decrease from an initial value to zero when G 22 increases. With further increasing G 22 above a critical value, it is shown that the quantum entanglement starts to increase again. This phenomenon is closely related to the quantum interference between the eigenmodes of the system [34].
We proceed to study the effect of the entangling-interaction strength J 12 on the entanglement in Figure 3, assuming J 21 = J 12 for simplicity. Figure 3a and Figure 3c plot, respectively, the dynamics of the cavity–cavity entanglement and mechanical–mechanical entanglement for ϕ = 0 , π / 2 , and π . It is shown that stationary entanglement can be achieved after t 2 μs, which depends on the phase difference ϕ . Moreover, Figure 3b and Figure 3d plot the stationary cavity–cavity entanglement and mechanical–mechanical entanglement as functions of the entangling-interaction strength J 12 / 2 π for different values of phase difference ϕ . At ϕ = 0 , both the cavity–cavity entanglement and the mechanical–mechanical entanglement are enhanced monotonically with the increase in the coupling strength J 12 in the parameter regimes. When the phase difference ϕ is tuned to be π / 2 or π , the stationary logarithmic negativity E N becomes larger and reaches the maximum with the increase of J 12 . With further increasing J 12 , the entanglement will decrease and be zero above a critical value of J 12 . On one hand, it is straightforward that the entanglement can be enhanced by increasing the coupling strengths J 12 ( J 21 ) , which are strengthened by the blue-detuned driving fields. On the other hand, increasing J 12 ( J 21 ) can suppress the cavity-cooling effect for given G 12 ( G 21 ) and increase the effective temperature of the cavity and mechanical modes, which will reduce the entanglement. In addition, when the coupling strengths J 12 and J 21 are too strong, the system will become unstable. Therefore, the maximum entanglement is achieved by balancing these two competing effects [35]. Throughout this work, we have checked numerically that the chosen parameters satisfy the stability condition.
Furthermore, we study the phase dependence of the optomechanical entanglement in Figure 4. The logarithmic negativity E N between the two cavity modes is plotted as a function of the phase difference ϕ for different values of J 12 in Figure 4a. At J 12 / ( 2 π ) = 0.1 MHz, the maximum cavity–cavity entanglement is achieved around ϕ = π and the minimum value locates at ϕ = 0 ( 2 π ) . However, when the value of J 12 is increased to be 0.2 MHz and 0.25 MHz, the peak around ϕ = π turns into a dip. Especially, the cavity–cavity entanglement in the vicinity of ϕ = π is absent for J 12 = 0.25 MHz. Therefore, the phase difference ϕ can be exploited as a switch of the cavity–cavity entanglement between zero and nonzero value. Similarly, Figure 4b plots the logarithmic negativity E N between the two mechanical modes versus the phase difference ϕ . It is shown that the mechanical–mechanical entanglement reaches the minimum at ϕ = π and the maximum at ϕ = 0 ( 2 π ) for J 12 / 2 π = 0.1 , 0.2, and 0.25 MHz. Moreover, the entanglement can be strengthened monotonically by increasing the coupling strength J 12 at ϕ = 0 ( 2 π ) . Consequently, the quantum entanglement can be controlled by the phase difference in this optomechanical system.
As mentioned above, the entanglement between the two cavity (mechanical) modes is realized via the entanglement transfer between the cavity and mechanical modes. In Figure 5, we study the optomechanical entanglement between the cavity mode a 1 ( a 2 ) and mechanical mode b 2 ( b 1 ) . Figure 5a plots the time evolution of the logarithmic negativity E N between the cavity mode a 1 and mechanical mode b 2 with J 12 / 2 π = 0.15 MHz for different values of phase difference ϕ . It is shown that stationary entanglement can be achieved after t 2 μs and it is strengthened when the the phase difference ϕ is increased from 0 to π . The logarithmic negativity E N between a 1 and b 2 versus the entangling-interaction strength J 12 is plotted in Figure 5b, which shows that stationary entanglement can be enhanced by increasing the strength J 12 and reaches the maximum at some critical value for ϕ = 0 and π / 2 . In the chosen parameter regime, the stationary entanglement increases monotonically with the increase in J 12 for ϕ = π . In addition, Figure 5c plots the time evolution of the logarithmic negativity E N between a 2 and b 1 for ϕ = 0 , π / 2 , and π , where stationary entanglement can also be obtained. The stationary entanglement between a 2 and b 1 versus J 12 is plotted in Figure 5d, where the maximum entanglement is achieved around J 12 / 2 π = 0.15 MHz for various phase differences. Figure 5 also shows that the stationary entanglement between a 2 and b 1 is weaker than the entanglement between a 1 and b 2 , which is mainly due to G 22 > G 11 under consideration.
The robustness of the quantum entanglement with respect to the environmental temperature is shown in Figure 6. The value of E N decreases with the increasing temperature for various phase differences. Figure 6a plots the logarithmic negativity E N between the cavity modes a 1 and a 2 versus the temperature T. At ϕ = 0 , the cavity–cavity entanglement persists up to 150 mK. The entanglement vanishes at higher temperature for ϕ = π / 2 and ϕ = π . The entanglement between the two mechanical modes shown in Figure 6b is stronger than the cavity–cavity entanglement and survives at higher temperature. In addition, the entanglement between the cavity mode a 1 and mechanical mode b 2 versus the temperature is plotted in Figure 6c, where E N 0.11 can be obtained for T = 400 mK and ϕ = π . The logarithmic negativity E N between the cavity mode a 2 and mechanical mode b 1 is much smaller than the value in Figure 6c and survives at much lower temperature, which results from G 22 > G 11 . The value of E N between a 2 and b 1 can be increased by reducing G 22 , but the entanglement between the two cavity (mechanical) modes will be absent in some parameter regimes shown in Figure 2.
Finally, the entanglement usually survives at low final phonon numbers. The beam-splitter interaction such as G 11 ( a 1 b 1 + a 1 b 1 ) can not only result in state transfer between a 1 and b 1 but also cool the mechanical mode b 1 . In Figure 7, we numerically study the final average phonon numbers n 1 f and n 2 f by choosing the parameters for observing the entanglement mentioned above. Figure 7a plots the phonon numbers n 1 f and n 2 f as a function of the optomechanical coupling strength J 12 / 2 π with ϕ = 0 , which shows that both n 1 f and n 2 f increase monotonically with the enhancement of J 12 . The parametric-amplifier interaction J 12 ( a 1 b 2 + a 1 b 2 ) is caused by driving the cavity mode a 1 on the blue sideband of mechanical mode b 2 , where a drive photon is scattered into a photon at frequency ω 1 and a phonon at frequency Ω 2 , resulting in the monotonic increase in phonon number at larger J 12 . It is noted that quantum phase transition in the two-mode Rabi model has been recently investigated by analyzing the eigenenergy of the system [75], but quantum phase transition of the phonon number of the mechanical mode in optomechanical systems has seldom been reported. In addition, here, we show that both n 1 f and n 2 f are smaller than 1 in the chosen parameter regime. Therefore, simultaneous cooling of the mechanical modes and quantum entanglement can be realized in this optomechanical system. In addition, the phase dependence of phonon numbers is plotted in Figure 7b. The final average phonon numbers of the two mechanical modes reach the maximum at ϕ = π and the minimum at ϕ = 0 ( 2 π ) , which is contrary to the entanglement between the two mechanical modes shown in Figure 4b.

4. Conclusions

In conclusion, we have explored the phase-controlled entanglement in a multimode optomechanical system, where two cavity modes are, respectively, coupled to two mechanical modes via radiation pressure. Four tones are applied to the cavity modes to generate the entanglement, where each cavity is driven on the red sideband and blue sideband of the two mechanical modes, respectively. The dynamics of the bipartite entanglement is studied, which shows that the stationary entanglement can be generated. We have shown that the entanglement can be effectively modulated by the phase difference between the optomechanical coupling strengths that form a loop interaction. By tuning the optomechanical coupling strengths and phase difference in some parameter regimes, we can realize simultaneously the entanglement between the two cavity mode, between the two mechanical modes, and between the cavity and mechanical modes. We also find that the entanglement is robust against temperature and the mechanical resonators can be cooled to the ground state when the entanglement is generated.

Author Contributions

Conceptualization and model, C.J.; numerical simulation, H.L. and Z.Z.; writing—original draft preparation, C.J.; writing—review and editing, G.C. and C.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by National Natural Science Foundation of China (No: 11874170), Qinglan Project of Jiangsu Province of China, and Natural Science Research of Jiangsu Higher Education Institutions of China (19KJA150011).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of the Linearized Hamiltonian

We can linearize the Hamiltonian Equation (1) by the following two steps:
(1) Expanding the cavity modes as the sum of their steady-state field amplitudes and the small fluctuations [64,65,67,68], i.e., a i a i + j = 1 2 α i j e i ω d , i j t , where α i j is the coherent state amplitude produced in cavity i due to the drive. Then, we can obtain the following approximation:
a i a i a i j = 1 2 α i j e i ω d , i j t + a i j = 1 2 α i j e i ω d , i j t .
(2) Moving into the interaction picture with respect to [64,65,67,68],
H 0 / = i = 1 2 ω i a i a i + ( Ω i + δ i ) b i b i ,
the linearized Hamiltonian becomes
H / = i = 1 2 δ i b i b i + i = 1 2 j = 1 2 g 0 , i j ( a i j = 1 2 α i j e i ω d , i j e i ω i t + a i j = 1 2 α i j e i ω d , i j t e i ω i t ) [ b j e i ( Ω j + δ j ) t + b j e i ( Ω j + δ j ) t ] = i = 1 2 δ i b i b i + i = 1 2 j = 1 2 g 0 , i j [ a i ( α i 1 e i Δ i 1 t + α i 2 e i Δ i 2 ) + a i ( α i 1 e i Δ i 1 t + α i 2 e i Δ i 2 t ) ] [ b j e i ( Ω j + δ j ) t + b j e i ( Ω j + δ j ) t ] = i = 1 2 δ i b i b i + g 0 , 11 α 11 a 1 b 1 + g 0 , 11 α 11 a 1 b 1 + g 0 , 12 α 12 a 1 b 2 + g 0 , 12 α 12 a 1 b 2 + g 0 , 21 α 21 a 2 b 1 + g 0 , 21 α 21 a 2 b 1 + g 0 , 22 α 22 a 2 b 2 + g 0 , 22 α 22 a 2 b 2 + H RWA + H off ,
where H RWA describes the rapidly oscillating terms at frequencies ± 2 [ Ω 1 ( 2 ) + δ 1 ( 2 ) ] and ± ( Ω 1 + δ 1 + Ω 2 + δ 2 ) , and H off denotes the off-resonant terms oscillating at ± ( Ω 1 + δ 1 Ω 2 δ 2 ) . Here, Δ i j = ω i ω d , i j is the detuning of the drive tone with respect to cavity a i , and we consider the situation where
Δ 11 = Ω 1 + δ 1 , Δ 12 = ( Ω 2 + δ 2 ) , Δ 21 = ( Ω 1 + δ 1 ) , Δ 22 = Ω 2 + δ 2 .
Under the condition Ω i , | Ω 1 Ω 2 | > { κ i , G 11 , G 22 , J 12 , J 21 } , the rotation wave approximation (RWA) can be applied and the off-resonant terms can be ignored. Defining g 0 , 11 α 11 = G 11 e i ϕ 11 , g 0 , 22 α 22 = G 22 e i ϕ 22 , g 0 , 12 α 12 = J 12 e i ϕ 12 , and g 0 , 21 α 21 = J 21 e i ϕ 21 as the field-enhanced optomechanical coupling strengths, the linearized Hamiltonian can be written as
H / = i = 1 2 δ i b i b i + ( G 11 e i ϕ 11 a 1 b 1 + J 12 e i ϕ 12 a 1 b 2 + J 21 e i ϕ 21 a 2 b 1 + G 22 e i ϕ 22 a 2 b 2 + H . c ) ,
which is Equation (2) in the main text.

References

  1. Berkley, A.J.; Xu, H.; Ramos, R.C.; Gubrud, M.A.; Strauch, F.W.; Johnson, P.R.; Anderson, J.R.; Dragt, A.J.; Lobb, C.J.; Wellstood, F.C. Entangled macroscopic quantum states in two superconducting qubits. Science 2003, 300, 1548–1550. [Google Scholar] [CrossRef] [PubMed]
  2. Julsgaard, B.; Kozhekin, A.; Polzik, E.S. Experimental long-lived entanglement of two macroscopic objects. Nature 2001, 413, 400–403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Chou, C.W.; Riedmatten, H.; Felinto, D.; Polyakov, S.V.; van Enk, S.J.; Kimble, H.J. Measurement-induced entanglement for excitation stored in remote atomic ensembles. Nature 2005, 438, 828–832. [Google Scholar] [CrossRef] [Green Version]
  4. Ritter, S.; Nölleke, C.; Hahn, C.; Reiserer, A.; Neuzner, A.; Uphoff, M.; Mücke, M.; Figueroa, E.; Bochmann, J.; Rempe, G. An elementary quantum network of single atoms in optical cavities. Nature 2012, 484, 195–200. [Google Scholar] [CrossRef] [Green Version]
  5. Moehring, D.L.; Maunz, P.; Olmschenk, S.; Younge, K.C.; Matsukevich, D.N.; Duan, L.-M.; Monroe, C. Entanglement of single-atom quantum bits at a distance. Nature 2007, 449, 68–71. [Google Scholar] [CrossRef] [Green Version]
  6. Hensen, B.; Bernien, H.; Dréau, A.E.; Reiserer, A.; Kalb, N.; Blok, M.S.; Ruitenberg, J.; Vermeulen, R.F.L.; Schouten, R.N.; Abellán, C.; et al. Loophole-free Bell inequality violation using electron spins separated by 1.3 kilometres. Nature 2015, 526, 682–686. [Google Scholar] [CrossRef]
  7. O’Connell, A.D.; Hofheinz, M.; Ansmann, M.; Bialczak, R.C.; Lenander, M.; Lucero, E.; Neeley, M.; Sank, D.; Wang, H.; Weides, M.; et al. Quantum ground state and single-phonon control of a mechanical resonator. Nature 2010, 464, 697–703. [Google Scholar] [CrossRef]
  8. Yeo, I.; Assis, P.L.; Gloppe, A.; Dupont-Ferrier, E.; Verlot, P.; Malik, N.S.; Dupuy, E.; Claudon, J.; Gérard, J.-M.; Auffêves, A.; et al. Strain-mediated coupling in a quantum dot-mechanical oscillator system. Nat. Nanotechnol. 2014, 9, 106–110. [Google Scholar] [CrossRef] [Green Version]
  9. Kolkowitz, S.; Jayich, A.C.B.; Unterreithmeier, Q.P.; Bennett, S.D.; Rabl, P.; Harris, J.G.E.; Lukin, M.D. Coherent sensing of a mechanical resonator with a single-spin qubit. Science 2012, 335, 1603–1606. [Google Scholar] [CrossRef] [Green Version]
  10. Hunger, D.; Camerer, S.; Hänsch, T.W.; König, D.; Kotthaus, J.P.; Reichel, J.; Treutlein, P. Resonant coupling of a Bose-Einstein Condensate to a micromechanical oscillator. Phys. Rev. Lett. 2010, 104, 143002. [Google Scholar] [CrossRef]
  11. Thomas, R.A.; Parniak, M.; Østfeldt, C.; Møller, C.B.; Bærentsen, C.; Tsaturyan, Y.; Schliesser, A.; Appel, J.; Zeuthen, E.; Polzik, E.S. Entanglement between distant macroscopic mechanical and spin systems. Nat. Phys. 2021, 17, 228–233. [Google Scholar] [CrossRef]
  12. Bienfait, A.; Satzinger, K.J.; Zhong, Y.P.; Chang, H.S.; Chou, M.H.; Conner, C.R.; Dumur, E.; Grebel, J.; Peairs, G.A.; Povey, R.G.; et al. Phonon-mediated quantum state transfer and remote qubit entanglement. Science 2019, 364, 368–371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Kimble, H.J. The quantum internet. Nature 2008, 453, 1023–1030. [Google Scholar] [CrossRef] [PubMed]
  14. Aspelmeyer, M.; Kippenberg, T.J.; Marquardt, F. Cavity optomechanics. Rev. Mod. Phys. 2014, 86, 1391–1452. [Google Scholar] [CrossRef] [Green Version]
  15. Stannigel, K.; Rabl, P.; Sørensen, A.S.; Zoller, P.; Lukin, M.D. Optomechanical transducers for long-distance quantum communication. Phys. Rev. Lett. 2010, 105, 220501. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Y.D.; Clerk, A.A. Using interference for high fidelity quantum state transfer in optomechanics. Phys. Rev. Lett. 2012, 108, 153603. [Google Scholar] [CrossRef] [Green Version]
  17. Tian, L. Adiabatic state conversion and pulse transmission in optomechanical systems. Phys. Rev. Lett. 2012, 108, 153604. [Google Scholar] [CrossRef] [Green Version]
  18. Barzanjeh, S.; Abdi, M.; Milburn, G.J.; Tombesi, P.; Vitali, D. Reversible optical-to-microwave quantum interface. Phys. Rev. Lett. 2012, 109, 130503. [Google Scholar] [CrossRef]
  19. Andrews, R.W.; Peterson, R.W.; Purdy, T.P.; Cicak, K.; Simmonds, R.W.; Regal, C.A.; Lehnert, K.W. Bidirectional and efficient conversion between microwave and optical light. Nat. Phys. 2014, 10, 321–326. [Google Scholar] [CrossRef] [Green Version]
  20. Dong, C.H.; Wang, Y.D.; Wang, H.L. Optomechanical interfaces for hybrid quantum networks. Natl. Sci. Rev. 2015, 2, 510–519. [Google Scholar] [CrossRef]
  21. Forsch, M.; Stockill, R.; Wallucks, A.; Marinković, I.; Gärtner, C.; Norte, R.A.; van Otten, F.; Fiore, A.; Srinivasan, K.; Gröblacher, S. Microwave-to-Optics conversion using a mechanical oscillator in its quantum ground state. Nat. Phys. 2020, 16, 69–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Jiang, W.T.; Sarabalis, C.J.; Dahmani, Y.D.; Patel, R.N.; Mayor, F.M.; McKenna, T.P.; Laer, R.V.; Safavi-Naeini, A.H. Efficient bidirectional piezo-optomechanical transduction between microwave and optical frequency. Nat. Commun. 2020, 11, 1166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Teufel, J.D.; Donner, T.; Li, D.; Harlow, J.W.; Allman, M.S.; Cicak, K.; Sirois, A.J.; Whittaker, J.D.; Lehnert, K.W.; Simmonds, R.W. Sideband cooling of micromechanical motion to the quantum ground state. Nature 2011, 475, 359–363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chan, J.; Alegre, T.P.M.; Safavi-Naeini, A.H.; Hill, J.T.; Krause, A.; Gröblacher, S.; Aspelmeyer, M.; Painter, O. Laser cooling of a nanomechanical oscillator into its quantum ground state. Nature 2011, 478, 89–92. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Delić, U.; Reisenbauer, M.; Dare, K.; Grass, D.; Vuletic, V.; Kiesel, N.; Aspelmeyer, M. Cooling of a levitated nanoparticle to the motional quantum ground state. Science 2020, 367, 892–895. [Google Scholar] [CrossRef] [Green Version]
  26. Safavi-Naeini, A.H.; Gröblacher, S.; Hill, J.T.; Chan, J.; Aspelmeyer, M.; Painter, O. Squeezed light from a silicon micromechanical resonator. Nature 2013, 500, 185–189. [Google Scholar] [CrossRef] [Green Version]
  27. Purdy, T.P.; Yu, P.L.; Peterson, R.W.; Kampel, N.S.; Regal, C.A. Strong optomechanical squeezing of light. Phys. Rev. X 2013, 3, 031012. [Google Scholar] [CrossRef] [Green Version]
  28. Wollman, E.E.; Lei, C.U.; Weinstein, A.J.; Suh, J.; Kronwald, A.; Marquardt, F.; Clerk, A.A.; Schwab, K.C. Quantum squeezing of motion in a mechanical resonator. Science 2015, 349, 952–955. [Google Scholar] [CrossRef] [Green Version]
  29. Pirkkalainen, J.M.; Damskägg, E.; Brandt, M.; Massel, F.; Sillanpää, M.A. Squeezing of quantum noise of motion in a micromechanical resonator. Phys. Rev. Lett. 2015, 115, 243601. [Google Scholar] [CrossRef] [Green Version]
  30. Vitali, D.; Gigan, S.; Ferreira, A.; Böhm, H.R.; Tombesi, P.; Guerreiro, A.; Vedral, V.; Zeilinger, A.; Aspelmeyer, M. Optomechanical entanglement between a movable mirror and a cavity field. Phys. Rev. Lett. 2007, 98, 030405. [Google Scholar] [CrossRef]
  31. Genes, C.; Mari, A.; Tombesi, P.; Vitali, D. Robust entanglement of a micromechanical resonator with output optical fields. Phys. Rev. A 2008, 78, 032316. [Google Scholar] [CrossRef] [Green Version]
  32. Barzanjeh, S.; Vitali, D.; Tombesi, P.; Milburn, G.J. Entangling optical and microwave cavity modes by means of a nanomechanical resonator. Phys. Rev. A 2011, 84, 042342. [Google Scholar] [CrossRef] [Green Version]
  33. Kuzyk, M.C.; Enk, S.J.; Wang, H.L. Generating robust optical entanglement in weak-coupling optomechanical systems. Phys. Rev. A 2013, 88, 062341. [Google Scholar] [CrossRef] [Green Version]
  34. Tian, L. Robust photon entanglement via quantum interference in optomechanical interfaces. Phys. Rev. Lett. 2013, 110, 233602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wang, Y.D.; Clerk, A.A. Reservoir-engineered entanglement in optomechanical systems. Phys. Rev. Lett. 2013, 110, 253601. [Google Scholar] [CrossRef] [Green Version]
  36. Wang, G.L.; Huang, L.; Lai, Y.C.; Grebogi, C. Nonlinear dynamics and quantum entanglement in optomechanical systems. Phys. Rev. Lett. 2014, 112, 110406. [Google Scholar] [CrossRef] [Green Version]
  37. Liao, J.Q.; Wu, Q.Q.; Nori, F. Entangling two macroscopic mechanical mirrors in a two-cavity optomechanical system. Phys. Rev. A 2014, 89, 014302. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, M.; Lu, X.Y.; Wang, Y.D.; You, J.Q.; Wu, Y. Macroscopic quantum entanglement in modulated optomechanics. Phys. Rev. A 2016, 94, 053807. [Google Scholar] [CrossRef] [Green Version]
  39. Sun, F.X.; Mao, D.; Dai, Y.T.; Ficek, Z.; He, Q.Y.; Gong, Q.H. Phase control of entanglement and quantum steering in a three-mode optomechanical system. New J. Phys. 2017, 19, 123039. [Google Scholar] [CrossRef]
  40. Li, J.; Li, G.; Zippilli, S.; Vitali, D.; Zhang, T. Enhanced entanglement of two different mechanical resonators via coherent feedback. Phys. Rev. A 2017, 95, 043819. [Google Scholar] [CrossRef]
  41. Qin, L.G.; Wang, Z.Y.; Ma, H.Y.; Zhang, C.M.; Ren, L.; Wang, L.L.; Gong, S.Q. Optomechanical entanglement switch in the hybrid opto-electromechanical device. JOSA B 2019, 36, 1544–1550. [Google Scholar] [CrossRef]
  42. Hu, C.S.; Liu, Z.Q.; Liu, Y.; Shen, L.T.; Wu, H.; Zheng, S.B. Entanglement beating in a cavity optomechanical system under two-field driving. Phys. Rev. A 2020, 101, 033810. [Google Scholar] [CrossRef]
  43. Zhang, Y.L.; Yang, C.S.; Shen, Z.; Dong, C.H.; Guo, G.C.; Zou, C.L.; Zou, X.B. Enhanced optomechanical entanglement and cooling via dissipation engineering. Phys. Rev. A 2020, 101, 063836. [Google Scholar] [CrossRef]
  44. Zhong, C.C.; Han, X.; Tang, H.X.; Jiang, L. Entanglement of microwave-optical modes in a strongly coupled electro-optomechanical system. Phys. Rev. A 2020, 101, 032345. [Google Scholar] [CrossRef] [Green Version]
  45. Li, G.Y.; Nie, W.J.; Wu, Y.; Liao, Q.; Chen, A.; Lan, Y. Manipulating the steady-state entanglement via three-level atoms in a hybrid levitated optomechanical system. Phys. Rev. A 2020, 102, 063501. [Google Scholar] [CrossRef]
  46. Jiao, Y.F.; Zhang, S.D.; Zhang, Y.L.; Miranowicz, A.; Kuang, L.M.; Jing, H. Nonreciprocal optomechanical entanglement against backscattering losses. Phys. Rev. Lett. 2020, 125, 143605. [Google Scholar] [CrossRef]
  47. Lin, Q.; He, B.; Xiao, M. Entangling two macroscopic mechanical resonators at high temperature. Phys. Rev. Appl. 2020, 13, 034030. [Google Scholar] [CrossRef] [Green Version]
  48. Jiang, C.; Tserkis, S.; Collins, K.; Onoe, S.; Li, Y.; Tian, L. Switchable bipartite and genuine tripartite entanglement via an optoelectromechanical interface. Phys. Rev. A 2020, 101, 042320. [Google Scholar] [CrossRef] [Green Version]
  49. Palomaki, T.A.; Teufel, J.D.; Simmonds, R.W.; Lehnert, K.W. Entangling mechanical motion with microwave fields. Science 2013, 342, 710–713. [Google Scholar] [CrossRef]
  50. Ockeloen-Korppi, C.F.; Damskägg, E.; Pirkkalainen, J.M.; Asjad, M.; Clerk, A.A.; Massel, F.; Woolley, M.J.; Sillanpää, M.A. Stabilized entanglement of massive mechanical oscillators. Nature 2018, 556, 478–482. [Google Scholar] [CrossRef]
  51. Riedinger, R.; Wallucks, A.; Marinković, I.; Löschnauzer, C.; Aspelmeyer, M.; Hong, S.; Gröblacher, S. Remote quantum entanglement between two micromechanical oscillators. Nature 2018, 556, 473–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. de Lépinay, L.M.; Ockeloen-Korppi, C.F.; Woolley, M.J.; Sillanpää, M.A. Quantum mechanics–free subsystem with mechanical oscillators. Science 2021, 372, 625. [Google Scholar] [CrossRef] [PubMed]
  53. Kotler, S.; Peterson, G.A.; Shojaee, E.; Lecocq, F.; Cicak, K.; Kwiatkowski, A.; Geller, S.; Glancy, S.; Knill, E.; Simmonds, R.W.; et al. Direct observation of deterministic macroscopic entanglement. Science 2021, 372, 622. [Google Scholar] [CrossRef]
  54. Chen, J.X.; Rossi, M.; Mason, D.; Schliesser, A. Entanglement of propagating optical modes via a mechanical interface. Nat. Commun. 2020, 11, 943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Barzanjeh, S.; Redchenko, E.S.; Peruzzo, M.; Wulf, M.; Lewis, D.P.; Arnold, G.; Fink, J.M. Stationary entangled radiation from micromechanical motion. Nature 2019, 570, 480–483. [Google Scholar] [CrossRef] [Green Version]
  56. Xu, X.W.; Li, Y.; Chen, A.X.; Liu, Y. Nonreciprocal conversion between microwave and optical photons in electro-optomechanical systems. Phys. Rev. A 2016, 93, 023827. [Google Scholar] [CrossRef] [Green Version]
  57. Li, Y.; Huang, Y.Y.; Zhang, X.Z.; Tian, L. Optical directional amplification in a three-mode optomechanical system. Opt. Express 2017, 25, 18907. [Google Scholar] [CrossRef] [Green Version]
  58. Zhang, X.Z.; Tian, L.; Li, Y. Optomechanical transistor with mechanical gain. Phys. Rev. A 2018, 97, 043818. [Google Scholar] [CrossRef] [Green Version]
  59. Tian, L.; Li, Z. Nonreciprocal quantum-state conversion between microwave and optical photons. Phys. Rev. A 2017, 96, 013808. [Google Scholar] [CrossRef] [Green Version]
  60. Jiang, C.; Song, L.N.; Li, Y. Directional amplifier in an optomechanical system with optical gain. Phys. Rev. A 2018, 97, 053812. [Google Scholar] [CrossRef]
  61. Ruesink, F.; Miri, M.A.; Alù, A.; Verhagen, E. Nonreciprocity and magnetic-free isolation based on optomechanical interactions. Nat. Commun. 2016, 7, 13662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Miri, M.A.; Ruesink, F.; Verhagen, E.; Alù, A. Optical nonreciprocity based on optomechanical coupling. Phys. Rev. Appl. 2017, 7, 064014. [Google Scholar] [CrossRef] [Green Version]
  63. Fang, K.J.; Luo, J.; Metelmann, A.; Matheny, M.H.; Marquardt, F.; Clerk, A.A.; Painter, O. Generalized non-reciprocity in an optomechanical circuit via synthetic magnetism and reservoir engineering. Nat. Phys. 2017, 13, 465. [Google Scholar] [CrossRef] [Green Version]
  64. Bernier, N.R.; Tóth, L.D.; Koottavida, A.; Ioannou, M.A.; Malz, D.; Nunnenkamp, A.; Feofanov, A.K.; Kippenberg, T.J. Nonreciprocal reconfigurable microwave optomechanical circuit. Nat. Commun. 2017, 8, 604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Barzanjeh, S.; Wulf, M.; Peruzzo, M.; Kalaee, M.; Dieterle, P.B.; Painter, O.; Fink, J.M. Mechanical on-chip microwave circulator. Nat. Commun. 2017, 8, 953. [Google Scholar] [CrossRef] [Green Version]
  66. Peterson, G.A.; Lecocq, F.; Cicak, K.; Simmonds, R.W.; Aumentado, J.; Teufel, J.D. Demonstration of efficient nonreciprocity in a microwave optomechanical circuit. Phys. Rev. X 2017, 7, 031001. [Google Scholar] [CrossRef] [Green Version]
  67. de Lépinay, L.M.; Damskägg, E.; Ockeloen-Korppi, C.F.; Sillanpää, M.A. Realization of directional amplification in a microwave optomechanical device. Phys. Rev. Appl. 2019, 11, 034027. [Google Scholar] [CrossRef] [Green Version]
  68. Malz, D.; Tóth, L.D.; Bernier, N.R.; Feofanov, A.K.; Kippenberg, T.J.; Nunnenkamp, A. Quantum-limited directional amplifiers with optomechanics. Phys. Rev. Lett. 2018, 120, 023601. [Google Scholar] [CrossRef] [Green Version]
  69. Shen, Z.; Zhang, Y.L.; Chen, Y.; Sun, F.W.; Zou, X.B.; Guo, G.C.; Zou, C.L.; Dong, C.H. Reconfigurable optomechanical circulator and directional amplifier. Nat. Commun. 2018, 9, 1797. [Google Scholar] [CrossRef] [Green Version]
  70. Ruesink, F.; Mathew, J.P.; Miri, M.A.; Alù, A.; Verhagen, E. Optical circulation in a multimode optomechanical resonator. Nat. Commun. 2018, 9, 1798. [Google Scholar] [CrossRef]
  71. Lai, D.G.; Liao, J.Q.; Miranowicz, A.; Nori, F. Noise-tolerant optomechanical entanglement via synthetic magnetism. Phys. Rev. Lett. 2022, 129, 063602. [Google Scholar] [CrossRef] [PubMed]
  72. Liu, J.X.; Jiao, Y.F.; Li, Y.; Xu, X.W.; He, Q.Y.; Jing, H. Phase-controlled asymmetric optomechanical entanglement against defect losses. arXiv 2022, arXiv:2209.12508. [Google Scholar]
  73. DeJesus, E.X.; Kaufman, C. Routh-Hurwitz criterion in the examination of eigenvalues of a system of nonlinear ordinary differential equations. Phys. Rev. A 1987, 35, 5288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Adesso, G.; Serafini, A.; Illuminati, F. Extremal entanglement and mixedness in continuous variable systems. Phys. Rev. A 2004, 70, 022318. [Google Scholar] [CrossRef] [Green Version]
  75. Shen, L.T.; Yang, J.W.; Zhong, Z.R.; Yang, Z.B.; Zheng, S.B. Quantum phase transition and quench dynamics in the two-mode Rabi model. Phys. Rev. A 2021, 104, 063703. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic of the optomechanical system. Two cavity modes a 1 and a 2 are, respectively, coupled to two mechanical modes b 1 and b 2 with the effective optomechanical coupling strengths G 11 , J 12 , J 21 , and G 22 . (b) Frequency domain schematic of the four-tone driving scheme. Cavity mode a 1 with resonance frequency ω 1 is driven by two tones at frequencies ω d , 11 = ω 1 ( Ω 1 + δ 1 ) and ω d , 12 = ω 1 + ( Ω 2 + δ 2 ) . Cavity mode a 2 with resonance frequency ω 2 is driven by two tones at frequencies ω d , 21 = ω 2 + ( Ω 1 + δ 1 ) and ω d , 22 = ω 2 ( Ω 2 + δ 2 ) .
Figure 1. (a) Schematic of the optomechanical system. Two cavity modes a 1 and a 2 are, respectively, coupled to two mechanical modes b 1 and b 2 with the effective optomechanical coupling strengths G 11 , J 12 , J 21 , and G 22 . (b) Frequency domain schematic of the four-tone driving scheme. Cavity mode a 1 with resonance frequency ω 1 is driven by two tones at frequencies ω d , 11 = ω 1 ( Ω 1 + δ 1 ) and ω d , 12 = ω 1 + ( Ω 2 + δ 2 ) . Cavity mode a 2 with resonance frequency ω 2 is driven by two tones at frequencies ω d , 21 = ω 2 + ( Ω 1 + δ 1 ) and ω d , 22 = ω 2 ( Ω 2 + δ 2 ) .
Photonics 09 00818 g001
Figure 2. Contour plots of the logarithmic negativity E N between (a) the two cavity modes and (b) the two mechanical modes as functions of the coupling strength G 22 / 2 π and the phase difference ϕ / π . Parameters are Ω 1 / 2 π = 4.34 MHz, Ω 2 / 2 π = 5.64 MHz, κ 1 / 2 π = 2.42 MHz, κ 2 / 2 π = 1.98 MHz, γ 1 / 2 π = 40 Hz, γ 2 / 2 π = 80 Hz, δ 1 = δ 2 = 0 , G 11 / 2 π = 0.5 MHz, J 21 = J 12 = 2 π × 0.15 MHz, and T = 20 mK.
Figure 2. Contour plots of the logarithmic negativity E N between (a) the two cavity modes and (b) the two mechanical modes as functions of the coupling strength G 22 / 2 π and the phase difference ϕ / π . Parameters are Ω 1 / 2 π = 4.34 MHz, Ω 2 / 2 π = 5.64 MHz, κ 1 / 2 π = 2.42 MHz, κ 2 / 2 π = 1.98 MHz, γ 1 / 2 π = 40 Hz, γ 2 / 2 π = 80 Hz, δ 1 = δ 2 = 0 , G 11 / 2 π = 0.5 MHz, J 21 = J 12 = 2 π × 0.15 MHz, and T = 20 mK.
Photonics 09 00818 g002
Figure 3. Logarithmic negativity E N affected by the phase difference ϕ . (a) Time evolution of the logarithmic negativity E N between the cavity modes a 1 and a 2 with J 12 / 2 π = 0.2 MHz for ϕ = 0 , π / 2 , and π . (b) Stationary logarithmic negativity E N between the cavity modes a 1 and a 2 as a function of J 12 / 2 π for ϕ = 0 , π / 2 , and π . (c) Time evolution of the logarithmic negativity E N between the mechanical modes b 1 and b 2 with J 12 / 2 π = 0.2 MHz for ϕ = 0 , ϕ = π / 2 , and ϕ = π . (d) Stationary logarithmic negativity E N between the mechanical modes b 1 and b 2 as a function of J 12 / 2 π for ϕ = 0 , π / 2 , and π . Here, ϕ = 0 , π / 2 , π correspond to the solid line, dashed line, and dash-dotted line, respectively. Other parameters are the same as Figure 2, except G 22 / 2 π = 1 MHz. Here, we keep J 21 = J 12 .
Figure 3. Logarithmic negativity E N affected by the phase difference ϕ . (a) Time evolution of the logarithmic negativity E N between the cavity modes a 1 and a 2 with J 12 / 2 π = 0.2 MHz for ϕ = 0 , π / 2 , and π . (b) Stationary logarithmic negativity E N between the cavity modes a 1 and a 2 as a function of J 12 / 2 π for ϕ = 0 , π / 2 , and π . (c) Time evolution of the logarithmic negativity E N between the mechanical modes b 1 and b 2 with J 12 / 2 π = 0.2 MHz for ϕ = 0 , ϕ = π / 2 , and ϕ = π . (d) Stationary logarithmic negativity E N between the mechanical modes b 1 and b 2 as a function of J 12 / 2 π for ϕ = 0 , π / 2 , and π . Here, ϕ = 0 , π / 2 , π correspond to the solid line, dashed line, and dash-dotted line, respectively. Other parameters are the same as Figure 2, except G 22 / 2 π = 1 MHz. Here, we keep J 21 = J 12 .
Photonics 09 00818 g003
Figure 4. Plots of the logarithmic negativity E N between (a) the two cavity modes and (b) the two mechanical modes versus the phase difference ϕ / π for J 12 / 2 π = 0.1 MHz (solid lines), 0.2 MHz (dashed lines), and 0.25 MHz (dash-dotted lines). The other parameters are the same as those in Figure 3.
Figure 4. Plots of the logarithmic negativity E N between (a) the two cavity modes and (b) the two mechanical modes versus the phase difference ϕ / π for J 12 / 2 π = 0.1 MHz (solid lines), 0.2 MHz (dashed lines), and 0.25 MHz (dash-dotted lines). The other parameters are the same as those in Figure 3.
Photonics 09 00818 g004
Figure 5. Logarithmic negativity E N between the cavity and mechanical modes affected by the phase difference ϕ . (a) Time evolution of the logarithmic negativity E N between the cavity mode a 1 and mechanical mode b 2 with J 12 / 2 π = 0.15 MHz for ϕ = 0 , π / 2 , and π . (b) Stationary logarithmic negativity E N between the cavity mode a 1 and mechanical mode b 2 as a function of J 12 / 2 π for various phase differences. (c) Time evolution of the logarithmic negativity E N between the cavity mode a 2 and mechanical mode b 1 with J 12 / 2 π = 0.15 MHz for ϕ = 0 , π / 2 , and π . (d) Stationary logarithmic negativity E N between the cavity mode a 2 and mechanical mode b 1 versus the J 12 / 2 π for various phase differences. Here, the solid line, dashed line, and dash-dotted line represent ϕ = 0 , π / 2 , and π , respectively. The other parameters are the same as those in Figure 3.
Figure 5. Logarithmic negativity E N between the cavity and mechanical modes affected by the phase difference ϕ . (a) Time evolution of the logarithmic negativity E N between the cavity mode a 1 and mechanical mode b 2 with J 12 / 2 π = 0.15 MHz for ϕ = 0 , π / 2 , and π . (b) Stationary logarithmic negativity E N between the cavity mode a 1 and mechanical mode b 2 as a function of J 12 / 2 π for various phase differences. (c) Time evolution of the logarithmic negativity E N between the cavity mode a 2 and mechanical mode b 1 with J 12 / 2 π = 0.15 MHz for ϕ = 0 , π / 2 , and π . (d) Stationary logarithmic negativity E N between the cavity mode a 2 and mechanical mode b 1 versus the J 12 / 2 π for various phase differences. Here, the solid line, dashed line, and dash-dotted line represent ϕ = 0 , π / 2 , and π , respectively. The other parameters are the same as those in Figure 3.
Photonics 09 00818 g005
Figure 6. Plots of logarithmic negativity E N as functions of the environment temperature T for ϕ = 0 (solid lines), π / 2 (dashed lines), and π (dash-dotted lines). Panel (a) represents the entanglement between the two cavity modes, panel (b) represents the entanglement between the two mechanical modes, panel (c) corresponds to the entanglement between the cavity mode a 1 and mechanical mode b 2 , and panel (d) corresponds to the entanglement between the cavity mode a 2 and mechanical mode b 1 . The other parameters are the same as those in Figure 3, except J 12 / 2 π = J 21 / 2 π = 0.15 MHz.
Figure 6. Plots of logarithmic negativity E N as functions of the environment temperature T for ϕ = 0 (solid lines), π / 2 (dashed lines), and π (dash-dotted lines). Panel (a) represents the entanglement between the two cavity modes, panel (b) represents the entanglement between the two mechanical modes, panel (c) corresponds to the entanglement between the cavity mode a 1 and mechanical mode b 2 , and panel (d) corresponds to the entanglement between the cavity mode a 2 and mechanical mode b 1 . The other parameters are the same as those in Figure 3, except J 12 / 2 π = J 21 / 2 π = 0.15 MHz.
Photonics 09 00818 g006
Figure 7. The final average phonon numbers n 1 f (solid line) and n 2 f (dashed line) versus (a) the optomechanical coupling strength J 12 / 2 π with ϕ = 0 and (b) the phase difference ϕ / π with J 12 / 2 π = J 21 / 2 π = 0.2 MHz. The other parameters are the same as those in Figure 3.
Figure 7. The final average phonon numbers n 1 f (solid line) and n 2 f (dashed line) versus (a) the optomechanical coupling strength J 12 / 2 π with ϕ = 0 and (b) the phase difference ϕ / π with J 12 / 2 π = J 21 / 2 π = 0.2 MHz. The other parameters are the same as those in Figure 3.
Photonics 09 00818 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, C.; Lu, H.; Zhai, Z.; Chen, G. Phase-Controlled Entanglement in a Four-Mode Optomechanical System. Photonics 2022, 9, 818. https://doi.org/10.3390/photonics9110818

AMA Style

Jiang C, Lu H, Zhai Z, Chen G. Phase-Controlled Entanglement in a Four-Mode Optomechanical System. Photonics. 2022; 9(11):818. https://doi.org/10.3390/photonics9110818

Chicago/Turabian Style

Jiang, Cheng, Hongxia Lu, Zhangyin Zhai, and Guibin Chen. 2022. "Phase-Controlled Entanglement in a Four-Mode Optomechanical System" Photonics 9, no. 11: 818. https://doi.org/10.3390/photonics9110818

APA Style

Jiang, C., Lu, H., Zhai, Z., & Chen, G. (2022). Phase-Controlled Entanglement in a Four-Mode Optomechanical System. Photonics, 9(11), 818. https://doi.org/10.3390/photonics9110818

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop