## Lindblad dynamics of the damped and forced quantum harmonic oscillator

## H. J. Korsch ∗

FB Physik, Technische Universit¨t Kaiserslautern a D-67653 Kaiserslautern, Germany

September 6, 2019

## Abstract

The quantum dynamics of a damped and forced harmonic oscillator is investigated in terms of a Lindblad master equation. Elementary algebraic techniques are employed allowing for example to analyze the long time behavior, i.e. the quantum limit cycle. The time evolution of various expectation values is obtained in closed form as well as the entropy and the Husimi phase space distribution. We also discuss the related description in terms of a non-Hermitian Hamiltonian.

## 1 Introduction

In classical mechanics, the damped and driven harmonic oscillator is one of the paradigmatic systems discussed in elementary physics lectures, modeling for example a mass attached to a linear spring in the gravitational field. The dynamics is linear and allows a closed form solution, which then serves as a basic model for analyzing more general and, of course, much more demanding general dynamical systems. The harmonic oscillator also constitutes a basic building block in elementary quantum mechanics, however almost exclusively in context with Hermitian dynamics, which is also in most cases restricted to systems with no explicit time-dependence, though also a time-depended force

∗ Electronic address: h.j.korsch@gmail.com

can be treated quite elementary, as for instance in the book by Louisell [1, sect. 3.11]. The damped quantum harmonic oscillator is, however, addressed very rarely in elementary texts. Exceptions are, e.g., the recent paper by Fujii [2] describing in much detail the dissipative quantum dynamics in terms of the Lindblad master equations, however for the force-free case. This system is also described in the article [3] as well as in lecture notes by Englert and Morigi [4], where also a quantum-optical system with periodic driving is studied, describing a decaying photon field excited by periodic kicking.

An elementary study of the Lindblad dynamics of a forced harmonic oscillator seems to be missing, however, and here we will try to fill this gap. We follow the analysis in [2], describing briefly these results for the non-driven system in section 3. In this case the general solution can be constructed and the density operator is expressed in closed form. This is, however, not fully achieved in the present analysis of the driven oscillator in section 4. We will derive a special class of solutions using similar techniques as in our previous study of a related model [5]. Nevertheless we will be able to describe some important features of the dynamics, as for example the density operator for the quantum limit cycle. Section 5 is devoted to a related description of quantum damping in terms of non-Hermitian Hamiltonians, in the present case a harmonic oscillator with a complex frequency. For convenience of the reader we begin with a brief repetition of the well-known classical dynamics. We will use units with glyph[planckover2pi1] = 1 throughout.

## 2 The classical oscillator

Let us briefly recall the damped and driven harmonic oscillator in classical mechanics. The equation of motion

$$\ddot { x } + 2 \gamma \dot { x } + \omega _ { 0 } ^ { 2 } x = \tilde { f } ( t )$$

models a harmonic oscillator with frequency ω 0 under a time-dependent force ˜ ( f t ), which will be assumed to be harmonic,

$$\tilde { f } ( t ) = \tilde { f } _ { 0 } \cos \Omega t \,,$$

in the following. The interaction with an environment is modeled by the term 2 γx ˙ where the constant γ &gt; 0 describes the friction. For a free motion with vanishing force the general solution is given by

$$x ( t ) = x _ { + } e ^ { + i \omega t - \gamma t } + x _ { - } e ^ { - i \omega t - \gamma t } \,,$$

where the constants x + and x -are determined by the initial conditions. The frequency

$$\omega = \sqrt { \omega _ { 0 } ^ { 2 } - \gamma ^ { 2 } }$$

is real valued in the underdamped case, γ &lt; ω 0 , resulting in an oscillatory motion, whereas ω is imaginary for strong damping, γ &gt; ω 0 , and both terms in (3) are exponentially decaying. In any case, the free solution finally approaches x = 0. For a harmonic force (2) one can construct a special solution, the limit cycle

$$x ^ { ( \ell c ) } ( t ) = A \cos ( \Omega t + \phi ) \,,$$

where the motion follows the external force, shifted by a phase

$$\phi = - \arctan \frac { 2 \gamma \Omega } { \omega _ { 0 } ^ { 2 } - \Omega ^ { 2 } }$$

with an amplitude

$$A = \frac { \tilde { f } _ { 0 } } { \sqrt { ( \omega _ { 0 } ^ { 2 } - \Omega ^ { 2 } ) ^ { 2 } + ( 2 \gamma \Omega ) ^ { 2 } } }.$$

The general solution is given by a sum of (3) and (5). Then for all initial conditions the solution x t ( ) converges to the limit cycle (5). As a function of the external frequency Ω, the limit cycle shows the celebrated resonance behavior, where the amplitude A shows a maximum of height A R = ˜ f / 0 (2 γ √ ω 2 0 -γ 2 ) at the resonance frequency Ω R = √ ω 2 0 -2 γ 2 for ω 2 0 &gt; γ 2 and the phase shift φ changes from φ = 0 for Ω = 0 to φ = -π for large Ω. If we identify ˙ x with the momentum p (assuming unit mass) the limit cycle appears in phase space as an ellipse, namely

$$p ^ { 2 } / \Omega ^ { 2 } + x ^ { 2 } = A ^ { 2 } \,.$$

## 3 The Lindblad master equation

The quantum dynamics of our system is described by the Lindblad master equation

$$\begin{array} { c } \cdot \\ \frac { \text{d} \hat { \rho } } { \text{d} t } = \mathcal { L } \hat { \rho } = - \text{i} \left [ \hat { H }, \hat { \rho } \right ] \\ & + \frac { \mu } { 2 } \left ( 2 \hat { a } \hat { \rho } \hat { a } ^ { \dagger } - \hat { a } ^ { \dagger } \hat { a } \hat { \rho } - \hat { \rho } \hat { a } ^ { \dagger } \hat { a } \right ) + \frac { \nu } { 2 } \left ( 2 \hat { a } ^ { \dagger } \hat { \rho } \hat { a } - \hat { a } \hat { a } ^ { \dagger } \hat { \rho } - \hat { \rho } \hat { a } \hat { a } ^ { \dagger } \right ), \end{array}$$

where the µ - and ν -terms represent (homogeneous) loss and (inhomogeneous) pumping of particles, respectively. Note that often these coefficients are written as µ = γ 0 (1 + n ) and ν = γ 0 n , where n = ν/ µ ( -ν ) is the thermal population of the environment. In the following it will be convenient to rewrite the damping coefficients as

$$\mu = \gamma ^ { \prime } + \gamma \, \quad \nu = \gamma ^ { \prime } - \gamma$$

where γ ′ and γ are denoted as diffusion and dissipation constants [6]. We will assume µ &gt; ν ≥ 0 so that γ ′ = ( µ + ) 2 and ν / γ = ( µ -ν / ) 2 are positive. A motivation and some applications of this evolution equation can be found in [4] or [7], where it is shown that this evolution equation can, for example, describe the damping of an electromagnetic field mode inside a cavity. The modes outside the cavity in a thermal reservoir cause a damping with rate γ .

It may be recalled that the density operator ˆ( ) ρ t is a positive operator with unit trace and the Lindblad evolution conserves the trace and the positivity of ˆ. ρ The expectation value of an observable ˆ is given by A

$$\langle \hat { A } \rangle _ { t } = \text{ trace } \left ( \hat { A } \hat { \rho } ( t ) \right ),$$

where we have assumed that ˆ A is not explicitly time-dependent. If ˆ H is Hermitian, we have

$$\langle \hat { A } ^ { \dagger } \rangle _ { t } = \langle \hat { A } \rangle _ { t } ^ { * } \,.$$

The Hamiltonian of the forced harmonic oscillator is

$$\hat { H } = \omega ( \hat { a } ^ { \dagger } \hat { a } + 1 / 2 ) - f ^ { * } ( t ) \, \hat { a } ^ { \dagger } - f ( t ) \, \hat { a }$$

with time-independent real ω , where ˆ a † and ˆ are the familiar bosonic crea ation and destruction operators with commutator | ˆ ˆ ] = 1. a, a † We will also use the number operator ˆ = ˆ ˆ, whose eigenvalues n a a † n = 0 , 1 , . . . are the number of oscillator quanta, sometimes also denoted as the number of photons or simply particles. This Hamiltonian models, for example, the coherent excitation of the cavity field by a monochromatic laser.

The Hamiltonian is Hermitian and can be rewritten as

$$\hat { H } = \omega ( \hat { a } ^ { \dagger } \hat { a } + 1 / 2 ) - \tilde { f } ( t ) \, \hat { x }$$

for a real valued force, using

$$\hat { x } = ( \hat { a } ^ { \dagger } + \hat { a } ) / \sqrt { 2 \omega }, \quad \hat { p } = \mathrm i ( \hat { a } ^ { \dagger } - \hat { a } ) \sqrt { \omega / 2 } \quad \text{and} \quad \tilde { f } ( t ) = \sqrt { 2 \omega } \, f ( t ) \,. \quad ( 1 5 )$$

In the following we will consider some special cases in more detail, namely a force-free system, f ( ) t = 0, which is time-independent and allows the construction of a general analytical solution ˆ( ), ρ t and a harmonic driving f ( ) = t f 0 cos Ω t with real valued amplitude f 0 . Here a closed form solution can be derived for a special set of initial distributions.

