Here n ¯ ¯ is the mean number of excitations in the reservoir damping the oscillator and γ is the decay rate. If we also add additional unitary evolution generated by the quantum harmonic oscillator Hamiltonian with frequency ω ω c , we obtain