Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities

Abstract

We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on $\sim \mathrm{AU}$ scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains drift super-Alfv{é}nically, with Lorentz forces strongly dominating over drag forces. Dust grains are unstable to the RDIs and form concentrated columns and sheets, whose scale grows until saturating at the simulation box size. Initially perfectly-streaming CRs are strongly scattered by RDI-excited Alfv{é}n waves, with the growth rate of the CR perpendicular velocity components equaling the growth rate of magnetic field perturbations. These rates are well-predicted by analytic linear theory. CRs finally become isotropized and drift at least at $\sim v_{\mathrm{A}}$ by unidirectional Alfvén waves excited by the RDIs, with a uniform distribution of the pitch angle cosine $\mu$ and a flat profile of the CR pitch angle diffusion coefficient $D_{\mu \mu }$ around $\mu =0$, without the "$90$ degree pitch angle problem." With CR feedback on the gas included, $D_{\mu \mu }$ decreases by a factor of a few, indicating a lower CR scattering rate, because the backreaction on the RDI from the CR pressure adds extra wave damping, leading to lower quasi-steady-state scattering rates. Our study demonstrates that the dust-induced CR confinement can be very important under certain conditions, e.g., the dusty circumgalactic medium around quasars or superluminous galaxies.

14 pages, 10 figures, published on MNRAS

References in corpus (9)

• The NumPy array: a structure for efficient numerical computation
• Feedback Heating by Cosmic Rays in Clusters of Galaxies
• Cosmic Ray Propagation: Nonlinear Diffusion Parallel and Perpendicular to Mean Magnetic Field
• Cosmic-ray propagation with DRAGON2: I. numerical solver and astrophysical ingredients
• Transverse Sizes of CIV Absorption Systems Measured from Multiple QSO Sightlines
• Magnetohydrodynamic-Particle-in-Cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution
• Kinematics of solid particles in a turbulent protoplanetary disc
• The impact of astrophysical dust grains on the confinement of cosmic rays
• Local simulations of MRI turbulence with meshless methods

Cited by in corpus (1)

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