Search papers, labs, and topics across Lattice.
We show that classical simulation of lossy Gaussian Boson Sampling using classical light replacement is limited to regimes where surviving photons grow as $\sqrt{N}$, proving this is a fundamental limit for non-negative P-representation states. To overcome this, we introduce a novel approach that approximates the conditional state after detection with a Gaussian, achieving $O(N\eta^3)$ total variation error at polynomial cost. This allows for simulating scenarios where the mean surviving photon number grows as $N^{2/3}$ (or $N^{3/4}$ with a specific construction) at fixed accuracy, significantly extending the classically simulable regime.
Classical simulation of lossy GBS breaks the $\sqrt{N}$ barrier by approximating conditional states post-detection, enabling simulations with $N^{2/3}$ surviving photons.
Photon loss makes Gaussian boson sampling easier to simulate classically. A standard approach replaces the lossy input by classical light, and it stays accurate while the mean number of surviving photons grows no faster than the square root of the number $N$ of squeezed inputs at fixed squeezing and accuracy. We prove that this scale is a limit of replacement by classical light. For explicit interferometers, no state with a nonnegative Glauber--Sudarshan $P$ representation, Gaussian or not and however correlated, keeps the error in the photon count distribution small beyond this scale. We then surpass this limit by approximating the detection instead of the input. The sampler evolves the state exactly between detections and, after each detection, replaces the conditional state by a Gaussian with the same covariance. For every interferometer its total variation error is $O(Nη^3)$ at polynomial expected cost, where $η$ is the photon survival probability. The mean surviving photon number can therefore grow as $N^{2/3}$ at fixed accuracy. A separate pair construction reaches $N^{3/4}$ when every input mode is equally squeezed.