Influence of Ion–Neutral Damping on the Cosmic-Ray Streaming Instability: Magnetohydrodynamic Particle-in-cell Simulations

• Published in: The Astrophysical journal Vol. 914; no. 1; pp. 3 - 19
• Main Authors: Plotnikov, Illya, Ostriker, Eve C., Bai, Xue-Ning
• Format: Journal Article
• Language: English
• Published: Philadelphia The American Astronomical Society 01.06.2021; IOP Publishing
• Subjects: Alfven waves; Amplitudes; Astrophysics; Collision rates; Computational fluid dynamics; Computational methods; Cosmic ray showers; Cosmic rays; Damping; Distribution functions; Fluid flow; Growth rate; Instability; Interstellar chemistry; Interstellar gas; Interstellar matter; Interstellar medium; Isotropy; Linear phase; Magnetic fields; Magnetohydrodynamics; Molecular gases; Particle in cell technique; Stability; Stability analysis; Wave spectra; Wavelengths
• ISSN: 0004-637X; 1538-4357
• DOI: 10.3847/1538-4357/abf7b3

Abstract We explore the physics of the gyroresonant cosmic-ray streaming instability (CRSI) including the effects of ion–neutral (IN) damping. This is the main damping mechanism in (partially ionized) atomic and molecular gas, which are the primary components of the interstellar medium (ISM) by mass. Limitation of CRSI by IN damping is important in setting the amplitude of Alfvén waves that scatter cosmic rays (CRs) and control galactic-scale transport. Our study employs the magnetohydrodynamic (MHD)–particle-in-cell hybrid fluid-kinetic numerical technique to follow linear growth as well as post-linear and saturation phases. During the linear phase of the instability—where simulations and analytical theory are in good agreement—IN damping prevents wave growth at small and large wavelengths, with the unstable bandwidth lower for higher IN collision rates ν in . Purely MHD effects during the post-linear phase extend the wave spectrum toward larger k . In the saturated state, the CR distribution evolves toward greater isotropy (lower streaming velocity) by scattering off of Alfvén waves excited by the instability. In the absence of low- k waves, CRs with sufficiently high momentum are not isotropized. The maximum wave amplitude and rate of isotropization of the distribution function decrease at higher ν in . When the IN damping rate approaches the maximum growth rate of CRSI, wave growth and isotropization are suppressed. Implications of our results for CR transport in partially ionized ISM phases are discussed.