skip to main content


Title: The physics of galactic winds driven by cosmic rays I: Diffusion
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.

 
more » « less
NSF-PAR ID:
10360865
Author(s) / Creator(s):
; ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
510
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
p. 1184-1203
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    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.

     
    more » « less
  2. ABSTRACT

    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
  3. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant. 
    more » « less
  4. 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
  5. Abstract We present the implementation and the first results of cosmic ray (CR) feedback in the Feedback In Realistic Environments (FIRE) simulations. We investigate CR feedback in non-cosmological simulations of dwarf, sub-L⋆ starburst, and L⋆ galaxies with different propagation models, including advection, isotropic and anisotropic diffusion, and streaming along field lines with different transport coefficients. We simulate CR diffusion and streaming simultaneously in galaxies with high resolution, using a two moment method. We forward-model and compare to observations of γ-ray emission from nearby and starburst galaxies. We reproduce the γ-ray observations of dwarf and L⋆ galaxies with constant isotropic diffusion coefficient κ ∼ 3 × 1029 cm2 s−1. Advection-only and streaming-only models produce order-of-magnitude too large γ-ray luminosities in dwarf and L⋆ galaxies. We show that in models that match the γ-ray observations, most CRs escape low-gas-density galaxies (e.g. dwarfs) before significant collisional losses, while starburst galaxies are CR proton calorimeters. While adiabatic losses can be significant, they occur only after CRs escape galaxies, so they are only of secondary importance for γ-ray emissivities. Models where CRs are “trapped” in the star-forming disk have lower star formation efficiency, but these models are ruled out by γ-ray observations. For models with constant κ that match the γ-ray observations, CRs form extended halos with scale heights of several kpc to several tens of kpc. 
    more » « less