A Particle Module for the PLUTO Code: I - an implementation of the MHD-PIC equations

Abstract

We describe an implementation of a particle physics module available for the PLUTO code, appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a non-thermal component represented by relativistic charged particles, or cosmic rays (CR). While the fluid is approached using standard numerical schemes for magnetohydrodynamics, CR particles are treated kinetically using conventional Particle-In-Cell (PIC) techniques. The module can be used to describe either test particles motion in the fluid electromagnetic field or to solve the fully coupled MHD-PIC system of equations with particle backreaction on the fluid as originally introduced by \cite{Bai_etal.2015}. Particle backreaction on the fluid is included in the form of momentum-energy feedback and by introducing the CR-induced Hall term in Ohm's law. The hybrid MHD-PIC module can be employed to study CR kinetic effects on scales larger than the (ion) skin depth provided the Larmor gyration scale is properly resolved. When applicable, this formulation avoids to resolve microscopic scales offering a substantial computational saving with respect to PIC simulations. We present a fully-conservative formulation which is second-order accurate in time and space and extends to either Runge-Kutta (RK) or corner-transport-upwind (CTU) time-stepping schemes (for the fluid) while a standard Boris integrator is employed for the particles. For highly-energetic relativistic CRs and in order to overcome the time step restriction a novel sub-cycling strategy that retains second-order accuracy in time is presented. Numerical benchmarks and applications including Bell instability, diffusive shock acceleration and test particle acceleration in reconnecting layers are discussed.

arXiv comment27 pages, 16 figures. Accepted for publication in ApJ Supplement series

References in corpus (39)

• PLUTO: a Numerical Code for Computational Astrophysics
• Relativistic Reconnection: an Efficient Source of Non-Thermal Particles
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport
• Simulations of Ion Acceleration at Non-relativistic Shocks. I. Acceleration Efficiency
• The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics
• The Maximum Energy of Accelerated Particles in Relativistic Collisionless Shocks
• Modeling the Emission Processes in Blazars
• Simulations of particle acceleration beyond the classical synchrotron burnoff limit in magnetic reconnection: An explanation of the Crab flares
• A Reconnection Switch to Trigger Gamma-Ray Burst Jet Dissipation
• Modeling the spectral evolution of PWNe inside SNRs
• Reconnection-driven plasmoids in blazars: fast flares on a slow envelope
• Prompt GRB emission from gradual energy dissipation
• Long Term Evolution of Magnetic Turbulence in Relativistic Collisionless Shocks: Electron-Positron Plasmas
• A kinetic approach to cosmic ray induced streaming instability at supernova shocks
• High Resolution 3D Relativistic MHD Simulations of Jets
• Simulations and Theory of Ion Injection at Non-relativistic Collisionless Shocks
• High-order conservative finite difference GLM-MHD schemes for cell-centered MHD
• Simulated synchrotron emission from Pulsar Wind Nebulae
• Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks
• Formation of dynamical structures in relativistic jets: the FRI case
• Simulations of Ion Acceleration at Non-relativistic Shocks. III. Particle Diffusion
• A Second-Order Unsplit Godunov Scheme for Cell-Centered MHD: the CTU-GLM scheme
• Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. III. Rotating Relativistic Jets
• Pulsar-Wind Nebulae: Recent Progress in Observations and Theory
• Three-dimensional Simulations of AGN Jets: Magnetic Kink Instability versus Conical Shocks
• Cosmic-ray-induced filamentation instability in collisionless shocks
• Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics
• dHybrid: a massively parallel code for hybrid simulations of space plasmas
• Numerical modelling of the lobes of radio galaxies in cluster environments - III. Powerful relativistic and nonrelativistic jets
• Prompt Gamma Ray Burst emission from gradual magnetic dissipation
• On magnetic field amplification and particle acceleration near non-relativistic astrophysical shocks: Particles in MHD Cells simulations
• Multi-D magnetohydrodynamic modelling of pulsar wind nebulae: recent progress and open questions
• The emission mechanism in magnetically dominated GRB outflows
• Modelling the Kinked Jet of the Crab Nebula
• Three-dimensional modeling of the asymmetric blast wave from the 2006 outburst of RS Ophiuchi: Early X-ray emission
• Collisionless shocks in a partially ionized medium: III. Efficient cosmic ray acceleration
• Three Dimensional Structure of Relativistic Jet Formation
• Photon-Plasma: a modern high-order particle-in-cell code
• Modeling the shock-cloud interaction in SN 1006: unveiling the origin of nonthermal X-ray and gamma-ray emission

