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

Quantum time‐marching algorithms for solving linear transport problems including boundary conditions

Sergio Bengoechea, Over, Paul, Thomas Rung

Published Apr 9, 2026• 0 citations
Publisher

Abstract

This article presents the first complete application of a quantum time‐marching algorithm for simulating multidimensional linear transport phenomena with arbitrary boundaries, whereby the success probabilities are problem intrinsic. The method adapts the linear combination of unitaries algorithm to block encode the diffusive dynamics, while arbitrary boundary conditions are enforced by the method of images only at the cost of one additional qubit per spatial dimension. As an alternative to the nonperiodic reflection, the direct encoding of Neumann conditions by the unitary decomposition of the discrete time‐marching operator is proposed. All presented algorithms indicate optimal success probabilities while maintaining linear time complexity, thereby securing the practical applicability of the quantum algorithm on fault‐tolerant quantum computers. The proposed time‐marching method is demonstrated through state‐vector simulations of the heat equation in combination with Neumann, Dirichlet, and mixed boundary conditions, showing excellent agreement with classical finite differences.

DOI
10.15480/882.16966

Related papers

Quantum time-marching algorithms for solving linear transport problems including boundary conditions

Sergio Bengoechea, Paul Over, Thomas Rung

This article presents the first complete application of a quantum time-marching algorithm for simulating multidimensional linear transport phenomena with arbitrary boundaries, whereby the success probabilities are problem intrinsic. The method adapts the linear combination of unitaries algorithm to block encode the diffusive dynamics, while arbitrary boundary conditions are enforced by the method of images only at the cost of one additional qubit per spatial dimension. As an alternative to the non-periodic reflection, the direct encoding of Neumann conditions by the unitary decomposition of the discrete time-marching operator is proposed. All presented algorithms indicate optimal success probabilities while maintaining linear time complexity, thereby securing the practical applicability of the quantum algorithm on fault-tolerant quantum computers. The proposed time-marching method is demonstrated through state-vector simulations of the heat equation in combination with Neumann, Dirichlet, and mixed boundary conditions.

quant-pharXiv:2511.04271•International Journal for Numerical Methods in Engineering•Nov 7, 2025•1 citation

Quantum algorithm for the advection-diffusion equation by direct block encoding of the time-marching operator

Paul Over, Sergio Bengoechea, Peter Brearley +2 more

A quantum algorithm for simulating multidimensional scalar transport problems using a time-marching strategy is presented. A direct unitary block encoding of the explicit time-marching operator is constructed, resulting in the intrinsic success probability of the squared solution norm without the need for amplitude amplification, thereby retaining a linear dependence on the simulation time. The algorithm separates the explicit time-marching operator into an advection-like component and a corrective shift operator. The advection-like component is mapped to a Hamiltonian simulation and combined with the shift operator through the linear combination of unitaries algorithm. State-vector simulations of a scalar transported in a steady two-dimensional Taylor-Green vortex support the theoretical findings.

arXiv:2410.07909•Physical Review A•Jun 4, 2025•9 citations

Unconditionally successful quantum Time-Marching algorithm via LCU for nonlinear Burgers equation

Niccolo Fonio, Giuseppe Di Molfetta, Pierre Sagaut

Most recently proposed quantum algorithms for solving linear and nonlinear partial differential equations rely on non-unitary operations. These operations are typically implemented probabilistically, requiring postselection and thus increasing the computational cost. We show that quantum lattice gas algorithms enable unconditionally successful quantum simulation of nonlinearities, yielding, to our knowledge, the first quantum algorithm for Burgers equation whose time steps can be concatenated without probabilistic failure. The key idea is to exploit the correspondence between the stochasticity of quantum measurement in the linear combination of unitaries framework and the intrinsic randomness of the classical lattice gas algorithm. In doing so, we identify general properties that characterize probabilistic classical algorithms amenable to this time-marching formulation, and illustrate the approach with an additional application.

quant-pharXiv:2608.02130•Aug 4, 2026•0 citations

Quantum algorithm for the advection-diffusion equation by direct block encoding of the time-marching operator

