Abstract:The state of art of charge-conserving electromagnetic finite element particle-in-cell has grown by leaps and bounds in the past few years. These advances have primarily been achieved for leap-frog time stepping schemes for Maxwell solvers, in large part, due to the method strictly following the proper space for representing fields, charges, and measuring currents. Unfortunately, leap-frog based solvers (and their other incarnations) are only conditionally stable. Recent advances have made Electromagnetic Finite Element Particle-in-Cell (EM-FEMPIC) methods built around unconditionally stable time stepping schemes were shown to conserve charge. Together with the use of a quasi-Helmholtz decomposition, these methods were both unconditionally stable and satisfied Gauss' Laws to machine precision. However, this architecture was developed for systems with explicit particle integrators where fields and velocities were off by a time step. While completely self-consistent methods exist in the literature, they follow the classic rubric: collect a system of first order differential equations (Maxwell and Newton equations) and use an integrator to solve the combined system. These methods suffer from the same side-effect as earlier--they are conditionally stable. Here we propose a different approach; we pair an unconditionally stable Maxwell solver to an exponential predictor-corrector method for Newton's equations. As we will show via numerical experiments, the proposed method conserves energy within a PIC scheme, has an unconditionally stable EM solve, solves Newton's equations to much higher accuracy than a traditional Boris solver and conserves charge to machine precision. We further demonstrate benefits compared to other polynomial methods to solve Newton's equations, like the well known Boris push.

Relativistic extension of a charge-conservative finite element solver for time-dependent Maxwell-Vlasov equations

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

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

Mitigation of Numerical Cerenkov Radiation and Instability Using a Hybrid Finite Difference-Fft Maxwell Solver and a Local Charge Conserving Current Deposit

A Charge Conserving Exponential Predictor Corrector FEMPIC Formulation for Relativistic Particle Simulations

A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

Relativistic reference frame PIC simulations for electron beam dynamics with meter-scale propagation inside plasma and under external fields

A fully implicit, scalable, conservative nonlinear relativistic Fokker-Planck 0D-2P solver for runaway electron

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

Relativistic Vlasov-Maxwell modelling using finite volumes and adaptive mesh refinement

General relativistic force-free electrodynamics with a discontinuous Galerkin-finite difference hybrid method

Loading Relativistic Maxwell Distributions in Particle Simulations

Explicit energy-conserving modification of relativistic PIC method

Particle-in-Cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Three-Dimensional Pic-Mc Modeling for Relativistic Electron Beam Transport Through Dense Plasma

An Unstructured Body-of-Revolution Electromagnetic Particle-in-Cell Algorithm with Radial Perfectly Matched Layers and Dual Polarizations

Higher-Order Accurate and Conservative Hybrid Numerical Scheme for Relativistic Time-Fractional Vlasov-Maxwell System