No Arabic abstract
The inflow of cosmological gas onto haloes, while challenging to directly observe and quantify, plays a fundamental role in the baryon cycle of galaxies. Using the EAGLE suite of hydrodynamical simulations, we present a thorough exploration of the physical properties of gas accreting onto haloes -- namely, its spatial characteristics, density, temperature, and metallicity. Classifying accretion as ``hot or `` cold based on a temperature cut of $10^{5.5}{rm K}$, we find that the covering fraction ($f_{rm cov}$) of cold-mode accreting gas is significantly lower than the hot-mode, with $z=0$ $f_{rm cov}$ values of $approx 50%$ and $approx 80%$ respectively. Active Galactic Nuclei (AGN) feedback in EAGLE reduces inflow $f_{rm cov}$ values by $approx 10%$, with outflows decreasing the solid angle available for accretion flows. Classifying inflow by particle history, we find that gas on first-infall onto a halo is metal-depleted by $approx 2$~dex compared to pre-processed gas, which we find to mimic the circum-galactic medium (CGM) in terms of metal content. We also show that high (low) halo-scale gas accretion rates are associated with metal-poor (rich) CGM in haloes below $10^{12}M_{odot}$, and that variation in halo-scale gas accretion rates may offer a physical explanation for the enhanced scatter in the star-forming main sequence at low ($lesssim10^{9}M_{odot}$) and high ($gtrsim10^{10}M_{odot}$) stellar masses. Our results highlight how gas inflow influences several halo- and galaxy-scale properties, and the need to combine kinematic and chemical data in order to confidently break the degeneracy between accreting and outgoing gas in CGM observations.
We use the EAGLE cosmological, hydrodynamical simulations to predict the column density and equivalent width distributions of intergalactic O VII ($E=574$ eV) and O VIII ($E=654$ eV) absorbers at low redshift. These two ions are predicted to account for 40% of the gas-phase oxygen, which implies that they are key tracers of cosmic metals. We find that their column density distributions evolve little at observable column densities from redshift 1 to 0, and that they are sensitive to AGN feedback, which strongly reduces the number of strong (column density $N gtrsim 10^{16} , mathrm{cm}^{-2})$ absorbers. The distributions have a break at $N sim 10^{16} , mathrm{cm}^{-2}$, corresponding to overdensities of $sim 10^{2}$, likely caused by the transition from sheet/filament to halo gas. Absorption systems with $N gtrsim 10^{16} mathrm{cm}^{-2}$ are dominated by collisionally ionized O VII and O VIII, while the ionization state of oxygen at lower column densities is also influenced by photoionization. At these high column densities, O VII and O VIII arising in the same structures probe systematically different gas temperatures, meaning their line ratio does not translate into a simple estimate of temperature. While O VII and O VIII column densities and covering fractions correlate poorly with the H I column density at $N_{mathrm{H , I}} gtrsim 10^{15} , mathrm{cm}^{-2}$, O VII and O VIII column densities are higher in this regime than at the more common, lower H I column densities. The column densities of O VI and especially Ne VIII, which have strong absorption lines in the UV, are good predictors of the strengths of O VII and O VIII absorption and can hence aid in the detection of the X-ray lines.
Determining the spatial distribution and intrinsic physical properties of neutral hydrogen on cosmological scales is one of the key goals of next-generation radio surveys. We use the EAGLE galaxy formation simulations to assess the properties of damped Lyman-alpha absorbers (DLAs) that are associated with galaxies and their underlying dark matter haloes between 0 $leq$ z $leq$ 2. We find that the covering fraction of DLAs increases at higher redshift; a significant fraction of neutral atomic hydrogen (HI) resides in the outskirts of galaxies with stellar mass greater than or equal to 10$^{10}$ M$_odot$; and the covering fraction of DLAs in the circumgalactic medium (CGM) is enhanced relative to that of the interstellar medium (ISM) with increasing halo mass. Moreover, we find that the mean density of the HI in galaxies increases with increasing stellar mass, while the DLAs in high- and low-halo-mass systems have higher column densities than those in galaxies with intermediate halo masses (~ 10$^{12}$ M$_odot$ at z = 0). These high-impact CGM DLAs in high-stellar-mass systems tend to be metal-poor, likely tracing smooth accretion. Overall, our results point to the CGM playing an important role in DLA studies at high redshift (z $geq$ 1). However, their properties are impacted both by numerical resolution and the detailed feedback prescriptions employed in cosmological simulations, particularly that of AGN.
We use the eagle simulations to study the connection between the quenching timescale, $tau_{rm Q}$, and the physical mechanisms that transform star-forming galaxies into passive galaxies. By quantifying $tau_{rm Q}$ in two complementary ways - as the time over which (i) galaxies traverse the green valley on the colour-mass diagram, or (ii) leave the main sequence of star formation and subsequently arrive on the passive cloud in specific star formation rate (SSFR)-mass space - we find that the $tau_{rm Q}$ distribution of high-mass centrals, low-mass centrals and satellites are divergent. In the low stellar mass regime where $M_{star}<10^{9.6}M_{odot}$, centrals exhibit systematically longer quenching timescales than satellites ($approx 4$~Gyr compared to $approx 2$~Gyr). Satellites with low stellar mass relative to their halo mass cause this disparity, with ram pressure stripping quenching these galaxies rapidly. Low mass centrals are quenched as a result of stellar feedback, associated with long $tau_{rm Q}gtrsim 3$~Gyr. At intermediate stellar masses where $10^{9.7},rm M_{odot}<M_{star}<10^{10.3},rm M_{odot}$, $tau_{rm Q}$ are the longest for both centrals and satellites, particularly for galaxies with higher gas fractions. At $M_{star}gtrsim 10^{10.3},rm M_{odot}$, galaxy merger counts and black hole activity increase steeply for all galaxies. Quenching timescales for centrals and satellites decrease with stellar mass in this regime to $tau_{rm Q}lesssim2$~Gyr. In anticipation of new intermediate redshift observational galaxy surveys, we analyse the passive and star-forming fractions of galaxies across redshift, and find that the $tau_{rm Q}$ peak at intermediate stellar masses is responsible for a peak (inflection point) in the fraction of green valley central (satellite) galaxies at $zapprox 0.5-0.7$.
We study the $zapprox3.5$ intergalactic medium (IGM) by comparing new, high-quality absorption spectra of eight QSOs with $langle z_{rm QSO} rangle=3.75$, to virtual observations of the EAGLE cosmological hydrodynamical simulations. We employ the pixel optical depth method and uncover strong correlations between various combinations of HI, CIII, CIV, SiIII, SiIV, and OVI. We find good agreement between many of the simulated and observed correlations, including OVI(HI). However, the observed median optical depths for the CIV(HI) and SiIV(HI) relations are higher than those measured from the mock spectra. The discrepancy increases from up to $approx0.1$ dex at $tau_{rm HI}=1$ to $approx1$ dex at $tau_{rm HI}=10^2$, where we are likely probing dense regions at small galactocentric distances. As possible solutions, we invoke (a) models of ionizing radiation softened above 4 Ryd to account for delayed completion of HeII reionization; (b) simulations run at a higher resolution; (c) the inclusion of additional line broadening due to unresolved turbulence; and (d) increased elemental abundancess; however, none of these factors can fully explain the observed differences. Enhanced photoionization of HI by local sources, which was not modelled, could offer a solution. However, the much better agreement with the observed OVI(HI) relation, which we find probes a hot and likely collisionally-ionized gas phase, indicates that the simulations are not in tension with the hot phase of the IGM, and suggests that the simulated outflows may entrain insufficient cool gas.
We investigate the abundance of galactic molecular hydrogen (H$_2$) in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological hydrodynamic simulations. We assign H$_2$ masses to gas particles in the simulations in post-processing using two different prescriptions that depend on the local dust-to-gas ratio and the interstellar radiation field. Both result in H$_2$ galaxy mass functions that agree well with observations in the local and high-redshift Universe. The simulations reproduce the observed scaling relations between the mass of H$_2$ and the stellar mass, star formation rate and stellar surface density. Towards high edshifts, galaxies in the simulations display larger H$_2$ mass fractions, and correspondingly lower H$_2$ depletion timescales, also in good agreement with observations. The comoving mass density of H$_2$ in units of the critical density, $Omega_{rm H_2}$, peaks at $zapprox 1.2-1.5$, later than the predicted peak of the cosmic star formation rate activity, at $zapprox 2$. This difference stems from the decrease in gas metallicity and increase in interstellar radiation field with redshift, both of which hamper H$_2$ formation. We find that the cosmic H$_2$ budget is dominated by galaxies with $M_{rm H_2}>10^9,rm M_{odot}$, star formation rates $>10,rm M_{odot},rm yr^{-1}$ and stellar masses $M_{rm stellar}>10^{10},rm M_{odot}$, which are readily observable in the optical and near-IR. The match between the H$_2$ properties of galaxies that emerge in the simulations and observations is remarkable, particularly since H$_2$ observations were not used to adjust parameters in EAGLE.