There are several definitions of a (bosonic) gaussian state (related question).
One of them is given in the paper "Gaussian quantum channels" by J. Eisert and M.M. Wolf on page 3, equation 4. Here it is:
Gaussian states are exactly those having a Gaussian characteristic function $$\chi_\rho(\xi) = e^{-\xi^T\Gamma\xi/4 + D^T\xi}.$$ Here, the $2n \times 2n$-matrix $\Gamma$ and the vector $D \in \mathbb{R}^{2n}$ are essentially the first and second moments
The definition of the characteristic function is given in the same paper - equation 3.
But I also found some slides of Dr. Alessio Serafini for ICPT Winter College on Optics. Here on page 4 is the following definition
Let us define the set of Gaussian states as all the ground and thermal states of quadratic Hamiltonians with positive definite Hamiltonian matrix $H > 0$.
Are they really equivalent? How to prove it?