Computing semi-classical quantum dynamics with Hagedorn wavepackets - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue SIAM Journal on Scientific Computing Année : 2009

Computing semi-classical quantum dynamics with Hagedorn wavepackets

Résumé

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.
Fichier non déposé

Dates et versions

hal-00777664 , version 1 (17-01-2013)

Identifiants

  • HAL Id : hal-00777664 , version 1

Citer

Erwan Faou, Christian Lubich, Vasile Gradinaru. Computing semi-classical quantum dynamics with Hagedorn wavepackets. SIAM Journal on Scientific Computing, 2009, 31 (4), pp.3027-3041. ⟨hal-00777664⟩
184 Consultations
0 Téléchargements

Partager

Gmail Facebook X LinkedIn More