Over, Paul, Sergio Bengoechea, Peter Brearley +2 more

A quantum algorithm for simulating multidimensional scalar transport problems using a time-marching strategy is presented. A direct unitary block encoding of the explicit time-marching operator is constructed, resulting in the intrinsic success probability of the squared solution norm without the need for amplitude amplification, thereby retaining a linear dependence on the simulation time. The algorithm separates the explicit time-marching operator into an advection-like component and a corrective shift operator. The advection-like component is mapped to a Hamiltonian simulation and combined with the shift operator through the linear combination of unitaries algorithm. State-vector simulations of a scalar transported in a steady two-dimensional Taylor-Green vortex support the theoretical findings.

•Jul 1, 2025•0 citations

Quantum Implementation of Numerical Methods for Convection-Diffusion Equations: Toward Computational Fluid Dynamics

Bofeng Liu, Lixing Zhu, Zixuan Yang +1 more

We present quantum numerical methods for the typical initial boundary value problems (IBVPs) of convection-diffusion equations in fluid dynamics. The IBVP is discretized into a series of linear systems via finite difference methods and explicit time marching schemes. To solve these discrete systems in quantum computers, we design a series of quantum circuits, including four stages of encoding, amplification, adding source terms, and incorporating boundary conditions. In the encoding stage, the initial condition is encoded in the amplitudes of quantum registers as a state vector to take advantage of quantum algorithms in space complexity. In the following three stages, the discrete differential operators in classical computing are converted into unitary evolutions to satisfy the postulate in quantum systems. The related arithmetic calculations in quantum amplitudes are also realized to sum up the increments from these stages. The proposed quantum algorithm is implemented within the open-source quantum computing framework Qiskit [2]. By simulating one-dimensional transient problems, including the Helmholtz equation, the Burgers’ equation, and Navier-Stokes equations, we demonstrate the capability of quantum computers in fluid dynamics.

Communications in Computational Physics•Mar 12, 2023•16 citations

Quantum algorithm for solving the advection equation using Hamiltonian simulation

Peter Brearley, Sylvain Laizet

A quantum algorithm for solving the advection equation by embedding the discrete time-marching operator into Hamiltonian simulations is presented. One-dimensional advection can be simulated directly since the central finite-difference operator for first-order derivatives is anti-Hermitian. Here this is extended to industrially relevant multidimensional flows with realistic boundary conditions and arbitrary finite-difference stencils. A single copy of the initial quantum state is required and the circuit depth grows linearly with the required number of time steps, the sparsity of the time-marching operator, and the inverse of the allowable error. State-vector simulations of a scalar transported in a two-dimensional channel flow and lid-driven cavity configuration are presented as a proof of concept of the proposed approach. Published by the American Physical Society 2024

