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.

arXiv comment20 pages, 19 figures, submitted to MNRAS

References in corpus (58)

• The Interstellar Environment of our Galaxy
• Athena: A New Code for Astrophysical MHD
• Theoretical Challenges in Galaxy Formation
• Relativistic Reconnection: an Efficient Source of Non-Thermal Particles
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport
• Simulations of Ion Acceleration at Non-relativistic Shocks. I. Acceleration Efficiency
• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• Particle Acceleration in Relativistic Magnetized Collisionless Electron-Ion Shocks
• 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
• Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos
• Simulations and Theory of Ion Injection at Non-relativistic Collisionless Shocks
• Magnetic reconnection in the era of exascale computing and multiscale experiments
• Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks
• Kinetic simulations of mildly relativistic shocks I: particle acceleration in high Mach number shocks
• Reconnection-driven particle acceleration in relativistic shear flows
• Efficient Nonthermal Particle Acceleration by the Kink Instability in Relativistic Jets
• Cosmic Ray Scattering in Compressible Turbulence
• A Particle Module for the PLUTO Code: I - an implementation of the MHD-PIC equations
• Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics
• Electron Heating by the Ion Cyclotron Instability in Collisionless Accretion Flows. I. Compression-Driven Instabilities and the Electron Heating Mechanism
• dHybrid: a massively parallel code for hybrid simulations of space plasmas
• Kinetic Simulations of Cosmic-Ray-Modified Shocks II: Particle Spectra
• Kinetic Simulations of Cosmic-Ray-Modified Shocks I: Hydrodynamics
• A filamentation instability for streaming cosmic-rays
• On magnetic field amplification and particle acceleration near non-relativistic astrophysical shocks: Particles in MHD Cells simulations
• A computational model for exploring particle acceleration during reconnection in macro-scale systems
• First Predicted Cosmic Ray Spectra, Primary-to-Secondary Ratios, and Ionization Rates from MHD Galaxy Formation Simulations
• Magnetohydrodynamic-Particle-in-Cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution
• VPIC 2.0: Next Generation Particle-in-Cell Simulations
• The Acceleration of Thermal Protons at Parallel Collisionless Shocks: Three-dimensional Hybrid Simulations
• An efficient and portable SIMD algorithm for charge/current deposition in Particle-In-Cell codes
• On the Growth and Saturation of the Gyroresonant Streaming Instabilities
• Particle acceleration in kink-unstable jets
• Cosmic Ray transport through gyroresonance instability in compressible turbulence
• Adaptive critical balance and firehose instability in an expanding, turbulent, collisionless plasma
• Bottom-up Acceleration of Ultra-High-Energy Cosmic Rays in the Jets of Active Galactic Nuclei
• Electron Heating by the Ion Cyclotron Instability in Collisionless Accretion Flows. II. Electron Heating Efficiency as a Function of Flow Conditions
• Propagation of Alfvén waves in the expanding solar wind with the fast-slow stream interaction
• Influence of Ion-Neutral Damping on the Cosmic-Ray Streaming Instability: Magnetohydrodynamic Particle-in-cell Simulations
• Large-Scale Parallel Electric Fields and Return Currents in a Global Simulation Model
• Kinetic-MHD simulations of gyroresonance instability driven by CR pressure anisotropy
• Adaptive SIMD optimizations in particle-in-cell codes with fine-grain particle sorting
• Towards First-principle Characterization of Cosmic-ray Transport Coefficients from Multi-scale Kinetic Simulations
• The Role of Pressure Anisotropy in Cosmic Ray Hydrodynamics
• Cosmic Ray Transport in the Ionized and Neutral ISM: MHD-PIC Simulations and Effective Fluid Treatments
• Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets
• A guiding center implementation for relativistic particle dynamics in the PLUTO code
• 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
• 3D PIC Simulations for Relativistic Jets with a Toroidal Magnetic Field
• Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities
• On the Conservation of Cross Helicity and Wave Action in Solar-Wind Models with Non-WKB Alfven Wave Reflection
• Mushroom-instability-driven Magnetic Reconnections in Collisionless Relativistic Jets
• A Reduced-Speed-of-Light Formulation of the Magnetohydrodynamic-Particle-in-Cell Method

Cited by in corpus (2)

• Cosmic Ray Feedback on Bi-stable ISM Turbulence
• Linear Analysis and Simulations of the Cosmic-Ray Streaming Instability: the Importance of Oblique Waves