Quantum correlation in a nano-electro-optomechanical system enhanced by an optical parametric amplifier and Coulomb-type interaction

In this paper, we investigated the quantum correlation of nano-electro-optomechanical system enhanced by an optical parametric amplifier (OPA) and Coulomb-type interaction. In particular, we consider a hybrid system consisting of a cavity and two charged mechanical oscillators with an OPA, where the optical cavity mode is coupled with a charged mechanical oscillator via radiation pressure, and the two charged mechanical oscillators are coupled through a Coulomb interaction. We use logarithmic negativity to quantify quantum entanglement, and quantum discord to measure the quantumness correlation between the two mechanical oscillators. We characterize quantum steering using the steerability between the two mechanical oscillators. Our results show that the presence of OPA and strong Coulomb coupling enhances the quantum correlations between the two mechanical oscillators. In addition, Coulomb interactions are more prominent in quantum correlations. Besides, in the presence of OPA, the maximum amount of quantum entanglement, quantum steering, and quantum discord were achieved between the two mechanical oscillators is greater than in the absence of OPA. Moreover, a proper phase choice of the optical field driving the OPA enhances quantum correlations under suitable conditions. We obtain quantum entanglement confines quantum steering and quantum discord beyond entanglement. Furthermore, quantum entanglement, quantum steering, and quantum discord decrease rapidly with increasing temperature as a result of decoherence. In addition, quantum discord persists at higher temperature values, although the quantum entanglement between the systems also vanishes completely. Our proposed scheme enhances quantum correlation and proves robust against fluctuations in the bath environment. We believe that the present scheme of quantum correlation provides a promising platform for the realization of continuous variable quantum information processing.


Quantum correlation in a nano-electro-optomechanical system enhanced by an optical parametric amplifier and Coulomb-type interaction
Habtamu Dagnaw Mekonnen 1,2 , Tesfay Gebremariam Tesfahannes 3* , Tewodros Yirgashewa Darge 1,4 & Alemayehu Getahun Kumela 1,4 In this paper, we investigated the quantum correlation of nano-electro-optomechanical system enhanced by an optical parametric amplifier (OPA) and Coulomb-type interaction.In particular, we consider a hybrid system consisting of a cavity and two charged mechanical oscillators with an OPA, where the optical cavity mode is coupled with a charged mechanical oscillator via radiation pressure, and the two charged mechanical oscillators are coupled through a Coulomb interaction.We use logarithmic negativity to quantify quantum entanglement, and quantum discord to measure the quantumness correlation between the two mechanical oscillators.We characterize quantum steering using the steerability between the two mechanical oscillators.Our results show that the presence of OPA and strong Coulomb coupling enhances the quantum correlations between the two mechanical oscillators.In addition, Coulomb interactions are more prominent in quantum correlations.Besides, in the presence of OPA, the maximum amount of quantum entanglement, quantum steering, and quantum discord were achieved between the two mechanical oscillators is greater than in the absence of OPA.Moreover, a proper phase choice of the optical field driving the OPA enhances quantum correlations under suitable conditions.We obtain quantum entanglement confines quantum steering and quantum discord beyond entanglement.Furthermore, quantum entanglement, quantum steering, and quantum discord decrease rapidly with increasing temperature as a result of decoherence.In addition, quantum discord persists at higher temperature values, although the quantum entanglement between the systems also vanishes completely.Our proposed scheme enhances quantum correlation and proves robust against fluctuations in the bath environment.We believe that the present scheme of quantum correlation provides a promising platform for the realization of continuous variable quantum information processing.
Quantum correlation and entanglement are special types of quantum coherence, and present many charming properties to realize quantum information processing 1-3 .Consequentlly, quantum correlations are widely used in numerous tasks of quantum information processing, for instance, in quantum communication 4 , computation 5 and metrology 6 .Specifically, quantum entanglement is a key ingredient of quantum information processing that characterize the non-classical property of multipartite quantum systems 7 .Accordingly, numerous researchers show that the nonzero entanglement assures the existence of quantum correlations but the zero entanglement does not assure the absence of quantum correlations in a bipartite quantum state.Unlike the entanglement, the quantum steering can be certifying and judge for creation and verification of optomechanical entanglement 8 .Subsequently, both theoretical and experimental researchers have precisely indicated that the observation of quantum steering is an essential resource in several fundamental applications [9][10][11][12][13] .Thus, the quantum steering is used to quantify how much the two entangled bipartite states are steerable.Such a quantifier exhibits the

