Section: New Results

Computing semi-classical quantum dynamics with Hagedorn wavepackets

Participant : Erwan Faou.

In Ref. [27] , we consider the approximation of multi-particle quantum dynamics in the semiclassical regime by Hagedorn wavepackets, which are products of complex Gaussians with polynomials that form an orthonormal L2 basis and preserve their type under propagation in Schrödinger equations with quadratic potentials. We build a time-reversible, fully explicit time-stepping algorithm to approximate the solution of the Hagedorn wavepacket dynamics. The algorithm is based on a splitting between the kinetic and potential part of the Hamiltonian operator, as well as on a splitting of the potential into its local quadratic approximation and the remainder. The algorithm is robust in the semi-classical limit. It reduces to the Strang splitting of the Schrödinger equation in the limit of the full basis set, and it advances positions and momenta by the Störmer-Verlet method for the classical equations of motion. The algorithm allows for the treatment of multi-particle problems by thinning out the basis according to a hyperbolic cross approximation, and of high-dimensional problems by Hartree-type approximations in a moving coordinate frame.