The Magnetohydrodynamic-Particle-In-Cell Module in Athena++: Implementation and Code Tests

Abstract

We present a new magnetohydrodynamic-particle-in-cell (MHD-PIC) code integrated into the Athena++ framework. It treats energetic particles as in conventional PIC codes while the rest of thermal plasmas are treated as background fluid described by MHD, thus primarily targeting at multi-scale astrophysical problems involving the kinetic physics of the cosmic-rays (CRs). The code is optimized toward efficient vectorization in interpolation and particle deposits, with excellent parallel scaling. The code is also compatible with static/adaptive mesh refinement, with dynamic load balancing to further enhance multi-scale simulations. In addition, we have implemented a compressing/expanding box framework which allows adiabatic driving of CR pressure anisotropy, as well as the $\delta f$ method that can dramatically reduce Poisson noise in problems where distribution function $f$ is only expected to slightly deviate from the background. The code performance is demonstrated over a series of benchmark test problems including particle acceleration in non-relativistic parallel shocks. In particular, we reproduce the linear growth of the CR gyro-resonant (streaming and pressure anisotropy) instabilities, under both the periodic and expanding/compressing box setting. We anticipate the code to open up the avenue for a wide range of astrophysical and plasma physics applications.

20 pages, 19 figures, submitted to MNRAS

References in corpus (22)

• Athena: A New Code for Astrophysical MHD
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• Simultaneous Acceleration of Protons and Electrons at Nonrelativistic Quasiparallel Collisionless Shocks
• Simulations of relativistic collisionless shocks: shock structure and particle acceleration
• In-situ switchback formation in the expanding solar wind
• Simulations and Theory of Ion Injection at Non-relativistic Collisionless Shocks
• Magnetic reconnection in the era of exascale computing and multiscale experiments
• Cosmic Ray Scattering in Compressible Turbulence
• Electron Heating by the Ion Cyclotron Instability in Collisionless Accretion Flows. I. Compression-Driven Instabilities and the Electron Heating Mechanism
• On magnetic field amplification and particle acceleration near non-relativistic astrophysical shocks: Particles in MHD Cells simulations
• Magnetohydrodynamic-Particle-in-Cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution
• Cosmic Ray transport through gyroresonance instability in compressible turbulence
• Propagation of Alfvén waves in the expanding solar wind with the fast-slow stream interaction
• Electron Heating by the Ion Cyclotron Instability in Collisionless Accretion Flows. II. Electron Heating Efficiency as a Function of Flow Conditions
• Towards First-principle Characterization of Cosmic-ray Transport Coefficients from Multi-scale Kinetic Simulations
• The Role of Pressure Anisotropy in Cosmic Ray Hydrodynamics
• Diffusive shock acceleration at oblique high Mach number shocks
• Importance of accurate consideration of the electron inertia in hybrid-kinetic simulations of collisionless plasma turbulence: 1. The 2D limit
• Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities
• Mushroom-instability-driven Magnetic Reconnections in Collisionless Relativistic Jets