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:2608.02130

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

Niccolo Fonio, Giuseppe Di Molfetta, Pierre Sagaut

quant-phSubmitted Aug 3, 2026• Announced Aug 4, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

Related papers

Quantum algorithm for nonlinear Burgers' equation for high-speed compressible flows

Esmaeil Esmaeilifar, Doyeol Ahn, Rho Shin Myong

Recent advances in quantum hardware and quantum computing algorithms promise significant breakthroughs in computational capabilities. Quantum computers can achieve exponential improvements in speed vs classical computers by employing principles of quantum mechanics like superposition and entanglement. However, designing quantum algorithms to solve the nonlinear partial differential equations governing fluid dynamics is challenging due to the inherent linearity of quantum mechanics, which requires unitary transformation. In this study, we first address in detail several challenges that arise when trying to deal with nonlinearity using quantum algorithms and then propose a novel pure quantum algorithm for solving a nonlinear Burgers' equation. We employed multiple copies of the state vector to calculate the nonlinear term, which is necessary due to the no-cloning theorem. By reusing qubits from the previous time steps, we significantly reduced the number of qubits required for multi-step simulations, from exponential/quadratic scaling in earlier studies to linear scaling in time in the current study. We also employed various advanced quantum techniques, including block-encoding, quantum Hadamard product, and the linear combination of unitaries, to design a quantum circuit for the proposed quantum algorithm. The quantum circuit was executed on quantum simulators, and the obtained results demonstrated excellent agreement with those from classical simulations.

Physics of Fluids•Oct 1, 2024•9 citations

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

Efficient Quantum Simulation for Nonlinear Stochastic Differential Equations

Xiangyu Li, Ahmet Burak Catli, Ho Kiat Lim +4 more

Nonlinear stochastic differential equations (NSDEs) are a pillar of mathematical modeling for scientific and engineering applications. Accurate and efficient simulation of large-scale NSDEs is prohibitive on classical computers due to the large number of degrees of freedom, and it is challenging on quantum computers due to the linear and unitary nature of quantum mechanics. We develop a quantum algorithm to tackle nonlinear differential equations driven by the Ornstein-Uhlenbeck (OU) stochastic process. The query complexity of our algorithm scales logarithmically with the error tolerance and nearly quadratically with the simulation time. Our algorithmic framework comprises probabilistic Carleman linearization (PCL) to tackle nonlinearity coupled with stochasticity, and stochastic linear combination of Hamiltonian simulations (SLCHS) to simulate stochastic non-unitary dynamics. We obtain probabilistic exponential convergence for the Carleman linearization of Liu et al. [1], provided the NSDE is stable and reaches a steady state. We extend deterministic LCHS to stochastic linear differential equations, retaining near-optimal parameter scaling from An et al. [2] except for the nearly quadratic time scaling. This is achieved by using Monte Carlo integration for time discretization of both the stochastic inhomogeneous term in LCHS and the truncated Dyson series for each Hamiltonian simulation.

quant-pharXiv:2603.12398•Mar 16, 2026•0 citations

A Tensor Network-Based Quantum Algorithm for the Nonlinear 1D Burgers' Equation

Muralikrishnan Gopalakrishnan Meena, Vincent Jones, Yu Zhang +1 more

In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.

AIAA SciTech Forum•Jan 8, 2026•0 citations

Quantum algorithm for differential equations via permutation matrix representation with application to the Burgers equation

Hriday Sabharwal, Amir Kalev, Itay Hen

We develop a quantum algorithm for solving the dynamics of the nonlinear viscous Burgers equation. We apply the Carleman linearization procedure on the spatially discretized equation, followed by a padding scheme that allows implementation on qubit registers. Existing Carleman-based quantum algorithms commonly formulate the lifted linear differential equation in an oracle model. Here we decompose the padded generator into diagonal masks and reversible arithmetic permutations using the Permutation Matrix Representation (PMR), which we show to be naturally compatible with the Linear Combination of Hamiltonian Simulations (LCHS) algorithm. Under the assumptions required by LCHS - most importantly positive semidefiniteness of the Hermitian part of the linear generator, possibly after a stabilizing shift - the algorithm prepares a normalized quantum state proportional to the solution of the truncated lifted system; the stabilizing shift introduces an exponential postselection overhead, which we quantify and mitigate through a rescaling scheme. We show that our algorithm scales with the off-diagonal norm of the Carleman generator instead of the matrix norm, which can be advantageous for other generators that are diagonally dominant. We also extend the PMR scheme to general fluid equations that may contain higher-order derivatives or nonlinear terms, or may involve multiple fluid variables or spatial dimensions. The construction illustrates how PMR can serve as a convenient Hamiltonian-simulation primitive for a broader class of LCU-based algorithms.

quant-pharXiv:2608.19508•Aug 21, 2026•0 citations

Quantum Lattice-Gas Model for the Burgers Equation. Revision 1.4

Jeffrey Yepez

