In this paper, we integrate neural networks and Gaussian wave packets to numerically solve the Schr & ouml;dinger equation with a smooth potential near the semi-classical limit. Our focus is not only on accurately obtaining solutions when the non-dimensional Planck's constant, epsilon, is small, but also on constructing an operator that maps initial values to solutions for the Schr & ouml;dinger equation with multiscale properties. Using Gaussian wave packets framework, we first reformulate the Schr & ouml;dinger equation as a system of ordinary differential equations. For a single initial condition, we solve the resulting system using PINNs or MscaleDNNs. Numerical simulations indicate that MscaleDNNs outperform PINNs, improving accuracy by one to two orders of magnitude. When dealing with a set of initial conditions, we adopt an operator-learning approach, such as physics-informed DeepONets. Numerical examples validate the effectiveness of physics-informed DeepONets with Gaussian wave packets in accurately mapping initial conditions to solutions.

Publication:

JOURNAL OF COMPUTATIONAL PHYSICS

http://dx.doi.org/10.1016/j.jcp.2026.114869

Author:Jizu Huang

Academy of Mathematics and Systems Science, SKLMS, Chinese Academy of Sciences, Beijing, 100190, PR China

School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing, 100190, PR Chin

Rukang You

Academy of Mathematics and Systems Science, SKLMS, Chinese Academy of Sciences, Beijing, 100190, PR China

School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing, 100190, PR Chin

Corresponding author.

E-mail address: yourukang@lsec.cc.ac.cn (R. You).

Tao Zhou

Academy of Mathematics and Systems Science, SKLMS, Chinese Academy of Sciences, Beijing, 100190, PR China