We use analytical calculations and time-dependent spherically symmetric simulations to study the properties of isothermal galactic winds driven by cosmic rays (CRs) streaming at the Alfvén velocity. The simulations produce time-dependent flows permeated by strong shocks; we identify a new linear instability of sound waves that sources these shocks. The shocks substantially modify the wind dynamics, invalidating previous steady state models: the CR pressure pc has a staircase-like structure with dpc/dr ≃ 0 in most of the volume, and the time-averaged CR energetics are in many cases better approximated by pc ∝ ρ1/2, rather than the canonical pc ∝ ρ2/3. Accounting for this change in CR energetics, we analytically derive new expressions for the mass-loss rate, momentum flux, wind speed, and wind kinetic power in galactic winds driven by CR streaming. We show that streaming CRs are ineffective at directly driving cold gas out of galaxies, though CR-driven winds in hotter ISM phases may entrain cool gas. For the same physical conditions, diffusive CR transport (Paper I) yields mass-loss rates that are a few-100 times larger than streaming transport, and asymptotic wind powers that are a factor of ≃4 larger. We discuss the implications of our results for galactic wind theory and observations; strong shocks driven by CR-streaming-induced instabilities produce gas with a wide range of densities and temperatures, consistent with the multiphase nature of observed winds. We also quantify the applicability of the isothermal gas approximation for modelling streaming CRs and highlight the need for calculations with more realistic thermodynamics.
Large-scale galactic winds driven by stellar feedback are one phenomenon that influences the dynamical and chemical evolution of a galaxy, redistributing material throughout the circumgalatic medium. Non-thermal feedback from galactic cosmic rays (CRs) – high-energy charged particles accelerated in supernovae and young stars – can impact the efficiency of wind driving. The streaming instability limits the speed at which they can escape. However, in the presence of turbulence, the streaming instability is subject to suppression that depends on the magnetization of turbulence given by its Alfvén Mach number. While previous simulations that relied on a simplified model of CR transport have shown that super-Alfvénic streaming of CRs enhances galactic winds, in this paper we take into account a realistic model of streaming suppression. We perform three-dimensional magnetohydrodynamic simulations of a section of a galactic disc and find that turbulent damping dependent on local magnetization of turbulent interstellar medium (ISM) leads to more spatially extended gas and CR distributions compared to the earlier streaming calculations, and that scale heights of these distributions increase for stronger turbulence. Our results indicate that the star formation rate increases with the level of turbulence in the ISM. We also find that the instantaneous wind mass loading is sensitive to local streaming physics with the mass loading dropping significantly as the strength of turbulence increases.
more » « less- Award ID(s):
- 1715140
- NSF-PAR ID:
- 10121091
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 490
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 1271-1282
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
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.
-
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
-
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
-
more » « less
ABSTRACT The physics of Cosmic ray (CR) transport remains a key uncertainty in assessing whether CRs can produce galaxy-scale outflows consistent with observations. In this paper, we elucidate the physics of CR-driven galactic winds for CR transport dominated by diffusion. A companion paper considers CR streaming. We use analytic estimates validated by time-dependent spherically symmetric simulations to derive expressions for the mass-loss rate, momentum flux, and speed of CR-driven galactic winds, suitable for cosmological-scale or semi-analytic models of galaxy formation. For CR diffusion coefficients κ ≳ r0ci, where r0 is the base radius of the wind and ci is the isothermal gas sound speed, the asymptotic wind energy flux is comparable to that supplied to CRs, and the outflow rapidly accelerates to supersonic speeds. By contrast, for κ ≲ r0ci, CR-driven winds accelerate more slowly and lose most of their energy to gravity, a CR analogue of photon-tired stellar winds. Given CR diffusion coefficients estimated using Fermi gamma-ray observations of pion decay, we predict mass-loss rates in CR-driven galactic winds of the order of the star formation rate for dwarf and disc galaxies. The dwarf galaxy mass-loss rates are small compared to the mass-loadings needed to reconcile the stellar and dark matter halo mass functions. For nuclear starbursts (e.g. M82, Arp 220), CR diffusion and pion losses suppress the CR pressure in the galaxy and the strength of CR-driven winds. We discuss the implications of our results for interpreting observations of galactic winds and for the role of CRs in galaxy formation.
