Section: New Results
Analysis and numerical simulation of the Schrödinger equation
The linear or nonlinear Schrödinger equation with potential is one of the basic equations of quantum mechanics and it arises in many areas of physical and technological interest, e.g. in quantum semiconductors, in electromagnetic wave propagation, and in seismic migration. The Schrödinger equation is the lowest order one-way approximation (paraxial wave equation) to the Helmholtz equation and is called Fresnel equation in optics, or standard parabolic equation in underwater acoustics. The solution of the equation is defined on an unbounded domain. If one wants to solve such a whole space evolution problem numerically, one has to restrict the computational domain by introducing artificial boundary conditions. So, the objective is to approximate the exact solution of the whole-space problem, restricted to a finite computational domain. A review article  was written this year to describe and compare the different current approaches of constructing and discretizing the transparent boundary conditions in one and two dimensions. However, these approaches are limited to the linear case (or nonlinear with the classical cubic nonlinearity: an article written was dedicated to this case this year  ) and constant potentials. Therefore, in collaboration with X. Antoine (IECN Nancy and INRIA Lorraine), we proposed to P. Klein to study, in her PhD thesis, the case of the Schrödinger equation with variable potentials. The study of the non-stationary one-dimensional case has already led to one publication  and some preliminary results in the stationary case are really promising. These cases are relevant since for example the equations appear in the Bose Einstein condensate with a quadratic potential.
This problem is obviously not limited to the Schrödinger equation and new developments are in progress on the Korteweg de Vries equation with M. Ehrhardt. This equation is more difficult to study due to its third order derivative in space.
The publications of Guillaume Dujardin for the last two years deal with the long time analysis of exact and numerical solutions of Schrodinger equations. In  , G.D. shows that a class of high order exponential integrators that are of common use for parabolic equations performs well in finite time for linear and nonlinear Schrodinger equations in periodic d -dimensional domains. In particular, G.D. provides sufficient conditions to achieve high orders with such methods applied to Schrodinger equations. In  , G.D. shows with F. Castella that the Lie-Trotter splitting method for the computation of numerical solutions of linear Schrodinger equations in a fully discrete setting preserves the regularity of the numerical solution over long times. In particular, the bounds on the regularity of the numerical solutions do not depend on the spatial discretization parameter. This completes a previous result obtained by G.D and E. Faou in a semi-discrete setting. In  , G.D. studies the long time asymptotics of the solutions of linear Schrodinger equations considered as initial-boundary value problems on the half-line and on bounded intervals when the boundary data are periodic functions of time. G.D. obtains theoretical results using a transformation method introduced by T. Fokas and provides several numerical experiments to support them.