An efficient quantum algorithm for simulating polynomial differential equations

We present an efficient quantum algorithm to simulate nonlinear differential equations with polynomial vector fields of arbitrary degree on quantum platforms. Models of physical systems that are governed by ordinary differential equations (ODEs) or partial differential equation (PDEs) can be challenging to solve on classical computers due to high dimensionality, stiffness, nonlinearities, and sensitive dependence to initial conditions. For sparse $n$-dimensional linear ODEs, quantum algorithms have been developed which can produce a quantum state proportional to the solution in poly(log(nx)) time using the quantum linear systems algorithm (QLSA). Recently, this framework was extended to systems of nonlinear ODEs with quadratic polynomial vector fields by applying Carleman linearization that enables the embedding of the quadratic system into an approximate linear form. A detailed complexity analysis was conducted which showed significant computational advantage under certain conditions. We present an extension of this algorithm to deal with systems of nonlinear ODEs with $k$-th degree polynomial vector fields for arbitrary (finite) values of $k$. The steps involve: 1) mapping the $k$-th degree polynomial ODE to a higher dimensional quadratic polynomial ODE; 2) applying Carleman linearization to transform the quadratic ODE to an infinite-dimensional system of linear ODEs; 3) truncating and discretizing the linear ODE and solving using the forward Euler method and QLSA. Alternatively, one could apply Carleman linearization directly to the $k$-th degree polynomial ODE, resulting in a system of infinite-dimensional linear ODEs, and then apply step 3. This solution route can be computationally more efficient. We present detailed complexity analysis of the proposed algorithms, prove polynomial scaling of runtime on $k$ and demonstrate the framework on an example.

[1]  Shi Jin,et al.  Time complexity analysis of quantum algorithms via linear representations for nonlinear ordinary and partial differential equations , 2022, J. Comput. Phys..

[2]  A. Amini,et al.  Carleman Linearization of Nonlinear Systems and Its Finite-Section Approximations , 2022, ArXiv.

[3]  R. Lowrie,et al.  Koopman von Neumann mechanics and the Koopman representation: A perspective on solving nonlinear dynamical systems with quantum computers , 2022, ArXiv.

[4]  H. Krovi Improved quantum algorithms for linear and nonlinear differential equations , 2022, Quantum.

[5]  A. Ourmazd,et al.  Embedding classical dynamics in a quantum computer , 2020, Physical Review A.

[6]  Andrew M. Childs,et al.  Efficient quantum algorithm for dissipative nonlinear differential equations , 2020, Proceedings of the National Academy of Sciences.

[7]  N. Linden,et al.  Quantum vs. Classical Algorithms for Solving the Heat Equation , 2020, Communications in Mathematical Physics.

[8]  I. Joseph Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics , 2020, 2003.09980.

[9]  Andrew M. Childs,et al.  High-precision quantum algorithms for partial differential equations , 2020, Quantum.

[10]  Andrew M. Childs,et al.  Quantum Spectral Methods for Differential Equations , 2019, Communications in Mathematical Physics.

[11]  Pedro C. S. Costa,et al.  Quantum algorithm for simulating the wave equation , 2017, Physical Review A.

[12]  Amaury Pouly,et al.  Explicit Error Bounds for Carleman Linearization , 2017, ArXiv.

[13]  Andrew M. Childs,et al.  Quantum Algorithm for Linear Differential Equations with Exponentially Improved Dependence on Precision , 2017, Communications in Mathematical Physics.

[14]  I. Bose,et al.  Allee dynamics: growth, extinction and range expansion , 2017, bioRxiv.

[15]  Ashley Montanaro,et al.  Quantum algorithms and the finite element method , 2015, 1512.05903.

[16]  Andrew M. Childs,et al.  Quantum Algorithm for Systems of Linear Equations with Exponentially Improved Dependence on Precision , 2015, SIAM J. Comput..

[17]  Vitaly Volpert,et al.  Elliptic Partial Differential Equations: Volume 2: Reaction-Diffusion Equations , 2014 .

[18]  Dominic W. Berry,et al.  High-order quantum algorithm for solving linear differential equations , 2010, ArXiv.

[19]  Tobias J. Osborne,et al.  A quantum algorithm to solve nonlinear differential equations , 2008, 0812.4423.

[20]  Thierry Paul,et al.  Quantum computation and quantum information , 2007, Mathematical Structures in Computer Science.

[21]  Robert Kersner,et al.  Travelling Waves in Nonlinear Diffusion-Convection Reaction , 2004 .

[22]  Charles R. Johnson,et al.  Topics in Matrix Analysis , 1991 .

[23]  W. Steeb,et al.  Nonlinear dynamical systems and Carleman linearization , 1991 .

[24]  L. Nirenberg,et al.  On elliptic partial differential equations , 1959 .