Simulating realistic non-Gaussian state preparation.

We consider conditional photonic non-Gaussian state preparation using multimode Gaussian states and photon-number-resolving detectors in the presence of photon loss. While simulation of such state preparation is often computationally challenging, we show that obtaining the required multimode Gaussian state Fock matrix elements can be reduced to the computation of matrix functions known as loop hafnians, and develop a tailored algorithm for their calculation that is faster than previously known methods. As an example of its utility, we use our algorithm to explore the loss parameter space for three specific non-Gaussian state preparation schemes: Fock state heralding, cat state heralding, and weak cubic-phase state heralding. We confirm that these schemes are fragile with respect to photon loss, yet find that there are regions in the loss parameter space that are potentially accessible in an experimental setting which correspond to heralded states with non-zero non-Gaussianity.