Quantum entanglement maintained by virtual excitations in an ultrastrongly-coupled-oscillator system.

We study the effect of quantum entanglement maintained by virtual excitations in an ultrastrongly-coupled harmonic-oscillator system. Here, the quantum entanglement is caused by the counterrotating interaction terms and hence it is maintained by the virtual excitations. We obtain the analytical expression for the ground state of the system and analyze the relationship between the average excitation numbers and the ground-state entanglement. We also study the entanglement dynamics between the two oscillators in both the closed- and open-system cases. In the latter case, the quantum master equation is microscopically derived in the normal-mode representation of the coupled-oscillator system. This work will open a route to the study of quantum information processing and quantum physics based on virtual excitations.

Introduction

The ultrastrong coupling (USC) physics[1,2] has recently attracted much attention from the communities of quantum physics, quantum optics, and condensed matter physics. Great advances have been made in both theory[3-15] and experiments in various physical platforms, including semiconductor cavity quantum electrodynamical (QED) systems[16-18], superconducting circuit-QED systems[19-24], coupled photon-2D-electron-gas[25-27], light-molecule[28, 29], and photon-magnon systems[30]. In the USC regime[1, 2], the coupling strength is comparable to the transition frequencies in the system, and then the rotating-wave approximation (RWA) is invalid, namely the counterrotating (CR) terms should be kept in the interactions. It has been demonstrated that the CR terms could produce some unpredictable physical phenomena[3] and have wide applications in quantum information processing[31, 32]. In particular, the development of the ultrastrong coupling field promotes various studies in quantum optics topics beyond the RWA such as the quantum Rabi model[33-39].

One of the interesting effects associated with the CR terms in the USC regime is the generation of virtual excitations. In the presence of the CR terms, the ground states of some typical quantum systems possess virtual excitations. For example, in the quantum Rabi model, it has been shown that virtual photons exist in the ground state[11]. These virtual photons cannot be detected directly even if this absorber is placed inside the cavity, except with very small probability at short times set by the time-energy uncertainty[40]. On the basis of these properties, the ground-state photons in the USC system are considered virtual photons[2]. However, even though these virtual photons cannot be detected directly, there are still ways to probe them. One proposal is to measure the change that they produce in the Lamb shift of an ancillary probe qubit coupled to the cavity[41]. Another proposal is to detect the radiation pressure that they give rise to if the cavity is an optomechanical system[42]. These proposals rely on the rapid modulation of either g (light–matter coupling strength) or the atomic frequency. Then the virtual photons can be converted into real ones and extracted from the system[3, 4, 9, 43–48].

In this paper, we propose to study another quantum effect, quantum entanglement, associated with the virtual excitations. Here the quantum entanglement is created by the CR terms and hence it is maintained by the virtual excitations. We note that the relationship between quantum entanglement and the CR terms has been previously considered in the quantum Rabi model[8]. In addition, the role of the CR terms in the creation of entanglement between two atoms has been investigated in Ref.[49]. We consider an ultrastrongly-coupled two-harmonic-oscillator system. We study the ground state entanglement of the two oscillators and analyze the average excitation numbers in the system. We also study the entanglement dynamics of the system when it is initially in the zero-excitation state and hence all the excitations are created by the CR terms. The influence of the environment dissipations on the system is analyzed based on a microscopically derived quantum master equation in the normal-mode representation.

The rest of this paper is organized as follows. Firstly, we present the physical model of two coupled harmonic oscillators and the Hamiltonian, we also analyze the property of the parity chain in this system. Secondly, we obtain the exact analytical eigensystem of the coupled two-oscillator system. Thirdly, the average virtual excitation numbers are calculated analytically and the quantum entanglement of the ground state is analyzed by calculating the logarithmic negativity. Fourthly, we study the dynamics of the average virtual excitation numbers and quantum entanglement between the two oscillators in both the closed- and open-system cases. Finally, we present a brief conclusion.

Results

Model and hamiltonian

We consider an ultrastrong coupling system, in which two harmonic oscillators are ultrastrongly coupled to each other through the so-called “position–position” type interaction (Fig. 1). This system is described by the Hamiltonianwhere () and () are, respectively, the coordinate and momentum operators of the first (second) oscillator with the resonance frequency () and mass , the parameter is the coupling strength between the two oscillators. By expanding the interaction term, Hamiltonian (1) can be expressed aswhere we introduce the renormalized frequencies and coupling strength asBy introducing the following creation and annihilation operatorsHamiltonian (2) becomeswith being a constant term. Here (a) and (b) are, respectively, the creation (annihilation) operators of the two oscillators with the corresponding resonance frequencies and . In Eq. (5), the first two terms and the constant term represent the free Hamiltonian of the two oscillators. The parameter denotes the coupling strength between the two oscillators. We note that this interaction includes both the rotating-wave and CR terms. In general, in the case of weak coupling and near resonance, the rotating-wave approximation can be made by discarding the CR terms. In this paper, we consider the ultrastrong-coupling case in which the CR terms cannot be discarded. In the presence of the CR terms, the ground state of the system will include excitations and hence quantum entanglement will exist in the ground state. Note that an ultrastrongly-coupled two-mode system has recently been realized in superconducting circuits[50].

