On Fourier time-splitting methods for nonlinear Schrodinger equations in the semi-classical limit
Provides rigorous numerical analysis for a known method in a specific asymptotic regime, benefiting computational quantum mechanics.
The paper proves error estimates for Lie-Trotter splitting methods applied to nonlinear Schrödinger equations in the semi-classical limit, showing that the time step can be independent of the Planck constant for quadratic observables.
We prove an error estimate for a Lie-Trotter splitting operator associated to the Schrodinger-Poisson equation in the semiclassical regime, when the WKB approximation is valid. In finite time, and so long as the solution to a compressible Euler-Poisson equation is smooth, the error between the numerical solution and the exact solution is controlled in Sobolev spaces, in a suitable phase/amplitude representation. As a corollary, we infer the numerical convergence of the quadratic observables with a time step independent of the Planck constant. A similar result is established for the nonlinear Schrodinger equation in the weakly nonlinear regime.