Abstract:In this paper, we present an efficient quantum algorithm to simulate nonlinear differential equations with polynomial vector fields of arbitrary (finite) degree on quantum platforms. Ordinary differential equations (ODEs) and partial differential equations (PDEs) arise extensively in science and engineering applications. Examples of ODE models include mechanics of rigid bodies, molecular dynamics, chemical kinetics, and epidemiology. Nonlinear PDEs arise in fluid dynamics, combustion, weather forecasting, structural mechanics, plasma dynamics, and finance to name a few. In practice, it is challenging to simulate such equations on classical computers due to high dimensionality, stiffness arising from multiple spatial/temporal scales, nonlinearities, and chaotic dynamics. Typically, high performance computing is used to mitigate computational challenges and involves approximations for tractability. For sparse n -dimensional linear ODEs, quantum algorithms have been developed which can produce a quantum state proportional to the solution in ) 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 and prove polynomial scaling of runtime on k . We demonstrate the computational framework on a numerical example.

Application of General Orthogonal Polynomials to Fast Simulation of Nonlinear Descriptor Systems Through Piecewise-Linear Approximation

A Kind of Improved Discrete Orthogonal Polynomials in Fourier Domain for Experimental Modal Analysis

Investigation and Application of Feature Extraction Based on the Best Uniform Approximation of Chebyshev′s Polynomial with Nonlinear Partial Orthogonal Analysis

Simulation of nonlinear thermomechanical waves with an empirical low dimensional model

A Novel Method of Polynomial Approximation for Parametric Problems in Power Systems

Adaptive Estimation Of Time-Varying Parameters In Linear Systems

Simultaneous Reconstruction of Phase and Amplitude for Wavefront Measurements Based on Nonlinear Optimization Algorithms

A NOVEL METHOD FOR IDENTIFYING NONLINEAR SYSTEMS USING ORTHOGONAL WAVELET NETWORKS

Nonlinear System Identification Using Orthogonal Wavelet Networks

Iterative Cut-Based PWA Approximation of Multi-Dimensional Nonlinear Systems

Machine Learning Optimized Orthogonal Basis Piecewise Polynomial Approximation

Multidimensional piecewise general orthogonal polynomials operator and its applications to parameter identification of distributed parameter systems

An Adaptive Wavelet Method for Nonlinear Circuit Simulation

Piecewise General Orthogonal Polynomial Operator and Its Properties

Anamnesic Neural Differential Equations with Orthogonal Polynomial Projections

An efficient quantum algorithm for simulating polynomial dynamical systems

Log Orthogonal Functions: Approximation Properties and Applications

Lattice Trajectory Piecewise Linear Method for the Simulation of Diode Circuits

Orthogonal Complete U-System and Its Application in CAGD

Orthogonal polynomial approximation and Extended Dynamic Mode Decomposition in chaos

Parameter Identification of a Class of Nonlinear Systems Via Piecewise Multiple Legendre Polynomials