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:2406.15821 · Advances in Applied Mathematics and Mechanics

HAM-Schrödingerisation: A Generic Framework of Quantum Simulation for Any Nonlinear PDEs

Shijun Liao

Published Sep 29, 2025• 1 citation
arXivHTMLPDF

Abstract

Recently, Jin et al. proposed a quantum simulation technique for any linear partial differential equations (PDEs), called Schrödingerisation [1-3]. In this paper, the Schrödingerisation technique for quantum simulation is expanded to any nonlinear PDEs by combining it with the homotopy analysis method (HAM) [4-6]. The HAM can transfer a nonlinear PDE into a series of linear PDEs with guaranteeing convergence of the series. In this way, any nonlinear PDEs can be solved by quantum simulation using a quantum computer. For simplicity, we call the procedure “HAM-Schrödingerisation quantum algorithm”. Quantum computing is a groundbreaking technique. Hopefully, the “HAM-Schrödingerisation quantum algorithm” can open a door to highly efficient simulation of complicated turbulent flows by means of quantum computing in future.

DOI
10.4208/aamm.oa-2024-0242

Related papers

Lindbladian Homotopy Analysis Method to Solve Nonlinear Partial Differential Equations

Eunsik Choi, Jungin E. Kim, Xueling Lu +1 more

Quantum scientific computing is to solve engineering and science problems such as simulation and optimization on quantum computers. Solving ordinary and partial differential equations (PDEs) is essential in simulations. However, existing quantum approaches to solve nonlinear PDEs suffer from the issues of curse of dimensionality and convergence during the linearization process. In this paper, a Lindbladian homotopy analysis method (LHAM) is proposed as a quantum differential equation solver to simulate nonlinear and nonunitary dynamics. The original nonlinear problem is first converted to a recursive sequence of linear PDEs with the homotopy analysis method and reformulated as a higher-dimensional lower block triangular linear homogeneous autonomous system. The solution is then embedded in the density matrix and obtained through the Lindblad dynamics simulation. Compared to other methods such as Carleman linearization and the Koopman-von Neumann approach where the dimension of Hilbert space increases polynomially with the inverse of truncation error, the Hilbert space dimension in LHAM increases only logarithmically. LHAM is demonstrated with nonlinear PDEs including Burgers' equation and reduced magnetohydrodynamics equations.

math.NACross-listedarXiv:2604.18924•Apr 22, 2026•0 citations

Quantum Homotopy Algorithm for Solving Nonlinear PDEs and Flow Problems

Sachin S. Bharadwaj, Balasubramanya Nadiga, Stephan Eidenbenz +1 more

Quantum algorithms to integrate nonlinear PDEs governing flow problems are challenging to discover but critical to enhancing the practical usefulness of quantum computing. We present a near-optimal, robust, and end-to-end quantum algorithm to solve time-dependent, dissipative, nonlinear PDEs. We embed the PDEs in a truncated, high-dimensional linear space on the basis of quantum homotopy analysis. The linearized system is discretized and integrated using finite-difference methods with a compact quantum algorithm. The present approach can adapt its input to the nature of nonlinearity and underlying physics. The complexity estimates improve existing approaches in terms of the time-marching system size, simulation time, accuracy parameters, and post-selection parameters. We provide a general embedding strategy, bounds on stability criteria, accuracy, gate counts, and query complexity. A physically motivated measure of nonlinearity is connected to a parameter similar to the flow Reynolds number $Re_{\textrm{H}}$, whose inverse marks the allowed integration window, for given accuracy and complexity. We illustrate the embedding scheme with numerical simulations of Burgers, Fisher--KPP, damped Kuramoto--Sivashinsky, and real Ginzburg--Landau/Allen--Cahn equations. Together, these examples encompass cubic nonlinearity and saturation, quadratic transport, and fourth-order dissipative stiffness. This work shows the potential of hybrid quantum algorithms for simulating nonlinear problems on near-term and fault-tolerant devices.

quant-pharXiv:2512.21033•Dec 25, 2025•0 citations

Quantum simulation of partial differential equations: Applications and detailed analysis

Shi Jin, Nana Liu, Yue Yu

