text:
# Decoherence of Quantum Damped Oscillators

## I Introduction

Physical systems are believed without doubt to obey quantum theory that has been tested in various areas. Microscopic systems are described most accurately by the quantum theory, whereas most of the macroscopic systems are described also by classical theory to any desired precision. The most prominent feature of quantum theory is quantum coherence, the interference among superposed states. Then the question may be raised: how quantum states of a system can have the characteristic features of classical theory such as quantum decoherence (loss of quantum coherence) and classical correlation along classical trajectory? This quantum-to-classical transition or classicality of quantum system has been an important problem since the advent of quantum theory. Though not completely settled, it is well known that the interaction of a quantum system with an environment can lead to classicality of quantum decoherence and classical correlation . In particular, dissipation is known to result in quantum decoherence . In most of literature, the quantum decoherence has been studied for a system coupled to an environment or thermal bath with many degrees of freedom.

In this paper we study quantum damped oscillators as a simple model for the dissipative system. There are different Hamiltonian representations for this damped oscillator. One representation is the Caldirola-Kanai (CK) oscillator, which is a one-dimensional system with an exponentially increasing mass . This oscillator is an open system because its parameters such as mass or frequency depend explicitly on time. The other representation is the Bateman or Feshbach-Tikochinsky (BFT) oscillator, which consists of a damped oscillator and an amplified oscillator . The second oscillator is a closed system as the total energy is conserved and the energy dissipated from the damped oscillator is transferred to the amplified one. These quantum damped oscillators have been studied intensively as a model to understand dissipation in quantum theory and the connection between them has been found . The CK oscillator has also been investigated to find the characteristic features of quantum states  (for review and references, see Ref. ). The BFT oscillator has also been used to study quantum dissipation .

The main purpose of this paper is to study the classicality of the damped CK and BFT oscillators as a quantum mechanical system with dissipation and few degrees of freedom. More concretely we investigate quantum decoherence, the necessary condition, and classical correlation, the sufficient condition for classicality. For that purpose we first find the density matrix of the Gaussian state for the CK oscillator and the reduced density matrix for the damped part of the BFT oscillator and then apply the criterion on measuring quantitatively quantum decoherence and classical correlation. The measure of quantum decoherence is defined as the ratio of the diagonal element to the off-diagonal element of the (reduced) density matrix . It is found that the density matrix for the CK oscillator shows no quantum decoherence, whereas the reduced density matrix for the BFT oscillator shows both quantum decoherence and classical correlation for the damped and amplified parts.

The organization of this paper is as follows. In Sec. II we find the quantum state and density matrix of CK oscillator and calculate the measure of quantum decoherence and classical correlation. In Sec. III we find the density matrix and thereby the reduced density matrix of the BFT oscillator. The measure of decoherence is evaluated. In Sec. IV we study the amplified oscillator, the opposite case of damped oscillator and investigate the condition for classicality.

## II CK Oscillator

The CK oscillator is an an open system with the variable mass $`m(t)=me^{\gamma t/m}`$

$$H_{\mathrm{CK}}=\frac{1}{2me^{\gamma t/m}}p_x^2+\frac{m\omega ^2e^{\gamma t/m}}{2}x^2,$$
(1)
and has the classical equation of motion for dissipative system

$$\ddot{x}+\frac{\gamma }{m}\dot{x}+\omega ^2x=0.$$
(2)
It is assumed that the quantum theory of the damped oscillator is prescribed by the time-dependent Schrödinger equation

$$i\mathrm{}\frac{}{t}\mathrm{\Psi }(x,t)=\widehat{H}_{\mathrm{CK}}(t)\mathrm{\Psi }(x,t).$$
(3)
The wave functions are found in various methods  (for review and references, see Ref. ). Lewis and Riesenfeld introduced an invariant operator for the general time-dependent oscillator, whose eigenstate is an exact quantum state up to a time-dependent phase factor .

Following Refs. , we introduce a pair of first order operators in position and momentum

$`\widehat{a}(t)`$ $`=`$ $`i[u^{}(t)\widehat{p}_x\dot{u}^{}(t)\widehat{x}],`$ (4)
$`\widehat{a}^{}(t)`$ $`=`$ $`i[u(t)\widehat{p}_x\dot{u}(t)\widehat{x}],`$ (5)
and require them to satisfy the quantum Liouville-von Neumann equation

