Scalable spectral solver in Galilean coordinates for eliminating the numerical Cherenkov instability in particle-in-cell simulations of streaming plasmas.

Discretizing Maxwell's equations in Galilean (comoving) coordinates allows the derivation of a pseudospectral solver that eliminates the numerical Cherenkov instability for electromagnetic particle-in-cell simulations of relativistic plasmas flowing at a uniform velocity. Here we generalize this solver by incorporating spatial derivatives of arbitrary order, thereby enabling efficient parallelization by domain decomposition. This allows scaling of the algorithm to many distributed compute units.