Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm

Abstract

Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole-Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.

Cited by in corpus (37)

• A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm
• Implementation of a hybrid particle code with a PIC description in r-z and a gridless description in $\varphi$ into OSIRIS
• High quality electron bunch generation using a longitudinal density-tailored plasma-based accelerator in the three-dimensional blowout regime
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• Physics of Weibel-mediated relativistic collisionless shocks
• Suppressing the Numerical Cherenkov Instability in FDTD PIC Codes
• Persistence of precursor waves in two-dimensional relativistic shocks
• Numerical stability analysis of the Pseudo-Spectral Analytical Time-Domain PIC algorithm
• Temporal resolution criterion for correctly simulating relativistic electron motion in a high-intensity laser field
• Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
• PIC Simulation Methods for Cosmic Radiation and Plasma Instabilities
• Elimination of Numerical Cherenkov Instability in flowing-plasma Particle-In-Cell simulations by using Galilean coordinates
• Stable discrete representation of relativistically drifting plasmas
• Physics and phenomenology of weakly magnetized, relativistic astrophysical shock waves
• Improved Numerical Cherenkov Instability Suppression in the Generalized PSTD PIC Algorithm
• Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
• Attosecond-Angstrom free-electron-laser towards the cold beam limit
• Modeling of Laser wakefield acceleration in Lorentz boosted frame using EM-PIC code with spectral solver
• Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
• Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
• 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
• Stability Property of Numerical Cherenkov Radiation and its Application to Relativistic Shock Simulations
• Relativistic Shear Flow Between Electron-Ion and Electron-Positron Plasmas and Astrophysical Applications
• Microphysics of Relativistic Collisionless Electron-ion-positron Shocks
• Numerical Stability Improvements for the Pseudo-Spectral EM PIC Algorithm
• On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability
• Lorentz boosted frame simulation of Laser wakefield acceleration in quasi-3D geometry
• On the generation of ultra-bright and low energy spread electron beams in laser wakefield acceleration in a uniform plasma
• 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
• PIC simulation study of the interaction between a relativistically moving leptonic micro-cloud and ambient electrons
• Stopping Power Enhancement From Discrete Particle-Wake Correlations in High Energy Density Plasmas
• A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
• EZ: An Efficient, Charge Conserving Current Deposition Algorithm for Electromagnetic Particle-In-Cell Simulations
• Symplectic Modeling of Beam Loading in Electromagnetic Cavities
• Modeling of Advanced Accelerator Concepts