Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2407.06290 · Physical Review Letters

Gate-Based Quantum Simulation of Gaussian Bosonic Circuits on Exponentially Many Modes

Alice Barthe, Marco Cerezo, Andrew T. Sornborger, Martín Larocca, Diego García-Martín

Published Feb 20, 2025• 7 citations
arXivHTMLPDF

Abstract

We introduce a framework for simulating, on an (n+1)-qubit quantum computer, the action of a Gaussian bosonic (GB) circuit on a state over 2^{n} modes. Specifically, we encode the initial bosonic state's expectation values over quadrature operators (and their covariance matrix) as an input qubit state. This is then evolved by a quantum circuit that effectively implements the symplectic propagators induced by the GB gates. We find families of GB circuits and initial states leading to efficient quantum simulations. For this purpose, we introduce a dictionary that maps between GB and qubit gates such that particle- (non-particle-) preserving GB gates lead to real- (imaginary-) time evolutions at the qubit level. For the special case of particle-preserving circuits, we present a bounded-error-quantum-polynomial time (BQP)-complete GB decision problem, indicating that GB evolutions of Gaussian states on exponentially many modes are as powerful as universal quantum computers. We also perform numerical simulations of an interferometer on ∼8×10^{9} modes, illustrating the power of our framework.

DOI
10.1103/physrevlett.134.070604

Related papers

Gate-based emulation of boson sampling using photonic qubits

Aastha P. Zalone, S. P. Dinesh, C. M. Chandrashekar

Boson sampling arising from multiphoton interference in linear-optical networks is a prominent non-universal model for quantum computation. Here, by encoding the multi-qubit state to bosonic Fock state, we present a scalable quantum-circuit framework for simulating boson sampling on a universal quantum computing platform. Beginning with balanced beam-splitter transformations on the single- and two-photon sectors, we derive equivalent quantum-circuit implementations and unify them within a common Hilbert-space representation using an ancilla-assisted encoding. This construction is then generalized to arbitrary interferometers by replacing each optical beam splitter with a repeating quantum-circuit unit that selectively acts only within the relevant local interference subspace, requiring $N+1$ qubits for a two-photon $N$-mode interferometer and a linear-overhead subspace-identification procedure. Using this framework, gate-based quantum circuit for a four-mode boson-sampling circuit is developed and experimentally implemented on a four-qubit gate-based photonic qubit system. The qubit framework for emulating boson sampling of $n-$photons in $m-$mode will be useful to solve a broad class of sampling complexity problem on a gate-based quantum computers.

quant-pharXiv:2608.09509•Aug 11, 2026•0 citations

Efficient Qubit Simulation of Hybrid Oscillator-Qubit Quantum Computation

Xi Lu, Bojko N. Bakalov, Yuan Liu

We introduce a framework for simulating hybrid oscillator-qubit quantum processors on qubit-only systems through position encoding. By encoding continuous-variable position and momentum wave functions into qubit amplitudes, our method efficiently simulates all Gaussian and conditional Gaussian operations -- encompassing the phase-space instruction set (beam splitter, single-qubit rotation, conditional displacement) and extending to squeezing, conditional squeezing, conditional rotation, and conditional beam splitter -- using $O\!\left(\log^2\!\left(Γ+ \log(1/ε)\right)\right)$ qubit gates per hybrid gate, where $Γ$ is the Fock-level bound and $ε$ is the target precision. This polylogarithmic per-gate complexity represents an exponential improvement over Fock basis encoding approaches, which require exponential quantum or classical resources in the number of qubits per mode. We provide rigorous numerical characterization of quantum Fourier transform errors for Fock-bounded states, enabling precise resource estimation for practical implementations. This work establishes that hybrid oscillator-qubit algorithms can be implemented on qubit processors with polynomial overhead, providing new insights into the computational power trade-offs between discrete-variable and hybrid continuous-discrete-variable quantum computing.

quant-pharXiv:2603.09233•Mar 11, 2026•0 citations

Coherent-State Propagation: A Computational Framework for Simulating Bosonic Quantum Systems

Nikita Guseynov, Zoë Holmes, Armando Angrisani

