Explicit energy-conserving modification of relativistic PIC method

Abstract

The use of explicit particle-in-cell (PIC) method for relativistic plasma simulations is restricted by numerical heating and instabilities that may significantly constrain the choice of time and space steps. To partially eliminate these limitations we consider a possibility to enforce exact energy conservation by altering the standard time step splitting. Instead of updating particles in a given field and then the field using the current they produce, we consider subsystems that describe the coupling of a single particle and the field at the nearby nodes and solve them with enforced energy conservation sequentially for all particles, which is completed by the field update with zero current. Such an approach is compatible with various advances, ranging from accounting for additional physical effects to the use of numerical-dispersion-free field solvers, high-order weighting shapes and particle push subcycling. To facilitate further considerations and use, we provide a basic implementation in a 3D, relativistic, spectral code $\pi$-PIC, which we make publicly available. The method and its implementations are verified using simulations of cold plasma oscillations, Landau damping and relativistic two-stream instability. The capabilities of the method to deal with large time and space steps are demonstrated in the problem of plasma heating by intense incident radiation.

References in corpus (8)

• QED cascades induced by circularly polarized laser fields
• Time-dependent pair cascades in magnetospheres of neutron stars I. Dynamics of the polar cap cascade with no particle supply from the neutron star surface
• VPIC 2.0: Next Generation Particle-in-Cell Simulations
• Temporal resolution criterion for correctly simulating relativistic electron motion in a high-intensity laser field
• Strategies for particle resampling in PIC simulations
• Efficient modeling of laser plasma interactions in high energy density scenarios
• Particle integrator for particle-in-cell simulations of ultra-high intensity laser-plasma interactions
• Optimized computation of tight focusing of short pulses using mapping to periodic space

Cited by in corpus (1)

• Grid instability growth rates for explicit, electrostatic momentum- and energy-conserving particle-in-cell algorithms