Thermal state of the general time-dependent harmonic oscillator

Taking advantage of dynamical invariant operator, we derived quantum mechanical solution of general time-dependent harmonic oscillator. The uncertainty relation of the system is always larger than ħ/2 not only in number but also in the thermal state as expected. We used the diagonal elements of density operator satisfying Leouville-von Neumann equation to calculate various expectation values in the thermal state. We applied our theory to a special case which is the forced Caldirola-Kanai oscillator.