skip to main content


Title: Thermal instability of halo gas heated by streaming cosmic rays
ABSTRACT Heating of virialized gas by streaming cosmic rays (CRs) may be energetically important in galaxy haloes, groups, and clusters. We present a linear thermal stability analysis of plasmas heated by streaming CRs. We separately treat equilibria with and without background gradients, and with and without gravity. We include both CR streaming and diffusion along the magnetic-field direction. Thermal stability depends strongly on the ratio of CR pressure to gas pressure, which determines whether modes are isobaric or isochoric. Modes with $\boldsymbol {k \cdot B }\ne 0$ are strongly affected by CR diffusion. When the streaming time is shorter than the CR diffusion time, thermally unstable modes (with $\boldsymbol {k \cdot B }\ne 0$) are waves propagating at a speed ∝ the Alfvén speed. Halo gas in photoionization equilibrium is thermally stable independent of CR pressure, while gas in collisional ionization equilibrium is unstable for physically realistic parameters. In gravitationally stratified plasmas, the oscillation frequency of thermally overstable modes can be higher in the presence of CR streaming than the buoyancy/free-fall frequency. This may modify the critical tcool/tff at which multiphase gas is present. The criterion for convective instability of a stratified, CR-heated medium can be written in the familiar Schwarzschild form dseff/dz < 0, where seff is an effective entropy involving the gas and CR pressures. We discuss the implications of our results for the thermal evolution and multiphase structure of galaxy haloes, groups, and clusters.  more » « less
Award ID(s):
1715054
NSF-PAR ID:
10302583
Author(s) / Creator(s):
 ;  
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
493
Issue:
2
ISSN:
0035-8711
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. ABSTRACT

    We investigate how cosmic rays (CRs) affect thermal and hydrostatic stability of circumgalactic (CGM) gas, in simulations with both CR streaming and diffusion. Local thermal instability can be suppressed by CR-driven entropy mode propagation, in accordance with previous analytic work. However, there is only a narrow parameter regime where this operates, before CRs overheat the background gas. As mass dropout from thermal instability causes the background density and hence plasma β ≡ Pg/PB to fall, the CGM becomes globally unstable. At the cool disc-to-hot−halo interface, a sharp drop in density boosts Alfven speeds and CR gradients, driving a transition from diffusive to streaming transport. CR forces and heating strengthen, while countervailing gravitational forces and radiative cooling weaken, resulting in a loss of both hydrostatic and thermal equilibrium. In lower β haloes, CR heating drives a hot, single-phase diffuse wind with velocities v ∝ (theat/tff)−1, which exceeds the escape velocity when theat/tff ≲ 0.4. In higher β haloes, where the Alfven Mach number is higher, CR forces drive multi-phase winds with cool, dense fountain flows and significant turbulence. These flows are CR dominated due to ‘trapping’ of CRs by weak transverse B-fields, and have the highest mass loading factors. Thus, local thermal instability can result in winds or fountain flows where either the heat or momentum input of CRs dominates.

     
    more » « less
  3. ABSTRACT

    Active Galactic Nuclei (AGN) are believed to provide the energy that prevents runaway cooling of gas in the cores of galaxy clusters. However, how this energy is transported and thermalized throughout the Intracluster Medium (ICM) remains unclear. In recent work, we showed that streaming cosmic rays (CRs) destabilize sound waves in dilute ICM plasmas. Here, we show that CR streaming in the presence of gravity also destabilizes a pressure-balanced wave. We term this new instability the CR buoyancy instability (CRBI). In stark contrast to standard results without CRs, the pressure-balanced mode is highly compressible at short wavelengths due to CR streaming. Maximal growth rates are of order (pc/pg)β1/2ωff, where pc/pg is the ratio of CR pressure to thermal gas pressure, β is the ratio of thermal to magnetic pressure, and ωff is the free-fall frequency. The CRBI operates alongside buoyancy instabilities driven by background heat fluxes, i.e. the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI). When the thermal mean free path lmfp is ≪ the gas scale height H, the HBI/MTI set the growth rate on large scales, while the CRBI sets the growth rate on small scales. Conversely, when lmfp ∼ H and (pc/pg)β1/2 ≳ 1, CRBI growth rates exceed HBI/MTI growth rates even on large scales. Our results suggest that CR-driven instabilities may be partially responsible for the sound waves/weak shocks and turbulence observed in galaxy clusters. CR-driven instabilities generated near radio bubbles may also play an important role redistributing AGN energy throughout clusters.

     
    more » « less
  4. null (Ed.)
    ABSTRACT We study the impact of cosmic rays (CRs) on the structure of virial shocks, using a large suite of high-resolution cosmological FIRE-2 simulations accounting for CR injection by supernovae. In Milky Way-mass, low-redshift (z ≲ 1−2) haloes, which are expected to form ‘hot haloes’ with slowly cooling gas in quasi-hydrostatic equilibrium (with a stable virial shock), our simulations without CRs do exhibit clear virial shocks. The cooler phase condensing out from inflows becomes pressure confined to overdense clumps, embedded in low-density, volume-filling hot gas with volume-weighted cooling time longer than inflow time. The gas thus transitions sharply from cool free-falling inflow, to hot and thermal-pressure supported at approximately the virial radius (≈Rvir), and the shock is quasi-spherical. With CRs, we previously argued that haloes in this particular mass and redshift range build up CR-pressure-dominated gaseous haloes. Here, we show that when CR pressure dominates over thermal pressure, there is no significant virial shock. Instead, inflowing gas is gradually decelerated by the CR pressure gradient and the gas is relatively subsonic out to and even beyond Rvir. Rapid cooling also maintains subvirial temperatures in the inflowing gas within ∼Rvir. 
    more » « less
  5. Abstract We study the propagation of mildly relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta et al. to postprocess three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfvén speed, and diffusion relative to the Alfvén waves, with the diffusion coefficient set by the balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the postprocessed simulations to make predictions regarding the potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extraplanar region (∣ z ∣ > 0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR. 
    more » « less