Phase control of entanglement and quantum steering in a three-mode optomechanical system

The theory of phase control of coherence, entanglement and quantum steering is developed for an optomechanical system composed of a single mode cavity containing a partially transmitting dielectric membrane and driven by short laser pulses. The closed loop in the coupling creates interfering channels which depend on the relative phase of the coupling strengths of the field modes to the mechanical mode. We show several interesting phase dependent effects such as reversible population transfer from one field mode to the other, creation of collective modes, and induced coherence without induced emission. These effects result from perfect mutual coherence between the field modes which is preserved even if one of the modes is not populated. Depending on the phase, the field modes can act on the mechanical mode collectively or individually resulting, respectively, in tripartite or bipartite entanglement. In addition, we examine the phase sensitivity of quantum steering of the mechanical mode by the field modes is investigated. Deterministic phase transfer of the steering from bipartite to collective is predicted and optimum steering corresponding to perfect EPR state can be achieved. These different types of quantum steering can be distinguished experimentally by measuring the coincidence rate between two detectors adjusted to collect photons of the output cavity modes. In particular, we find that the minima of the interference pattern of the coincidence rate signal the bipartite steering, while the maxima signal the collective steering.


Introduction
It is well know that optical coherence and quantum entanglement between two systems can occur when the source systems are prepared in a superposition state [1,2]. Especially, the mutual coherence is described by the first-order correlation between two undistinguishable pathways or channels and can determine the phase of a radiating unknown system relative to the phase of the other known system. On the other hand, quantum entanglement determines inseparability of quantum systems, and has been recognized as one of the most intrinsic features of quantum mechanics with many useful applications ranging from quantum cryptography, quantum metrology, to quantum computation [3][4][5]. A particular interest is devoted to a special type of entanglement called quantum steering [6][7][8]. The concept of quantum steering was originally introduced by Schrödinger [9] to explore the fact that entanglement would allow one to remotely steer or pilot the state of a distant system, as considered in the original Einstein-Podolsky-Rosen (EPR) paradox [10]. The EPR steering allows two parties to verify the shared entanglement even if one measurement device is untrusted, which makes it an essential resource for one-sided device independent quantum cryptography [11][12][13][14][15][16], one-way quantum computing [17,18], secure quantum teleportation [19,20], and subchannel discrimination [21].
Recent studies have revealed that coherence is closely related to entanglement and quantum steering. It was pointed out by Suzuki et. al [22] that entanglement can be detected from interference fringes in atom-photon systems. It has also been shown that the coherence in a system and entanglement between that system and another initially incoherent one are quantitatively, or operationally, equivalent [23]. The power of quantum steering for the generation of coherence has also been demonstrated [24,25]. It has been shown that the presence of mutual coherence among two systems may have a distractive effect on the creation of entanglement between these systems [26]. On the other hand, it has been shown that in a tripartite system the mutual coherence between two parties may help to collectively steer the third party [27].
The hope to demonstrate entanglement and quantum steering in macroscopic systems has encouraged research on optomechanical systems.
Apart from the extensive efforts toward the creation of entangled states, it is crucial for quantum information processing to be able to control the creation and evolution of the entangled states. Recent research on phase dependent systems has addressed the problem of controlling the entanglement in two qubits systems [58,59], a triple spin qubit [60], and quantum parametric oscillators [61]. Other simple systems have demonstrated the optical nonreciprocal behaviour induced by the phase difference between the coupling constants in a fully coupled tripartite optomechanical system [62], or in a multimode on-chip electromechanical system [63], and novel quantum interference effects and correlations arising from the interplay and competition between different excitation channels [64].
In this paper we explore the possibility of the phase controlled generation and transfer of coherence, entanglement and quantum steering in a closed loop multimode system. We propose to consider a close loop coupling between modes in optomechanics. In particular, we examine phase dependent dynamics an optomechanical system composed of a single mode cavity containing a dielectric membrane in its interior. We assume that the membrane can behave as a partly transmitting and reflecting mirror which divides the cavity into two mutually coupled optomechanical cavities. We examine phase properties of the coherence functions, the first-order coherence and anomalous correlation functions involved in the generation of entanglement between modes. We then propose to use such phase dependent functions in the controlled generation and transfer of entanglement among the modes. Specifically, we show that the presence of the mutual coherence between the field modes may result in perfect entanglement between the mechanical and one of the field modes. By varying the relative phase of the coupling constants of the field modes to the mechanical mode, one can achieve the perfectly entangled tripartite state involving the mechanical mode and a linear superposition of the field modes. In addition, we demonstrate the crucial role of the phase in the quantum steering of the mechanical mode by the field modes, the achievement of perfectly steerable EPR state and the phase controlled distinguishability of the collective and bipartite steering.
The paper is organized as follows. In section 2 we describe in detail our model and discuss the method of the normalized temporal pulse-shape amplitudes used to solve the appropriate equations of motion describing the evolution of the mechanical mode and the cavity field operators. These equations are then used to find the analytical expressions for the population of the modes and for the first-order correlations and anomalous correlation functions. We apply the expressions in section 3 to determine the mutual coherence between the modes and the phase dependent transfer of the population between the modes. We also present there results concerning the degree of coherence, visibility of the interference fringes, distinguishability of the modes and the phase dependence of the coincidence rate between two detectors adjusted to collect photons of the output modes of the two cavities. Section 4 is devoted to the discussion of the inseparability criterion. We present the analytical form of the separability parameter and discuss in details its phase properties. In section 5 we focus on steering properties of the modes, discuss the possibility of achieving a perfectly steerable EPR state between the modes, and suggest a measurement technique to distinguish between the bipartite and collective steerings. We summarize our results in section 6. Finally, in the Appendix, we give results for the separability parameters evaluated for excitation of the cavity modes with a noisy laser.

