On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability

Abstract

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct represenation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.

25 pages, 6 figures

References in corpus (4)

• On the structure of relativistic collisionless shocks in electron-ion plasmas
• Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• Particle-In-Cell Modeling of Plasma-Based Accelerators in Two and Three Dimensions

Cited by in corpus (9)

• Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
• Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
• Accurate simulation of direct laser acceleration in a laser wakefield accelerator
• Fundamental form of the electrostatic $\delta f$-PIC algorithm and discovery of a converged numerical instability
• 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
• Finite spatial-grid effects in energy-conserving particle-in-cell algorithms
• Particle Acceleration in Relativistic Electron-Ion Outlfows
• An \mathcal{O}(N) Maxwell solver with improved numerical dispersion properties