Gaussian quantum operator representation for bosons

We introduce the Gaussian quantum operator representation, using the most general multi-mode Gaussian operator basis. The representation unifies and substantially extends existing phase-space representations of density matrices for Bose systems, and also includes generalized squeezed-state and thermal bases. It enables first-principles dynamical or equilibrium calculations in quantum many-body systems, with quantum uncertainties appearing as dynamical objects. Any quadratic Liouville equation for the density operator results in a purely deterministic time evolution. Any cubic or quartic master equation can be treated using stochastic methods.