Uniformly accurate time-splitting methods for the semiclassical linear Schr{\"o}dinger equation

This article is devoted to the construction of numerical methods which remain insensitive to the smallness of the semiclassical parameter for the linear Schr{\"o}dinger equation in the semiclassical limit. We specifically analyse the convergence behavior of the first-order splitting. Our main result is a proof of uniform accuracy. We illustrate the properties of our methods with simulations

AN ASYMPTOTIC PRESERVING SCHEME BASED ON A NEW FORMULATION FOR NLS IN THE SEMICLASSICAL LIMIT

Abstract. We consider the semiclassical limit for the nonlinear Schrödinger equation. We introduce a phase/amplitude representation given by a system similar to the hydrodynamical formulation, whose novelty consists in including some asymptotically vanishing viscosity. We prove that the system is always locally well-posed in a class of Sobolev spaces, and globally well-posed for a fixed positive Planck constant in the one-dimensional case. We propose a second order numerical scheme which is asymptotic preserving. Before singularities appear in the limiting Euler equation, we recover the quadratic physical observables as well as the wave function with mesh size and time step independent of the Planck constant. This approach is also well suited to the linear Schrödinger equation