Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

TL;DR

This paper introduces a novel integration of a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code that uses azimuthal mode decomposition, significantly enhancing simulation efficiency for plasma-based acceleration modeling.

Contribution

It presents the implementation of a PGC algorithm compatible with azimuthal mode expansion in QSA-PIC codes, enabling larger time steps and improved computational performance.

Findings

Enables larger time steps in plasma simulations

Reduces computational complexity of 3D PIC modeling

Achieves significant speedup in plasma acceleration simulations

Abstract

High-fidelity modeling of plasma-based acceleration (PBA) requires the use of 3D fully nonlinear and kinetic descriptions based on the particle-in-cell (PIC) method. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the beam-plasma interaction. In a QSA-PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA-PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA-PIC algorithm that combines another speedup technique known as azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can significantly reduce the algorithmic complexity. Modeling the laser-plasma interaction in a full 3D PIC algorithm is very computationally expensive due to the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible with the QSA PIC algorithms based on the azimuthal mode expansion. This algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described.