More Like this

1. Tracing the origin of low diffusivity and CR bubbles around sources

   https://doi.org/10.22323/1.395.0163  

   Schroer, Benedikt; Pezzi, Oreste; Caprioli, Damiano; Haggerty, Colby; Blasi, Pasquale (July 2021, 37th International Cosmic Ray Conference (ICRC2021))

   null (Ed.)

   Cosmic rays leave their sources mainly along the local magnetic field present in the region around the source and in doing so they excite both resonant and non-resonant modes through streaming instabilities. The excitation of these modes leads to enhanced particle scattering and in turn to a large pressure gradient that causes the formation of expanding bubbles of gas and self-generated magnetic fields. By means of hybrid particle-in-cell simulations, we demonstrate that, by exciting this instability, cosmic rays excavate a cavity around their source where the diffusivity is strongly suppressed. This phenomenon is general and is expected to occur around any sufficiently powerful cosmic ray source in the Galaxy. Our results are consistent with recent γ-ray observations where emission from the region around supernova remnants and stellar clusters have been used to infer that the diffusion coefficient around these sources is ∼10−100 times smaller than the typical Galactic one. 

   more » « less
2. Effects of Different Cosmic Ray Transport Models on Galaxy Formation

   https://doi.org/10.1093/mnras/staa3692  

   Hopkins, Philip F; Chan, T K; Squire, Jonathan; Quataert, Eliot; Ji, Suoqing; Kereš, Dušan; Faucher-Giguére, Claude-André (January 2021, Monthly Notices of the Royal Astronomical Society)

   null (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 to lie “in between” models which ignore CRs entirely and models which treat CRs with constant diffusivity. 

   more » « less
3. Numerical study of cosmic ray confinement through dust resonant drag instabilities

   https://doi.org/10.1093/mnras/stac895  

   Ji, Suoqing; Squire, Jonathan; Hopkins, Philip F (April 2022, Monthly Notices of the Royal Astronomical Society)

   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 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. 

   more » « less
4. Properties of the circumgalactic medium in cosmic ray-dominated galaxy haloes

   https://doi.org/10.1093/mnras/staa1849  

   Ji, Suoqing; Chan, T K; Hummels, Cameron B; Hopkins, Philip F; Stern, Jonathan; Kereš, Dušan; Quataert, Eliot; Faucher-Giguère, Claude-André; Murray, Norman (August 2020, Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass haloes hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, and collisional and streaming losses, with constant parallel diffusivity $$\kappa \sim 3\times 10^{29}\, \mathrm{cm^2\ s^{-1}}$$ chosen to match γ-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass haloes at z ≲ 1–2. The gas in these ‘CR-dominated’ haloes differs significantly from runs without CRs: the gas is primarily cool (a few $${\sim}10^{4}\,$$ K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ‘low’ and ‘mid’ ions in this diffuse cool gas is dominated by photoionization, with O vi columns $${\gtrsim}10^{14.5}\, \mathrm{cm^{-2}}$$ at distances $${\gtrsim}150\, \mathrm{kpc}$$. CR and thermal gas pressure are locally anticorrelated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same haloes are primarily warm/hot ($$T\gtrsim 10^{5}\,$$K) with thermal pressure balancing gravity, collisional ionization dominates, O vi columns are lower and Ne viii higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase. 

   more » « less
5. A simple sub-grid model for cosmic ray effects on galactic scales

   https://doi.org/10.1093/mnras/stad976  

   Hopkins, Philip F.; Butsky, Iryna S.; Ji, Suoqing; Kereš, Dušan (April 2023, Monthly Notices of the Royal Astronomical Society)

   ABSTRACT Many recent numerical studies have argued that cosmic rays (CRs) from supernovae (SNe) or active galactic nuclei (AGNs) could play a crucial role in galaxy formation, in particular by establishing a CR-pressure-dominated circumgalactic medium (CGM). But explicit CR-magnetohydrodynamics (CR-MHD) remains computationally expensive, and it is not clear whether those results can be applied to simulations that do not explicitly treat magnetic fields or resolved interstellar medium phase structure. We therefore present an intentionally extremely simplified ‘sub-grid’ model for CRs, which attempts to capture the key qualitative behaviors of greatest interest for those interested in simulations or semi-analytical models including some approximate CR effects on galactic (≳ kpc) scales, while imposing negligible computational overhead. The model is numerically akin to some recently developed sub-grid models for radiative feedback, and allows for a simple constant parametrization of the CR diffusivity and/or streaming speed; it allows for an arbitrary distribution of sources (proportional to black hole accretion rates or star–particle SNe rates or gas/galaxy star formation rates), and interpolates between the limits where CRs escape the galaxies with negligible losses and those where CRs lose most of their energy catastrophically before escape (relevant in e.g. starburst galaxies). The numerical equations are solved trivially alongside gravity in most codes. We compare this to explicit CR-MHD simulations and discuss where the (many) sub-grid approximations break down, and what drives the major sources of uncertainty. 

   more » « less