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

Abstract

In this paper we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift ($\hat{1}$ direction). We show that this eliminates the main NCI modes with moderate $|k_1|$, while keeps additional main NCI modes well outside the range of physical interest with higher $|k_1|$. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along $\hat{1}$ which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss' Law is satisfied. We present simulation examples of one relativistically drifting plasmas, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.

References in corpus (14)

• Generating multi-GeV electron bunches using single stage laser wakefield acceleration in a 3D nonlinear regime
• Particle acceleration in relativistic collisionless shocks: Fermi process at last?
• The Maximum Energy of Accelerated Particles in Relativistic Collisionless Shocks
• Ion dynamics and acceleration in relativistic shocks
• Weibel-instability-mediated collisionless shocks in laboratory with ultraintense lasers
• Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
• Implementation of a hybrid particle code with a PIC description in r-z and a gridless description in $\varphi$ into OSIRIS
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• Suppressing the Numerical Cherenkov Instability in FDTD PIC Codes
• Numerical stability analysis of the Pseudo-Spectral Analytical Time-Domain PIC algorithm
• Improved Numerical Cherenkov Instability Suppression in the Generalized PSTD PIC Algorithm
• Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
• Modeling of Laser wakefield acceleration in Lorentz boosted frame using EM-PIC code with spectral solver
• Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit

Cited by in corpus (24)

• High quality electron bunch generation using a longitudinal density-tailored plasma-based accelerator in the three-dimensional blowout regime
• PIC Simulation Methods for Cosmic Radiation and Plasma Instabilities
• Accurate modeling of plasma acceleration with arbitrary order pseudo-spectral particle-in-cell methods
• Stable Positron Acceleration in Thin, Warm, Hollow Plasma Channels
• Physics and phenomenology of weakly magnetized, relativistic astrophysical shock waves
• X-Dispersionless Maxwell solver for plasma-based particle acceleration
• A multi-sheath model for highly nonlinear plasma wakefields
• Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
• On numerical errors to the fields surrounding a relativistically moving particle in PIC codes
• Ultra-Bright Electron Bunch Injection in a Plasma Wakefield Driven by a Superluminal Flying Focus Electron Beam
• On the generation of ultra-bright and low energy spread electron beams in laser wakefield acceleration in a uniform plasma
• Generating high quality ultra-relativistic electron beams using an evolving electron beam driver
• Direct laser acceleration: A model for the electron injection from the walls of a cylindrical guiding structure
• Accurate simulation of direct laser acceleration in a laser wakefield accelerator
• Overcoming timestep limitations in boosted-frame Particle-In-Cell simulations of plasma-based acceleration
• Positron acceleration in plasma waves driven by non-neutral fireball beams
• Particle-in-Cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms
• Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition
• Multi-scale simulations of particle acceleration in astrophysical systems
• A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
• Stopping Power Enhancement From Discrete Particle-Wake Correlations in High Energy Density Plasmas
• An \mathcal{O}(N) Maxwell solver with improved numerical dispersion properties
• Optimized matching conditions for self-guided laser wakefield accelerators
• Modeling of Advanced Accelerator Concepts