The stellar decomposition of Gaussian quantum states
Abstract
We introduce the stellar decomposition, a novel method for characterizing non-Gaussian states produced by photon-counting measurements on Gaussian states. Given an (m+n)-mode Gaussian state G, we express it as an (m+n)-mode "Gaussian core state" G_core followed by a fixed m-mode Gaussian transformation T that only acts on the first m modes. The defining property of the Gaussian core state G_core is that measuring the last n of its modes in the photon-number basis leaves the first m modes on a finite Fock support, i.e. a core state. Since T is measurement-independent and G_core has an exact and finite Fock representation, this decomposition exactly describes all non-Gaussian states obtainable by projecting n modes of G onto the Fock basis. For pure states we prove that a physical pair (G_core, T) always exists with G_core pure and T unitary. For mixed states, we establish necessary and sufficient conditions for (G_core, T) to be a Gaussian mixed state and a Gaussian channel. Finally, we develop a semidefinite program to extract the "largest" possible Gaussian channel when these conditions fail. The stellar decomposition leads to practical bounds on achievable state quality in photonic circuits and for GKP state generation in particular. Our results are based on a new characterization of Gaussian completely positive maps in the Bargmann picture, which may be of independent interest. As a result, this work provides novel tools for improved simulations of quantum optical systems, and for understanding the generation of non-Gaussian resource states.