$`i\mathrm{}{\displaystyle \frac{}{t}}\widehat{a}(t)+[\widehat{a}(t),\widehat{H}_{\mathrm{CK}}(t)]=0,`$ (6)
$`i\mathrm{}{\displaystyle \frac{}{t}}\widehat{a}^{}(t)+[\widehat{a}^{}(t),\widehat{H}_{\mathrm{CK}}(t)]=0.`$ (7)
Then $`u`$ satisfies the classical equation of motion (2). The Wronskian condition

$$\mathrm{}me^{\gamma t/m}\left[\dot{u}^{}(t)u(t)\dot{u}(t)u^{}(t)\right]=i,$$
(8)
guarantees the standard commutation relation for all times

$$[\widehat{a}(t),\widehat{a}^{}(t)]=1.$$
(9)
We are interested in the underdamped motion given by

$$u(t)=\frac{1}{\sqrt{2\mathrm{}m\mathrm{\Omega }}}e^{\gamma t/(2m)}e^{i\mathrm{\Omega }t},$$
(10)
where

$$\mathrm{\Omega }=\sqrt{\omega ^2\left(\frac{\gamma }{2m}\right)^2},(\gamma 2m\omega ).$$
(11)
The number operator defined by

$$\widehat{N}(t)=\widehat{a}^{}(t)\widehat{a}(t)$$
(12)
also satisfies Eq. (7) and yields the number state as an exact quantum state

$$\widehat{N}(t)|n,t=n|n,t.$$
(13)
The wave function for the number state that satisfies Eq. (3) is given by

$$\mathrm{\Psi }_n(x,t)=\left(\frac{m\mathrm{\Omega }e^{\gamma t/m}}{\pi \mathrm{}}\right)^{1/4}\frac{e^{i(n+1/2)\mathrm{\Omega }t}}{\sqrt{2^nn!}}H_n\left(\sqrt{\frac{m\mathrm{\Omega }e^{\gamma t/m}}{\mathrm{}}}x\right)\mathrm{exp}\left[e^{\gamma t/m}\left(\frac{m\mathrm{\Omega }}{2\mathrm{}}+i\frac{\gamma }{4\mathrm{}}\right)x^2\right],$$
(14)
where $`H_n`$ is the Hermite polynomial. The wave function has the dispersion relations

$`\widehat{x}^2`$ $`=`$ $`\mathrm{}^2u^{}(t)u(t)={\displaystyle \frac{\mathrm{}}{2m\mathrm{\Omega }}}e^{\gamma t/m},`$ (15)
$`\widehat{p}_x^2`$ $`=`$ $`\mathrm{}^2m^2(t)\dot{u}^{}\dot{u}={\displaystyle \frac{\mathrm{}m\omega ^2}{2\mathrm{\Omega }}}e^{\gamma t/m}.`$ (16)
The more the system dissipates in time, the more the wave function becomes sharply peaked around $`x=0`$, whereas the wave function is more dispersed in momentum space. The uncertainty is a constant

$$(\mathrm{\Delta }x)(\mathrm{\Delta }P_x)=\frac{\mathrm{}\omega }{2\mathrm{\Omega }}\frac{\mathrm{}}{2},$$
(17)
and the Hamiltonian expectation value is also a constant.

$$0,t|\widehat{H}_{\mathrm{CK}}(t)|0,t=\frac{\mathrm{}\omega ^2}{\mathrm{\Omega }}\frac{\mathrm{}\omega }{2}.$$
(18)
The density matrix of the Gaussian wave function (14), the ground state with $`n=0`$, has the form

$`\rho _{\mathrm{CK}}(x^{},x,t)`$ $`=`$ $`\mathrm{\Psi }_0(x^{})\mathrm{\Psi }_0^{}(x)`$ (19)
$`=`$ $`\left({\displaystyle \frac{m\mathrm{\Omega }e^{\gamma t/m}}{\pi \mathrm{}}}\right)^{1/2}\mathrm{exp}\left[\mathrm{\Gamma }_cx_c^2\mathrm{\Gamma }_\delta x_\delta ^2\mathrm{\Gamma }_\mu x_cx_\delta \right],`$ (20)
where

$$x_c=\frac{1}{2}(x^{}+x),x_\delta =\frac{1}{2}(x^{}x),$$
(21)
and

