Dynamics of two levitated nanospheres nonlinearly coupling with non-Markovian environment
Li Xun, Xiong Biao, Chao Shilei, Jin Jiasen, Zhou Ling
School of Physics, Dalian University of Technology, Dalian 116026, China

 

† Corresponding author. E-mail: jsjin@dlut.edu.cn zhlhxn@dlut.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11874099, 11605022, 11775040, 11747317, and 11474044).

Abstract
Abstract

The dynamics of two nanospheres nonlinearly coupling with non-Markovian reservoir is investigated. A master equation of the two nanospheres is derived by employing quantum state diffusion method. It is shown that the nonlinear coupling can improve the non-Markovianity. Due to the sharing of the common non-Markovian environment, the state transfer between the two nanospheres can be realized. The entanglement and the squeezing of the individual mode, as well as the jointed two-mode are analyzed. The present system can be realized by trapping two nanospheres in a wideband cavity, which might provide a method to study adjustable non-Markovian dynamics of mechanical motion.

1. Introduction

Trapping particles via an electromagnetic field is widely researched theoretically and experimentally because of its potential application in precise instrument and macroscopic quantum effects detection. For example, an optical tweezer is applied to capture and control cells in biological and medical research. The realization of Bose–Einstein condensation (BEC) and atom cooling have strict requirements for optical traps.[14] With the development of quantum techniques, the optical trap has been exploited in several research domains, such as generating squeezing states[5] and macroscopic superposition states of trapped matters,[68] entangling the trapped particles,[9,10] and cooling the trapped matter to ultra-low temperature.[1114] In addition, the nonlinearity of the levitated system has been reported in Refs.[1217]. The dynamics of a nanoparticle that nonlinearly couples to its reservoir has been theoretically studied in the Markovian regime.[18,19]

Non-Markovian theory plays an important role in treating a realistic system. The non-Markovian environment usually can keep the coherence and suppress the dissipation. The theory of open quantum system is researched through several perspectives, such as the projection operator theory,[20,21] Green's functions method,[2224] perturbation theory,[25] and stochastic Schrödinger SSEs techniques.[2630] However, these researches mostly investigate the system linearly coupling with its non-Markovian environment. However, the system coupled with its environment in nonlinear form remains under-studied.

In this paper we consider two nanospheres trapped by a continuous frequency field and we write the Hamiltonian of the system with a nonlinear coupling to its environment. By employing quantum state diffusion (QSD) method,[2830] we derive a nonlinear master equation. We investigate the dynamics of the two levitated nanospheres with open quantum system techniques, and study the entanglement, squeezing of the two nanospheres, as well as the state transfer. We find that the nonlinear coupling can improve non-Markovianity.

2. Model and Hamiltonian

We consider two nanospheres trapped by a wideband cavity. The wideband cavity field can be considered as a non-Markovian environment of the two levitated nanospheres. As found in Refs. [5] and [17], the vibration frequency of the trapped nanospheres is on the order of kHz or MHz, which is smaller than the width of the trapping field. This means that we can treat the trapping fields as a wideband cavity. Recently, the crystal cavity with engineering modes, whose frequency is continuous, has been employed for trapping nanoparticles.[31] When acted upon by the optical forces, named the gradient force Fgrad and the scattering force Fscatt, the two nanospheres can be described as harmonic oscillators that linearly and nonlinearly couple to the wideband cavity field (the detail illustration can be seen in Ref. [6]). The Hamiltonian can be written as

where ak ( ) is the annihilation (creation) operator of the k-th mode with detuning , is the levitated frequency, qi is position (momentum) operator of the i-th nanosphere, and Gk is the coupling between the k-th cavity field and nanospheres.

The Hamiltonian (1) describes the two nanospheres wrapped up by photon gas with a continuous spectrum distribution, which is analogous to Brownian motion. The system should be solved with non-Markovian theory. However, the third summation in Eq. (1) contains the nonlinear term in order of qi which differs from the usual linear non-Markovian interaction. The system can not be solved with the general non-Markovian theory as done in Refs. [23] and [32].

To treat the nonlinear and non-Markovian dynamics, we employ QSD method. For convenience, we rewrite the Hamiltonian into three parts, the free energy of the environment , the energy of the two nanospheres , and the interaction part , where . It is convenient to change into the interaction picture rotating with HB, then we have the time-dependent Hamiltonian