Model and dynamical equations
The hybrid optomechanical system shown in Fig. 1 is composed of a fixed partially transmitting mirror, two charged nano-mechanical oscillators, and an OPA embedded in the cavity.The cavity mode couples the charged nano-mechanical oscillator through the radiation-pressure interaction, and also the first charged nano-mechanical oscillator can be connected to the second spatially separated charged nano-mechanical oscillator by the Coulomb interaction 46 .The cavity is coherently driven by an external laser with frequency ω l and amplitude from the left side of a mirror.The total Hamiltonian of the system can be written as where ĉ(ĉ † ) is the annihilation (creation) operator of the cavity optical mode with cavity frequency ω c .The first term describes the free energy of the cavity field.The second term denotes the energy of the mechanical modes with frequency ω i , effective mass of mechanical oscillators m i and the momentum (position) pi (q i ) .We define the left mechanical oscillator as mechanical oscillator − 1 and the right mechanical oscillator as mechanical oscil- lator − 2 .The third term is the energy of the photon-phonon interaction between the cavity mode and the mechan- ical oscillator − 1 , with the single photon optomechanical coupling constant g = ω c L 2m 1 ω 1 .The fourth term describes the cavity field driven by an input field with frequency ω l and amplitude � = 2κP ω l , where κ and P are the cavity damping rate and the input laser power, respectively.The fifth term denotes the energy between the OPA and the cavity field, G a is the nonlinear gain of the OPA, and θ is the phase of the optical field driving the OPA.The last term shows the Coulomb interaction potential of a charged mechanical oscillators is given as , where C i is the gate capacitance, V i is the voltage of the bias gate, ǫ 0 is the vacuum dielectric constant and r 0 is the separation of the equilibrium positions of the two mechanical oscillators.
We assume the distance between the two charged mechanical oscillators is much greater than the small oscillations of the charged mechanical oscillators ( q 1 , q 2 ≪ r 0 ) , the term describing the interaction between two charged mechanical oscillators can be expanded to the second order as 47 where the linear terms are absorbed in equilibrium positions, and quadratic terms are incorporated in the renormalization of oscillation frequencies.After omitting the constant term, the Coulomb interaction can be reduced to this form 48,49 where is Coulomb coupling strength.The Hamiltonian of the system in a rotating frame at the frequency ω l takes the form (1) where � = ω c − ω l is optical cavity detuning.Using Eq. ( 4) and considering the damping and noise terms into Heisenberg equations, we can obtain the quantum Langevin equations of the hybrid system as where γ 1 (γ 2 ) is the damping rate of the two mechanical oscillators and ĉin is the input vacuum noise operator with zero mean value and nonzero correlation function 26,27 .The quantum Brownian noise operators ξ 1 (t) (ξ 2 (t) ) a rise from the coupling between the two mechanical oscillators with the environment, and their mean values are zero and correlation function where κ B is the Boltzmann constant and T is the temperature of the reservoir related to the mechanical oscilla- tors.However, quantum effects are revealed just for the oscillators with a large quality factor Q i = ω i γ i ≫ 1 .In this limit, the correlation function of the noise ξ i (t) can be written as 50 In which, n = (exp( ω i κ B T − 1)) −1 is the mean thermal phonon number.We utilize the nonlinear quantum Lan- gevin equations for the optical mode and the two mechanical oscillators by taking into account the condition of intense laser driving and weak coupling, i.e., ω l ≫ κ ≫ g 51 .We use the linearization approach, expanding each field operator as the sum of its steady-state mean values and fluctuation operator, which can be treated separately as ĉ = c s + δĉ, qi = q s i + δ qi , qi = q s i + δ pi , where c s , q s i and p s i are the mean values for operators ĉ, qi and pi .The steady-state mean values of the operators can easily be obtained by setting all time derivatives equal to zero in Eqs.(5-9) such that the steady-state mean value for the hybrid system can be obtained where ′ = − gq s 1 is the effective cavity detuning from the frequency of the input laser in the presence of the radiation pressure.In order to examine the quantum correlation between the two mechanical oscillators, we need to calculate the fluctuations of their corresponding operators.The cavity is intensively driven with a very large input laser power, which means that at the steady state, the intracavity field has a large amplitude, i.e., |c s | ≫ 1 .Under, the strong driving limit, here we have neglected the high-order small terms of the fluctuation part, the linearized quantum Langevin equations can be written as We choose the phase reference of the cavity field c s is real.In order to study the quantum statistical properties of the system through the small fluctuations of the system around the steady-state regime.Specifically, the full information about the correlations and entanglement can be characterized by analyzing the variance between the quadrature components of the fields.For this purpose, we define dimensionless quadrature operators.We define the quadrature operators for each mode as Vol.:(0123456789) Thus, we obtain the linearized quantum fluctuations equations where G 0 = √ 2gc s is the effective optomechanical coupling.The equations of motion for the quantum fluctua- tions from Eqs. ( 19)-( 24) can be written in a compact form where ŷin ) T are the column vector of the fluctuation and the column vector of the noises sources, respectively.The drift matrix A can be defined by Therefore, the drift matrix A of Eq. ( 26) contains all the information about the system.It is worth noting that the system can achieve a stable steady-state condition when all of the real parts of the eigenvalues of the drift matrix A are negative.The stability condition can be obtained by using the Routh-Hurwitz criterion 52 .Therefore, the steady state of the quantum fluctuations is a continuous variable Gaussian state.This state is fully characterized by a 6 × 6 covariance matrix (CM) with corresponding components defined as Accordingly, we can express the above as where f (t) = exp(At) , and as t → ∞ , the system is stable.
is the matrix of the stationary noise correlation functions.As a consequence, and using the fact that the components of n(t) are uncorrelated, and using Eq. ( 11), then we get , where D = diag(0, γ 1 (2n + 1), 0, γ 2 (2n + 1), κ, κ) is the noise correlation matrix.The solution of Eq. ( 28) becomes R = ∞ 0 dtf (t)Df (t) T .The stability conditions of the systems are satisfied, then the steady-state correlation matrix can be derived by considering the Lyapunov equation 53 The correlation matrix R can be written in the form of a block matrix where each block represents a 2 × 2 matrix.The blocks on the diagonal represent variance within each subsys- tem (the optical cavity mode, the two mechanical modes ), whereas the blocks off the diagonal represent the correlations between subsystems.

