ABSTRACT In linear optics, photons are scattered in a network through passive optical elements including beam splitters and phase shifters, leading to many intriguing applications in physics, such as Mach–Zehnder interferometry, the Hong–Ou–Mandel effect, and tests of fundamental quantum mechanics. Here we present the fundamental limit in the transition amplitudes of bosons, applicable to all physical linear optical networks. Apart from boson sampling, this transition bound results in many other interesting applications, including behaviors of Bose–Einstein condensates (BEC) in optical networks, counterparts of Hong–Ou–Mandel effects for multiple photons, and approximating permanents of matrices. In addition, this general bound implies the existence of a polynomial-time randomized algorithm for estimating the transition amplitudes of bosons, which represents a solution to an open problem raised by Aaronson and Hance (Quantum Inf Comput 2012; 14: 541–59). Consequently, this bound implies that computational decision problems encoded in linear optics, prepared and detected in the Fock basis, can be solved efficiently by classical computers within additive errors. Furthermore, our result also leads to a classical sampling algorithm that can be applied to calculate the many-body wave functions and the S-matrix of bosonic particles.
ABSTRACT In linear optics, photons are scattered in a network through passive optical elements including beam splitters and phase shifters, leading to many intriguing applications in physics, such as Mach–Zehnder interferometry, the Hong–Ou–Mandel effect, and tests of fundamental quantum mechanics. Here we present the fundamental limit in the transition amplitudes of bosons, applicable to all physical linear optical networks. Apart from boson sampling, this transition bound results in many other interesting applications, including behaviors of Bose–Einstein condensates (BEC) in optical networks, counterparts of Hong–Ou–Mandel effects for multiple photons, and approximating permanents of matrices. In addition, this general bound implies the existence of a polynomial-time randomized algorithm for estimating the transition amplitudes of bosons, which represents a solution to an open problem raised by Aaronson and Hance (Quantum Inf Comput 2012; 14: 541–59). Consequently, this bound implies that computational decision problems encoded in linear optics, prepared and detected in the Fock basis, can be solved efficiently by classical computers within additive errors. Furthermore, our result also leads to a classical sampling algorithm that can be applied to calculate the many-body wave functions and the S-matrix of bosonic particles.
In linear optics, photons are scattered in a network through passive optical elements including beamsplitters and phase shifters, leading to many intriguing applications in physics, such as Mach-Zehnder interferometry, Hong-Ou-Mandel effect, and tests of fundamental quantum mechanics. Here we present a general analytic expression governing the upper limit of the transition amplitudes in sampling bosons, through all realizable linear optics. Apart from boson sampling, this transition bound results in many other interesting applications, including behaviors of Bose-Einstein Condensates (BEC) in optical networks, counterparts of Hong-Ou-Mandel effects for multiple photons, and approximating permanents of matrices. Also, this general bound implies the existence of a polynomial-time randomized algorithm for estimating transition amplitudes of bosons, which represents a solution to an open problem raised by Aaronson and Hance in 2012.
Boson Sampling is a prominent candidate for the demonstration of quantum computational advantage but it remains unclear whether a large-scale boson sampler can find useful computational applications. The challenge is that to use a boson sampler to estimate, for example, a physical observable, it is necessary to coarse grain the outcome distribution, owing to the exponential size of the outcome space and the anti-concentration properties of the Boson Sampling distribution. In this work, we analyse the complexity of a specific type of coarse-graining of Boson Sampling distributions based on linear functions of the output photon occupation numbers, which we refer to as linear statistics. We present an efficient classical algorithm for approximating linear statistics of boson samplers within additive error for different kinds of input states, such as Fock states and squeezed states. This allows us to unify in the same framework recent results on efficient quantum-inspired classical algorithms for simulating molecular vibronic spectra, or efficient classical approximations of coarse-grained distributions based on detector binning. Additionally, we show how our algorithm can be used to classically evaluate a proposed one-way function based on Boson Sampling, while other proposals for cryptographic applications escape our classical simulation techniques. We leave open the question of classical simulability of non-linear statistics and connect it to the problem of computing transition amplitudes of linear-optical circuits with one layer of interactions.
Abstract A fundamental question in linear optical quantum computing is to understand the origin of the quantum supremacy in the physical system. It is found that the multimode linear optical transition amplitudes are calculated through the permanents of transition operator matrices, which is a hard problem for classical simulations (boson sampling problem). We can understand this problem by considering a quantum measure that directly determines the runtime for computing the transition amplitudes. In this paper, we suggest a quantum measure named “Fock state concurrence sum” C S , which is the summation over all the members of “the generalized Fock state concurrence” (a measure analogous to the generalized concurrences of entanglement and coherence). By introducing generalized algorithms for computing the transition amplitudes of the Fock state boson sampling with an arbitrary number of photons per mode, we show that the minimal classical runtime for all the known algorithms directly depends on C S . Therefore, we can state that the Fock state concurrence sum C S behaves as a collective measure that controls the computational complexity of Fock state BS. We expect that our observation on the role of the Fock state concurrence in the generalized algorithm for permanents would provide a unified viewpoint to interpret the quantum computing power of linear optics.
We give new evidence that quantum computers -- moreover, rudimentary quantum computers built entirely out of linear-optical elements -- cannot be efficiently simulated by classical computers. In particular, we define a model of computation in which identical photons are generated, sent through a linear-optical network, then nonadaptively measured to count the number of photons in each mode. This model is not known or believed to be universal for quantum computation, and indeed, we discuss the prospects for realizing the model using current technology. On the other hand, we prove that the model is able to solve sampling problems and search problems that are classically intractable under plausible assumptions. Our first result says that, if there exists a polynomial-time classical algorithm that samples from the same probability distribution as a linear-optical network, then P^#P=BPP^NP, and hence the polynomial hierarchy collapses to the third level. Unfortunately, this result assumes an extremely accurate simulation. Our main result suggests that even an approximate or noisy classical simulation would already imply a collapse of the polynomial hierarchy. For this, we need two unproven conjectures: the "Permanent-of-Gaussians Conjecture", which says that it is #P-hard to approximate the permanent of a matrix A of independent N(0,1) Gaussian entries, with high probability over A; and the "Permanent Anti-Concentration Conjecture", which says that |Per(A)|>=sqrt(n!)/poly(n) with high probability over A. We present evidence for these conjectures, both of which seem interesting even apart from our application. This paper does not assume knowledge of quantum optics. Indeed, part of its goal is to develop the beautiful theory of noninteracting bosons underlying our model, and its connection to the permanent function, in a self-contained way accessible to theoretical computer scientists.
We present a classical algorithm for approximating the expectation values of observables in linear optical circuits with arbitrary product input states, achieving additive-error accuracy. This result indicates that current applications of photonic systems aimed at demonstrating practical quantum advantage through expectation value estimation, such as photonic variational algorithms, may face challenges in attaining the computational advantage. It also implies the output probabilities of boson sampling with arbitrary product input states can be efficiently approximated by our method, so that boson sampling becomes efficiently simulable when its output probability distribution is polynomially sparse. We also develop an efficient classical algorithm for estimating transition amplitudes of arbitrary product states in linear optical circuits. This provides additive-error approximation algorithms for matrix functions associated with linear optical circuits, such as the (loop-)Hafnian, which are of independent interest. As an application, it solves the generalized molecular vibronic spectra problem [C. Oh et al., Nat. Phys. 20, 225 (2024)10.1038/s41567-023-02308-9], previously suggested as a candidate for practical quantum advantage. Finally, we extend our framework to near-Clifford circuits, enabling classical approximation of their expectation values.
In his famous 1982 lecture, Nobel laureate Richard Feynman seeded the idea of quantum computation. He reasoned that the amount of classical information required to describe a quantum system scaled exponentially with its size, and thus could not be efficiently simulated by a classical computer. Indeed, a fairly modest system of 100 quantum bits would already require more information to describe than all the present data storage on Earth [1]. Given this, there would likely be computational tasks that achieve quantum supremacy—impossible to perform practically without the aid of a quantum information processor. In the 37 intervening years, these initial ideas have solidified, as the discovery of a string of quantum algorithms—from efficient factoring of large numbers to simulating chemical reactions—fueled a race for quantum supremacy. Yet, even theoretically, conclusively deciding—if quantum supremacy has been achieved—is non-trivial. Judging what problems can be efficiently solved by classical computers turns out to be exceedingly hard in itself. Take factoring, for instance. It is really classical intractable, or is it merely that we have not yet found an efficient classical means to factor? As we better understood how quantum computers process information, we also got better at finding classical methods to simulate them. This in turn has inspired new efficient classical algorithms in surprising situations [2]. Prof. Man-Hong Yung, associate professor at the Southern University of Science and Technology, and colleagues provided a beautiful illustration of such a discovery. Their work has been published in the July 2019 issue of National Science Review [3]. It considers a specialized type of linear optical quantum processor. The input and output are encoded within an array of optical modes by adjusting their initial photon number. Computation is achieved by interfering with these modes via a network of beam splitters, such that photons in one mode are transferred to another. Such processors have captured scientific attention due to their ability to perform a task known as boson sampling. When each mode is initialized with a single photon, the output statistics allows efficient evaluation of the matrix permanent—a computational problem with very strong evidence for being classically intractable [4]. This, coupled with their relative ease of engineering, made this specialized type of quantum processor an ideal holy grail for demonstrating quantum supremacy. Given this, one may expect that classical computers can say little about what such processors will output. After all, the number of paths that a photon can take to reach one mode from another scales exponentially with the size of the beam-splitter network. One may therefore assume that it would take exponential resources to track these paths to determine the probability of a photon moving from mode A to mode B. However, Yung and colleagues showed that a classical algorithm can efficiently estimate the transition probabilities—as long as one allows for additive error. This result represents a significant generalization of prior results by Scott Aaronson, which considers the case where there exists at most one photon in each mode [5]. In contrast, Yung's solution works regardless of the number of photons encoded in each mode. The result not only benefits the simulation of bosonic interference, but also has strong consequences in computing. Suppose that we use such specialized quantum processors to make binary decisions. That is, depending on input, we feed in particular distributions of photons—and then decide to take a particular action depending on whether the probability of seeing x photons in a particular output is above some threshold. In such contexts, a classical algorithm that estimates this probability up to the bounded additive error is sufficient for making identical decisions as the quantum processer itself. As such, Yung's algorithm enables a classical computer to be able to efficiently mimic the decisions of its specialized linear optical counterpart. This will no doubt have important ramifications in the present race for quantum supremacy. Much present effort is spent on demonstrating quantum supremacy through boson sampling, and the results here indicate that we will need to introduce more sophisticated quantum architectures if we wish to demonstrate such supremacy in making the right decisions.