The model
We consider a system composed of a single mode cavity containing a dielectric membrane in its interior, as illustrated in figure 1(a). The membrane divides the cavity into two cavities resulting in a system effectively behaving as a three-mode optomechanical system: two field modes of the two cavities plus a mechanical mode representing the vibrating membrane. The field modes and the mechanical mode are treated as quantized: a j (a † j ), j = 1, 2, is the annihilation (creation) operator for the mode of cavity j, and c (c † ) is the annihilation (creation) operator for the mechanical mode. The modes of these two cavities have frequencies ω 1 and ω 2 depending on the position of the membrane relative to the position of the cavity mirrors. The mechanical mode has a frequency ω m . The field modes are driven by external detuned short laser pulses of equal duration time τ , equal frequencies ω L , and electric field amplitudes The laser pulses are derived from the same laser to ensue a constant phase difference δϕ 0 = ϕ 01 − ϕ 02 . The thickness of the membrane is very small that it can be treated as a nonabsorptive partly reflecting and transmitting mirror. The possibility of the field transmission through the membrane results in a direct linear coupling of the modes a 1 and a 2 with a strength J. This coupling together with the coupling to the mechanical mode creates a three-mode closed linkage or "loop" coupling of the modes of the system, as illustrated in figure 1(b). The presence of this closed loop will give rise to interference effects and phase dependence of the dynamics of the modes. The Hamiltonian of this system may be written as where g 0,j is the single-photon coupling of the mode j to the mechanical mode. The Hamiltonian (2.1) may be met in typical experimental optomechanical systems with a membrane in a Fabry-Perot cavity [37].
Using the Heisenberg equation of motion and Hamiltonian (2.1), we obtain the following equations of motion for the annihilation operatorṡ along with the corresponding equations for the creation operators. Here, a r j = a j exp{i [ω L t + ϕ 0j ]}, a in,r j = a in j exp{i [ω L t + ϕ 0j ]}, and ∆ j = ω j − ω L is the detuning of the laser frequency from the frequency of the mode j. In writing equation (2.2) we have included relaxation terms of the modes: κ j is the damping rate of the field mode j, and γ is the damping rate of the mechanical mode. We have also included the noises of the input modes a in j and c in arising from the coupling of the modes to their surrounding environments. Note that the dependence of the equations of motion (2.2) on the phase difference δϕ 0 only arises from the linear coupling of the two cavity modes. For the statistics of the input modes we assume that the input modes of the cavities are in the ordinary vacuum state characterized by the correlation function a in j (t)(a in j ) † (t ) = δ(t − t ), and the input mode to the vibrating membrane is in a thermal vacuum states characterized by the correlation functions c in (t)(c in ) is the average number of thermal phonons at the frequency of the mechanical mode, k B is the Boltzmann constant and T is the temperature of the environment.
Assuming E 0j κ j , γ, the equations of motion (2.2) may be solved by the linearization approach [65]. In this approach, we write the operators of the system as composed of their steady-state mean values and a small fluctuation around the steadystate. The steady-state mean values of the operators, a j ≡ α j and c ≡ χ, are obtained by setting the left-hand sides of equation (2.2) to zero yielding where ∆ j = ω j − ω L , with ω j = ω j + g 0,j (χ + χ * ). Using the expansions a j → α j + δa j and c → χ + δc, we obtain the equations of motion for the fluctuation parts of the operators. Keeping the linear terms only, the equations of motion are of the form where g j = g 0,j α j is the effective optomechanical coupling strength of the mode j to the mechanical mode. It is seen from equation (2.4) that the fluctuation operators oscillate at frequencies ±ω m . Therefore, it is convenient to introduce slowly varying fluctuation operators δc m = δce iωmt , δc in,m = δc in e iωmt , δa r,m j = δa r j e −i(ωmt+ϕ 0j ) , δa in,r,m j = δa in,r j e −i(ωmt+ϕ 0j ) , and substitute them into equation (2.4). Discarding all terms oscillating at 2ω m based on the rotating-wave approximation, we then obtaiṅ Here we have assumed that the laser pulses are tuned to the blue sideband of the average frequency of the field modes, i.e., ω L = ω 0 + ω m , where ω 0 = (ω 1 + ω 2 )/2 and ∆ = (ω 1 − ω 2 )/2. For clarity of the notation, we have omitted in equation (2.5) the symbol δ and the superscripts r and m on the displacement operators. From equation (2.5) we can see that the modes a 1 and a 2 are directly coupled to each other with the strength J, and are also indirectly coupled to each other through the coupling to the mechanical mode with strengths g 1 and g 2 , respectively. This coupling configuration forms a closed loop, as illustrated in figure 1(b), and therefore the dynamics of the system can exhibit phase dependent effects.
Let us first consider the result of the coupling J on the dynamics of the field modes. By introducing column vectors Y a = (a 1 , a 2 ) T , Y in a = (a in 1 , a in 2 ) T , and Y c = (c † , c † ) T , we can put the equations of motion for a 1 and a 2 into a matrix forṁ where the matrixM describes the effects of the linear coupling J between the modes and the detuning ∆, whereas the matrixN governs the influence of the nonlinear coupling to the mechanical mode, Diagonalization of the matrixM results in orthogonal superposition modes represented by the superposition operators a w = a 1 cos θ + a 2 sin θ, a u = a 1 sin θ − a 2 cos θ, (2.8) where The angle θ belongs to the interval [0, π/2].
Using equations (2.5) and (2.8) one can easily find that the annihilation operators of the superposition modes satisfy the following equations of motioṅ where g w = g 1 cos θ + g 2 sin θ, g u = g 1 sin θ − g 2 cos θ are effective coupling strengths of the superposition modes to the mirror mode. Note that the superposition modes a w and a u only couple to the mechanical mode and not to each other. This coupling configuration contains no closed loops, it resembles a chain coupling sequence. However, the dynamics of this system can exhibit phase dependence as it involves phase-dependent coupling constants g w and g u . Furthermore, the dynamics reduce to those involving two modes only if the conditions ∆ = 0 and g 1 = g 2 ≡ g are fulfilled. In this case, g w = √ 2g and g u = 0 that the mode u completely decouples from the remaining modes.
Equation (2.9) is the basic equation for calculating the dynamics of the modes. In the following we always concentrate on transient effects resulting from the excitation of the modes by short laser pulses.
The solution of equation (2.9) is in general complicated. A simple analytical solution arises, however, in two cases, the bad cavity limit, κ |g w,u |, or at a large difference between the frequencies of the superposition modes, w |g w,u |. Under such conditions, we can approximate a w and a u as slowly varying in time, and putȧ w ≈ 0 andȧ u ≈ 0, which leads to where φ = arctan (w/κ). The equation of motion for the mechanical mode becomeṡ where δ = (|g w | 2 − |g u | 2 )w/(κ 2 + w 2 ), G w = |g w | 2 κ/(κ 2 + w 2 ), G u = |g u | 2 κ/(κ 2 + w 2 ), and G = G w +G u −γ. Note that the following results are all obtained in this approximation. Before proceeding further, we briefly comment about the possible setups and the parameters of the proposed model could be realized with the current experiments. A relevant for the possible realisation of our proposal could be, for example, the experimental system reported in reference [66], with a silicon nitride membrane inside of a Fabry-Perot optical cavity. The parameters achieved were ω m /2π = 1.55 MHz, κ/2π = 0.89 MHz, γ/2π = 0.47 Hz, single-photon optomechanical coupling rate g 0 /2π = 16 Hz, and the photon number N a = 3.6 × 10 8 , such that the effective coupling rate was estimated as g/2π = 0.3 MHz. In a similar system reported in reference [67], the parameters values achieved were ω m /2π = 1.48 MHz, κ/2π = 2.6 MHz, γ/2π = 0.18 Hz, and g/2π ≤ 30 kHz. Furthermore, this system was prepared in a dilution refrigerator where n m ∼ 10 3 and could be cooled down to the ground state with n 0 = 0.2. Besides, our proposal could be realized in microwave optomechanics, where LC resonators are used [63,68]. In those experiments, parameters achieved were ω m /2π ≈ 5 MHz, κ/2π ≈ 2MHz, γ/2π ≈ 4 Hz, g 0 /2π ≈ 13 ∼ 34 Hz and the photon number N a ≈ 2350 ∼ 4 × 10 7 such that g/2π = 1.5 ∼ 200kHz. The LC circuit at milliKelvin temperatures in a dilution refrigerator can be cooled down to an average phonon occupancy of n m = 0.32. The linear coupling strength J tuneable in a large range from J ≈ 0.2 MHz to J ≈ ω m (ω m /2π ≈ 5 MHz) have been reported [69,70]. We point out that the computed theoretical results will be graphically presented for the feasible parameters chosen from the above mentioned experiments.
We now proceed to evaluate the populations of the modes and correlations between them. We assume that G > 0 that G is the gain associated with the system. The input noise to the mechanical mode will be amplified during the evolution. Therefore, we shall concentrate on short-time effects. Using the standard cavity input-output relation a out = a in + √ 2κa, we define annihilation operators of normalized temporal field modes [46,52]: where τ is the interaction time with the laser pulses. It is easily verified that the operators obey the canonical commutation relations In what follows, we assume that the input fields to the field modes are in the ordinary vacuum state, the mechanical mode is initially in a thermal state with the mean number of phonons n 0 . The mode is subjected to the damping γ and to the Brownian noise n m . Under these assumptions, we find the populations of the output modes to be and first-order correlation functions Here, r = Gτ represents an effective squeezing parameter, φ gw and φ gu are phase angles associated with the coupling strengths of the superposition modes, g w = |g w | exp(iφ gw ) and g u = |g u | exp(iφ gu ), respectively. Two remarks must be made about the solutions given by equations (2.13) and (2.14). Firstly, there is no first-order coherence between the mechanical and superposition modes, but there is a nonzero anomalous coherence between these modes. Inversely, the first-order coherence is observed between the superposition modes but there is no anomalous coherence between these modes. This means that the first-order coherence between modes rules out the anomalous coherence between these modes and vice versa, the anomalous coherences between modes results in no first-order coherence. Secondly, the phases φ gw and φ gu do not affect measurable quantities such as populations of the modes, the first-order coherence between the modes or variances of the quadrature components of the modes. However, these measurable quantities can still be phase dependent since the absolute values |g w | 2 and |g u | 2 are not constant parameters. In the presence of the coupling J these parameters turn out to be dependent on the phases of g 1 and g 2 . It is easy to see. Writing g 1 = |g 1 | exp(iφ g 1 ) and g 2 = |g 2 | exp(iφ g 2 ), the phase-sensitive contribution to |g w | 2 and |g u | 2 can be identified, and the parameters |g w | 2 and |g u | 2 may be expressed as where 2ψ = φ g 1 − φ g 2 is the relative phase of the g 1 and g 2 coupling strengths. Clearly, the parameters |g w | and |g u | vary periodically with 2ψ, but only if sin 2θ = 0. For J = 0 we generally have sin 2θ = 0, and then phase dependent effects are to be expected. On the other hand, for J = 0, we have sin 2θ = 0 and then there is no phase dependence. Referring to equation (2.3) the phases φ g 1 and φ g 2 can be controlled through the phases of the driving lasers and a constant phases difference 2ψ can be established. In an experimental situation, we would envisage a single laser providing the pump for both modes.
We may summarize that the presence of the linear coupling J between the field modes results in the phase dependence of the populations of the modes and coherences between them. In other words, the dependence of the solutions (2.14) on the relative phase 2ψ results from the presence of the three-mode loop in the coupling between the modes, as illustrated in figure 1(b). This close loop coupling leads to interesting new effects, which will be discussed in details in the following section.