Cited by in corpus (46)

• The effects of cosmic rays on the formation of Milky Way-mass galaxies in a cosmological context
• A Particle Module for the PLUTO code: II - Hybrid Framework for Modeling Non-thermal emission from Relativistic Magnetized flows
• First Predicted Cosmic Ray Spectra, Primary-to-Secondary Ratios, and Ionization Rates from MHD Galaxy Formation Simulations
• Evolution of cosmic ray electron spectra in magnetohydrodynamical simulations
• Simulating the dynamics and synchrotron emission from relativistic jets II. Evolution of non-thermal electrons
• PIC Simulation Methods for Cosmic Radiation and Plasma Instabilities
• Amplitudes and energy fluxes of simulated decayless kink oscillations
• Wave heating in gravitationally stratified coronal loops in the presence of resistivity and viscosity
• Resistive relativistic MHD simulations of astrophysical jets
• dHybridR: a Hybrid--Particle-in-Cell Code Including Relativistic Ion Dynamics
• A Particle Module for the PLUTO Code: III -- Dust
• Influence of Ion-Neutral Damping on the Cosmic-Ray Streaming Instability: Magnetohydrodynamic Particle-in-cell Simulations
• Three-dimensional simulations of non-resonant streaming instability and particle acceleration near non-relativistic astrophysical shocks
• Dust dynamics in RAMSES -- I. Methods and turbulent acceleration
• The Magnetohydrodynamic-Particle-In-Cell Module in Athena++: Implementation and Code Tests
• A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks
• A comparison study of extrapolation models and empirical relations in forecasting solar wind
• Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse
• On the growth of the thermally modified non-resonant streaming instability
• Modeling kilonova afterglows: Effects of the thermal electron population and interaction with GRB outflows
• Alfvénic motions arising from asymmetric acoustic wave drivers in solar magnetic structures
• Particle acceleration by magnetic reconnection in relativistic jets: the transition from small to large scales
• PeV proton acceleration in Gamma-ray Binaries
• Modelling the connection between propagating disturbances and solar spicules
• Radiative pulsar magnetospheres: oblique rotators
• Magnetic reconnection-driven turbulence and turbulent reconnection acceleration
• A guiding center implementation for relativistic particle dynamics in the PLUTO code
• Acceleration of Solar Energetic Particles by the shock of Interplanetary Coronal Mass Ejection
• Effects of a Velocity Shear on Double Current Sheet Systems: Explosive Reconnection and Particle Acceleration
• PyPLUTO: a data analysis Python package for the PLUTO code
• On the influence of supra-thermal particle acceleration on the morphology of low-Mach, high-$\beta$ shocks
• Role of magnetic arcades in explaining the puzzle of the gamma-ray emission from the solar disk
• Enhancement of the non-resonant streaming instability by particle collisions
• The impact of resistive electric fields on particle acceleration in reconnection layers
• A Reduced-Speed-of-Light Formulation of the Magnetohydrodynamic-Particle-in-Cell Method
• The PLUTO Code on GPUs: Offloading Lagrangian Particle Methods
• Studying the properties of reconnection-driven turbulence
• Multi-scale simulations of particle acceleration in astrophysical systems
• On the Impact of the Numerical Method on Magnetic Reconnection and Particle Acceleration -- I. The MHD case
• Interplay between the non-resonant streaming instability and self-generated pressure anisotropies
• How many VHE gamma-ray binaries with young pulsars can be observed?
• From electrons to Janskys: Full stokes polarized radiative transfer in 3D relativistic particle-in-cell jet simulations
• Acoustic wave properties in footpoints of coronal loops in 3D MHD simulations
• Planetesimal gravitational collapse in a gaseous environment: Thermal and dynamic evolution
• Numerical modeling and physical interplay of stochastic turbulent acceleration for non-thermal emission processes
• Direct simulations of very high energy cosmic ray acceleration in 3D MHD model of a compact star cluster