© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2610.16719

Realistic simulation of universal multiport interferometers with losses within a Gaussian formalism

Sebastiano Corli, Dario Chemoli, Samuele Altilia, Edoardo Suerra, Dario Tamascelli, Simone Cialdi, Enrico Prati

quant-phSubmitted Oct 12, 2026• Announced Oct 13, 2026• 0 citations

Abstract

Photonic technologies are among the leading platform contenders for demonstrating quantum advantage. Nevertheless, photon losses represent the major hurdle for the realization of a fault-tolerant optical quantum computer, and the need for an accurate characterization of this phenomenon is crucial for classical simulations of quantum algorithms on noisy backends. The noisy regime would require the simulation of the algorithm on a Fock backend or rather assuming the a priori hypothesis on the uniformity of the optical losses, which may lead to wrong predictions about the outcome from the algorithm on a cerain hardware. We propose a Gaussian model that can reproduce a lossy dynamics through the experimental characterization of imperfect interferometers, yielding the corresponding mean vector and the covariance matrix. Such model can be adapted to any continuous-variable architecture (spatial or time-bin encoded) and for both the Gaussian boson sampling or general-dyne based algorithms. Eventually, we quantify the drift of optical losses from being uniform throughout the multiport interferometer by a numerical simulation. This work proposes a blueprint for future characterizations, tailor-suited for specific devices and architectures, as well as for potential error correction and mitigation techniques, and it paves the way for realistic simulations of quantum algorithms on NISQ optical platforms.

Related papers

Regimes of Classical Simulability for Noisy Gaussian Boson Sampling

Haoyu Qi, Daniel J. Brod, Nicolás Quesada +1 more

As a promising candidate for exhibiting quantum computational supremacy, Gaussian boson sampling (GBS) is designed to exploit the ease of experimental preparation of Gaussian states. However, sufficiently large and inevitable experimental noise might render GBS classically simulable. In this work, we formalize this intuition by establishing a sufficient condition for approximate polynomial-time classical simulation of noisy GBS-in the form of an inequality between the input squeezing parameter, the overall transmission rate, and the quality of photon detectors. Our result serves as a nonclassicality test that must be passed by any quantum computational supremacy demonstration based on GBS. We show that, for most linear-optical architectures, where photon loss increases exponentially with the circuit depth, noisy GBS loses its quantum advantage in the asymptotic limit. Our results thus delineate intermediate-sized regimes where GBS devices might considerably outperform classical computers for modest noise levels. Finally, we find that increasing the amount of input squeezing is helpful to evade our classical simulation algorithm, which suggests a potential route to mitigate photon loss.

arXiv:1905.12075•Physical Review Letters•Mar 13, 2020•87 citations

Boosting Gaussian Boson Sampling using Optical Parametric Amplification Networks

Yukuan Zhao, Xiao-Ye Xu, Chuan-Feng Li +1 more

Gaussian Boson Sampling (GBS) provides a route toward demonstrating quantum computational advantage. However, optical loss, which reduces the entanglement in the system, can render GBS results classically simulable. We propose a nonlinear photonic architecture based on optical parametric amplifiers (OPAs) arranged in an interferometer network. This active configuration amplifies quantum correlations within the circuit while preserving the #P-hard Hafnian structure of the output probabilities. Using logarithmic negativity, we numerically show that entanglement scales linearly with both the OPA gain and network depth in the lossless limit, and maintains linear scaling with the number of modes under realistic loss rate. These scaling behaviors suggest that classical simulation in lossy scenarios remains computationally intractable. Our results demonstrate that OPA-boosted GBS preserves computational hardness in noisy environments, offering a more effective implementations of near-term photonic quantum computers.

quant-pharXiv:2512.00753•Physical Review A•Dec 2, 2025•0 citations

Sufficient conditions for hardness of lossy Gaussian boson sampling

Byeongseon Go, Changhun Oh, Hyunseok Jeong

Gaussian boson sampling (GBS) is a prominent candidate for the experimental demonstration of quantum advantage. However, while the current implementations of GBS are unavoidably subject to noise, the robustness of the classical intractability of GBS against noise remains largely unexplored. In this work, we establish the complexity-theoretic foundations for the classical intractability of noisy GBS under photon loss, which is a dominant source of imperfection in current implementations. We identify the loss threshold below which lossy GBS maintains the same complexity-theoretic level as ideal GBS, and show that this holds when at most a logarithmic fraction of photons is lost. We additionally derive an intractability criterion for the loss rate through a direct quantification of the statistical distance between ideal and lossy GBS. This work presents the first rigorous characterization of classically intractable regimes of lossy GBS, thereby serving as a crucial step toward demonstrating quantum advantage with near-term implementations.