Figure 1

Schematic diagram of the coupled two-harmonic-oscillator system. Two harmonic oscillators with resonance frequencies and are coupled to each other via a “position–position” type interaction with the coupling strength g. The parameters and are the decay rates associated with the heat baths in contacted with the oscillators a and b, respectively.

In this two-oscillator system, we introduce the parity operator as , which has the standard properties of a parity operator, such as , , and [7, 51]. The Hamiltonian H in Eq. (5) remains invariant under the transformation , based on the relations , , , and . The Hilbert space of the system can be divided into two subspaces with different parities: odd and even. The basis states of the odd- and even-parity subspaces are, respectively, given byandThe eigenvalues of the parity operator P corresponding to the odd and even parity states are and 1, respectively.

Eigensystem of the coupled two-oscillator system

To study the quantum entanglement of the eigenstates, we need to diagonalize the Hamiltonian H in Eq. (2). To this end, we introduce the transformation operator[52]where the mixing angle is defined byIn terms of the transformation, the Hamiltonian in Eq. (2) can be diagonalized aswhere the resonance frequencies are defined bywith . By introducing the annihilation and creation operatorsHamiltonian (10) can be expressed asThe relations between the operators A (), B (), a (), and b () can be obtained asHere the concrete forms of coefficients () have been given by

Based on Eq. (13), we know the eigenstates of the system in the representation associated with and aswhere the eigenvalues are given byIt is obvious that the ground state of the two-oscillator system is . To study the virtual excitations in the system, we need to know the eigenstates which are expressed in the representation associated with and . It implies that we need to diagonalize the Hamiltonian in the representation of a and b. To this end, we express Hamiltonian (10) with the bosonic creation (annihilation) operators (a) and (b) aswhere we introduce the coupling strengthsTo diagonalize Hamiltonian (18), we introduce the squeezing operatorswhere the two real squeezing parameters are defined byThe transformed Hamiltonian can be written aswhere the unitary operator U can be expressed with the operators a and b asThe eigenstates of the Hamiltonian can be obtained aswhere the eigenvalues are defined in Eq. (17). The eigensystem of the Hamiltonian H can be obtained asAs a result, the ground state of the system can be expressed asIn general, it is hard to write out the ground state in the number state representation. However, we can obtain a number-state expansion of the ground state in the degenerate two-oscillator case[53], i.e., . In this case, we have and the ground state becomesBy expanding the squeezing operators, we then haveIn terms of the relations we then obtainIt can be seen that the superposition components in the ground state are even parity states. This property can be confirmed because the transform U conserves the excitation number and the squeezing operators change the excitation number two by two, without changing the parity.

Ground-state entanglement and quadrature squeezing

We study the ground-state entanglement in this system by calculating the logarithmic negativity. For the two-oscillator system, if the coupling is sufficiently weak, i.e., , the interaction Hamiltonian between the two oscillators can be reduced by the RWA as , which conserves the number of excitations. In this case, the ground state of the system is a trivial direct product of two vacuum states , which does not contain excitations. In the presence of the CR terms, the is not an eigenstate of the system and the ground state will possess excitations. Below, we use numerical method to obtain the ground state of Hamiltonian (5) and calculate the ground state entanglement between the two oscillators. In the presence of the CR terms, the ground state of the two-oscillator system can be expressed aswhere these superposition coefficients are given by , which should be solved numerically. The and reveal that the ground state of the system contains excitations. These excitations in the ground state are called virtual excitations because these excitations cannot be extracted from the system.

The effect of the virtual excitations can be seen from the probability amplitudes in the ground state. The distribution of these probability amplitudes can also exhibit the parity of the ground state. As the ground state is an even parity state, and hence these probability amplitudes associated with the odd parity basis states will disappear. In Fig. 2, we show the absolute values of these probability amplitudes . Here we can see that the values of decrease with the increase of m and n and that there is a symmetric relation . In addition, the values of these odd-parity probability amplitudes with being an odd number are zero, which is a consequence of the fact that the ground state is an even-parity state.

Figure 2

