Modeling of laser wakefield acceleration in the Lorentz boosted frame using OSIRIS and UPIC framework

Abstract

We examine numerical issues related to modeling laser wakefield acceleration (LWFA) in a Lorentz boosted frame. We present some theoretical results on a numerical instability that results from a relativistic plasma drift on a grid. We compare these predictions with results obtained from a finite difference time domain (FDTD) code OSIRIS and a spectral (FFT) time domain code UPIC. Consistent with the theoretical predictions, it is found that the spectral solver has advantages for mitigating the numerical instability when compared against some prevailing (FDTD) solvers. We also compared using finite difference or spectral field solvers in simplified 2D LWFA simulations using a plane wave for the laser. We find that the use of the spectral solver has advantages for mitigating the numerical instability induced by the relativistically drifting plasma.