$$\mathrm{\Gamma }_c=\mathrm{\Gamma }_\delta =\frac{m\mathrm{\Omega }}{\mathrm{}}e^{\gamma t/m},\mathrm{\Gamma }_\mu =i\frac{\gamma }{\mathrm{}}e^{\gamma t/m}.$$
(22)
The off-diagonal element, the coefficient of $`x_\delta ^2`$, measures the degree of quantum coherence, i.e., the interference between two different trajectories. The representation-independent measure of quantum decoherence  is now given by

$$\delta _{\mathrm{QD}}=\frac{1}{2}\sqrt{\frac{\mathrm{\Gamma }_c}{\mathrm{\Gamma }_\delta }}=\frac{1}{2}.$$
(23)
The condition for quantum decoherence $`(\delta _{\mathrm{QD}}1)`$ shows no decoherence for the CK oscillator. Likewise, the measure of classical correlation

$$\delta _{\mathrm{CC}}=\sqrt{\frac{\mathrm{\Gamma }_c^2\mathrm{\Gamma }_\delta ^2}{\mathrm{\Gamma }_\mu ^{}\mathrm{\Gamma }_\mu }}=\frac{(m\mathrm{\Omega })^2}{\mathrm{}\gamma }e^{\gamma t/m},$$
(24)
shows no classical correlation conditioned by $`\delta _{\mathrm{CC}}1`$, except for the case of $`\mathrm{\Omega }0`$ in the large-damping limit $`\gamma 2m\omega `$. Even in this case the exponentially growing factor dominates at later times and classical correlation is lost. The zero-damping limit $`(\gamma =0)`$ of a pure harmonic oscillator does not lead to any classical correlation with the infinite $`\delta _{\mathrm{CC}}`$ as expected. The CK oscillator does achieve neither quantum decoherence nor classical correlation.

In summary, the CK oscillator does not have the genuine properties of dissipative systems, though its equation of motion does show such a damping effect. First, the energy defined by the expectation value of the Hamiltonian operator does not have any damping factor that implies the dissipation of energy. Second, the uncertainty does not grow as the evolution proceeds. Third, there is neither quantum decoherence nor classical correlation regardless of the magnitude of damping factor.

## III BFT Oscillator

The one-dimensional CK oscillator has a constant expectation value of the Hamiltonian. To be a genuine dissipative system, the energy of the damped subsystem of the system must be dissipated away and transferred to another subsystem. This means that the damped oscillator may be properly described by a two-dimensional system, one subsystem of which dissipates the energy and the transferred energy amplifies the other subsystem. Such a model has been suggested long ago by Bateman  and later by Feshbach and Tikochinsky . The BFT oscillator is described by the Lagrangian

$$L=m\dot{x}\dot{y}+\frac{\gamma }{2}(x\dot{y}\dot{x}y)kxy.$$
(25)
The one subsystem with $`x`$ variable obeys the damped equation of motion (2) where $`\omega =\sqrt{k/m}`$. The other subsystem obeys the equation for an amplified oscillator

$$\ddot{y}\frac{\gamma }{m}\dot{y}+\omega ^2y=0.$$
(26)
The energy of $`y`$ increases as $`e^{\gamma t/m}`$. The Hamiltonian is given by

$$H_\gamma =\frac{1}{m}p_xp_y+\frac{\gamma }{2m}(yp_yxp_x)+\mathrm{\Omega }^2xy,$$
(27)
where $`p_x=m\dot{y}`$ and $`p_y=m\dot{x}`$.

In the limit of zero-dissipation $`(\gamma =0)`$, the BFT oscillator is the sum of two decoupled oscillators with opposite signs

$$H_0=\frac{1}{2m}p_\xi ^2+\frac{k}{2}\xi ^2\frac{1}{2m}p_\zeta ^2\frac{k}{2}\zeta ^2,$$
(28)
where

$$\xi =\frac{1}{\sqrt{2}}(x+y),\zeta =\frac{1}{\sqrt{2}}(x+y).$$
(29)
The ground state of each oscillator leads to the zero energy and a density matrix

$`\rho (x^{},y^{},x,y)`$ $`=`$ $`\left({\displaystyle \frac{m\omega }{\pi \mathrm{}}}\right)\mathrm{exp}\left[{\displaystyle \frac{m\omega }{2\mathrm{}}}(\xi ^2+\zeta ^2+\xi ^2+\zeta ^2)\right]`$ (30)
$`=`$ $`\left({\displaystyle \frac{m\omega }{\pi \mathrm{}}}\right)\mathrm{exp}\left[{\displaystyle \frac{m\omega }{2\mathrm{}}}(x^2+y^2+x^2+y^2)\right].`$ (31)
For the dissipation case $`(\gamma 0)`$, the density matrix satisfies the quantum Liouville-von Neumann equation, whose coordinate representation is given by

