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
The large-scale distribution of ionized metals in IllustrisTNG
We study the intrinsic large-scale distribution and evolution of seven ionized metals in the IllustrisTNG magnetohydrodynamical cosmological simulation. We focus on the fractions of C ii, C iv, Mg ii, N v, Ne viii, O vi, and Si iv in different cosmic web structures (filaments, haloes, and voids) and gas phases (warm–hot intergalactic medium, hot, diffuse, and condensed gas) from z = 6 to z = 0. Our analysis provides a new perspective to the study of the distribution and evolution of baryons across cosmic time while offering new hints in the context of the well-known missing baryons problem. The cosmic web components are here identified using the local comoving dark matter density, which provides a simple but effective way of mapping baryons on large scales. Our results show that C ii and Mg ii are mostly located in condensed gas inside haloes in high-density and low-temperature star-forming regions ($\rho _{\rm gas}/\bar{\rho }_{\rm bar}\gtrsim 10^3$, and T ≲ 105 K). C iv and Si iv present similar evolution of their mass fractions in haloes and filaments across cosmic time. In particular, their mass budgets in haloes in condensed phase ($\rho _{\rm gas}/\bar{\rho }_{\rm bar}\gtrsim 10^3$, and T ≲ 105 K) are driven by gas cooling and star formation with a peak at z ∼ 2. Finally, our results confirm that O vi, Ne viii, and N v are good tracers of warm/hot and low-density gas at low redshift ($\rho _{\rm gas}/\bar{\rho }_{\rm bar}\lesssim 10^3$, and T ≳ 105 K), regions that are likely to contain most of the missing baryons in the local Universe.
more » « less- NSF-PAR ID:
- 10360977
- 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. 399-412
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Intervening C iv absorbers are key tracers of metal-enriched gas in galaxy haloes over cosmic time. Previous studies suggest that the C iv cosmic mass density ($\Omega _{\rm C \, \small {IV}}$) decreases slowly over 1.5 $\lesssim \, z\lesssim$ 5 before declining rapidly at z ≳ 5, but the cause of this downturn is poorly understood. We characterize the $\Omega _{\rm C \, \small {IV}}$ evolution over 4.3 ≲ z ≲ 6.3 using 260 absorbers found in 42 XSHOOTER spectra of z ∼ 6 quasars, of which 30 come from the ESO Large Program XQR-30. The large sample enables us to robustly constrain the rate and timing of the downturn. We find that $\Omega _{\rm C \, \small {IV}}$ decreases by a factor of 4.8 ± 2.0 over the ∼300 Myr interval between z ∼ 4.7 and ∼5.8. The slope of the column density (log N) distribution function does not change, suggesting that C iv absorption is suppressed approximately uniformly across 13.2 ≤ log N/cm−2 < 15.0. Assuming that the carbon content of galaxy haloes evolves as the integral of the cosmic star formation rate density (with some delay due to stellar lifetimes and outflow travel times), we show that chemical evolution alone could plausibly explain the fast decline in $\Omega _{\rm C \, \small {IV}}$ over 4.3 ≲ z ≲ 6.3. However, the C iv/C ii ratio decreases at the highest redshifts, so the accelerated decline in $\Omega _{\rm C \, \small {IV}}$ at z ≳ 5 may be more naturally explained by rapid changes in the gas ionization state driven by evolution of the UV background towards the end of hydrogen reionization.more » « less
-
ABSTRACT We analyse the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z = 0.25 and 2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range $M_{\rm h}(z=0) \sim 10^{10} \!-\! 10^{12} \rm {\,M}_\odot$. We analyse both where the gas ends up after first leaving the CGM (its ‘proximate’ fate) and its location at z = 0 (its ‘ultimate’ fate). Of the CGM at z = 2, about half is found in the ISM or stars of the central galaxy by z = 0 in Mh(z = 2) ∼ 5 × 1011 M⊙ haloes, but most of the CGM in lower mass haloes ends up in the intergalactic medium (IGM). This is so even though most of the CGM in Mh(z = 2) ∼ 5 × 1010 M⊙ haloes first accretes on to the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z = 0.25 remains in the CGM by z = 0 at all halo masses analysed. Of the CGM gas that subsequently accretes on to the central galaxy in the progenitors of Mh(z = 0) ∼ 1012 M⊙ haloes, most of it is cool (T ∼ 104 K) at z = 2 but hot (∼Tvir) at z ∼ 0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z = 0.25 and $2 \, {\gtrsim} 80{{\ \rm per\ cent}}$ of the cool gas in $M_{\rm h} \gtrsim 10^{11} \rm {M}_\odot$ haloes will accrete on to a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM.more » « less
-
null (Ed.)ABSTRACT Observations of emission lines in active galactic nuclei (AGNs) often find fast (∼1000 km s−1) outflows extending to kiloparsec scales, seen in ionized, neutral atomic and molecular gas. In this work we present radiative transfer calculations of emission lines in hydrodynamic simulations of AGN outflows driven by a hot wind bubble, including non-equilibrium chemistry, to explore how these lines trace the physical properties of the multiphase outflow. We find that the hot bubble compresses the line-emitting gas, resulting in higher pressures than in the ambient interstellar medium or that would be produced by the AGN radiation pressure. This implies that observed emission line ratios such as [O iv]$_{25 \, \rm {\mu m}}$ / [Ne ii]$_{12 \, \rm {\mu m}}$, [Ne v]$_{14 \, \rm {\mu m}}$ / [Ne ii]$_{12 \, \rm {\mu m}}$, and [N iii]$_{57 \, \rm {\mu m}}$ / [N ii]$_{122 \, \rm {\mu m}}$ constrain the presence of the bubble and hence the outflow driving mechanism. However, the line-emitting gas is under-pressurized compared to the hot bubble itself, and much of the line emission arises from gas that is out of pressure, thermal and/or chemical equilibrium. Our results thus suggest that assuming equilibrium conditions, as commonly done in AGN line emission models, is not justified if a hot wind bubble is present. We also find that ≳50 per cent of the mass outflow rate, momentum flux, and kinetic energy flux of the outflow are traced by lines such as [N ii]$_{122 \, \rm {\mu m}}$ and [Ne iii]$_{15 \, \rm {\mu m}}$ (produced in the 10$^{4} \, \rm {K}$ phase) and [C ii]$_{158 \, \rm {\mu m}}$ (produced in the transition from 10$^{4} \, \rm {K}$ to 100 K).more » « less
-
more » « less
ABSTRACT We analyse the internal structure and dynamics of cosmic-web filaments connecting massive high-z haloes. Our analysis is based on a high-resolution arepo cosmological simulation zooming-in on three Mpc-scale filaments feeding three massive haloes of $\sim 10^{12}\, \text{M}_\odot$ at z ∼ 4, embedded in a large-scale sheet. Each filament is surrounded by a cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \, {\rm kpc}$. The post-shock gas is in virial equilibrium within the potential well set by an isothermal dark-matter filament. The filament line-mass is $\sim 9\times 10^8\, \text{M}_\odot \, {\rm kpc}^{-1}$, the gas fraction within rshock is the universal baryon fraction, and the virial temperature is ∼7 × 105 K. These all match expectations from analytical models for filament properties as a function of halo mass and redshift. The filament cross-section has three radial zones. In the outer ‘thermal’ (T) zone, $r \ge 0.65 \, r_{\rm shock}$, inward gravity, and ram-pressure forces are overbalanced by outward thermal pressure forces, decelerating the inflowing gas and expanding the shock outwards. In the intermediate ‘vortex’ (V) zone, 0.25 ≤ r/rshock ≤ 0.65, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outward force is dominated by centrifugal forces associated with these vortices, with additional contributions from global rotation and thermal pressure. Shear and turbulent forces associated with the vortices act inwards. The inner ‘stream’ (S) zone, $r \lt 0.25 \, r_{\rm shock}$, is a dense isothermal core, $T\sim 3 \times 10^4 \, {\rm K}$ and $n_{\rm H}\sim 0.01 \, {\rm cm^{-3}}$, defining the cold streams that feed galaxies. The core is formed by an isobaric cooling flow and is associated with a decrease in outward forces, though exhibiting both inflows and outflows.
