Simulating Gaussian boson sampling on graphs in polynomial time
Konrad Anand, Zongchen Chen, Mary Cryan, Graham Freifeld, Leslie Ann Goldberg, Heng Guo, Xinyuan Zhang
TL;DR
The paper addresses whether Gaussian Boson Sampling on graphs, and a related non-negative Boson Sampling distribution, can be classically simulated in polynomial time. It achieves this by reducing the sampling task to weighted perfect-matchings on the Cartesian product graph $G\square K_2$ and applying the Jerrum–Sinclair chain with a carefully chosen weight scaling, yielding a runtime of $O(\overline{c} m n^4 \log^2(n\overline{c}/\varepsilon))$. For standard Boson Sampling with non-negative input matrices, it provides a similar polynomial-time sampler with runtime $O\left(\frac{m^7 n^{14}}{\varepsilon^7} \log^4\left(\frac{mn}{\varepsilon}\right)\right)$. Together, these results challenge the expectation of exponential quantum advantage for graph-based applications of BS/GBS in the non-negative regime and motivate further exploration of relaxing assumptions and linking to permanents in more general settings.
Abstract
We show that a distribution related to Gaussian Boson Sampling (GBS) on graphs can be sampled classically in polynomial time. Graphical applications of GBS typically sample from this distribution, and thus quantum algorithms do not provide exponential speedup for these applications. We also show that another distribution related to Boson sampling can be sampled classically in polynomial time.
