@article{70, author = {Matthew Kunz and James Stone and Xue-Ning Bai}, title = {Pegasus: A new hybrid-kinetic particle-in-cell code for astrophysical plasma dynamics}, abstract = {

We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order{\textendash}accurate, three-stage predictor{\textendash}predictor{\textendash}corrector integration algorithm. The constrained transport method is used to enforce the divergence-free constraint on the magnetic field. A δf scheme is included to facilitate a reduced-noise study of systems in which only small departures from an initial distribution function are anticipated. The effects of rotation and shear are implemented through the shearing-sheet formalism with orbital advection. These algorithms are embedded within an architecture similar to that used in the popular astrophysical magnetohydrodynamics code Athena, one that is modular, well-documented, easy to use, and efficiently parallelized for use on thousands of processors. We present a series of tests in one, two, and three spatial dimensions that demonstrate the fidelity and versatility of the code.

}, year = {2014}, journal = {J. Computational Physics}, volume = {259}, chapter = {154}, url = {https://doi.org/10.1016/j.jcp.2013.11.035}, language = {eng}, }