Phase control of mutual coherence and populations of the optical fields
We now turn to analyse the effect of the correlation between the superposition modes on the correlation and mutual coherence of the output field modes A out 1 and A out 2 . To do this, we invert the transformations (2.6) and find Then, defining annihilation operators of the normalized temporal output modes we readily find that the relationship between A out 1 , A out Then the populations of the modes and correlation functions are found to be and we have extracted the phases φ A j of the complex amplitude A j (ψ) = |A j (ψ)| exp(iφ A j ). The solutions (3.4) are of a simple form with the terms A 1 (ψ) and A 2 (ψ) representing the interference effects, while the terms Υ(r) and Λ(r) give the time-dependent effect of the driving laser pulses. It is interesting to note that the interference terms factorize from the temporal terms that the interference effects are independent of the duration of the laser pulses. The solutions (3.4) have the same type correlation properties as in the case of the superposition modes, equations (2.13) and (2.14). However, the dependence on the phases θ and ψ is more complicated than that for the superposition modes and, in addition, it involves the phase φ.
It should be noted here that for a given r the total population of the modes is constant, independent of the phases θ, φ, and ψ, In other words, the variation of the populations with the phase ψ is due to the transfer of the population from one mode to the other not due to the generation or losses of photons.

Mutual coherence and population transfer
Let us discuss properties of the solutions (3.4) for some special cases. We are particularly interested in the properties of the mutual coherence between the output field modes. Firstly, we note that for J = 0 we have sin θ = 0, cos θ = 1, and then the coherence Clearly, in the absence of the coupling between the field modes (J = 0), the mutual coherence | is independent of the phases. When the modes are degenerate (∆ = 0), we have sin 2 θ = cos 2 θ = 1/2 and then we obtain Further simplification to the case of the symmetric coupling |g 1 | = |g 2 | = g gives | vanishes only for sin 2 2φ = 1, i.e., for φ = π/4 corresponding to J = κ. For J = κ, the mutual coherence never vanishes.
It is interesting to noted that (A out 1 ) † A out 2 can vanish even for an asymmetric coupling |g 1 | = |g 2 |. The expression (3.8) is a sum of two positive numbers and therefore it would vanish when simultaneously both numbers are equal to zero. It is seen that both numbers are simultaneously equal to zero when cos 2ψ = 0, sin 2ψ = ±1, and then tan 2φ = ∓2|g 1 ||g 2 |/(|g 1 | 2 − |g 2 | 2 ). Thus, in the case of unequal coupling strengths the coherence function can vanish but only for a specific value of the ratio J/κ that can satisfy this condition.
Apart from the mutual coherence, the populations of the superposition modes A out w and A out u , as well as the modes A out 1 and A out 2 can depend on the phases. It is not difficult to see from equations (2.13) and (2.16) that the populations of the superposition modes in the case |g 1 | = |g 2 | never become zero as the phase ψ is varied. For example, in the case ψ = 0, the minimum of the population is simply proportional to (|g 1 | sin θ−|g 2 | cos θ) 2 /2, and the maximum is proportional to (|g 1 | cos θ + |g 2 | sin θ) 2 /2. Only in the limit of |g 1 | = |g 2 | = g, the population may be completely transferred between the modes. In this case, the populations are given by (1 + sin 2θ cos 2ψ) , We see that in the absence of the linear coupling (sin 2θ = 0) the modes are equally populated and the populations are independent of the phase ψ. The presence of the coupling J clearly results in the populations becoming dependent on the phase ψ and thus to allow the transfer of the population between the modes. However, the total transfer of the population takes place only when sin 2θ = 1, i.e., in the degenerate case of ∆ = 0. Otherwise the populations can vary with the phase ψ but cannot be totally transferred between the modes. For θ = π/4 and ψ = nπ, (n = 0, 1, 2, . . .), we see from equation (3.10) that the mode u is unpopulated for all values of r, whereas the mode w is maximally populated. On the other hand, for ψ = (n + 1/2)π this relationship is reversed and u is the mode which is maximally populated. Consider now the populations of the output field modes, A out 1 and A out 2 . From equations (3.4) and (3.5) we see that apart from the dependence on the phase angles θ and ψ, the populations depend on the phase angle φ. This may lead to some limits in the transfer of the populations not present for the superposition modes. It is easy to see from equation (3.5) that for J = 0 at which sin 2θ = 0, sin 4θ = 0, the populations are independent of the phases φ and ψ. On the other hand, for ∆ = 0 and J = 0, we have sin 2θ = 1, sin 4θ = 0 and then In this case the populations depend on the phases and are also strongly dependent on the relationship between the coupling constants |g 1 | and |g 2 |. The involvement of the phase φ allows to achieve the complete transfer of the population between the states even if |g 1 | = |g 2 |. In particular, for the symmetric coupling |g 1 | = |g 2 | = g, we get Viewed as a function of ψ, the population can be transferred between the modes reversibly. However, there is a notable limit in the transfer of the populations provided by phase angle φ. Only for φ = π/4 the population can be totally transferred from one mode to the other.