$`i\mathrm{}{\displaystyle \frac{t}{t}}\rho (x^{},y^{},x,y)=[{\displaystyle \frac{\mathrm{}^2}{m}}({\displaystyle \frac{^2}{x^{}y^{}}}{\displaystyle \frac{^2}{xy}})+i{\displaystyle \frac{\mathrm{}\gamma }{2m}}(x^{}{\displaystyle \frac{}{x^{}}}y^{}{\displaystyle \frac{}{y^{}}}+x{\displaystyle \frac{}{x}}y{\displaystyle \frac{}{y}})`$ (32)
$`+m\mathrm{\Omega }^2(x^{}y^{}xy)]\rho (x^{},y^{},x,y).`$ (33)
The density matrix quadratic in $`x^{},y^{},x,y`$ has the general form

$`\rho (x^{},y^{},x,y)=N\mathrm{exp}[e^{\gamma t/m}(A^{}x^2+A_1x^{}x+Ax^2)e^{\gamma t/m}(B^{}y^2+B_1y^{}y+By^2)`$ (34)
$`C(x^{}y^{}+xy)(Dx^{}y+D^{}xy^{})],`$ (35)
where the coefficients satisfy the set of equations

$`\dot{A}`$ $`=`$ $`i{\displaystyle \frac{\mathrm{}}{m}}(2ACA_1D^{}),`$ (36)
$`\dot{B}`$ $`=`$ $`i{\displaystyle \frac{\mathrm{}}{m}}(2BCB_1D),`$ (37)
$`\dot{A}_1`$ $`=`$ $`i{\displaystyle \frac{2\mathrm{}}{m}}(ADA^{}D^{}),`$ (38)
$`\dot{B}_1`$ $`=`$ $`i{\displaystyle \frac{2\mathrm{}}{m}}(BD^{}B^{}D),`$ (39)
$`\dot{C}`$ $`=`$ $`i{\displaystyle \frac{\mathrm{}}{m}}(4AB+C^2A_1B_1D^{}D)i{\displaystyle \frac{m\mathrm{\Omega }^2}{\mathrm{}}},`$ (40)
$`\dot{D}`$ $`=`$ $`i{\displaystyle \frac{2\mathrm{}}{m}}(A_1BA^{}B_1).`$ (41)
There is a symmetry under $`xy`$ and $`\gamma \gamma `$, which is in fact the time-reversal symmetry of the Hamiltonian (27). Hence we can set $`A=B^{}`$ and $`A_1=B_1^{}`$. Also in the zero-dissipation limit, by comparing the density matrices (31) and (35) we find that $`A_1,C,D`$ approach zero. A particular solution is found to be

$`A(\gamma )`$ $`=`$ $`B^{}(\gamma )=\left[\left({\displaystyle \frac{m\mathrm{\Omega }}{2\mathrm{}}}\right)^2+{\displaystyle \frac{D^{}(\gamma )D(\gamma )}{4}}\right]^{1/2}e^{i(\pi \theta )},`$ (42)
$`D(\gamma )`$ $`=`$ $`|D(\gamma )|e^{i\theta },`$ (43)
$`A_1`$ $`=`$ $`B_1=C=0,`$ (44)
where

$$D(\gamma 0)0.$$
(45)
Now the density matrix from the particular solution (44) becomes

$$\rho (x^{},y^{},x,y)=N\mathrm{exp}\left[e^{\gamma t/m}(A^{}x^2+Ax^2)e^{\gamma t/m}(Ay^2+A^{}y^2)\left(Dx^{}y+D^{}xy^{}\right)\right],$$
(46)
where

$$N=\frac{1}{\pi }\left[(A+A^{})^2\frac{1}{4}(D+D^{})^2\right]^{1/2}.$$
(47)
Letting $`y^{}=y`$ and integrating over $`y`$, one obtains the reduced density matrix for the damped subsystem of $`x^{}`$ and $`x`$:

$$\rho _{\mathrm{red}}(x^{},x)=N_1\mathrm{exp}\left[e^{\gamma t/m}\left\{(A^{}x^2+Ax^2)\frac{(Dx^{}+D^{}x)^2}{4(A+A^{})}\right\}\right],$$
(48)
where