with
The open quantum system obeys the non-Markovian quantum state diffusion equation
We assume the initial state of the system and environment as . The system state is arbitrary, and the initial states of the environment are their ground states . The total state can be expanded as
where , and with identity for any Bargmann coherent state
The stochastic state of system obeys QSD equation
where
is the correlation function of the bath. The variation in Eq. (7) can be taken place by some time-dependent operator O satisfying . Therefore, the non-Markovian QSD equation (7) is rewritten in a pseudo-time-local form
where . The operators and O satisfy
The reduced density operator of the open system can be obtained after partial trace operating over total state , that is
We denote the reduced density operator with operator ensemble mean over all stochastic states. The master equation is obtained as
However, it is still difficult to find the exact solution of O and with Eq. (10), especially when the system contains nonlinear or multi-particle interaction. An alternative method is to expand it in terms of the weak coupling strength, i.e., , where G represents the coupling strength. The lowest order of this operator can be written as
by employing post-Markov approximation, which is valid when the square of the ratio of the environment memory time scale to system time scale is negligible.[33,34] The integrated operator should be expressed as
where the coefficients can be found in Appendix B. From Eq. (12), the master equation with nonlinear terms is obtained as
where is the label of nonlinear coupling. The contributions of the linear and nonlinear couplings can be discussed by choosing equal to zero or not. In addition, the method of the master equation deduction (15) can be easily extended to that with nonlinear coupling. If the two nanospheres interact with other fields except the non-Markovian reservoir, then one only need to add their interaction into the Hamiltonian and perform the above procedure again.

3. Squeezing of levitated nanospheres induced by the wideband cavity field

From the first line of the summation in Eq. (15), we can obtain an effective interaction between the two nanospheres resulting from the linear coupling, i.e., , which can induce a self-squeezing and two-mode squeezing of the nanospheres.

Meanwhile, the nonlinear coupling between photon and phonon with induces another effective interaction

which describes a Kerr effect.

We would like to investigate the effect of nonlinear coupling under Lorentz spectrum , which is a typical example to analyze an open system thanks to its realism in experiment, where is the spectrum strength and d0 is the spectrum width.

Squeezing can be quantified by the fluctuation of operators. The position fluctuation is defined as . The two-mode fluctuation can be defined as , where . We choose and assume that the initial states of the nanospheres are coherent states with and . As shown in Fig. 1(a), the fluctuation with linear coupling attains a relative “stable” value (with weak amplitude oscillation). When , the single-mode squeezing finally disappears after long time evolution. For two-mode squeezing, there is little difference for or not. Since the two-mode squeezing exhibits the correlation between the two nanospheres, the nonlinear coupling slightly affects the correlation. In addition, we can observe that the larger , the stronger the two-mode squeezing, because is the spectrum strength. The spectrum width d0 also has a little bit influence on the two-mode squeezing. One can see that the squeezing with is stronger than that with as shown in Fig. 1(b), but it does not mean the bigger d0 the stronger squeezing. In our simulation, we find that for , there is no squeezing. If , , then it is Markovian and the squeezing is obviously worse than that of non-Markovian.

Fig. 1. The evolution of (a) single-mode fluctuation and (b) two-mode fluctuation of the levitated nanospheres. In panel (a), , , and . In panel (b), for the green line, and for other dashed lines, . We let .
4. Entanglement and state transfer

The two-mode squeezing represents the correlation between two nanospheres. Quantum entanglement, as an alternative correlation between two subsystems, can also be induced by effective interaction. The entanglement can be measured by negativity, which is defined as

where is the minimum symplectic eigenvalue of 2 × 2 block covariance matrix
with element ( ) and .

We plot the negativity in Fig. 2. The entanglement does not appear when then increases fast(sudden birth) and stabilizes without nonlinear coupling. When the nonlinear coupling is included, the entanglement is always less than that with , finally drops to zero, which means that the nonlinear coupling is not good for preserving the existence entanglement for a long time. This happens because the nonlinear coupling also induces an addition damping. We will explain this point later.

Fig. 2. The evolution of negativity where and . Other parameters are the same as those in fig. 1.