Degree of coherence and visibility
We have already seen that the coherence functions (A out w ) † A out u and (A out 1 ) † A out 2 vary with the phase angles θ, φ and ψ, and may vanish for particular choices of the phases. However, if we calculate the the first-order coherence γ (1) ij between the output modes, defined as we find that the first-order coherences γ wu and γ 12 are always unity, i.e., γ wu = 1 and γ (1) 12 = 1. This means that the superposition modes as well as the output cavity modes , with the phase angle ψ for equal (|g 1 | = |g 2 |) and unequal (|g 2 | > |g 1 |) couplings of the field modes to the mechanical mode. It is clearly seen that the firstorder coherence function and the populations vary periodically with the phase ψ. The coherence function is greatest half-way between the zeros of the populations and there are zeros in the coherence for phases at which the population is completely transferred to one of the modes. While figure 2(a) shows that in the case of |g 1 | = |g 2 | the population is periodically transferred from one mode to the other, figure 2(b) shows that in the case of |g 1 | = |g 2 | the complete transfer of the population occurs only to mode with larger coupling strength.
It is interesting that regardless of the distribution of the population between the modes the degree of coherence γ (1) jk = 1. The constant phase relation between the modes is preserved even if one of the modes is not populated. This surprising behaviour is an example of induced coherence between the modes without induced emission, a phenomenon discussed and observed by several authors [71][72][73][74][75][76][77][78]. In fact, the system considered here of two cavity modes coupled parametrically to the mechanical mode is analogous to the system of Wang, Zou and Mandel [71,72] composed of two coupled parametric downconverters in which the induced coherence without induced emission was demonstrated experimentally. In their experiment two downconverters were pumped by the same laser and arranged in cascade, that the idler field emitted from one of the downconverters was used as the input field of the other downconverter. Under this arrangement interference effects were observed between the signal fields of the two downconverters. It was further shown that in their system the degree of coherence γ (1) jk = 1. In our system the two cavity modes are parametrically coupled to the mechanical mode and therefore can be treated as two downconverters. Since the modes are linearly coupled to each other through the coupling to the mechanical mode, they are locked in phase also and therefore exhibit perfect first-order coherence.
Although the modes A out 1 and A out 2 are mutually perfectly coherent, the visibility of the interference pattern can be zero that the modes can be completely distinguishable. The degree of visibility V of the interference pattern is related to distinguishability D through the expression |V| 2 + |D| 2 ≤ 1 of complementarity. When |V| = 1 then |D| = 0 which means that the modes are indistinguishable. On the other hand, when |V| = 0 then |D| = 1 that the modes are perfectly distinguished. The visibility |V| is given by (3.14) Similarly, we may define a distinguishability Note that |V| 2 +|D| 2 = 1, which is a special case of the general expression |V| 2 +|D| 2 ≤ 1.
If we take ∆ = 0, we then find after substituting equation (3.5) into equation (3.15) that |D| = | (|g 1 | 2 − |g 2 | 2 ) cos 2φ − 2|g 1 ||g 2 | sin 2φ sin 2ψ| It is not difficult to check that the conditions for the modes to be perfectly indistinguishable (D = 0) are different from the conditions for the modes to be perfectly distinguishable (|D| = 1). For example, the modes are perfectly indistinguishable for the phase ψ such that sin 2ψ = ±1 if tan 2φ = ± (|g 1 | 2 − |g 2 | 2 ) /2|g 1 ||g 2 |. On the other hand, the modes can be perfectly distinguishable at the same phase ψ only if |g 1 | = |g 2 | and sin 2φ = 1, the latter happens when J = κ. Thus, for J = κ the modes are always at least partly indistinguishable. The phase dependent transfer of the population between the output cavity modes can be experimentally observed by measuring the coincidence rate R 12 between two detectors D 1 and D 2 adjusted to collect photons of the output modes A out 1 and A out 2 , respectively. The coincidence rate is proportional to the second order correlation function. In the case of the two output cavity fields, the rate is given by The choice of ∆ = 0 and |g 1 | = |g 2 | ≡ g leads to the result is zero. In other words, R 12 = 0 signals the complete transfer of the population to one of the output cavity modes.
In closing this section, we would like to point out that nonzero anomalous correlations are present between the mechanical and field modes, i.e., A out m A out j , (j = 1, 2, w, u) are all different from zero. These correlations are known to be responsible for entanglement between two modes. At the same time, all the first-order coherence functions (A out m ) † A out j are zero. It then follows that the mechanical and the field modes are mutually incoherent. Therefore, we may conclude that interference effects between two modes signal the complete separability of the modes and vice versa, entanglement signals the modes are mutually incoherent. In this sense, pairs of modes which are mutually coherent are clearly distinguishable from those which are entangled.

