The family of generalized Pseudo-Spectral Time Domain (including the Pseudo-Spectral Analytical Time Domain) Particle-in-Cell algorithms offers substantial versatility for simulating particle beams and plasmas, and well written codes using these algo
rithms run reasonably fast. When simulating relativistic beams and streaming plasmas in multiple dimensions, they are, however, subject to the numerical Cherenkov instability. Previous studies have shown that instability growth rates can be reduced substantially by modifying slightly the transverse fields as seen by the streaming particles . Here, we offer an approach which completely eliminates the fundamental mode of the numerical Cherenkov instability while minimizing the transverse field corrections. The procedure, numerically computed residual growth rates (from weaker, higher order instability aliases), and comparisons with WARP simulations are presented. In some instances, there are no numerical instabilities whatsoever, at least in the linear regime.
A procedure for largely suppressing the numerical Cherenkov instability in finite difference time-domain (FDTD) particle-in-cell (PIC) simulations of cold, relativistic beams is derived, and residual growth rates computed and compared with WARP code
simulation results. Sample laser-plasma acceleration simulation output is provided to further validate the new procedure.
The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwells equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle be
am simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper presents several approaches that, when combined with digital filtering, almost completely eliminate the numerical Cherenkov instability. The paper also investigates the numerical stability of the PSATD algorithm at low beam energies.
The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwells equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle be
am simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper derives and solves the numerical dispersion relation for the PSATD algorithm and compares the results with corresponding behavior of the more conventional pseudo-spectral time-domain (PSTD) and finite difference time-domain (FDTD) algorithms. In general, PSATD offers superior stability properties over a reasonable range of time steps. More importantly, one version of the PSATD algorithm, when combined with digital filtering, is almost completely free of the numerical Cherenkov instability for time steps (scaled to the speed of light) comparable to or smaller than the axial cell size.
A method is proposed to generate low emittance electron bunches from two color laser pulses in a laser-plasma accelerator. A two-region gas structure is used, containing a short region of a high-Z gas (e.g., krypton) for ionization injection, followe
d by a longer region of a low-Z gas for post-acceleration. A long-laser-wavelength (e.g., 5 micron) pump pulse excites plasma wake without triggering the inner-shell electron ionization of the high-Z gas due to low electric fields. A short-laser-wavelength (e.g., 0.4 micron) injection pulse, located at a trapping phase of the wake, ionizes the inner-shell electrons of the high-Z gas, resulting in ionization-induced trapping. Compared with a single-pulse ionization injection, this scheme offers an order of magnitude smaller residual transverse momentum of the electron bunch, which is a result of the smaller vector potential amplitude of the injection pulse.
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.
