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.
Type de document :
Article dans une revue
SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2009, 31 (4), pp.3027-3041
Liste complète des métadonnées

https://hal.inria.fr/hal-00777664
Contributeur : Erwan Faou <>
Soumis le : jeudi 17 janvier 2013 - 17:21:23
Dernière modification le : mardi 19 juin 2018 - 11:12:07

Identifiants

  • HAL Id : hal-00777664, version 1

Citation

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〉

Partager

Métriques

Consultations de la notice

395