$$N_1=N\times \sqrt{\frac{\pi e^{\gamma t/m}}{(A+A^{})}}.$$
(49)
Similarly, the reduced density matrix for the amplified subsystem is obtained by using the symmetry $`xy`$ and $`\gamma \gamma `$:

$$\rho _{\mathrm{red}}(y^{},y)=N_2\mathrm{exp}\left[e^{\gamma t/m}\left\{(Ay^2+A^{}y^2)\frac{(Dy+D^{}y^{})^2}{4(A+A^{})}\right\}\right],$$
(50)
where

$$N_2=N\times \sqrt{\frac{\pi e^{\gamma t/m}}{(A+A^{})}}$$
(51)
Then the reduced density matrix (48) is written in the form

$$\rho _{\mathrm{red}}(x^{},x)=N_1\mathrm{exp}\left[\mathrm{\Gamma }_cx_c^2\mathrm{\Gamma }_\delta x_\delta ^2\mathrm{\Gamma }_\mu x_cx_\delta \right],$$
(52)
where

$`\mathrm{\Gamma }_c(|D|,\theta )`$ $`=`$ $`e^{\gamma t/m}\left[(A+A^{}){\displaystyle \frac{(D+D^{})^2}{4(A+A^{})}}\right],`$ (53)
$`\mathrm{\Gamma }_\delta (|D|,\theta )`$ $`=`$ $`e^{\gamma t/m}\left[(A+A^{}){\displaystyle \frac{(DD^{})^2}{4(A+A^{})}}\right],`$ (54)
$`\mathrm{\Gamma }_\mu (|D|,\theta )`$ $`=`$ $`e^{\gamma t/m}(2)\left[(AA^{})+{\displaystyle \frac{D^2D^2}{4(A+A^{})}}\right].`$ (55)
As the density matrix (46) and the reduced one (48) depend only two real parameters $`|D|`$ and $`\theta `$, the measure of quantum decoherence is given by