Without the forcing term, the density ˆ( ) ρ t approaches for ν &gt; 0 in the long time limit the steady state distribution [4]

$$\hat { \rho } ^ { ( \text{ss} ) } = ( 1 - \nu / \mu ) \, \text{e} ^ { \ln ( \nu / \mu ) \, \hat { n } }$$

and the mean value 〈 ˆ n 〉 is equal to n = ν/ µ ( -ν ) = ν/ γ 2 , which agrees with the mean number of photons in the environment if the system models a photon field in a cavity. For ν = 0, where we have no gain, the steady state distribution is given by the oscillator ground state ˆ ρ (ss) = 0 | 〉〈 0 . | In the following we will be interested in the Lindblad dynamics of a forced harmonic oscillator, where the force is time-dependent with emphasis on a harmonic driving.

The Lindblad equation (9) can, of course, be solved numerically (a simple Matlab program can be found in [8]), but the dynamics shows various features, which are not directly obvious and a closed form exact solution is in many cases preferable, in particular for such a quite elementary system. In the following we will discuss this in some detail. A general solution describing the density operator ˆ( ) ρ t for any given ˆ ρ 0 has been presented recently by Fujii [2] for the non-driven case and some of these results will be given below. Solutions, however, can also be found for a driven system, as for instance for the expectation values of position, momentum and energy, or for a special class of density operators as shown in the following sections.

## 3.1 The force-free oscillator

For a free evolution, f ( ) = 0, the Hamiltonian is time-independent, which t considerably simplifies the analysis. A full solution has been presented recently by Fujii [2] based on previous work [9, 10]. Here the dynamics of the density operator is mapped on a vector evolution using the Kronecker

product (see also the book by Steeb and Hardy [11] as well as [12]). 1 This derivation is heavily based on techniques from operator algebra and here we only state the final result obtained in [2], namely

$$\dots \nolimits _ { \nu } \dots \dots \dots \dots \dots \dots \dots \dots \nolimits _ { \nu } \dots \nolimits _ { \nu } \\ \hat { \rho } ( t ) \ = \ \frac { e ^ { \gamma t } } { F ( t ) } \sum _ { j = 0 } ^ { \infty } \frac { G ( t ) ^ { j } } { j! } \hat { a } ^ { \dagger j } \, e ^ { \{ - \i \omega t - \ln F ( t ) \} \hat { n } } \, \times \\ \left \{ \sum _ { k = 0 } ^ { \infty } \frac { E ( t ) ^ { k } } { k! } \hat { a } ^ { k } \hat { \rho } _ { 0 } \, \hat { a } ^ { \dagger k } \right \} \, \text{e} ^ { \{ + \i \omega t - \ln F ( t ) \} \hat { n } } \, \hat { a } ^ { j } \,, \quad \ \ ( 1 7 ) \\ \cdots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolimits _ { \nu } \dots \ nolibration$$

where the three functions E t ( ), F t ( ) and G t ( ) are solutions of some first order differential equations with F (0) = 1 and E (0) = G (0) = 0. Let us note, for example, that G t ( ) solves the Riccati equation

$$\dot { u } + ( \mu + \nu ) u - \mu u ^ { 2 } = \nu \,,$$

for an initial condition u (0) = 0. As shown in [2] these functions are given by

$$F ( t ) \ = \ \cosh \gamma t + ( \gamma ^ { \prime } / \gamma ) \sinh \gamma t$$

$$E ( t ) \ \ = \ \frac { \mu } { \gamma F ( t ) } \sinh \gamma t \, \ \ G ( t ) = \frac { \nu } { \gamma F ( t ) } \sinh \gamma t \,.$$

The general solution (17) simplifies for special cases (see [2] for details):

## ( · ) For ν = 0 one obtains

$$\hat { \rho } ( t ) = e ^ { - ( \mu / 2 + i \omega ) t \hat { n } } \Big \{ \sum _ { k = 0 } ^ { \infty } \frac { E ( t ) ^ { k } } { k! } \hat { a } ^ { k } \hat { \rho } _ { 0 } \, \hat { a } ^ { \dagger k } \Big \} e ^ { - ( \mu / 2 - i \omega ) t \hat { n } } \,.$$

glyph[negationslash]

( · ) If ν = 0 and the initial state is the oscillator ground state ˆ ρ 0 = 0 | 〉〈 0 , we | have

$$\hat { \rho } ( t ) = ( 1 - G ( t ) ) \, \text{e} ^ { \ln G ( t ) \, \hat { n } }$$

and the mean value of ˆ is given by n

$$\langle \hat { n } \rangle _ { t } = \text{trace} \left ( \hat { n } \hat { \rho } ( t ) \right ) = \frac { \nu } { 2 \gamma } \left ( 1 - \text{e} ^ { - 2 \gamma t } \right ).$$

1 Note that such a technique has also been developed as a Liouville space formulation by Ban [13].

Here the system is for all times in a thermal state, i.e. a canonical distribution. In the long time limit we clearly see with G t ( ) → ν/µ the convergence toward the steady state distribution (16).

glyph[negationslash]

( · ) For ν = 0 and a coherent initial state ˆ ρ 0 = | α 0 〉〈 α 0 | the solution is given by

$$\hat { \rho } ( t ) = e ^ { c ( t ) - \log G ( t ) \left ( \alpha ( t ) \hat { a } ^ { \dagger } + \alpha ^ { * } ( t ) \hat { a } - \hat { n } \right ) }$$

with

$$c ( t ) = | \alpha ( t ) | ^ { 2 } \log G ( t ) + \log ( 1 - G ( t ) ) \, \quad \alpha ( t ) = \alpha _ { 0 } \text{e} ^ { - ( \gamma + \imath \omega ) t } \,. \quad \ ( 2 5 )$$

(compare eq. (4.3) in [2]), or, written in product form using the disentangling relation given in eq. (56) below,

$$\hat { \rho } ( t ) = \mathrm e ^ { \tilde { c } ( t ) } \mathrm e ^ { ( 1 - G ( t ) ) \alpha ( t ) \hat { a } ^ { \dagger } } \mathrm e ^ { \log G ( t ) \hat { n } } \mathrm e ^ { ( 1 - G ( t ) ) \alpha ^ { * } ( t ) \hat { a } }$$

with

$$\tilde { c } ( t ) = \log ( 1 - G ( t ) ) - | \alpha ( t ) | ^ { 2 } ( 1 - G ( t ) ) \,.$$

( · ) Finally, for ˆ ρ 0 = | α 0 〉〈 α 0 | and ν = 0, we find from (21)

$$\hat { \rho } ( t ) \ & = \ e ^ { ( 1 - e ^ { - \mu t } ) | \alpha _ { 0 } | ^ { 2 } } \, e ^ { - ( \mu / 2 + i \omega ) t \, \hat { n } } \, | \alpha _ { 0 } \rangle \langle \alpha _ { 0 } | \, e ^ { - ( \mu / 2 - i \omega ) t \, \hat { n } } \\ & = \ | \alpha ( t ) \rangle \langle \alpha ( t ) | \quad \text{with} \quad \alpha ( t ) = \alpha _ { 0 } \, e ^ { - ( \mu / 2 + i \omega ) t }$$

(compare eq. (127) in appendix A), i.e. the system remains in a coherent state for all times and the mean values of ˆ and ˆ decay as a n

$$\langle \hat { a } \rangle _ { t } = \alpha ( t ) \, \quad \langle \hat { n } \rangle _ { t } = | \alpha ( t ) | ^ { 2 } = | \alpha _ { 0 } \, e ^ { - ( \mu / 2 + i \omega ) t } | ^ { 2 } = n _ { 0 } \, e ^ { - \mu t } \, \quad \ \quad ( 2 9 )$$

with n 0 = | α 0 | 2 .

## 3.2 Mean value evolution

Returning to the forced oscillation, we will first consider the time evolution of some expectation values. For the operators ˆ, ˆ a a † and ˆ this can be achieved n

quite easily. Let us demonstrate this in some detail for ˆ. The time derivative a of 〈 ˆ a 〉 t = trace (ˆˆ( )) yields aρ t

$$& \frac { \mathrm d \langle \hat { a } \rangle } { \mathrm d t } = \, \text{trace} \, ( \hat { a } \, \frac { \mathrm d \hat { \rho } } { \mathrm d t } ) \\ & = - \mathrm i \omega \, \text{trace} \, ( \hat { a } [ \hat { a } ^ { \dagger } \hat { a }, \hat { \rho } ] ) + \mathrm i \, \text{trace} \, ( \hat { a } \, [ f ^ { * } ( t ) \hat { a } ^ { \dagger } + f ( t ) \hat { a }, \hat { \rho } ] )$$

$$& \stackrel { \mathcal { W } } { = } - \mathrm i \omega \text{ trace } ( \hat { a } [ \hat { a } ^ { \dagger } \hat { a }, \hat { \rho } ] ) + \mathrm i \text{ trace } ( \hat { a } \, [ f ^ { * } ( t ) \hat { a } ^ { \dagger } + f ( t ) \hat { a }, \hat { \rho } ] ) \\ & + \frac { \mu } { 2 } \text{ trace } ( 2 \hat { a } \hat { a } \hat { \rho } \hat { a } ^ { \dagger } - \hat { a } \hat { a } ^ { \dagger } \hat { a } \hat { \rho } - \hat { a } \hat { \rho } \hat { a } ^ { \dagger } \hat { a } ) + \frac { \nu } { 2 } \text{ trace } ( 2 \hat { a } \hat { a } ^ { \dagger } \hat { \rho } \hat { a } - \hat { a } \hat { a } \hat { a } ^ { \dagger } \hat { \rho } - \hat { a } \hat { \rho } \hat { a } \hat { a } ^ { \dagger } ).$$

