Diagnosing symplecticity in simulations of high-dimensional Hamiltonian systems

Abstract

Integrals of the Liouville $1$-form, known as the first Poincar\'e integral invariant, provide a computable figure of merit for monitoring the conservation of symplecticity in the numerical integration of Hamiltonian systems. These integrals may be approximated with spectral convergence in the number of sample points, limited only by the regularity of the Hamiltonian. We devise a numerical integral invariant diagnostic for checking preservation of symplecticity in particle-in-cell (PIC) kinetic plasma simulation codes. As a first application of this diagnostic tool, we check the preservation of symplecticity in symplectic electrostatic particle-in-cell (PIC) methods. Surprisingly, such PIC methods fail to have symplectic time-advance maps if the charge is interpolated to the grid using linear shape functions, as is commonly done in practice. It is found that at least quadratic interpolation is needed for a structure-preserving PIC method to truly be symplectic.