$`\delta _{\mathrm{QD}}`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{{\displaystyle \frac{\mathrm{\Gamma }_c(|D|,\theta )}{\mathrm{\Gamma }_\delta (|D|,\theta )}}}`$ (56)
$`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{{\displaystyle \frac{(m\mathrm{\Omega }/\mathrm{})^2}{(m\mathrm{\Omega }/\mathrm{})^2\mathrm{cos}^2\theta +|D|^2}}}.`$ (57)
and that of classical correlation by

$`\delta _{\mathrm{CC}}`$ $`=`$ $`\sqrt{{\displaystyle \frac{\mathrm{\Gamma }_c^2(|D|,\theta )\mathrm{\Gamma }_\delta ^2(|D|,\theta )}{\mathrm{\Gamma }_\mu ^{}(|D|,\theta )\mathrm{\Gamma }_\mu (|D|,\theta )}}}`$ (58)
$`=`$ $`{\displaystyle \frac{1}{2}}|\mathrm{cot}\theta |\left[\left({\displaystyle \frac{m\mathrm{\Omega }}{\mathrm{}}}\right)^2\mathrm{cos}^2\theta +|D|^2\right].`$ (59)
If the condition

$$\theta \frac{\pi }{2},|D|>\frac{m\mathrm{\Omega }}{\mathrm{}}$$
(60)
is satisfied, the density matrices (46) and (48) achieve a significant degree of quantum decoherence $`\delta _{\mathrm{QD}}<1/2`$ and a sufficient degree of classical correlation $`\delta _{\mathrm{CC}}1`$. This is the case of the large dissipation.

In summary, the BFT oscillator has the density matrix which achieves quantum decoherence as well as almost complete classical correlation. This means that the BFT oscillator may be the quantum analog of a classical dissipative oscillator.

## IV Amplified Oscillator

We now consider the opposite case of the damped oscillator, that is, the amplified oscillator, which is physically motivated by an unstable system, for instance, the second order phase transition during the spinodal instability. Quantum decoherence has not been observed for an unstable, exponentially growing, single oscillator . This oscillator has the Hamiltonian

$$H_y=\frac{1}{2me^{\gamma t/m}}p_y^2+\frac{m\omega ^2e^{\gamma t/m}}{2}y^2.$$
(61)
The time-dependent annihilation and creation operators are given by

$`\widehat{b}(t)`$ $`=`$ $`i[v^{}(t)\widehat{p}_y\dot{v}^{}(t)\widehat{y}],`$ (62)
$`\widehat{b}^{}(t)`$ $`=`$ $`i[v(t)\widehat{p}_y\dot{v}(t)\widehat{y}],`$ (63)
where

$$\ddot{v}\frac{\gamma }{m}\dot{v}+\omega ^2v=0.$$
(64)
Then the exponentially growing solution is given by

$$v(t)=\frac{1}{\sqrt{2\mathrm{}m\mathrm{\Omega }_y}}e^{\gamma t/(2m)}e^{i\mathrm{\Omega }t},$$
(65)
where $`\mathrm{\Omega }`$ is given by Eq. (11). The harmonic oscillator wave functions are obtained from Eq. (14) by replacing $`\gamma `$ by $`(\gamma )`$:

$$\mathrm{\Psi }_n(y,t)=\left(\frac{m\mathrm{\Omega }e^{\gamma t/m}}{\pi \mathrm{}}\right)^{1/4}\frac{e^{i(n+1/2)\mathrm{\Omega }t}}{\sqrt{2^nn!}}H_n\left(\sqrt{\frac{m\mathrm{\Omega }e^{\gamma t/m}}{\mathrm{}}}y\right)\mathrm{exp}\left[e^{\gamma t/m}\left(\frac{m\mathrm{\Omega }}{2\mathrm{}}i\frac{\gamma }{4\mathrm{}}\right)y^2\right].$$
(66)
Repeating the steps in Sec. II, the dispersion relations are given by

$`\widehat{y}^2`$ $`=`$ $`{\displaystyle \frac{\mathrm{}}{2m\mathrm{\Omega }_y}}e^{\gamma t/m},`$ (67)
$`\widehat{p}_y^2`$ $`=`$ $`{\displaystyle \frac{\mathrm{}m\omega ^2}{2\mathrm{\Omega }_y}}e^{\gamma t/m}.`$ (68)
So the uncertainty relation and the Hamiltonian expectation value have the same values as Eqs. (17) and (18), respectively. As the ground state wave function $`(n=0)`$ in Eq. (66) is obtained by replacing $`\gamma `$ by $`(\gamma )`$ and all the steps are the same in Sec. II, we find the measures for quantum decoherence and classical correlation

$`\delta _{\mathrm{QD}}`$ $`=`$ $`{\displaystyle \frac{1}{2}},`$ (69)
$`\delta _{\mathrm{CC}}`$ $`=`$ $`{\displaystyle \frac{(m\mathrm{\Omega })^2}{\mathrm{}\gamma }}e^{\gamma t/m}.`$ (70)
The wave functions for the amplified oscillator are certainly classically correlated as shown in Ref. . But this does not mean that the wave functions achieve quantum decoherence because $`\delta _{\mathrm{QD}}=1/2`$.

Now we turn to the amplified subsystem of the BFT oscillator in Sec. III. The reduced density matrix for the amplified subsystem can be written in the form

$$\rho _{\mathrm{red}}(y^{},y)=N_2\mathrm{exp}\left[\mathrm{\Gamma }_c^yy_c^2\mathrm{\Gamma }_\delta ^yy_\delta ^2\mathrm{\Gamma }_\mu ^yy_cy_\delta \right],$$
(71)
where

$`\mathrm{\Gamma }_c^y(|D|,\theta )`$ $`=`$ $`e^{\gamma t/m}\left[(A+A^{}){\displaystyle \frac{(D+D^{})^2}{4(A+A^{})}}\right],`$ (72)
$`\mathrm{\Gamma }_\delta ^y(|D|,\theta )`$ $`=`$ $`e^{\gamma t/m}\left[(A+A^{}){\displaystyle \frac{(DD^{})^2}{4(A+A^{})}}\right],`$ (73)
$`\mathrm{\Gamma }_\mu ^y(|D|,\theta )`$ $`=`$ $`e^{\gamma t/m}(2)\left[(AA^{})+{\displaystyle \frac{D^2D^2}{4(A+A^{})}}\right].`$ (74)
Therefore, for the particular solution (44) the measures for quantum decoherence and classical correlation for the amplified subsystem are given by the same Eqs. (57) and (59) for the damped subsystem. This may be expected from the symmetry of the Hamiltonian under $`xy`$ and $`\gamma \gamma `$. There is also the symmetry of reduced density matrices of $`x`$ and $`y`$ under $`\gamma \gamma `$ and $`D(\gamma )=D^{}(\gamma )`$. When the condition (60) is satisfied, the amplified subsystem also achieves both quantum decoherence and classical correlation. Such classicality of wave functions has been observed for the system of an unstable amplified oscillator coupled to a stable oscillator in a quantum phase transition model . The BFT oscillator provides an exactly solvable model for both quantum decoherence and classical correlation and may shed some light in understanding how the classicality of quantum systems can be achieved.

## V Discussion

We have studied the CK and BFT oscillators as the quantum analogs of a classical dissipative oscillator that obeys the classical damped equation of motion. The CK oscillator is a one-dimensional oscillator with an exponentially increasing mass in time. The wave function of the CK oscillator is sharply peaked around the origin with the exponentially decreasing position-dispersion. However, the momentum-dispersion increases exponentially so that the uncertainty of both position and momentum is constant and has the value greater than the minimum value of Heisenberg uncertainty relation. Similarly, the Hamiltonian expectation value has a constant value greater than the minimum value of a harmonic oscillator without the damping factor. Further, CK oscillator does not achieve both quantum decoherence and classical correlation.

On the other hand, the BFT oscillator is a two-dimensional system, one subsystem describing the damped oscillator and the other describing an amplified oscillator. In this oscillator the energy dissipated away from the damped oscillator is transferred to the other so that the total energy is conserved. We have found two parameter-dependent density matrix (46), whose reduced density matrix for the damped oscillator or the amplified oscillator shows not only quantum decoherence but also classical correlation. The quantum decoherence of the BFT oscillator may be understood as a consequence of the coupling between the damped and amplified modes.

The fact that the interaction or coupling of the system is indispensable for quantum decoherence has been observed in the two-oscillator model with the Hamiltonian

$$H=\frac{1}{2m}(p_1^2+p_2^2)+\frac{m}{2}(\omega _1^2x_1^2+\omega _2^2x_2^2)+\lambda x_1x_2.$$
(75)
The reduced density matrix for $`x_1`$ is given by

$$\rho _{\mathrm{red}}(x_1^{},x_1)=\left(\frac{1}{\pi D}\right)^{1/2}\mathrm{exp}\left[\frac{1}{D}(x_c^2+\mathrm{\Gamma }_\delta x_\delta ^2)\right],$$
(76)
where

$`D`$ $`=`$ $`\mathrm{cosh}\eta \mathrm{sinh}\eta \mathrm{cos}(2\vartheta ),`$ (77)
$`\mathrm{\Gamma }_\delta `$ $`=`$ $`\mathrm{cosh}^2\eta \mathrm{sinh}^2\eta \mathrm{cos}^2(2\vartheta ),`$ (78)
where

$`e^\eta `$ $`=`$ $`{\displaystyle \frac{m(\omega _1^2+\omega _2^2)+\sqrt{m^2(\omega _1^2\omega _2^2)^2+4\lambda ^2}}{2\sqrt{m^2\omega _1^2\omega _2^2\lambda ^2}}},`$ (79)
$`\mathrm{tan}(2\vartheta )`$ $`=`$ $`{\displaystyle \frac{2\lambda }{m(\omega _2^2\omega _1^2)}}.`$ (80)
The measure of decoherence

$$\delta _{\mathrm{QD}}=\frac{1}{2\sqrt{\mathrm{cosh}^2\eta \mathrm{sinh}^2\eta \mathrm{cos}^2(2\vartheta )}}$$
(81)
has a value $`\delta _{\mathrm{QD}}1/2`$, the equality corresponding to the zero mode-mixing, and has small values $`\delta _{\mathrm{QD}}1/2)`$ for large mixing angles $`(\vartheta \pi /4`$ and $`\eta 1)`$. However, the measure of classical correlation is infinite, the same as the quantum state of a simple harmonic oscillator. Therefore we may conclude that a system with dissipation (damping) interacting with another system can achieve both quantum decoherence and classical correlation.

###### Acknowledgements.

S.P.K. and A.E.S. would like to appreciate the warm hospitality of the Theoretical Physics Institute, Univ. of Alberta. S.P.K. also would like to thank D.N. Page for helpful discussions on quantum decoherence and Prof. Y.S. Kim for useful information. F.C.K. was supported by NSERC of Canada, S.P.K. by KRF under Grant No. 2000-015-DP0080, and A.E.S. by CNPq of Brazil.