EZ: An Efficient, Charge Conserving Current Deposition Algorithm for Electromagnetic Particle-In-Cell Simulations

Abstract

We present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as ``EZ'' as shorthand for ``Esirkepov meets ZigZag''. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

References in corpus (10)

• SMILEI: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation
• A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm
• One-to-one direct modeling of experiments and astrophysical scenarios: pushing the envelope on kinetic plasma simulations
• Demonstration of a compact plasma accelerator powered by laser-accelerated electron beams
• Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
• DC magnetic field generation in unmagnetized shear flows
• Probing Ultrafast Magnetic-Field Generation by Current Filamentation Instability in Femtosecond Relativistic Laser-Matter Interactions
• Photon-Plasma: a modern high-order particle-in-cell code
• On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability
• Radiation from relativistic jets in turbulent magnetic fields