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 nonzero non-Gaussianity.
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 nonzero non-Gaussianity.
We present a detailed analytic framework for studying multimode non-Gaussian states that are conditionally generated when few modes of a multimode Gaussian state are subject to photon-number-resolving detectors. From the output state Wigner function, we deduce that the state factorizes into a Gaussian gate applied to a finite Fock-superposition non-Gaussian state. The framework provides an approach to find the optimal strategy to generate a given target non-Gaussian state. We explore examples, such as the generation of cat states, weak cubic phase states, and bosonic code states, and achieve improvements of success probability over other schemes. Our framework also applies to the case in which the measured Gaussian state is mixed which is very important for the analysis of experimental imperfections such as photon loss. The framework has potential far-reaching implications to the generation of bosonic error-correcting codes and for the implementation of non-Gaussian gates using resource states, among other applications requiring non-Gaussianity
Generation of highly non-classical quantum states of light is essential for optical quantum information processing and quantum metrology. Given the lack of sufficiently strong nonlinear interactions between optical fields, the commonly employed optical quantum-state preparation schemes are conditional, based on nonlinearity induced by heralding photon number measurement on a part of a multimode squeezed Gaussian state. Development and optimization of such probabilistic quantum-state engineering schemes represents one of the central challenges in current quantum optics. As technology advances and experiments progress to detection of higher numbers of photons, the maximization of the heralding probability becomes essential to ensure sufficiently high state-preparation rates. Here, we show that for the conditional quantum state preparation schemes based on Gaussian states and photon number measurements the maximization of the heralding probability can be formulated as finding solution to a system of polynomial equations, which offers an efficient way to find the optimal configuration and allows us to apply techniques dedicated specifically to solving such systems of equations. Our approach can seamlessly incorporate bounds on the available single-mode quadrature squeezing, which is highly experimentally relevant. We mainly consider generation of finite superpositions of Fock states but show that the approach can be straightforwardly extended to generation of squeezed superpositions of Fock states. We focus on Gaussian states with vanishing coherent displacements, hence the conditionally generated states have well defined photon number parity. We illustrate our general methodology through examples of generation of single-mode states with two and three heralding modes, and generation of two-mode states with two heralding modes.
Generation of high-fidelity photonic non-Gaussian states is a crucial ingredient for universal quantum computation using continuous-variable platforms, yet it remains a challenge to do this efficiently. We present a general framework for a probabilistic production of multimode non-Gaussian states by measuring a few modes of multimode Gaussian states via photon-number-resolving detectors. We use Gaussian elements consisting of squeezed displaced vacuum states and interferometers, the only non-Gaussian elements consisting of photon-number-resolving detectors. We derive analytic expressions for the output Wigner function, and the probability of generating the states in terms of the mean and the covariance matrix of the Gaussian state and the photon detection pattern. We find that the output states can be written as a Fock-basis superposition state followed by a Gaussian gate, and we derive explicit expressions for these parameters. These analytic expressions show exactly what non-Gaussian states can be generated by this probabilistic scheme. Further, it provides a method to search for the Gaussian circuit and measurement pattern that produce a target non-Gaussian state with optimal fidelity and success probability. We present specific examples such as the generation of cat states, ON states, Gottesman-Kitaev-Preskill states, NOON states, and bosonic-code states. The proposed framework has potentially far-reaching implications for the generation of bosonic error-correction codes that require non-Gaussian states and resource states for the implementation of non-Gaussian gates needed for universal quantum computation, among other applications requiring non-Gaussianity. The tools developed here could also prove useful for the quantum resource theory of non-Gaussianity.
To understand quantum optics experiments, we must perform calculations that consider the principal sources of noise, such as losses, spectral impurity, and partial distinguishability. In both discrete and continuous variable systems, these can be modeled as mixed Gaussian states over multiple modes. The modes are not all resolved by photon-number measurements and so require calculations on coarse-grained photon-number distributions. Existing methods can lead to a combinatorial explosion in the time complexity, making this task unfeasible for even moderately sized experiments of interest. We prove that the computation of this type of distribution can be done in exponential time, providing a combinatorial speed-up. We develop numerical techniques that allow us to determine coarse-grained photon number distributions of Gaussian states, as well as density matrix elements of heralded non-Gaussian states prepared in the presence of spectral impurity and partial distinguishability. These results offer significant speed-up and accuracy improvements to validation tests of both Fock and Gaussian boson samplers that rely on binned probability distributions. Moreover, they pave the way to a more efficient simulation of realistic photonic circuits, unlocking the ability to perform exact calculations at scales that were previously out of reach. In addition, our results, including loop Hafnian master theorems, may be of interest to the fields of combinatorics and graph theory.
To understand quantum optics experiments, we must perform calculations that consider the principal sources of noise, such as losses, spectral impurity and partial distinguishability. In both discrete and continuous variable systems, these can be modeled as mixed Gaussian states over multiple modes. The modes are not all resolved by photon-number measurements and so require calculations on coarse-grained photon-number distribution. Existing methods can lead to a combinatorial explosion in the time complexity, making this task unfeasible for even moderate sized experiments of interest. In this work, we prove that the computation of this type of distributions can be done in exponential time, providing a combinatorial speedup. We develop numerical techniques that allow us to determine coarse-grained photon number distributions of Gaussian states, as well as density matrix elements of heralded non-Gaussian states prepared in the presence of spectral impurity and partial distinguishability. These results offer significant speed-up and accuracy improvements to validation tests of both Fock and Gaussian boson samplers that rely on binned probability distributions.
We analyze the impact of photon loss on the photon-number statistics of Gaussian states. Specifically, we propose and carefully evaluate several methods to mitigate deviations in the photon-number distributions of lossy (displaced) squeezed vacuum states from those of their lossless counterparts. These methods rely on appropriately redefining the parameters of Gaussian states when the loss budget is known in order to recover, as closely as possible, the desired photon-number distribution associated with each target state. While it is intrinsically hard to directly optimize the photon-number distribution of high-dimensional, correlated multimode Gaussian states, the proposed methods are instead based on optimizing specific key properties such as fidelity, phase-space functions, low-order moments of the underlying photon-number statistics, or overlap with the vacuum state. In particular, our results show that optimizing the fidelity between a pure Gaussian target state and a modified Gaussian state that has passed through a loss channel does typically not result in closeness of the corresponding photon-number distributions. Furthermore, we show that correcting for the vacuum overlap minimizes the deviation in the photon-number distribution for large parameter ranges which we explicitly prove for single-mode squeezed vacuum and provide numerical evidence for general (displaced) squeezed vacuum states. As photon loss is a key limitation for Gaussian boson sampling, our results provide insights into the feasibility and limitations of such photonic quantum simulations in lossy environments and offer guidelines for mitigating these imperfections.