This can be simplified using trace ( ˆ ˆ ) = AB trace ( ˆ ˆ ) with the result BA

$$\frac { \mathrm d \langle \hat { a } \rangle } { \mathrm d t } = - ( \mathrm i \omega + \gamma ) \langle \hat { a } \rangle + \mathrm i f ^ { * } ( t ) \,,$$

which can be integrated as

$$\langle \hat { a } \rangle _ { t } = \text{e} ^ { - ( i \omega + \gamma ) t } \langle \hat { a } \rangle _ { 0 } + \mathcal { F } ( t )$$

with

$$\mathcal { F } ( t ) = \mathfrak { i } \int _ { 0 } ^ { t } \mathrm d t ^ { \prime } f ^ { * } ( t ^ { \prime } ) \, \mathrm e ^ { - ( \mathrm i \omega + \gamma ) ( t - t ^ { \prime } ) } \,.$$

For a time-periodic force f ( ) = t f ( t +Ω 2 ) one can evaluate the integral / π (33) by Fourier expansion

$$f ( t ) = \sum _ { k } c _ { k } \, \text{e} ^ { i k \Omega t }$$

with the result

$$\mathcal { F } ( t ) = \sum _ { k } c _ { k } ^ { * } \, \frac { \mathrm e ^ { i k \Omega t } - \mathrm e ^ { - ( i \omega + \gamma ) t } } { \omega + k \Omega - \mathrm i \gamma } \,.$$

In the same way one can derive an equation for 〈 ˆ a † 〉 t , which agrees with 〈 ˆ a 〉 ∗ t (compare eq. (12)).

For vanishing driving we have, of course,

$$\langle \hat { a } \rangle _ { t } = \text{e} ^ { - ( \text{i} \omega + \gamma ) t } \langle \hat { a } \rangle _ { 0 }$$

and for a harmonic driving f ( ) = t f 0 cos Ω t and one obtains

$$\mathcal { F } ( t ) = \frac { f _ { 0 } } { 2 } \left ( \frac { e ^ { i \Omega t } - e ^ { - ( i \omega + \gamma ) t } } { \omega + \Omega - i \gamma } + \frac { e ^ { - i \Omega t } - e ^ { - ( i \omega + \gamma ) t } } { \omega - \Omega - i \gamma } \right ),$$

which then provides the solution for the expectation values of position ˆ = x (ˆ a † + ˆ) a / √ 2 ω and momentum ˆ = i p √ ω/ 2 (ˆ a † -ˆ). a In the limit t → ∞ we have

$$\mathcal { F } ( t ) \rightarrow \mathcal { F } ^ { ( \ell c ) } ( t ) = \frac { f _ { 0 } } { 2 } \left ( \frac { e ^ { i \Omega t } } { \omega + \Omega - i \gamma } + \frac { e ^ { - i \Omega t } } { \omega - \Omega - i \gamma } \right ),$$

and from (32) we see that

$$\langle \hat { a } \rangle _ { t } \rightarrow \langle \hat { a } \rangle _ { t } ^ { ( \ell c ) } = \mathcal { F } ^ { ( \ell c ) } ( t )$$

the expectation value for the quantum limit cycle, which leads to the asymptotic expressions for 〈 x 〉 ( glyph[lscript] c) ( ) and t 〈 p 〉 ( glyph[lscript] c) ( ). t

Alternatively, for real ˜ ( f t ), one can can start from the equations of motion

$$\frac { \mathrm d \langle \hat { x } \rangle } { \mathrm d t } \ = \ \langle \hat { p } \rangle - \gamma \, \langle \hat { x } \rangle$$

$$\frac { \mathrm d \langle \hat { p } \rangle } { \mathrm d t } \ = \ - \omega ^ { 2 } \left \langle \hat { x } \right \rangle - \gamma \left \langle \hat { p } \right \rangle + \tilde { f } ( t )$$

with ˜ ( f t ) = √ 2 ωf t ( ). A second time derivative of (40) yields, after inserting (41) and eliminating 〈 ˆ , p 〉

$$\frac { d ^ { 2 } \langle \hat { x } \rangle } { d t ^ { 2 } } + 2 \gamma \frac { \mathrm d \langle \hat { x } \rangle } { \mathrm d t } + \left ( \omega ^ { 2 } + \gamma ^ { 2 } \right ) \langle \hat { x } \rangle = \tilde { f } ( t ) \,.$$

We observe, that this differential equation agrees with the classical one in (1), where the oscillator frequency ω 0 of in the classical equation is replaced by √ ω 2 + γ 2 . Comparing with section 2 we see that here only the underdamped case is realized.

For a cosine-driving, the quantum limit cycle can be copied from the classical one in (5) as

$$\langle \hat { x } \rangle _ { t } = A _ { q } \cos ( \Omega t + \phi _ { q } ),$$

where the motion follows the external force, shifted by a phase

$$\phi _ { q } = - \arctan \frac { 2 \gamma \Omega } { \omega ^ { 2 } + \gamma ^ { 2 } - \Omega ^ { 2 } }$$

with an amplitude

$$A _ { q } = \frac { \tilde { f } _ { 0 } } { \sqrt { ( \omega ^ { 2 } + \gamma ^ { 2 } - \Omega ^ { 2 } ) ^ { 2 } + ( 2 \gamma \Omega ) ^ { 2 } } } \,.$$

In the quantum system, the resonance amplitude is A R = ˜ f / 0 (2 γω ) at the resonance frequency Ω R = √ ω 2 -γ 2 . We also note, that in this limit we have

$$\langle \hat { p } \rangle ^ { ( \ell c ) } ( t ) \ & = \ \text{d} \langle \hat { x } \rangle / \text{d} t + \gamma \langle \hat { x } \rangle \\ & = \ - \Omega A _ { q } \sin ( \Omega t + \phi _ { q } ) + \gamma A _ { q } \cos ( \Omega t + \phi _ { q } ) \,,$$

and in phase space we again obtain an ellipse, namely

$$\left ( \left \langle \hat { p } \right \rangle - \gamma \left \langle \hat { x } \right \rangle \right ) ^ { 2 } / \Omega ^ { 2 } + \left \langle \hat { x } \right \rangle ^ { 2 } = A _ { q } ^ { 2 } \,.$$

Finally we will evaluate the mean value of ˆ by means of n

$$\frac { d \langle \hat { n } \rangle } { d t } = \, \text{trace} \, ( \hat { n } \, \frac { d \hat { \rho } } { d t } )$$

$$+ \frac { \mu } { 2 } \text{ trace } ( 2 \hat { n } \hat { a } \hat { \rho } \hat { a } ^ { \dagger } - \hat { n } \hat { a } ^ { \dagger } \hat { a } \hat { \rho } - \hat { n } \hat { \rho } \hat { a } ^ { \dagger } \hat { a } ) + \frac { \nu } { 2 } \text{ trace } ( 2 \hat { n } \hat { a } ^ { \dagger } \hat { \rho } \hat { a } - \hat { n } \hat { a } \hat { a } ^ { \dagger } \hat { \rho } - \hat { n } \hat { \rho } \hat { a } \hat { a } ^ { \dagger } ),$$

$$& \frac { \mathrm d \langle \hat { n } \rangle } { \mathrm d t } = \, \text{trace} \, ( \hat { n } \, \frac { \mathrm d \hat { \rho } } { \mathrm d t } ) \\ & = \, \text{trace} \, ( - \, \mathrm i \omega \hat { n } [ \hat { a } ^ { \dagger } \hat { a }, \hat { \rho } ] + \mathrm i ( \hat { n } [ f ^ { * } ( t ) \hat { a } ^ { \dagger } + f ( t ) \hat { a }, \hat { \rho } ] ) ) \\ & + \frac { \mu } { \mathrm n } \, \text{trace} \, ( 2 \hat { n } \hat { a } \hat { \rho } \hat { a } ^ { \dagger } - \hat { n } \hat { a } ^ { \dagger } \hat { a } \hat { \rho } - \hat { n } \hat { \rho } \hat { a } ^ { \dagger } ) + \frac { \nu } { \mathrm n } \, \text{trace} \, ( 2 \hat { n } \hat { a } ^ { \dagger } \hat { \rho } \hat { a } - \hat { n } \hat { a } \hat { a } ^ { \dagger } \hat { \rho } - \hat { n } \hat { \rho } \hat { a } \hat { a } ^ { \dagger } ),$$

which can be reduced to

$$\frac { \mathrm d \langle \hat { n } \rangle } { \mathrm d t } = \nu - 2 \gamma \langle \hat { n } \rangle + \mathrm i f ^ { * } ( t ) \langle \hat { a } ^ { \dagger } \rangle - \mathrm i f ( t ) \langle \hat { a } \rangle \,.$$

For the force-free oscillator the solution is