The effective interaction connects the two nanospheres, thus if the nanospheres are in different initial states , then the states ( ) will transfer to ( ) through effective interaction. We employ the conversion fidelity to measure the efficiency of state transfer.[35] We plot the fidelity in Fig. 3. We can observe that the fidelity oscillates periodically with decreasing amplitude, whether or not the nonlinear coupling exists. When , the fidelity is obviously worse than that with . The maximum fidelity can reach 0.9, which is relative ideal, while the entanglement is weak (see Fig. 2) and the squeezing is not so deep (see Fig. 1). This happens because the equal frequency of the two nanospheres makes the effective interaction induced by linear coupling Heff with beam split (BS) form, therefore the state transfer can be achieved while the squeezing and the entanglement are not so ideal.

Fig. 3. The evolution of fidelity. Other parameters are the same as those in fig. 1.
5. Non-Markovianity measure

We have shown that squeezing, entanglement, and state transfer between the two nanospheres can be achieved in our system. These properties result from their common coupling with the same non-Markovian environment; i.e., relating with the non-Markovianity in the present system. We now investigate the non-Markovianity of the system. For the non-Markovianity measure, several proposals[3639] based on violation of dynamical semigroups or information backflow have been put forward. Here, we employ the measure of non-Markovianity proposed in Ref. [36]

where is the derivative of trace distance between different initial states and
with . Under this measure, implies that the dynamical map is non-Markovian. We calculate the non-Markovianity measure by solving Eq. (15). As shown in Fig. 4(a), with the increase of , the non-Markovianity is improved whenever or . Unfortunately, as we found in the previous section, the nonlinear coupling damages the squeezing, entanglement, and state transfer to some degree. This result conflicts with previous works, which declare that the information backflow could preserve those properties rather than damage them.[4043] To uncover the reason for this, we plot the trace distance between different initial states in Fig. 4(b) with a random initial state pair where the initial state is the same for red and blue lines. After a long time evolution, the quasi-steady state may be achieved no matter what the initial state is. Therefore, the decline of for different initial state indicates the dissipation of the system. As shown in Fig. 4(b), when , a larger value of nonlinearity will lead to a faster decline of , which means that a larger value of will lead to a larger dissipation rate. Thus, on the one hand the nonlinearity enhances the non-Markovianity, and on the other hand increases the loss rate. This is the reason why the squeezing and entanglement become bad after a long time evolution.

Fig. 4. (a) The non-Markovianity measure as a function of where . (b) The trace distance between two different initial states for and . Other parameters are the same as those in fig. 11.
6. Discussion and conclusion

We investigate two nanospheres that are trapped in a wideband cavity field, where the two nanospheres nonlinearly couple with the environment. The main contribution of this work is that we put forward a method to treat non-Markovian dynamics with nonlinear coupling. A master equation of the nonlinear non-Markovian system is derived, in which the effective linear and nonlinear interactions can be induced, which can generate the squeezing, entanglement, and realize state transfer. We study the influence of this nonlinear coupling on the dynamics of the system. We find that the nonlinear coupling can improve non-Markovianity. After a long enough evolution, the nonlinear coupling will damage the squeezing and entanglement because the nonlinearity enlarges the loss rate. To compensate for the loss of the system, one can introduce another pumping field. In addition, one can see that the entanglement and the squeezing are not as good as the state transfer because the equal frequency of the two nanospheres results in the main interaction with beam split form. If one would like to obtain large entanglement and squeezing, then one could manipulate the two nanospheres with red and blue detunings, respectively. The spectrum of the environment depends on the distribution of the cavity field, which can be manipulated. Therefore, it might be a way to simulate the spectrum of the environment to study the non-Markovian system.

Acknowledgement

We would like to thank Xinyu Zhao and W. L. Li for helpful discussions.

