General solution of the time evolution of two interacting harmonic oscillators.

We study the time evolution of an ideal system composed of two harmonic oscillators coupled through a quadratic Hamiltonian with arbitrary interaction strength. We solve the dynamics analytically by employing Lie algebraic tools that allow to decouple the time-evolution operator induced by quadratic Hamiltonians. In particular, we use this result to completely chracterize the dynamics of the two oscillators interacting in the ultrastrong coupling regime. Furthermore, we compute quantities of interest, such as the average number of excitations and the correlations that are established between the two subsystems due to the evolution. We also provide an exact decoupling of the time evolution in terms of simple quantum optical operations, which can be used for practical implementations and studies. Finally, we show how our techniques can be extended to include more oscillators and higher order interactions.