$$\lambda \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \stackrel { \sim } { \dots } \, } \, \right.$$

with 〈 n 〉 t → ν/ γ 2 for t →∞ , which agrees with eq. (23) for n 0 = 0.

For the driven system with f ( ) = t f 0 cos Ω t eq. (49) reads

$$\frac { \mathrm d \langle \hat { n } \rangle } { \mathrm d t } = \nu - 2 \gamma \langle \hat { n } \rangle + \sqrt { \frac { 2 } { \omega } } \, f _ { 0 } \cos \Omega t \, \langle \hat { p } \rangle \,,$$

and, using (46), one can quite easily see that in the long time limit we have

$$\langle \hat { n } \rangle _ { t } ^ { ( \ell c ) } \ = \ \overline { n } + \frac { \tilde { f } _ { 0 } A _ { q } } { 4 \omega } \, \cos ( 2 \Omega t + \phi _ { q } )$$

$$\overline { n } \ = \ \frac { \nu } { 2 \gamma } + \frac { \stackrel { \lambda } { \tilde { f } _ { 0 } } A _ { q } } { 4 \gamma \omega } \left ( \gamma \cos \phi _ { q } - \Omega \sin \phi _ { q } \right ),$$

which solves the differential equation (49), i.e. the mean value of ˆ oscillates n with twice the forcing frequency, again shifted by the same phase φ q , with an average value n .

## 4 A special class of solutions

The fact that the mean values 〈 ˆ n 〉 t , 〈 ˆ x 〉 t and 〈 ˆ p 〉 t can be described by a closed set of differential equations, which can be solved without the necessity to evaluate the full density matrix ˆ( ), tremendously simplifies the analysis. ρ t Note that this is no longer true for other operators as for instance 〈 ˆ a 2 〉 t . This simplification is due to the fact that the set of operators ˆ , I ˆ, a ˆ a † and ˆ ˆ a a † closes under the Lindblad evolution and the dynamics in this reduced subset simplifies. We can therefore conjecture that one can find special solutions generated solely by these operators expressing the density operator as

$$\hat { \rho } ( t ) = \text{e} ^ { z ( t ) + v ( t ) \hat { a } ^ { \dagger } \hat { a } + \delta ( t ) \hat { a } ^ { \dagger } + \delta ^ { * } ( t ) \hat { a } } \,,$$

which is Hermitian for real parameters z and v . This exponential representation can be rewritten as an exponential product

$$\hat { \rho } ( t ) = e ^ { c ( t ) } \, e ^ { \beta ( t ) \hat { a } ^ { \dagger } } e ^ { \sigma ( t ) \hat { a } ^ { \dagger } \hat { a } } e ^ { \beta ^ { * } ( t ) \hat { a } } \,,$$

which is often more convenient. The coefficients of both expressions satisfy the disentangling relations

$$v = \sigma \, \quad \delta = \frac { \beta \sigma } { e ^ { \sigma } - 1 } \, \quad z = c + | \beta | ^ { 2 } \, \frac { 1 + \sigma - e ^ { \sigma } } { ( 1 - e ^ { \sigma } ) ^ { 2 } }$$

or

$$\sigma = v \, \ \beta = \frac { \delta } { v } \left ( e ^ { v } - 1 \right ) \, \ \ c = z - \left | \, \frac { \delta } { v } \, \right | ^ { 2 } \left ( 1 + v - e ^ { v } \right )$$

(see, e.g., [2], eqs. (4.4), (4.5)).

Clearly this class of density operators is not general. It covers, however, some interesting cases. For β = 0 one obtains the canonical distribution, in the limit σ →-∞ the pure state distribution ˆ = ρ | β 〉〈 β | for a coherent state | β 〉 and, as will be shown below, it describes the limit cycle distribution of a forced oscillator.

In order to verify that such a ˆ( ) is actually a solution of the Lindblad ρ t equation we start from the well known operator identities (see, e.g., [1])

$$\mathrm e ^ { - \beta \hat { a } ^ { \dagger } } \hat { a } \, \mathrm e ^ { \beta \hat { a } ^ { \dagger } } = \hat { a } + \beta \, \quad \mathrm e ^ { - \sigma \hat { a } ^ { \dagger } \hat { a } } \hat { a } \, \mathrm e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } = \mathrm e ^ { \sigma } \hat { a } \,,$$

$$e ^ { \beta \hat { a } } \hat { a } ^ { \dagger } \, e ^ { - \beta \hat { a } } = \hat { a } ^ { \dagger } + \beta \, \quad e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \hat { a } ^ { \dagger } \, e ^ { - \sigma \hat { a } ^ { \dagger } \hat { a } } = e ^ { \sigma } \hat { a } ^ { \dagger } \,,$$

which imply

$$\hat { a } \, e ^ { \beta \hat { a } ^ { \dagger } } = e ^ { \beta \hat { a } ^ { \dagger } } ( \hat { a } + \beta ) \, \quad \ \hat { a } \, e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } = e ^ { \sigma } \, e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \hat { a } \,,$$

$$e ^ { \beta \hat { a } } \, \hat { a } ^ { \dagger } = ( \hat { a } ^ { \dagger } + \beta ) e ^ { \beta \hat { a } } \, \quad e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \hat { a } ^ { \dagger } = e ^ { \sigma } \, \hat { a } ^ { \dagger } \, e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \,.$$

A short calculation using these relations and ˆ as given in (55) yields ρ

$$\hat { a } \hat { \rho } = \mathrm e ^ { \sigma } \hat { \rho } \hat { a } + \beta \hat { \rho } \, \quad \hat { \rho } \hat { a } ^ { \dagger } = \mathrm e ^ { \sigma } \hat { a } ^ { \dagger } \hat { \rho } + \beta ^ { * } \hat { \rho } \,,$$

and the time derivative is equal to

$$\frac { \text{d} \hat { \rho } } { \text{d} t } \ & = \ \dot { c } \, \hat { \rho } + \dot { \beta } \hat { a } ^ { \dagger } \hat { \rho } + \dot { \sigma } \text{e} ^ { c } \text{e} ^ { \beta \hat { a } ^ { \dagger } } \hat { a } ^ { \dagger } \hat { a } \, \text{e} ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \text{e} ^ { \beta ^ { * } \hat { a } } + \dot { \beta } ^ { * } \hat { \rho } \hat { a } \\ & = \ \dot { c } \, \hat { \rho } + \dot { \beta } \hat { a } ^ { \dagger } \hat { \rho } + \dot { \sigma } \hat { a } ^ { \dagger } \hat { \rho } \hat { a } + \dot { \beta } ^ { * } \hat { \rho } \hat { a } \,,$$

where we have used

$$e ^ { c } e ^ { \beta \hat { a } ^ { \dagger } } \hat { a } ^ { \dagger } \hat { a } \, e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } e ^ { \beta ^ { * } \hat { a } } = e ^ { c } \hat { a } ^ { \dagger } \, e ^ { \beta \hat { a } ^ { \dagger } } e ^ { \sigma } \, e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \, e ^ { \beta ^ { * } \hat { a } } = e ^ { \sigma } \hat { a } ^ { \dagger } \hat { \rho } \, \hat { a } \,.$$

As a next step we evaluate L ˆ for the Hamiltonian (13) and reorder the terms ρ starting with the Lindblad ones using ˆˆ aa † = ˆ ˆ + 1: a a †

$$\frac { \mu } { 2 } \left ( 2 \hat { a } \hat { \rho } \hat { a } ^ { \dagger } & ^ { \dagger } \hat { a } \hat { \rho } - \hat { \rho } \hat { a } ^ { \dagger } \hat { a } \right ) + \frac { \nu } { 2 } \left ( 2 \hat { a } ^ { \dagger } \hat { \rho } \hat { a } - \hat { a } \hat { a } ^ { \dagger } \hat { \rho } - \hat { \rho } \hat { a } \hat { a } ^ { \dagger } \right ) \\ & = \mu \hat { a } \hat { \rho } \hat { a } ^ { \dagger } + \nu \hat { a } ^ { \dagger } \hat { \rho } \hat { a } - \gamma ^ { \prime } \hat { a } ^ { \dagger } \hat { a } \hat { \rho } - \gamma ^ { \prime } \hat { \rho } \hat { a } ^ { \dagger } \hat { a } - \nu \hat { \rho } \,.$$

We then have

