Recent theoretical studies predict that the circumgalactic medium (CGM) around low-redshift, ∼L* galaxies could have substantial non-thermal pressure support in the form of cosmic rays. However, these predictions are sensitive to the specific model of cosmic ray transport employed, which is theoretically and observationally underconstrained. In this work, we propose a novel observational constraint for calculating the lower limit of the radially averaged, effective cosmic ray transport rate, ${\kappa _{\rm eff}^{\rm min}}$. Under a wide range of assumptions (so long as cosmic rays do not lose a significant fraction of their energy in the galactic disc, regardless of whether the cosmic ray pressure is important or not in the CGM), we demonstrate a well-defined relationship between ${\kappa _{\rm eff}^{\rm min}}$ and three observable galaxy properties: the total hydrogen column density, the average star formation rate, and the gas circular velocity. We use a suite of Feedback in Realistic Environments 2 galaxy simulations with a variety of cosmic ray transport physics to demonstrate that our analytical model of ${\kappa _{\rm eff}^{\rm min}}$ is a robust lower limit of the true cosmic ray transport rate. We then apply our new model to calculate ${\kappa _{\rm eff}^{\rm min}}$ for galaxies in the COS-Halos sample, and confirm this already reveals strong evidence for an effective transport rate that rises rapidly away from the interstellar medium to values ${\kappa _{\rm eff}^{\rm min}}\gtrsim 10^{30\!-\!31}\, {\rm cm}^2\, {\rm s}^{-1}$ (corresponding to anisotropic streaming velocities of $v^{\rm stream}_{\rm eff} \gtrsim 1000\, {\rm km}\, {\rm s}^{-1}$) in the diffuse CGM, at impact parameters larger than 50–100 kpc. We discuss how future observations can provide qualitatively new constraints in our understanding of cosmic rays in the CGM and intergalactic medium.
This content will become publicly available on February 15, 2025
The relevance of some galactic feedback mechanisms, in particular cosmic-ray (CR) feedback and the hydrogen ionizing radiation field, has been challenging to definitively describe in a galactic context, especially far outside the galaxy in the circumgalactic medium (CGM). Theoretical and observational uncertainties prevent conclusive interpretations of multiphase CGM properties derived from ultraviolet (UV) diagnostics. We conduct three-dimensional magnetohydrodynamic simulations of a section of a galactic disc with star formation and feedback, including radiative heating from stars, a UV background, and CR feedback. We utilize the temperature phases present in our simulations to generate Cloudy models to derive spatially and temporally varying synthetic UV diagnostics. We find that radiative effects without additional heating mechanisms are not able to produce synthetic diagnostics in the observed ranges. For low CR diffusivity $\kappa _{\rm {cr}}=10^{28} \rm {cm}^2 \rm {s}^{-1}$, CR streaming heating in the outflow helps our synthetic line ratios roughly match observed ranges by producing transitional temperature gas (T ∼ 105–106 K). High CR diffusivity $\kappa _{\rm {cr}}=10^{29} \rm {cm}^2 \rm {s}^{-1}$, with or without CR streaming heating, produced transitional temperature gas. The key parameter controlling the production of this gas phase remains unclear, as the different star formation history and outflow evolution itself influences these diagnostics. Our work demonstrates the use of UV plasma diagnostics to differentiate between galactic/circumgalactic feedback models.
more » « less- Award ID(s):
- 2044303
- NSF-PAR ID:
- 10516756
- Publisher / Repository:
- MNRAS
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 528
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 7543 to 7563
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT We study the properties of cosmic-ray (CR) driven galactic winds from the warm interstellar medium using idealized spherically symmetric time-dependent simulations. The key ingredients in the model are radiative cooling and CR-streaming-mediated heating of the gas. Cooling and CR heating balance near the base of the wind, but this equilibrium is thermally unstable, leading to a multiphase wind with large fluctuations in density and temperature. In most of our simulations, the heating eventually overwhelms cooling, leading to a rapid increase in temperature and a thermally driven wind; the exception to this is in galaxies with the shallowest potentials, which produce nearly isothermal $T \approx 10^4\,$ K winds driven by CR pressure. Many of the time-averaged wind solutions found here have a remarkable critical point structure, with two critical points. Scaled to real galaxies, we find mass outflow rates $\dot{M}$ somewhat larger than the observed star-formation rate in low-mass galaxies, and an approximately ‘energy-like’ scaling $\dot{M} \propto v_{\rm esc}^{-2}$. The winds accelerate slowly and reach asymptotic wind speeds of only ∼0.4vesc. The total wind power is $\sim 1~{{\ \rm per\ cent}}$ of the power from supernovae, suggesting inefficient preventive CR feedback for the physical conditions modelled here. We predict significant spatially extended emission and absorption lines from 104–105.5 K gas; this may correspond to extraplanar diffuse ionized gas seen in star-forming galaxies.
-
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
-
more » « less
ABSTRACT Many recent numerical studies have argued that cosmic rays (CRs) from supernovae (SNe) or active galactic nuclei (AGNs) could play a crucial role in galaxy formation, in particular by establishing a CR-pressure-dominated circumgalactic medium (CGM). But explicit CR-magnetohydrodynamics (CR-MHD) remains computationally expensive, and it is not clear whether those results can be applied to simulations that do not explicitly treat magnetic fields or resolved interstellar medium phase structure. We therefore present an intentionally extremely simplified ‘sub-grid’ model for CRs, which attempts to capture the key qualitative behaviors of greatest interest for those interested in simulations or semi-analytical models including some approximate CR effects on galactic (≳ kpc) scales, while imposing negligible computational overhead. The model is numerically akin to some recently developed sub-grid models for radiative feedback, and allows for a simple constant parametrization of the CR diffusivity and/or streaming speed; it allows for an arbitrary distribution of sources (proportional to black hole accretion rates or star–particle SNe rates or gas/galaxy star formation rates), and interpolates between the limits where CRs escape the galaxies with negligible losses and those where CRs lose most of their energy catastrophically before escape (relevant in e.g. starburst galaxies). The numerical equations are solved trivially alongside gravity in most codes. We compare this to explicit CR-MHD simulations and discuss where the (many) sub-grid approximations break down, and what drives the major sources of uncertainty.
-
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
