Extended phase-space symplectic integration for electron dynamics

Abstract

We investigate the use of extended phase-space symplectic integration for simulating two different classes of electron dynamics. The first one, with one and a half degrees of freedom, comes from plasma physics and describes the classical dynamics of a charged particle in a strong, constant, and uniform magnetic field perturbed by a turbulent electrostatic potential. The second one, with an infinite number of degrees of freedom, comes from physical chemistry and corresponds to Kohn-Sham time-dependent density-functional theory. For both we lay out the extension procedure and stability condition for numerical integration of the dynamics using high-order symplectic split-operator schemes. We also identify a computationally inexpensive metric that can be used for on-the-fly estimation of the accuracy of simulations. Our work paves the way for broad application of symplectic split-operator integration of classical and quantum Hamiltonian systems with finite and infinite number of degrees of freedom by comparing different modes of implementation of extended phase space integration.

16 pages, 7 figures

References in corpus (14)

• GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
• Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
• New families of symplectic splitting methods for numerical integration in dynamical astronomy
• Strong field double ionization : The phase space perspective
• Explicit methods in extended phase space for inseparable Hamiltonian problems
• Optimized Forest-Ruth- and Suzuki-like algorithms for integration of motion in many-body systems
• Capturing photoelectron motion with guiding fictitious particles
• Semiexplicit Symplectic Integrators for Non-separable Hamiltonian Systems
• Integrability and chaos: the classical uncertainty
• High order regularised symplectic integrator for collisional planetary systems
• Explicit K-symplectic methods for nonseparable non-canonical Hamiltonian systems
• QMol-grid: A MATLAB package for quantum-mechanical simulations in atomic and molecular systems
• Beatification: Flattening the Poisson Bracket for Two-Dimensional Fluid and Plasma Theories
• Hamiltonian formulation and symplectic split-operator schemes for time-dependent density-functional-theory equations of electron dynamics in molecules