$$\mathcal { L } \hat { \rho } \ & = \ ( - \mathrm i \omega + \gamma ^ { \prime } ) \hat { a } ^ { \dagger } \hat { a } \hat { \rho } + ( \mathrm i \omega - \gamma ^ { \prime } ) \hat { \rho } \hat { a } ^ { \dagger } \hat { a } - \nu \hat { \rho } \\ & \quad + \mathrm i f ^ { * } \hat { a } ^ { \dagger } \hat { \rho } + \mathrm i f \hat { a } \hat { \rho } - \mathrm i f ^ { * } \hat { \rho } \hat { a } ^ { \dagger } - \mathrm i f \hat { \rho } \hat { a } + \mu \hat { a } \hat { \rho } \hat { a } ^ { \dagger } + \nu \hat { a } ^ { \dagger } \hat { \rho } \hat { a } \\ & = \ ( - \mathrm i \omega + \gamma ^ { \prime } ) \hat { a } ^ { \dagger } ( \mathrm e ^ { \sigma } \hat { \rho } \hat { a } + \beta \hat { \rho } ) + ( \mathrm i \omega - \gamma ^ { \prime } ) ( \mathrm e ^ { \sigma } \hat { a } ^ { \dagger } \hat { \rho } + \beta ^ { * } \hat { \rho } ) \hat { a } - \nu \hat { \rho } \\ & \quad + \mathrm i f ^ { * } \hat { a } ^ { \dagger } \hat { \rho } + \mathrm i f ( \mathrm e ^ { \sigma } \hat { \rho } \hat { a } + \beta \hat { \rho } ) - \mathrm i f ^ { * } ( \mathrm e ^ { \sigma } \hat { a } ^ { \dagger } \hat { \rho } + \beta ^ { * } \hat { \rho } ) - \mathrm i f \hat { \rho } \hat { a } \\ & \quad + \mu \hat { a } \hat { \rho } \hat { a } ^ { \dagger } + \nu \hat { a } ^ { \dagger } \hat { \rho } \hat { a } \,,$$

and with relations (60) - (62) we find

$$\hat { a } \hat { \rho } \hat { a } ^ { \dagger } & \ = \ \ e ^ { \sigma } \hat { \rho } \hat { a } \hat { a } ^ { \dagger } + \beta \hat { \rho } \hat { a } ^ { \dagger } = e ^ { \sigma } \hat { \rho } ( \hat { a } ^ { \dagger } \hat { a } + 1 ) + \beta e ^ { \sigma } \hat { a } ^ { \dagger } \hat { \rho } + | \beta | ^ { 2 } \hat { \rho } \\ & = \ \ e ^ { 2 \sigma } \hat { a } ^ { \dagger } \hat { \rho } \hat { a } + e ^ { \sigma } \beta ^ { * } \hat { \rho } \hat { a } + \beta e ^ { \sigma } \hat { a } ^ { \dagger } \hat { \rho } + ( e ^ { \sigma } + | \beta | ^ { 2 } ) \hat { \rho } \,.$$

We finally arrive at

$$\mathcal { L } \hat { \rho } \ = \ \text{if} \, & \beta - \text{if} ^ { * } \beta ^ { * } - \nu + \mu ( \text{e} ^ { \sigma } + | \beta | ^ { 2 } ) ) \hat { \rho } \\ & + ( \text{if} ^ { * } ( 1 - \text{e} ^ { \sigma } ) + ( - \text{i} \omega - \gamma ^ { \prime } ) \beta + \mu \beta \text{e} ^ { \sigma } ) \hat { a } ^ { \dagger } \hat { \rho } \\ & + ( \text{if} ( 1 - \text{e} ^ { \sigma } ) + ( \text{i} \omega - \gamma ^ { \prime } ) \beta ^ { * } + \mu \beta ^ { * } \text{e} ^ { y } ) \hat { \rho } \hat { a } \\ & + ( \mu \text{e} ^ { 2 y } - 2 \gamma ^ { \prime } \text{e} ^ { \sigma } + \nu ) \hat { a } ^ { \dagger } \hat { \rho } \hat { a } \,,$$

and in comparison with (63) we see that the coefficients satisfy the differential equations

$$\dot { c } \ = \ \text{if} \beta - \mathrm i \beta ^ { * } - \nu + \mu ( \mathrm e ^ { \sigma } + | \beta | ^ { 2 } )$$

$$\dot { \beta } \ = \ \text{if} ^ { * } ( 1 - \text{e} ^ { \sigma } ) + ( - \text{i} \omega - \gamma ^ { \prime } ) \beta + \mu \beta \text{e} ^ { \sigma }$$

$$\dot { \sigma } \ = \ - 2 \gamma ^ { \prime } + \mu \mathrm e ^ { \sigma } + \nu \mathrm e ^ { - \sigma } \,.$$

The substitution u = e σ transforms the last equation into

$$\dot { u } = \nu - 2 \gamma ^ { \prime } u + \mu u ^ { 2 } \,,$$

which agrees with eq. (18). A solution for an initial condition 0 &lt; u 0 = u (0) = e σ (0) &lt; 1 is given by

$$u ( t ) = \frac { \mu u _ { 0 } - \nu + \nu ( 1 - u _ { 0 } ) e ^ { 2 \gamma t } } { \mu u _ { 0 } - \nu + \mu ( 1 - u _ { 0 } ) e ^ { 2 \gamma t } }$$

with u t ( ) &lt; 1. The function u t ( ) is monotonic, decreasing for u 0 &gt; ν/µ and increasing for u 0 &lt; ν/µ , with u t ( ) → ν/µ for t →∞ . For ν = 0 we have

$$u ( t ) = \frac { u _ { 0 } } { u _ { 0 } + ( 1 - u _ { 0 } ) \, e ^ { 2 \gamma t } } \,.$$

Equation (70) for β t ( ) can be simplified by means of the transformation β = α (1 -u ), which leads to

$$\dot { \alpha } = - ( \text{i} \omega + \gamma ) \alpha + \text{i} f ^ { * } ( t )$$

with initial condition α (0) = β (0) / (1 -u (0)). This agrees with (31) for 〈 ˆ a 〉 t and is solved by

$$\alpha ( t ) = \text{e} ^ { - ( i \omega + \gamma ) t } \alpha ( 0 ) + \mathcal { F } ( t )$$

with F ( ) given in (33). t

Integration of (69) yields c t ( ), i.e. the normalization factor Z t ( ) = e c t ( ) of

$$\hat { \rho } = Z \, \hat { \rho } ^ { \prime } \quad \text{with} \quad \hat { \rho } ^ { \prime } = e ^ { \beta \hat { a } ^ { \dagger } } \, e ^ { \sigma \, \hat { a } ^ { \dagger } \hat { a } } \, e ^ { \beta ^ { * } \hat { a } } \,,$$

but preferably the normalization factor can be obtained in closed form by evaluating the trace in coherent states | α 〉 :

$$\vec { \sim } \quad & \quad \cdot \quad \\ \text{trace} \, \hat { \rho } ^ { \prime } \ & = \ \int \frac { d ^ { 2 } \alpha } { \pi } \, \langle \alpha | e ^ { \beta \hat { a } ^ { \dagger } } \, e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \, e ^ { \beta ^ { * } \hat { a } } | \alpha \rangle \\ & = \ \int \frac { d ^ { 2 } \alpha } { \pi } \, e ^ { \beta \alpha ^ { * } - b ( t ) | \alpha | ^ { 2 } + \beta ^ { * } \alpha } = ( 1 / b ) \, e ^ { \beta ^ { * } \beta / b } \quad \quad \quad ( 7 8 ) \,.$$

with

$$b ( t ) = 1 - e ^ { \sigma ( t ) } = 1 - u ( t ) \,,$$

which leads to

$$Z = b \, \mathrm e ^ { - \beta ^ { * } \beta / b } \,.$$

We see that ˆ is normalizable for ρ b &gt; 0, i.e. σ &lt; 0.

It should be noted that the normalization factor Z can be employed in the same way the partition function in statistical mechanics to generate mean values by parameter differentiation. By means of

$$\frac { \partial \hat { \rho } ^ { \prime } } { \partial \beta ^ { * } } = \hat { \rho } ^ { \prime } \, \hat { a } \, \quad \frac { \partial \hat { \rho } ^ { \prime } } { \partial \beta } = \hat { a } ^ { \dagger } \, \hat { \rho } ^ { \prime }$$

one obtains

$$\langle \hat { a } \rangle = - \frac { 1 } { Z } \, \frac { \partial Z } { \partial \beta ^ { * } } = \frac { \beta } { b } \quad \text{and} \quad \langle \hat { a } ^ { \dagger } \rangle = - \frac { 1 } { Z } \, \frac { \partial Z } { \partial \beta } = \frac { \beta ^ { * } } { b }$$

and therefore

$$\langle \hat { x } \rangle = \frac { 1 } { \sqrt { 2 \omega } } \, \frac { \beta ^ { * } + \beta } { b } \ \, \ \ \langle \hat { p } \rangle = \mathrm i \sqrt { \frac { \omega } { 2 } } \, \frac { \beta ^ { * } - \beta } { b }.$$

One can also derive expressions for the expectation values of ˆ = ˆ ˆ or ˆ n a a † a † 2 and ˆ a 2 etc. using this technique. As an example, we consider ˆ = ˆ ˆ. Using n a a †

$$\frac { \partial \hat { \rho } ^ { \prime } } { \partial \sigma } \ & = \ e ^ { x \hat { a } ^ { \dagger } } \, \hat { a } ^ { \dagger } \hat { a } e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } = \hat { a } ^ { \dagger } e ^ { \beta \hat { a } ^ { \dagger } } \, \hat { a } e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } \\ & = \ \hat { a } ^ { \dagger } ( \hat { a } - \beta ) e ^ { \beta \hat { a } ^ { \dagger } } \, e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } = \hat { a } ^ { \dagger } \hat { a } \, \hat { \rho } ^ { \prime } - \beta \hat { a } ^ { \dagger } \hat { a } \, \hat { \rho } ^ { \prime }$$

one gets

$$\langle \hat { n } \rangle \ = \ - \frac { 1 } { Z } \, \frac { \partial Z } { \partial \sigma } + \beta \langle \hat { a } ^ { \dagger } \rangle = \frac { 1 } { b } - 1 + \frac { | \beta | ^ { 2 } } { b ^ { 2 } } \,.$$

