Classical simulation and quantum resource theory of non-Gaussian optics
Oliver Hahn, Ryuji Takagi, Giulia Ferrini, Hayata Yamasaki
TL;DR
This work addresses the classical simulation of continuous-variable systems with non-Gaussian inputs under Gaussian operations by decomposing non-Gaussian states into Gaussian superpositions and tracking relative phases with an extended covariance formalism. It introduces two simulators—an exact one costing $\mathcal{O}(\chi^2)$ and an approximate, sparsified one with linear scaling in $\|\mathbf{c}\|_1$—and defines the Gaussian rank $\chi$ and Gaussian extent $\xi$ as operational measures of non-Gaussianity. The authors connect these measures to a quantum-resource framework via robustness, prove monotonicity under Gaussian operations, and derive optimal decompositions for states relevant to CV quantum computing, including bounds for Gaussian boson sampling and cat-state breeding. The results provide both practical tools for classically bounding CV quantum advantages and fundamental insights into the resource cost of non-Gaussianity in continuous-variable computation.
Abstract
We propose efficient algorithms for classically simulating Gaussian unitaries and measurements applied to non-Gaussian initial states. The constructions are based on decomposing the non-Gaussian states into linear combinations of Gaussian states. We use an extension of the covariance matrix formalism to efficiently track relative phases in the superpositions of Gaussian states. We get an exact simulation algorithm, which costs quadratically with the number of Gaussian states required to represent the initial state, and an approximate simulation algorithm, which costs linearly with the $l_1$ norm of the coefficients associated with the superposition. We define measures of non-Gaussianity quantifying this simulation cost, which we call the Gaussian rank and the Gaussian extent. From the perspective of quantum resource theories, we investigate the properties of this type of non-Gaussianity measure and compute optimal decompositions for states relevant to continuous-variable quantum computing.
