English

RelSIM: A Relativistic Semi-implicit Method for Particle-in-Cell Simulations

High Energy Astrophysical Phenomena 2023-08-29 v2 Computational Physics Plasma Physics

Abstract

We present a novel Relativistic Semi-Implicit Method (RelSIM) for particle-in-cell (PIC) simulations of astrophysical plasmas, implemented in a code framework ready for production runs. While explicit PIC methods have gained widespread recognition in the astrophysical community as a reliable tool to simulate plasma phenomena, implicit methods have been seldom explored. This is partly due to the lack of a reliable relativistic implicit PIC formulation that is applicable to state-of-the-art simulations. We propose the RelSIM to fill this gap: our new method is relatively simple, being free of nonlinear iterations and only requiring a global linear solve of the field equations. With a set of one- and two-dimensional tests, we demonstrate that the RelSIM produces more accurate results with much smaller numerical errors in the total energy than standard explicit PIC, particularly when characteristic plasma scales (skin depth and plasma frequency) are heavily underresolved on the numerical grid. By construction, the RelSIM also performs much better than the Relativistic Implicit-Moment Method (RelIMM), originally proposed for semi-implicit PIC simulations in the relativistic regime. Our results are promising to conduct large-scale (in terms of duration and domain size) PIC simulations of astrophysical plasmas, potentially reaching physical regimes inaccessible by standard explicit PIC codes.

Keywords

Cite

@article{arxiv.2306.04685,
  title  = {RelSIM: A Relativistic Semi-implicit Method for Particle-in-Cell Simulations},
  author = {Fabio Bacchini},
  journal= {arXiv preprint arXiv:2306.04685},
  year   = {2023}
}

Comments

Accepted in ApJS

R2 v1 2026-06-28T10:59:14.860Z