Reference
[1] Jochim S Bartenstein M Altmeyer A Hendl G Riedl S Chin C Hecker Denschlag J Grimm R 2003 Science 302 2101
[2] Barrett M D Sauer J A Chapman M S 2001 Phys. Rev. Lett. 87 010404
[3] Zippilli S Morigi G 2005 Phys. Rev. Lett. 95 143001
[4] Yi Z Li G X Yang Y P 2013 Phys. Rev. 87 053408
[5] Rashid M Tufarelli T Bateman J Vovrosh J Hempston D Kim M S Ulbricht H 2016 Phys. Rev. Lett. 117 273601
[6] Monroe C Meekhof D M King B E Wineland D J 1996 Science 272 1131
[7] Abdi M Degenfeld-Schonburg P Sameti M Navarrete-Benlloch C Hartmann M J 2016 Phys. Rev. Lett. 116 233604
[8] Zheng S B 2004 Phys. Rev. 69 055801
[9] Welte S Hacker B Daiss S Ritter S Rempe G 2017 Phys. Rev. Lett. 118 210503
[10] Roghani M Helm H Breuer H P 2011 Phys. Rev. Lett. 106 040502
[11] Li T Kheifets S Raizen M G 2011 Nat. Phys. 7 527
[12] Fonseca P Z G Aranas E B Millen J Monteiro T S Barker P F 2016 Phys. Rev. Lett. 117 173602
[13] Millen J Fonseca P Z Mavrogordatos T Monteiro T S Barker P F 2015 Phys. Rev. Lett. 114 243601
[14] Zhang S Wu W Wu C W Li F G Li T Wang X Bao W S 2017 Chin. Phys. 26 074205
[15] Jain V Gieseler J Moritz C Dellago C Quidant R Novotny L 2016 Phys. Rev. Lett. 116 243601
[16] Monteiro T S Millen J Pender G A T Marquardt F Chang D Barker P F 2013 New J. Phys. 15 015001
[17] Gieseler J Novotny L Quidant R 2013 Nat. Phys. 9 806
[18] Voje A Croy A Isacsson A 2015 Phys. Rev. 92 012313
[19] Duarte O S Caldeira A O 2009 Phys. Rev. 80 032110
[20] Nakajima S 1958 Prog. Theor. Phys. 20 948
[21] Zwanzig R 1960 J. Chem. Phys. 33 1338
[22] Hu B L Paz J P Zhang Y 1992 Phys. Rev. 45 2843
[23] Zhang W M Lo P Y Xiong H N Tu M W Y Nori F 2012 Phys. Rev. Lett. 109 170402
[24] Cheng J Liang X T Zhang W Z Duan X 2018 Chin. Phys. 27 120302
[25] Cavina V Mari A Giovannetti V 2017 Phys. Rev. Lett. 119 050601
[26] Piilo J Maniscalco S Härkönen K Suominen K A 2008 Phys. Rev. Lett. 100 180402
[27] Strunz W T Diosi L Gisin N Yu T 1999 Phys. Rev. Lett. 83 4909
[28] Megier N Strunz W T Viviescas C Luoma K 2018 Phys. Rev. Lett. 120 150402
[29] Zhao X Hedemann S R Yu T 2013 Phys. Rev. 88 022321
[30] Jing J Yu T 2010 Phys. Rev. Lett. 105 240403
[31] Milord L Gerelli E Jamois C Harouri A Chevalier C Viktorovitch P Letartre X Benyattou T 2015 Appl. Phys. Lett. 106 121110
[32] Tan H T Zhang W M Li G X 2011 Phys. Rev. 83 062310
[33] Strunz W T Yu T 2004 Phys. Rev. 69 052115
[34] Yu T Diósi L Gisin N Strunz W T 2000 Phys. Lett. 265 331
[35] Fuchs C A Caves C M 1994 Phys. Rev. Lett. 73 3047
[36] Breuer H P Laine E M Piilo J 2009 Phys. Rev. Lett. 103 210401
[37] Breuer H P Laine E M Piilo J Vacchini B 2016 Rev. Mod. Phys. 88 021002
[38] Hou S C Liang S L Yi X X 2015 Phys. Rev. 91 012109
[39] Ali M M Lo P Y Tu M W Y Zhang W M 2015 Phys. Rev. 92 062306
[40] Cheng J Zhang W Z Zhou L Zhang W 2016 Sci. Rep. 6 23678
[41] Mu Q Zhao X Yu T 2016 Phys. Rev. 94 012334
[42] Golkar M K T 2018 Chin. Phys. 27 40303
[43] Wen J 2018 Chin. Phys. Lett. 35 060301
[44] Ashkin A 1970 Phys. Rev. Lett. 24 156