Cosmic rays (CRs) are an important component in the interstellar medium, but their effect on the dynamics of the disc–halo interface (<10 kpc from the disc) is still unclear. We study the influence of CRs on the gas above the disc with high-resolution FIRE-2 cosmological simulations of late-type L⋆ galaxies at redshift z ∼ 0. We compare runs with and without CR feedback (with constant anisotropic diffusion κ∥ ∼ 3 × 1029 cm2 s−1 and streaming). Our simulations capture the relevant disc–halo interactions, including outflows, inflows, and galactic fountains. Extra-planar gas in all of the runs satisfies dynamical balance, where total pressure balances the weight of the overlying gas. While the kinetic pressure from non-uniform motion (≳1 kpc scale) dominates in the mid-plane, thermal and bulk pressures (or CR pressure if included) take over at large heights. We find that with CR feedback, (1) the warm (∼104 K) gas is slowly accelerated by CRs; (2) the hot (>5 × 105 K) gas scale height is suppressed; (3) the warm-hot (2 × 104–5 × 105 K) medium becomes the most volume-filling phase in the disc–halo interface. We develop a novel conceptual model of the near-disc gas dynamics in low-redshift L⋆ galaxies: with CRs, the disc–halo interface is filled with CR-driven warm winds and hot superbubbles that are propagating into the circumgalactic medium with a small fraction falling back to the disc. Without CRs, most outflows from hot superbubbles are trapped by the existing hot halo and gravity, so typically they form galactic fountains.
Cosmic-Ray Driven Outflows to Mpc Scales from L * Galaxies
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
- NSF-PAR ID:
- 10278900
- 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
-
more » « less
ABSTRACT -
null (Ed.)ABSTRACT Many phenomenologically successful cosmological simulations employ kinetic winds to model galactic outflows. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ gadget-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed and the scaling of the mass-loading factor with velocity dispersion. We design our fiducial model to reproduce the measured wind properties at 25 per cent of the virial radius from the Feedback In Realistic Environments simulations. We find that a strong dependence of η ∼ σ5 in low-mass haloes with $\sigma \lt 106\mathrm{\, km\, s^{-1}}$ is required to match the faint end of the stellar mass functions at $z$ > 1. In addition, faster winds significantly reduce wind recycling and heat more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z$ = 2 and 0 within a single model, suggesting that an additional feedback source such as active galactic nucleus might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc}\, h^{-1}$, 2 × 5763 simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z$ ∼ 0–2.more » « less
-
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
-
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
-
ABSTRACT The discovery of giant quasar Ly α nebulae at z > 2 has opened up the possibility to directly study in emission the Circumgalactic and Intergalactic Medium (CGM/IGM). However, the resonant nature of the Ly α line and its different emission mechanisms hamper the ability to constrain both the kinematics and physical properties of the CGM/IGM. Here, we present results of a pilot project aiming at the detection of CGM H α emission, a line which does not suffer from these limitations. To this end, we first used KCWI to detect Ly α emission around three bright quasars with 2.25 < z < 2.27, a range which is free from bright IR sky lines for H α, and then selected the most extended nebula for H α follow-up with MOSFIRE. Within the MOSFIRE slit, we detected H α emission extending up to 20 physical kpc with a total H α flux of FH α = (9.5 ± 0.9) × 10$^{-18}~\mathrm{erg\, s^{-1}\, cm^{-2}}$. Considering the Ly α flux in the same region, we found FLy α/FH α = 3.7 ± 0.3 consistent with that obtained for the Slug Nebula at z = 2.275 and with recombination radiation. This implies high densities or a very broad density distribution within the CGM of high-redshift quasars. Moreover, the H α line profile suggests the presence of multiple emitting components overlapping along our line of sight and relatively quiescent kinematics, which seems incompatible with either quasar outflows capable of escaping the potential well of the host halo or disc-like rotation in a massive halo (>1012 M⊙).more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa3690