The absolute values of the probability amplitudes for the ground state G in the degenerate two-oscillator case when the coupling strength takes (a) and (b) .

We also calculate the average excitation numbers and in the ground state aswhere we have used the formula,

In Fig. 3a, we show the average excitation numbers and in the ground state as functions of the scaled coupling strength in the degenerate oscillator case . These results show that the average excitation numbers of the two modes are identical (two curves overlap each other). This is because the corotating terms conserve the excitations and the CR terms create simultaneously the excitations in the two modes. The average excitation numbers increase with the coupling strength since a larger coupling strength corresponds to a faster excitation creation.

Figure 3

(a) The average excitation numbers , and (b) the logarithmic negativity in the ground state of the degenerate two-oscillator system as functions of the ratio .

The degree of entanglement between the two oscillators a and b in the ground state of the system can be obtained by calculating the logarithmic negativity. Combining with Eq. (31), the density matrix of the ground state can be written asThe degree of entanglement of the ground state can be quantized by calculating the logarithmic negativity[54, 55]. For a bipartite system described by the density matrix , the logarithmic negativity can be defined bywhere denotes the partial transpose of the density matrix of the system with respect to the oscillator b, and the trace norm is defined byUsing Eqs. (34), (35), and (36), the logarithmic negativity of ground state of the two coupled oscillators can be obtained. In Fig. 3b, we show the logarithmic negativity N as a function of the coupling parameter . The curve shows that the degree of entanglement between the two oscillators in the ground state monotonically increases over the entire range of g. This is because the CR terms in Hamiltonian (5) cause the virtual excitations in the ground state of the system and maintain the quantum entanglement between the two oscillators. If the CR terms are discarded, then the ground state of the system becomes a separate state .

We also study the quadrature squeezing in the ground state by calculating the fluctuations of the rotated quadrature operators. We introduce the rotated quadrature operators for the two modes asThe commutation relation of the above two rotated quadrature operators isAccording to the uncertainty relation, we haveThen the quadrature squeezing appears along the angle if the variances of the rotated quadrature operators satisfy the relation[56]For the ground state given in Eq. (27), the variances of the rotated quadrature operators can be obtained asWhen we exchange the subscripts a and b in Eq. (41), the expression does not change for a given rotating angle. This means that, in the resonance case , the squeezing is the same for the two bosonic modes in the ground state. This point can also be seen from Hamiltonian (5), which is symmetric under the exchange of the subscripts and operators for the two modes in the resonance case.

In Fig. 4a, we show the variance as a function of the rotating angle in the resonance case . Here we can see that the variance is periodic function of and that the minimum of is obtained at and . Note that in the present case . We also show the variance as a function of the coupling strength  in the resonance case , as shown in Fig. 4b. We observe that the squeezing increases with the scaled coupling strength . This is because the quadrature squeezing is caused by the CR interaction terms.

Figure 4

(a) The variance of the rotated quadrature operators as a function of the angle in the resonant case when , 0.2, and 0.4. (b) The variance as a function of the coupling strength in the resonant case .

Dynamics of quantum entanglement

The phenomenon of quantum entanglement accompanied with virtual excitations can also be seen by analyzing the entanglement dynamics of the system. We consider the case in which the system is initially in the zero-excitation state . In the closed-system case, a general state of the system can be written asBy substituting Eqs. (5) and (42) into the Schrödinger equation, the equations of motion for these probability amplitudes are obtained asFor the initial state , the initial condition of these probability amplitudes reads . By numerically solving Eq. (43) under this initial condition, the evolution of these probability amplitudes can be obtained. Using Eqs. (35), (36), and (42), we can calculate numerically the average excitation numbers and and the logarithmic negativity of the state .

In Fig. 5a, we show the time evolution of the average excitation numbers and in modes a and b. Here we can see that, similar to the ground state case, the average excitation numbers in the two modes are identical (the two curves overlap each other). In addition, the average excitation numbers experience a periodic oscillation. In Fig. 5b, we show the time dependence of the logarithmic negativity N(t) of the state . The curve shows that logarithmic negativity between the two oscillators also experiences a periodic oscillation. Here we choose the initial state of the system as , the existence of the CR terms still causes the appearance of virtual excitations, which leads to entanglement between the two oscillators. This result is different from that in the RWA case in which the CR terms are discarded in the two oscillators under the same initial state. When we discard the CR terms and choose the initial state as , which is the eigenstate of the corotating interaction term , the system will stay this state. Then there are no virtual excitations in the system and no quantum entanglement between the two oscillators.

Figure 5

Dynamics of (a) the average excitation numbers , and (b) the logarithmic negativity when the degenerate two-oscillator system is initially in the state . The used parameter is .

