Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems

• Published in: Physics of plasmas Vol. 22; no. 11
• Main Authors: Xiao, Jianyuan, Qin, Hong, Liu, Jian, He, Yang, Zhang, Ruili, Sun, Yajuan
• Format: Journal Article
• Language: English
• Published: Melville American Institute of Physics 01.11.2015
• Subjects: Algorithms; Bernstein waves; Computer simulation; Differential calculus; Electromagnetic fields; Landau damping; Maxwell's equations; Particle in cell technique; Plasma physics; Supercomputers
• ISSN: 1070-664X; 1089-7674
• DOI: 10.1063/1.4935904

Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithms conserve a discrete non-canonical symplectic structure derived from the Lagrangian of the particle-field system, which is naturally discrete in particles. The electromagnetic field is spatially discretized using the method of discrete exterior calculus with high-order interpolating differential forms for a cubic grid. The resulting time-domain Lagrangian assumes a non-canonical symplectic structure. It is also gauge invariant and conserves charge. The system is then solved using a structure-preserving splitting method discovered by He et al. [preprint arXiv:1505.06076 (2015)], which produces five exactly soluble sub-systems, and high-order structure-preserving algorithms follow by combinations. The explicit, high-order, and conservative nature of the algorithms is especially suitable for long-term simulations of particle-field systems with extremely large number of degrees of freedom on massively parallel supercomputers. The algorithms have been tested and verified by the two physics problems, i.e., the nonlinear Landau damping and the electron Bernstein wave.