Quantification of quantum correlations
In this section, we measure the quantum correlation between the subsystems.To this aim, it is commonly accepted that this quantum correlation transfer provides a potential tool to exploit the quantum information encoded in mechanical modes that can be more resilient against decoherence effects.The enhancement of this transfer is crucial.In this sense, we measure the quantum correlation between the subsystems through three quantum quantifiers such as quantum entanglement, quantum steering, and quantum discord.Particularly, the (19)   www.nature.com/scientificreports/logarithmic negativity E N is a witness of the entanglement between the bipartite subsystems in a continuous variable system.It is defined by 54,55 which ν s is the smallest partial transposed symplectic eigenvalue χ and given by where = detK m 1 + detK m 2 − 2detL m 1 m 2 , and the correlation matrix is associated with the selected bipartite, and by neglecting the rows and columns, we obtained the interesting mode from Eq. ( 30) in the form of 2 × 2 , block matrix Furthermore, we can see that the necessary and sufficient condition for the Gaussian state being entangled if ν s < 0.5 , which is entirely identical to Simon's criterion, which states that the necessary and sufficient condition for entanglement of non positive partial transpose condition for Gaussian states 56 .Consequently, we numerically describe the results of those calculations via plotted in Figs. 2, 3, 4 and 5. Furthermore, quantum steering is another quantum correlation quantifier that is an essential resource in several fundamental applications 12,13 .
It is a measure of asymmetric property between two entangled observers (between the mechanical oscillator − 1 and mechanical oscillator − 2 in this case).Moreover, it provides a way to quantify how much steerability the two separate mechanical oscillators possess.If we consider the information transfer between the mechanical oscillator − 1 and mechanical oscillator − 2 due to their correlation and label Alice (A: mechanical oscillator − 1 ) and Bob (B: mechanical oscillator − 2 ), we use the covariance matrix of mechanical oscillators of Eq. ( 33) and the steerability of mechanical oscillator − 2 by mechanical oscillator − 1 A → B is defined as 14 where vB is the symplectic eigenvalues of derived from the Schur com- plement of K m 1 in the covariance matrix χ .The measure of Gaussian quantum steering of mechanical oscillator − 2 by mechanical oscillator − 1 is given by The corresponding measure of the Gaussian steerability B → A can be found by From the above, they are two possibilities for quantum steering between A and B: If S A→B = S B→A = 0 , there is no-way steering, which means Alice cannot steer Bob and vice versa even if they are not separable, and two-way steering if S A→B = S B→A > 0 .In actuality, a non-separable state is not always a steerable state, while a steerable state is always not separable.www.nature.com/scientificreports/Finally, we use another quantum correlation quantifier which is a fundamental notion allowing for the description of the quantumness of the correlations present in the state of a quantum system.In our case, the Gaussian quantum discord denotes non-classical correlations if the mechanical oscillator − 1 and mechanical oscillator − 2 are separable or not.Employing Eq. ( 33), the Gaussian quantum discord for the mechanical oscil- lator − 1 defined as 57 The simplectic eigenvalues given by with � ′ = det(K m 1 ) + det(K m 2 ) + 2det(L m 1 m 2 ) and ǫ is defined by Similarly, the quantum discord for the mechanical oscillator − 2 can be found as where (37) Other parameters are the same as Fig. 2.
where the function f is defined by The quantum state of the mechanical oscillator − 1 and mechanical oscillator − 2 are not separable if

