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 its dynamics analytically by employing tools from symplectic geometry. In particular, we use this result to completely characterize the dynamics of the two oscillators interacting in the ultrastrong-coupling regime with additional single-mode squeezing on both oscillators, as...