Exploiting higher-order correlation functions for photon-statistics-based characterization and reconstruction of arbitrary Gaussian states

Gaussian states are an essential building block for various applications in quantum optics and quantum information science, yet the precise relation between their second- and third-order correlation functions remains not fully explored. We discuss connections between these correlation functions by constructing an explicit decomposition formula for arbitrary sixth-order moments of ladder operators for general Gaussian states and demonstrate how the derived relations enable state classification from correlation data alone. Whereas violating these relations certifies non-Gaussianity, satisfying them provides evidence for a Gaussian-state description and allows a direct distinction among non-displaced, non-squeezed, and displaced-squeezed sectors of the Gaussian state space. Further, we show that it is not possible to uniquely extract state parameters solely from correlation-function measurements without prior assumptions about the Gaussian state. Resolving this ambiguity requires additional loss-sensitive information, e.g., measuring the mean intensity or the vacuum overlap of each mode. In particular, we show under which circumstances these measurements can be used to reconstruct a generic Gaussian state.