Results and discussion
In this section, we investigate the quantum correlation through an optical parametric amplifier and coulombtype interaction in the hybrid system.This can be understood analytically by inspecting the structure of the drift matrix A of Eq. ( 26) which contains all the information about the system.Specifically, we have employed a three-quantum correlation quantifier's expression of Eqs. ( 31), (35), and (37).Moreover, to obtain the covariance matrix, we numerically solve Eq. ( 29) and numerically calculate the logarithmic negativity, which is used as the witness of quantum entanglement, and Gaussian quantum discord that gives the measure of all non-classical correlations.Furthermore, we numerically calculate the quantum steering to characterize the steerability between the two mechanical modes.We now investigate the properties of the quantum correlation of the hybrid system.
( The remaining set of parameters is the same as Fig. 2.
In Fig. 2, we report the quantum correlations measuring between the two mechanical oscillators using logarithmic negativity E N , quantum steering ( S A→B : S B→A ), and quantum discord ( D A : D B ) as a function of the normalized detuning �/ω m for different values of nonlinear gain of the OPA (a) G a = 0 , (b) G a = 0.26κ , (c) G a = 0.35κ , and (d) G a = 0.4κ .In Fig. 2a-d, we numerically display the effect of the nonlinear gain of the OPA on the quantum entanglement, quantum steering, and quantum discord between the two mechanical oscillators.We generate the bipartite entanglement, quantum steering, and quantum discord between the two mechanical modes as functions of normalized detuning for different values of the nonlinear gain of the OPA.In particular, using realistic parameters for which a significant amount of quantum correlation is achievable.Specifically, one can observe that a correlation exists between the two oscillators, implying that there is a quantum correlation between them, even though they are separated.Furthermore, Fig. 2a-d, shows one can clearly understand that the higher the nonlinear gain the coupling parameter is the stronger the oscillators entangle and the broader the range of the correlation between the subsystems can be realized.This shows that increasing the nonlinear gain of the OPA enhances quantum entanglement, quantum steering, and quantum discord compared to the absence of the OPA, i,e., under ordinary light driving.Our result shows that if there is no Coulomb coupling, it is impossible to entangle the two separated oscillators.This is because increasing the nonlinear gain of the OPA corresponds to the increase in the photon number in the optical cavity, which leads to a stronger radiation pressure acting on the mechanical oscillator − 1 and the Coulomb coupling between the mechanical oscillators.Furthermore, our results are consistent with those reported in 45 .
Next, we explore the crucial role of Coulomb coupling strength η on the quantum correlations, e.g., quantum entanglement, quantum steering, and quantum discord between the two mechanical oscillators separated in space.Figure 3, we plot the logarithmic negativity E N , quantum steering ( S A→B : S B→A ), and quantum discord ( D A : D B ) as function of the normalized detuning �/ω m for different value of Coulomb coupling strength (a) η = 0.95ω m , (b) η = 0.97ω m , (c) η = 0.98ω m , and (d ) η = 0.99ω m .As illustrated in the previous section, as long as the logarithmic negativity which characterizes the entanglement remains positive, there is an entanglement between the oscillators, meaning that there is a quantum correlation between two mechanical oscillators, even though they are separated in space.As can be seen from Fig. 3a-d that the presence of strong Coulomb coupling enhances the entanglement, quantum steering, and quantum discord of mechanical oscillators separated by space.It is worthwhile to point out that the larger the coulomb coupling results the more strongly entangled the mechanical oscillators.Furthermore, if there is no Coulomb coupling strength, it is impossible to entangle the two oscillators separately.Therefore, the Coulomb interaction between the two oscillators is the most essential parameter to realize the state transfer between the subsystems.Furthermore, our results are consistent with those reported in 40,45 .Now, we consider the feasibility of the choice of the numerical value of the coupling strength η in the experiment.If we apply the reported experimental parameters, i.e., the gate voltage V 1 = V 2 = 200V , the capacitance of the gate C 1 = C 2 = 2.4nF and the separation between mechanical oscillators without coulomb and optomechanical interaction r 0 = 160µm 40,47 , in this situation η ≈ 0.33ω m .If we compare the numerical values used in our coupled optomechanical system, it is obvious that our choice of the numerical value of coulomb coupling strength is easily executable in experiments.Thus, our hybrid optomechanical system can be realized by choosing appropriate experimental parameters from 58,59 .The results show that the presence of OPA and strong Coulomb coupling enhances the quantum correlations between the two mechanical oscillators, and Coulomb interactions are more prominent in quantum correlations.
We next examine the effect of the phase of the optical field driving the OPA on quantum correlations between the two mechanical oscillators separated in space.In Fig. 4, we plot the quantum correlations, e.g., logarithmic negativity E N , quantum steering ( S A→B : S B→A ), and quantum discord (D A : D B ) as a function of the normalized detuning �/ω m for different values of a phase of the optical field driving the OPA (a) θ = π/16 , (b) θ = π/8 , (c) θ = π/4 , and (d) θ = π/3 .As shown in Fig. 4a-d, it can be seen that the entanglement, quantum steering, and quantum discord larger with the phase θ increases with fixed values of the nonlinear gain of the OPA (G a ) and Coulomb coupling strength η .Because the proper choice of a phase of the optical field driving the OPA may lead to maximum noise suppression, thereby resulting in maximum entanglement, quantum steering, and quantum discord between the two mechanical oscillators.Besides, the optimum values of quantum steering and quantum discord are occurred at �/ω m = 0.13 , and �/ω m = 0.05 , respectively, for the case θ = π/3 .Further- more, the quantum correlation between the two mechanical oscillators increases as the phase of the optical field driving the OPA ( θ ) increases.Thus, we deduced that the presence of phase fluctuation of driving fields affects the quantum correlations.
It is also significant to study the robustness of quantum correlations between the two mechanical oscillators against temperature.Figure 5 shows plots of quantum correlations, e.g., logarithmic negativity ( E N ), steering ( S A→B : S B→A ), and quantum discord (D A : D B ) as a function of the temperature for different values drive power (a) P = 50 mW, (b) P = 70 mW , (c) P = 80 mW , and (d) P = 100 mW with � = 0.6ω m .It is obvious that quantum correlations, such as entanglement, quantum steering, and quantum discord, become stronger as the laser power increases.The reason for this is that increasing the laser power causes a stronger coupling between the mechanical oscillator − 1 and the cavity field due to an increase in the photon number in the cav- ity.The amount of entanglement, quantum steering, and quantum discord also decreases monotonically as the temperature increases because thermal noise in the environment induces decoherence.To be more specific, the entanglement degrades to zero, this is the phenomenon of entanglement of sudden death 60 .According to these results, an increase in temperature leads to the transition from quantum to classical regimes, which is caused by thermal fluctuations.As a result, the hybrid system does not exhibit any entanglement in classical regimes, despite the presence of quite strong laser power.We deduced that the stronger the thermal noise, the higher the temperature of the environment.The entanglement between two separated oscillators is then submerged by the strong thermal noise 50 .We note that the two-way steerable state vanishes with higher values of the environmental temperature.In addition, we obtain better quantum discord at low temperatures, and even at higher values of temperature, quantum discord is found to persist, although entanglement E N vanishes completely, indicating that quantum discord extends beyond entanglement and confirms the robustness of this measure against the fluctuations of the bath environment.Therefore, thermal fluctuations affect quantum correlation, while laser power enhances quantum correlation.We believe that these are appropriate measures to quantify quantum entanglement, quantum steering, and quantum discord and show that our proposed scheme enhances the quantum correlation and proves robust against fluctuations in the bath environment.
Furthermore, we can summarize Figs. 2, 3, 4 and 5 show that entanglement, quantum steering, and quantum discord all behave in the same way.As can be seen, quantum steering is bounded by quantum entanglement.We have S A→B = S B→A > 0 and logarithmic negativity E N > 0 is the witness of two-way steering, while for S A→B = S B→A = 0 and E N > 0 , the two mechanical oscillators are not steerable (i.e.no-way steering).Thus, the quantum discord is more robust than entanglement as shown in Fig. 2, 3, 4 and 5, because when 0 ≤ D A < 1 : 0 ≤ D B < 1 (i.e. the two mechanical oscillators are separable (if E N = 0 ) or entangled (if E N > 0 ).On the other hand when D A > 1 : D B > 1 (i.e. the two mechanical oscillators must be entangled) and we have E N > 0 .Moreover, we can also see from Figs. 2, 3, 4 and 5, that the quantum discord is more dominant than entanglement and also is a good quantifier of quantum correlation.The presence of OPA and strong Coulomb coupling enhances the quantum correlations between the two mechanical oscillators, and Coulomb interactions are more prominent in the nano-electro-optomechanical system.Such a phenomenon is because the enhancement of the effective coupling accelerates the quantum correlation in hybrid systems.Thus, entanglement detection is still a challenge experimentally, but quantum correlation detection is relatively easy.At present, some promising schemes have been suggested in 26,45 , so we can employ homodyne measurement techniques indirectly to detect quantum correlation especially the quantum entanglement 41,61 .Therefore, we believe that our scheme will be experimentally feasible with the quantum information experimental applications, and demand for the development of skills in quantum state manipulation.

Conclusions
In conclusion, we have studied the quantum correlation in a nano-electro-optomechanical system enhanced by an optical parametric amplifier and Coulomb-type interaction.We consider a hybrid system that comprises a cavity and two charged mechanical oscillators with an OPA, in which the cavity mode is coupled to a charged mechanical oscillator via radiation pressure, and the two charged mechanical oscillators are coupled via Coulomb interaction.We showed that the Coulomb interaction between mechanical oscillators is the primary reason for the existence of a quantum correlation between the two mechanical oscillators.Our result shows that the presence of OPA and strong Coulomb coupling enhances the quantum correlations between the two mechanical oscillators.In addition, Coulomb interactions are more prominent in quantum correlations.Besides, in the presence of OPA, the maximum amount of quantum entanglement, quantum steering, and quantum discord achieved between the two mechanical oscillators than in the absence of OPA.This is because increasing the nonlinear gain of the OPA increases the photon number in the cavity, which leads to a stronger radiation pressure acting on the mechanical oscillator − 1 .Additionally, we show that a proper phase choice of the optical field driving the OPA enhances quantum correlations under suitable conditions.This is because the proper choice of a phase of the optical field driving the OPA may lead to maximum noise suppression, thereby resulting in maximum quantum correlations.Furthermore, the quantum correlations decline rapidly with increasing temperature as a result of decoherence.Specifically, we noted that the quantum entanglement degrades to zero; this is the phenomenon of the quantum entanglement of sudden death.According to these results, an increase in temperature leads to the transition from quantum to classical regimes, which is caused by thermal fluctuations.Furthermore, when compared to quantum entanglement, the two-way steerable state vanishes with higher values of the environmental temperature.Finally, we have shown that quantum discord persists at higher temperature values, even though quantum entanglement disappears completely at higher temperature values.In this regard, quantum discord extends beyond entanglement and confirms the robustness of this measure against the fluctuations of the bath environment.Our proposed scheme of quantum correlation provides a promising platform for realizing continuous variable quantum information processing.

Figure 1 .
Figure 1.Schematic representation of the system.Hybrid OMS consists of the optical cavity, two charged nanomechanical resonators and a degenerate OPA placed inside the cavity, and the pump of the OPA is not shown.

Figure 2 .
Figure 2. Plots of quantum correlations, e.g., logarithmic negativity E N , quantum steering(S A→B : S B→A ), and quantum discord(D A : D B ) as function of the normalized detuning �/ω m for different value of the nonlinear gain of the OPA (a) G a = 0 , (b) G a = 0.26κ , (c) G a = 0.35κ , and (d) G a = 0.4κ with phase of the optical field driving the OPA ( θ = 0 ), temperature ( T = 4 mK), laser power (P = 10 mW) , and Coulomb coupling strength (η = 0.95ω m ) .Other parameters are listed in the main text.

( 34 )Figure 3 .
Figure 3. Plots of quantum correlations, e.g., logarithmic negativity E N , quantum steering ( S A→B : S B→A ), and quantum discord (D A : D B ) as a function of the normalized detuning �/ω m for different value of Coulomb coupling strength (a) η = 0.95ω m , (b) η = 0.97ω m , (c) η = 0.98ω m , and (d) η = 0.99ω m with the nonlinear gain of the OPA (G a = 0.26κ) .Other parameters are the same as Fig. 2.

Figure 5 .
Figure 5. Plots of quantum correlations, e.g., logarithmic negativity E N , quantum steering ( S A→B : S B→A ), and quantum discord ( D A : D B ) as a function of the temperature for different values drive laser power (a) P = 50mW , (b) P = 70mW , (c) P = 80 mW , and (d) P = 100 mW with � = 0.6ω m with the nonlinear gain of the OPA (G a = 0.26κ) .The remaining set of parameters is the same as Fig. 2.
oscillator − 1 and mechanical oscillator − 2 can be in a separable state or an entangled state.Thus, we numerically describe the results of those calculations and plotted in Figs. 2, 3, 4 and 5.