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.

arXiv comment14 pages, 10 figures, published on MNRAS

References in corpus (51)

• The NumPy array: a structure for efficient numerical computation
• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Particle Clumping and Planetesimal Formation Depend Strongly on Metallicity
• Scattering of Cosmic Rays by Magnetohydrodynamic Interstellar Turbulence
• 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
• Measuring the galaxy-mass and galaxy-dust correlations through magnification and reddening
• But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation
• STARFORGE: Toward a comprehensive numerical model of star cluster formation and feedback
• Wave damping by MHD turbulence and its effect upon cosmic ray propagation in the ISM
• Global simulations of galactic winds including cosmic ray streaming
• Accurate, Meshless Methods for Magneto-Hydrodynamics
• The Effect of the Radial Pressure Gradient in Protoplanetary Disks on Planetesimal Formation
• Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion
• Cosmic Ray Transport in the Galaxy: a Review
• Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos
• Cosmic ray feedback in the FIRE simulations: constraining cosmic ray propagation with GeV gamma ray emission
• Diffusive propagation of cosmic rays from supernova remnants in the Galaxy. I: spectrum and chemical composition
• Binary Quasars at High Redshift I: 24 New Quasar Pairs at z ~ 3-4
• Models for Galactic cosmic-ray propagation
• Impact of Cosmic Ray Transport on Galactic Winds
• The Role of Cosmic Ray Transport in Shaping the Simulated Circumgalactic Medium
• Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks
• Simulating galactic dust grain evolution on a moving mesh
• Anisotropic Diffusion in Mesh-Free Numerical Magnetohydrodynamics
• A Constrained-Gradient Method to Control Divergence Errors in Numerical MHD
• The effects of cosmic rays on the formation of Milky Way-mass galaxies in a cosmological context
• Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies
• Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation
• Role of Cosmic Rays in the Circumgalactic Medium
• Cosmic Ray Parallel and Perpendicular Transport in Turbulent Magnetic Fields
• The Fundamentally Different Dynamics of Dust and Gas in Molecular Clouds
• Dust in the Circumgalactic Medium of Low-Redshift Galaxies
• Transverse Sizes of CIV Absorption Systems Measured from Multiple QSO Sightlines
• Resonant drag instability of grains streaming in fluids
• Numerical study of Cosmic Ray Diffusion in MHD turbulence
• Cosmic Rays or Turbulence can Suppress Cooling Flows (Where Thermal Heating or Momentum Injection Fail)
• Magnetohydrodynamic-Particle-in-Cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution
• Ubiquitous Instabilities of Dust Moving in Magnetized Gas
• On the Growth and Saturation of the Gyroresonant Streaming Instabilities
• The Dynamics of Charged Dust in Magnetized Molecular Clouds
• Kinematics of solid particles in a turbulent protoplanetary disc
• Virial shocks are suppressed in cosmic ray-dominated galaxy halos
• The impact of astrophysical dust grains on the confinement of cosmic rays
• Local simulations of MRI turbulence with meshless methods
• Nonlinear Evolution of the Resonant Drag Instability in Magnetized Gas
• Nonlinear Evolution of Instabilities Between Dust and Sound Waves
• Simulating the Diverse Instabilities of Dust in Magnetized Gas
• Three-dimensional simulations of non-resonant streaming instability and particle acceleration near non-relativistic astrophysical shocks

Cited by in corpus (4)

• First Predicted Cosmic Ray Spectra, Primary-to-Secondary Ratios, and Ionization Rates from MHD Galaxy Formation Simulations
• The Magnetohydrodynamic-Particle-In-Cell Module in Athena++: Implementation and Code Tests
• A Reduced-Speed-of-Light Formulation of the Magnetohydrodynamic-Particle-in-Cell Method
• Velocity decorrelation functions of high-energy cosmic rays propagating in magnetic fields