  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.
