Skip to Main content Skip to Navigation
Journal articles

Computing semi-classical quantum dynamics with Hagedorn wavepackets

Abstract : We consider the approximation of multi-particle quantum dynamics in the semi-classical 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.
Document type :
Journal articles
Complete list of metadatas
Contributor : Erwan Faou <>
Submitted on : Thursday, January 17, 2013 - 5:21:23 PM
Last modification on : Monday, August 17, 2020 - 5:06:01 PM


  • HAL Id : hal-00777664, version 1


Erwan Faou, Christian Lubich, Vasile Gradinaru. Computing semi-classical quantum dynamics with Hagedorn wavepackets. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2009, 31 (4), pp.3027-3041. ⟨hal-00777664⟩



Record views