Section: Research Program
Geometric schemes for the Schrödinger equation
Participants : François Castella, Philippe Chartier, Erwan Faou, Florian Méhats.
Schrödinger equation, variational splitting, energy conservation.
Given the Hamiltonian structure of the Schrödinger equation, we are led to consider the question of energy preservation for time-discretization schemes.
At a higher level, the Schrödinger equation is a partial differential equation which may exhibit Hamiltonian structures. This is the case of the time-dependent Schrödinger equation, which we may write as
where is the wave function depending on the spatial variables with (e.g., with or 3 in the partition) and the time . Here, is a (small) positive number representing the scaled Planck constant and is the complex imaginary unit. The Hamiltonian operator is written
with the kinetic and potential energy operators
where is a particle mass and the Laplacian in the variable , and where the real-valued potential acts as a multiplication operator on .
The multiplication by in (8 ) plays the role of the multiplication by in classical mechanics, and the energy is conserved along the solution of (8 ), using the physicists' notations where denotes the Hermitian -product over the phase space. In quantum mechanics, the number of particles is very large making the direct approximation of (8 ) very difficult.
The numerical approximation of (8 ) can be obtained using projections onto submanifolds of the phase space, leading to various PDEs or ODEs: see  ,  for reviews. However the long-time behavior of these approximated solutions is well understood only in this latter case, where the dynamics turns out to be finite dimensional. In the general case, it is very difficult to prove the preservation of qualitative properties of (8 ) such as energy conservation or growth in time of Sobolev norms. The reason for this is that backward error analysis is not directly applicable for PDEs. Overwhelming these difficulties is thus a very interesting challenge.
A particularly interesting case of study is given by symmetric splitting methods, such as the Strang splitting:
where is the time increment (we have set all the parameters to 1 in the equation). As the Laplace operator is unbounded, we cannot apply the standard methods used in ODEs to derive long-time properties of these schemes. However, its projection onto finite dimensional submanifolds (such as Gaussian wave packets space or FEM finite dimensional space of functions in ) may exhibit Hamiltonian or Poisson structure, whose long-time properties turn out to be more tractable.