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-pro
cessing 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.
We introduce the Virgo Consortiums EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolutio
n and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_odot lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* gtrsim 10^9 M_odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.
X-ray observations of circumgalactic coronae provide a valuable means by which to test galaxy formation theories. Two primary mechanisms are thought to be responsible for the establishment of such coronae: accretion of intergalactic gas (IGM) and/or
galactic feedback. In this paper, we first compare our Chandra sample of galactic coronae of 53 nearby highly-inclined disc galaxies to an analytical model considering only the accretion of IGM. We confirm the existing conclusion that this pure accretion model substantially over-predicts the coronal emission. We then select 30 field galaxies from our original sample, and correct their coronal luminosities to uniformly compare them to deep X-ray measurements of several massive disc galaxies from the literature, as well as to a comparable sample of simulated galaxies drawn from the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC). These simulations explicitly model both accretion and SNe feedback and yield galaxies exhibit X-ray properties in broad agreement with our observational sample. However, notable and potentially instructive discrepancies exist between the slope and scatter of the Lx-M200 and Lx-SFR relations, highlighting some known shortcomings of GIMIC, e.g., the absence of AGN feedback, and possibly the adoption of constant stellar feedback parameters. The simulated galaxies exhibit a tight Lx-M200 correlation with little scatter. Having inferred M200 for our observational sample via the Tully-Fisher relation, we find a weaker and more scattered correlation. In the simulated and observed samples alike, massive non-starburst galaxies above a typical transition mass of M*~2e11Msun or M200~1e13Msun tend to have higher Lx/M* and Lx/M200 than low-mass counterparts, indicating that the accretion of IGM plays an increasingly important role in establishing the observable hot circumgalactic medium with increasing galaxy mass.
We study the effect of baryons on the abundance of structures and substructures in a Lambda-CDM cosmology, using a pair of high resolution cosmological simulations from the GIMIC project. Both simulations use identical initial conditions, but while o
ne contains only dark matter, the other also includes baryons. We find that gas pressure, reionisation, supernova feedback, stripping, and truncated accretion systematically reduce the total mass and the abundance of structures below ~10^12 solar masses compared to the pure dark matter simulation. Taking this into account and adopting an appropriate detection threshold lowers the abundance of observed galaxies with maximum circular velocities below 100 km/s, significantly reducing the reported discrepancy between Lambda-CDM and the measured HI velocity function of the ALFALFA survey. We also show that the stellar-to-total mass ratios of galaxies with stellar masses of ~10^5 - 10^7 solar masses inferred from abundance matching of the (sub)halo mass function to the observed galaxy mass function increase by a factor of ~2. In addition, we find that an important fraction of low-mass subhaloes are completely devoid of stars. Accounting for the presence of dark subhaloes below 10^10 solar masses further reduces the abundance of observable objects, and leads to an additional increase in the inferred stellar-to-total mass ratio by factors of 2 - 10 for galaxies in haloes of 10^9 - 10^10 solar masses. This largely reconciles the abundance matching results with the kinematics of individual dwarf galaxies in Lambda-CDM. We propose approximate corrections to the masses of objects derived from pure dark matter calculations to account for baryonic effects.
[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1,
0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of downsizing (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.
