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.

Large-scale Parallel Particle-in-cell Code CHIPIC

Three-dimensional Electromagnetic Particle-in-cell Simulation by Parallel Computing

Parallel 3D electromagnetic PIC simulation based on OpenMP

Implementation of parallel plasma particle-in-cell codes on PC cluster

Three-dimensional Parallelized Fully-Electromagnetic Particle Simulation Code UNIPIC-3D

PIFE-PIC: Parallel Immersed-Finite-Element Particle-In-Cell For 3-D Kinetic Simulations of Plasma-Material Interactions

Particle-in-Cell Code for GPU Systems

UNIPIC Code for Simulations of High Power Microwave Devices

QuickPIC: a highly efficient fully parallelized PIC code for plasma-based acceleration

Characterizing the Performance of the Implicit Massively Parallel Particle-in-Cell iPIC3D Code

Parallelization And Optimization Of Electrostatic Particle-In-Cell/Monte-Carlo Coupled Codes As Applied To Rf Discharges

Electromagnetic Field Algorithms of CHIPIC Code

Unipic: A Conformal Particle-In-Cell Code for High Power Micrwowave and Pulse Discharge

2.5-Dimensional Electromagnetic Particle-in-cell Code-Unipic for High Power Microwave Simulations

CPIC: a Curvilinear Particle-In-Cell code for plasma-material interaction studies

The magnetohydrodynamic-particle-in-cell module in athena plus plus : implementation and code tests

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

A NVIDA kepler based acceleration of PIC method

Performance Analysis and Application of 2-D Plasma Parallel Particle-in-cell Simulations

PHARE : Parallel hybrid particle-in-cell code with patch-based adaptive mesh refinement

Parallel Algorithm Design and Implementation of Electromagnetic Simulation in EM PIC Software