Abstract:In a previous paper, we developed a new particle-in-cell method for the Vlasov-Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell's equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in an non-staggered PIC formulation.

Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

Geometric particle-in-cell methods for Vlasov--Poisson equations with Maxwell--Boltzmann electrons

Universal Vector–Scalar Potential Framework for Inhomogeneous Electromagnetic System and Its Application in Semiclassical Quantum Electromagnetics

A Unified Hamiltonian Solution to Maxwell–Schrödinger Equations for Modeling Electromagnetic Field–particle Interaction

Geometric Particle-in-Cell Simulations of the Vlasov--Maxwell System in Curvilinear Coordinates

GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods

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

Variational Framework for Structure-Preserving Electromagnetic Particle-in-Cell Methods

An efficient implicit boundary integral solver for the Vlasov-Maxwell system

Discontinuous Galerkin Methods for the Vlasov-Maxwell Equations

A Particle-in-cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

A Tutorial on Variational Integrators

Variational Integration for Ideal Magnetohydrodynamics with Built-in Advection Equations

Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

A numerical integration scheme for vectorised phase-space of one-dimensional collision-free, electrostatic systems

Variational Symplectic Integrator for Long-Time Simulations of the Guiding-Center Motion of Charged Particles in General Magnetic Fields.

A Particle-in-cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

Variational formulation of macro-particle plasma simulation algorithms

The Hamiltonian Structure and Euler-Poincaré Formulation of the Vlasov-Maxwell and Gyrokinetic Systems

Conservative semi-Lagrangian kinetic scheme coupled with implicit finite element field solver for multidimensional Vlasov Maxwell system