Efficient quantum algorithm for nonlinear reaction-diffusion equations and energy estimation
Nonlinear differential equations exhibit rich phenomena in many fields but are notoriously challenging to solve. Recently, Liu et al. [1] demonstrated the first efficient quantum algorithm for dissipative quadratic differential equations under the condition R < 1, where R measures the ratio of nonlinearity to dissipation using the ℓ_2 norm. Here we develop an efficient quantum algorithm based on [1] for reaction-diffusion equations, a class of nonlinear partial differential equations (PDEs). To achieve this, we improve upon the Carleman linearization approach introduced in [1] to obtain a faster convergence rate under the condition R_D < 1, where R_D measures the ratio of nonlinearity to dissipation using the ℓ_∞ norm. Since R_D is independent of the number of spatial grid points n while R increases with n, the criterion R_D<1 is significantly milder than 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.
READ FULL TEXT