quant-pharXiv:2511.07853•Nov 12, 2025•0 citations

Noise-tolerant tomography of multimode linear optical interferometers with single photons

Yu. A. Biriukov, R. D. Morozov, I. V. Dyakonov +11 more

Linear optical networks are fundamental to the advancement of quantum technologies, including quantum computing, communication, and sensing. The accurate characterization of these networks, described by unitary matrices, is crucial to their effective utilization and scalability. In this work, we present the method for reconstructing the transfer matrix of a linear optical interferometer based on the analysis of cross-correlation functions of photon counts between pairs of output modes. Our approach accounts for losses and photon indistinguishability, making it robust to experimental imperfections. By minimizing the requirements for the input states, the method simplifies the experimental implementation. We demonstrate the effectiveness of our technique through theoretical modeling and experimental validation in a 4-mode programmable integrated optical interferometer. The results show high fidelity in matrix reconstruction and successful application in boson sampling experiments. In addition, we provide a comprehensive formalism for correlation functions and discuss the robustness of the method to measurement errors. This work offers a practical and efficient solution for characterizing linear-optical networks, paving the way for scaling up photonic quantum technologies.

arXiv:2506.20490•Jun 25, 2025•0 citations

Simulating lossy and partially distinguishable quantum optical circuits: theory, algorithms, and applications to experiment validation and state preparation

Jacob F. F. Bulmer, Javier Martínez-Cifuentes, Bryn A. Bell +1 more

To understand quantum optics experiments, we must perform calculations that consider the principal sources of noise, such as losses, spectral impurity, and partial distinguishability. In both discrete and continuous variable systems, these can be modeled as mixed Gaussian states over multiple modes. The modes are not all resolved by photon-number measurements and so require calculations on coarse-grained photon-number distributions. Existing methods can lead to a combinatorial explosion in the time complexity, making this task unfeasible for even moderately sized experiments of interest. We prove that the computation of this type of distribution can be done in exponential time, providing a combinatorial speed-up. We develop numerical techniques that allow us to determine coarse-grained photon number distributions of Gaussian states, as well as density matrix elements of heralded non-Gaussian states prepared in the presence of spectral impurity and partial distinguishability. These results offer significant speed-up and accuracy improvements to validation tests of both Fock and Gaussian boson samplers that rely on binned probability distributions. Moreover, they pave the way to a more efficient simulation of realistic photonic circuits, unlocking the ability to perform exact calculations at scales that were previously out of reach. In addition, our results, including loop Hafnian master theorems, may be of interest to the fields of combinatorics and graph theory.

arXiv:2412.17742•Advanced Photonics•Feb 18, 2026•1 citation

A Recurrent Multiport Interferometer Architecture for Fault-Tolerant Programmable Photonics

Mihai Crisan, David A. Carpenter, Shamsul Arafin

We propose a novel multiport interferometer architecture that uses optical recursion to mitigate hardware errors and scales as O(N). We also show that it is possible to construct perfect transformations with imperfect beamsplitters meshes.

•Jan 1, 2024•0 citations

Photonic technologies are among the leading platform contenders for demonstrating quantum advantage. Nevertheless, photon losses represent the major hurdle for the realization of a fault-tolerant optical quantum computer, and the need for an accurate characterization of this phenomenon is crucial for classical simulations of quantum algorithms on noisy backends. The noisy regime would require the simulation of the algorithm on a Fock backend or rather assuming the a priori hypothesis on the uniformity of the optical losses, which may lead to wrong predictions about the outcome from the algorithm on a cerain hardware. We propose a Gaussian model that can reproduce a lossy dynamics through the experimental characterization of imperfect interferometers, yielding the corresponding mean vector and the covariance matrix. Such model can be adapted to any continuous-variable architecture (spatial or time-bin encoded) and for both the Gaussian boson sampling or general-dyne based algorithms. Eventually, we quantify the drift of optical losses from being uniform throughout the multiport interferometer by a numerical simulation. This work proposes a blueprint for future characterizations, tailor-suited for specific devices and architectures, as well as for potential error correction and mitigation techniques, and it paves the way for realistic simulations of quantum algorithms on NISQ optical platforms.

