05 52 2 v 1 [ qu an tph ] 1 5 M ay 2 01 8 Output Entanglement Saturation and Optimization in Cavity Optomechanics

Xiao-Bo Yan, ∗ Zhi-Jiao Deng, Xue-Dong Tian, and Jin-Hui Wu † College of Electronic Science, Northeast Petroleum University, Daqing 163318, China Department of Physics, National University of Defense Technology, Changsha 410073, China College of Physics Science and Technology, Guangxi Normal University, Guilin 541004, China School of Physics, Northeast Normal University, Changchun 130024, China (Dated: May 16, 2018)

PACS numbers: 42.50.Ex, 42.50.Wk, 07. 10.Cm Introduction.-Entanglement is the distinguishing feature of quantum mechanics because it is responsible for nonlocal correlations between observables and now it has become a basic resource for many quantum information processing schemes [1]. For example, entanglement is required in quantum teleportation in which quantum information can be transmitted from one location to another with security. Extending entanglement into macroscopic systems has become a prominent experimental objective, and would allow us to explore the quantum-classical boundary. So far, a number of theoretical and experimental works on entanglement between macroscopic objects have been studied, such as between atomic ensembles [2,3], between superconducting qubits [4][5][6][7], and between mechanical oscillator and microwave fields [8].
Recently, quantum entanglement in cavity optomechanics has received increasing attention for the potential to use radiation pressure to generate various entanglement between subsystems . Due to the small mechanical decay rate in optomechanical systems, the information can be stored in mechanical modes for long time. Hence, entangled optomechanical systems could be profitably used for the realization of quantum communication networks, in which the mechanical modes play the role of local nodes where quantum information can be stored and retrieved, and optical modes carry this information between the nodes [9][10][11][12].
In fact, any quantum communication application involves traveling output modes rather than intracavity ones. Therefore, it is very important to study how to obtain optimal output entanglement in cavity optomechanics. Here, we use a three-mode optomechanical system (see Fig. 1) to generate output entanglement between two filtered output optical fields. This setup has been realized in several recent experiments [39][40][41]. Because in such a system the parametric-amplifier interaction and the beam-splitter interaction can entangle the two intracavity modes, the output cavity ones are also entangled with each other. In previous theory works [34][35][36][37][38], the output entanglement in this model has been studied. But in most of them, the output entanglement were studied with equal-coupling, zero filter bandwidth or without time delay between the two output fields.
In this paper, we mainly focus on how output entanglement saturation emerges, and how to obtain optimal output entanglement with optimal coupling and optimal time delay for large bandwidth. First, we obtain the important expression of optimal time delay between the two filtered output fields. With large bandwidth and no time delay, the optimal output entanglement appears at the point where the optimal time delay equals zero, from which we obtain the expression of optimal coupling. It is very interesting that the optimal output entanglement will saturate if the driving field is strong enough, and the saturation value is sensitive to filter bandwidth. As far as we know, the phenomenon of entanglement saturation has not been studied before. Finally, we obtain the optimal output entanglement with optimal time delay and optimal couplings. We believe the results of this paper are very important to experimental and theoretical physicists who work on entanglement or quantum information processing.
System model.-We consider a three-mode optomechanical system in which two cavities are coupled to a common mechanical resonator (see Fig. 1). The Hamiltonian of the system reads Here,â i is the annihilation operator for cavity i with frequency ω i and damping rate κ i ,b is the annihilation operator for mechanics resonator with frequency ω m and damping rate γ, and g i is the optomechanical coupling strength. In order to generate the steady entanglement between the two output fields, we drive cavity 1 (2) at the red (blue) sideband with respect to mechanical resonator: ω d1 = ω 1 −ω m and ω d2 = ω 2 +ω m . If we work in a rotating frame with respect to the free Hamiltonian, following the standard linearization procedure, and making the rotating-wave approximation (in this paper, we focus on the resolved-sideband regime ω m ≫ κ 1 , κ 2 ), hence, the Hamiltonian of the system can be written aŝ Here,d i =â i −ā i ,ā i being the classical cavity amplitude. G i is the effective coupling strength which can be easily controlled by adjusting the strength of driving fields. Here, without loss of generality, we take G i as real number. Based on Eq. (2), the dynamics of the system is described by the following quantum Langevin equations for relevant operators of mechanical and optical modes Here,b in ,d in i are the input noise operators of mechanical resonator and cavity i(i = 1, 2), whose correlation functions are bin, respectively. N m and N i are the average thermal populations of mechanical mode and cavity i, respectively. In the following discussion, we mainly study how output entanglement saturation emerges, so we assume these average thermal populations are zero (zero temperature). According to the Routh-Hurwitz stability conditions [42] and we focus on the regime of strong cooperativities C i ≡ 4G 2 i /(γκ i ) ≫ 1 and κ i ≫ γ in this paper, the stability condition of our system can be obtained as G 2 1 /G 2 2 > max(κ 1 /κ 2 , κ 2 /κ 1 ) for κ 1 = κ 2 , and the system is always stable if κ 1 = κ 2 and G 2 ≤ G 1 [24,37].
Output fields and optimal time delay.-In a quantum network, entangled output photon pairs are a useful resource for quantum information processing. Here, the output entanglement can be generated when the output optical fields pass through the filter functions f (ω) which satisfy +∞ −∞ |f (ω)| 2 dω = 1. For simplicity, we adopt a rectangle filter with a bandwidth σ centered about the frequency ω to generate the output temporal modes, i.e., the Heaviside step function. Then, the filtered optical output fields can be written aŝ Here, ω ± = ω ± σ 2 , and τ i is the absolute time at which the wave packet of interest is emitted from cavity i. The time delay between this two output fields is defined as τ = τ 1 − τ 2 . Without loss of generality, we set τ 1 = τ and τ 2 = 0. Using the Langevin equation Eq. (3) and the input-output relation [43], we can derive the output operators in terms of the input operators. And we use the logarithmic negativity [44,45] to quantify the entanglement between the filtered output cavity modeŝ D out 1 [ω, σ, τ ] andD out 2 [−ω, σ, 0]. We writeD out i [ω, σ, τ i ] asD i and set equal cavity damping rate κ 1 = κ 2 = κ for simplicity, and the mechanical damping rate γ is much less than other parameters.
It can be proven that our system is equivalent with a two-mode squeezed thermal state [38]. From the equivalent relation, it can be concluded that the optimal time delay for output entanglement between the two filtered output fields is the time delay that makes the modulus of the correlator D 1D2 reach a maximum [38]. Through complex calculation, we obtain the expression of the optimal time delay for resonant frequency (ω = 0 in the rotating frame) as follows For the special case of equal coupling (G 1 = G 2 ) and σ ≪ κ, the expression of τ opt will become κ/4G 2 1 which is consistent with the result in Ref. [36,37]. It can be seen from Eq. (5) the optimal time delay τ opt will be positive with the increase of G 2 . We will see that the coupling G 2 at which the optimal time delay equals zero is a very special coupling. In the following, we study the output entanglement in two cases, without and with time delay respectively.
Optimal couplings and output entanglement saturation.-In this part, we study how to obtain the optimal coupling and the output entanglement saturation in the case of without time delay (τ = 0). The entangling interaction G 2 (bd 2 +b †d † 2 ) inĤ int entangles the mechanical resonatorb and cavity moded 2 , and the beam splitter interaction G 1 (b †d 1 +bd † 1 ) swaps theb andd 1 states, these two combined interactions yield the net entanglement betweend 1 andd 2 . The entanglement generated within the cavities can be transferred to the filtered output fields. It is obvious that the output entanglement will disappear as coupling G 2 = 0. Moreover, the output entanglement for resonant frequency almost equals zero in the case of equal-coupling G 2 = G 1 [38]. Therefore, it is certain that there will be an optimal coupling G opt 2 that makes the output entanglement E N (ω = 0) reach a maximum.
Before getting the expression of optimal coupling G opt 2 , we do some analysis first. It is certain that the output entanglement with optimal time delay will be larger than that without time delay. In addition, it can be seen from Eq. (5) that the optimal time delay may equal zero. It means that the curve of output entanglement with optimal time delay will be tangent to that without time delay at the point where the optimal time delay equals zero. Hence, we guess the output entanglement without time delay will take the maximum at the tangent point. In Fig. 2(a), the output entanglement E N (ω = 0) is plotted vs G 2 /G 1 with large filter bandwidth σ = κ (red line), σ = κ/10 (black line), and the dashed lines are their corresponding lines with numerical optimal time delay. The other parameters are γ = 1, κ = 10 5 , and G 1 = 10κ. It can be seen from Fig. 2(a) that the output entanglement actually takes the maximum at the tangent point for large bandwidth. Hence, according to Eq. (5), the optimal coupling for large bandwidth can be obtained as In Fig. 2(b), we plot the output entanglement E N (ω = 0) vs G 2 /G 1 for small bandwidth σ = κ/10 4 (red line), σ = 2κ/10 3 (yellow line) and their corresponding common line with numerical optimal time delay (blue dashed line). It can be seen from Fig. 2(b) that the maximum of output entanglement will not appear at the tangent point. This is because time delay has no significantly effect on the output entanglement in the case of small bandwidth. For small bandwidth, the optimal coupling can be obtained as From Eq. (6) and Eq. (7), we obtain the boundary between small bandwidth and large bandwidth as . We also plot the output entanglement E N (ω = 0) with boundary bandwidth σ = σ b (violet line) in Fig.  2(b) from which it can be seen the maximum value of output entanglement will not appear at the tangent point anymore as σ < σ b .
Based on Eq. (6) and Eq. (7), we can obtain the expressions about maximum value of output entanglement. But most of the expressions are lengthy. Here, we just give an interesting result for large bandwidth, i.e. output entanglement saturation. In general, entanglement will increase with the increase of driving strength. But we find the optimal output entanglement E N (ω = 0) (with optimal coupling Eq. (6)) will saturate as G 1 ≫ κ 5 √ 3σ 3 . In Fig. 2(c), we plot optimal output entanglement E N (ω = 0) vs G 1 /κ for bandwidth σ = κ (red line), σ = κ/2 (blue line), σ = κ/10 (black line). It can be seen clearly from Fig. 2(c) that the optimal output entanglement will approach a constant (saturation value) with the increase of coupling G 1 . The saturation value can be obtained as with α = 5κ 2 + 3σ 2 , β = κ arctan σ κ . The saturation value plotted according to Eq. (8) in Fig. 2(c) (three green dashed lines) perfectly fits the numerical result. Due to entanglement saturation, the coupling G 1 is not the stronger the better. It can also be seen from Fig.  2(c) that for larger bandwidth, such as σ = κ (red line), the coupling G 1 slightly greater than κ can make output entanglement saturation occur. For the case of σ ≪ κ, the saturation value Eq. (8) can be simplified as ln[ 175κ 6 4σ 6 ] which is sensitive to filter bandwidth for high power.
Optimal output entanglement.-Now, we study how to find the optimal output entanglement with optimal time delay (τ = τ opt ). First, we plot the optimal time delay vs. G 2 /G 1 according to Eq. (5) (blue solid line) and the numerical result (red dashed line) in Fig. 3(a). The optimal time delay will become positive with the increase of G 2 (see the inset of Fig. 3(a)). In Fig. 3(b), we plot the corresponding output entanglement with these two optimal time delay with parameters γ = 1, κ = 10 5 , σ = κ, G 1 = 10κ. The two curves of output entanglement fit very well except a very small area (see the inset of Fig.  3(b)). It can be seen from Fig. 3(b) that the maximum of output entanglement does not appear at the point where the optimal time delay equals zero. If G 1 ≫ κ, σ, we can obtain the optimal coupling and the optimal output entanglement as and ]. (10) We plot the highest point according to Eq. (9) and Eq. (10) in Fig. 3(b) (see the green dot), and it fits the numerical result very well. Until now, we study the output entanglement just for resonant frequency (ω = 0). In Fig. 3(c), we plot output entanglement E N (ω) vs. ω/κ with optimal coupling Eq. (9) and with numerical optimal time delay (red line). It is obviously seen from Fig. 3(c) that the optimal output entanglement Eq. (10) is just the optimal one in the whole center frequency domain of output fields. For contrast, we also plot the output entanglement with equal-coupling (G 2 = G 1 ) and without time delay (τ = 0) (blue dashed line). It can be seen from Fig. 3(c) the output entanglement in the vicinity of resonant frequency can be largely enhanced by adopting optimal coupling Eq. (9) and optimal time delay. In experiment, it just need to adjust the coupling strength G 2 according to the parameters of the systems to obtain the optimal output entanglement.
Conclusions.-In summary, we have studied theoretically the output entanglement between two filtered output fields in a three-mode cavity optomechanical system. We obtain the important expression of optimal time delay Eq. (5) between the two filtered output fields. Based on the analysis about this expression, we obtain the optimal coupling Eq. (6) for output entanglement in the case of without time delay. We give a reasonable boundary between large bandwidth and small bandwidth. For large bandwidth, the optimal output entanglement will saturate if the coupling strength satisfies G 1 ≫ κ 5 √ 3σ 3 , and the saturation value is obtained as Eq. (8). In addition, using the optimal coupling Eq. (9) and optimal time delay, we find the optimal output entanglement Eq. (10) for resonant frequency (ω = 0) is just the optimal one in the whole center frequency domain of output fields. Our results can also be applied to other parametrically coupled three-mode bosonic systems, and may be useful to experimentalists to obtain large entanglement.