arXiv:2312.09784•Physical Review A•Jul 9, 2024•28 citations
This article presents the first complete application of a quantum time‐marching algorithm for simulating multidimensional linear transport phenomena with arbitrary boundaries, whereby the success probabilities are problem intrinsic. The method adapts the linear combination of unitaries algorithm to block encode the diffusive dynamics, while arbitrary boundary conditions are enforced by the method of images only at the cost of one additional qubit per spatial dimension. As an alternative to the nonperiodic reflection, the direct encoding of Neumann conditions by the unitary decomposition of the discrete time‐marching operator is proposed. All presented algorithms indicate optimal success probabilities while maintaining linear time complexity, thereby securing the practical applicability of the quantum algorithm on fault‐tolerant quantum computers. The proposed time‐marching method is demonstrated through state‐vector simulations of the heat equation in combination with Neumann, Dirichlet, and mixed boundary conditions, showing excellent agreement with classical finite differences.
This article presents the first complete application of a quantum time-marching algorithm for simulating multidimensional linear transport phenomena with arbitrary boundaries, whereby the success probabilities are problem intrinsic. The method adapts the linear combination of unitaries algorithm to block encode the diffusive dynamics, while arbitrary boundary conditions are enforced by the method of images only at the cost of one additional qubit per spatial dimension. As an alternative to the non-periodic reflection, the direct encoding of Neumann conditions by the unitary decomposition of the discrete time-marching operator is proposed. All presented algorithms indicate optimal success probabilities while maintaining linear time complexity, thereby securing the practical applicability of the quantum algorithm on fault-tolerant quantum computers. The proposed time-marching method is demonstrated through state-vector simulations of the heat equation in combination with Neumann, Dirichlet, and mixed boundary conditions.
A quantum algorithm for simulating multidimensional scalar transport problems using a time-marching strategy is presented. A direct unitary block encoding of the explicit time-marching operator is constructed, resulting in the intrinsic success probability of the squared solution norm without the need for amplitude amplification, thereby retaining a linear dependence on the simulation time. The algorithm separates the explicit time-marching operator into an advection-like component and a corrective shift operator. The advection-like component is mapped to a Hamiltonian simulation and combined with the shift operator through the linear combination of unitaries algorithm. State-vector simulations of a scalar transported in a steady two-dimensional Taylor-Green vortex support the theoretical findings.
Most recently proposed quantum algorithms for solving linear and nonlinear partial differential equations rely on non-unitary operations. These operations are typically implemented probabilistically, requiring postselection and thus increasing the computational cost. We show that quantum lattice gas algorithms enable unconditionally successful quantum simulation of nonlinearities, yielding, to our knowledge, the first quantum algorithm for Burgers equation whose time steps can be concatenated without probabilistic failure. The key idea is to exploit the correspondence between the stochasticity of quantum measurement in the linear combination of unitaries framework and the intrinsic randomness of the classical lattice gas algorithm. In doing so, we identify general properties that characterize probabilistic classical algorithms amenable to this time-marching formulation, and illustrate the approach with an additional application.
A quantum algorithm for simulating multidimensional scalar transport problems using a time-marching strategy is presented. A direct unitary block encoding of the explicit time-marching operator is constructed, resulting in the intrinsic success probability of the squared solution norm without the need for amplitude amplification, thereby retaining a linear dependence on the simulation time. The algorithm separates the explicit time-marching operator into an advection-like component and a corrective shift operator. The advection-like component is mapped to a Hamiltonian simulation and combined with the shift operator through the linear combination of unitaries algorithm. State-vector simulations of a scalar transported in a steady two-dimensional Taylor-Green vortex support the theoretical findings.
We present quantum numerical methods for the typical initial boundary value problems (IBVPs) of convection-diffusion equations in fluid dynamics. The IBVP is discretized into a series of linear systems via finite difference methods and explicit time marching schemes. To solve these discrete systems in quantum computers, we design a series of quantum circuits, including four stages of encoding, amplification, adding source terms, and incorporating boundary conditions. In the encoding stage, the initial condition is encoded in the amplitudes of quantum registers as a state vector to take advantage of quantum algorithms in space complexity. In the following three stages, the discrete differential operators in classical computing are converted into unitary evolutions to satisfy the postulate in quantum systems. The related arithmetic calculations in quantum amplitudes are also realized to sum up the increments from these stages. The proposed quantum algorithm is implemented within the open-source quantum computing framework Qiskit [2]. By simulating one-dimensional transient problems, including the Helmholtz equation, the Burgers’ equation, and Navier-Stokes equations, we demonstrate the capability of quantum computers in fluid dynamics.
A quantum algorithm for solving the advection equation by embedding the discrete time-marching operator into Hamiltonian simulations is presented. One-dimensional advection can be simulated directly since the central finite-difference operator for first-order derivatives is anti-Hermitian. Here this is extended to industrially relevant multidimensional flows with realistic boundary conditions and arbitrary finite-difference stencils. A single copy of the initial quantum state is required and the circuit depth grows linearly with the required number of time steps, the sparsity of the time-marching operator, and the inverse of the allowable error. State-vector simulations of a scalar transported in a two-dimensional channel flow and lid-driven cavity configuration are presented as a proof of concept of the proposed approach. Published by the American Physical Society 2024