As a promising candidate for exhibiting quantum computational supremacy, Gaussian boson sampling (GBS) is designed to exploit the ease of experimental preparation of Gaussian states. However, sufficiently large and inevitable experimental noise might render GBS classically simulable. In this work, we formalize this intuition by establishing a sufficient condition for approximate polynomial-time classical simulation of noisy GBS-in the form of an inequality between the input squeezing parameter, the overall transmission rate, and the quality of photon detectors. Our result serves as a nonclassicality test that must be passed by any quantum computational supremacy demonstration based on GBS. We show that, for most linear-optical architectures, where photon loss increases exponentially with the circuit depth, noisy GBS loses its quantum advantage in the asymptotic limit. Our results thus delineate intermediate-sized regimes where GBS devices might considerably outperform classical computers for modest noise levels. Finally, we find that increasing the amount of input squeezing is helpful to evade our classical simulation algorithm, which suggests a potential route to mitigate photon loss.

Gaussian Boson Sampling (GBS) provides a route toward demonstrating quantum computational advantage. However, optical loss, which reduces the entanglement in the system, can render GBS results classically simulable. We propose a nonlinear photonic architecture based on optical parametric amplifiers (OPAs) arranged in an interferometer network. This active configuration amplifies quantum correlations within the circuit while preserving the #P-hard Hafnian structure of the output probabilities. Using logarithmic negativity, we numerically show that entanglement scales linearly with both the OPA gain and network depth in the lossless limit, and maintains linear scaling with the number of modes under realistic loss rate. These scaling behaviors suggest that classical simulation in lossy scenarios remains computationally intractable. Our results demonstrate that OPA-boosted GBS preserves computational hardness in noisy environments, offering a more effective implementations of near-term photonic quantum computers.

Gaussian boson sampling (GBS) is a prominent candidate for the experimental demonstration of quantum advantage. However, while the current implementations of GBS are unavoidably subject to noise, the robustness of the classical intractability of GBS against noise remains largely unexplored. In this work, we establish the complexity-theoretic foundations for the classical intractability of noisy GBS under photon loss, which is a dominant source of imperfection in current implementations. We identify the loss threshold below which lossy GBS maintains the same complexity-theoretic level as ideal GBS, and show that this holds when at most a logarithmic fraction of photons is lost. We additionally derive an intractability criterion for the loss rate through a direct quantification of the statistical distance between ideal and lossy GBS. This work presents the first rigorous characterization of classically intractable regimes of lossy GBS, thereby serving as a crucial step toward demonstrating quantum advantage with near-term implementations.

Linear optical networks are fundamental to the advancement of quantum technologies, including quantum computing, communication, and sensing. The accurate characterization of these networks, described by unitary matrices, is crucial to their effective utilization and scalability. In this work, we present the method for reconstructing the transfer matrix of a linear optical interferometer based on the analysis of cross-correlation functions of photon counts between pairs of output modes. Our approach accounts for losses and photon indistinguishability, making it robust to experimental imperfections. By minimizing the requirements for the input states, the method simplifies the experimental implementation. We demonstrate the effectiveness of our technique through theoretical modeling and experimental validation in a 4-mode programmable integrated optical interferometer. The results show high fidelity in matrix reconstruction and successful application in boson sampling experiments. In addition, we provide a comprehensive formalism for correlation functions and discuss the robustness of the method to measurement errors. This work offers a practical and efficient solution for characterizing linear-optical networks, paving the way for scaling up photonic quantum technologies.

To understand quantum optics experiments, we must perform calculations that consider the principal sources of noise, such as losses, spectral impurity, and partial distinguishability. In both discrete and continuous variable systems, these can be modeled as mixed Gaussian states over multiple modes. The modes are not all resolved by photon-number measurements and so require calculations on coarse-grained photon-number distributions. Existing methods can lead to a combinatorial explosion in the time complexity, making this task unfeasible for even moderately sized experiments of interest. We prove that the computation of this type of distribution can be done in exponential time, providing a combinatorial speed-up. We develop numerical techniques that allow us to determine coarse-grained photon number distributions of Gaussian states, as well as density matrix elements of heralded non-Gaussian states prepared in the presence of spectral impurity and partial distinguishability. These results offer significant speed-up and accuracy improvements to validation tests of both Fock and Gaussian boson samplers that rely on binned probability distributions. Moreover, they pave the way to a more efficient simulation of realistic photonic circuits, unlocking the ability to perform exact calculations at scales that were previously out of reach. In addition, our results, including loop Hafnian master theorems, may be of interest to the fields of combinatorics and graph theory.