Efficient Quantum Algorithm for Nonlinear Reaction–Diffusion Equations and Energy Estimation

Abstract Nonlinear differential equations exhibit rich phenomena in many fields but are notoriously challenging to solve. Recently, Liu et al. (in: Proceedings of the National Academy of Sciences 118(35), 2021) demonstrated the first efficient quantum algorithm for dissipative quadrat...

Full description

Bibliographic Details
Main Authors: Liu, Jin-Peng, An, Dong, Fang, Di, Wang, Jiasu, Low, Guang H., Jordan, Stephen
Other Authors: Massachusetts Institute of Technology. Center for Theoretical Physics
Format: Article
Language:English
Published: Springer Berlin Heidelberg 2023
Online Access:https://hdl.handle.net/1721.1/152927
Description
Summary:Abstract Nonlinear differential equations exhibit rich phenomena in many fields but are notoriously challenging to solve. Recently, Liu et al. (in: Proceedings of the National Academy of Sciences 118(35), 2021) demonstrated the first efficient quantum algorithm for dissipative quadratic differential equations under the condition $$R < 1$$ R < 1 , where R measures the ratio of nonlinearity to dissipation using the $$\ell _2$$ ℓ 2 norm. Here we develop an efficient quantum algorithm based on Liu et al. (2021) for reaction–diffusion equations, a class of nonlinear partial differential equations (PDEs). To achieve this, we improve upon the Carleman linearization approach introduced in Liu et al. (2021) to obtain a faster convergence rate under the condition $$R_D < 1$$ R D < 1 , where $$R_D$$ R D measures the ratio of nonlinearity to dissipation using the $$\ell _{\infty }$$ ℓ ∞ norm. Since $$R_D$$ R D is independent of the number of spatial grid points n while R increases with n, the criterion $$R_D<1$$ R D < 1 is significantly milder than $$R<1$$ R < 1 for high-dimensional systems and can stay convergent under grid refinement for approximating PDEs. As applications of our quantum algorithm we consider the Fisher-KPP and Allen-Cahn equations, which have interpretations in classical physics. In particular, we show how to estimate the mean square kinetic energy in the solution by postprocessing the quantum state that encodes it to extract derivative information.