We introduce coherent-state propagation, a computational framework for simulating bosonic systems. We focus on bosonic circuits composed of displaced linear optics augmented by Kerr nonlinearities, a universal model of bosonic quantum computation that is also physically motivated by driven Bose-Hubbard dynamics. The method works in the Schrödinger picture representing the evolving state as a sparse superposition of coherent states. We develop approximation strategies that keep the simulation cost tractable in physically relevant regimes, notably when the number of Kerr gates is small or the Kerr nonlinearities are weak, and prove rigorous guarantees for both observable estimation and sampling. In particular, bosonic circuits with logarithmically many Kerr gates admit quasi-polynomial-time classical simulation at exponentially small error in trace distance. We further identify a weak-nonlinearity regime in which the runtime is polynomial for arbitrarily small constant precision. We complement these results with numerical benchmarks on the Bose-Hubbard model with all-to-all connectivity. The method reproduces Fock-basis and matrix-product-state reference data, suggesting that it offers a useful route to the classical simulation of bosonic systems.

quant-pharXiv:2604.19625•Apr 22, 2026•0 citations

Exponentially-improved effective descriptions of physical bosonic systems

Varun Upreti, Nicolás Quesada, Ulysse Chabaud

The effective description of a bosonic quantum system identifies the minimum finite dimension required to capture its essential dynamics. This effective dimension plays an important role in the complexity of classical and quantum algorithms for learning and simulating bosonic systems. While generic bosonic states require a dimension scaling as $1/ε^2$ for a precision of approximation $ε$, here we identify a natural energy condition which allows us to improve this scaling exponentially to $\log(1/ε)$. We then prove that most bosonic quantum states satisfy this condition, and in particular those produced by combining Gaussian dynamics with generic energy-preserving dynamics, which include the output states of universal bosonic quantum circuits. We apply this finding to enhance learning algorithms for bosonic quantum states and we further obtain new classical simulation algorithms for a large class of bosonic systems. Finally, using efficient decompositions of Kerr gates as sums of Gaussian gates, we significantly refine these classical simulation algorithms for universal bosonic quantum circuits. Our results demonstrate that physical bosonic systems are significantly more well-behaved than previously assumed, allowing for efficient descriptions even at high precision.

quant-pharXiv:2604.18720•Apr 22, 2026•0 citations

Unified framework for bosonic quantum information encoding, resources, and universality from superselection rules

Éloi Descamps, Astghik Saharyan, Adrien Chivet +2 more

Since the early days of quantum information science, diverse strategies have been proposed to encode qubits and implement quantum gates across various physical platforms. Modes and bosonic states offer a flexible framework for quantum information processing, from single-photon optics to superconducting circuits, combining massless and massive particles. Universal gate sets have been established for specific encodings—such as those based on Fock states or the field’s quadratures. However, a general framework clarifying how different bosonic resources—such as mode structure, particle-number statistics, Gaussian and non-Gaussian states and operations—combine to enable universality and potential quantum advantage irrespectively of the encoding remains lacking. A main difficulty lies in the coexistence of distinct paradigms that appear to combine loosely in bosonic quantum information: discrete-variable schemes—where information is encoded in qubits or effective finite-dimensional systems—and continuous-variable schemes—which rely on infinite-dimensional Hilbert spaces associated with the field’s quadratures. Despite the success of each approach, they remain seemingly incompatible in many aspects—in particular in terms of physical resource analysis for universality—due to their underlying mathematical differences. In this work, we apply a superselection-rule-compliant framework for bosonic quantum information that unifies these paradigms. This representation explicitly incorporates the phase reference and enforces total particle number conservation, allowing for a consistent treatment of bosonic quantum information protocols and resources across both discrete and continuous regimes. We derive necessary and sufficient conditions expressed in terms of Gaussian and non-Gaussian resources for universality that apply to any encoding, without resorting to specific mappings. Our results show that continuous variable encodings are not fundamentally distinct from discrete ones, but rather emerge as subspaces within a general bosonic structure. Importantly, in such a widely studied framework, non-Gaussian properties of bosonic states and operations can be washed out, providing evidence that the continuous variables regime is not the most suitable to fully assess the quantum properties of bosons.

arXiv:2501.03943•Optica Quantum•Feb 6, 2026•2 citations