We also study the influence of the environment dissipations on the dynamics of the system. As we consider the ultrastrong-coupling regime of the coupled system, we derive the quantum master equation in the normal-mode representation of these two coupled oscillators. We employ the standard Born–Markov approximation under the condition of weak system-bath couplings and short bath correlation times to derive the quantum master equation. The secular approximation is made by discarding these high-frequency oscillating terms including , , and . The quantum master equation in the normal-mode representation of Hamiltonian (5) can be written aswith the effective decay ratesIn Eq. (44), the dissipatoris the standard Lindblad superoperator that describes the dampings of the oscillators. The parameters and are the decay rates relating to the heat bath in contact with the oscillators a and b, respectively. Here we consider the zero temperature environments such that the thermal excitation effect can be excluded.

In our simulations, we need to calculate the evolution of the system in the bare-mode representation, i.e., . The relationship between the density matrix in the normal-mode representation ( ) and the density matrix in the bare-mode representation () is determined by the transformation . Combining with Eq. (10), the quantum master equation in the bare-mode representation can be obtained asThe transformation relations between the operators and arewhere the forms of these coefficients () are given by By substituting Eqs. (46), (48), and (49) into Eq. (47), we obtainwhere we introduce the superoperator asThe coefficients introduced in Eq. (50) are defined by Note that the cross terms between the two modes a and b in Eq. (50) are induced by the interaction between the two oscillators.

For below calculations, we express the density matrix of the two-oscillator system in the bare-mode representation aswith the density matrix elements . For an initial state , the nonzero density matrix element is . By numerically solving Eq. (47) under the initial condition, the time evolution of the density matrix can be obtained.

Below we study the dynamics of the average excitation numbers and quantum entanglement in this system. Based on Eq. (47), the expressions of the average excitation numbers and can be expanded asTherefore, the average excitation numbers and can be obtained by solving the equations of motion for these density matrix elements in the number-state representation.

In Fig. 6a, the dynamics of the average excitation numbers and is shown in the open-system case with different time t. We observe that the two excitation numbers and overlap each other and initially experience a large oscillation. With the increase of time t, the oscillation amplitudes of the average excitation numbers decrease gradually. In the long-time limit , the average excitation numbers will reach steady values due to the dissipations.

Figure 6

Dynamics of (a, c) the average excitation numbers , and (b, d) the logarithmic negativity as functions of the evolution time t when the system is initially in the state . The parameters used are and . The results in panels (a) and (b) [(c) and (d)] are calculated with the microscopic quantum master equation (the phenomenological quantum master equation).

The entanglement of the density matrix can be quantified by calculating the logarithmic negativity. In terms of Eqs. (35), (47), and (53), the logarithmic negativity of the state can be obtained numerically. In Fig. 6b, we show the logarithmic negativity N(t) of the density matrix versus the time t. The result shows that the logarithmic negativity oscillates very fast due to the free evolution of the system. We also find that the envelope of the logarithmic negativity converges gradually with the evolution time t and eventually reaches a stable value due to the dissipations. The time scale of the oscillation-pattern decay for the logarithmic negativity is very similar to that of the excitations created by the CR interaction terms. In particular, we find that there exists steady-state entanglement due to the presence of the CR interaction terms in this system.

In this work, we consider the ultrastrong-coupling regime and hence the quantum master equation is derived in the normal mode representation. For comparison, we show in Fig. 6c,d the evolution of the average excitation numbers and the logarithmic negativity calculated by solving the phenomenological quantum master equation, which is obtained by adding the dissipators of two free bosonic modes into the Liouville equation,The initial state of the system is the same as that considered in the microscopic quantum master equation. We see from Fig. 6 that, for the average excitation numbers, though these results can approach steady-state values, the envelop and the oscillation amplitude are different for the results obtained with two different quantum master equations. However, for the logarithmic negativity, we find that the difference between the two results exists but is small when . We checked the fact that the difference will increase as the increase of the ratio . Therefore, the microscopic quantum master equation should be used in the ultrastrongly-coupled-oscillator system.

Conclusion

In conclusion, we have studied quantum entanglement in an ultrastrongly-coupled two-harmonic-oscillator system. Concretely, we have studied the ground-state entanglement by calculating the logarithmic negativity of the ground state. Here, the quantum entanglement is maintained by the virtual excitations generated by the CR terms and bounded in the ground state. We have also studied the dynamics of quantum entanglement of the system. By microscopically deriving a quantum master equation in the normal-mode representation of the two oscillators, we analyzed the influence of the dissipations on the entanglement dynamics and found that there exists steady-state entanglement in this system.