Testing physical models for cosmic ray transport coefficients on galactic scales: self-confinement and extrinsic turbulence at ∼GeV energies
ABSTRACT The microphysics of ∼ GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant diffusivity). We explore different physically motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological Feedback In Realistic Environment (FIRE) simulations of dwarf and ∼L* galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfvén-wave damping, etc. We self-consistently predict observables including γ-rays (Lγ), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multiphase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or Lγ), as transport is rate-limited by the ionized WIM and ‘inner CGM’ gaseous halo (104–106 K gas within $\lesssim 10\!-\!30\,$ kpc), and Lγ can be dominated by trapping in small ‘patches’. Most physical ET models contribute negligible scattering of ∼1–10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely more »
Authors:
; ; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10278897
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
501
Issue:
3
Page Range or eLocation-ID:
4184 to 4213
ISSN:
0035-8711
National Science Foundation
##### More Like this
1. ABSTRACT

Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ∼MeV-TeV, these fluctuations are usually attributed either to ‘extrinsic turbulence’ (ET) – a cascade from larger scales – or ‘self-confinement’ (SC) – self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed ‘live’ numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally required scattering rates disagree with ET predictions by orders of magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvénic streaming are unstable. In live simulations, they revert to either arbitrarily rapid CR escape with zero secondary production, ormore »

2. (Ed.)
Abstract Cosmic rays (CRs) with ∼ GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic γ-ray constraints. Here, we compare the effects of these models relative to simpler “diffusion+streaming” models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend tomore »
3. 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 ∼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 ∼vA by unidirectional Alfvén waves excited by the RDIs, with a uniform distribution of the pitch angle cosine μ and a flat profile of the CR pitch angle diffusion coefficient Dμμ around μ = 0, without the ‘90○ pitch angle problem.’ With CR feedback on the gas included, Dμμ decreases by a factor of a few, indicating a lower CRmore »
4. ABSTRACT

Phenomenological models of cosmic ray (CR) transport in the Milky Way can reproduce a wide range of observations assuming that CRs scatter off of magnetic-field fluctuations with spectrum ∝ k−δ and δ ∼ [1.4, 1.67]. We study the extent to which such models can be reconciled with current microphysical theories of CR transport, specifically self-confinement due to the streaming instability and/or extrinsic turbulence due to a cascade of magnetohydrodynamic (MHD) fast modes. We first review why it is that on their own neither theory is compatible with observations. We then highlight that CR transport is a strong function of local plasma conditions in the multiphase interstellar medium, and may be diffusive due to turbulence in some regions and streaming due to self-confinement in others. A multiphase combination of scattering mechanisms can in principle reproduce the main trends in the proton spectrum and the boron-to-carbon ratio. However, models with a combination of scattering by self-excited waves and fast-mode turbulence require significant fine-tuning due to fast-mode damping, unlike phenomenological models that assume undamped Kolmogorov turbulence. The assumption that fast modes follow a weak cascade is also not well justified theoretically, as the weak cascade is suppressed by wave steepening and weak-shock dissipation evenmore »

5. Abstract

Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as nonresonant reacceleration. It is often invoked to explain the observed ratios of primary to secondary CRs at ∼GeV energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasmaβ. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-βenvironments like the interstellar medium but remain unchanged in high-βenvironments like the intracluster medium. We also quantify the streaming energy-loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution dependent (hence unconverged in large-scale simulations) and heavily suppressed compared to the isotropic loss ratevA· ∇PCR/PCRvA/L0, due to misalignment between the mean field and isotropic CR gradients. Unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength overβ∼ 1–100 and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarilymore »