English

Neural Q-learning for solving PDEs

Numerical Analysis 2023-06-27 v2 Machine Learning Numerical Analysis Analysis of PDEs Probability Machine Learning

Abstract

Solving high-dimensional partial differential equations (PDEs) is a major challenge in scientific computing. We develop a new numerical method for solving elliptic-type PDEs by adapting the Q-learning algorithm in reinforcement learning. Our "Q-PDE" algorithm is mesh-free and therefore has the potential to overcome the curse of dimensionality. Using a neural tangent kernel (NTK) approach, we prove that the neural network approximator for the PDE solution, trained with the Q-PDE algorithm, converges to the trajectory of an infinite-dimensional ordinary differential equation (ODE) as the number of hidden units \rightarrow \infty. For monotone PDE (i.e. those given by monotone operators, which may be nonlinear), despite the lack of a spectral gap in the NTK, we then prove that the limit neural network, which satisfies the infinite-dimensional ODE, converges in L2L^2 to the PDE solution as the training time \rightarrow \infty. More generally, we can prove that any fixed point of the wide-network limit for the Q-PDE algorithm is a solution of the PDE (not necessarily under the monotone condition). The numerical performance of the Q-PDE algorithm is studied for several elliptic PDEs.

Keywords

Cite

@article{arxiv.2203.17128,
  title  = {Neural Q-learning for solving PDEs},
  author = {Samuel N. Cohen and Deqing Jiang and Justin Sirignano},
  journal= {arXiv preprint arXiv:2203.17128},
  year   = {2023}
}
R2 v1 2026-06-24T10:33:32.425Z