Logarithmic-depth quantum simulation of boson sampling

Changhun Oh

We show that boson sampling with an arbitrary $m$-mode interferometer and $n\le m$ single-photon inputs can be simulated to inverse-polynomial total-variation error by a logarithmic-depth qubit circuit with polynomially many qubits. The circuit uses Clifford+$T$ gates, arbitrary qubit connectivity, and a single final measurement, and its family is logspace uniform. The key idea is to enlarge the optical system, decompose the resulting transformation into six quadratic shears, and distribute each mode over many submodes. This redistribution permits a fixed local occupation cutoff, after which local basis changes and parallel phase gates give the qubit circuit. Consequently, our result places boson sampling within shallow quantum computation.

quant-pharXiv:2609.18907•Sep 17, 2026•0 citations
We introduce a framework for simulating, on an (n+1)-qubit quantum computer, the action of a Gaussian bosonic (GB) circuit on a state over 2^{n} modes. Specifically, we encode the initial bosonic state's expectation values over quadrature operators (and their covariance matrix) as an input qubit state. This is then evolved by a quantum circuit that effectively implements the symplectic propagators induced by the GB gates. We find families of GB circuits and initial states leading to efficient quantum simulations. For this purpose, we introduce a dictionary that maps between GB and qubit gates such that particle- (non-particle-) preserving GB gates lead to real- (imaginary-) time evolutions at the qubit level. For the special case of particle-preserving circuits, we present a bounded-error-quantum-polynomial time (BQP)-complete GB decision problem, indicating that GB evolutions of Gaussian states on exponentially many modes are as powerful as universal quantum computers. We also perform numerical simulations of an interferometer on ∼8×10^{9} modes, illustrating the power of our framework.
Boson sampling arising from multiphoton interference in linear-optical networks is a prominent non-universal model for quantum computation. Here, by encoding the multi-qubit state to bosonic Fock state, we present a scalable quantum-circuit framework for simulating boson sampling on a universal quantum computing platform. Beginning with balanced beam-splitter transformations on the single- and two-photon sectors, we derive equivalent quantum-circuit implementations and unify them within a common Hilbert-space representation using an ancilla-assisted encoding. This construction is then generalized to arbitrary interferometers by replacing each optical beam splitter with a repeating quantum-circuit unit that selectively acts only within the relevant local interference subspace, requiring N+1N+1N+1 qubits for a two-photon NNN-mode interferometer and a linear-overhead subspace-identification procedure. Using this framework, gate-based quantum circuit for a four-mode boson-sampling circuit is developed and experimentally implemented on a four-qubit gate-based photonic qubit system. The qubit framework for emulating boson sampling of n−n-n−photons in m−m-m−mode will be useful to solve a broad class of sampling complexity problem on a gate-based quantum computers.
We introduce a framework for simulating hybrid oscillator-qubit quantum processors on qubit-only systems through position encoding. By encoding continuous-variable position and momentum wave functions into qubit amplitudes, our method efficiently simulates all Gaussian and conditional Gaussian operations -- encompassing the phase-space instruction set (beam splitter, single-qubit rotation, conditional displacement) and extending to squeezing, conditional squeezing, conditional rotation, and conditional beam splitter -- using O ⁣(log⁡2 ⁣(Γ+log⁡(1/ε)))O\!\left(\log^2\!\left(Γ+ \log(1/ε)\right)\right)O(log2(Γ+log(1/ε))) qubit gates per hybrid gate, where ΓΓΓ is the Fock-level bound and εεε is the target precision. This polylogarithmic per-gate complexity represents an exponential improvement over Fock basis encoding approaches, which require exponential quantum or classical resources in the number of qubits per mode. We provide rigorous numerical characterization of quantum Fourier transform errors for Fock-bounded states, enabling precise resource estimation for practical implementations. This work establishes that hybrid oscillator-qubit algorithms can be implemented on qubit processors with polynomial overhead, providing new insights into the computational power trade-offs between discrete-variable and hybrid continuous-discrete-variable quantum computing.
We introduce coherent-state propagation, a computational framework for simulating bosonic systems. We focus on bosonic circuits composed of displaced linear optics augmented by Kerr nonlinearities, a universal model of bosonic quantum computation that is also physically motivated by driven Bose-Hubbard dynamics. The method works in the Schrödinger picture representing the evolving state as a sparse superposition of coherent states. We develop approximation strategies that keep the simulation cost tractable in physically relevant regimes, notably when the number of Kerr gates is small or the Kerr nonlinearities are weak, and prove rigorous guarantees for both observable estimation and sampling. In particular, bosonic circuits with logarithmically many Kerr gates admit quasi-polynomial-time classical simulation at exponentially small error in trace distance. We further identify a weak-nonlinearity regime in which the runtime is polynomial for arbitrarily small constant precision. We complement these results with numerical benchmarks on the Bose-Hubbard model with all-to-all connectivity. The method reproduces Fock-basis and matrix-product-state reference data, suggesting that it offers a useful route to the classical simulation of bosonic systems.
The effective description of a bosonic quantum system identifies the minimum finite dimension required to capture its essential dynamics. This effective dimension plays an important role in the complexity of classical and quantum algorithms for learning and simulating bosonic systems. While generic bosonic states require a dimension scaling as 1/ε21/ε^21/ε2 for a precision of approximation εεε, here we identify a natural energy condition which allows us to improve this scaling exponentially to log⁡(1/ε)\log(1/ε)log(1/ε). We then prove that most bosonic quantum states satisfy this condition, and in particular those produced by combining Gaussian dynamics with generic energy-preserving dynamics, which include the output states of universal bosonic quantum circuits. We apply this finding to enhance learning algorithms for bosonic quantum states and we further obtain new classical simulation algorithms for a large class of bosonic systems. Finally, using efficient decompositions of Kerr gates as sums of Gaussian gates, we significantly refine these classical simulation algorithms for universal bosonic quantum circuits. Our results demonstrate that physical bosonic systems are significantly more well-behaved than previously assumed, allowing for efficient descriptions even at high precision.
Since the early days of quantum information science, diverse strategies have been proposed to encode qubits and implement quantum gates across various physical platforms. Modes and bosonic states offer a flexible framework for quantum information processing, from single-photon optics to superconducting circuits, combining massless and massive particles. Universal gate sets have been established for specific encodings—such as those based on Fock states or the field’s quadratures. However, a general framework clarifying how different bosonic resources—such as mode structure, particle-number statistics, Gaussian and non-Gaussian states and operations—combine to enable universality and potential quantum advantage irrespectively of the encoding remains lacking. A main difficulty lies in the coexistence of distinct paradigms that appear to combine loosely in bosonic quantum information: discrete-variable schemes—where information is encoded in qubits or effective finite-dimensional systems—and continuous-variable schemes—which rely on infinite-dimensional Hilbert spaces associated with the field’s quadratures. Despite the success of each approach, they remain seemingly incompatible in many aspects—in particular in terms of physical resource analysis for universality—due to their underlying mathematical differences. In this work, we apply a superselection-rule-compliant framework for bosonic quantum information that unifies these paradigms. This representation explicitly incorporates the phase reference and enforces total particle number conservation, allowing for a consistent treatment of bosonic quantum information protocols and resources across both discrete and continuous regimes. We derive necessary and sufficient conditions expressed in terms of Gaussian and non-Gaussian resources for universality that apply to any encoding, without resorting to specific mappings. Our results show that continuous variable encodings are not fundamentally distinct from discrete ones, but rather emerge as subspaces within a general bosonic structure. Importantly, in such a widely studied framework, non-Gaussian properties of bosonic states and operations can be washed out, providing evidence that the continuous variables regime is not the most suitable to fully assess the quantum properties of bosons.
We show that boson sampling with an arbitrary mmm-mode interferometer and n≤mn\le mn≤m single-photon inputs can be simulated to inverse-polynomial total-variation error by a logarithmic-depth qubit circuit with polynomially many qubits. The circuit uses Clifford+TTT gates, arbitrary qubit connectivity, and a single final measurement, and its family is logspace uniform. The key idea is to enlarge the optical system, decompose the resulting transformation into six quadratic shears, and distribute each mode over many submodes. This redistribution permits a fixed local occupation cutoff, after which local basis changes and parallel phase gates give the qubit circuit. Consequently, our result places boson sampling within shallow quantum computation.