A quantum algorithm is presented for modeling the time evolution of a continuous field governed by the nonlinear Burgers equation in one spatial dimension. It is a microscopic-scale algorithm for a type-II quantum computer, a large lattice of small quantum computers interconnected in nearest neighbor fashion by classical communication channels. A formula for quantum state preparation is presented. The unitary evolution is governed by a conservative quantum gate applied to each node of the lattice independently. Following each quantum gate operation, ensemble measurements over independent microscopic realizations are made resulting in a finite-difference Boltzmann equation at the mesoscopic scale. The measured values are then used to re-prepare the quantum state and one time step is completed. The procedure of state preparation, quantum gate application, and ensemble measurement is continued ad infinitum. The Burgers equation is derived as an effective field theory governing the behavior of the quantum computer at its macroscopic scale where both the lattice cell site and the time step interval become infinitesimal. A numerical simulation of shock formation is carried out and agrees with the exact analytical solution.

•Apr 1, 2002•0 citations
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.
Recent advances in quantum hardware and quantum computing algorithms promise significant breakthroughs in computational capabilities. Quantum computers can achieve exponential improvements in speed vs classical computers by employing principles of quantum mechanics like superposition and entanglement. However, designing quantum algorithms to solve the nonlinear partial differential equations governing fluid dynamics is challenging due to the inherent linearity of quantum mechanics, which requires unitary transformation. In this study, we first address in detail several challenges that arise when trying to deal with nonlinearity using quantum algorithms and then propose a novel pure quantum algorithm for solving a nonlinear Burgers' equation. We employed multiple copies of the state vector to calculate the nonlinear term, which is necessary due to the no-cloning theorem. By reusing qubits from the previous time steps, we significantly reduced the number of qubits required for multi-step simulations, from exponential/quadratic scaling in earlier studies to linear scaling in time in the current study. We also employed various advanced quantum techniques, including block-encoding, quantum Hadamard product, and the linear combination of unitaries, to design a quantum circuit for the proposed quantum algorithm. The quantum circuit was executed on quantum simulators, and the obtained results demonstrated excellent agreement with those from classical simulations.
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.
Nonlinear stochastic differential equations (NSDEs) are a pillar of mathematical modeling for scientific and engineering applications. Accurate and efficient simulation of large-scale NSDEs is prohibitive on classical computers due to the large number of degrees of freedom, and it is challenging on quantum computers due to the linear and unitary nature of quantum mechanics. We develop a quantum algorithm to tackle nonlinear differential equations driven by the Ornstein-Uhlenbeck (OU) stochastic process. The query complexity of our algorithm scales logarithmically with the error tolerance and nearly quadratically with the simulation time. Our algorithmic framework comprises probabilistic Carleman linearization (PCL) to tackle nonlinearity coupled with stochasticity, and stochastic linear combination of Hamiltonian simulations (SLCHS) to simulate stochastic non-unitary dynamics. We obtain probabilistic exponential convergence for the Carleman linearization of Liu et al. [1], provided the NSDE is stable and reaches a steady state. We extend deterministic LCHS to stochastic linear differential equations, retaining near-optimal parameter scaling from An et al. [2] except for the nearly quadratic time scaling. This is achieved by using Monte Carlo integration for time discretization of both the stochastic inhomogeneous term in LCHS and the truncated Dyson series for each Hamiltonian simulation.
In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.
We develop a quantum algorithm for solving the dynamics of the nonlinear viscous Burgers equation. We apply the Carleman linearization procedure on the spatially discretized equation, followed by a padding scheme that allows implementation on qubit registers. Existing Carleman-based quantum algorithms commonly formulate the lifted linear differential equation in an oracle model. Here we decompose the padded generator into diagonal masks and reversible arithmetic permutations using the Permutation Matrix Representation (PMR), which we show to be naturally compatible with the Linear Combination of Hamiltonian Simulations (LCHS) algorithm. Under the assumptions required by LCHS - most importantly positive semidefiniteness of the Hermitian part of the linear generator, possibly after a stabilizing shift - the algorithm prepares a normalized quantum state proportional to the solution of the truncated lifted system; the stabilizing shift introduces an exponential postselection overhead, which we quantify and mitigate through a rescaling scheme. We show that our algorithm scales with the off-diagonal norm of the Carleman generator instead of the matrix norm, which can be advantageous for other generators that are diagonally dominant. We also extend the PMR scheme to general fluid equations that may contain higher-order derivatives or nonlinear terms, or may involve multiple fluid variables or spatial dimensions. The construction illustrates how PMR can serve as a convenient Hamiltonian-simulation primitive for a broader class of LCU-based algorithms.
A quantum algorithm is presented for modeling the time evolution of a continuous field governed by the nonlinear Burgers equation in one spatial dimension. It is a microscopic-scale algorithm for a type-II quantum computer, a large lattice of small quantum computers interconnected in nearest neighbor fashion by classical communication channels. A formula for quantum state preparation is presented. The unitary evolution is governed by a conservative quantum gate applied to each node of the lattice independently. Following each quantum gate operation, ensemble measurements over independent microscopic realizations are made resulting in a finite-difference Boltzmann equation at the mesoscopic scale. The measured values are then used to re-prepare the quantum state and one time step is completed. The procedure of state preparation, quantum gate application, and ensemble measurement is continued ad infinitum. The Burgers equation is derived as an effective field theory governing the behavior of the quantum computer at its macroscopic scale where both the lattice cell site and the time step interval become infinitesimal. A numerical simulation of shock formation is carried out and agrees with the exact analytical solution.