We have therefore shown that the exponential product (55),

$$\hat { \rho } ( t ) = Z ( t ) \, e ^ { \beta ( t ) \hat { a } ^ { \dagger } } e ^ { \sigma ( t ) \hat { a } ^ { \dagger } \hat { a } } e ^ { \beta ^ { * } ( t ) \hat { a } } \,,$$

is form-invariant under Lindblad dynamics. Moreover we have constructed an explicit solution, which requires only the numerical solution of a single differential equation, namely (75) for α t ( ):

$$\dot { \alpha } = \text{i} f ( t ) - ( \text{i} \omega + \gamma ) \alpha$$

with α (0) = β (0) /b (0). This determines the parameter

$$\beta ( t ) = \alpha ( t ) b ( t ) \,,$$

with b t ( ) = 1 -u t ( ), where u t ( ) given in (73) does not depend on the force, and the remaining parameter in (86) is given by σ t ( ) = ln u t ( ). With the normalization Z = b e -| α | 2 b from (80) eq. (86) can be written as

$$\hat { \rho } ( t ) = b ( t ) \, \text{e} ^ { - | \beta ( t ) | ^ { 2 } / b ( t ) } \, \text{e} ^ { \beta ( t ) \hat { a } ^ { \dagger } } \text{e} ^ { \sigma ( t ) \hat { n } } \text{e} ^ { \beta ^ { * } ( t ) \hat { a } }$$

or in pure exponential form by means of the disentangling relation (56) as

$$\hat { \rho } ( t ) = b ( t ) \, e ^ { \sigma ( t ) \left ( | \alpha ( t ) | ^ { 2 } + \hat { n } - \alpha ( t ) \hat { a } ^ { \dagger } - \alpha ^ { * } ( t ) \hat { a } \right ) } \,.$$

The most important expectation values are then known analytically, as for example

$$\langle \hat { a } \rangle _ { t } = \beta ( t ) / b ( t ) = \alpha ( t )$$

given in eq. (82)).

From eq. (89) one can also derive a simple equation for the diagonal matrix elements of the density operator in coherent states,

$$\rho ( \alpha, t ) \ & = \ \langle \alpha | \hat { \rho } ( t ) | \alpha \rangle = b \, e ^ { - | \beta | ^ { 2 } / b } \langle \alpha | e ^ { \beta \hat { a } ^ { \dagger } } e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } e ^ { \beta ^ { \dagger } \hat { a } } | \alpha \rangle \\ & = \ b \, e ^ { - | \beta | ^ { 2 } / b } e ^ { \beta \hat { a } ^ { \dagger } } \langle \alpha | e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } | \alpha \rangle e ^ { \beta ^ { \ast } \alpha } \\ & = \ b \, e ^ { - | \beta | ^ { 2 } / b + \beta \alpha ^ { \ast } + b | \alpha | ^ { 2 } + \beta ^ { \ast } \alpha } = b \, e ^ { \beta \alpha ^ { \ast } + \beta ^ { \ast } \alpha },$$

where we have used

$$\langle \alpha | e ^ { \sigma \hat { a } ^ { \dagger } \hat { a } } | \alpha \rangle = \mathrm e ^ { - ( 1 - e ^ { \sigma } ) \alpha ^ { * } \alpha }$$

(see [1, eq. (3.3.24)]). ˆ( ρ α, t ) is known as the Husimi distribution, a (quasi) density in phase space ( x, p ) with x = ( α ∗ + ) α / √ 2 ω and p = i √ ω/ 2 ( α ∗ -α ). The Husimi distribution (92) can be rewritten in Gaussian form as

$$\hat { \rho } ( \alpha, t ) = b ( t ) \, \text{e} ^ { - b ( t ) \, | \alpha - \alpha ( t ) \, | ^ { 2 } },$$

where α t ( ) describes the motion of the center of the density in phase space, i.e. the mean value of ˆ. a The width is given by 1 / √ 2 ( ). b t For the special case ν = 0 we see, using u t ( ) from (74), that the width parameter is given by

$$b ( t ) = \frac { b _ { 0 } } { b _ { 0 } + ( 1 - b _ { 0 } ) \, e ^ { - 2 \gamma t } }$$

with b t ( ) → 1 for t →∞ .

Let us consider some examples in more detail. First, for a force-free system, eq. (75) is solved by

$$\alpha ( t ) = \alpha _ { 0 } \, \text{e} ^ { - ( \text{i} \omega + \gamma ) t }$$

with α 0 = α (0), and we see that the center of the Gaussian Husimi density distribution (92) traces out the path α t ( ) toward the center at α = 0, whereas the width parameter b t ( ) converges to 2 γ/µ . This limit agrees, of course, with the steady state distribution (16). For the harmonically driven system with f ( ) = t f 0 cos Ω t we can easily verify that α t ( ) is given by

$$\alpha ( t ) = \alpha _ { 0 } \, e ^ { - ( i \omega + \gamma ) t } - \frac { f _ { 0 } } { 2 } \left ( \frac { e ^ { - ( i \omega + \gamma ) t } - e ^ { i \Omega t } } { \omega + \Omega - i \gamma } + \frac { e ^ { - ( i \omega + \gamma ) t } - e ^ { - i \Omega t } } { \omega - \Omega - i \gamma } \right ) \quad ( 9 7 )$$

and, according to (85), the expectation value of ˆ is equal to n

$$\langle \hat { n } \rangle _ { t } = \frac { 1 } { b ( t ) } - 1 + | \alpha ( t ) | ^ { 2 } \,.$$

In the long time limit we have

$$\alpha ( t ) \rightarrow \alpha ^ { ( \ell c ) } ( t ) = \frac { f _ { 0 } } { 2 } \Big { ( } \frac { e ^ { i \Omega t } } { \omega + \Omega - i \gamma } + \frac { e ^ { - i \Omega t } } { \omega - \Omega - i \gamma } \Big { ) }$$

and therefore, because of b → 2 γ/µ ,

$$\langle \hat { n } \rangle _ { t } & = \frac { \nu } { 2 \gamma } + \left | \alpha ^ { ( \ell c ) } \right | ^ { 2 }$$

in agreement with the previous result (52). The density operator converges for long times to the limit cycle distribution

$$\hat { \rho } ^ { ( \ell c ) } ( t ) = \frac { 2 \gamma } { \mu } \, e ^ { | \alpha ^ { ( \ell c ) } ( t ) | ^ { 2 } } e ^ { \beta ( t ) \hat { a } ^ { \dagger } } e ^ { \log ( \nu / \mu ) \hat { n } } e ^ { \beta ^ { * } ( t ) \hat { a } }$$

with β t ( ) = 2 γ µ α ( glyph[lscript] c) ( ), where we have used t b → 2 γ/µ and σ = log(1 -b ) → log( ν/µ ). Furthermore, the phase space distribution of the limit cycle is simply

$$\rho ^ { ( \ell c ) } ( \alpha, t ) = \frac { 2 \gamma } { \mu } \, \text{e} ^ { - \frac { 2 \gamma } { \mu } | \alpha - \alpha ^ { ( \ell c ) } ( t ) | ^ { 2 } },$$

which is a periodic solution with period T = 2 π/ Ω, a cyclic steady state distribution. Note that this also implies 〈 ˆ a 〉 ( glyph[lscript] c) ( ) t = α ( glyph[lscript] c) ( ) t in agreement with the result (39). For the special case ν = 0 the limit cycle distribution is a coherent state distribution

$$\rho ^ { ( \ell c ) } ( \alpha, t ) = \, e ^ { - | \alpha - \alpha ^ { ( \ell c ) } ( t ) | ^ { 2 } } \,.$$

Figure 1 illustrates the long time dynamics for a system with both loss and gain, µ = 0 6 and . ν = 0 4, for . ω = 1 1, and resonant harmonic driving . Ω = Ω R = √ ω 2 0 -γ 2 = 1 0954 with an amplitude . f 0 = 1 4. . The Husimi phase space distribution ˆ ρ ( glyph[lscript] c) ( ) for the limit cycle given in eq. (102) is shown t in the ( x, p )-plane ( α = ( ωx + i p / ) √ 2 ω ) for six equidistant times in the period 2 π/ Ω. The center of the distribution follows the ellipse (47) where 〈 ˆ x 〉 ( glyph[lscript] c) ( ) and t 〈 ˆ p 〉 ( glyph[lscript] c) ( ) are given in (43) and (46). t

Let us finally remark that for ν = 0 the limit cycle distribution turns out to be a pure coherent state, which implies that any initial distribution is driven into a coherent state. However we know that a coherent state remains coherent under the evolution, so that its limiting distribution must be necessarily coherent, with the unavoidable consequence that the limit state must be coherent.

Closing this section, we briefly look at the entropy

$$S ( t ) = - \, \text{trace} \left ( \hat { \rho } ( t ) \log \hat { \rho } ( t ) \right ),$$

Figure 1: Husimi phase space distribution ˆ ρ ( glyph[lscript] c) ( α, t ) for the limit cycle given in eq. (102) shown for six equidistant times. The center of the distribution follows the ellipse (47) shown as a red curve.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01187v2-figure-c5d943da-b3eb-412c-8067-29a77d386b08.jpg)

