Structure-Preserving Integration for Magnetic Gaussian Wave Packet Dynamics
Fuente:
arXiv
Guardado en:
| Autores principales: | , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2026
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866918411198529536 |
|---|---|
| author | Merk, Sebastian Lasser, Caroline |
| author_facet | Merk, Sebastian Lasser, Caroline |
| contents | We develop structure-preserving time integration schemes for Gaussian wave packet dynamics associated with the magnetic Schrödinger equation. The variational Dirac--Frenkel formulation yields a finite-dimensional Hamiltonian system for the wave packet parameters, where the presence of a magnetic vector potential leads to a non-separable structure and a modified symplectic geometry. By introducing kinetic momenta through a minimal substitution, we reformulate the averaged dynamics as a Poisson system that closely parallels the classical equations of charged particle motion. This representation enables the construction of Boris-type integrators adapted to the variational setting. In addition, we propose explicit high-order symplectic schemes based on splitting methods and partitioned Runge--Kutta integrators. The proposed methods conserve the quadratic invariants characterizing the Hagedorn parametrization, preserve linear and angular momentum under symmetry assumptions, and exhibit near-conservation of the averaged Hamiltonian over long time intervals. Rigorous error estimates are derived for both the wave packet parameters and observable quantities, with bounds uniform in the semiclassical parameter. Numerical experiments demonstrate the favorable long-time behavior and structure preservation of the integrators. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2603_25596 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Structure-Preserving Integration for Magnetic Gaussian Wave Packet Dynamics Merk, Sebastian Lasser, Caroline Numerical Analysis Mathematical Physics 65P10, 81Q05, 37M15, 81S30, 78A35 We develop structure-preserving time integration schemes for Gaussian wave packet dynamics associated with the magnetic Schrödinger equation. The variational Dirac--Frenkel formulation yields a finite-dimensional Hamiltonian system for the wave packet parameters, where the presence of a magnetic vector potential leads to a non-separable structure and a modified symplectic geometry. By introducing kinetic momenta through a minimal substitution, we reformulate the averaged dynamics as a Poisson system that closely parallels the classical equations of charged particle motion. This representation enables the construction of Boris-type integrators adapted to the variational setting. In addition, we propose explicit high-order symplectic schemes based on splitting methods and partitioned Runge--Kutta integrators. The proposed methods conserve the quadratic invariants characterizing the Hagedorn parametrization, preserve linear and angular momentum under symmetry assumptions, and exhibit near-conservation of the averaged Hamiltonian over long time intervals. Rigorous error estimates are derived for both the wave packet parameters and observable quantities, with bounds uniform in the semiclassical parameter. Numerical experiments demonstrate the favorable long-time behavior and structure preservation of the integrators. |
| title | Structure-Preserving Integration for Magnetic Gaussian Wave Packet Dynamics |
| topic | Numerical Analysis Mathematical Physics 65P10, 81Q05, 37M15, 81S30, 78A35 |
| url | https://arxiv.org/abs/2603.25596 |