Phase control of entanglement
The presence of phase dependent populations of the modes and anomalous correlations between the mechanical and the field modes suggests that it could be possible to control and transfer of entanglement between these modes by manipulating the phase. To see the effects of the phase on entanglement between the modes, we examine the Duan-Simon inseparability criterion which involves linear combinations X − X, or P − P or X − P of the quadrature components of the modes.
The solutions given in equation (2.10) show that the X j quadrature component is linearly related to the P m quadrature component, and vice versa, the P j quadrature component is linearly related to the X m quadrature component. Therefore, to quantify entanglement we will use the inseparability criterion for asymmetric X −P combinations of the quadrature components of the output fields, which is determined by the separability parameter ( = 1) [47,[79][80][81] where h j is a weight factor which is chosen to minimize the variances. The optimal value of h j which minimizes ∆ i,j is found using the variational method, Here the variances is defined as The modes A out m and A out j are said to be entangled if ∆ m,j < 1. The inequality ∆ m,j < 1 is necessary and sufficient condition to confirm entanglement between modes in Gaussian states under Gaussian measurements. The entanglement will be said to be maximal or perfect EPR state is created between the modes A out m and A out j if ∆ m,j can become zero.
To determine the inseparability criterion we need the variances of the quadrature components of the modes and correlations between them. An obvious question is how to choose in our phase dependent system proper phase dependent quadratures which would correspond to the maximal reduction of fluctuations. The structure of the solutions given in equation (2.10) suggests to define quadrature components of the optical fields that are in phase and out of phase with the coupling strengths of the modes to the mechanical mode. With this identification, we define the following quadrature components for the output field of the mechanical mode, for the output field of the superposition modes (j = u, w), and for the output field of the cavity modes (j = 1, 2). For the input field modes, we define the following quadrature components For the field modes initially in the ordinary vacuum state the variances of the quadrature components of the input fields are ∆ 2 X in 1 = ∆ 2 P in 1 = ∆ 2 X in 2 = ∆ 2 P in 2 = 1/2. For the mechanical mode initially in a thermal state, ∆ 2 X in m = ∆ 2 P in m = n 0 + 1/2. To determine variances of the output field quadratures and correlations between them, which are needed in equation (4.1), we make use of equations (2.14) and (3.4), and find where |U j (ψ)| =    κ G(κ 2 +w 2 ) |g j (ψ)|, for j = u, w, κ G(κ 2 +w 2 ) |A j (ψ)|, for j = 1, 2. (4.8) Note that the variances of the quadratures of the output fields are all greater than 1/2. This means that the fluctuations in the output modes are not squeezed and therefore correlations between the modes are necessary to produce entanglement.
To illustrate the analytic structure of the parameter ∆ m,j , we consider two simpler but standard cases where one case is the minimized ∆ m,j by the optimal gain factor with the limit γ/G 1 [66] and the other case is ∆ m,j | h j =1 for a fixed gain factor h j = 1. After making use of equation (4.7) in equation (4.1), we find that in the limit of γ/G 1, the parameter ∆ m,j is given by It is clear from the form of the two terms on the right-hand side of equation (4.9) that the separability parameter ∆ m,j attains a minimum value, corresponding to maximum entanglement, when the term inside the curl brackets is equal to zero. It happens when e −2r ≈ 0 and h j |U j (ψ)| = 1, in which case ∆ m,j = (h 2 j − 1)/(h 2 j + 1). It follows that optimal entanglement (∆ m,j = 0) could be observable in principle for the symmetric (h j = 1) combination of the X m and P j quadrature components and sufficiently long laser pulses.
However, the question whether ∆ m,j < 1 and whether it is possible to achieve ∆ m,j = 0 depends strongly on the phase dependent factors |U j (ψ)| and the optimal gain factor h j . The variation of the parameters ∆ m,j with the relative phase ψ for large r is plotted in figure 3. The horizontal dashed-dotted line represents the degree of entanglement in the absence of the coupling J. In this case both field modes are equally entangled with the mechanical mode and the entanglement is independent of ψ. In fact, ∆ m,j can never be reduced to zero even for large τ due to the presence of initial thermal noise in the mechanical mode. In the presence of the coupling J, the separability parameters oscillate periodically with the phase. The amplitude of the oscillations is equal to one. The most interesting is that optimal entanglement (∆ m,j = 0), resulting in perfect EPR state, is achieved at sin 2ψ = ±1 for j = 1, 2 and cos 2ψ = ±1 for j = w, u. For example, the parameter ∆ m,1 = 0 happens periodically for the phase satisfying ψ = (n + 3/4)π, (n = 0, 1, 2, . . .) where the optimal gain h 1 = 1 and the phase dependent factor |U 1 (ψ)| = 1. Meanwhile, ∆ m,2 = 1 at these values of the phase since there is no population in the mode A out 2 . On the other hand, ∆ m,w = 0 happens for the phase ψ satisfying ψ = nπ, (n = 0, 1, 2, . . .) where the optimal gain h w = 1 and |U w (ψ)| = 1. Meanwhile, ∆ m,u = 1 at these phases since there is no population in the mode A out u . Note that the case ∆ m,j = 0 occurs for phases ψ at which either A out 1 and A out 2 or A out w and A out u are distinguishable. In addition, the system behaves deterministically that the entanglement can be periodically transferred between the modes by varying the relative phase ψ. This is of course a reflection of the fact that the field modes are perfectly coherent and the coherence is preserved independent of the redistribution of the population between the modes. Note that the location of the zeros of ∆ m,j corresponds to the zeros of the population of the other field mode (see figure 2). Thus, a constant phase relation between the modes allows for coherent transfer of the population and correlations between the modes.     Figure 4 shows the corresponding situation for unequal coupling strengths, |g 1 | = |g 2 | with the degenerate modes ∆ = 0 and γ = 0. We see that the amplitudes of the periodic maxima and minima of ∆ m,j occur at the same phases as in figure 3, but are noticeable reduced in magnitude. Figure 4a shows that the amplitude of ∆ m,2 can be reduced to zero at ψ = (n + 1/4)π, (n = 0, 1, 2, . . .) where h 2 = 1, |U 2 (ψ)| = 1, while ∆ m,1 cannot due to the required condition cannot be satisfied. This is easy to understand if one refers to the fact that in the case |g 2 | > |g 1 | only the population of the mode A out 1 can be completely transferred to the mode A out 2 . When the population is completely transferred ∆ m,1 = 1 and simultaneously ∆ m,2 = 0. Figure 4b shows that the behavior of the parameters ∆ m,w and ∆ m,u is much different than ∆ m,1 and ∆ m,2 . The amplitudes of their oscillation are equal but not one anymore, and become very small for greatly unequal coupling strengths. Thus we conclude that the result of unequal coupling strengths is marked, always reducing the amplitudes of the variation of the separability parameters, and optimal entanglement occurs only between the mechanical mode and one of the cavity modes whose the coupling strength is larger than the other mode.

.(4.10)
Similarly, by substituting for |U j (ψ)| from equation (2.16) we obtain an explicit solution for the phase dependence of the separability parameters (j = w, u), in which the upper sign "+" at cos 2ψ is for j = w and the lower sign "−" is for j = u. It is seen from equations (4.10) and (4.11) that ∆ m,j | h j =1 varies periodically with the phase 2ψ. It is not difficult to verified that the periodicity is the same as the periodicity of the oscillation of the populations of the mode j. For a sufficiently large squeezing parameter r, such that e −2r ≈ 0, choosing sin 2φ = 1, one finds from equation (4.10) that the choice of the relative phase sin 2ψ = ±1 leads to a reduction of ∆ m,1 | h 1 =1 and ∆ m,2 | h 2 =1 to zero. Similarly, for a particular choice of ψ at which cos 2ψ = ±1, the term inside the curl brackets of equation (4.11) is equal to zero. Thus, ∆ m,u | hu=1 and ∆ m,w | hw=1 can be reduced to zero. It follows that perfect entanglement is achievable. Hence, with the presence of the linear coupling J, optimal entanglement can be created between the mechanical and one of the field modes. Note that the optimal entanglement in this case is same with the first case given in equation (4.9), where the optimal gain factor h j = 1 at same specific choices of phase for mode j. Moreover, we have also investigated the effects of phase fluctuations on the optimal entanglement by involving the phase noises of the two driving lasers ϕ 1 (t) and ϕ 2 (t) in the Hamiltonian, the details can be found in Appendix A. There we find that the optimal entanglement is reduced for a finite correlation time of the phase noise process. However, we can still achieve entanglement when the phase noise is in a reasonable size.
We have so far ignored the effect of the damping γ of the mechanical mode as this has been secondary to our goal of demonstrating how entanglement varies with the phase. However, to see how the damping γ, if not ignored could affect the entanglement between the modes, we plot in figure 5 the separability parameter ∆ m,j as a function of r for several different values of γ/G. It is seen that as soon as γ/G = 0, there is no optimal entanglement, but the modes can be significantly entangled over entire range of r when γ(n + 1) < G. Even for γ(n + 1) > G the mode still can be entangled but the entanglement is restricted to small r. It is seen that for large γ, the parameter ∆ m,j passes through a minimum, which occurs at small values of r, and then tends to a nonzero value γ(n + 1)/G as r increases. Notice that the value at which ∆ m,j saturates is smaller than unity when γ(n + 1) < G.  The behavior of ∆ m,j can be understood as follows. For ∆ = 0, |g 1 | = |g 2 |, and the particular values of ψ at which |U j (ψ)| = 1, the separability parameter ∆ m,j can be expresses as the sum of two terms (n ≡ n 0 = n m ) ∆ m,j = (n + 1) e 2r 1 + (4.12) Notice that the first term on the right-hand side (RHS) of equation (4.12) exhibits the same variation with r as in the absence of the damping. This term tends to zero as r increases. Therefore, the second term on RHS represents the major contribution of the damping. This term tends to reduce the entanglement but only for large r. It is easily found that in the case of small r, the second term can be approximated by − 1 − 2 √ 2r γ/G, from which we see that there is a threshold value of r = 1/8 below which this term is negative. This has the effect that below threshold (r < 1/8) the parameter ∆ m,j is not much affected by the damping γ, the negative term compensates a part of the positive contribution of γ to the amplitude of the first term. Above threshold r > 1/8, the second term is positive resulting in an increased value of ∆ m,j . Nevertheless, even in this case ∆ m,j can still be smaller than unity. It is easy to see, for large values of r, such that e −2r ≈ 0, the first term on RHS of equation (4.12) vanishes, but the second does not and takes a simple asymptotic value ∆ m,j = (n + 1) γ/G. Thus, in the presence of the damping γ, the separability parameters saturate at a nonzero value, which is smaller than unity whenever γ(n + 1) < G. Hence we may conclude that ∆ m,j is not affected much by the damping γ when r < 1/8. The damping γ tends to affect ∆ m,j at large values of r and causes the separability parameter to saturate at the value γ(n + 1)/G.

Phase control of quantum steering
In Section 4 we discussed the control and transfer of entanglement between the mechanical and field modes by manipulating the relative phase of their coupling strengths. In this section, we discuss how such phase sensitive coupling may be used to control quantum steering of the modes.
Quantum steering takes place between entangled modes and provides the information as to how a given mode "steers" the other mode to be entangled [6,7]. The quantum steering of mode i by mode j can be identified by E i|j = ∆ inf,j X i ∆ inf,j P i < 1/2 ( = 1), where the inferred quadrature variances are defined as ∆ are arbitrary observables of the system j, and u j , u j are gain factors. The quadrature components and the gain factor are selected such that they minimize the inference (conditional standard deviation) product. In our case, ∆ 2 X j = ∆ 2 P j as shown in equation (4.7), and we choose O j = P j and O j = X j , the optimal gain factor is obtained as u j = u j = − X i , P j /∆ 2 P j via variation method ∂E i|j /∂u j = 0. As a result, we have We see from equation (5.1) that a nonzero correlation X i , P j is required to achieve the steering condition. In Section 4, we found that the nonzero correlations exist only between the mechanical mode and the output cavity modes as well as between the optomechanical mode and the output superposition modes. Therefore, we will look at the steering of the mechanical mode by either one of the output cavity modes or one of the output superposition modes. Substituting the variances and correlations given in equation (4.7) into equation (5.1), we can calculate the steering parameters analytically. In the limit of γ G, we get a simple analytical expression of the steering parameter where |U j (ψ)| 2 is given in equation (4.8). We see that the mechanical mode can be steered in a similar manner by the two field modes (j = 1, 2) and their superposition modes (j = u, w). The phase dependent function |U j (ψ)| 2 plays the major role in the steering properties of the modes. If the relative phase ψ is such that |U j (ψ)| 2 = 1, then the steering parameter E m|j falls below its vacuum level 1/2 by a minimum value which dependent directly on n 0 and the degree of squeezing r. If n 0 = 0 and if |U j (ψ)| 2 = 1, then quantum steering (E m|j < 1/2) is seen to occur over the entire range of r. When n 0 = 0, quantum steering is possible for r > r th , where r th = ln[(2n 0 + 1)/(n 0 + 1)]/2. The minimum value of E m|j is reached periodically when |U j (ψ)| 2 = 1 and r 1, in which case E m|j = 0 and the steering is perfect. Thus, we may speak about perfect EPR state between the mechanical mode m and mode j. The above considerations are illustrated in figure 6, which shows the variation of the steering parameters E m|1 and E m|2 with the phase difference ψ for |g 1 | = |g 2 | and |g 1 | = |g 2 |. It is seen that for r 1 and the phases ψ at which |U j (ψ)| 2 = 1, the steering parameter E m|j approaches zero that optimum steering is achieved. For the asymmetric case of |g 1 | = |g 2 | the optimum steering occurs only between the mods A out m and A out 2 . Thus, we can achieve perfect steering of the mechanical mode by one of the field modes simply through varying the phase difference ψ. We can also see that for a given phase ψ only one of E m|1 and E m|2 can be smaller than 1/2. In other words, when mode 1 can steer the mechanical mode, mode 2 cannot. There is a simple physical interpretation of this feature by referring to the monogamy condition that two or more modes cannot simultaneously steer another mode [82]. It is interesting to note that no such rigorous constrain exists for entanglement. As one can see from figures 3 and 4, two field modes can be simultaneously entangled with the mechanical mode.
Another interesting feature of the quantum steering present in this system concerns the steering of the mechanical mode by the superposition modes appearing as collective modes. In accordance with the monogamy condition, the presence of the collective steering of the mechanical mode, either E m|u < 1/2 or E m|w < 1/2 is accompanied by the absence of the bipartite steering, both E m|1 > 1/2 and E m|2 > 1/2. In other words, collective steering is present in regimes that do not show any bipartite steering. Therefore, the phase dependence of the steering parameters makes the collective steering clearly distinguishable from the bipartite steering. In effect, we can observe perfect one-sided devices independent quantum secret sharing, which is achieved by creating collective steering without creating the bipartite steering [27].
Finally, we consider the effect of the damping rate γ and the thermal noise n at the mechanical mode on the steering properties. Figure 7 shows the effects of the rate γ and the thermal photons n on the optimum value of the steering parameter by taking |U j (ψ)| 2 = 1. The damping γ and the thermal noise n affect the optimum steering in a similar way. Interestingly, however, a weak sensitivity to γ and n occurs in a region of the squeezing parameter r significantly different from that which maximizes steering in the absence of the damping. Namely, for short pulses (r < 1), the effect of γ and n is seen to be negligibly small, that E m,j remains quite close to zero over a large range of γ and n. For r 1, the minimum of E m|j degrades more rapidly with γ and n. In concluding this section, we comment about how to test for the presence of bipartite steering of the mechanical mode generated by the cavity modes 1 and 2, and the collective steering generated by the modes u and w. The generation of the bipartite and collective steerings is ultimately connected to the redistribution of the population among the cavity modes. The bipartite steering E m|1 < 1/2 or E m|2 < 1/2 is generated by the population transfer from mode 2 to 1 or mode 1 to 2 induced by phase difference, respectively. The collective steering, on the other hand, is generated by the equal distribution of the population between the modes 1 and 2. One can confirm bipartite and tripartite steerings simply by measuring the coincidence rate R 12 between two detectors D 1 and D 2 adjusted to collect photons of the output modes A out 1 and A out 2 , respectively. According to equation (3.18), the coincidence rate exhibits a cosine modulation with the phase ψ, with the depth of modulation varying between 0 when one of the output modes, either A out 1 or A out 2 is not populated, and 2[κg 2 Υ(r)/G(κ 2 + J 2 )] 2 when they are equally populated. Thus, the minima of the interference pattern of R 12 signal the bipartite steering, either E m|1 < 1/2 or E m|2 < 1/2, and the maxima signal the collective steering, either E m|u < 1/2 or E m|w < 1/2.

Conclusions
We have developed a theory of the phase control of mutual coherence, entanglement, and quantum steering in a system composed of a single mode cavity containing a dielectric membrane in its interior. The membrane was treated as a partly transmitting and reflecting mirror which divides the cavity into two mutually coupled optomechanical cavities. The resulting system has been found to behave effectively as a three-mode linked coupling or closed loop system, the two cavity (optical) modes and the mechanical mode representing the oscillating membrane. The effect of this closed loop is to create interfering channels which depend on the relative phase of the coupling strengths of the field modes to the mechanical mode. The influence of this phase on the population distribution and correlations between the modes has been calculated for the case of short laser pulse excitations of the cavity modes. Manipulation of the phases, which can be done by varying the laser phases, can give rise to a transfer of the population between the cavity modes and to create superposition (collective) modes.
Our calculations demonstrate that the phase controlled transfer of the population may result in an interesting phenomenon of the induced coherence without induced emission. This phenomenon allows the cavity modes to remain perfectly coherent even if one of the modes is not populated. It has been found that the perfect coherence among the field modes excludes the possibility of the modes to be entangled. On the other hand, it makes possible of the modes to be maximally entangled with the mechanical mode and to behave deterministically that the entanglement can be periodically transferred among the modes by varying the relative phase. Furthermore, we have investigated the effect of the relative phase on quantum steering of the mechanical mode by the field modes. It has been found that for a given phase difference only one of the four steering parameters can be smaller than the vacuum level limit. In other words, the mechanical mode can be steered by one of the field modes only. Moreover, we have demonstrated that the phase dependence makes the collective steering clearly distinguishable from the bipartite steering. Such effects can be observed by measuring the coincidence rate between two detectors adjusted to collect photons of the output cavity modes.
The phase dependent dynamics of this close loop system suggest a measurement technique to distinguish between the bipartite and collective steerings. These different types of quantum steering can be observed by measuring the coincidence rate between two detectors adjusted to collect photons of the output cavity modes. In particular, we have found that the minima of the interference pattern of the coincidence rate signal the bipartite steering, while the maxima signal the collective steering. This makes the system suitable to observe perfect one-sided devices independent quantum secret sharing. Finally, we have considered the effect of the damping rate of the mechanical mode on the optimal entanglement and quantum steering and have found that for short pulses the effect of the damping is negligibly small, that the inseparability parameter and the degree of quantum steering remain almost unchanged over a large range of the damping rate. In addition, we have found that the effect of the thermal noise at the mechanical mode on the entanglement and quantum steering is similar to that of the damping of the mode.
Here we have used the same rotating farm for equation (2.5), δc m = δce iωmt , δc in,m = δc in e iωmt , δa r,m j = δa r j e −i(ωmt+ϕ 0j ) , δa in,r,m j = δa in,r j e −i(ωmt+ϕ 0j ) . And for clarity of the notation, we have omitted the symbol δ and the superscripts r and m on the displacement operators, as done for equation (2.5).
In the bad cavity limit κ |g w,u |, we setȧ w =ȧ u = 0, and find Incorporating the input-output relation a out = a in + √ 2κa, we can get the following expressions for the populations and the correlation functions of the normalized temporal field modes (A out m ) † A out m = n 0 + Γ(r) + |β wu | 2 e 2r D † (t)D(t ) , in which β wu = g * w α w e −iφ + g * u α u e iφ √ κ 2 + w 2 , D = We see that the contribution of the terms resulting from the phase fluctuations are proportional to |β wu | and to the noise correlation functions D † (t)D(t ) etc. Thus, the contribution of these terms could be negligibly small if either |β wu | 1 or D † (t)D(t ) 1. Since |β wu | ∼ |α i |, (i = 1, 2) and |α i | 1, the condition of |β wu | 1 is not met in our optomechanical system. Therefore, to determine the correlation functions D † (t)D(t ) etc., appearing in equation (A.6), we consider a colour-noise model of the excitation laser for which the correlation function σ(t)σ(t ) is given by [84] σ(t)σ(t ) = Γ l γ c e −γc|t−t | , (A.8) where Γ l is the laser bandwidth, and the parameter γ −1 c characterizes a finite correlation time of the phase noise process. The limit γ c → ∞ corresponds to the case of a whitenoise laser.
Note that in the white-noise limit of γ c → ∞, the phase fluctuation terms in equation (A.6) become proportional to |β wu |Γ l , which for a finite linewidth Γ l could be negligibly small only if |β wu | 1. As we have already mentioned the condition |β wu | 1 is not met in the system, so one can expect a significant contribution of these terms to the populations and the correlation functions of the normalized temporal field modes in the white-noise limit.
For a colour-noise laser γ c could be very small resulting in a negligibly small contribution of the phase fluctuations terms. In order to show this, we evaluate the separability parameters, as given in equation (4.9), but with the correlation functions (A.6) which include the effect of the fluctuations of the excitation lasers. In particular, we evaluate the effect of the phase fluctuations on the optimal entanglement previously obtained for the monochromatic case and illustrated in figure 3(b), i.e., ∆ m,w = 0 for the phase ψ satisfying ψ = nπ, (n = 0, 1, 2, ...) where the optimal gain h w = 1, and ∆ m,u = 0 for the phase ψ satisfying ψ = (n + 1/2)π, (n = 0, 1, 2, ...) where the optimal gain h u = 1, respectively. Figure A1 shows the variation of the optimal entanglement with γ c as determined by the separability parameters ∆ m,w and ∆ m,u for different relative phases, frame (a) ψ = 0, and frame (b) ψ = π/2. One can see that the optimal entanglement between field mode A out w and the mechanical mode A out m rapidly disappears at γ c ∼ 10 −5 Hz, while  Figure A1. The effect of the laser fluctuations on the minimum values of the inseparability parameters ∆ m,w and ∆ m,u obtained for the case of the monochromatic excitation. The parameters are g 0,1 = g 0,2 = 10 Hz, |g 1 | = |g 2 | ≡ g = 0.1 MHz, J = κ = 10g, ∆ = 0, γ = 0, n 0 = n m = 200, r = Gτ = 5, Γ l = 1 kHz. The black solid line is for ∆ m,w , and the blue dashed line is for ∆ m,u . In frame (a), ψ = 0 at which ∆ m,w approaches 0 in the case of the excitation with a monochromatic laser. In frame (b), ψ = π/2 at which ∆ m,u → 0 in the monochromatic case.
the optimal entanglement between the modes A out u and A out m disappears at γ c ∼ 10 2 Hz. The reason why the separability parameters ∆ m,w < 1 and ∆ m,u < 1 have different thresholds for γ c is that α u α w determined by the considered parameters |g 1 | = |g 2 |, ∆ = 0, κ = J. We may conclude that similar to the laser phase fluctuation effects on the optomechanical cavity cooling [83,84] and on the entanglement [85][86][87], the optimal entanglement generated in this scheme, corresponding the perfect EPR state, can be experimentally observed if a colour rather than white-noise laser is used to excite the cavity modes.