The Classical Complexity of Boson Sampling

We study the classical complexity of the exact Boson Sampling problem where the objective is to produce provably correct random samples from a particular quantum mechanical distribution. The computational framework was proposed by Aaronson and Arkhipov in 2011 as an attainable demonstration of `quantum supremacy', that is a practical quantum computing experiment able to produce output at a speed beyond the reach of classical (that is non-quantum) computer hardware. Since its introduction Boson Sampling has been the subject of intense international research in the world of quantum computing. On the face of it, the problem is challenging for classical computation. Aaronson and Arkhipov show that exact Boson Sampling is not efficiently solvable by a classical computer unless $P^{\#P} = BPP^{NP}$ and the polynomial hierarchy collapses to the third level. The fastest known exact classical algorithm for the standard Boson Sampling problem takes $O({m + n -1 \choose n} n 2^n )$ time to produce samples for a system with input size $n$ and $m$ output modes, making it infeasible for anything but the smallest values of $n$ and $m$. We give an algorithm that is much faster, running in $O(n 2^n + \operatorname{poly}(m,n))$ time and $O(m)$ additional space. The algorithm is simple to implement and has low constant factor overheads. As a consequence our classical algorithm is able to solve the exact Boson Sampling problem for system sizes far beyond current photonic quantum computing experimentation, thereby significantly reducing the likelihood of achieving near-term quantum supremacy in the context of Boson Sampling.

[1]  Leslie G. Valiant,et al.  The Complexity of Computing the Permanent , 1979, Theor. Comput. Sci..

[2]  Andrew G. White,et al.  Photonic Boson Sampling in a Tunable Circuit , 2012, Science.

[3]  Christian P. Robert,et al.  Bayesian computation: a summary of the current state, and samples backwards and forwards , 2015, Statistics and Computing.

[4]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[5]  H. Neven,et al.  Characterizing quantum supremacy in near-term devices , 2016, Nature Physics.

[6]  Nicolò Spagnolo,et al.  Experimental scattershot boson sampling , 2015, Science Advances.

[7]  Andreas Buchleitner,et al.  Stringent and efficient assessment of boson-sampling devices. , 2013, Physical review letters.

[8]  Michael J. Bremner,et al.  Quantum sampling problems, BosonSampling and quantum supremacy , 2017, npj Quantum Information.

[9]  Albert Nijenhuis,et al.  Combinatorial Algorithms for Computers and Calculators , 1978 .

[10]  David Bruce Wilson,et al.  Exact sampling with coupled Markov chains and applications to statistical mechanics , 1996, Random Struct. Algorithms.

[11]  William Feller,et al.  An Introduction to Probability Theory and Its Applications , 1967 .

[12]  Aram W. Harrow,et al.  Quantum computational supremacy , 2017, Nature.

[13]  Scott Aaronson,et al.  Bosonsampling is far from uniform , 2013, Quantum Inf. Comput..

[14]  Andreas Buchleitner,et al.  Statistical benchmark for BosonSampling , 2014, 1410.8547.

[15]  Nicolò Spagnolo,et al.  Experimental validation of photonic boson sampling , 2014, Nature Photonics.

[16]  Yang Wang,et al.  Computing Permanents for Boson Sampling on Tianhe-2 Supercomputer , 2016, National Science Review.

[17]  V. Shchesnovich,et al.  Asymptotic evaluation of bosonic probability amplitudes in linear unitary networks in the case of large number of bosons , 2013, 1304.6675.

[18]  Peter W. Shor,et al.  Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer , 1995, SIAM Rev..

[19]  Christian Schneider,et al.  Multi-photon boson-sampling machines beating early classical computers , 2016 .

[20]  David P. DiVincenzo,et al.  Adaptive quantum computation, constant depth quantum circuits and arthur-merlin games , 2002, Quantum Inf. Comput..

[21]  Andreas Björklund Below All Subsets for Some Permutational Counting Problems , 2016, SWAT.

[22]  Pavlos S. Efraimidis,et al.  Weighted Random Sampling over Data Streams , 2010, Algorithms, Probability, Networks, and Games.

[23]  Raphaël Clifford,et al.  Classical boson sampling algorithms with superior performance to near-term experiments , 2017, Nature Physics.

[24]  B. J. Metcalf,et al.  Boson Sampling on a Photonic Chip , 2012, Science.

[25]  A. Crespi,et al.  Integrated multimode interferometers with arbitrary designs for photonic boson sampling , 2013, Nature Photonics.

[26]  Fabio Sciarrino,et al.  Towards quantum supremacy with lossy scattershot boson sampling , 2016, 1610.02279.

[27]  David G. Glynn,et al.  The permanent of a square matrix , 2010, Eur. J. Comb..

[28]  Feller William,et al.  An Introduction To Probability Theory And Its Applications , 1950 .

[29]  Philip Walther,et al.  Experimental boson sampling , 2012, Nature Photonics.

[30]  Tim Hesterberg,et al.  Monte Carlo Strategies in Scientific Computing , 2002, Technometrics.

[31]  Greg Kuperberg,et al.  The bosonic birthday paradox , 2011, 1106.0849.

[32]  Lov K. Grover A fast quantum mechanical algorithm for database search , 1996, STOC '96.

[33]  Ashley Montanaro,et al.  No imminent quantum supremacy by boson sampling , 2017, 1705.00686.

[34]  Sheng-Tao Wang,et al.  Certification of Boson Sampling Devices with Coarse-Grained Measurements , 2016, 1601.02627.

[35]  Scott Aaronson,et al.  The computational complexity of linear optics , 2010, STOC '11.

[36]  A. J. Walker New fast method for generating discrete random numbers with arbitrary frequency distributions , 1974 .

[37]  R. Kronmal,et al.  On the Alias Method for Generating Random Variables From a Discrete Distribution , 1979 .