which can be most easily evaluated by expressing the density operator in pure exponential form ˆ = ρ b e σ ( | α | 2 +ˆ n -αa ˆ † -α a ∗ ˆ) as in (90). We then have

$$S ( t ) \ & = \ - \text{trace} \left ( \hat { \rho } \left ( \log b + \sigma ( | \alpha | ^ { 2 } + \hat { n } - \alpha \hat { a } ^ { \dagger } - \alpha ^ { * } \hat { a } ) \right ) \right ) \\ & = \ - \log b - \sigma ( | \alpha | ^ { 2 } + \langle \hat { n } \rangle - \alpha \langle \hat { a } ^ { \dagger } \rangle - \alpha ^ { * } \langle \hat { a } \rangle ) \,.$$

Inserting the expectation values 〈 ˆ a 〉 = α , 〈 ˆ a † 〉 = α ∗ (see eq. (91)) and 〈 ˆ n 〉 = 1 /b -1 + | α | 2 from eqs. (82) and (85), one arrives, with b = 1 -e σ and σ = log u , at

$$S ( t ) = - \log ( 1 - u ( t ) ) - \frac { u ( t ) } { 1 - u ( t ) } \, \log u ( t ) \,,$$

where u t ( ) is given in closed form in eq. (73). Because of d S/ u &gt; d 0 and the monotonicity of u t ( ), the entropy is monotonically decreasing for u 0 &gt; ν/µ and increasing for u 0 &lt; ν/µ . It should be pointed out, that the parameter u t ( ), and hence the entropy S t ( ), does only depend on the Lindblad parameters µ and ν and not at all on the driving force. In the long time limit u t ( ) approaches ν/µ and the resulting value S ∞ is

$$S _ { \infty } = - \frac { 1 } { 2 \gamma } \left ( \nu \log \nu - \mu \log \mu + 2 \gamma \log 2 \gamma \right ) \,.$$

In particular we find S ∞ = 0 for ν = 0, in accordance with the fact that the limit cycle is a coherent state. Let us also note that the entropy of the (time-dependent) limit cycle distribution ˆ( ) is ρ t constant in time, which is, of course, obvious in view of the fact of the dynamics of the phase density discussed above. More results concerning the concept of entropy in context with the force-free damped harmonic oscillator can be found in [14].

## 5 Non-Hermitian Hamiltonian dynamics

In the preceding sections the quantum damped harmonic oscillator has been described by a Lindblad master equation. It is, however, also possible to account for damping in terms of non-Hermitian Hamiltonians. As an example, the Hamiltonian

$$\hat { H } = \tilde { \omega } ( \hat { a } ^ { \dagger } \hat { a } + 1 / 2 ) - f ( t ) \left ( \hat { a } ^ { \dagger } + \hat { a } \right )$$

with a complex-valued frequency

$$\tilde { \omega } = \omega - \mathrm i \gamma \ \, \ \ \omega, \gamma > 0$$

has been considered in [15] in context with an analysis of a semiclassical limit of non-Hermitian quantum dynamics. Such a non-Hermitian system can be used to describe a pure loss of particles and in the following we will compare non-Hermitian and Lindblad dynamics, in particular for the case of non gain ( ν = 0).

Let us first consider the simple force-free case f ( ) = 0. t Then the solution of the Schr¨dinger equation i o d d t | ψ 〉 = ˆ H ψ | 〉 is

$$| \psi ( t ) \rangle = e ^ { - i \hat { H } t } | \psi _ { 0 } \rangle = e ^ { - i \tilde { \omega } ( \hat { a } ^ { t } \hat { a } + 1 / 2 ) t } | \psi _ { 0 } \rangle \,.$$

The non-Hermitian time evolution does not conserve the norm, which decays as

$$\langle \psi ( t ) | \psi ( t ) \rangle = \langle \psi _ { 0 } | e ^ { i \bar { \omega } ^ { * } ( \hat { a } ^ { \dagger } \hat { a } + \frac { 1 } { 2 } ) t } e ^ { - i \bar { \omega } ( \hat { a } ^ { \dagger } \hat { a } + \frac { 1 } { 2 } ) t } | \psi _ { 0 } \rangle = \langle \psi _ { 0 } | e ^ { - 2 \gamma ( \hat { a } ^ { \dagger } \hat { a } + \frac { 1 } { 2 } ) t } | \psi _ { 0 } \rangle \,. \ ( 1 1 1 )$$

Using (58) we have

$$e ^ { i \tilde { \omega } ^ { * } ( \hat { a } ^ { \dagger } \hat { a } + \frac { 1 } { 2 } ) t } \hat { a } e ^ { - i \tilde { \omega } ( \hat { a } ^ { \dagger } \hat { a } + \frac { 1 } { 2 } ) t } = e ^ { - \gamma t } e ^ { i \tilde { \omega } ^ { * } \hat { n } t } \hat { a } e ^ { - i \tilde { \omega } \hat { a } ^ { \dagger } \hat { a } t } \\ = e ^ { - \gamma t - i \tilde { \omega } t } e ^ { i \tilde { \omega } ^ { * } \hat { a } ^ { \dagger } \hat { a } t } e ^ { - i \tilde { \omega } \hat { a } ^ { \dagger } \hat { a } t } \hat { a } = e ^ { - \gamma t - i \tilde { \omega } t } e ^ { - 2 \gamma \hat { a } ^ { \dagger } \hat { a } t } \hat { a } \,, \quad ( 1 1 2 )$$

and for a coherent initial state | α 0 〉 the expectation value of ˆ is given by a

$$\langle \hat { a } \rangle _ { t } = \frac { \langle \psi ( t ) | \hat { a } | \psi ( t ) \rangle } { \langle \psi ( t ) | \psi ( t ) \rangle } = \alpha _ { 0 } \, \text{e} ^ { - \mathrm i \tilde { \omega } t } = \alpha ( t ) \,,$$

and the expectation value of ˆ is n

$$\langle \hat { n } \rangle _ { t } = | \alpha ( t ) | ^ { 2 } = n _ { 0 } \, \text{e} ^ { - 2 \gamma t }.$$

with n 0 = | α 0 | 2 . These results are in perfect agreement with the corresponding Lindblad expectation values in eq. (29) for ν = 0.

For the time-dependent driven system, the analysis is somewhat more involved. A closed form solution of the Schr¨dinger equation can be found o in [15] and here we will simply copy those results. For a coherent initial state | α 0 〉 , the solution is given by

$$| \psi ( t ) \rangle = \mathrm e ^ { - \mathrm i \tilde { \omega } t / 2 + A ( t ) + B ( t ) \alpha _ { 0 } - | \alpha _ { 0 } | ^ { 2 } / 2 + | \alpha ( t ) | ^ { 2 } / 2 } \, | \alpha ( t ) \rangle,$$

where the parameter of the coherent state | α t ( ) 〉 is given by

$$\alpha ( t ) = C ( t ) + \alpha _ { 0 } \, \mathrm e ^ { - i \tilde { \omega } t }$$

and the coefficients A t ( ), B t ( ) and C t ( ) can be found in appendix B. These terms are, however, only important if one needs to calculate explicitly the norm 〈 ψ t ( ) | ψ t ( ) 〉 . If one is only interested in the expectation values

$$\langle \hat { A } \rangle _ { t } = \frac { \langle \psi ( t ) | \hat { A } | \psi ( t ) \rangle } { \langle \psi ( t ) | \psi ( t ) \rangle }$$

most of these terms cancel, and we obtain the results

$$\langle \hat { a } \rangle _ { t } = \alpha ( t ) = C ( t ) + \alpha _ { 0 } \, e ^ { - i \tilde { \omega } t } \, \quad \langle \hat { n } \rangle _ { t } = | \alpha ( t ) | ^ { 2 },$$

which reduce, of course, for f = 0 to the expressions (113) and (114). These findings are immediately clear if one considers the renormalized density operator

$$\hat { \rho } ( t ) = \frac { | \psi ( t ) \rangle \langle \psi ( t ) | } { \langle \psi ( t ) | \psi ( t ) \rangle } = | \alpha ( t ) \rangle \langle \alpha ( t ) | \,.$$

The Husimi phase space distribution can be also calculated easily from the wave function (115):

$$\hat { \rho } ( \alpha, t ) \ = \ \frac { | \langle \alpha | \psi ( t ) \rangle | ^ { 2 } } { \langle \psi ( t ) | \psi ( t ) \rangle } = | \langle \alpha | \alpha ( t ) \rangle | ^ { 2 } = e ^ { - | \alpha - \alpha ( t ) | ^ { 2 } }. \quad \ \ ( 1 2 0 )$$

If these results are compared with the Lindblad evolution, we find precise agreement if the Lindblad parameter ν is chosen as ν = 0 and therefore γ = µ . It should be stressed, however, that this agreement is only valid for an initial coherent state, which remains coherent under the Lindblad evolution and, after renormalization, also in the non-Hermitian description.

## 6 Concluding remarks

In conclusion, we have described the dynamics of an elementary dynamical system, a damped and driven harmonic oscillator. This model is almost unavoidable in teaching physics, however only rarely described in quantum mechanics. Here we have presented a quite simple analysis using algebraic techniques, which allow an analytical solution of the Lindblad master equation and a calculation of important expectation values in closed form. In particular we have derived an analytic expression for the asymptotic density operator for harmonic driving, the limit cycle. Remarkably, for ν = 0 this is a coherent state moving along the classical limit cycle, a pure state.