We study a recently introduced simple method [S. Jin, N. Liu, and Y. Yu, Quantum simulation of partial differential equations via Schr\"odingerisation, arXiv:2212.13969] for solving general linear partial differential equations with quantum simulation. This method converts linear partial differential equations into a Hamiltonian system, using a simple transformation called the warped phase transformation. Here we provide a more-in-depth technical discussion and expand on this approach in a more detailed and pedagogical way. We apply this to examples of partial differential equations, including heat, convection, Fokker-Planck, linear Boltzmann, and Black-Scholes equations. This approach can also be extended to general linear partial differential equations, including the Vlasov-Fokker-Planck equation and the Liouville representation equation for nonlinear ordinary differential equations. Extension to higher-order time derivatives is also possible.

arXiv:2212.14703•Physical Review A•Sep 12, 2023•95 citations

Quantum simulation of discrete linear dynamical systems and simple iterative methods in linear algebra

Shi Jin, Nana Liu

Quantum simulation is capable of simulating certain dynamical systems in continuous time—Schrödinger’s equations being the most direct and well known—more efficiently than classical simulation. Any linear dynamical system can in fact be transformed into a system of Schrödinger’s equations via a method called Schrödingerisation (Jin et al. 2022. Quantum simulation of partial differential equations via Schrödingerisation. ( https://arxiv.org/abs/2212.13969 ) and Jin et al. 2023. Phys. Rev. A 108 , 032603. ( doi:10.1103/PhysRevA.108.032603 )). We show how Schrödingerisation allows quantum simulation to be directly used for the simulation of continuous-time versions of general (explicit) iterative schemes or discrete linear dynamical systems. In particular, we use this new method to solve linear systems of equations and for estimating the maximum eigenvector and eigenvalue of a matrix, respectively. This method is applicable using either discrete-variable quantum systems or on hybrid continuous-variable and discrete-variable quantum systems. This framework provides an interesting alternative to solve linear algebra problems using quantum simulation.

arXiv:2304.02865•Proceedings of the Royal Society A Mathematical Physical and Engineering Sciences•Jun 1, 2024•13 citations

Quantum Circuits for partial differential equations via Schrödingerisation

Junpeng Hu, Shi Jin, Nana Liu +1 more

Quantum computing has emerged as a promising avenue for achieving significant speedup, particularly in large-scale PDE simulations, compared to classical computing. One of the main quantum approaches involves utilizing Hamiltonian simulation, which is directly applicable only to Schrödinger-type equations. To address this limitation, Schrödingerisation techniques have been developed, employing the warped transformation to convert general linear PDEs into Schrödinger-type equations. However, despite the development of Schrödingerisation techniques, the explicit implementation of the corresponding quantum circuit for solving general PDEs remains to be designed. In this paper, we present detailed implementation of a quantum algorithm for general PDEs using Schrödingerisation techniques. We provide examples of the heat equation, and the advection equation approximated by the upwind scheme, to demonstrate the effectiveness of our approach. Complexity analysis is also carried out to demonstrate the quantum advantages of these algorithms in high dimensions over their classical counterparts.

arXiv:2403.10032•Quantum•Dec 12, 2024•27 citations

Variational Quantum Homotopy Perturbation Method to Solve Nonlinear Partial Differential Equations

Jungin E. Kim, Eunsik Choi, Yan Wang

Solving nonlinear partial differential equations is important in various scientific and engineering applications. Recently, quantum computing was introduced as an alternative computational paradigm for solving nonlinear partial differential equations. In this paper, a new method called the variational quantum homotopy perturbation method is proposed to improve the scalability of solving nonlinear partial differential equations through two aspects. First, the dimension of the Hilbert space remains the same after the nonlinear partial differential equation is linearized through the homotopy perturbation. Second, the solutions are obtained with a variational quantum simulation framework, where the number of qubits is decreased with functional encoding and the depth of parametrized circuits is reduced. The additional contribution of this paper is the introduction of new criteria for selecting the homotopy series truncation order and circuit depth for cost-effective variational quantum homotopy perturbation method. The proposed approach is demonstrated with several examples, including the vorticity transport equation and the reduced magnetohydrodynamics equations.

quant-pharXiv:2609.06357•Sep 8, 2026•0 citations
Recently, Jin et al. proposed a quantum simulation technique for any linear partial differential equations (PDEs), called Schrödingerisation [1-3]. In this paper, the Schrödingerisation technique for quantum simulation is expanded to any nonlinear PDEs by combining it with the homotopy analysis method (HAM) [4-6]. The HAM can transfer a nonlinear PDE into a series of linear PDEs with guaranteeing convergence of the series. In this way, any nonlinear PDEs can be solved by quantum simulation using a quantum computer. For simplicity, we call the procedure “HAM-Schrödingerisation quantum algorithm”. Quantum computing is a groundbreaking technique. Hopefully, the “HAM-Schrödingerisation quantum algorithm” can open a door to highly efficient simulation of complicated turbulent flows by means of quantum computing in future.
Quantum scientific computing is to solve engineering and science problems such as simulation and optimization on quantum computers. Solving ordinary and partial differential equations (PDEs) is essential in simulations. However, existing quantum approaches to solve nonlinear PDEs suffer from the issues of curse of dimensionality and convergence during the linearization process. In this paper, a Lindbladian homotopy analysis method (LHAM) is proposed as a quantum differential equation solver to simulate nonlinear and nonunitary dynamics. The original nonlinear problem is first converted to a recursive sequence of linear PDEs with the homotopy analysis method and reformulated as a higher-dimensional lower block triangular linear homogeneous autonomous system. The solution is then embedded in the density matrix and obtained through the Lindblad dynamics simulation. Compared to other methods such as Carleman linearization and the Koopman-von Neumann approach where the dimension of Hilbert space increases polynomially with the inverse of truncation error, the Hilbert space dimension in LHAM increases only logarithmically. LHAM is demonstrated with nonlinear PDEs including Burgers' equation and reduced magnetohydrodynamics equations.
Quantum algorithms to integrate nonlinear PDEs governing flow problems are challenging to discover but critical to enhancing the practical usefulness of quantum computing. We present a near-optimal, robust, and end-to-end quantum algorithm to solve time-dependent, dissipative, nonlinear PDEs. We embed the PDEs in a truncated, high-dimensional linear space on the basis of quantum homotopy analysis. The linearized system is discretized and integrated using finite-difference methods with a compact quantum algorithm. The present approach can adapt its input to the nature of nonlinearity and underlying physics. The complexity estimates improve existing approaches in terms of the time-marching system size, simulation time, accuracy parameters, and post-selection parameters. We provide a general embedding strategy, bounds on stability criteria, accuracy, gate counts, and query complexity. A physically motivated measure of nonlinearity is connected to a parameter similar to the flow Reynolds number ReHRe_{\textrm{H}}ReH​, whose inverse marks the allowed integration window, for given accuracy and complexity. We illustrate the embedding scheme with numerical simulations of Burgers, Fisher--KPP, damped Kuramoto--Sivashinsky, and real Ginzburg--Landau/Allen--Cahn equations. Together, these examples encompass cubic nonlinearity and saturation, quadratic transport, and fourth-order dissipative stiffness. This work shows the potential of hybrid quantum algorithms for simulating nonlinear problems on near-term and fault-tolerant devices.
We study a recently introduced simple method [S. Jin, N. Liu, and Y. Yu, Quantum simulation of partial differential equations via Schr\"odingerisation, arXiv:2212.13969] for solving general linear partial differential equations with quantum simulation. This method converts linear partial differential equations into a Hamiltonian system, using a simple transformation called the warped phase transformation. Here we provide a more-in-depth technical discussion and expand on this approach in a more detailed and pedagogical way. We apply this to examples of partial differential equations, including heat, convection, Fokker-Planck, linear Boltzmann, and Black-Scholes equations. This approach can also be extended to general linear partial differential equations, including the Vlasov-Fokker-Planck equation and the Liouville representation equation for nonlinear ordinary differential equations. Extension to higher-order time derivatives is also possible.
Quantum simulation is capable of simulating certain dynamical systems in continuous time—Schrödinger’s equations being the most direct and well known—more efficiently than classical simulation. Any linear dynamical system can in fact be transformed into a system of Schrödinger’s equations via a method called Schrödingerisation (Jin et al. 2022. Quantum simulation of partial differential equations via Schrödingerisation. ( https://arxiv.org/abs/2212.13969 ) and Jin et al. 2023. Phys. Rev. A 108 , 032603. ( doi:10.1103/PhysRevA.108.032603 )). We show how Schrödingerisation allows quantum simulation to be directly used for the simulation of continuous-time versions of general (explicit) iterative schemes or discrete linear dynamical systems. In particular, we use this new method to solve linear systems of equations and for estimating the maximum eigenvector and eigenvalue of a matrix, respectively. This method is applicable using either discrete-variable quantum systems or on hybrid continuous-variable and discrete-variable quantum systems. This framework provides an interesting alternative to solve linear algebra problems using quantum simulation.
Quantum computing has emerged as a promising avenue for achieving significant speedup, particularly in large-scale PDE simulations, compared to classical computing. One of the main quantum approaches involves utilizing Hamiltonian simulation, which is directly applicable only to Schrödinger-type equations. To address this limitation, Schrödingerisation techniques have been developed, employing the warped transformation to convert general linear PDEs into Schrödinger-type equations. However, despite the development of Schrödingerisation techniques, the explicit implementation of the corresponding quantum circuit for solving general PDEs remains to be designed. In this paper, we present detailed implementation of a quantum algorithm for general PDEs using Schrödingerisation techniques. We provide examples of the heat equation, and the advection equation approximated by the upwind scheme, to demonstrate the effectiveness of our approach. Complexity analysis is also carried out to demonstrate the quantum advantages of these algorithms in high dimensions over their classical counterparts.
Solving nonlinear partial differential equations is important in various scientific and engineering applications. Recently, quantum computing was introduced as an alternative computational paradigm for solving nonlinear partial differential equations. In this paper, a new method called the variational quantum homotopy perturbation method is proposed to improve the scalability of solving nonlinear partial differential equations through two aspects. First, the dimension of the Hilbert space remains the same after the nonlinear partial differential equation is linearized through the homotopy perturbation. Second, the solutions are obtained with a variational quantum simulation framework, where the number of qubits is decreased with functional encoding and the depth of parametrized circuits is reduced. The additional contribution of this paper is the introduction of new criteria for selecting the homotopy series truncation order and circuit depth for cost-effective variational quantum homotopy perturbation method. The proposed approach is demonstrated with several examples, including the vorticity transport equation and the reduced magnetohydrodynamics equations.