“…The Hamiltonian in Eq. (1) is quadratic, so it can be solved using many different methods; such as Heisenberg equations of motion, 9,12,[22][23][24]26,27,29,31,40,47,49 Feynman-Vernon path integrals, 37 extended quantum Langevin equations, 38 Green's functions, 8 and Keldysh Green's functions. 27,34,36,39,50 Here we use the Heisenberg equations of motion, which consist of a set of linear first-order differential equation, 9,12,[22][23][24]26,27,29,31,40,47,49 solved using a Laplace transform.…”