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" Gcore followed by an m-mode Gaussian transformation T that only acts on the first m modes. The defining property of the Gaussian core state Gcore 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 Gcore 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 (Gcore, T) always exists with Gcore pure and T unitary. For mixed states, we establish necessary and sufficient conditions for (Gcore, T) to be a Gaussian mixed state and a Gaussian channel. We also develop a semidefinite program to extract the "largest" possible Gaussian channel when these conditions fail. Finally, we present a formal stellar decomposition for generic operators, which is useful in simulations where the only requirement is that the two parts contract back to the original operator. 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.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.