Sampling experiments provide a viable route to show quantum advantages of quantum devices over classical computers in well-defined computational tasks. However, quantum devices such as boson samplers are susceptible to various errors, including stochastic errors due to fabrication imperfections. These cause the implemented unitary operations to deviate randomly from their intended targets, following distributions with finite variance. Whilst full-scale quantum error correction remains challenging in the near term, quantum error mitigation schemes have been devised to estimate expectation values, but it is unclear how these schemes would work for sampling experiments. In this work, we demonstrate that, given access to multiple stochastic unitaries, it is possible to mitigate the effect of these errors in sampling experiments. We adopt the unitary averaging protocol which employs multiple stochastic boson samplers to generate a distribution that approximates the ideal boson sampler distribution as the number of samplers increases. We derive a rigorous upper bound on the trace distance between the output probability distributions induced by invertible vacuum-heralded networks based on the Schur-Weyl duality. This result can be seen concretely as an error mitigation scheme in sampling experiments against stochastic errors. On a broader level, it suggests a path towards understanding error mitigation for sampling experiments and developing analysis tools for photonic circuits incorporating measurements and feed-forward. We further provide other applications of unitary averaging, including its use in implementing the linear combination of unitaries and benchmarking fabrication repeatability in linear optics.
Sampling experiments provide a viable route to show quantum advantages of quantum devices over classical computers in well-defined computational tasks. However, quantum devices such as boson samplers are susceptible to various errors, including stochastic errors due to fabrication imperfections. These cause the implemented unitary operations to deviate randomly from their intended targets, following distributions with finite variance. Whilst full-scale quantum error correction remains challenging in the near term, quantum error mitigation schemes have been devised to estimate expectation values, but it is unclear how these schemes would work for sampling experiments. In this work, we demonstrate that, given access to multiple stochastic unitaries, it is possible to mitigate the effect of these errors in sampling experiments. We adopt the unitary averaging protocol which employs multiple stochastic boson samplers to generate a distribution that approximates the ideal boson sampler distribution as the number of samplers increases. We derive a rigorous upper bound on the trace distance between the output probability distributions induced by invertible vacuum-heralded networks based on the Schur-Weyl duality. This result can be seen concretely as an error mitigation scheme in sampling experiments against stochastic errors. On a broader level, it suggests a path towards understanding error mitigation for sampling experiments and developing analysis tools for photonic circuits incorporating measurements and feed-forward. We further provide other applications of unitary averaging, including its use in implementing the linear combination of unitaries and benchmarking fabrication repeatability in linear optics.
Recent experimental breakthroughs have signalled the imminent arrival of the early fault-tolerant era. However, for a considerable period in the foreseeable future, relying solely on quantum error correction for full error suppression will remain extremely challenging due to its substantial hardware overhead. Additional help from quantum error mitigation (QEM) is essential for bridging this gap towards achieving quantum advantage. The application of QEM has so far been restricted to expectation value estimation, leaving its extension to sampling-based algorithms -- which is expected to play a pivotal role in the early fault-tolerant era -- an unresolved challenge. In this work, we present a framework for applying any QEM techniques to obtain the error-mitigated output distribution, showing that this incurs no greater cost than estimating a single observable. We also devised a way to sample from this distribution and constructed an explicit scheme for applying any QEM methods to quantum phase estimation, which can be generalised to other sampling algorithms. Numerical experiments were conducted to validate the efficacy of these methods. We believe our methods significantly broaden the scope of QEM, extending its applicability to most algorithms of practical interest and forming a crucial step towards realising quantum advantage.
Abstract We experimentally demonstrate a testing strategy for boson samplers that is based on efficiently computable expressions for the output photon counting distributions binned over multiple optical modes. We apply this method to validate boson sampling experiments with three photons on a reconfigurable photonic chip, which implements a four-mode interferometer, analyzing 50 Haar-random unitary transformations while tuning photon indistinguishability via controlled delays. We show that for (nearly) indistinguishable photons, the experiment accurately reproduces the ideal boson sampling binned-mode distributions, which exhibit variations that depend both on the specific interferometer implemented as well as the choice of bin, confirming the usefulness of the method to diagnose imperfections such as partial distinguishability or imperfect chip control. Finally, we analyze the behavior of Haar-averaged binned-mode distributions in relation to indistinguishability and demonstrate analytically that its variance is proportional to the average of the square of the photons’ indistinguishability parameter. These findings highlight the central role of binning in boson sampling validation, offering a scalable and efficient framework for assessing multiphoton interference and experimental performance.
Sampling from probability distributions of quantum circuits is a fundamentally and practically important task which can be used to demonstrate quantum supremacy using noisy intermediate-scale quantum devices. In the present work, we examine classical simulability of sampling from the output photon-number distribution of linear-optical circuits composed of random beam splitters with equally distributed squeezed vacuum states and single-photon states input. We provide efficient classical algorithms to simulate linear-optical random circuits and show that the algorithms' error is exponentially small up to a depth less than quadratic in the distance between sources using a classical random walk behavior of random linear-optical circuits. Notably, the average-case depth allowing an efficient classical simulation is larger than the worst-case depth limit, which is linear in the distance. Besides, our results together with the hardness of boson sampling give a lower-bound on the depth for constituting global Haar-random unitary circuits.
Despite rapid advances in quantum hardware, noise remains a central obstacle to deploying quantum algorithms on near-term devices. In particular, random coherent errors that accumulate during circuit execution constitute a dominant and fundamentally challenging noise source. We introduce a noise-resilient framework for Quantum Signal Processing (QSP) that mitigates such coherent errors without increasing circuit depth or ancillary qubit requirements. Our approach uses ensembles of noisy QSP circuits combined with measurement-level averaging to suppress random phase errors in Z rotations. Building on this framework, we develop robust QSP algorithms for implementing polynomial functions of Hermitian matrices and for estimating observables, with applications to Hamiltonian simulation, quantum linear systems, and ground-state preparation. We analyze the trade-off between approximation error and hardware noise, which is essential for practical implementation under the stringent depth and coherence constraints of current quantum hardware. Our results establish a practical pathway for integrating error mitigation seamlessly into algorithmic design, advancing the development of robust quantum computing, and enabling the discovery of scientific applications with near- and mid-term quantum devices.
Abstract Measurement errors limit the fidelity of sampling-based quantum tasks, yet extending error mitigation to sampling remains challenging. Here, we introduce a simple and scalable repetition-based scheme that mitigates measurement errors directly at the sampling level. By encoding each data qubit into multiple physical qubits prior to readout and applying majority-vote decoding, physical errors are converted into a higher-order process, enabling exponential suppression that persists despite realistic gate imperfections during encoding. Using experimental noise profiles from leading superconducting platforms, we demonstrate that a modest linear qubit overhead suppresses measurement errors by several orders of magnitude, significantly boosting fidelities for typical sampling tasks, such as measurement-based quantum computation and random circuit sampling. These results establish a practical, hardware-compatible, and scalable approach for improving the reliability of near-term quantum sampling experiments.
Boson sampling is a candidate for quantum protocols to truly realize the quantum computation advantage and to be used in advanced fields where complex computations are needed, such as quantum chemistry. However, this proposal is hard to achieve due to the existence of noise sources such as photon losses. In order to quantificationally analyze the influences of photon losses in optical networks, boson sampling is classically simulated based on the equivalent beam splitter mechanism, where the photon loss happening in optical units is equivalent to the photon transmission into the environmental paths through a virtual beam splitter. In our simulation, networks corresponding to random unitary matrices are made up, considering both the Reck structure and the Clements structure. The photon loss probability in an optical unit is well controlled by adjusting the parameters of the virtual beam splitter. Therefore, to simulate boson sampling with photon losses in optical networks is actually to simulate ideal boson sampling with more modes. It is found that when the photon loss probability is constant, boson sampling with Clements structures distinctly performs much better than that with Reck structures. Furthermore, the photon loss probability is also set to follow the normal distribution, which is thought to be closer to the situation in reality. It is found that when the mean value of photon loss probability is constant, for both network structures, errors of outputs become more obvious with the increase of standard deviation. It can be inferred that the increase of error rate can be explained by the network depth and the conclusion is suitable for larger-scale boson sampling. Finally, the number of output photons is taken into consideration, which is directly related to the classical computation complexity. It is found that with the photon loss probability, the ratio of output combinations without photon losses decreases sharply, implying that photon losses can obviously affect the quantum computation advantage of boson sampling. Our results indicate that photon losses can result in serious errors for boson sampling, even with a stable network structure such as that of Clements. This work is helpful for boson sampling experiments in reality and it is desired to develop a better protocol, for example, a well-designed network or excellent optical units, to well suppress photon losses.