Abstract:Gaussian boson sampling (GBS) is considered a candidate problem for demonstrating quantum advantage. We propose an algorithm for the approximate classical simulation of a lossy GBS instance. The algorithm relies on the Taylor series expansion, and increasing the number of terms of the expansion that are used in the calculation yields greater accuracy. The complexity of the algorithm is polynomial in the number of modes given the number of terms is fixed. We describe conditions for the input state squeezing parameter and loss level that provide the best efficiency for this algorithm (by efficient, we mean that the Taylor series converges quickly). In recent experiments that claim to have demonstrated quantum advantage, these conditions are satisfied; thus, this algorithm can be used to classically simulate these experiments.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How to simulate the lossy Gaussian Boson Sampling (GBS) problem through classical computers in order to verify the results in quantum advantage experiments**. ### Problem Background Gaussian Boson Sampling (GBS) is a candidate problem widely considered to be able to demonstrate quantum advantage. However, actual experimental devices always have various imperfections, such as photon loss, uncertainty of interferometer parameters, etc. These imperfect factors make the classical simulation more complex and challenging. ### Core Contributions of the Paper The author proposes a classical algorithm based on Taylor series expansion for approximately simulating lossy Gaussian Boson Sampling instances. The main features of this algorithm are as follows: 1. **Taylor Series Expansion**: The algorithm depends on Taylor series expansion. As the number of terms used increases, the computational accuracy also improves. 2. **Polynomial Complexity**: Given the number of modes, the time complexity of the algorithm is polynomial, provided that the number of terms in the Taylor series is fixed. 3. **Optimal Condition Analysis**: The author analyzes the optimal conditions of the squeezing parameters and loss levels of the input state, so that the Taylor series can converge quickly. These conditions are consistent with recent experiments claiming to demonstrate quantum advantage, so this algorithm can be used for classical simulation of these experiments. ### Specific Presentation of Formulas The key formulas involved in the paper include: - **Squeezed - state Representation of the Input State**: \[ |\psi\rangle = e^{\sum_i \frac{\alpha_i}{2} (\hat{a}_i^\dagger)^2} |0\rangle \] - **Density Matrix after Considering Loss**: \[ \hat{\rho} = \text{Tr}_{\text{loss}} \left[ e^{\sum_i \frac{\alpha_i}{2} (c \hat{a}_i^\dagger + s \hat{b}_i^\dagger)^2} |0_a 0_b\rangle \langle 0_a 0_b| e^{\sum_i \frac{\alpha_i}{2} (c \hat{a}_i + s \hat{b}_i)^2} \right] \] - **Taylor Series Expansion**: \[ e^{(ξ_0 + χ)c \hat{a}^\dagger}|0\rangle \langle 0| e^{(ξ_0 + \tilde{χ})c \hat{a}} \] ### Experimental Verification The author verifies the effectiveness of his algorithm by analyzing recent Gaussian Boson Sampling experimental data. The results show that under the conditions of high loss and low squeezing parameters, the ε value is small, the Taylor series can converge quickly, and thus the lost contribution can be controlled within 1% at the fourth - order expansion. ### Summary The classical algorithm proposed in this paper can efficiently simulate the lossy Gaussian Boson Sampling problem under certain conditions and verifies the experimental results demonstrating quantum advantage recently. This provides a powerful tool for further research on quantum advantage and helps to understand the boundary between quantum computing and classical computing.

Classical Modeling of a Lossy Gaussian Bosonic Sampler

Classical modelling of a lossy Gaussian bosonic sampler

Classical algorithm for simulating experimental Gaussian boson sampling

Classically simulating near-term partially-distinguishable and lossy boson sampling

The boundary for quantum advantage in Gaussian boson sampling

Gaussian boson sampling at finite temperature

Simulating Gaussian boson sampling quantum computers

Simulating lossy Gaussian boson sampling with matrix product operators

Gaussian boson sampling validation via detector binning

The Complexity of Bipartite Gaussian Boson Sampling

Effect of partial distinguishability on quantum supremacy in Gaussian Boson sampling

Classical models may be a better explanation of the Jiuzhang 1.0 Gaussian Boson Sampler than its targeted squeezed light model

Classical benchmarking of Gaussian Boson Sampling on the Titan supercomputer

Variational Tensor Network Simulation of Gaussian Boson Sampling and Beyond

Classical simulation of non-Gaussian bosonic circuits

Speedup in Classical Simulation of Gaussian Boson Sampling

Speeding up the classical simulation of Gaussian boson sampling with limited connectivity

Experimental demonstration of Gaussian boson sampling with displacement

Linear cross-entropy certification of quantum computational advantage in Gaussian Boson Sampling

Classical Simulation of Boson Sampling Based on Graph Structure