Classical Simulation of Lossy Gaussian Boson Sampling beyond Square-Root Regime
Abstract
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 N2/3 at fixed accuracy. A separate pair construction reaches N3/4 when every input mode is equally squeezed.
Related
- Classical Algorithms for Function Computation in Gaussian Boson Sampling
- A Topology-Aware Quantum-Classical Framework for Multi-Target Active Debris Removal using Gaussian Boson Sampling
- Improved classical optimization of fermionic Hamiltonians beyond Gaussian states
- Classical Simulation and Design Frontiers for IBM's Doped Clifford Sampling Experiment
- A Chip-scale Space-time Multiplexed Gaussian Boson Sampling Processor Beyond 10,000 Photons
BibTeX
@misc{lim2026classical,
title = {{Classical Simulation of Lossy Gaussian Boson Sampling beyond Square-Root Regime}},
author = {Youngrong Lim},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.06526},
eprint = {2610.06526},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.06526},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Classical Simulation of Lossy Gaussian Boson Sampling beyond Square-Root Regime
AU - Youngrong Lim
PY - 2026
DA - 2026/10/05
PB - arXiv
DO - 10.48550/arxiv.2610.06526
UR - https://arxiv.org/abs/2610.06526
AB - 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.
ER -