ABSTRACT Recently, cosmic rays (CRs) have emerged as a leading candidate for driving galactic winds. Small-scale processes can dramatically affect global wind properties. We run two-moment simulations of CR streaming to study how sound waves are driven unstable by phase-shifted CR forces and CR heating. We verify linear theory growth rates. As the sound waves grow non-linear, they steepen into a quasi-periodic series of propagating shocks; the density jumps at shocks create CR bottlenecks. The depth of a propagating bottleneck depends on both the density jump and its velocity; ΔPc is smaller for rapidly moving bottlenecks. A series of bottlenecks creates a CR staircase structure, which can be understood from a convex hull construction. The system reaches a steady state between growth of new perturbations, and stair mergers. CRs are decoupled at plateaus, but exert intense forces and heating at stair jumps. The absence of CR heating at plateaus leads to cooling, strong gas pressure gradients and further shocks. If bottlenecks are stationary, they can drastically modify global flows; if their propagation times are comparable to dynamical times, their effects on global momentum and energy transfer are modest. The CR acoustic instability is likely relevant in thermal interfaces between cold and hot gas, as well as galactic winds. Similar to increased opacity in radiative flows, the build-up of CR pressure due to bottlenecks can significantly increase mass outflow rates, by up to an order of magnitude. It seeds unusual forms of thermal instability, and the shocks could have distinct observational signatures, on ∼kpc scales.
more »
« less
The impact of Cosmic Rays on thermal and hydrostatic stability in galactic halos
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 disk 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 β halos, 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 β halos, 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
- Award ID(s):
- 1911198
- PAR ID:
- 10465358
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- ISSN:
- 0035-8711
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
null (Ed.)Abstract We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive (Mhalo ≳ 1011 M⊙), low-redshift (z ≲ 1 − 2) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos “traps” galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows “escape” and CR pressure gradients continuously accelerate this material well into the IGM in “fast” outflows, while lower-density gas at large radii is accelerated in-situ into “slow” outflows that extend to >Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to >Mpc. The CR-driven outflows are primarily cool (T ∼ 105 K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses (≲ 1011 M⊙) or higher redshifts (z ≳ 1 − 2), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena. We caution that these predictions may depend on uncertain CR transport physics.more » « less
-
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
-
null (Ed.)ABSTRACT We examine the thermodynamic state and cooling of the low-z circumgalactic medium (CGM) in five FIRE-2 galaxy formation simulations of Milky Way-mass galaxies. We find that the CGM in these simulations is generally multiphase and dynamic, with a wide spectrum of largely non-linear density perturbations sourced by the accretion of gas from the intergalactic medium (IGM) and outflows from both the central and satellite galaxies. We investigate the origin of the multiphase structure of the CGM with a particle-tracking analysis and find that most of the low-entropy gas has cooled from the hot halo as a result of thermal instability triggered by these perturbations. The ratio of cooling to free-fall time-scales tcool/tff in the hot component of the CGM spans a wide range of ∼1−100 at a given radius but exhibits approximately constant median values of ∼5−20 at all radii 0.1Rvir < r < Rvir. These are similar to the ≈10−20 value typically adopted as the thermal instability threshold in ‘precipitation’ models of the ICM. Consequently, a one-dimensional model based on the assumption of a constant tcool/tff and hydrostatic equilibrium approximately reproduces the number density and entropy profiles of each simulation but only if it assumes the metallicity profile and temperature boundary condition taken directly from the simulation. We explicitly show that the tcool/tff value of a gas parcel in the hot component of the CGM does not predict its probability of subsequently accreting on to the central galaxy. This suggests that the value of tcool/tff is a poor predictor of thermal stability in gaseous haloes in which large-amplitude density perturbations are prevalent.more » « less
An official website of the United States government