The present analysis covers, however, only a special class of solutions and the derivation of a general solution remains a challenging problem for future studies as well as a deeper study of the interrelation between the Lindblad description and non-Hermitian Hamiltonians.

## Appendix

## A Coherent initial states

Some remarks concerning coherent state distributions may be helpful. As noted in the beginning of section 4, the exponential product

$$\hat { \rho } = Z \, e ^ { \beta \hat { a } ^ { \dagger } } \mathrm e ^ { \sigma \hat { n } } e ^ { \beta ^ { * } \hat { a } } \,,$$

with real σ and Z = b e -| β | 2 /b with b = 1 -e σ describes in the limit σ →-∞ the pure state distribution | β 〉〈 β | of a coherent state | β 〉 . This is can be easily understood because of

$$\lim _ { \sigma \to - \infty } e ^ { \sigma \hat { n } } = \lim _ { \sigma \to - \infty } \sum _ { n = 0 } ^ { \infty } e ^ { \sigma n } | n \rangle \langle n | = | 0 \rangle \langle 0 |$$

and therefore with Z → e -| β | 2 /b for σ →-∞

$$\lim _ { \sigma \to - \infty } \hat { \rho } = e ^ { - | \beta | ^ { 2 } } e ^ { \beta \hat { a } ^ { \dagger } } | 0 \rangle \langle 0 | e ^ { \beta ^ { * } \hat { a } } = | \beta \rangle \langle \beta |$$

(for the last equality see, e.g., [1, eq. (3.3.42)]).

An example of such a density operator is found in eq. (24) for the forcefree oscillator, namely

$$\hat { \rho } ( t ) = e ^ { c ( t ) } \, \text{e} ^ { - \log G ( t ) \left ( \beta ( t ) \hat { a } ^ { t } + \beta ^ { * } ( t ) \hat { a } - \hat { n } \right ) }$$

with

$$c ( t ) = | \alpha _ { 0 } | ^ { 2 } e ^ { - 2 \gamma t } \log G ( t ) + \log ( 1 - G ( t ) ) \, \quad \beta ( t ) = \alpha _ { 0 } e ^ { - ( \gamma + 1 \omega ) t } \,, \quad ( 1 2 5 )$$

where G t ( ) is given in (20), and ˆ(0) = ρ | α 0 〉〈 α 0 | . By means of the disentangling relations (56) this can be rewritten as an exponential product

$$\hat { \rho } ( t ) = \mathrm e ^ { c ^ { \prime } ( t ) } \, \mathrm e ^ { \beta ( t ) ( 1 - G ( t ) ) \hat { a } ^ { \dagger } } \, \mathrm e ^ { \log G ( t ) \, \hat { n } } \, \mathrm e ^ { \beta ^ { ^ { * } } ( t ) ( 1 - G ( t ) ) \hat { a } }$$

with c ′ ( ) = t -| β t ( ) | 2 (1 -G t ( )) + log(1 -G t ( )). Initially we have G (0) = 0 and β (0) = α 0 and from the arguments above we see ˆ(0) = ρ | α 0 〉〈 α 0 | . For ν = 0, where we have G t ( ) = 0 for all times, we find in the same way

$$\hat { \rho } ( t ) = | \beta ( t ) \rangle \langle \beta ( t ) |$$

with β t ( ) = α 0 e -( µ/ 2+i ω t ) , i.e. the system stays in a coherent state for all times.

## B Details of non-Hermitian dynamics

For convenience of the reader we reproduce here some details of the nonHermitian dynamics from [15]. These results can be derived, e.g., by epressing the time-evolution operator as an exponential product as in section 4 or by extending the analysis of the Hermitian forced harmonic oscillator in [1, sect. 3.11] to complex frequencies. For a coherent initial state | α 0 〉 this leads to the the wave function (115), where the coefficients A t ( ), B t ( ) and C t ( ) are solutions of the differential equations

$$i \dot { B } = - f \mathrm e ^ { - i \tilde { \omega } t } \, \ \ i \dot { C } = \tilde { \omega } C - f ^ { * } \, \ \ i \dot { A } = - f C$$

with initial conditions A (0) = B (0) = C (0) = 0, namely

$$B ( t ) & = - \mathrm i \int _ { 0 } ^ { t } \mathrm d t ^ { \prime } \mathrm e ^ { \mathrm i \tilde { \omega } t ^ { \prime } } f ( t ^ { \prime } ) \, \quad C ( t ) = - \mathrm i \int _ { 0 } ^ { t } \mathrm d t ^ { \prime } \mathrm e ^ { \mathrm i \tilde { \omega } ( t - t ) } f ^ { * } ( t ^ { \prime } ) \\ A ( t ) & = \mathrm i \int _ { 0 } ^ { t } \mathrm d t ^ { \prime } \, f ( t ^ { \prime } ) C ( t ^ { \prime } ). \\ \Big. &$$

For harmonic driving force f ( ) = t ˜ cos Ω f 0 t these solutions are given by

$$B ( t ) = \frac { - f _ { 0 } } { 2 ( \tilde { \omega } ^ { 2 } - \Omega ^ { 2 } ) } \Big \{ ( \tilde { \omega } + \Omega ) e ^ { - i ( \tilde { \omega } - \Omega ) t } + ( \tilde { \omega } - \Omega ) e ^ { - i ( \tilde { \omega } + \Omega ) t } - 2 \tilde { \omega } \Big \} ( 1 3 0 )$$

$$C ( t ) = \frac { \cdot } { 2 ( \tilde { \omega } ^ { 2 } - \Omega ^ { 2 } ) } \Big \{ ( \tilde { \omega } - \Omega ) e ^ { i \Omega t } + ( \tilde { \omega } + \Omega ) e ^ { - i \Omega t } - 2 \tilde { \omega } e ^ { - i \tilde { \omega } t } \Big \} \quad ( 1 3 1 )$$

$$A ( t ) & = \frac { \tau ( \omega ^ { - } - \mathfrak { s } t ^ { - } ) \, \mathbb { C } } { 4 ( \tilde { \omega } ^ { 2 } - \Omega ^ { 2 } ) } \Big \{ 1 + 2 \mathfrak { i } \tilde { \omega } t + \frac { \tilde { \omega } + \Omega } { 2 \Omega } \mathrm e ^ { - i ( \tilde { \omega } - \Omega ) t } - \frac { \tilde { \omega } - \Omega } { 2 \Omega } \mathrm e ^ { - i ( \tilde { \omega } + \Omega ) t } \\ & \quad + \frac { 2 \tilde { \omega } } { \tilde { \omega } ^ { 2 } - \Omega ^ { 2 } } \Big ( ( \tilde { \omega } + \Omega ) \mathrm e ^ { - i ( \tilde { \omega } - \Omega ) t } + ( \tilde { \omega } - \Omega ) \mathrm e ^ { - i ( \tilde { \omega } + \Omega ) t } - 2 \tilde { \omega } \Big ) \Big \}. ( 1 3 2 ) \\ & \quad \,. \.$$

(Note that the factor 2 in the denominator in eq. (130) is missing in [15].)

## References

- [1] W. H. Louisell, Quantum Statistical Properties of Radiation , John Wiley, New York, 1973
- [2] K. Fujii, arXiv:1209.1437 (2012)
- [3] A. Lampo, S. H. Lim, J. Wehr, P. Massignan, and M. Lewenstein, Phys. Rev. A 94 (2016) 042123
- [4] B.-G. Englert and G. Morigi, arXiv:quant-ph/0206116 (2012)
- [5] M. Hensel and H. J. Korsch, J. Phys. A 25 (1992) 2043
- [6] J. Piilo and S. Maniscalco, Phys. Rev. A 74 (2006) 032303
- [7] H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems , Oxford University Press, Oxford, 2002
- [8] H. J. Korsch and K. Rapedius, Europ. J. Phys. 37 (2016) 055410

- [9] R. Endo, K. Fujii, and T. Suzuki, arXiv:0710.2724 (2012)
- [10] K. Fujii and T. Suzuki, arXiv:0806.2169 (2012)
- [11] W.-H. Steeb and Y. Hardy, Matrix Calculus and Kronecker Product , World Scientific, Singapore, 2011
- [12] M. Am-Shallem, A. Levy, I. Schaefer, and R. Kosloff, arXiv:1510.08634v2 (2015)
- [13] M. Ban, J. Math. Phys. 33 (1992) 3213
- [14] A. Isar, Fortschr. Phys. 47 (1999) 855
- [15] E. M. Graefe, M. H¨ning, and H. J. Korsch, J. Phys. A o 43 (2010) 075306
H. J. Korsch
Lindblad dynamics of the damped and forced quantum harmonic oscillator
The quantum dynamics of a damped and forced harmonic oscillator is
investigated in terms of a Lindblad master equation. Elementary algebraic
techniques are employed allowing for example to analyze the long time behavior,
i.e. the quantum limit cycle. The time evolution of various expectation values
is obtained in closed form as well as the entropy and the Husimi phase space
distribution. We also discuss the related description in terms of a
non-Hermitian Hamiltonian.