The metallicity of star-forming gas in galaxies from the EAGLE simulations increases with stellar mass. Here we investigate whether the scatter around this relation correlates with morphology and/or stellar kinematics. At redshift $z=0$, galaxies with more rotational support have lower metallicities on average when the stellar mass is below $M_\star\approx 10^{10}~{\rm M}_\odot$. This trend inverts at higher values of $M_\star$, when prolate galaxies show typically lower metallicity. At increasing redshifts, the trend between rotational support and metallicity becomes weaker at low stellar mass but more pronounced at high stellar mass. We argue that the secondary dependence of metallicity on stellar kinematics is another manifestation of the observed anti-correlation between metallicity and star formation rate at a given stellar mass. At low masses, such trends seem to be driven by the different star-formation histories of galaxies and stellar feedback. At high masses, feedback from active galactic nuclei and galaxy mergers play a dominant role.

The experiment to detect the global epoch of reionization signature (EDGES) collaboration reported the detection of a line at 78MHz in the sky-averaged spectrum due to neutral hydrogen (HI) 21-cm hyperfine absorption of cosmic microwave background photons at $z\sim 17$. This requires that the spin temperature of HI be coupled to the kinetic temperature of the gas at this $z$ through the scattering of Ly-$\alpha$ photons emitted by massive stars. To explain the experimental result, star formation needs to be sufficiently efficient at $z\sim 17$ and this can be used to constrain models in which small-scale structure formation is suppressed (DMF models), either due to dark matter free-streaming or non-standard inflationary dynamics. We combine simulations of structure formation with a simple recipe for star formation to investigate whether these models emit enough Ly-$\alpha$ photons to reproduce the experimental signal for reasonable values of the star formation efficiency, $f_\star$. We find that a thermal warm dark matter (WDM) model with mass $m_\mathrm{WDM}\sim 4.3\,\mathrm{keV}$ is consistent with the timing of the signal for $f_\star\lesssim 2\%$. The exponential growth of structure around $z\sim 17$ in such a model naturally generates a sharp onset of the absorption. A warmer model with $m_\mathrm{WDM}\sim 3\,\mathrm{keV}$ requires a higher star formation efficiency, $f_\star\sim 6\%$, which is a factor of few above predictions of current star formation models and observations of satellites in the Milky Way. However, uncertainties in the process of star formation at these $z$ do not allow to derive strong constrains on such models using 21-cm absorption line. The onset of the 21-cm absorption is generally faster in DMF compared to cold dark matter (CDM), unless some process significantly suppresses star formation in halos with masses below $\sim 10^8\,h^{-1}\,\mathrm{M}_\odot$.

We present infrared luminosity functions and dust mass functions for the EAGLE cosmological simulation, based on synthetic multi-wavelength observations generated with the SKIRT radiative transfer code. In the local Universe, we reproduce the observed infrared luminosity and dust mass functions very well. Some minor discrepancies are encountered, mainly in the high luminosity regime, where the EAGLE-SKIRT luminosity functions mildly but systematically underestimate the observed ones. The agreement between the EAGLE-SKIRT infrared luminosity functions and the observed ones gradually worsens with increasing lookback time. Fitting modified Schechter functions to the EAGLE-SKIRT luminosity and dust mass functions at different redshifts up to $z=1$, we find that the evolution is compatible with pure luminosity/mass evolution. The evolution is relatively mild: within this redshift range, we find an evolution of $L_{\star,250}\propto(1+z)^{1.68}$, $L_{\star,\text{TIR}}\propto(1+z)^{2.51}$ and $M_{\star,\text{dust}}\propto(1+z)^{0.83}$ for the characteristic luminosity/mass. For the luminosity/mass density we find $\varepsilon_{250}\propto(1+z)^{1.62}$, $\varepsilon_{\text{TIR}}\propto(1+z)^{2.35}$ and $\rho_{\text{dust}}\propto(1+z)^{0.80}$, respectively. The mild evolution of the dust mass density is in relatively good agreement with observations, but the slow evolution of the infrared luminosity underestimates the observed luminosity evolution significantly. We argue that these differences can be attributed to increasing limitations in the radiative transfer treatment due to increasingly poorer resolution, combined with a slower than observed evolution of the SFR density in the EAGLE simulation and the lack of AGN emission in our EAGLE-SKIRT post-processing recipe.

We introduce the ARTEMIS simulations, a new set of 42 zoomed-in, high-resolution (baryon particle mass of ~ 2x10^4 Msun/h), hydrodynamical simulations of galaxies residing in haloes of Milky Way mass, simulated with the EAGLE galaxy formation code with re-calibrated stellar feedback. In this study, we analyse the structure of stellar haloes, specifically the mass density, surface brightness, metallicity, colour and age radial profiles, finding generally very good agreement with recent observations of local galaxies. The stellar density profiles are well fitted by broken power laws, with inner slopes of ~ -3, outer slopes of ~ -4 and break radii that are typically ~ 20-40 kpc. The break radii generally mark the transition between in situ formation and accretion-driven formation of the halo. The metallicity, colour and age profiles show mild large-scale gradients, particularly when spherically-averaged or viewed along the major axes. Along the minor axes, however, the profiles are nearly flat, in agreement with observations. Overall, the structural properties can be understood by two factors: that in situ stars dominate the inner regions and that they reside in a spatially-flattened distribution that is aligned with the disc. Observations targeting both the major and minor axes of galaxies are thus required to obtain a complete picture of stellar haloes.

We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$\alpha$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$\alpha$ star formation rate surface density $\Sigma_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. S\'ersic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium.

There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \mu m}$, MIPS - 24, 70 and 160 ${\rm \mu m}$ or SPIRE - 250, 350, 500 ${\rm \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.

Observational measures of the primordial helium mass fraction, $Y_{\rm P}$, are of interest for cosmology and fundamental particle physics. Current measures obtained from H II regions agree with the Standard Model prediction to approximately 1% precision, although these determinations may be affected by systematic uncertainties. This possibility can only be tested by independently measuring the helium abundance in new ways. Here, we propose a novel method to obtain a measurement of $Y_{\rm P}$ using hydrogen and helium recombination line emission from RELHICs: pristine, gas-rich but star-free low-mass dark matter haloes whose existence is predicted by hydrodynamical simulations. Although expected to be uncommon and intrinsically faint in emission, the primordial composition and simple physical properties of these objects make them an ideal laboratory to determine $Y_{\rm P}$. We present radiative transfer simulations to demonstrate the effectiveness of this approach, finding that comparing the emission in H and He lines, either via their volumetric emissivities, or integrated properties such as the surface brightness and total flux, may be used to infer $Y_{\rm P}$. Furthermore, we show that RELHICs can be used to provide an entirely novel constraint on the spectral slope of the ultraviolet background, and discuss the possibility of measuring this slope and the primordial helium abundance simultaneously.

The distribution of the absorption line broadening observed in the Ly$\alpha$ forest carries information about the temperature and widths of the filaments in the intergalactic medium (IGM). In this work, we present and test a new method for inferring the temperature of the IGM, the size of the absorbing filaments and the photo-ionization rate of hydrogen from the properties of absorption line broadening in the Ly$\alpha$ forest. We apply our method to mock spectra from the reference model of the EAGLE cosmological simulation, and we demonstrate that we are able to reconstruct the IGM properties. Our method explicitly takes into account the degeneracy between IGM temperature, the photo-ionization rate of hydrogen and the widths of the intergalactic filaments.

We introduce a simple analytic model of galaxy formation that links the growth of dark matter haloes in a cosmological background to the build-up of stellar mass within them. The model aims to identify the physical processes that drive the galaxy-halo co-evolution through cosmic time. The model restricts the role of baryonic astrophysics to setting the relation between galaxies and their haloes. Using this approach, galaxy properties can be directly predicted from the growth of their host dark matter haloes. We explore models in which the effective star formation efficiency within haloes is a function of mass (or virial temperature) and independent of time. Despite its simplicity, the model reproduces self-consistently the shape and evolution of the cosmic star formation rate density, the specific star formation rate of galaxies, and the galaxy stellar mass function, both at the present time and at high redshifts. By systematically varying the effective star formation efficiency in the model, we explore the emergence of the characteristic shape of the galaxy stellar mass function. The origin of the observed double Schechter function at low redshifts is naturally explained by two efficiency regimes in the stellar to halo mass relation, namely, a stellar feedback regulated stage, and a supermassive black hole regulated stage. By providing a set of analytic differential equations, the model can be easily extended and inverted, allowing the roles and impact of astrophysics and cosmology to be explored and understood.

We present the $I\kappa\epsilon\alpha$ model of galaxy formation, in which a galaxy's star formation rate is set by the balance between energy injected by feedback from massive stars and energy lost by the deepening of the potential of its host dark matter halo due to cosmological accretion. Such a balance is secularly stable provided that the star formation rate increases with the pressure in the star forming gas. The $I\kappa\epsilon\alpha$ model has four parameters that together control the feedback from star formation and the cosmological accretion rate onto a halo. $I\kappa\epsilon\alpha$ reproduces accurately the star formation rate as a function of halo mass and redshift in the EAGLE hydrodynamical simulation, even when all four parameters are held constant. It predicts the emergence of a star forming main sequence along which the specific star formation rate depends weakly on stellar mass with an amplitude that increases rapidly with redshift. We briefly discuss the emerging mass-metallicity relation, the evolution of the galaxy stellar mass function, and an extension of the model that includes feedback from active galactic nuclei (AGN). These self-regulation results are independent of the star formation law and the galaxy's gas content. Instead, star forming galaxies are shaped by the balance between stellar feedback and cosmological accretion, with accurately accounting for energy losses associated with feedback a crucial ingredient.

We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight, including two bright z ~ 3.2 quasars with close separation (~ 500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, log(Mh/Msun) ~ 11 - 13.5. For four of the groups, we find associated MgII absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ~ 480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionised gas seen in emission in lower-redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross section of cool gas and leading to the high fraction of strong MgII absorbers that we detect.

The observed Lyman-$\alpha$ flux power spectrum (FPS) is suppressed on scales below $\sim~ 30~{\rm km~s}^{-1}$. This cutoff could be due to the high temperature, $T_0$, and pressure, $p_0$, of the absorbing gas or, alternatively, it could reflect the free streaming of dark matter particles in the early universe. We perform a set of very high resolution cosmological hydrodynamic simulations in which we vary $T_0$, $p_0$ and the amplitude of the dark matter free streaming, and compare the FPS of mock spectra to the data. We show that the location of the dark matter free-streaming cutoff scales differently with redshift than the cutoff produced by thermal effects and is more pronounced at higher redshift. We, therefore, focus on a comparison to the observed FPS at $z>5$. We demonstrate that the FPS cutoff can be fit assuming cold dark matter, but it can be equally well fit assuming that the dark matter consists of $\sim 7$ keV sterile neutrinos in which case the cutoff is due primarily to the dark matter free streaming.

We present adaptive optics assisted integral field spectroscopy of 34 star-forming galaxies at $z$ = 0.8-3.3 selected from the HiZELS narrow-band survey. We measure the kinematics of the ionised interstellar medium on $\sim$1 kpc scales, and show that the galaxies are turbulent, with a median ratio of rotational to dispersion support of $v$/$\sigma$=0.82$\pm$0.13. We combine the dynamics with high-resolution rest-frame optical imaging and extract emission line rotation curves. We show that high-redshift star-forming galaxies follow a similar power-law trend in specific angular momentum with stellar mass as that of local late type galaxies. We exploit the high resolution of our data and examine the radial distribution of angular momentum within each galaxy by constructing total angular momentum profiles. Although the stellar mass of a typical star-forming galaxy is expected to grow by a factor $\sim$8 in the $\sim$5 Gyrs between $z$$\sim$3.3 and $z$$\sim$0.8, we show that the internal distribution of angular momentum becomes less centrally concentrated in this period i.e the angular momentum grows outwards. To interpret our observations, we exploit the EAGLE simulation and trace the angular momentum evolution of star forming galaxies from $z$$\sim$3 to $z$$\sim$0, identifying a similar trend of decreasing angular momentum concentration. This change is attributed to a combination of gas accretion in the outer disk, and feedback that preferentially arises from the central regions of the galaxy. We discuss how the combination of the growing bulge and angular momentum stabilises the disk and gives rise to the Hubble sequence.

The scatter (${\rm\sigma_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm \sigma_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in $ {\rm z \sim 0-4}$. We find that the relation evolves with time, with the dispersion depending on both stellar mass and redshift. The models point to an evolving U-shape form for the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation with the scatter being minimal at a characteristic mass $M^{\star}$ of ${\rm 10^{9.5}}$ ${\rm M_{\odot}}$ and increasing both at lower and higher masses. This implication is that the diversity of SFHs increases towards both at the low- and high-mass ends. We find that active galactic nuclei feedback is important for increasing the ${\rm \sigma_{\text{sSFR}}}$ for high mass objects. On the other hand, we suggest that SNe feedback increases the ${\rm \sigma_{\text{sSFR}}}$ of galaxies at the low-mass end. We also find that excluding galaxies that have experienced recent mergers does not significantly affect the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. Furthermore, we employ the combination of the EAGLE simulations with the radiative transfer code SKIRT to evaluate the effect of SFR/stellar mass diagnostics in the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation and find that the ${\rm SFR/M_{\star}}$ methodologies (e.g. SED fitting, UV+IR, UV+IRX-$\beta$) widely used in the literature to obtain intrinsic properties of galaxies have a large effect on the derived shape and normalization of the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation.

Photoheating of the gas in low-mass dark matter (DM) haloes prevents baryons from cooling, leaving the haloes free of stars. Gas in these 'dark' haloes remains exposed to the ultraviolet background (UVB), and so is expected to emit via fluorescent recombination lines. We present a set of radiative transfer simulations, which model dark haloes as spherical gas clouds in hydrostatic equilibrium with a DM halo potential, and in thermal equilibrium with the UVB at redshift z = 0. We use these simulations to predict surface brightnesses in H-alpha, which we show to have a characteristic ring-shaped morphology for haloes in a narrow mass range between 10^9.5 and 10^9.6 M_sun. We explore how this emission depends on physical parameters such as the DM density profile and the UVB spectrum. We predict the abundance of fluorescent haloes on the sky, and discuss possible strategies for their detection. We demonstrate how detailed observations of fluorescent rings can be used to infer the properties of the haloes which host them, such as their density profiles and the mass-concentration relation, as well as to directly measure the UVB amplitude.

We present results from a survey of galaxies in the fields of six z>3 Damped Lyman alpha systems (DLAs) using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). We report a high detection rate of up to ~80% of galaxies within 1000 km/s from DLAs and with impact parameters between 25 and 280 kpc. In particular, we discovered 5 high-confidence Lyman alpha emitters associated with three DLAs, plus up to 9 additional detections across five of the six fields. The majority of the detections are at relatively large impact parameters (>50 kpc) with two detections being plausible host galaxies. Among our detections, we report four galaxies associated with the most metal-poor DLA in our sample (Z/Z_sun = -2.33), which trace an overdense structure resembling a filament. By comparing our detections with predictions from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological simulations and a semi-analytic model designed to reproduce the observed bias of DLAs at z>2, we conclude that our observations are consistent with a scenario in which a significant fraction of DLAs trace the neutral regions within halos with a characteristic mass of 10^11-10^12 M_sun, in agreement with the inference made from the large-scale clustering of DLAs. We finally show how larger surveys targeting ~25 absorbers have the potential of constraining the characteristic masses of halos hosting high-redshift DLAs with sufficient accuracy to discriminate between different models.

The nearby Perseus galaxy cluster is a key target for indirect detection searches for decaying dark matter. We use the C-EAGLE simulations of galaxy clusters to predict the flux, width and shape of a dark matter decay line, paying particular attention to the unexplained 3.55keV line detected in the spectra of some galaxies and clusters, and the upcoming XRISM X-ray observatory mission. We show that the line width in C-EAGLE clusters similar to Perseus is typically [600-800]$\mathrm{kms^{-1}}$, and therefore narrower than the amplitude of the velocity dispersion of galaxies in the cluster. Haloes that are significantly disturbed can, however, exhibit galaxy velocity dispersions higher than $1000\mathrm{kms^{-1}}$, and in this case will show a large difference between the line profiles of on- and off-center observations. We show that the line profile is likely to be slightly asymmetric, but still well approximated by a Gaussian at the 10% level, and that the halo asymmetry can lead to fluxes that vary by a factor of two. In summary, we predict that, if the previously reported 3.55keV line detections do originate from dark matter decay, the XRISM mission will detect a line with a roughly Gaussian profile at a rest frame energy of 3.55keV, with a width $>600\mathrm{kms^{-1}}$ and flux approximately in the range $[4-9]\times10^{-8}\mathrm{counts/sec/cm^{2}}$.

We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations.

We investigate the connection between the morphology and internal kinematics of the stellar component of central galaxies with mass $M_\star > {10}^{9.5} {\rm M}_\odot$ in the EAGLE simulations. We compare several kinematic diagnostics commonly used to describe simulated galaxies, and find good consistency between them. We model the structure of galaxies as ellipsoids and quantify their morphology via the ratios of their principal axes. We show that the differentiation of blue star-forming and red quiescent galaxies using morphological diagnostics can be achieved with similar efficacy to the use of kinematical diagnostics, but only if one is able to measure both the flattening and the triaxiality of the galaxy. Flattened oblate galaxies exhibit greater rotational support than their spheroidal counterparts, but there is significant scatter in the relationship between morphological and kinematical diagnostics, such that kinematically-similar galaxies can exhibit a broad range of morphologies. The scatter in the relationship between the flattening and the ratio of the rotation and dispersion velocities ($v/\sigma$) correlates strongly with the anisotropy of the stellar velocity dispersion: at fixed $v/\sigma$, flatter galaxies exhibit greater dispersion in the plane defined by the intermediate and major axes than along the minor axis, indicating that the morphology of simulated galaxies is influenced significantly by the structure of their velocity dispersion. The simulations reveal that this anisotropy correlates with the intrinsic morphology of the galaxy's inner dark matter halo, i.e. the halo's morphology that emerges in the absence of dissipative baryonic physics. This implies the existence of a causal relationship between the morphologies of galaxies and that of their host dark matter haloes.

We analyse maps of the spatially-resolved nebular emission of $\approx$1500 star-forming galaxies at $z\approx0.6$-$2.2$ from deep KMOS and MUSE observations to measure the average shape of their rotation curves. We use these to test claims for declining rotation curves at large radii in galaxies at $z\approx1$-$2$ that have been interpreted as evidence for an absence of dark matter. We show that the shape of the average rotation curves, and the extent to which they decline beyond their peak velocities, depends upon the normalisation prescription used to construct the average curve. Normalising in size by the galaxy stellar disk-scale length after accounting for seeing effects ($R_{\rm{d}}^{\prime}$), we construct stacked position-velocity diagrams that trace the average galaxy rotation curve out to $6R_{\rm{d}}^{\prime}$ ($\approx$13 kpc, on average). Combining these curves with average HI rotation curves for local systems, we investigate how the shapes of galaxy rotation curves evolve over $\approx$10 Gyr. The average rotation curve for galaxies binned in stellar mass, stellar surface mass density and/or redshift is approximately flat, or continues to rise, out to at least $6R_{\rm{d}}^{\prime}$. We find a trend between the outer slopes of galaxies' rotation curves and their stellar mass surface densities, with the higher surface density systems exhibiting flatter rotation curves. Drawing comparisons with hydrodynamical simulations, we show that the average shapes of the rotation curves for our sample of massive, star-forming galaxies at $z\approx0$-$2.2$ are consistent with those expected from $\Lambda$CDM theory and imply dark matter fractions within $6R_{\rm{d}}$ of at least $\approx60$ percent.

We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes.

The cosmic spectral energy distribution (CSED) is the total emissivity as a function of wavelength of galaxies in a given cosmic volume. We compare the observed CSED from the UV to the submm to that computed from the EAGLE cosmological hydrodynamical simulation, post-processed with stellar population synthesis models and including dust radiative transfer using the SKIRT code. The agreement with the data is better than 0.15 dex over the entire wavelength range at redshift $z=0$, except at UV wavelengths where the EAGLE model overestimates the observed CSED by up to a factor 2. Global properties of the CSED as inferred from CIGALE fits, such as the stellar mass density, mean star formation density, and mean dust-to-stellar-mass ratio, agree to within better than 20 per cent. At higher redshift, EAGLE increasingly underestimates the CSED at optical-NIR wavelengths with the FIR/submm emissivity underestimated by more than a factor of 5 by redshift $z=1$. We believe that these differences are due to a combination of incompleteness of the EAGLE-SKIRT database, the small simulation volume and the consequent lack of luminous galaxies, and our lack of knowledge on the evolution of the characteristics of the interstellar dust in galaxies. The impressive agreement between the simulated and observed CSED at lower $z$ confirms that the combination of EAGLE and SKIRT dust processing yields a fairly realistic representation of the local Universe.

We exploit EAGLE, a cosmological hydrodynamical simulation, to reproduce the selection of the observed sub-millimeter (submm) galaxy population by selecting the model galaxies at $z \geq 1$ with mock submm fluxes $S_{850} \geq 1$ mJy. There is a reasonable agreement between the galaxies within this sample and the properties of the observed submm population, such as their star formation rates (SFRs) at $z<3$, redshift distribution and many integrated galaxy properties. We find that the bulk of the $S_{850} \geq 1$ mJy model population is at $z = 2.5$, and that they are massive galaxies ($M_* \sim 10^{11}$ Msol) with high dust masses ($M_{\mathrm{dust}} \sim 10^{8}$ Msol), gas fractions ($f_{\mathrm{gas}} \approx 50$%) and SFRs ($\dot M_* \approx 100$ Msol/yr). They have major and minor merger fractions similar to the general population, suggesting that mergers are not the primary driver of the model submm galaxies. Instead, the $S_{850} \geq 1$ mJy model galaxies yield high SFRs primarily because they maintain a significant gas reservoir as a result of hosting an undermassive black hole. In addition, we find that not all highly star-forming EAGLE galaxies have submm fluxes $S_{850} > 1$ mJy. Thus, we investigate the nature of $z \geq 1$ highly star-forming Submm-Faint galaxies (i.e., $\dot M_* \geq 80$ Msol/yr but $S_{850}< 1$ mJy). We find they are similar to the model submm galaxies; being gas rich and hosting undermassive black holes, however they are typically lower mass ($M_* \sim 10^{10}$ Msol) and are at higher redshifts ($z>4$). These typically higher-$z$ galaxies show stronger evidence for having been triggered by major mergers, and critically, they are likely missed by current submm surveys due to their higher dust temperatures. This suggests a potentially even larger contribution to the SFR density at $z > 3$ from dust-obscured systems than implied by current observations.

Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations, designed to enable the study of high-mass galaxy formation and evolution in cluster environments. We find 69 $M_{\rm BH}(M_{\star})$ outliers at $z=0$, defined as those with $ \rm M_{BH} >10^{7} M_{\odot}$ and $\rm M_{BH}/\rm M_{\star}> 0.01$. This paper focusses on a sample of 5 extreme outliers, that have been selected based on their $\rm M_{BH}$ and $\rm M_{\star}$ values, which are comparable to the most recent estimates of observed positive outliers. This sample of 5 outliers, classified as `Black hole monster galaxies' (BMGs), was traced back in time to study their origin and evolution. In agreement with the results of previous simulations for lower-mass $\rm M_{BH}(\rm M_{\star})$ outliers, we find that these galaxies became outliers due to a combination of their early formation times and tidal stripping. For BMGs with $\rm M_{BH} > 10^9 M_{\odot}$, major mergers (with a stellar mass ratio of $\mu > 0.25$) at early times ($z>2$) precede the rapid growth of their supermassive BHs. Furthermore, the scatter in the relation between $\rm M_{BH}$ and stellar velocity dispersion, $\sigma$, correlates positively with the scatter in [Mg/Fe]($\sigma$). This indicates that the alpha enhancement of these galaxies, which is closely related to their star formation history, is related to the growth of their central BHs.

The EAGLE simulation suite has previously been used to investigate the relationship between the stellar mass of galaxies, $M_{*}$, and the properties of dark matter halos, using the hydrodynamical reference simulation combined with a dark matter only (DMO) simulation having identical initial conditions. The stellar masses of central galaxies in halos with $M_{\mathrm{200c}} > 10^{11} \mathrm{M_{\odot}}$ were shown to correlate with the DMO halo maximum circular velocity, with $\approx 0.2$ dex of scatter that is uncorrelated with other DMO halo properties. Here we revisit the origin of the scatter in the $M_{*}-V_{\mathrm{max, DMO}}$ relation in EAGLE at $z = 0.1$. We find that the scatter in $M_{*}$ correlates with the mean age of the galaxy stellar population such that more massive galaxies at fixed $V_{\mathrm{max, DMO}}$ are younger. The scatter in the stellar mass and mean stellar population age results from variation in the baryonic mass, $M_{\mathrm{bary}} = M_{\mathrm{gas}} + M_{*}$, of the galaxies' progenitors at fixed halo mass and concentration. At the redshift of peak correlation ($z \approx 1$), the progenitor baryonic mass accounts for $75\%$ of the variance in the $z=0.1$ $M_{*}-V_{\mathrm{max, DMO}}$ relation. The scatter in the baryonic mass, in turn, is primarily set by differences in feedback strength and gas accretion over the course of the evolution of each halo.

Carbon enhanced metal poor (CEMP) stars with low abundances of neutron capture elements (CEMP-no stars) are ubiquitous among metal poor stars in the Milky Way. Recent observations have uncovered their two subgroups that differ in the carbon to magnesium ([C/Mg]) abundance ratio. Here we demonstrate that similar abundance patterns are also present in Milky Way-like galaxies in the EAGLE cosmological hydrodynamical simulation, where these patterns originate from the fact that stars may form from gas enriched predominantly by asymptotic giant branch (AGB) stars or by type-II supernovae. This occurs when stars form in the poorly mixed interstellar medium of Milky Way progenitor galaxies that are undergoing bursty star formation. The CEMP-no stars with lower [C/Mg] form at the onset of a starburst from gas enriched by low metallicity type-II supernovae that power a strong outflow, quenching further star formation. When star formation resumes following cosmological gas accretion, the CEMP-no stars with higher [C/Mg] form, with enrichment by AGB ejecta evident in their higher abundance of barium and lower abundance of magnesium. This suggests that bursty star formation in the progenitors of the Galaxy leaves a permanent imprint in the abundance patterns of CEMP stars.

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter.

We investigate the rapid growth phase of supermassive black holes (BHs) within the hydrodynamical cosmological \eagle simulation. This non-linear phase of BH growth occurs within $\sim$$L_{*}$ galaxies, embedded between two regulatory states of the galaxy host: in sub $L_{*}$ galaxies efficient stellar feedback regulates the gas inflow onto the galaxy and significantly reduces the growth of the central BH, while in galaxies more massive than $L_{*}$ efficient AGN feedback regulates the gas inflow onto the galaxy and curbs further non-linear BH growth. We find evolving critical galaxy and halo mass scales at which rapid BH growth begins. Galaxies in the low-redshift Universe transition into the rapid BH growth phase in haloes that are approximately an order of magnitude more massive than their high-redshift counterparts (\M{200} $\approx 10^{12.4}$~\Msol at $z \approx 0$ decreasing to \M{200} $\approx 10^{11.2}$~\Msol at $z \approx 6$). Instead, BHs enter the rapid growth phase at a fixed critical halo virial temperature ($T_{\mathrm{vir}} \approx 10^{5.6}$~K). We additionally show that major galaxy--galaxy interactions ($\mu \geq \frac{1}{4}$, where $\mu$ is the stellar mass ratio) play a substantial role in triggering the rapid growth phase of BHs in the low-redshift Universe, whilst potentially having a lower influence at high redshift. Approximately 40\% of BHs that initiate the rapid BH growth phase at $z \approx 0$ do so within $\pm 0.5$ dynamical times of a major galaxy--galaxy merger, a fourfold increase above what is expected from the background merger rate. We find that minor mergers ($\frac{1}{10} \leq \mu < \frac{1}{4}$) have a substantially lower influence in triggering the rapid growth phase at all epochs.

In this paper we introduce the SEAGLE (i.e. Simulating EAGLE LEnses) program, that approaches the study of galaxy formation through strong gravitational lensing, using a suite of high-resolution hydrodynamic simulations, Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. We introduce the simulation and analysis pipeline and present the first set of results from our analysis of early-type galaxies. We identify and extract an ensemble of simulated lens galaxies and use the GLAMER ray-tracing lensing code to create mock lenses similar to those observed in the SLACS and SL2S surveys, using a range of source parameters and galaxy orientations, including observational effects such as the Point-Spread-Function (PSF), pixelization and noise levels, representative of single-orbit observations with the Hubble Space Telescope (HST) using the ACS-F814W filter. We subsequently model these mock lenses using the code LENSED, treating them in the same way as observed lenses. We also estimate the mass model parameters directly from the projected surface mass density of the simulated galaxy, using an identical mass model family. We perform a three-way comparison of all the measured quantities with real lenses. We find the average total density slope of EAGLE lenses, $t=2.26\; (0.25\; \rm{rms})$ to be higher than SL2S, $t=2.16$ or SLACS, $t=2.08$. We find a very strong correlation between the external shear ($\gamma$) and the complex ellipticity ($\epsilon$), with $\gamma \sim \epsilon/4$. This correlation indicates a degeneracy in the lens mass modeling. We also see a dispersion between lens modeling and direct fitting results, indicating systematical biases.

We investigate the population of dwarf galaxies with stellar masses similar to the Large Magellanic Cloud (LMC) and M33 in the EAGLE galaxy formation simulation. In the field, galaxies reside in haloes with stellar-to-halo mass ratios of $1.03^{+0.50}_{-0.31}\times10^{-2}$ (68% confidence level); systems like the LMC, which have an SMC-mass satellite, reside in haloes about 1.3 times more massive, which suggests an LMC halo mass at infall, $M_{200}=3.4^{+1.8}_{-1.2}\times10^{11}M_\odot$ (68% confidence level). The colour distribution of dwarfs is bimodal, with the red galaxies ($g-r>0.6$) being mostly satellites. The fraction of red LMC-mass dwarfs is 15% for centrals, and for satellites this fraction increases rapidly with host mass: from 10% for satellites of Milky Way (MW)-mass haloes to nearly 90% for satellites of groups and clusters. The quenching timescale, defined as the time after infall when half of the satellites have acquired red colours, decreases with host mass from ${>}5$ Gyrs for MW-mass hosts to $2.5$ Gyrs for cluster mass hosts. The satellites of MW-mass haloes have higher star formation rates and bluer colours than field galaxies. This is due to enhanced star formation triggered by gas compression shortly after accretion. Both the LMC and M33 have enhanced recent star formation that could be a manifestation of this process. After infall into their MW-mass hosts, the $g-r$ colours of LMC-mass dwarfs become bluer for the first 2 Gyrs, after which they rapidly redden. LMC-mass dwarfs fell into their MW-mass hosts only relatively recently, with more than half having an infall time of less than 3.5 Gyrs.

We use the EAGLE simulations to study the oxygen abundance gradients of gas discs in galaxies within the stellar mass range [10^9.5, 10^10.8]Mo at z=0. The estimated median oxygen gradient is -0.011 (0.002) dex kpc^-1, which is shallower than observed. No clear trend between simulated disc oxygen gradient and galaxy stellar mass is found when all galaxies are considered. However, the oxygen gradient shows a clear correlation with gas disc size so that shallower abundance slopes are found for increasing gas disc sizes. Positive oxygen gradients are detected for ~40 per cent of the analysed gas discs, with a slight higher frequency in low mass galaxies. Galaxies that have quiet merger histories show a positive correlation between oxygen gradient and stellar mass, so that more massive galaxies tend to have shallower metallicity gradients. At high stellar mass, there is a larger fraction of rotational-dominated galaxies in low density regions. At low stellar mass, non-merger galaxies show a large variety of oxygen gradients and morphologies. The normalization of the disc oxygen gradients in non-merger galaxies by the effective radius removes the trend with stellar mass. Conversely, galaxies that experienced mergers show a weak relation between oxygen gradient and stellar mass. Additionally, the analysed EAGLE discs show no clear dependence of the oxygen gradients on local environment, in agreement with current observational findings.

We present the first results from our survey of intervening and proximate Lyman limit systems (LLSs) at $z$$\sim$2.0-2.5 using the Wide Field Camera 3 on-board the Hubble Space Telescope. The quasars in our sample are projected pairs with proper transverse separations $R_\perp$$\leq$150 kpc and line of sight velocity separations $\lesssim$11,000 km/s. We construct a stacked ultraviolet (rest-frame wavelengths 700-2000\AA) spectrum of pairs corrected for the intervening Lyman forest and Lyman continuum absorption. The observed spectral composite presents a moderate flux excess for the most prominent broad emission lines, a $\sim$30% decrease in flux at $\lambda$=800-900\AA\ compared to a stack of brighter quasars not in pairs at similar redshifts, and lower values of the mean free path of the HI ionizing radiation for pairs ($\lambda_{\rm mfp}^{912}=140.7\pm20.2~h_{70}^{-1}$Mpc) compared to single quasars ($\lambda_{\rm mfp}^{912}=213.8\pm28~h_{70}^{-1}$Mpc) at the average redshift $z\simeq2.44$. From the modelling of LLS absorption in these pairs, we find a higher ($\sim$20%) incidence of proximate LLSs with $\log N_{\rm HI}\geq17.2$ at $\delta v$$<$5,000 km/s compared to single quasars ($\sim$6%). These two rates are different at the 5$\sigma$ level. Moreover, we find that optically-thick absorbers are equally shared between foreground and background quasars. Based on these pieces of evidence, we conclude that there is a moderate excess of gas absorbing Lyman continuum photons in our closely-projected quasar pairs compared to single quasars. We argue that this gas arises mostly within large-scale structures or partially neutral regions inside the dark matter haloes where these close pairs reside.

The Hubble sequence provides a useful classification of galaxy morphology at low redshift. However, morphologies are not static, but rather evolve as the growth of structure proceeds through mergers, accretion and secular processes. We investigate how kinematically defined disc and spheroidal structures form and evolve in the EAGLE hydrodynamic simulation of galaxy formation. At high redshift most galaxies of all masses are asymmetric. By redshift $z\simeq 1.5$ the Hubble sequence is established and after this time most of the stellar mass is in spheroids whose contribution to the stellar mass budget continues to rise to the present day. The stellar mass fraction in discs peaks at $z\simeq 0.5$ but overall remains subdominant at all times although discs contribute most of the stellar mass in systems of mass $M^*\sim 10^{10.5}{\rm M_\odot}$ at $z \leq 1.5$. Star formation occurs predominantly in disc structures throughout most of cosmic time but morphological transformations rearrange stars, thus establishing the low-redshift morphological mix. Morphological transformations are common and we quantify the rates at which they occur. The rate of growth of spheroids decreases at $z < 2$ while the rate of decay of discs remains roughly constant at $z < 1$. Finally, we find that the prograde component of galaxies becomes increasingly dynamically cold with time.

We study the correlation between stellar mass and stellar metallicity in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological hydrodynamical simulations. At a given stellar mass, simulated galaxies with lower stellar metallicities show, on average, higher gas fractions, higher specific star formation rates and younger stellar populations. Active galactic nuclei feedback seems to play an important role on the determination of the stellar metallicity at high stellar masses. In general, simulated systems follow a well-defined anticorrelation between stellar metallicity and gas fraction, which does not evolve significantly with redshift. All these trends are consistent with previous findings regarding the metallicity of the star-forming gas in EAGLE.

We investigate the effect of the accelerated expansion of the Universe due to a cosmological constant, $\Lambda$, on the cosmic star formation rate. We utilise hydrodynamical simulations from the EAGLE suite, comparing a $\Lambda$CDM Universe to an Einstein-de Sitter model with $\Lambda=0$. Despite the differences in the rate of growth of structure, we find that dark energy, at its observed value, has negligible impact on star formation in the Universe. We study these effects beyond the present day by allowing the simulations to run forward into the future ($t>13.8$ Gyr). We show that the impact of $\Lambda$ becomes significant only when the Universe has already produced most of its stellar mass, only decreasing the total co-moving density of stars ever formed by ${\approx}15\%$. We develop a simple analytic model for the cosmic star formation rate that captures the suppression due to a cosmological constant. The main reason for the similarity between the models is that feedback from accreting black holes dramatically reduces the cosmic star formation at late times. Interestingly, simulations without feedback from accreting black holes predict an upturn in the cosmic star formation rate for $t>15$ Gyr due to the rejuvenation of massive ($ > 10^{11} \mathrm{M}_{\odot}$) galaxies. We briefly discuss the implication of the weak dependence of the cosmic star formation on $\Lambda$ in the context of the anthropic principle.

The massive stars that ionised the Universe have short lifetimes and can only be studied near the time of formation, but any low mass stars that formed contemporaneously might be observable in the local Universe today. We study the abundance pattern and spatial distribution of these `siblings of reionizers' (SoRs) in the EAGLE cosmological hydrodynamical simulation. SoRs tend to be enriched to super-solar levels in $\alpha$-elements compared to iron. In {\sc eagle} galaxies resembling the Milky Way, $\sim 40$ percent of carbon-enhanced metal poor (CEMP) stars are SoRs. Conversely, $\sim 10$ percent of all SoRs are CEMP stars. This fraction increases to $\gtrsim 50$ percent for SoRs of metallicity [Fe/H]$<-4$, and at such low metallicities, most of the CEMP stars are of CEMP-no subtype that are lacking neutron capture elements. Although these numbers may well depend on the details of the physical models implemented in EAGLE, the trends we describe are robust as they result from the strong feedback from star formation in early galaxies, itself a key ingredient of most current models of galaxy formation. We further find that most SoRs today reside in halos with mass $M_h\gtrapprox 10^{12}$ M$_\odot$, and 50 percent of them are in the halo of their central galaxy (distance $>10$ kpc), mainly because they were accreted onto their current host rather than formed in-situ. To a good approximation, the SoRs are CEMP-no stars that reside in the stellar halos of massive galaxies, with nearly half of them contributing to the intracluster light in groups and clusters.

Models of the very early universe, including inflationary models, are argued to produce varying universe domains with different values of fundamental constants and cosmic parameters. Using the cosmological hydrodynamical simulation code from the eagle collaboration, we investigate the effect of the cosmological constant on the formation of galaxies and stars. We simulate universes with values of the cosmological constant ranging from Lambda = 0 to Lambda_0 = 300, where Lambda_0 is the value of the cosmological constant in our Universe. Because the global star formation rate in our Universe peaks at t = 3.5 Gyr, before the onset of accelerating expansion, increases in Lambda of even an order of magnitude have only a small effect on the star formation history and efficiency of the universe. We use our simulations to predict the observed value of the cosmological constant, given a measure of the multiverse. Whether the cosmological constant is successfully predicted depends crucially on the measure. The impact of the cosmological constant on the formation of structure in the universe does not seem to be a sharp enough function of Lambda to explain its observed value alone.

We investigate how the nature of the galaxies that reionized the Universe affects the duration of reionization. We contrast two sets of models: one in which galaxies on the faint side of the luminosity function dominate the ionizing emissivity, and a second in which the galaxies on the bright side of the luminosity function dominate. The faint-end of the luminosity function evolves slowly, therefore the transition from mostly neutral to mostly ionized state takes a much longer time in the first set of models compared to the second. Existing observational constraints on the duration of this transition are relatively weak, but taken at face value prefer the model in which galaxies on the bright side play a major role. Measurements of the kinetic Sunyaev Zeldovich effect in the cosmic microwave background from the epoch of reionization also point in the same direction.

We examine the origin of radial and vertical gradients in the age/metallicity of the stellar component of a galaxy disc formed in the APOSTLE cosmological hydrody- namical simulations. Some of these gradients resemble those in the Milky Way, where they have sometimes been interpreted as due to internal evolution, such as scattering off giant molecular clouds, radial migration driven by spiral patterns, or orbital reso- nances with a bar. Secular processes play a minor role in the simulated galaxy, which lacks strong spiral or bar patterns, and where such gradients arise as a result of the gradual enrichment of a gaseous disc that is born thick but thins as it turns into stars and settles into centrifugal equilibrium. The settling is controlled by the feedback of young stars; which links the star formation, enrichment, and equilibration timescales, inducing radial and vertical gradients in the gaseous disc and its descendent stars. The kinematics of coeval stars evolve little after birth and provide a faithful snapshot of the gaseous disc structure at the time of their formation. In this interpretation, the age-velocity dispersion relation would reflect the gradual thinning of the disc rather than the importance of secular orbit scattering; the outward flaring of stars would result from the gas disc flare rather than from radial migration; and vertical gradients would arise because the gas disc gradually thinned as it enriched. Such radial and vertical trends might just reflect the evolving properties of the parent gaseous disc, and are not necessarily the result of secular evolutionary processes.

Spectroscopic surveys of the Galaxy reveal that its disc stars exhibit a spread in $\mathrm{[\alpha/Fe]}$ at fixed $\mathrm{[Fe/H]}$, manifest at some locations as a bimodality. The origin of these diverse, and possibly distinct, stellar populations in the Galactic disc is not well understood. We examine the Fe and $\alpha$-element evolution of 133 Milky Way-like galaxies from the EAGLE simulation, to investigate the origin and diversity of their $\mathrm{[\alpha/Fe]}$-$\mathrm{[Fe/H]}$ distributions. We find that bimodal $\mathrm{[\alpha/Fe]}$ distributions arise in galaxies whose gas accretion histories exhibit episodes of significant infall at both early and late times, with the former fostering more intense star formation than the latter. The shorter characteristic consumption timescale of gas accreted in the earlier episode suppresses its enrichment with iron synthesised by Type Ia SNe, resulting in the formation of a high-$\mathrm{[\alpha/Fe]}$ sequence. We find that bimodality in $\mathrm{[\alpha/Fe]}$ similar to that seen in the Galaxy is rare, appearing in approximately 5 percent of galaxies in our sample. We posit that this is a consequence of an early gas accretion episode requiring the mass accretion history of a galaxy's dark matter halo to exhibit a phase of atypically-rapid growth at early epochs. The scarcity of EAGLE galaxies exhibiting distinct sequences in the $\mathrm{[\alpha/Fe]}$-$\mathrm{[Fe/H]}$ plane may therefore indicate that the Milky Way's elemental abundance patterns, and its accretion history, are not representative of the broader population of $\sim L^\star$ disc galaxies.

We present dust-attenuated and dust emission fluxes for sufficiently resolved galaxies in the EAGLE suite of cosmological hydrodynamical simulations, calculated with the SKIRT radiative transfer code. The post-processing procedure includes specific components for star formation regions, stellar sources, and diffuse dust, and takes into account stochastic heating of dust grains to obtain realistic broad-band fluxes in the wavelength range from ultraviolet to sub-millimeter. The mock survey includes nearly half a million simulated galaxies with stellar masses above 10^8.5 solar masses across six EAGLE models. About two thirds of these galaxies, residing in 23 redshift bins up to z=6, have a sufficiently resolved metallic gas distribution to derive meaningful dust attenuation and emission, with the important caveat that the same dust properties were used at all redshifts. These newly released data complement the already publicly available information about the EAGLE galaxies, which includes intrinsic properties derived by aggregating the properties of the smoothed particles representing matter in the simulation. We further provide an open source framework of Python procedures for post-processing simulated galaxies with the radiative transfer code SKIRT. The framework allows any third party to calculate synthetic images, SEDs, and broadband fluxes for EAGLE galaxies, taking into account the effects of dust attenuation and emission.

We investigate the nature of carbon-enhanced metal poor (CEMP) stars in Milky Way (MW) analogues selected from the EAGLE cosmological hydrodynamical simulation. The stellar evolution model in EAGLE includes the physics of enrichment by asymptotic giant branch (AGB) stars, winds from massive stars, and type I and type II supernovae (SNe). In the simulation, star formation in young MW progenitors is bursty due to efficient stellar feedback, which causes poor metal mixing leading to the formation of CEMP stars with extreme abundance patterns. In this scenario, two classes of CEMP stars emerge: those mostly enriched by low-metallicity type II SNe with low Fe yields that drive galactic outflows, and those mostly enriched by AGB stars when a gas-poor progenitor accretes pristine gas. The first class resembles CEMP-no stars with high [C/Fe] and low [C/O], the second class resembles CEMP-s stars overabundant in s-process elements and high values of [C/O]. This scenario explains several trends seen in data: (i) the increase in the scatter and median of [C/O] at low and decreasing [O/H], (ii) the trend of stars with very low [Fe/H] or [C/H] to be of type CEMP-no, and (iii) the reduction in the scatter of [{\alpha}/Fe] with atomic number in metal poor stars. In this interpretation CEMP-no stars were enriched by the stars that enabled galaxies to reionise the Universe.

The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies are consistent with constant-density `cores' in their inner regions and rule out `cuspy' Navarro-Frenk-White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Lambda cold dark matter cosmological hydrodynamics simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 53 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed.

We present sensitive 870$\mu$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distributions in the infrared band (8 - 1000$\mu$m) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the EAGLE hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above $\approx 10^{10}$M$_{\odot}$.

We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range $14.0 \le \log_{10}(M_{\rm 200c}/\mathrm{M_\odot}) \le 15.4$, with their sophisticated sub-grid physics modelling and high numerical resolution (sub-kpc gravitational softening) making them ideal for this purpose. We find that selecting galaxies by their total mass results in a velocity dispersion that is 5-10 per cent higher than the dark matter particles. However, selecting galaxies by their stellar mass results in an almost unbiased ($<5$ per cent) estimator of the velocity dispersion. This result holds out to $z=1.5$ and is relatively insensitive to the choice of cluster aperture, varying by less than 5 per cent between $r_{\rm 500c}$ and $r_{\rm 200m}$. We show that the velocity bias is a function of the time spent by a galaxy inside the cluster environment. Selecting galaxies by their total mass results in a larger bias because a larger fraction of objects have only recently entered the cluster and these have a velocity bias above unity. Galaxies that entered more than $4 \, \mathrm{Gyr}$ ago become progressively colder with time, as expected from dynamical friction. We conclude that velocity bias should not be a major issue when estimating cluster masses from kinematic methods.

We use EAGLE to quantify the effect galaxy mergers have on the stellar specific angular momentum of galaxies, $j_{\rm stars}$. We split mergers into: dry (gas-poor)/wet (gas-rich), major/minor, and different spin alignments and orbital parameters. Wet (dry) mergers have an average neutral gas-to-stellar mass ratio of $1.1$ ($0.02$), while major (minor) mergers are those with stellar mass ratios $\ge 0.3$ ($0.1-0.3$). We correlate the positions of galaxies in the $j_{\rm stars}$-stellar mass plane at $z=0$ with their merger history, and find that galaxies of low spins suffered dry mergers, while galaxies of normal/high spins suffered predominantly wet mergers, if any. The radial $j_{\rm stars}$ profiles of galaxies that went through dry mergers are deficient by $\approx 0.3$~dex at $r\lesssim 10\,r_{50}$ (with $r_{50}$ being the half-stellar mass radius), compared to galaxies that went through wet mergers. Studying the merger remnants reveals that dry mergers reduce $j_{\rm stars}$ by $\approx 30$\%, while wet mergers increase it by $\approx 10$\%, on average. The latter is connected to the build-up of the bulge by newly formed stars of high rotational speed. Moving from minor to major mergers accentuates these effects. When the spin vectors of the galaxies prior to the dry merger are misaligned, $j_{\rm stars}$ decreases to a greater magnitude, while in wet mergers co-rotation and high orbital angular momentum efficiently spun-up galaxies. We predict what would be the observational signatures in the $j_{\rm stars}$ profiles driven by dry mergers: (i) shallow radial profiles and (ii) profiles that rise beyond $\approx 10\,r_{50}$, both of which are significantly different from spiral galaxies.

It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.

We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo gas provides a clear signature of the formation of a hot corona. Haloes of total mass $10^{11.5}-10^{12}{\rm{M}}_{\odot}$ develop a hot corona independent of redshift, at least in the interval $z=0-4$ where the simulation has sufficiently good statistics. We analyse the build up of the hot gas mass in the halo, $M_{\rm{hot}}$, as a function of halo mass and redshift and find that while more energetic galactic winds powered by SNe increases $M_{\rm{hot}}$, AGN feedback reduces it by ejecting gas from the halo. We also study the thermal properties of gas accreting onto haloes and measure the fraction of shock-heated gas as a function of redshift and halo mass. We develop analytic and semianalytic approaches to estimate a `critical halo mass', $M_{\rm{crit}}$, for hot halo formation. We find that the mass for which the heating rate produced by accretion shocks equals the radiative cooling rate, reproduces the mass above which haloes develop a significant hot atmosphere. This yields a mass estimate of $M_{\rm{crit}} \approx 10^{11.7}{\rm{M}}_{\odot}$ at $z=0$, which agrees with the simulation results. The value of $M_{\rm{crit}}$ depends more strongly on the cooling rate than on any of the feedback parameters.

This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K.

We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhalos compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or `simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low mass haloes is required to form galaxies than is predicted by the EAGLE simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models.

We quantify the correlations between gas-phase and stellar metallicities and global properties of galaxies, such as stellar mass, halo mass, age and gas fraction, in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological hydrodynamical simulations. The slope of the correlation between stellar mass and metallicity of star-forming (SF) gas ($M_*-Z_{\rm SF,gas}$ relation) depends somewhat on resolution, with the higher-resolution run reproducing a steeper slope. This simulation predicts a non-zero metallicity evolution, increasing by $\approx 0.5$ dex at $\sim 10^9 {\rm M}_{\odot}$ since $z = 3$. The simulated relation between stellar mass, metallicity and star formation rate at $z \lesssim 5$ agrees remarkably well with the observed fundamental metallicity relation. At $M_* \lesssim 10^{10.3} {\rm M}_{\odot}$ and fixed stellar mass, higher metallicities are associated with lower specific star formation rates, lower gas fractions and older stellar populations. On the other hand, at higher $M_*$, there is a hint of an inversion of the dependence of metallicity on these parameters. The fundamental parameter that best correlates with the metal content, in the simulations, is the gas fraction. The simulated gas fraction-metallicity relation exhibits small scatter and does not evolve significantly since $z = 3$. In order to better understand the origin of these correlations, we analyse a set of lower resolution simulations in which feedback parameters are varied. We find that the slope of the simulated $M_*-Z_{\rm SF,gas}$ relation is mostly determined by stellar feedback at low stellar masses ($M_* \lesssim 10^{10} {\rm M}_{\odot}$), and at high masses ($M_* \gtrsim 10^{10} {\rm M}_{\odot}$) by the feedback from active galactic nuclei.

We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8 $ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$\alpha$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations. We find that the constraints implied by different indicators are inconsistent with each other for the highest star-forming objects at z < 2, a problem that is possibly related to selection biases and the uncertainties of dust attenuation effects. EAGLE's feedback parameters were calibrated to reproduce realistic galaxy sizes and stellar masses at z = 0.1. In this work we test if and why those choices yield realistic Star Formation Rates (SFRs) for $z \sim 0-8$ as well. We demonstrate that SNe feedback plays a major role at setting the abundance of galaxies at all star-forming regimes, especially at high redshifts. On the contrary, Active Galactic Nuclei (AGN) feedback becomes more prominent at lower redshifts and is a major mechanism that affects only the highest star-forming systems. Furthermore, we find that galaxies with SFR $\sim 1-10 \, {\rm M_{\odot} \, yr^{-1}}$ dominate the CSFRD at redshifts z < 5, while rare high star-forming galaxies (SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$) contribute significantly only briefly around the peak era ($z \sim 2$) and then are quenched by AGN feedback. In the absence of this prescription objects with SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$ would dominate the CSFRD, while the cosmic budget of star formation would be extremely high. Finally, we demonstrate that the majority of the cosmic star formation occurs in relatively rare high mass halos ($ {\rm M_{Halo}} \sim 10^{11-13} \, {\rm M_{\odot}}$) even at the earliest epochs.

We present new observations acquired with the Multi Unit Spectroscopic Explorer instrument on the Very Large Telescope in a quasar field that hosts a high column-density damped Ly{\alpha} absorber (DLA) at z~3.25. We detect Ly{\alpha} emission from a nebula at the redshift of the DLA with line luminosity (27+/-1)x1e41 erg/s, which extends over 37+/-1 kpc above a surface brightness limit of 6x1e-19 erg/s/cm2/arcsec2 at a projected distance of 30.5+/-0.5 kpc from the quasar sightline. Two clumps lie inside this nebula, both with Ly{\alpha} rest-frame equivalent width > 50 A and with relative line-of-sight velocities aligned with two main absorption components seen in the DLA spectrum. In addition, we identify a compact galaxy at a projected distance of 19.1+/-0.5 kpc from the quasar sightline. The galaxy spectrum is noisy but consistent with that of a star-forming galaxy at the DLA redshift. We argue that the Ly{\alpha} nebula is ionized by radiation from star formation inside the two clumps, or by radiation from the compact galaxy. In either case, these data imply the presence of a structure with size >>50 kpc inside which galaxies are assembling, a picture consistent with galaxy formation in groups and filaments as predicted by cosmological simulations such as the EAGLE simulations.

We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bah\'e et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of $1.8\times10^{6}\,\mathrm{M}_{\odot}$ and physical softening length of $0.7\,\mathrm{kpc}$. In this paper, we introduce the sample and present the low-redshift global properties of the clusters. We calculate the X-ray properties in a manner consistent with observational techniques, demonstrating the bias and scatter introduced by using estimated masses. We find the total stellar content and black hole masses of the clusters to be in good agreement with the observed relations. However, the clusters are too gas rich, suggesting that the AGN feedback model is not efficient enough at expelling gas from the high-redshift progenitors of the clusters. The X-ray properties, such as the spectroscopic temperature and the soft-band luminosity, and the Sunyaev-Zel'dovich properties are in reasonable agreement with the observed relations. However, the clusters have too high central temperatures and larger-than-observed entropy cores, which is likely driven by the AGN feedback after the cluster core has formed. The total metal content and its distribution throughout the ICM are a good match to the observations.

We study HI and metal-line absorption around $z\approx2$ star-forming galaxies by comparing an analysis of data from the Keck Baryonic Structure Survey to mock spectra generated from the EAGLE cosmological, hydrodynamical simulations. We extract sightlines from the simulations and compare the properties of the absorption by HI, CIV and SiIV around simulated and observed galaxies using pixel optical depths. We mimic the resolution, pixel size, and signal-to-noise ratio of the observations, as well as the distributions of impact parameters and galaxy redshift errors. We find that the EAGLE reference model is in excellent agreement with the observations. In particular, the simulation reproduces the high metal-line optical depths found at small galactocentric distances, the optical depth enhancements out to impact parameters of 2 proper Mpc, and the prominent redshift-space distortions which we find are due to peculiar velocities rather than redshift errors. The agreement is best for halo masses $\sim10^{12.0}$ M$_\odot$, for which the observed and simulated stellar masses also agree most closely. We examine the median ion mass-weighted radial gas velocities around the galaxies, and find that most of the gas is infalling, with the infall velocity depending on halo rather than stellar mass. From this we conclude that the observed redshift-space distortions are predominantly caused by infall rather than outflows.

We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $\Lambda$CDM. This is because: (i) disk galaxies in $\Lambda$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $\Lambda$CDM galaxy formation model. Explaining the MDAR in $\Lambda$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process.

Cosmological parameter constraints from observations of time-delay lenses are becoming increasingly precise. However, there may be significant bias and scatter in these measurements due to, among other things, the so-called mass-sheet degeneracy. To estimate these uncertainties, we analyze strong lenses from the largest EAGLE hydrodynamical simulation. We apply a mass-sheet transformation to the radial density profiles of lenses, and by selecting lenses near isothermality, we find that the bias on H0 can be reduced to 5% with an intrinsic scatter of 10%, confirming previous results performed on a different simulation data set. We further investigate whether combining lensing observables with kinematic constraints helps to minimize this bias. We do not detect any significant dependence of the bias on lens model parameters or observational properties of the galaxy, but depending on the source--lens configuration, a bias may still exist. Cross lenses provide an accurate estimate of the Hubble constant, while fold (double) lenses tend to be biased low (high). With kinematic constraints, double lenses show bias and intrinsic scatter of 6% and 10%, respectively, while quad lenses show bias and intrinsic scatter of 0.5% and 10%, respectively. For lenses with a reduced $\chi^2 > 1$, a power-law dependence of the $\chi^2$ on the lens environment (number of nearby galaxies) is seen. Lastly, we model, in greater detail, the cases of two double lenses that are significantly biased. We are able to remove the bias, suggesting that the remaining biases could also be reduced by carefully taking into account additional sources of systematic uncertainty.

We study present-day galaxy clustering in the EAGLE cosmological hydrodynamical simulation. EAGLE's galaxy formation parameters were calibrated to reproduce the redshift $z=0.1$ galaxy stellar mass function, and the simulation also reproduces galaxy colours well. The simulation volume is too small to correctly sample large-scale fluctuations and we therefore concentrate on scales smaller than a few mega parsecs. We find very good agreement with observed clustering measurements from the Galaxy And Mass Assembly (GAMA) survey, when galaxies are binned by stellar mass, color, or luminosity. However, low-mass red-galaxies are clustered too strongly, which is at least partly due to limited numerical resolution. Apart from this limitation, we conclude that EAGLE galaxies inhabit similar dark matter haloes as observed GAMA galaxies, and that the radial distribution of satellite galaxies as function of stellar mass and colour is similar to that observed as well.

We present mock optical images, broad-band and H$\alpha$ fluxes, and D4000 spectral indices for 30,145 galaxies from the EAGLE hydrodynamical simulation at redshift $z=0.1$, modelling dust with the SKIRT Monte Carlo radiative transfer code. The modelling includes a subgrid prescription for dusty star-forming regions, with both the subgrid obscuration of these regions and the fraction of metals in diffuse interstellar dust calibrated against far-infrared fluxes of local galaxies. The predicted optical colours as a function of stellar mass agree well with observation, with the SKIRT model showing marked improvement over a simple dust screen model. The orientation dependence of attenuation is weaker than observed because EAGLE galaxies are generally puffier than real galaxies, due to the pressure floor imposed on the interstellar medium. The mock H$\alpha$ luminosity function agrees reasonably well with the data, and we quantify the extent to which dust obscuration affects observed H$\alpha$ fluxes. The distribution of D4000 break values is bimodal, as observed. In the simulation, 20$\%$ of galaxies deemed `passive' for the SKIRT model, i.e. exhibiting D4000 $> 1.8$, are classified `active' when ISM dust attenuation is not included. The fraction of galaxies with stellar mass greater than $10^{10}$ M$_\odot$ that are deemed passive is slightly smaller than observed, which is due to low levels of residual star formation in these simulated galaxies. Colour images, fluxes and spectra of EAGLE galaxies are to be made available through the public EAGLE database.

The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make use of the APOSTLE cosmological hydrodynamical simulations of the Local Group to test the validity and accuracy of such mass estimators for both dispersion and rotation-supported galaxies, for field and satellite galaxies, and for galaxies of varying masses, shapes, and velocity dispersion anisotropies. We find that the mass estimators of Walker et al. and Wolf et al. are able to recover the masses of dispersion-dominated systems with little systematic bias, but with a 1-sigma scatter of 25 and 23 percent, respectively. The error on the estimated mass is dominated by the impact of the 3D shape of the stellar mass distribution, which is difficult to constrain observationally. This intrinsic scatter becomes the dominant source of uncertainty in the masses estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky Way, where the observational errors in their sizes and velocity dispersions are small. Such scatter may also affect the inner density slopes of dSphs derived from multiple stellar populations, relaxing the significance with which Navarro-Frenk-White profiles may be excluded, depending on the degree to which the relevant properties of the different stellar populations are correlated. Finally, we derive a new optimal mass estimator that removes the residual biases and achieves a statistically significant reduction in the scatter to 20 percent overall for dispersion-dominated galaxies, allowing more precise and accurate mass estimates.

We investigate the relation between kinematic morphology, intrinsic colour and stellar mass of galaxies in the EAGLE cosmological hydrodynamical simulation. We calculate the intrinsic u-r colours and measure the fraction of kinetic energy invested in ordered corotation of 3562 galaxies at z=0 with stellar masses larger than $10^{10}M_{\odot}$. We perform a visual inspection of gri-composite images and find that our kinematic morphology correlates strongly with visual morphology. EAGLE produces a galaxy population for which morphology is tightly correlated with the location in the colour- mass diagram, with the red sequence mostly populated by elliptical galaxies and the blue cloud by disc galaxies. Satellite galaxies are more likely to be on the red sequence than centrals, and for satellites the red sequence is morphologically more diverse. These results show that the connection between mass, intrinsic colour and morphology arises from galaxy formation models that reproduce the observed galaxy mass function and sizes.

We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different prescriptions for unresolved stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: different feedback implementations -- which produce different galaxy populations -- mainly shift galaxies along the relation, rather than perpendicular to it. Furthermore, galaxies exhibit a characteristic acceleration, $g_{\dagger}$, above which baryons dominate the mass budget, as observed. These observations, consistent with simple modified Newtonian dynamics, can be accommodated within the standard cold dark matter paradigm.

We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses Mstar > 10^10 M_Sun this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (Mstar < 10^10 M_Sun) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching. Conversely, the simulations predict an overabundance of massive galaxies in clusters compared to the field that persists to their far outskirts (> 5r_200c). This is caused by a significantly increased stellar mass fraction of (sub-)haloes in the cluster environment, by up to ~0.3 dex even well beyond r_200c. Haloes near clusters are also more concentrated than equally massive field haloes, but these two effects are largely uncorrelated.

We investigate the connection between the star formation rate (SFR) of galaxies and their central black hole accretion rate (BHAR) using the EAGLE cosmological hydrodynamical simulation. We find, in striking concurrence with recent observational studies, that the <SFR>--BHAR relation for an AGN selected sample produces a relatively flat trend, whilst the <BHAR>--SFR relation for a SFR selected sample yields an approximately linear trend. These trends remain consistent with their instantaneous equivalents even when both SFR and BHAR are time-averaged over a period of 100~Myr. There is no universal relationship between the two growth rates. Instead, SFR and BHAR evolve through distinct paths that depend strongly on the mass of the host dark matter halo. The galaxies hosted by haloes of mass M200 $\lesssim 10^{11.5}$Msol grow steadily, yet black holes (BHs) in these systems hardly grow, yielding a lack of correlation between SFR and BHAR. As haloes grow through the mass range $10^{11.5} \lesssim$ M200 $\lesssim 10^{12.5 }$Msol BHs undergo a rapid phase of non-linear growth. These systems yield a highly non-linear correlation between the SFR and BHAR, which are non-causally connected via the mass of the host halo. In massive haloes (M200 $\gtrsim 10^{12.5}$Msol) both SFR and BHAR decline on average with a roughly constant scaling of SFR/BHAR $\sim 10^{3}$. Given the complexity of the full SFR--BHAR plane built from multiple behaviours, and from the large dynamic range of BHARs, we find the primary driver of the different observed trends in the <SFR>--BHAR and <BHAR>--SFR relationships are due to sampling considerably different regions of this plane.

We investigate the properties of the galaxies that reionized the Universe and the history of cosmic reionization using the "Evolution and Assembly of GaLaxies and their environments" (EAGLE) cosmological hydrodynamical simulations. We obtain the evolution of the escape fraction of ionizing photons in galaxies assuming that galactic winds create channels through which 20~percent of photons escape when the local surface density of star formation is greater than $0.1$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$. Such threshold behaviour for the generation of winds is observed, and the rare local objects which have such high star formation surface densities exhibit high escape fractions of $\sim 10$ percent. In our model the luminosity-weighted mean escape fraction increases with redshift as $\bar f_{\rm esc}=0.045~((1+z)/4)^{1.1}$ at $z>3$, and the galaxy number weighted mean as $\langle f_{\rm esc} \rangle=2.2\times10^{-3}~((1+z)/4)^4$, and becomes constant $\approx0.2$ at redshift $z>10$. The escape fraction evolves as an increasingly large fraction of stars forms above the critical surface density of star formation at earlier times. This evolution of the escape fraction, combined with that of the star formation rate density from EAGLE, reproduces the inferred evolution of the filling factor of ionized regions during the reionization epoch ($6<z<8$), the evolution of the post-reionization ($0\leq z<6$) hydrogen photoionisation rate, and the optical depth due to Thomson scattering of the cosmic microwave background photons measured by the Planck satellite.

We report the detection of extended Halpha emission from the tip of the HI disk of the nearby edge-on galaxy UGC 7321, observed with the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope. The Halpha surface brightness fades rapidly where the HI column density drops below N(HI) = 10^19 cm^-2 , consistent with fluorescence arising at the ionisation front from gas that is photoionized by the extragalactic ultraviolet background (UVB). The surface brightness measured at this location is (1.2 +/- 0.5)x10^-19 erg/s/cm^2/arcsec^2, where the error is mostly systematic and results from the proximity of the signal to the edge of the MUSE field of view, and from the presence of a sky line next to the redshifted Halpha wavelength. By combining the Halpha and the HI 21 cm maps with a radiative transfer calculation of an exponential disk illuminated by the UVB, we derive a value for the HI photoionization rate of Gamma ~ (6-8)x10^-14 1/s . This value is consistent with transmission statistics of the Lyalpha forest and with recent models of a UVB which is dominated by quasars.

We present a MUSE and KMOS dynamical study 405 star-forming galaxies at redshift z=0.28-1.65 (median redshift z=0.84). Our sample are representative of star-forming, main-sequence galaxies, with star-formation rates of SFR=0.1-30Mo/yr and stellar masses M=10^8-10^11Mo. For 49+/-4% of our sample, the dynamics suggest rotational support, 24+/-3% are unresolved systems and 5+/-2% appear to be early-stage major mergers with components on 8-30kpc scales. The remaining 22+/-5% appear to be dynamically complex, irregular (or face-on systems). For galaxies whose dynamics suggest rotational support, we derive inclination corrected rotational velocities and show these systems lie on a similar scaling between stellar mass and specific angular momentum as local spirals with j*=J/M*\propto M^(2/3) but with a redshift evolution that scales as j*\propto M^{2/3}(1+z)^(-1). We identify a correlation between specific angular momentum and disk stability such that galaxies with the highest specific angular momentum, log(j*/M^(2/3))>2.5, are the most stable, with Toomre Q=1.10+/-0.18, compared to Q=0.53+/-0.22 for galaxies with log(j*/M^(2/3))<2.5. At a fixed mass, the HST morphologies of galaxies with the highest specific angular momentum resemble spiral galaxies, whilst those with low specific angular momentum are morphologically complex and dominated by several bright star-forming regions. This suggests that angular momentum plays a major role in defining the stability of gas disks: at z~1, massive galaxies that have disks with low specific angular momentum, appear to be globally unstable, clumpy and turbulent systems. In contrast, galaxies with high specific angular have evolved in to stable disks with spiral structures.

We use the hydrodynamic, cosmological EAGLE simulations to investigate how hot gas in haloes condenses to form and grow galaxies. We select haloes from the simulations that are actively cooling and study the temperature, distribution, and metallicity of their hot, cold, and transitioning `cooling' gas, placing these in context of semi-analytic models. Our selection criteria lead us to focus on Milky Way-like haloes. We find the hot-gas density profiles of the haloes form a progressively stronger core over time, the nature of which can be captured by a beta profile that has a simple dependence on redshift. In contrast, the hot gas that will cool over a time-step is broadly consistent with a singular isothermal sphere. We find that cooling gas carries a few times the specific angular momentum of the halo and is offset in spin direction from the rest of the hot gas. The gas loses ~60% of its specific angular momentum during the cooling process, generally remaining greater than that of the halo, and it precesses to become aligned with the cold gas already in the disc. We find tentative evidence that angular-momentum losses are slightly larger when gas cools onto dispersion-supported galaxies. We show that an exponential surface density profile for gas arriving on a disc remains a reasonable approximation, but a cusp containing ~20% of the mass is always present, and disc scale radii are larger than predicted by a vanilla Fall & Efstathiou model. These scale radii are still closely correlated with the halo spin parameter, for which we suggest an updated prescription for galaxy formation models.

Galaxies that have fallen into massive haloes may no longer be able to accrete gas from their surroundings, a process referred to as 'starvation' or 'strangulation' of satellites. We study the environmental dependence of gas accretion onto galaxies using the cosmological, hydrodynamical EAGLE simulation. We quantify the dependence of gas accretion on stellar mass, redshift, and environment, using halo mass and galaxy overdensity as environmental indicators. We find a strong suppression, by many orders of magnitude, of the gas accretion rate in dense environments, primarily for satellite galaxies. This suppression becomes stronger at lower redshift. However, the scatter in accretion rates is very large for satellites. This is (at least in part) due to the variation in halocentric radius, since gas accretion is more suppressed at smaller radii. Central galaxies are influenced less strongly by their environment and exhibit less scatter in their gas accretion rates. The star formation rates of both centrals and satellites show similar behaviour to their gas accretion rates. The relatively small differences between gas accretion and star formation rates demonstrate that galaxies generally exhaust their gas reservoir somewhat faster at higher stellar mass, lower redshift, and in denser environments. We conclude that the environmental suppression of gas accretion could directly result in the quenching of star formation.

We examine the baryon content of low-mass {\Lambda}CDM halos $(10^8<M_{200}/{\rm M_\odot}<5\times 10^{9})$ using the APOSTLE cosmological hydrodynamical simulations. Most of these systems are free of stars and have a gaseous content set by the combined effects of cosmic reionization, which imposes a mass-dependent upper limit, and of ram pressure stripping, which reduces it further in high-density regions. Halos mainly affected by reionization (RELHICs; REionization-Limited HI Clouds) inhabit preferentially low-density regions and make up a population where the gas is in hydrostatic equilibrium with the dark matter potential and in thermal equilibrium with the ionizing UV background. Their thermodynamic properties are well specified, and their gas density and temperature profiles may be predicted in detail. Gas in RELHICs is nearly fully ionized but with neutral cores that span a large range of HI masses and column densities and have negligible non-thermal broadening. We present predictions for their characteristic sizes and central column densities: the massive tail of the distribution should be within reach of future blind HI surveys. Local Group RELHICs (LGRs) have some properties consistent with observed Ultra Compact High Velocity Clouds (UCHVCs) but the sheer number of the latter suggests that most UCHVCs are not RELHICs. Our results suggest that LGRs (i) should typically be beyond 500 kpc from the Milky Way or M31; (ii) have positive Galactocentric radial velocities; (iii) HI sizes not exceeding 1 kpc, and (iv) should be nearly round. The detection and characterization of RELHICs would offer a unique probe of the small-scale clustering of cold dark matter.

We use the hydrodynamical EAGLE simulation to study the magnitude and origin of the scatter in the stellar mass - halo mass relation for central galaxies. We separate cause and effect by correlating stellar masses in the baryonic simulation with halo properties in a matched dark matter only (DMO) simulation. The scatter in stellar mass increases with redshift and decreases with halo mass. At $z = 0.1$ it declines from 0.25 dex at $M_{200, \rm DMO} \approx 10^{11}$ M$_{\odot}$ to 0.12 dex at $M_{200, \rm DMO} \approx 10^{13}$ M$_{\odot}$, but the trend is weak above $10^{12}$ M$_{\odot}$. For $M_{200, \rm DMO} < 10^{12.5}$ M$_{\odot}$ up to 0.04 dex of the scatter is due to scatter in the halo concentration. At fixed halo mass, a larger stellar mass corresponds to a more concentrated halo. This is likely because higher concentrations imply earlier formation times and hence more time for accretion and star formation, and/or because feedback is less efficient in haloes with higher binding energies. The maximum circular velocity, $V_{\rm max, DMO}$, and binding energy are therefore more fundamental properties than halo mass, meaning that they are more accurate predictors of stellar mass, and we provide fitting formulae for their relations with stellar mass. However, concentration alone cannot explain the total scatter in the $M_{\rm star} - M_{200, \rm DMO}$ relation, and it does not explain the scatter in $M_{\rm star} -V_{\rm max, DMO}$. Halo spin, sphericity, triaxiality, substructure and environment are also not responsible for the remaining scatter, which thus could be due to more complex halo properties or non-linear/stochastic baryonic effects.

Galaxies fall into two clearly distinct types: `blue-sequence' galaxies that are rapidly forming young stars, and `red-sequence' galaxies in which star formation has almost completely ceased. Most galaxies more massive than $3\times10^{10} M_\odot$ follow the red-sequence while less massive central galaxies lie on the blue sequence. We show that these sequences are created by a competition between star formation-driven outflows and gas accretion on to the supermassive black hole at the galaxy's center. We develop a simple analytic model for this interaction. In galaxies less massive than $3\times10^{10} M_\odot$, young stars and supernovae drive a high entropy outflow that is more buoyant that any diffuse corona. The outflow balances the rate of gas inflow, preventing high gas densities building up in the central regions. More massive galaxies, however, are surrounded by a hot corona. We argue that above a halo mass of $\sim 10^{12} M_\odot$, the supernova-driven outflow is no longer buoyant and star formation is unable to prevent the build up of gas in the central regions. This triggers a strongly non-linear response from the black hole. Its accretion rate rises rapidly, heating the galaxy's corona, disrupting the incoming supply of cool gas and starving the galaxy of the fuel for star formation. The host galaxy makes a transition to the red sequence, and further growth predominantly occurs through galaxy mergers. We show that the analytic model provides a good description of galaxy evolution in the EAGLE hydrodynamic simulations, and demonstrate that, so long as star formation-driven outflows are present, the transition mass scale is almost independent of subgrid parameter choice. The transition mass disappears entirely, however, if star formation driven outflows are absent.

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

We use the EAGLE cosmological hydrodynamic simulation suite to study the specific angular momentum of galaxies, $j$, with the aims of (i) investigating the physical causes behind the wide range of $j$ at fixed mass and (ii) examining whether simple, theoretical models can explain the seemingly complex and non-linear nature of the evolution of $j$. We find that $j$ of the stars, $j_{\rm stars}$, and baryons, $j_{\rm bar}$, are strongly correlated with stellar and baryon mass, respectively, with the scatter being highly correlated with morphological proxies such as gas fraction, stellar concentration, (u-r) intrinsic colour, stellar age and the ratio of circular velocity to velocity dispersion. We compare with available observations at $z=0$ and find excellent agreement. We find that $j_{\rm bar}$ follows the theoretical expectation of an isothermal collapsing halo under conservation of specific angular momentum to within $\approx 50$%, while the subsample of rotation-supported galaxies are equally well described by a simple model in which the disk angular momentum is just enough to maintain marginally stable disks. We extracted evolutionary tracks of the stellar spin parameter of EAGLE galaxies and found that the fate of their $j_{\rm stars}$ at $z=0$ depends sensitively on their star formation and merger histories. From these tracks, we identified two distinct physical channels behind low $j_{\rm stars}$ galaxies at $z=0$: (i) galaxy mergers, and (ii) early star formation quenching. The latter can produce galaxies with low $j_{\rm stars}$ and early-type morphologies even in the absence of mergers.

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ' star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor's mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM.

Observations of galaxies in the local Universe have shown that both the ionized gas and the stars of satellites are more metal-rich than of equally massive centrals. To gain insight into the connection between this metallicity enhancement and other differences between centrals and satellites, such as their star formation rates, gas content, and growth history, we study the metallicities of >3600 galaxies with M_star > 10^10 M_sun in the cosmological hydrodynamical EAGLE 100 Mpc `Reference' simulation, including ~1500 in the vicinity of galaxy groups and clusters (M_200 >= 10^13 M_sun). The simulation predicts excess gas and stellar metallicities in satellites consistent with observations, except for stellar metallicities at M_star <~ 10^10.2 M_sun where the predicted excess is smaller than observed. The exact magnitude of the effect depends on galaxy selection, aperture, and on whether the metallicity is weighted by stellar mass or luminosity. The stellar metallicity excess in clusters is also sensitive to the efficiency scaling of star formation feedback. We identify stripping of low-metallicity gas from the galaxy outskirts, as well as suppression of metal-poor inflows towards the galaxy centre, as key drivers of the enhancement of gas metallicity. Stellar metallicities in satellites are higher than in the field as a direct consequence of the more metal-rich star forming gas, whereas stripping of stars and suppressed stellar mass growth, as well as differences in accreted vs. in-situ star formation between satellites and the field, are of secondary importance.

Stellar population synthesis models can now reproduce the photometry of old stellar systems (age $>$ 2 Gyr) in the near-infrared (NIR) bands at 3.4 and 4.6$\mu$m (WISE W1 $\&$ W2 or IRAC 1 $\&$ 2). In this paper we derive stellar mass-to-light ratios for these and optical bands, and confirm that the NIR M/L shows dramatically reduced sensitivity to both age and metallicity compared to optical bands, and further, that this behavior leads to significantly more robust stellar masses for quiescent galaxies with [Fe/H] > -0.5 regardless of star formation history (SFH). We then use realistic early-type galaxy SFHs and metallicity distributions from the EAGLE simulations of galaxy formation to investigate two methods to determine the appropriate M/L for a galaxy: 1) We show that the uncertainties introduced by an unknown SFH can be largely removed using a spectroscopically inferred luminosity-weighted age and metallicity for the population to select the appropriate single stellar population (SSP) equivalent M/L. Using this method, the maximum systematic error due to SFH on the M/L of an early-type galaxy is $<$ 4$\%$ at 3.4 $\mu$m and typical uncertainties due to errors in the age and metallicity create scatter of $\lesssim$13$\%$. The equivalent values for optical bands are more than 2-3 times greater, even before considering uncertainties associated with internal dust extinction. 2) We demonstrate that if the EAGLE SFHs and metallicities accurately reproduce the true properties of early-type galaxies, the use of an iterative approach to select a mass dependent M/L can provide even more accurate stellar masses for early-type galaxies, with typical uncertainties $<$ 9$\%$.

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

We study the $z\approx3.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 report on the search for galaxies in the proximity of two very metal-poor gas clouds at z~3 towards the quasar Q0956+122. With a 5-hour MUSE integration in a ~500x500 kpc^2 region centred at the quasar position, we achieve a >80% complete spectroscopic survey of continuum-detected galaxies with m<25 mag and Ly{\alpha} emitters with luminosity L>3e41 erg/s. We do not identify galaxies at the redshift of a z~3.2 Lyman limit system (LLS) with log Z/Zsun = -3.35 +/- 0.05, placing this gas cloud in the intergalactic medium or circumgalactic medium of a galaxy below our sensitivity limits. Conversely, we detect five Ly{\alpha} emitters at the redshift of a pristine z~3.1 LLS with log Z/Zsun < -3.8, while ~0.4 sources were expected given the z~3 Ly{\alpha} luminosity function. Both this high detection rate and the fact that at least three emitters appear aligned in projection with the LLS suggest that this pristine cloud is tracing a gas filament that is feeding one or multiple galaxies. Our observations uncover two different environments for metal-poor LLSs, implying a complex link between these absorbers and galaxy halos, which ongoing MUSE surveys will soon explore in detail. Moreover, in agreement with recent MUSE observations, we detected a ~90 kpc Ly{\alpha} nebula at the quasar redshift and three Ly{\alpha} emitters reminiscent of a "dark galaxy" population.

Subhalo abundance matching (SHAM) is a widely-used method to connect galaxies with dark matter structures in numerical simulations. SHAM predictions agree remarkably well with observations, yet they still lack strong theoretical support. We examine the performance, implementation, and assumptions of SHAM using the EAGLE project simulations. We find that $V_{\rm relax}$, the highest value of the circular velocity attained by a subhalo while it satisfies a relaxation criterion, is the subhalo property that correlates most strongly with galaxy stellar mass ($M_{\rm star}$). Using this parameter in SHAM, we retrieve the real-space clustering of EAGLE to within our statistical uncertainties on scales greater than $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<10.77$. Conversely, clustering is overestimated by $30\%$ on scales below $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<9.77$ because SHAM slightly overpredicts the fraction of satellites in massive haloes compared to EAGLE. The agreement is even better in redshift-space, where the clustering is recovered to within our statistical uncertainties for all masses and separations. Additionally, we analyse the dependence of galaxy clustering on properties other than halo mass, i.e. the assembly bias. We demonstrate assembly bias alters the clustering in EAGLE by $20\%$ and that $V_{\rm relax}$ captures its effect to within $15\%$. We trace small differences in the clustering to the failure of SHAM as typically implemented, i.e. the $M_{\rm star}$ assigned to a subhalo does not depend on i) its host halo mass, ii) whether it is a central or a satellite. In EAGLE we find that these assumptions are not completely satisfied.

We estimate the expected event rate of gravitational wave signals from mergers of supermassive black holes that could be resolved by a space-based interferometer, such as the Evolved Laser Interferometer Space Antenna (eLISA), utilising the reference cosmological hydrodynamical simulation from the EAGLE suite. These simulations assume a $\Lambda$CDM cosmogony with state-of-the-art subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. They have been shown to reproduce the observed galaxy population with unprecedented fidelity. We combine the merger rates of supermassive black holes in EAGLE with the latest phenomenological waveform models to calculate the gravitational waves signals from the intrinsic parameters of the merging black holes. The EAGLE models predict $\sim2$ detections per year by a gravitational wave detector such as eLISA. We find that these signals are largely dominated by mergers between seed mass black holes merging at redshifts between $z\sim2$ and $z\sim1$. In order to investigate the dependence on the assumed black hole seed mass, we introduce an additional model with a black hole seed mass an order of magnitude smaller than in our reference model. We also consider a variation of the reference model where a prescription for the expected delays in the black hole merger timescale has been included after their host galaxies merge. We find that the merger rate is similar in all models, but that the initial black hole seed mass could be distinguished through their detected gravitational waveforms. Hence, the characteristic gravitational wave signals detected by eLISA will provide profound insight into the origin of supermassive black holes and the initial mass distribution of black hole seeds.

The EAGLE cosmological simulations reproduce the observed galaxy stellar mass function and many galaxy properties. In this work, we study the dust-related properties of present-day EAGLE galaxies through mock observations in the far-infrared and submm wavelength ranges obtained with the 3D dust radiative transfer code SKIRT. To prepare an EAGLE galaxy for radiative transfer processing, we derive a diffuse dust distribution from the gas particles and we re-sample the star-forming gas particles and the youngest star particles into star-forming regions that are assigned dedicated emission templates. We select a set of redshift-zero EAGLE galaxies that matches the K-band luminosity distribution of the galaxies in the Herschel Reference Survey (HRS), a volume-limited sample of about 300 normal galaxies in the Local Universe. We find overall agreement of the EAGLE dust scaling relations with those observed in the HRS, such as the dust-to-stellar mass ratio versus stellar mass and versus NUV-r colour relations. A discrepancy in the f_250/f_350 versus f_350/f_500 submm colour-colour relation implies that part of the simulated dust is insufficiently heated, likely because of limitations in our sub-grid model for star-forming regions. We also investigate the effect of adjusting the metal-to-dust ratio and the covering factor of the photodissociation regions surrounding the star-forming cores. We are able to constrain the important dust-related parameters in our method, informing the calculation of dust attenuation for EAGLE galaxies in the UV and optical domain.

The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $\Lambda$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $\Lambda$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution.

We investigate the evolution of supermassive black holes in the `Evolution and Assembly of GaLaxies and their Environments' (EAGLE) cosmological hydrodynamic simulations. The largest of the EAGLE volumes covers a $(100 \,\rm cMpc)^3$ and includes state-of-the-art physical models for star formation and black hole growth that depend only on local gas properties. We focus on the black hole mass function, Eddington ratio distribution and the implied duty cycle of nuclear activity. The simulation is broadly consistent with observational constraints on these quantities. In order to make a more direct comparison with observational data, we calculate the soft and hard X-ray luminosity functions of the active galactic nuclei (AGN). Between redshifts $0$ and $1$, the simulation is in agreement with data. At higher redshifts, the simulation tends to underpredict the luminosities of the brightest observed AGN. This may be due to the limited volume of the simulation, or a fundamental deficiency of the underlying model. It seems unlikely that additional unresolved variability can account for this difference. The simulation shows a similar `downsizing' of the AGN population as seen in observational surveys.

The spatial distribution of the satellite populations of the Milky Way and Andromeda are puzzling in that they are nearly perpendicular to the disks of their central galaxies. To understand the origin of such configurations we study the alignment of the central galaxy, satellite system and dark matter halo in the largest of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation. We find that centrals and their satellite systems tend to be well aligned with their haloes, with a median misalignment angle of $33^{\circ}$ in both cases. While the centrals are better aligned with the inner $10$ kpc halo, the satellite systems are better aligned with the entire halo indicating that satellites preferentially trace the outer halo. The central - satellite alignment is weak (median misalignment angle of $52^{\circ}$) and we find that around $20\%$ of systems have a misalignment angle larger than $78^{\circ}$, which is the value for the Milky Way. The central - satellite alignment is a consequence of the tendency of both components to align with the dark matter halo. As a consequence, when the central is parallel to the satellite system, it also tends to be parallel to the halo. In contrast, if the central is perpendicular to the satellite system, as in the case of the Milky Way and Andromeda, then the central - halo alignment is much weaker. Dispersion-dominated (spheroidal) centrals have a stronger alignment with both their halo and their satellites than rotation-dominated (disk) centrals. We also found that the halo, the central galaxy and the satellite system tend to be aligned with the surrounding large-scale distribution of matter, with the halo being the better aligned of the three.

We examine the evolution of intrinsic u-r colours of galaxies in the EAGLE cosmological hydrodynamical simulations, which has been shown to reproduce the observed redshift z=0.1 colour-magnitude distribution well. The median u-r of star-forming ('blue cloud') galaxies reddens by 1 mag from z=2 to 0 at fixed stellar mass, as their specific star formation rates decrease with time. A red sequence starts to build-up around z=1, due to the quenching of low-mass satellite galaxies at the faint end, and due to the quenching of more massive central galaxies by their active galactic nuclei (AGN) at the bright end. This leaves a dearth of intermediate-mass red sequence galaxies at z=1, which is mostly filled in by z=0. We quantify the time-scales of colour transition due to satellite and AGN quenching, finding that most galaxies spend less than 2 Gyr in the 'green valley'. On examining the trajectories of galaxies in a colour-stellar mass diagram, we identify three characteristic tracks that galaxies follow (quiescently star-forming, quenching and rejuvenating galaxies) and quantify the fraction of galaxies that follow each track.

We use the EAGLE suite of cosmological hydrodynamical simulations to study how the HI content of present-day galaxies depends on their environment. We show that EAGLE reproduces observed HI mass-environment trends very well, while semi-analytic models typically overpredict the average HI masses in dense environments. The environmental processes act primarily as an on/off switch for the HI content of satellites with stellar mass Mstar>10^9 Msun. At a fixed Mstar, the fraction of HI-depleted satellites increases with increasing host halo mass M200 in response to stronger environmental effects, while at a fixed M200 it decreases with increasing satellite Mstar as the gas is confined by deeper gravitational potentials. HI-depleted satellites reside mostly, but not exclusively, within the virial radius r200 of their host halo. We investigate the origin of these trends by focussing on three environmental mechanisms: ram pressure stripping by the intra-group medium, tidal stripping by the host halo, and satellite-satellite encounters. By tracking back in time the evolution of the HI-depleted satellites, we find that the most common cause of HI removal is satellite encounters. The timescale for HI removal is typically less than 0.5 Gyr. Tidal stripping occurs in halos of M200<10^{14} Msun within half r200, while the other processes act also in more massive halos, generally within r200. Conversely, we find that ram pressure stripping is the most common mechanism that disturbs the HI morphology of galaxies at redshift z=0. This implies that HI removal due to satellite-satellite interactions occurs on shorter timescales than the other processes.

We study the origin of the stellar $\alpha$-element-to-iron abundance ratio, $[\alpha/\mathrm{Fe}]_{\ast}$, of present-day central galaxies, using cosmological, hydrodynamical simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. For galaxies with stellar masses of $M_{\ast} > 10^{10.5}$ M$_{\odot}$, $[\alpha/\mathrm{Fe}]_{\ast}$ increases with increasing galaxy stellar mass and age. These trends are in good agreement with observations of early-type galaxies, and are consistent with a `downsizing' galaxy formation scenario: more massive galaxies have formed the bulk of their stars earlier and more rapidly, hence from an interstellar medium that was mostly $\alpha$-enriched by massive stars. In the absence of feedback from active galactic nuclei (AGN), however, $[\alpha/\mathrm{Fe}]_{\ast}$ in $M_{\ast} > 10^{10.5}$ M$_{\odot}$ galaxies is roughly constant with stellar mass and decreases with mean stellar age, extending the trends found for lower-mass galaxies in both simulations with and without AGN. We conclude that AGN feedback can account for the $\alpha$-enhancement of massive galaxies, as it suppresses their star formation, quenching more massive galaxies at earlier times, thereby preventing the iron from longer-lived intermediate-mass stars (supernova Type Ia) from being incorporated into younger stars.

We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

Observations of local galaxies harbouring supermassive black holes (BHs) of anomalously high mass, M_BH, relative to their stellar mass, M_star, appear to be at odds with simple models of the co-evolution between galaxies and their central BHs. We study the origin of such outliers in a Lambda cold dark matter context using the EAGLE cosmological, hydrodynamical simulation. We find 15 'M_BH(M_star)-outlier' galaxies, defined as having M_BH more than 1.5 dex above the median M_BH(M_star) relation in the simulation, M_{BH,med}. All M_BH(M_star)-outliers are satellite galaxies, typically with M_star ~ 10^10 M_sun and M_BH ~ 10^8 M_sun. They have all become outliers due to a combination of tidal stripping of their outer stellar component acting over several Gyr and early formation times leading to rapid BH growth at high redshift, with the former mechanism being most important for 67 per cent of these outliers. The same mechanisms also cause the M_BH(M_star)-outlier satellites to be amongst the most compact galaxies in the simulation, making them ideal candidates for ultracompact dwarf galaxy progenitors. The 10 most extreme central galaxies found at z=0 (with log_{10}(M_BH/M_{BH,med}) in [1.2, 1.5]) grow rapidly in M_BH to lie well above the present-day M_BH-M_star relation at early times (z > 2), and either continue to evolve parallel to the z=0 relation or remain unchanged until the present day, making them 'relics' of the high-redshift universe. This high-z formation mechanism may help to explain the origin of observed M_BH(M_star)-outliers with extended dark matter haloes and undisturbed morphologies.

We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the Solar neighbourhood for this set of Milky Way analogues, and use them to analyse the results of current direct detection experiments. For most Milky Way analogues, the event rates in direct detection experiments obtained from the best fit Maxwellian distribution (with peak speed of 223 - 289 km/s) are similar to those obtained directly from the simulations. As a consequence, the allowed regions and exclusion limits set by direct detection experiments in the dark matter mass and spin-independent cross section plane shift by a few GeV compared to the Standard Halo Model, at low dark matter masses. For each dark matter mass, the halo-to-halo variation of the local dark matter density results in an overall shift of the allowed regions and exclusion limits for the cross section. However, the compatibility of the possible hints for a dark matter signal from DAMA and CDMS-Si and null results from LUX and SuperCDMS is not improved.

We use the EAGLE galaxy formation simulation to study the effects of baryons on the power spectrum of the total matter and dark matter distributions and on the velocity fields of dark matter and galaxies. On scales $k{\stackrel{>}{{}_\sim}} 4{h\,{\rm Mpc}^{-1}}$ the effect of baryons on the amplitude of the total matter power spectrum is greater than $1\%$. The back-reaction of baryons affects the density field of the dark matter at the level of $\sim3\%$ on scales of $1\leq k/({h\,{\rm Mpc}^{-1}})\leq 5$. The dark matter velocity divergence power spectrum at $k{\stackrel{<}{{}_\sim}}0.5{h\,{\rm Mpc}^{-1}}$ is changed by less than $1\%$. The 2D redshift-space power spectrum is affected at the level of $\sim6\%$ at $|\vec{k}|{\stackrel{>}{{}_\sim}} 1{h\,{\rm Mpc}^{-1}}$ (for $\mu>0.5$), but for $|\vec{k}|\leq 0.4{h\,{\rm Mpc}^{-1}}$ it differs by less than $1\%$. We report vanishingly small baryonic velocity bias for haloes: the peculiar velocities of haloes with $M_{200}>3\times10^{11}{{\rm M}_{\odot}}$ (hosting galaxies with $M_{*}>10^9{{\rm M}_{\odot}}$) are affected at the level of at most $1~$km/s, which is negligible for $1\%$-precision cosmology. We caution that since EAGLE overestimates cluster gas fractions it may also underestimate the impact of baryons, particularly for the total matter power spectrum. Nevertheless, our findings suggest that for theoretical modelling of redshift space distortions and galaxy velocity-based statistics, baryons and their back-reaction can be safely ignored at the current level of observational accuracy. However, we confirm that the modelling of the total matter power spectrum in weak lensing studies needs to include realistic galaxy formation physics in order to achieve the accuracy required in the precision cosmology era.

We present new results on the auto- and cross-correlation functions of galaxies and OVI absorbers in a $\sim 18~\textrm{Gpc}^3$ comoving volume at $z < 1$. We use a sample of 51,296 galaxies and 140 OVI absorbers in the column density range $13 \lesssim \log N \lesssim 15$ to measure two-point correlation functions in the two dimensions transverse and orthogonal to the line-of-sight $\xi(r_{\perp}, r_{\parallel})$. We furthermore infer the corresponding 'real-space' correlation functions, $\xi(r)$, by projecting $\xi(r_{\perp}, r_{\parallel})$ along $r_{\parallel}$, and assuming a power-law form, $\xi(r) = (r / r_0)^{-\gamma}$. Comparing the results from the absorber-galaxy cross-correlation function, $\xi_{\textrm{ag}}$, the galaxy auto-correlation function, $\xi_{\textrm{gg}}$, and the absorber auto-correlation function, $\xi_{\textrm{aa}}$, we constrain the statistical connection between galaxies and the metal-enriched intergalactic medium as a function of star-formation activity. We also compare these results to predictions from the EAGLE cosmological hydrodynamical simulation and find a reasonable agreement. We find that: (i) OVI absorbers show very little velocity dispersion with respect to galaxies on $\sim$ Mpc scales, likely $\lesssim$ 100 \kms; (ii) OVI absorbers and galaxies may not linearly trace the same underlying distribution of matter in general. In particular, our results demonstrate that OVI absorbers are less clustered, and potentially more extended around galaxies than galaxies are around themselves; (iii) On $\gtrsim 100$ kpc scales, the likelihood of finding OVI absorbers around star-forming galaxies is similar to the likelihood of finding OVI absorbers around non star-forming galaxies (abridged)

We study the distribution and evolution of highly ionised intergalactic metals in the Evolution and Assembly of Galaxies and their Environment (EAGLE) cosmological, hydrodynamical simulations. EAGLE has been shown to reproduce a wide range of galaxy properties while its subgrid feedback was calibrated without considering gas properties. We compare the predictions for the column density distribution functions (CDDFs) and cosmic densities of SiIV, CIV, NV, OVI and NeVIII absorbers with observations at redshift z = 0 to ~ 6 and find reasonable agreement, although there are some differences. We show that the typical physical densities of the absorbing gas increase with column density and redshift, but decrease with the ionization energy of the absorbing ion. The typical metallicity increases with both column density and time. The fraction of collisionally ionized metal absorbers increases with time and ionization energy. While our results show little sensitivity to the presence or absence of AGN feedback, increasing/decreasing the efficiency of stellar feedback by a factor of two substantially decreases/increases the CDDFs and the cosmic densities of the metal ions. We show that the impact of the efficiency of stellar feedback on the CDDFs and cosmic densities is largely due to its effect on the metal production rate. However, the temperatures of the metal absorbers, particularly those of strong OVI, are directly sensitive to the strength of the feedback.

We investigate correlations between different physical properties of star-forming galaxies in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamical simulation suite over the redshift range $0\le z\le 4.5$. A principal component analysis reveals that neutral gas fraction ($f_{\rm gas, neutral}$), stellar mass ($M_{\rm stellar}$) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of $f_{\rm gas, neutral}-M_{\rm stellar}-\rm SFR$ with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their $f_{\rm gas, neutral}$ and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from $f_{\rm gas, neutral}$, or from the $M_{\rm stellar}$ and SFR. We argue that the appearance of this "fundamental plane of star formation" is a consequence of self-regulation, with the plane's curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range $0\lesssim \rm z\lesssim 2.5$, and find that such a plane is also present in the data. The properties of the observed fundamental plane of star formation are in good agreement with EAGLE's predictions.

Hydrogen in the Universe was (re)ionised between redshifts $z \approx 10$ and $z \approx 6$. The nature of the sources of the ionising radiation is hotly debated, with faint galaxies below current detection limits regarded as prime candidates. Here we consider a scenario in which ionising photons escape through channels punctured in the interstellar medium by outflows powered by starbursts. We take account of the observation that strong outflows occur only when the star formation density is sufficiently high, and estimate the galaxy-averaged escape fraction as a function of redshift and luminosity from the resolved star formation surface densities in the EAGLE cosmological hydrodynamical simulation. We find that the fraction of ionising photons that escape from galaxies increases rapidly with redshift, reaching values of 5-20 percent at $z > 6$, with the brighter galaxies having higher escape fractions. Combining the dependence of escape fraction on luminosity and redshift with the observed luminosity function, we demonstrate that galaxies emit enough ionising photons to match the existing constraints on reionisation while also matching the observed UV-background post-reionisation. Our findings suggest that galaxies above the current Hubble Space Telescope detection limit emit half of the ionising radiation required to reionise the Universe.

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.

We explore the co-evolution of the specific angular momentum of dark matter haloes and the cold baryons that comprise the galaxies within. We study over two thousand central galaxies within the reference cosmological hydrodynamical simulation of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project. We employ a methodology within which the evolutionary history of a system is specified by the time-evolving properties of the Lagrangian particles that define it at z=0. We find a strong correlation between the evolution of the specific angular momentum of today's stars (cold gas) and that of the inner (whole) dark matter halo they are associated with. This link is particularly strong for the stars formed before the epoch of maximum expansion and subsequent collapse of the central dark matter halo (turnaround). Spheroids are typically assembled primarily from stars formed prior to turnaround, and are therefore destined to suffer a net loss of angular momentum associated with the strong merging activity during the assembly of the inner dark matter halo. Stellar discs retain their specific angular momentum since they are comprised of stars formed mainly after turnaround, from gas that mostly preserves the high specific angular momentum it acquired by tidal torques during the linear growth of the halo. Since the specific angular momentum loss of the stars is tied to the galaxy's morphology today, it may be possible to use our results to predict, statistically, the assembly history of a halo given the morphology of the galaxy it hosts.

We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $r\approx1.5$ kpc and $r\approx10$ kpc. At smaller radii, $r\lesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $\gamma$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.

We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical prescription from the literature. The resulting neutral hydrogen fractions, M_(HI+H2)/M_star, agree with observations to better than 0.1 dex for galaxies with M_star between 10^10 and 10^11 M_Sun. Our fiducial, empirical H2 model based on gas pressure results in galactic HI mass fractions, M_HI/M_star, that agree with observations from the GASS survey to better than 0.3 dex, but the alternative theoretical H2 formula leads to a negative offset in M_HI/M_star of up to 0.5 dex. Visual inspection reveals that most HI disks in simulated HI-rich galaxies are vertically disturbed, plausibly due to recent accretion events. Many galaxies (up to 80 per cent) contain spuriously large HI holes, which are likely formed as a consequence of the feedback implementation in EAGLE. The HI mass-size relation of all simulated galaxies is close to (but 16 per cent steeper than) observed, and when only galaxies without large holes in the HI disc are considered, the agreement becomes excellent (better than 0.1 dex). The presence of large HI holes also makes the radial HI surface density profiles somewhat too low in the centre, at \Sigma_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (\Sigma_HI < 1 M_Sun pc^-2), the simulated profiles agree quantitatively with observations. Scaled by HI size, the simulated profiles of HI-rich (M_HI > 10^9.8 M_Sun) and control galaxies (10^9.1 M_Sun > M_HI > 10^9.8 M_Sun) follow each other closely, as observed. (Abridged)

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" ({\Lambda}CDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in {\Lambda}CDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

We use cosmological, hydrodynamical simulations from the EAGLE and OWLS projects to assess the significance of recycled stellar ejecta as fuel for star formation. The fractional contributions of stellar mass loss to the cosmic star formation rate (SFR) and stellar mass densities increase with time, reaching $35 \%$ and $19 \%$, respectively, at $z=0$. The importance of recycling increases steeply with galaxy stellar mass for $M_{\ast} < 10^{10.5}$ M$_{\odot}$, and decreases mildly at higher mass. This trend arises from the mass dependence of feedback associated with star formation and AGN, which preferentially suppresses star formation fuelled by recycling. Recycling is more important for satellites than centrals and its contribution decreases with galactocentric radius. The relative contribution of AGB stars increases with time and towards galaxy centers. This is a consequence of the more gradual release of AGB ejecta compared to that of massive stars, and the preferential removal of the latter by star formation-driven outflows and by lock up in stellar remnants. Recycling-fuelled star formation exhibits a tight, positive correlation with galaxy metallicity, with a secondary dependence on the relative abundance of alpha elements (which are predominantly synthesized in massive stars), that is insensitive to the subgrid models for feedback. Hence, our conclusions are directly relevant for the origin of the mass-metallicity relation and metallicity gradients. Applying the relation between recycling and metallicity to the observed mass-metallicity relation yields our best estimate of the mass-dependent contribution of recycling. For centrals with a mass similar to that of the Milky Way, we infer the contributions of recycled stellar ejecta to the SFR and stellar mass to be $35 \%$ and $20 \%$, respectively.

We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 1003 cMpc3 simulation of the EAGLE project. We find that the sizes increase with stellar mass , but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size - mass relation is of order 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies to that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size - mass relation due to the evolving number density of passive galaxies. Observations of z \approx 2 galaxies have revealed an abundance of massive red compact galaxies, that depletes below z \approx 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies to their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60% of the compact galaxies increase in size further due to renewed star formation and/or mergers.

Modern analyses of structure formation predict a universe tangled in a 'cosmic web' of dark matter and diffuse baryons. These theories further predict that at low-z, a significant fraction of the baryons will be shock-heated to $T \sim 10^{5}-10^{7}$K yielding a warm-hot intergalactic medium (WHIM), but whose actual existence has eluded a firm observational confirmation. We present a novel experiment to detect the WHIM, by targeting the putative filaments connecting galaxy clusters. We use HST/COS to observe a remarkable QSO sightline that passes within $\Delta d = 3$ Mpc from the 7 inter-cluster axes connecting 7 independent cluster-pairs at redshifts $0.1 \le z \le 0.5$. We find tentative excesses of total HI, narrow HI (NLA; Doppler parameters $b<50$ km/s), broad HI (BLA; $b \ge 50$ km/s) and OVI absorption lines within rest-frame velocities of $\Delta v \lesssim 1000$ km/s from the cluster-pairs redshifts, corresponding to $\sim 2$, $\sim 1.7$, $\sim 6$ and $\sim 4$ times their field expectations, respectively. Although the excess of OVI likely comes from gas close to individual galaxies, we conclude that most of the excesses of NLAs and BLAs are truly intergalactic. We find that the covering fractions, $f_c$, of BLAs close to cluster-pairs are $\sim 4-7$ times higher than the random expectation (at the $\sim 2 \sigma$ c.l.), whereas the $f_c$ of NLAs and OVI are not significantly enhanced. We argue that a larger relative excess of BLAs compared to those of NLAs close to cluster-pairs may be a signature of the WHIM in inter-cluster filaments. By extending the present analysis to tens of sightlines our experiment offers a promising route to detect the WHIM.

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7<log_{10}(M_{200}/[h^{-1}M_\odot])<15$) and a large range of separations ($-1<log_{10}(r/[h^{-1}Mpc])< 2$). We focus on two classes of alignments: the orientations of galaxies with respect to either the directions to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is a strongly decreasing function of the distance between galaxies. The orientation-direction alignment can remain significant up to ~100 Mpc, for galaxies hosted by the most massive haloes in our simulations. Galaxies hosted by more massive subhaloes show stronger alignment. At a fixed halo mass, more aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites is anisotropic and significantly aligned with the major axis of the main host halo. The major axis of satellite galaxies, when all stars are considered, are preferentially aligned towards the centre of the main host halo. The predicted projected direction-orientation alignment, $\epsilon_{g+}(r_{p})$, is in broad agreement with recent observations when only stars within the typical observable extent of a galaxy are used to define galaxy orientations. We find that the orientation-orientation alignment is weaker than the orientation-direction alignment on all scales. Overall, the strength of galaxy alignments depends strongly on the subset of stars that are used to measure the orientations of galaxies and it is always weaker than the alignment of the dark matter haloes. Thus, alignment models that use halo orientation as a direct proxy for galaxy orientation will overestimate the impact of intrinsic alignments on weak lensing analyses.

We present results from a subset of simulations from the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) suite in which the formulation of the hydrodynamics scheme is varied. We compare simulations that use the same subgrid models without re-calibration of the parameters but employing the standard GADGET flavour of smoothed particle hydrodynamics (SPH) instead of the more recent state-of-the-art ANARCHY formulation of SPH that was used in the fiducial EAGLE runs. We find that the properties of most galaxies, including their masses and sizes, are not significantly affected by the details of the hydrodynamics solver. However, the star formation rates of the most massive objects are affected by the lack of phase mixing due to spurious surface tension in the simulation using standard SPH. This affects the efficiency with which AGN activity can quench star formation in these galaxies and it also leads to differences in the intragroup medium that affect the X-ray emission from these objects. The differences that can be attributed to the hydrodynamics solver are, however, likely to be less important at lower resolution. We also find that the use of a time step limiter is important for achieving the feedback efficiency required to match observations of the low-mass end of the galaxy stellar mass function.

Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a "subgrid" model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($\lsim 10^{11.5} \msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations.

Simulations of galaxy formation follow the gravitational and hydrodynamical interactions between gas, stars and dark matter through cosmic time. The huge dynamic range of such calculations severely limits strong scaling behaviour of the community codes in use, with load-imbalance, cache inefficiencies and poor vectorisation limiting performance. The new swift code exploits task-based parallelism designed for many-core compute nodes interacting via MPI using asynchronous communication to improve speed and scaling. A graph-based domain decomposition schedules interdependent tasks over available resources. Strong scaling tests on realistic particle distributions yield excellent parallel efficiency, and efficient cache usage provides a large speed-up compared to current codes even on a single core. SWIFT is designed to be easy to use by shielding the astronomer from computational details such as the construction of the tasks or MPI communication. The techniques and algorithms used in SWIFT may benefit other computational physics areas as well, for example that of compressible hydrodynamics. For details of this open-source project, see www.swiftsim.com

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with the local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations.

We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of ~ 5*10^12 M_sun. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of ~2 times smaller than suggested by the timing argument, but with large dispersion, spanning more than a decade in mass. Guided by these results, we have selected 12 halo pairs with total mass in the range 1.6-3.6 *10^12 M_sun for the APOSTLE project (A Project Of Simulations of The Local Environment), a suite of resimulations at various numerical resolution levels (reaching up to ~10^4 M_sun per gas particle) that use the hydrodynamical code and subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. These diagnostics are sensitive to the total mass assumed for the MW-M31 pair; indeed, the LG satellite population would be quite difficult to reproduce for pair masses as high as indicated by the timing argument. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the Lambda-CDM cosmology with observations of our local Universe.

We examine the circular velocity profiles of galaxies in {\Lambda}CDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect.

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 $z\approx 1.2-1.5$, later than the predicted peak of the cosmic star formation rate activity, at $z\approx 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.

The observed high covering fractions of neutral hydrogen (HI) with column densities above $\sim 10^{17} \rm{cm}^{-2}$ around Lyman-Break Galaxies (LBGs) and bright quasars at redshifts z ~ 2-3 has been identified as a challenge for simulations of galaxy formation. We use the EAGLE cosmological, hydrodynamical simulation, which has been shown to reproduce a wide range of galaxy properties and for which the subgrid feedback was calibrated without considering gas properties, to study the distribution of HI around high-redshift galaxies. We predict the covering fractions of strong HI absorbers ($N_{\rm{HI}} \gtrsim 10^{17} \rm{cm}^{-2}$) inside haloes to increase rapidly with redshift but to depend only weakly on halo mass. For massive ($M_{200} \gtrsim 10^{12} {\rm M_{\odot}}$) halos the covering fraction profiles are nearly scale-invariant and we provide fitting functions that reproduce the simulation results. While efficient feedback is required to increase the HI covering fractions to the high observed values, the distribution of strong absorbers in and around halos of a fixed mass is insensitive to factor of two variations in the strength of the stellar feedback. In contrast, at fixed stellar mass the predicted HI distribution is highly sensitive to the feedback efficiency. The fiducial EAGLE simulation reproduces both the observed global column density distribution function of HI and the observed radial covering fraction profiles of strong HI absorbers around LBGs and bright quasars.

We use the "Evolution and assembly of galaxies and their environments" (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $\rho_*(r) \propto r^{-3}$. Stars dominate the mass density for $r<10~\rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($r\gtrsim3.0~\rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy.

The evolution of the metal content of galaxies and its relations to other global properties [such as total stellar mass (M*), circular velocity, star formation rate (SFR), halo mass, etc.] provides important constraints on models of galaxy formation. Here we examine the evolution of metallicity scaling relations of simulated galaxies in the Galaxies-Intergalactic Medium Interaction Calculation suite of cosmological simulations. We make comparisons to observations of the correlation of gas-phase abundances with M* (the mass-metallicity relation, MZR), as well as with both M* and SFR or gas mass fraction (the so-called 3D fundamental metallicity relations, FMRs). The simulated galaxies follow the observed local MZR and FMRs over an order of magnitude in M*, but overpredict the metallicity of massive galaxies (log M* > 10.5), plausibly due to inefficient feedback in this regime. We discuss the origin of the MZR and FMRs in the context of galactic outflows and gas accretion. We examine the evolution of mass-metallicity relations defined using different elements that probe the three enrichment channels (SNII, SNIa, and AGB stars). Relations based on elements produced mainly by SNII evolve weakly, whereas those based on elements produced preferentially in SNIa/AGB exhibit stronger evolution, due to the longer timescales associated with these channels. Finally, we compare the relations of central and satellite galaxies, finding systematically higher metallicities for satellites, as observed. We show this is due to the removal of the metal poor gas reservoir that normally surrounds galaxies and acts to dilute their gas-phase metallicity (via cooling/accretion onto the disk), but is lost due to ram pressure stripping for satellites.

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $\Lambda$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass.

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation.

We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EAGLE simulation suite using the GALAXEV population synthesis models. We take into account obscuration by dust in birth clouds and diffuse ISM using a two-component screen model, following the prescription of Charlot and Fall. We compare models in which the dust optical depth is constant to models where it depends on gas metallicity, gas fraction and orientation. The colours of EAGLE galaxies for the more sophisticated models are in broad agreement with those of observed galaxies. In particular, EAGLE produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, with approximately the correct fraction of galaxies in each population and with g-r colours within 0.1 magnitudes of those observed. Luminosity functions from UV to NIR wavelengths differ from observations at a level comparable to systematic shifts resulting from a choice between Petrosian and Kron photometric apertures. Despite the generally good agreement there are clear discrepancies with observations. The blue cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, consistent with the modest underestimate of the passive fraction in massive EAGLE galaxies. There is also a moderate excess of bright blue galaxies compared to observations. The overall level of agreement with the observed colour distribution suggests that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that are comparable to those of observed galaxies.

We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

We performed a hydrodynamical cosmological simulation of the formation of a Milky Way-like galaxy in a warm dark matter (WDM) cosmology. Smooth and dense filaments, several co-moving mega parsec long, form generically above z 2 in this model. Atomic line cooling allows gas in the centres of these filaments to cool to the base of the cooling function, resulting in a very striking pattern of extended Lyman-limit systems (LLSs). Observations of the correlation function of LLSs might hence provide useful limits on the nature of the dark matter. We argue that the self-shielding of filaments may lead to a thermal instability resulting in star formation. We implement a sub-grid model for this, and find that filaments rather than haloes dominate star formation until z 6. Reionisation decreases the gas density in filaments, and the more usual star formation in haloes dominates below z 6, although star formation in filaments continues until z=2. Fifteen per cent of the stars of the z=0 galaxy formed in filaments. At higher redshift, these stars give galaxies a stringy appearance, which, if observed, might be a strong indication that the dark matter is warm.

We provide an analytical description of the line broadening of HI absorbers in the Lyman-alpha forest resulting from Doppler broadening and Jeans smoothing. We demonstrate that our relation captures the dependence of the line-width on column density for narrow lines in z~3 mock spectra remarkably well. Broad lines at a given column density arise when the underlying density structure is more complex, and such clustering is not captured by our model. Our understanding of the line broadening opens the way to a new method to characterise the thermal state of the intergalactic medium and to determine the sizes of the absorbing structures.

We seek to understand the relationship between galaxy properties and their local environment, which calls for a proper formulation of the notion of environment. We analyse the GIMIC suite of cosmological hydrodynamical simulations within the framework of the cosmic web as formulated by Hoffman et al., focusing on properties of simulated DM haloes and luminous galaxies with respect to voids, sheets, filaments and knots - the four elements of the cosmic web. We find that the mass functions of haloes depend on environment, which drives other environmental dependence of galaxy formation. The web shapes the halo mass function, and through the strong dependence of the galaxy properties on the mass of their host haloes, it also shapes the galaxy-(web) environment dependence.

The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm.

The abundance evolution of galaxies depends critically on the balance between the mixing of metals in their interstellar medium, the inflow of new gas and the outflow of enriched gas. We study these processes in gas columns perpendicular to a galactic disk using sub-parsec resolution simulations that track stellar ejecta with the Flash code. We model a simplified interstellar medium stirred and enriched by supernovae and their progenitors. We vary the density distribution of the gas column and integrate our results over an exponential disk to predict wind and ISM enrichment properties for disk galaxies. We find that winds from more massive galaxies are hotter and more highly enriched, in stark contrast to that which is often assumed in galaxy formation models. We use these findings in a simple model of galactic enrichment evolution, in which the metallicity of forming galaxies is the result of accretion of nearly pristine gas and outflow of enriched gas along an equilibrium sequence. We compare these predictions to the observed mass-metallicity relation, and demonstrate how the galaxy's gas fraction is a key controlling parameter. This explains the observed flattening of the mass-metallicity relation at higher stellar masses.

We introduce the Virgo Consortium's 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 resolution 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.

We use observations and simulation to study the relationship between star-forming galaxies and the intergalactic medium at z~3. The observed galaxy sample is based on spectroscopic redshift data from a combination of the VLT LBG Redshift Survey and Keck observations in fields centred on bright z>3 QSOs, whilst the simulation data is taken from GIMIC. In the simulation, we find that the dominant peculiar velocities are in the form of large-scale coherent motions of gas and galaxies. Gravitational infall of galaxies towards one another is also seen. At smaller scales, the peculiar velocities in the simulation over-predict the difference between the simulated real- and z-space galaxy correlation functions. Peculiar velocity pairs separated by <1Mpc/h have a smaller dispersion and explain the z-space correlation function better. The Ly{\alpha} auto- and cross-correlation functions in the GIMIC simulation show infall smaller than implied by previous work. This reduced infall may be due to the galaxy wide outflows implemented in the simulation. The main challenge in comparing these simulated results with the observed correlation functions comes from the presence of velocity errors for the observed LBGs which dominate at ~1Mpc/h scales. When these are taken into account, the observed LBG correlation function is well matched by a simulated $M_*>10^9M_\odot$ galaxy sample. The simulated cross-correlation shows similar neutral gas densities around galaxies as are seen in the observations. The simulated and observed Ly{\alpha} z-space autocorrelation functions agree well with each other. Our overall conclusion is that gas and galaxy peculiar velocities are towards the low end of expectation. Finally, little direct evidence is seen in either simulation or observations for high transmission near galaxies due to feedback. (Abridged)

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.

Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmological N-body simulations of galactic haloes of mass similar to the Milky Way's. On scales larger than the free-streaming cutoff, the initial conditions have the same power spectrum and phases as one of the cold dark matter (CDM) haloes previously simulated by Springel et al as part of the Virgo consortium Aquarius project. We have simulated four haloes with WDM particle masses in the range 1.4-2.3keV and, for one case, we have carried out further simulations at varying resolution. N-body simulations in which the power spectrum cutoff is resolved are known to undergo artificial fragmentation in filaments producing spurious clumps which, for small masses (<10^7Msun in our case) outnumber genuine haloes. We have developed a robust algorithm to identify these spurious objects and remove them from our halo catalogues. We find that the WDM subhalo mass function is suppressed by well over an order magnitude relative to the CDM case for masses <10^9Msun. Requiring that there should be at least as many subhaloes as there are observed satellites in the Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density distributions that are well described by NFW or Einasto profiles. Their central densities are lower for lower WDM particle masses and none of the models we have considered suffer from the "too big to fail" problem recently highlighted by Boylan-Kolchin et al.

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

We present HST/COS observations of highly ionized absorption lines associated with a radio-loud QSO at $z=1.1319$. The absorption system has multiple velocity components, tracing gas that is largely outflowing from the QSO at velocities of a few 100 km s$^{-1}$. There is an unprecedented range in ionization, with detections of HI, NIII, NIV, NV, OIV, OIV*, OV, OVI, NeVIII, MgX, SV and ArVIII. We estimate the total hydrogen number density from the column density ratio N(OIV*)/N(OIV) to be $\log(n_{\textrm{H}}/\textrm{cm}^3)\sim 3$. Assuming photoionization equilibrium, we derive a distance to the absorbing complex of $2.3<R<6.0$ kpc from the centre of the QSO. A range in ionization parameter, covering $\sim 2$ orders of magnitude, suggest absorption path lengths in the range $10^{-4.5}<l_{\textrm{abs}}<1$ pc. In addition, the absorbing gas only partially covers the background emission from the QSO continuum, which suggests clouds with transverse sizes $l_{\textrm{trans}}<10^{-2.5}$ pc. Widely differing absorption path lengths, combined with covering fractions less than unity across all ions pose a challenge to models involving simple cloud geometries. These issues may be mitigated by the presence of non-equilibrium effects, together with the possibility of multiple gas temperatures. The dynamics and expected lifetimes of the gas clouds suggest that they do not originate from close to the AGN, but are instead formed close to their observed location. Their inferred distance, outflow velocities and gas densities are broadly consistent with scenarios involving gas entrainment or condensations in winds driven by either supernovae, or the supermassive black hole accretion disc. In the case of the latter, the present data most likely does not trace the bulk of the outflow by mass, which could instead manifest itself as an accompanying warm absorber, detectable in X-rays.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) cosmological hydrodynamic simulation at z=0 to study the distribution and environmental dependence of neutral hydrogen (HI) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly$\alpha$ absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse HI in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a re-simulation of the average cosmic density GIMIC region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average HI radial profiles by stacking the individual profiles according to both mass and environment. We find high HI column densities at large impact parameters in group environments and markedly lower HI densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

We compute the z = 3 neutral hydrogen column density distribution function f(NHI) for 19 simulations drawn from the OWLS project using a post-processing correction for self-shielding calculated with full radiative transfer of the ionising background radiation. We investigate how different physical processes and parameters affect the abundance of Lyman-limit systems (LLSs) and damped Lyman-alpha absorbers (DLAs) including: i) metal-line cooling; ii) the efficiency of feedback from SNe and AGN; iii) the effective equation of state for the ISM; iv) cosmological parameters; v) the assumed star formation law and; vi) the timing of hydrogen reionization . We find that the normalisation and slope, D = d log10 f /d log10 NHI, of f(NHI) in the LLS regime are robust to changes in these physical processes. Among physically plausible models, f(NHI) varies by less than 0.2 dex and D varies by less than 0.18 for LLSs. This is primarily due to the fact that these uncertain physical processes mostly affect star-forming gas which contributes less than 10% to f(NHI) in the the LLS column density range. At higher column densities, variations in f(NHI) become larger (approximately 0.5 dex at NHI = 10^22 cm^-2 and 1.0 dex at NHI = 10^23 cm^-2) and molecular hydrogen formation also becomes important. Many of these changes can be explained in the context of self-regulated star formation in which the amount of star forming gas in a galaxy will adjust such that outflows driven by feedback balance inflows due to accretion. Data and code to reproduce all figures can be found at the following url: https://bitbucket.org/galtay/hi-cddf-owls-1

We present a new optical spectroscopic survey of 1777 'star-forming' ('SF') and 366 'non-star-forming' ('non-SF') galaxies at redshifts z < 1 (2143 in total), 22 AGN and 423 stars, observed by instruments such as DEIMOS, VIMOS and GMOS, in 3 fields containing 5 quasi-stellar objects (QSOs) with HST UV spectroscopy. We also present a new spectroscopic survey of 165 'strong' (10^14 < NHI < 10^17 cm^-2), and 489 'weak' (10^13 < NHI < 10^14 cm^-2) intervening HI absorption line systems at z < 1 (654 in total), observed in the spectra of 8 QSOs by COS and FOS on the HST. Combining these new data with previously published galaxy catalogs such as VVDS and GDDS, we have gathered a sample of 654 HI absorption systems and 17509 galaxies at transverse scales < 50 Mpc. We present observational results on the HI-galaxy and galaxy-galaxy correlations at transverse scales r < 10 Mpc, and the HI-HI auto-correlation at transverse scales r < 2 Mpc. The two-point correlation functions are measured both along and transverse to the line-of-sight. We constrain the HI-galaxy statistical connection, as a function of both HI column density and galaxy star-forming activity. Our results are consistent with the following conclusions: (1) the bulk of HI systems on Mpc scales have little velocity dispersion (<120 km/s) with respect to the bulk of galaxies; (2) the vast majority of strong HI systems and SF galaxies are distributed in the same locations, together with 75+-15% of non-SF galaxies, all of which typically reside in dark matter haloes of similar masses; (3) 25+-15% of non-SF galaxies reside in galaxy clusters and are not correlated with strong HI systems at scales < 2 Mpc; and (4) 50% of weak HI systems reside within galaxy voids (hence not correlated with galaxies), and are confined in dark matter haloes of masses smaller than those hosting... [abridged]

This paper describes a novel approach to neighbour-finding in Smoothed Particle Hydrodynamics (SPH) simulations with large dynamic range in smoothing length. This approach is based on hierarchical cell decompositions, sorted interactions, and a task-based formulation. It is shown to be faster than traditional tree-based codes, and to scale better than domain decomposition-based approaches on shared-memory parallel architectures such as multi-cores.

We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on metal-line cooling and feedback from star formation and active galactic nuclei (AGN). We find that if the sub-grid feedback from star formation is implemented kinetically, the feedback is only efficient if the initial wind velocity exceeds a critical value. This critical velocity increases with galaxy mass and also if metal-line cooling is included. This suggests that radiative losses quench the winds if their initial velocity is too low. If the feedback is efficient, then the star formation rate is inversely proportional to the amount of energy injected per unit stellar mass formed (which is proportional to the initial mass loading for a fixed wind velocity). This can be understood if the star formation is self-regulating, i.e. if the star formation rate (and thus the gas fraction) increase until the outflow rate balances the inflow rate. Feedback from AGN is efficient at high masses, while increasing the initial wind velocity with gas pressure or halo mass allows one to generate galaxy-wide outflows at all masses. Matching the observed galaxy mass function requires efficient feedback. In particular, the predicted faint-end slope is too steep unless we resort to highly mass loaded winds for low-mass objects. Such efficient feedback from low-mass galaxies (M_* << 10^10 Msun) also reduces the discrepancy with the observed specific star formation rates, which are higher than predicted unless the feedback transitions from highly efficient to inefficient just below the observed stellar mass range.

We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on the assumed star formation law, the equation of state imposed on the unresolved interstellar medium, the stellar initial mass function, the reionization history, and the assumed cosmology. This work complements that of Paper I, where we studied the effects of varying models for galactic winds driven by star formation and AGN. The normalisation of the matter power spectrum strongly affects the galaxy mass function, but has a relatively small effect on the physical properties of galaxies residing in haloes of a fixed mass. Reionization suppresses the stellar masses and gas fractions of low-mass galaxies, but by z = 2 the results are insensitive to the timing of reionization. The stellar initial mass function mainly determines the physical properties of galaxies through its effect on the efficiency of the feedback, while changes in the recycled mass and metal fractions play a smaller role. If we use a recipe for star formation that reproduces the observed star formation law independently of the assumed equation of state of the unresolved ISM, then the latter is unimportant. The star formation law, i.e. the gas consumption time scale as a function of surface density, determines the mass of dense, star-forming gas in galaxies, but affects neither the star formation rate nor the stellar mass. This can be understood in terms of self-regulation: the gas fraction adjusts until the outflow rate balances the inflow rate.

We describe URCHIN, a reverse ray tracing radiative transfer scheme optimised to model self-shielding from the post-reionisation ultraviolet (UV) background in cosmological simulations. The reverse ray tracing strategy provides several benefits over forward ray tracing codes including: (1) the preservation of adaptive density field resolution (2) completely uniform sampling of gas elements by rays; (3) the preservation of galilean invariance; (4) the ability to sample the UV background spectrum with hundreds of frequency bins; and (5) exact preservation of the input UV background spectrum and amplitude in optically thin gas. The implementation described here focuses on Smoothed Particle Hydrodynamics (SPH). However, the method can be applied to any density field representation in which resolution elements admit ray intersection tests and can be associated with optical depths. We characterise the errors in our implementation in stages beginning with comparison to known analytic solutions and ending with a realistic model of the z = 3 cosmological UV background incident onto a suite of spherically symmetric models of gaseous galactic halos.

We present a new parallel implementation of the PINpointing Orbit Crossing-Collapsed HIerarchical Objects (PINOCCHIO) algorithm, a quick tool, based on Lagrangian Perturbation Theory, for the hierarchical build-up of Dark Matter halos in cosmological volumes. To assess its ability to predict halo correlations on large scales, we compare its results with those of an N-body simulation of a 3 Gpc/h box sampled with 2048^3 particles taken from the MICE suite, matching the same seeds for the initial conditions. Thanks to the FFTW libraries and to the relatively simple design, the code shows very good scaling properties. The CPU time required by PINOCCHIO is a tiny fraction (~1/2000) of that required by the MICE simulation. Varying some of PINOCCHIO numerical parameters allows one to produce a universal mass function that lies in the range allowed by published fits, although it underestimates the MICE mass function of FoF halos in the high mass tail. We compare the matter-halo and the halo-halo power spectra with those of the MICE simulation and find that these 2-point statistics are well recovered on large scales. In particular, when catalogs are matched in number density, agreement within ten per cent is achieved for the halo power spectrum. At scales k>0.1 h/Mpc, the inaccuracy of the Zel'dovich approximation in locating halo positions causes an underestimate of the power spectrum that can be modeled as a Gaussian factor with a damping scale of d=3 Mpc/h at z=0, decreasing at higher redshift. Finally, a remarkable match is obtained for the reduced halo bispectrum, showing a good description of nonlinear halo bias. Our results demonstrate the potential of PINOCCHIO as an accurate and flexible tool for generating large ensembles of mock galaxy surveys, with interesting applications for the analysis of large galaxy redshift surveys.

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 one 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.

Models of galaxy formation in a CDM universe predict that massive galaxies are surrounded by a hot, quasi-hydrostatic circumgalactic corona of slowly cooling gas, predominantly accreted from the IGM. This prediction is borne out by the cosmological hydrodynamical simulations of Crain et al., which reproduce scaling relations between the X-ray and optical properties of nearby disc galaxies. Such coronae are metal poor, but observations of the X-ray emitting circumgalactic medium (CGM) of local galaxies typically indicate enrichment to near-solar iron abundance, potentially signalling a shortcoming in galaxy formation models. We show here that, while the hot CGM of galaxies formed in the simulations is metal poor in a mass-weighted sense, its X-ray luminosity-weighted metallicity is often close to solar. This bias arises because the soft X-ray emissivity of a typical 0.1 keV corona is dominated by collisionally-excited metal ions that are synthesised in stars and recycled into the hot CGM. We find that these metals are ejected primarily by stars that form in-situ to the main progenitor of the galaxy, rather than in satellites or external galaxies. The enrichment of the hot CGM therefore proceeds in an inside-out fashion throughout the assembly of the galaxy: metals are transported from the central galaxy by SNe-driven winds and convection over several gigayears, establishing a strong negative radial metallicity gradient. Whilst metal ions synthesised by stars are necessary to produce the X-ray emissivity that enables the hot CGM of isolated galaxies to be detected, the electrons that collisionally excite them are equally important. Since our simulations indicate that the electron density of hot coronae is dominated by the metal-poor gas accreted from the IGM, we infer that the hot CGM observed via X-ray emission is the outcome of both hierarchical accretion and stellar recycling.

We present a catalogue of 2135 galaxy redshifts from the VLT LBG Redshift Survey (VLRS), a spectroscopic survey of z ~ 3 galaxies in wide fields centred on background quasi-stellar objects. We have used deep optical imaging to select galaxies via the Lyman-break technique. Spectroscopy of the Lyman-break galaxies (LBGs) was then made using the Visible Multi-Object Spectrograph (VIMOS), giving a mean redshift of z=2.79. We analyse the clustering properties of the VLRS sample and also of the VLRS sample combined with the smaller area Keck-based survey of Steidel et al. From the semiprojected correlation function, wp({\sigma}) we find that the results are well fit with a single power-law model, with clustering scale lengths of r0=3.46+-0.41 and 3.83+-0.24 Mpc/h, respectively. We note that the corresponding combined {\xi}(r) slope is flatter than for local galaxies at {\gamma} = 1.5-1.6 rather than {\gamma}=1.8. This flat slope is confirmed by the z-space correlation function, {\xi}(s), and in the range 10<s<100 Mpc/h the VLRS shows ~2.5{\sigma} excess over the {\Lambda} cold dark matter. This excess may be consistent with recent evidence for non-Gaussianity in clustering results at z~1. We then analyse the LBG z-space distortions using the 2D correlation function, {\xi}({\sigma}, {\pi}), finding for the combined sample a large-scale infall parameter of $\beta$ = 0.38+-0.19 and a velocity dispersion of 420km/s. Based on our measured {\beta}, we are able to determine the gravitational growth rate, finding a value of f(z = 3)=0.99+-0.50 (or f{\sigma}8 = 0.26+-0.13), which is the highest redshift measurement of the growth rate via galaxy clustering and is consistent with {\Lambda}CDM. Finally, we constrain the mean halo mass for the LBG population, finding that the VLRS and combined sample suggest mean halo masses of log(MDM/Msun) = 11.57+-0.15 and 11.73+-0.07, respectively.

We have performed a series of numerical experiments to investigate how the primordial thermal velocities of fermionic dark matter particles affect the physical and phase space density profiles of the dark matter haloes into which they collect. The initial particle velocities induce central cores in both profiles, which can be understood in the framework of phase space density theory. We find that the maximum coarse-grained phase space density of the simulated haloes (computed in 6 dimensional phase space using the EnBid code) is very close to the theoretical fine-grained upper bound, while the pseudo phase space density, Q ~ {\rho}/{\sigma}^3, overestimates the maximum phase space density by up to an order of magnitude. The density in the inner regions of the simulated haloes is well described by a 'pseudo-isothermal' profile with a core. We have developed a simple model based on this profile which, given the observed surface brightness profile of a galaxy and its central velocity dispersion, accurately predicts its central phase space density. Applying this model to the dwarf spheroidal satellites of the Milky Way yields values close to 0.5 keV for the mass of a hypothetical thermal warm dark matter particle, assuming the satellite haloes have cores produced by warm dark matter free streaming. Such a small value is in conflict with the lower limit of 1.2 keV set by observations of the Lyman-{\alpha} forest. Thus, if the Milky Way dwarf spheroidal satellites have cores, these are likely due to baryonic processes associated with the forming galaxy, perhaps of the kind proposed by Navarro, Eke and Frenk and seen in recent simulations of galaxy formation in the cold dark matter model.

We present adaptive optics assisted integral field spectroscopy of nine Halpha-selected galaxies at z=0.84--2.23 selected from the HiZELS narrow-band survey. Our observations map the star-formation and kinematics of these representative star-forming galaxies on ~kpc-scales. We demonstrate that within the ISM of these galaxies, the velocity dispersion of the star-forming gas (\sigma) follows a scaling relation \sigma\propto\Sigma_SFR^(1/n)+constant (where \Sigma_SFR is the star formation surface density and the constant includes the stellar surface density). Assuming the disks are marginally stable (Toomre Q=1), we show that this follows from the Kennicutt-Schmidt relation (\Sigma_SFR = A\Sigma_gas^n), and we use the data to derive best fit parameters of n=1.34+/-0.15 and A=3.4_(-1.6)^(+2.5)x10^(-4)Mo/yr/kpc^2, consistent with the local relation and implying cold molecular gas masses of M_gas=10^(9-10)Mo and molecular gas fractions M_gas/(M_gas+Mstars)=0.3+/-0.1, with a range of 10-75%. These values confirm the high gas fractions for high-redshift star-forming galaxies, independent of CO-H_2 conversion factor. We also identify eleven ~kpc-scale star-forming regions (clumps) within our sample and show that their sizes are comparable to the wavelength of the fastest growing unstable mode. The luminosities and velocity dispersions of these clumps follow the same scaling relations as local HII regions, although their star formation densities are a factor 15+/-5x higher than typically found locally. We discuss how the clump properties are related to the disk, and show that their high masses and luminosities are a consequence of the high disk surface density.

Feedback from supernovae is an essential aspect of galaxy formation. In order to improve subgrid models of feedback we perform a series of numerical experiments to investigate how supernova explosions power galactic winds. We use the Flash hydrodynamic code to model a simplified ISM, including gravity, hydrodynamics, radiative cooling above 10,000 K, and star formation that reproduces the Kennicutt-Schmidt relation. By simulating a small patch of the ISM in a tall box perpendicular to the disk, we obtain sub-parsec resolution allowing us to resolve individual supernova events and we investigate how the wind properties depend on those of the ISM and the galaxy. We find that outflows are more efficient in disks with lower surface densities or gas fractions. A simple model in which the warm cloudy medium is the barrier that limits the expansion of blast waves reproduces the scaling of outflow properties with disk parameters at high star formation rates. The scaling we find sets the investigation of galaxy winds on a new footing, providing a physically motivated sub-grid description of winds that can be implemented in cosmological hydrodynamic simulations and phenomenological models. [Abridged]

We compare the observed probability distribution function of the transmission in the \HI\ Lyman-alpha forest, measured from the UVES 'Large Programme' sample at redshifts z=[2,2.5,3], to results from the GIMIC cosmological simulations. Our measured values for the mean transmission and its PDF are in good agreement with published results. Errors on statistics measured from high-resolution data are typically estimated using bootstrap or jack-knife resampling techniques after splitting the spectra into chunks. We demonstrate that these methods tend to underestimate the sample variance unless the chunk size is much larger than is commonly the case. We therefore estimate the sample variance from the simulations. We conclude that observed and simulated transmission statistics are in good agreement, in particular, we do not require the temperature-density relation to be 'inverted'.

We present adaptive optics assisted, spatially resolved spectroscopy of a sample of nine H-alpha-selected galaxies at z=0.84--2.23 drawn from the HiZELS narrow-band survey. These galaxies have star-formation rates of 1-27Mo/yr and are therefore representative of the typical high-redshift star-forming population. Our ~kpc-scale resolution observations show that approximately half of the sample have dynamics suggesting that the ionised gas is in large, rotating disks. We model their velocity fields to infer the inclination-corrected, asymptotic rotational velocities. We use the absolute B-band magnitudes and stellar masses to investigate the evolution of the B-band and stellar mass Tully-Fisher relationships. By combining our sample with a number of similar measurements from the literature, we show that, at fixed circular velocity, the stellar mass of star-forming galaxies has increased by a factor 2.5 between z=2 and z=0, whilst the rest-frame B-band luminosity has decreased by a factor ~6 over the same period. Together, these demonstrate a change in mass-to-light ratio in the B-band of Delta(M/L_B)/(M/L_B)_(z=0) \sim 3.5 between z=1.5 and z=0, with most of the evolution occurring below z=1. We also use the spatial variation of [NII]/Halpha to show that the metallicity of the ionised gas in these galaxies declines monotonically with galacto-centric radius, with an average Delta(log O/H)/DeltaR=-0.027+/-0.005dex/kpc. This gradient is consistent with predictions for high-redshift disk galaxies from cosmologically based hydrodynamic simulations.

We examine the rotation rates, sizes, and star formation (SF) efficiencies of a representative population of simulated disc galaxies extracted from the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamic simulations. These simulations include efficient, but energetically feasible supernova feedback, but have not been tuned in any way to produce 'realistic' disc galaxies. Yet, they generate a large number of discs, without requiring extremely high resolution. Over the wide galaxy stellar mass range, 9.0 < log10[Mstar (Msun)] < 10.5, the simulations reproduce the observed Tully-Fisher relation, the rotation curves of disc galaxies in bins of stellar mass, the mass-size relation of disc galaxies, the optical rotation to virial circular velocity ratio (Vopt/Vvir), and the SF efficiencies of disc galaxies as inferred from stacked weak lensing and stacked satellite kinematics observations. They also reproduce the specific star formation rates of ~L* galaxies but predict too low levels of star formation for low-mass galaxies, which is plausibly due to the finite resolution of the simulations. At higher stellar masses, log10[Mstar (\Msun)] > 10.6, the simulated galaxies are too concentrated and have too high SF efficiencies. We conjecture that this shortcoming reflects the neglect of feedback from accreting supermassive black holes in these simulations. We conclude that it is possible to generate a representative population of disc galaxies that reproduces many of the observed trends of local disc galaxies using standard numerical hydrodynamic techniques and a plausible implementation of the "subgrid" astrophysical processes thought to be relevant to galaxy formation.

In the simplest scenario, disk galaxies form predominantly in halos with high angular momentum and quiet recent assembly history, whereas spheroids are the slowly-rotating remnants of repeated merging events. We explore these assumptions using one hundred systems with halo masses similar to that of the Milky Way, identified in a series of cosmological gasdynamical simulations GIMIC. At z=0, the simulated galaxies exhibit a wide variety of morphologies, from dispersion-dominated spheroids to pure disk galaxies. Surprisingly, these morphological features are very poorly correlated with their halo properties: disks form in halos with high and low net spin, and mergers play a negligible role in the formation of spheroid stars, most of which form in-situ. More important to morphology is the coherent alignment of the angular momentum of baryons that accrete over time to form a galaxy. Spheroids tend to form when the spin of newly-accreted gas is misaligned with that of the extant galaxy, leading to the episodic formation of stars with different kinematics that cancel out the net rotation of the system. Disks, on the other hand, form out of gas that flows in with similar angular momentum to that of earlier-accreted material. Gas accretion from a hot corona thus favours disk formation, whereas gas that flows "cold", often along separate, misaligned filaments, favours the formation of spheroids. In this scenario, most spheroids consist of superpositions of stellar components with distinct kinematics, age, and metallicity, an arrangement that might survive to the present day given the paucity of major mergers. Since angular momentum is acquired largely at turnaround, morphology is imprinted early by the interplay of the tidal field and the shape of the material destined to form the galaxy.

We investigate the physical state of HI absorbing gas at low redshift (z = 0.25) using a subset of cosmological, hydrodynamic simulations from the OWLS project, focusing in particular on broad (b_HI > 40 km/s) HI Lyman-Alpha absorbers (BLAs), which are believed to originate in shock-heated gas in the warm-hot intergalactic medium (WHIM). Our fiducial model, which includes radiative cooling by heavy elements and feedback by supernovae and active galactic nuclei, predicts that by z = 0.25 nearly 60 per cent of the gas mass ends up at densities and temperatures characteristic of the WHIM and we find that half of this fraction is due to outflows. The standard HI observables (distribution of HI column densities N_HI, distribution of Doppler parameters b_HI, b_HI - N_HI correlation) and the BLA line number density predicted by our simulations are in remarkably good agreement with observations. BLAs arise in gas that is hotter, more highly ionised and more enriched than the gas giving rise to typical Lyman-Alpha forest absorbers. The majority of the BLAs arise in warm-hot (log (T/K) ~ 5) gas at low (log Delta < 1.5) overdensities. On average, thermal broadening accounts for at least 60 per cent of the BLA line width, which in turn can be used as a rough indicator of the thermal state of the gas. Detectable BLAs account for only a small fraction of the true baryon content of the WHIM at low redshift. In order to detect the bulk of the mass in this gas phase, a sensitivity at least one order of magnitude better than achieved by current ultraviolet spectrographs is required. We argue that BLAs mostly trace gas that has been shock-heated and enriched by outflows and that they therefore provide an important window on a poorly understood feedback process.

Ram pressure stripping of the hot gas that surrounds normal galaxies as they fall into groups and clusters (also referred to as `strangulation' or `starvation') is generally thought to shut down star formation on a time scale of a few Gyr. However, it has recently been suggested, on the basis of X-ray-optical scaling relations of galaxies in the field and the group/cluster environment, that confinement pressure by the intra-cluster medium can actually lead to an increase in the mass of hot gas surrounding these galaxies. We investigate the competition between pressure confinement and ram pressure stripping for satellite galaxies in orbit about galaxy groups and clusters using simple analytic models and detailed cosmological hydrodynamic simulations. It is found that, independent of host mass, ram pressure is generally dominant over confinement pressure --- only ~16 % of galaxies find themselves in the reverse situation. Furthermore, these galaxies have, on average, less hot gas than ram-pressure dominated ones, contrary to simple expectations. This is explained by the fact that the small number of galaxies which are confinement dominated are typically at first or second apocentre and have therefore already been maximally affected by ram pressure stripping around first pericentre. Our results are shown to be insensitive to host halo mass; we argue that the same is true for uncertain sub-grid processes, such as feedback.

The formation of galaxies is regulated by a balance between the supply of gas and the rate at which it is ejected. Traditional explanations of gas ejection equate the energy required to escape the galaxy or host halo to an estimate for the energy yield from supernovae. This yield is usually assumed to be a constant fraction of the total available from the supernova, or is derived from the assumption of a consistent momentum yield. By applying these ideas in the context of a cold dark matter cosmogony, we derive a 1st-order analytic connection between these working assumptions and the expected relationship between baryon content and galaxy circular velocity, and find that these quick predictions straddle recent observational estimates. To examine the premises behind these theories in more detail, we then explore their applicability to a set of gasdynamical simulations of idealised galaxies. We show that different premises dominate to differing degrees in the simulated outflow, depending on the mass of the system and the resolution with which it is simulated. Using this study to anticipate the emergent behaviour at arbitrarily high resolution, we motivate more comprehensive analytic model which allows for the range of velocities with which the gas may exit the system, and incorporates both momentum and energy-based constraints on the outflow. Using a trial exit velocity distribution, this is shown to be compatible with the observed baryon fractions in intermediate-mass systems, but implies that current estimates for low-mass systems can not be solely accounted for by supernova winds under commonly-held assumptions.

We investigate the properties of the HI Ly-a absorption systems (Ly-a forest) within and around galaxy voids at z<0.1. We find a significant excess (>99% c.l.) of Ly-a systems at the edges of galaxy voids with respect to a random distribution, on ~5 h^{-1} Mpc scales. We find no significant difference in the number of systems inside voids with respect to the random expectation. We report differences between both column density (N_{HI}) and Doppler parameter (b_{HI}) distributions of Ly-a systems found inside and at the edge of galaxy voids at the >98% and >90% c.l. respectively. Low density environments (voids) have smaller values for both N_{HI} and b_{HI} than higher density ones (edges of voids). These trends are theoretically expected and also found in GIMIC, a state-of-the-art hydrodynamical simulation. Our findings are consistent with a scenario of at least three types of Ly-alpha: (1) containing embedded galaxies and so directly correlated with galaxies (referred as `halo-like'), (2) correlated with galaxies only because they lie in the same over-dense LSS, and (3) associated with under-dense LSS with a very low auto-correlation amplitude (~ random) that are not correlated with luminous galaxies. We argue the latter arise in structures still growing linearly from the primordial density fluctuations inside galaxy voids that have not formed galaxies because of their low densities. We estimate that these under-dense LSS absorbers account for 25-30% +- 6% of the current Ly-a population (N_{HI} > 10^{12.5} cm^{-2}) while the other two types account for the remaining 70-75% +- 12%. Assuming that only N_{HI} > 10^{14} cm^{-2} systems have embedded galaxies nearby, we have estimated the contribution of the `halo-like' Ly-a to be ~12-15% +- 4% and consequently ~55-60% +- 13% of the Ly-a systems to be associated with the over-dense LSS.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.

Simulations predict that galaxies grow primarily through the accretion of gas that has not gone through an accretion shock near the virial radius and that this cold gas flows towards the central galaxy along dense filaments and streams. There is, however, little observational evidence for the existence of these cold flows. We use a large, cosmological, hydrodynamical simulation that has been post-processed with radiative transfer to study the contribution of cold flows to the observed z=3 column density distribution of neutral hydrogen, which our simulation reproduces. We find that nearly all of the HI absorption arises in gas that has remained colder than 10^5.5 K, at least while it was extragalactic. In addition, the majority of the HI is rapidly falling towards a nearby galaxy, with non-negligible contributions from outflowing and static gas. Above a column density of N_HI = 10^17 cm^-2, most of the absorbers reside inside haloes, but the interstellar medium only dominates for N_HI > 10^21 cm^-2. Haloes with total mass below 10^10 Msun dominate the absorption for 10^17<N_HI < 10^21 cm^-2, but the average halo mass increases sharply for higher column densities. Although very little of the HI in absorbers with N_HI <~ 10^20 cm^-2 resides inside galaxies, systems with N_HI > 10^17 cm^-2 are closely related to star formation: most of their HI either will become part of the interstellar medium before z=2 or has been ejected from a galaxy at z>3. Cold accretion flows are critical for the success of our simulation in reproducing the observed rate of incidence of damped Lyman-alpha and particularly that of Lyman limit systems. We therefore conclude that cold accretion flows exist and have already been detected in the form of high column density HI absorbers.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r < ~30 kpc. In the present study we show that this component was primarily formed in a proto-disc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, that is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axis ratio of ~0.6 at radii of < ~30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ vs. accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a 'dual stellar halo'.

We investigate the column density distribution function of neutral hydrogen at redshift z = 3 using a cosmological simulation of galaxy formation from the OverWhelmingly Large Simulations (OWLS) project. The base simulation includes gravity, hydrodynamics, star formation, supernovae feedback, stellar winds, chemodynamics, and element-by-element cooling in the presence of a uniform UV background. Self-shielding and formation of molecular hydrogen are treated in post-processing, without introducing any free parameters, using an accurate reverse ray-tracing algorithm and an empirical relation between gas pressure and molecular mass fraction. The simulation reproduces the observed z = 3 abundance of Ly-A forest, Lyman Limit and Damped Ly-A HI absorption systems probed by quasar sight lines over ten orders of magnitude in column density. Self-shielding flattens the column density distribution for NHI > 10^18 cm-2, while the conversion to fully neutral gas and conversion of HI to H2 steepen it around column densities of NHI = 10^20.3 cm-2 and NHI = 10^21.5 cm-2, respectively.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the formation of stellar spheroids of Milky Way-mass disc galaxies. The simulations contain accurate treatments of metal-dependent radiative cooling, star formation, supernova feedback, and chemodynamics, and the large volumes that have been simulated yield an unprecedentedly large sample of ~400 simulated L_* disc galaxies. The simulated galaxies are surrounded by low-mass, low-surface brightness stellar haloes that extend out to ~100 kpc and beyond. The diffuse stellar distributions bear a remarkable resemblance to those observed around the Milky Way, M31 and other nearby galaxies, in terms of mass density, surface brightness, and metallicity profiles. We show that in situ star formation typically dominates the stellar spheroids by mass at radii of r < 30 kpc, whereas accretion of stars dominates at larger radii and this change in origin induces a change in slope of the surface brightness and metallicity profiles, which is also present in the observational data. The system-to-system scatter in the in situ mass fractions of the spheroid, however, is large and spans over a factor of 4. Consequently, there is a large degree of scatter in the shape and normalisation of the spheroid density profile within r < 30 kpc (e.g., when fit by a spherical powerlaw profile the indices range from -2.6 to -3.4). We show that the in situ mass fraction of the spheroid is linked to the formation epoch of the system. Dynamically older systems have, on average, larger contributions from in situ star formation, although there is significant system-to-system scatter in this relationship. Thus, in situ star formation likely represents the solution to the longstanding failure of pure accretion-based models to reproduce the observed properties of the inner spheroid.

The intergalactic medium was not completely reionized until approximately a billion years after the Big Bang, as revealed by observations of quasars with redshifts of less than 6.5. It has been difficult to probe to higher redshifts, however, because quasars have historically been identified in optical surveys, which are insensitive to sources at redshifts exceeding 6.5. Here we report observations of a quasar (ULAS J112001.48+064124.3) at a redshift of 7.085, which is 0.77 billion years after the Big Bang. ULAS J1120+0461 had a luminosity of 6.3x10^13 L_Sun and hosted a black hole with a mass of 2x10^9 M_Sun (where L_Sun and M_Sun are the luminosity and mass of the Sun). The measured radius of the ionized near zone around ULAS J1120+0641 is 1.9 megaparsecs, a factor of three smaller than typical for quasars at redshifts between 6.0 and 6.4. The near zone transmission profile is consistent with a Ly alpha damping wing, suggesting that the neutral fraction of the intergalactic medium in front of ULAS J1120+0641 exceeded 0.1.

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.

We present a study of cooling in radiative shocks simulated with smoothed particle hydrodynamics (SPH) and adaptive mesh refinement codes. We obtain a similarity solution for a shock-tube problem in the presence of radiative cooling, and test how well the solution is reproduced in Gadget and Flash. Shock broadening governed by the details of the numerical scheme (artificial viscosity or Riemann solvers) leads to potentially significant overcooling in both codes. We interpret our findings in terms of a resolution criterion, and apply it to realistic simulations of cosmological accretion shocks onto galaxy haloes, cold accretion and thermal feedback from supernovae or active galactic nuclei. To avoid numerical overcooling of accretion shocks onto haloes that should develop a hot corona requires a particle or cell mass resolution of 10^6 M_sun, which is within reach of current state-of-the-art simulations. At this mass resolution, thermal feedback in the interstellar medium of a galaxy requires temperatures of supernova or AGN driven bubbles to be in excess of 10^7 K at densities of n_H=1.0 cm^-3, in order to avoid spurious suppression of the feedback by numerical overcooling.

We present the initial imaging and spectroscopic data acquired as part of the VLT VIMOS Lyman-break galaxy Survey. UBR (or UBVI) imaging covers five 36'x36' fields centred on bright z>3 QSOs, allowing ~21,000 2<z<3.5 galaxy candidates to be selected using the Lyman-break technique. We performed spectroscopic follow-up using VIMOS, measuring redshifts for 1020 z>2 LBGs and 10 z>2 QSOs from a total of 19 VIMOS pointings. From the galaxy spectra, we observe a 625+-510 km/s velocity offset between the ISM absorption and Ly-alpha emission line redshifts. Using the photometric and spectroscopic catalogues, we have analysed the galaxy clustering at z~3. In the photometric case, the angular correlation function, w(theta), is well fit by a double power-law with clustering scale-length, r_0 = 3.19+0.32-0.54 Mpc/h for r < 1 Mpc/h and r_0 = 4.59+0.31-0.33 Mpc/h at larger scales. Using the redshift sample we estimate the semi-projected correlation function, w_p(sigma) and find r_0 = 3.67+0.23-0.24 Mpc/h for the VLT sample and r_0 = 3.98+0.14-0.15 Mpc/h for a combined VLT+Keck sample. From the z-space correlation functions and assuming the above xi(r) models, we find that the combined VLT and Keck surveys require a galaxy velocity dispersion, <w_z^2>^1/2 ~ 700 km/s, higher than the ~400 km/s found by previous authors. We also measure a value for the gravitational growth rate parameter of beta(z=3) = 0.48+-0.17, implying a low value for the bias of b = 2.06+1.1-0.5. This value is consistent with the galaxy clustering amplitude which gives b = 2.22+-0.16, assuming the standard cosmology, implying that the evolution of the gravitational growth rate is also consistent with Einstein gravity. We have compared our LBG clustering amplitudes with lower redshift measurements and find that the clustering strength is not inconsistent with that of low-redshift L* spirals for simple 'long-lived' galaxy models.

[Abridged] We investigate how a range of physical processes affect the cosmic metal distribution using a suite of cosmological, hydrodynamical simulations. Focusing on z = 0 and 2, we study the metallicities and metal mass fractions for stars as well as for the ISM, and several more diffuse gas phases. We vary the cooling rates, star formation law, structure of the ISM, galactic winds, feedback from AGN, reionization history, stellar IMF, and cosmology. In all models stars and the warm-hot IGM (WHIM) constitute the dominant repository of metals, while for z > 2 the ISM is also important. In models with galactic winds, predictions for the metallicities of the various phases vary at the factor of two level and are broadly consistent with observations. The exception is the cold-warm IGM, whose metallicity varies at the order of magnitude level if the prescription for galactic winds is varied, even for a fixed wind energy per unit stellar mass formed, and falls far below the observed values if winds are not included. At the other extreme, the metallicity of the intracluster medium (ICM) is largely insensitive to the presence of galactic winds, indicating that its enrichment is regulated by other processes. The mean metallicities of stars (~ Z_sun), the ICM (~ 0.1 Z_sun), and the WHIM (~ 0.1 Z_sun) evolve only slowly, while those of the cold halo gas and the IGM increase by more than an order of magnitude from z = 5 to 0. Higher velocity outflows are more efficient at transporting metals to low densities, but actually predict lower metallicities for the cold-warm IGM since the winds shock-heat the gas to high temperatures, thereby increasing the fraction of the metals residing in, but not the metallicity of, the WHIM. Besides galactic winds driven by feedback from star formation, the metal distribution is most sensitive to the inclusion of metal-line cooling and feedback from AGN.

We analyse the stellar and hot gas content of 18 nearby, low-mass galaxy clusters, detected in redshift space and selected to have a dynamical mass 3E14<M/Msun<6E14, as measured from the 2dF Galaxy Redshift Survey. We combine X-ray measurements from both Chandra and XMM with ground-based near-infrared observations from CTIO, AAT and CFHT to compare the mass in hot gas and stars to the dynamical mass and state of the clusters. Only 13 of the clusters are detected in X-ray emission, and for these systems we find that a range of 7-20 per cent of their baryonic mass, and <3 per cent of their dynamical mass, is detected in starlight, similar to what is observed in more massive clusters. In contrast, the five undetected clusters are underluminous in X-ray emission, by up to a factor 10, given their stellar mass. Although the velocity distribution of cluster members in these systems is indistinguishable from a Gaussian, all show subtle signs of being unrelaxed: either they lack a central, dominant galaxy, or the bright galaxy distribution is less concentrated and/or more elongated than the rest of the sample. Thus we conclude that low-mass clusters and groups selected from the velocity distribution of their galaxies exhibit a dichotomy in their hot gas properties. Either they are detected in X-ray, in which case they generally lie on the usual scaling relations, or they are completely undetected in X-ray emission. The non-detections may be partly related to the apparently young dynamical state of the clusters, but it remains a distinct possibility that some of these systems are exceptionally devoid of hot emitting gas as the result of its expulsion or rarefaction.

We investigate the origin and physical properties of OVI absorbers at low redshift (z = 0.25) using a subset of cosmological, hydrodynamical simulations from the OverWhelmingly Large Simulations (OWLS) project. Intervening OVI absorbers are believed to trace shock-heated gas in the Warm-Hot Intergalactic Medium (WHIM) and may thus play a key role in the search for the missing baryons in the present-day Universe. When compared to observations, the predicted distributions of the different OVI line parameters (column density, Doppler parameter, rest equivalent width) from our simulations exhibit a lack of strong OVI absorbers. This suggests that physical processes on sub-grid scales (e.g. turbulence) may strongly influence the observed properties of OVI systems. We find that the intervening OVI absorption arises mainly in highly metal-enriched (0.1 << Z/Z_sun < 1) gas at typical overdensities of 1 << rho/<rho> < 100. One third of the OVI absorbers in our simulation are found to trace gas at temperatures T < 10^5 K, while the rest arises in gas at higher temperatures around T =10^5.3 K. The OVI resides in a similar region of (rho,T)-space as much of the shock-heated baryonic matter, but the vast majority of this gas has a lower metal content and does not give rise to detectable OVI absorption As a consequence of the patchy metal distribution, OVI absorbers in our simulations trace only a very small fraction of the cosmic baryons (<2 percent) and the cosmic metals. Instead, these systems presumably trace previously shock-heated, metal-rich material from galactic winds that is now cooling. The common approach of comparing OVI and HI column densities to estimate the physical conditions in intervening absorbers from QSO observations may be misleading, as most of the HI (and most of the gas mass) is not physically connected with the high-metallicity patches that give rise to the OVI absorption.

An ionization front expanding into a neutral medium can be slowed-down significantly by recombinations. In cosmological numerical simulations the recombination rate is often computed using a 'clumping factor', that takes into account that not all scales in the simulated density field are resolved. Here we demonstrate that using a single value of the clumping factor significantly overestimates the recombination rate, and how a local estimate of the clumping factor is both easy to compute, and gives significantly better numerical convergence. We argue that this lower value of the recombination rate is more relevant during the reionization process and hence that the importance of recombinations during reionization has been overestimated.

The existence of hot, X-ray luminous gaseous coronae surrounding present day L* galaxies is a generic prediction of galaxy formation theory in the cold dark matter cosmogony. While extended X-ray emission has been known to exist around elliptical galaxies for a long time, diffuse extra-planar emission has only recently been detected around disc galaxies. We compile samples of elliptical and disc galaxies that have Chandra and XMM-Newton measurements, and compare the scaling of the coronal X-ray luminosity (L_X) with both the K-band luminosity (L_K) and the coronal X-ray temperature (T_X). The X-ray flux measurements are corrected for non-thermal point source contamination by spatial excision and spectral subtraction for resolved and unresolved sources respectively. We find that the properties of the extended X-ray emission from galaxies of different morphological types are similar: for both elliptical and disc galaxies, the L_X - L_K and L_X - T_X relations have similar slope, normalisation and scatter. The observed universality of coronal X-ray properties suggests that the bulk of this emission originates from gas that has been accreted, shock-heated and compressed during the assembly of the galaxy and that outflows triggered by stellar processes make only a minor contribution to the X-ray emission. This reservoir of cooling gas is a potential source of morphological transformation; it provides a fresh supply of material for discs to grow around galaxies of all morphological types.

The Durham GALFORM semi-analytical galaxy formation model has been shown to reproduce the observed rest-frame 1500\AA\ luminosity function of galaxies well over the whole redshift range z=5-10. We show that in this model, this galaxy population also emits enough ionizing photons to reionize the Universe by redshift z=10, assuming a modest escape fraction of 20 per cent. The bulk of the ionizing photons is produced in faint galaxies during starbursts triggered by galaxy mergers. The bursts introduce a dispersion up to ~ 5 dex in galaxy ionizing luminosity at a given halo mass. Almost 90 per cent of the ionizing photons emitted at z=10 are from galaxies below the current observational detection limit at that redshift. Photo-ionization suppression of star formation in these galaxies is unlikely to affect this conclusion significantly, because the gas that fuels the starbursts has already cooled out of their host halos. The galaxies that dominate the ionizing emissivity at z=10 are faint, with M_{1500, AB} ~ -16, have low star formation rates, \dot{M_{*}} ~ 0.06 h^{-1} M_sun yr^{-1}, and reside in halos of mass M ~ 10^9 h^{-1} M_sun.

We use a suite of cosmological, hydrodynamical simulations to investigate the chemical enrichment history of the Universe. Specifically, we trace the origin of the metals back in time to investigate when various gas phases were enriched and by what halo masses. We find that the age of the metals decreases strongly with the density of the gas in which they end up. At least half of the metals that reside in the diffuse intergalactic medium (IGM) at redshift zero (two) were ejected from galaxies above redshift two (three). The mass of the haloes that last contained the metals increases rapidly with the gas density. More than half of the mass in intergalactic metals was ejected by haloes with total masses less than 1e11 solar masses and stellar masses less than 1e9 solar masses. The range of halo masses that contributes to the enrichment is wider for the hotter part of the IGM. By combining the `when' and `by what' aspects of the enrichment history, we show that metals residing in lower density gas were typically ejected earlier and by lower mass haloes.

We have measured redshifts for 243 z ~3 quasars in nine VLT VIMOS LBG redshift survey areas, each of which is centred on a known bright quasar. Using spectra of these quasars, we measure the cross-correlation between neutral hydrogen gas causing the Lya forest and 1020 Lyman-break galaxies at z ~3. We find an increase in neutral hydrogen absorption within 5 h^-1 Mpc of a galaxy in agreement with the results of Adelberger et al. (2003, 2005). The Lya-LBG cross-correlation can be described by a power-law on scales larger than 3 h^-1 Mpc. When galaxy velocity dispersions are taken into account our results at smaller scales (<2 h^-1 Mpc) are also in good agreement with the results of Adelberger et al. (2005). There is little immediate indication of a region with a transmission spike above the mean IGM value which might indicate the presence of star-formation feedback. To measure the galaxy velocity dispersions, which include both intrinsic LBG velocity dispersion and redshift errors, we have used the LBG-LBG redshift space distortion measurements of Bielby et al. (2010). We find that the redshift-space transmission spike implied in the results of Adelberger et al. (2003) is too narrow to be physical in the presence of the likely LBG velocity dispersion and is likely to be a statistical fluke. Nevertheless, neither our nor previous data can rule out the presence of a narrow, real-space transmission spike, given the evidence of the increased Lya absorption surrounding LBGs which can mask the spike's presence when convolved with a realistic LBG velocity dispersion. Finally, we identify 176 CIV systems in the quasar spectra and find an LBG-CIV correlation strength on scales of 10 h^-1 Mpc consistent with the relation measured at ~Mpc scales.

Approximately half the baryons in the local Universe are thought to reside in the warm-hot intergalactic medium (WHIM). Emission lines from metals in the UV band are excellent tracers of the cooler fraction of this gas. We present predictions for the surface brightness of a sample of UV lines that could potentially be observed by the next generation of UV telescopes at z<1. We use a subset of simulations from the OWLS project to create emission maps and to investigate the effect of varying the physical prescriptions for star formation, supernova and AGN feedback, chemodynamics and radiative cooling. Most models produce results in agreement within a factor of a few, indicating that the predictions are robust. Of the lines we consider, C III is the strongest line, but it typically traces gas colder than 10^5 K. The same is true for Si IV. The second strongest line, C IV, traces circum-galactic gas with T~10^5 K. O VI and Ne VIII probe the warmer (T~10^5.5 K and T~10^6 K, respectively) and more diffuse gas that may be a better tracer of the large scale structure. N V emission is intermediate between C IV and O VI. The intensity of all emission lines increases strongly with gas density and metallicity, and for the bright emission it is tightly correlated with the temperature for which the line emissivity is highest. In particular, the C III, C IV, Si IV and O VI emission that is sufficiently bright to be potentially detectable in the near future (>10^3 photon/s/cm^2/sr), comes from relatively dense (rho>10^2 rho_mean) and metal rich (Z>0.1 Z_sun) gas. As such, emission lines are highly biased tracers of the missing baryons and are not an optimal tool to close the baryon budget. However, they do provide a powerful means to detect the gas cooling onto or flowing out of galaxies and groups. (Abridged)

The existence of X-ray luminous gaseous coronae around massive disc galaxies is a long-standing prediction of galaxy formation theory in the cold dark matter cosmogony. This prediction has garnered little observational support, with non-detections commonplace and detections for only a relatively small number of galaxies which are much less luminous than expected. We investigate the coronal properties of a large sample of bright, disc-dominated galaxies extracted from the GIMIC suite of cosmological hydrodynamic simulations recently presented by Crain et al. Remarkably, the simulations reproduce the observed scalings of X-ray luminosity with K-band luminosity and star formation rate and, when account is taken of the density structure of the halo, with disc rotation velocity as well. Most of the star formation in the simulated galaxies (which have realistic stellar mass fractions) is fuelled by gas cooling from a quasi-hydrostatic hot corona. However, these coronae are more diffuse, and of a lower luminosity, than predicted by the analytic models of White & Frenk because of a substantial increase in entropy at z ~ 1-3. Both the removal of low entropy gas by star formation and energy injection from supernovae contribute to this increase in entropy, but the latter is dominant for halo masses M_200 <~ 10^(12.5) Msun. Only a small fraction of the mass of the hot gas is outflowing as a wind but, because of its high density and metallicity, it contributes disproportionally to the X-ray emission. The bulk of the X-ray emission, however, comes from the diffuse quasi-hydrostatic corona which supplies the fuel for ongoing star formation in discs today. Future deep X-ray observations with high spectral resolution (e.g. with NeXT/ASTRO-H or IXO) should be able to map the velocity structure of the hot gas and test this fundamental prediction of current galaxy formation theory.

Emission lines from metals offer one of the most promising ways to detect the elusive warm-hot intergalactic medium (WHIM; 10^5 K<T<10^7 K), which is thought to contain a substantial fraction of the baryons in the low-redshift Universe. We present predictions for the soft X-ray line emission from the WHIM using a subset of cosmological simulations from the OverWhelmingly Large Simulations (OWLS) project. We use the OWLS models to test the dependence of the predicted emission on a range of physical prescriptions, such as cosmology, gas cooling and feedback from star formation and accreting black holes. Provided that metal-line cooling is taken into account, the models give surprisingly similar results, indicating that the predictions are robust. Soft X-ray lines trace the hotter part of the WHIM (T>10^6 K). We find that the OVIII 18.97A is the strongest emission line, with a predicted maximum surface brightness of ~10^2 photon/s/cm^2/sr, but a number of other lines are only slightly weaker. All lines show a strong correlation between the intensity of the observed flux and the density and metallicity of the gas responsible for the emission. On the other hand, the potentially detectable emission consistently corresponds to the temperature at which the emissivity of the electronic transition peaks. The emission traces neither the baryonic nor the metal mass. In particular, the emission that is potentially detectable with proposed missions, traces overdense (rho>10^2rho_mean) and metal-rich (Z>0.1Z_sun) gas in and around galaxies and groups. While soft X-ray line emission is therefore not a promising route to close the baryon budget, it does offer the exciting possibility to image the gas accreting onto and flowing out of galaxies.

We investigate the properties of satellite galaxies in cosmological N-body/SPH simulations of galaxy formation in Milky Way-sized haloes. Because of their shallow potential wells, satellite galaxies are very sensitive to heating processes which affect their gas content. Their properties can therefore be used to constrain the nature of feedback processes that regulate galaxy formation. In our simulations, we assume that all the energy produced by supernovae is used as kinetic energy to drive galactic winds. Several of our simulations produce bright, disc-dominated galaxies. We find that wind models in which the wind speed, v_w, is proportional to local velocity dispersion of dark matter, sigma, (and thus the wind mass-loading, eta_w \propto sigma^{-2}) have episodic star formation histories, reproduce the observed satellite luminosity function reasonably well (down to M_v=-7) and match the luminosity-metallicity relation observed in the Local Group satellites. By contrast, models that assume a constant wind speed overproduce faint satellites and predict an incorrect luminosity-metallicity relation. Our simulations therefore suggest that the feedback processes that operate on the scale of satellite galaxies should generate galactic outflows whose intensity varies inversely with the depth of the potential.

[Abridged] The relatively recent insight that energy input from supermassive black holes (BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies motivates us to examine its effects on the scale of galaxy groups. At present, groups contain most of the galaxies and a significant fraction of the overall baryon content of the universe. To explore the effects of BH feedback on groups, we analyse two high resolution cosmological hydro simulations from the OverWhelmingly Large Simulations project. While both include galactic winds driven by supernovae, only one includes feedback from BHs. We compare the properties of the simulated groups to a wide range of observational data, including hot gas radial profiles and gas mass fractions (fgas), luminosity-mass-temperature (L-M-T) scaling relations, K-band luminosity of the group and its central brightest galaxy (CBG), SFRs and ages of the CBG, and gas/stellar metallicities. Both runs yield entropy profiles similar to the data, while the run without AGN feedback yields highly peaked temperature profiles, in discord with the observations. Energy input from BHs significantly reduces fgas for groups with masses less than ~10^14 Msun, yielding fgas-T and L-T relations that are in agreement with the data. The run without AGN feedback suffers from the well known overcooling problem; the resulting K-band luminosities are much larger than observed. By contrast, the run that includes BH feedback yields K-band luminosities and CBG SFRs and ages in agreement with current estimates. Both runs yield very similar gas-phase metallicities that match X-ray data, but they predict very different stellar metallicities. Based on the above, galaxy groups provide a compelling case that BH feedback is a crucial ingredient in the formation of massive galaxies.

We have implemented a simple model to identify the likely sites of the first stars and galaxies in the high-resolution simulations of the formation of galactic dark matter halos of the Aquarius Project. The first star in a galaxy like the Milky Way formed around redshift z=35; by z=10, the young galaxy contained up to 30000 dark matter haloes capable of forming stars by molecular hydrogen cooling. These minihaloes were strongly clustered and feedback may have severely limited the actual number of Population III stars that formed. By the present day, the remnants of the first stars would be strongly concentrated to the centre of the main halo. If a second generation of long-lived stars formed near the first (the first star relics), we would expect to find half of them within 30kpc/h of the galactic centre and a significant fraction in satellites where they may be singled out by their anomalous metallicity patterns. The first halo in which gas could cool by atomic hydrogen line radiation formed at z=25; by z=10, the number of such first galaxies had increased to ~300. Feedback might have decreased the number of first galaxies at the time they undergo significant star formation, but not the number that survive to the present because near neighbours merge. Half of all the "first galaxies" that form before z=10 merge with the main halo before z~3 and most lose a significant fraction of their mass. However, today there should still be more than 20 remnants orbiting within the central ~30kpc/h of the Milky Way. These satellites have circular velocities of a few kilometers per second or more, comparable to those of known Milky Way dwarfs. They are a promising hunting ground for the remnants of the earliest epoch of star formation.

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.

We examine the relationship between galaxies and the intergalactic medium at z < 1 using a group of three closely spaced background QSOs with z_em ~1 observed with the Hubble Space Telescope. Using a new grouping algorithm, we identify groups of galaxies and absorbers across the three QSO sightlines that may be physically linked. There is an excess number of such groups compared to the number we expect from a random distribution of absorbers at a confidence level of 99.9%. The same search is performed with mock spectra generated using a hydrodynamic simulation, and we find the vast majority of such groups arise in dense regions of the simulation. We find that at z<0.5, groups in the simulation generally trace the large-scale filamentary structure as seen in the projected 2-d distribution of the HI column density in a ~30 h^-1 Mpc region. We discover a probable sub-damped Lyman-alpha system at z=0.557 showing strong, low-ionisation metal absorption lines. Previous analyses of absorption across the three sightlines attributed these metal lines to HI. We show that even when the new line identifications are taken into account, evidence remains for planar structures with scales of ~1 Mpc absorbing across the three sightlines. We identify a galaxy at z=0.2272 with associated metal absorption in two sightlines, each 200 kpc away. By constraining the star formation history of the galaxy, we show the gas causing this metal absorption may have been enriched and ejected by the galaxy during a burst of star formation 2 Gyr ago.

(Abridged) We present an implementation of stellar evolution and chemical feedback for smoothed particle hydrodynamics (SPH) simulations. We consider the timed release of individual elements by both massive (Type II supernovae and stellar winds) and intermediate mass stars (Type Ia supernovae and asymptotic giant branch stars). We illustrate the results of our method using a suite of cosmological simulations that include new prescriptions for radiative cooling, star formation, and galactic winds. Radiative cooling is implemented element-by-element, in the presence of an ionizing radiation background, and we track all 11 elements that contribute significantly to the radiative cooling. We contrast two reasonable definitions of the metallicity of a resolution element and find that while they agree for high metallicities, there are large differences at low metallicities. We argue the discrepancy is indicative of the lack of metal mixing caused by the fact that metals are stuck to particles. We argue that since this is a (numerical) sampling problem, solving it using a poorly constrained physical process such as diffusion could have undesired consequences. We demonstrate that the two metallicity definitions result in redshift z = 0 stellar masses that can differ by up to a factor of two, because of the sensitivity of the cooling rates to the elemental abundances. We find that by z = 0 most of the metals are locked up in stars. The gaseous metals are distributed over a very wide range of gas densities and temperatures. The shock-heated warm-hot intergalactic medium has a relatively high metallicity of ~ 10^-1 Z_sun that evolves only weakly and is therefore an important reservoir of metals.

[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.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

(Abridged) The thermal state of the intracluster medium results from a competition between gas cooling and heating. The heating comes from two distinct sources: gravitational heating from the collapse of the dark matter halo and thermal input from galaxy/black hole formation. However, a long standing problem has been that cosmological simulations based on smoothed particle hydrodynamics (SPH) and Eulerian mesh codes predict different results even when cooling and galaxy/black hole heating are switched off. Clusters formed in SPH simulations show near powerlaw entropy profiles, while those formed in mesh simulations develop a core and do not allow gas to reach such low entropies. Since the cooling rate is closely connected to the minimum entropy of the gas, the differences are of potentially key importance. In this paper, we investigate the origin of this discrepancy. By comparing simulations run using the GADGET-2 SPH code and the FLASH adaptive Eulerian mesh code, we show that the discrepancy arises during the idealised merger of two clusters. The difference is not sensitive to the resolution of our simulations, nor is it is due differences in the gravity solvers, Galilean non-invariance of the mesh code, or an effect of unsuitable artificial viscosity in the SPH code. Instead, we find that the difference is inherent to the treatment of eddies and fluid instabilities. These are suppressed in the SPH simulations, while the cluster mergers generate strong vortices in the mesh simulations that efficiently mix the fluid and erase the low entropy gas. Consequently, particles in the SPH simulations retain a close connection to their initial entropy, while this connection is much weaker in the mesh simulations. We discuss the potentially profound implications of these results.

We explore the clustering properties of high redshift dark matter halos, focusing on halos massive enough to host early generations of stars or galaxies at redshift 10 and greater. Halos are extracted from an array of dark matter simulations able to resolve down to the "mini-halo" mass scale at redshifts as high as 30, thus encompassing the expected full mass range of halos capable of hosting luminous objects and sources of reionization. Halo clustering on large-scales agrees with the Sheth, Mo & Tormen halo bias relation within all our simulations, greatly extending the regime where large-scale clustering is confirmed to be "universal" at the 10-20% level (which means, for example, that 3sigma halos of cluster mass at z=0 have the same large-scale bias with respect to the mass distribution as 3sigma halos of galaxy mass at z=10). However, on small-scales, the clustering of our massive halos (> ~10^9 Msun/h) at these high redshifts is stronger than expected from comparisons with small-scale halo clustering extrapolated from lower redshifts. This implies "non-universality" in the scale-dependence of halo clustering, at least for the commonly used parameterizations of the scale-dependence of bias that we consider. We provide a fit for the scale-dependence of bias in our results. This study provides a basis for using extraordinarily high redshift galaxies (redshift ~10) as a probe of cosmology and galaxy formation at its earliest stages. We show also that mass and halo kinematics are strongly affected by finite simulation volumes. This suggests the potential for adverse affects on gas dynamics in hydrodynamic simulations of limited volumes, such as is typical in simulations of the formation of the "first stars", though further study is warranted.

We test four commonly used astrophysical simulation codes; Enzo, Flash, Gadget and Hydra, using a suite of numerical problems with analytic initial and final states. Situations similar to the conditions of these tests, a Sod shock, a Sedov blast and both a static and translating King sphere occur commonly in astrophysics, where the accurate treatment of shocks, sound waves, supernovae explosions and collapsed haloes is a key condition for obtaining reliable validated simulations. We demonstrate that comparable results can be obtained for Lagrangian and Eulerian codes by requiring that approximately one particle exists per grid cell in the region of interest. We conclude that adaptive Eulerian codes, with their ability to place refinements in regions of rapidly changing density, are well suited to problems where physical processes are related to such changes. Lagrangian methods, on the other hand, are well suited to problems where large density contrasts occur and the physics is related to the local density itself rather than the local density gradient.

We perform cosmological hydrodynamic simulations to determine to what extent galaxies lose their gas due to photoheating from an ionizing background. We find that the characteristic mass at which haloes on average have lost half of their baryons is Mc ~ 6.5 x 10^9 Msun/h at z = 0, which corresponds to a circular velocity of 25 km/s. This is significantly lower than the filtering mass obtained by the linear theory, which is often used in semianalytical models of galaxy formation. We demonstrate it is the gas temperature at the virial radius which determines whether a halo can accrete gas. A simple model that follows the merger history of the dark matter progenitors, and where gas accretion is not allowed when this temperature is higher than the virial temperature of the halo, reproduces the results from the simulation remarkably well. This model can be applied to any reionization history, and is easy to incorporate in semianalytical models.

We have studied the abundance of oxygen in the IGM by analyzing OVI, CIV, SiIV, and HI pixel optical depths derived from a set of high-quality VLT and Keck spectra of 17 QSOs at 2.1 < z < 3.6. Comparing OVI and CIV optical depth ratios to those in realistic, synthetic spectra drawn from a hydrodynamical simulation and comparing to existing constraints on [Si/C] places strong constraints on the ultraviolet background (UVB) model using weak priors on allowed values of [Si/O]: for example, a quasar-only background yields [Si/O] ~ 1.4, highly inconsistent with the [Si/O] ~ 0 expected from nucleosynthetic yields and with observations of metal-poor stars. Assuming a fiducial quasar+galaxy UVB consistent with these constraints yields a primary result that [O/C] = 0.66 +/- 0.06 +/- 0.2; this result pertains to gas with overdensity >~ 0.2. Consistent results are obtained by similarly comparing OVI to HI or OVI to SiIV optical depth ratios to simulation values, and also by directly ionization-correcting OVI optical depths as function of HI optical depths into [O/H] as a function of density. Subdividing the sample reveals no evidence for evolution, but low- and high-density samples are inconsistent, suggesting either density-dependence of [O/C] or -- more likely -- prevalence of collisionally-ionized gas at high density.

We have conducted a long slit search for low surface brightness Lyman-alpha emitters at redshift 2.67 < z < 3.75. A 92 hour long exposure with VLT/FORS2 down to a 1-sigma surface brightness detection limit of 8x10^-20 erg/cm2/s/sqarcsec yielded a sample of 27 single line emitters with fluxes of a few times 10^-18 erg/s/cm2. We present arguments that most objects are indeed Lyman-alpha. The large comoving number density, the large covering factor, dN/dz ~ 0.2-1, and the often extended Lyman-alpha emission suggest that the emitters be identified with the elusive host population of damped Lyman-alpha systems (DLAS) and high column density Lyman limit systems. A small inferred star formation rate, perhaps supplanted by cooling radiation, appears to energetically dominate the Lyman-alpha emission, and is consistent with the low metallicity, low dust content, and theoretically inferred low masses of DLAS, and with the relative lack of success of earlier searches for their optical counterparts. (abridged)

We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous halos of galaxies as they fall into massive groups and clusters. The sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ~30% of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott (1972), except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disk) and takes into account that stripping is not instantaneous but occurs on a characteristic timescale. The model reproduces the results of the simulations to within approximately 10% at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10% of the group or cluster into which it is falling (in which case the model under-predicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically-simulated galaxies and as a way to improve current semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.

The first stars in the Universe form when chemically pristine gas heats as it falls into dark matter potential wells, cools radiatively due to the formation of molecular hydrogen, and becomes self-gravitating. We demonstrate with super-computer simulations that their properties depend critically on the currently unknown nature of the dark matter. If the dark matter particles have intrinsic velocities that wipe-out small-scale structure, then the first stars form in filaments with lengths of order the free-streaming scale, which can be about 10^20m (~3kpc, baryonic masses 10^7 solar masses) for realistic "warm dark matter" candidates. Fragmentation of the filaments forms stars with a range of masses which may explain the observed peculiar element abundance pattern of extremely metal-poor stars, while coalescence of fragments and stars during the filament's ultimate collapse may seed the super massive black holes that lurk in the centres of most massive galaxies.

Simulations of an isolated Milky Way-like galaxy, in which supernovae power a galactic fountain, reproduce the observed velocity and 21cm brightness statistics of galactic neutral hydrogen (HI). The simulated galaxy consists of a thin HI disk, similar in extent and brightness to that observed in the Milky Way, and extra-planar neutral gas at a range of velocities due to the galactic fountain. Mock observations of the neutral gas resemble the HI flux measurements from the Leiden-Argentine-Bonn (LAB) HI, survey, including a high-velocity tail which matches well with observations of high-velocity clouds. The simulated high-velocity clouds are typically found close to the galactic disk, with a typical line-of-sight distance of 13kpc from observers on the solar circle. The fountain efficiently cycles matter from the centre of the galaxy to its outskirts at a rate of around 0.5 M_sun/yr

We describe a numerical implementation of star formation in disk galaxies, in which the conversion of cooling gas to stars in the multiphase interstellar medium is governed by the rate at which molecular clouds are formed and destroyed. In the model, clouds form from thermally unstable ambient gas and get destroyed by feedback from massive stars and thermal conduction. Feedback in the ambient phase cycles gas into a hot galactic fountain or wind. We model the ambient gas hydrodynamically using smoothed particle hydrodynamics (SPH). However, we cannot resolve the Jeans mass in the cold and dense molecular gas and, therefore, represent the cloud phase with ballistic particles that coagulate when colliding. We show that this naturally produces a multiphase medium with cold clouds, a warm disk, hot supernova bubbles and a hot, tenuous halo. Our implementation of this model is based on the Gadget N-Body code. We illustrate the model by evolving an isolated Milky Way-like galaxy and study the properties of a disk formed in a rotating spherical collapse. Many observed properties of disk galaxies are reproduced well, including the molecular cloud mass spectrum, the molecular fraction as a function of radius, the Schmidt law, the stellar density profile and the appearance of a galactic fountain.

(Abridged) The thermal history of the intracluster medium (ICM) is complex. Heat input from cluster mergers, AGN, and galaxy winds offsets and may even halt the cooling of the ICM. Consequently, the processes that set the properties of the ICM play a key role in determining how galaxies form. In this paper we focus on the shock heating of the ICM during cluster mergers, with the eventual aim of incorporating this mechanism into semi-analytic models of galaxy formation. We use a suite of hydrodynamic simulations to track the evolution of the ICM in idealised two-body mergers. We find the heating of the ICM can be understood relatively simply by considering the evolution of the gas entropy during the mergers. We examine the processes that generate the entropy in order to understand why previous analytic shock heating models failed. We find that: (1) The energy that is thermalised in the collision greatly exceeds the kinetic energy available when the systems first touch. The smaller system penetrates deep into the potential well before it is disrupted. (2) For unequal mass mergers, most of the energy is thermalised in the more massive component. The heating of the smaller system is minor and its gas sinks to the centre of the final system. (3) The bulk of the entropy generation occurs in two distinct episodes. The first episode occurs following core collision, when a shock wave is generated that propagates outwards from the centre. This causes the combined system to expand rapidly and overshoot hydrostatic equilibrium. The second episode occurs as this material is shock heated as it re-collapses. This revised model for entropy generation significantly improves our physical understanding of cosmological gas simulations.

We use an array of high-resolution N-body simulations to determine the mass function of dark matter haloes at redshifts 10-30. We develop a new method for compensating for the effects of finite simulation volume that allows us to find an approximation to the true ``global'' mass function. By simulating a wide range of volumes at different mass resolution, we calculate the abundance of haloes of mass 10^{5-12} Msun/h. This enables us to predict accurately the abundance of the haloes that host the sources that reionize the universe. In particular, we focus on the small mass haloes (>~10^{5.5-6} Msun/h) likely to harbour population III stars where gas cools by molecular hydrogen emission, early galaxies in which baryons cool by atomic hydrogen emission at a virial temperature of ~10^4K (10^{7.5-8} Msun/h), and massive galaxies that may be observable at redshift ~10. When we combine our data with simulations that include high mass halos at low redshift, we find that the best fit to the halo mass function depends not only on linear overdensity, as is commonly assumed in analytic models, but also upon the slope of the linear power spectrum at the scale of the halo mass. The Press-Schechter model gives a poor fit to the halo mass function in the simulations at all epochs; the Sheth-Tormen model gives a better match, but still overpredicts the abundance of rare objects at all times by up to 50%. Finally, we consider the consequences of the recently released WMAP 3-year cosmological parameters. These lead to much less structure at high redshift, reducing the number of z=10 ``mini-haloes'' by more than a factor of two and the number of z=30 galaxy hosts by more than four orders of magnitude. Code to generate our best-fit halo mass function may be downloaded from http://icc.dur.ac.uk/Research/PublicDownloads/genmf_readme.html

We summarise the science cases for an ELT that were presented in the parallel session on the intergalactic medium, and the open discussion that followed the formal presentations. Observations of the IGM with an ELT provides tremendous potential for dramatic improvements in current programmes in a very wide variety of subjects. These range from fundamental physics (expansion of the Universe, nature of the dark matter, variation of physical constants), cosmology (geometry of the Universe, large-scale structure), reionisation (ionisation state of the IGM at high redshift>6, to more traditional astronomy, such as the interactions between galaxies and the IGM (metal enrichment, galactic winds and other forms of feedback), and the study of the interstellar medium in high redshift galaxies through molecules. The requirements on ELTs and their instruments for fulfilling this potential are discussed.

We model gas cooling in high-resolution N-body simulations in order to investigate the formation of the first generation of stars. We follow a region of a LCDM universe especially selected to contain a rich cluster by the present day. The properties of the dark haloes that form in these sub-solar mass resolution simulations are presented in a companion paper by Gao et al. The first gas clouds able to cool by molecular hydrogen line emission collapse at extremely high redshift, z ~ 47, when the mass of the dark halo is 2.4 x 10^5 Msun/h. By z ~ 30, a substantial population of haloes are capable of undergoing molecular hydrogen cooling although their ability to form stars is dependent on the efficiency of feedback processes such as dissociating Lyman-Werner radiation. The mass of the main halo grows extremely rapidly and, by z ~ 36, its virial temperature has reached 10^4K, at which point gas cooling becomes dominated by more effective atomic line processes. By z ~ 30, a small ``group'' of such potential galaxies will have formed unless prevented from doing so by feedback processes. By this redshift, massive (>~ 100 Msun) population III stars are able to ionise gas well beyond their own host halo and neighbouring HII regions can percolate to form an ionized superbubble. Such patches would be too widely separated to contribute significantly to reionisation at this time. The large number density of early cooling haloes in the pre-reionised universe raises the exciting prospect that this ultra-early generation of stars may be observable as gamma-ray bursts or supernovae.

The physical state of the intergalactic medium can be probed in great detail with the intervening absorption systems seen in quasar spectra. The properties of the Hydrogen absorbers depend on many cosmological parameters, such as the matter-power spectrum, reionisation history, ionising background and the nature of the dark matter. The spectra also contain metal lines, which can be used to constrain the star formation history and the feedback processes acting in large and small galaxies. Simulations have been instrumental in investigating to what extent these parameters can be unambiguously constrained with current and future data. This paper is meant as an introduction to this subject, and reviews techniques and methods for simulating the intergalactic medium.

We investigate the effect of peculiar velocities on the redshift space distribution of z>~2 galaxies, and we focus in particular on Ly-alpha emitters. We generate catalogues of dark matter (DM) halos and identify emitters with halos of the same co-moving space density (M(Ly-alpha emitters) ~ 3x10^11 M_sun). We decompose the peculiar velocity field of halos into streaming, gradient and random components, and compute and analyse these as a function of scale. Streaming velocities are determined by fluctuations on very large scales, strongly affected by sample variance, but have a modest impact on the interpretation of observations. Gradient velocities are the most important as they distort structures in redshift space, changing the thickness and orientation of sheets and filaments. Random velocities are typically below or of the same order as the typical observational uncertainty on the redshift. We discuss the importance of these effects for the interpretation of data on the large-scale structure as traced by Ly-alpha emitters (or similar kinds of astrophysical high-redshift objects), focusing on the induced errors in the viewing angles of filaments. We compare our predictions of velocity patterns for Ly-alpha emitters to observations and find that redshift clumping of Ly-alpha emitters, as reported for instance in the fields of high-redshift radio galaxies, does not allow to infer whether an observed field is sampling an early galaxy overdensity.

We present a method for studying the proximity effect and the density structure around redshift z=2-3 quasars. It is based on the probability distribution of Lyman-alpha pixel optical depths and its evolution with redshift. We validate the method using mock spectra obtained from hydrodynamical simulations, and then apply it to a sample of 12 bright quasars at redshifts 2-3 observed with UVES at the VLT-UT2 Kueyen ESO telescope. These quasars do not show signatures of associated absorption and have a mean monochromatic luminosity of 5.4 10^{31} h^{-2} erg/s/Hz at the Lyman limit. The observed distribution of optical depth within 10 Mpc/h from the QSO is statistically different from that measured in the general intergalactic medium at the same redshift. Such a change will result from the combined effects of the increase in photoionisation rate above the mean UV-background due to the extra ionizing photons from the quasar radiation (proximity effect), and the higher density of the IGM if the quasars reside in overdense regions (as expected from biased galaxy formation). The first factor decreases the optical depth whereas the second one increases the optical depth, but our measurement cannot distinguish a high background from a low overdensity. An overdensity of the order of a few is required if we use the amplitude of the UV-background inferred from the mean Lyman-$\alpha$ opacity. If no overdensity is present, then we require the UV-background to be higher, and consistent with the existing measurements based on standard analysis of the proximity effect.

Using the statistics of pixel optical depths, we compare HI, CIV and CIII absorption in a set of six high quality z ~ 3-4 quasar absorption spectra to that in spectra drawn from two different state-of-the-art cosmological simulations that include galactic outflows. We find that the simulations predict far too little CIV absorption unless the UVB is extremely soft, and always predict far too small CIII/CIV ratios. We note, however, that much of the enriched gas is in a phase (T ~ 10^5-10^7 K, overdensity ~ 0.1-10, Z >~0.1 Z_sol) that should cool by metal line emission -- which was not included in our simulations. When the effect of cooling is modeled, the predicted CIV absorption increases substantially, but the CIII/CIV ratios are still far too small because the density of the enriched gas is too low. Finally, we find that the predicted metal distribution is much too inhomogeneous to reproduce the observed probability distribution of CIV absorption. These findings suggest that strong z <~ 6 winds cannot fully explain the observed enrichment, and that an additional (perhaps higher-z) contribution is required.

The observed cooling rate of hot gas in clusters is much lower than that inferred from the gas density profiles. This suggests that the gas is being heated by some source. We use an adaptive-mesh refinement code (FLASH) to simulate the effect of multiple, randomly positioned, injections of thermal energy within 50 kpc of the centre of an initially isothermal cluster with mass M_200=3x10^(14) Msol and kT=3.1 keV. We have performed eight simulations with spherical bubbles of energy generated every 10^8 years, over a total of 1.5 Gyr. Each bubble is created by injecting thermal energy steadily for 10^7 years; the total energy of each bubble ranges from 0.1--3x10^(60) erg, depending on the simulation. We find that 2x10^(60) erg per bubble (corresponding to a average power of 6.3x10^(44) erg/s) effectively balances energy loss in the cluster and prevents the accumulation of gas below kT=1 keV from exceeding the observational limits of 30 Msol/yr. This injection rate is comparable to the radiated luminosity of the cluster, and the required energy and periodic timescale of events are consistent with observations of bubbles produced by central AGN in clusters. The effectiveness of this process depends primarily on the total amount of injected energy and the initial location of the bubbles, but is relatively insensitive to the exact duty cycle of events.

(Modified) We study the abundance of silicon in the intergalactic medium by analyzing the statistics of SiIV, CIV, and HI pixel optical depths in a sample of 19 high-quality quasar absorption spectra spanning redshifts z ~ 2 - 4, which we compare to realistic spectra drawn from a hydrodynamical simulation. We find that silicon is highly overabundant relative to carbon and that the scatter in Si/C is much smaller than that in C/H, implying a common origin for Si and C. The inferred [Si/C] depends upon the shape of the UV background (UVB) (harder backgrounds leading to higher [Si/C]), ranging from [Si/C] ~ 1.5 for a quasar-only UVB, to [Si/C] ~ 0.25 for a UVB including both galaxies and an artificial softening. For our fiducial UVB, which includes both quasars and galaxies, we find [Si/C]=0.77 +/- 0.05, with a possible systematic error of ~ 0.1 dex. There is no evidence for evolution in [Si/C] and the data are inconsistent with previous claims of a sharp change in the SiIV/CIV ratio (or the UVB) at z ~ 3. Comparisons with low-metallicity halo stars and nucleosynthetic yields suggest that either our fiducial UVB is too hard or that supermassive Pop III stars might have to be included. The inferred [Si/C] corresponds to a contribution to the cosmic Si abundance of [Si/H] = -2.0, a significant fraction of all Si production expected by z ~ 3.

(Modified) We measure the distribution of carbon in the intergalactic medium as a function of redshift z and overdensity delta. Using a hydrodynamical simulation to link the H I absorption to the density and temperature of the absorbing gas, and a model for the UV background radiation, we convert ratios of C IV to H I pixel optical depths into carbon abundances. A careful analysis of 19 high-quality quasar absorption spectra reveals that the carbon abundance is spatially highly inhomogeneous and is well-described by a lognormal distribution for fixed delta and z. Using data in the range log(delta) = -0.5 - 1.8 and z = 1.8 - 4.1, and a renormalized version of the Haardt & Madau (2001) model for the UV background radiation from galaxies and quasars, we measure a median metallicity of [C/H] = -3.47 + 0.08(z-3) + 0.65[log(delta)-0.5] and a lognormal scatter of sigma([C/H]) = 0.76 + 0.02(z-3) - 0.23[log(delta)-0.5]. Thus, we find significant trends with overdensity, but no evidence for evolution. These measurements imply that gas in this density range accounts for a cosmic carbon abundance of [C/H] = -2.80 +/- 0.13 (Omega_C = 2e-7), with no evidence for evolution. The dominant source of systematic error is the spectral shape of the UV background, with harder spectra yielding higher carbon abundances. While the systematic errors due to uncertainties in the spectral hardness may exceed the quoted statistical errors for delta < 10, we stress that UV backgrounds that differ significantly from our fiducial model give unphysical results. The measured lognormal scatter is strictly independent of the spectral shape, provided the background radiation is uniform. We also present measurements of the C III/C IV ratio (which rule out collisional ionization) and of the evolution of the effective Lyman-alpha optical depth.

Using hydrodynamical simulations of the LCDM cosmology, that include both radiative cooling and a phenomenological model for star formation and supernovae feedback, we investigate the impact of galaxy formation on the X-ray properties of groups at redshift zero. Motivated by the observed ``break'' in the L_x-T_x relation at kT~1-2keV, our feedback model is based on the assumption that supernovae imprint a temperature scale on the hot gas, with the star formation rate and corresponding reheated gas mass then depending only on the available energy budget. We demonstrate that a strong feedback model with a heating temperature comparable to this break (kT_=2keV), and an energy budget twice that available from supernovae (epsilon=2), raises the core entropy of groups sufficiently to produce an adequate match to their observed X-ray properties. A lower value of epsilon increases the star formation rate without significantly affecting the X-ray properties of groups, and a model with epsilon~0.1 reproduces the observed fraction of baryons in stars. However, a heating temperature which is lower than the virial temperatures of the groups leads to an excess of cooling gas that boosts their X-ray luminosities, due to the failure of the reheated material to escape from the gravitational potential. A limited study of numerical resolution effects reveals that the temperature of poorly-resolved objects is underestimated, therefore (in our case) a fully-resolved group population would lead to a steeper L_x-T_x relation, bringing our results into even better agreement with the observations.

In quintessence models, the dark energy content of the universe is described by a slowly rolling scalar field whose pressure and energy density obey an equation of state of the form p=w $\rho$; w is in general a function of time such that w<-1/3, in order to drive the observed acceleration of the Universe today. The cosmological constant model LCDM corresponds to the limiting case w=-1. In this paper, we explore the prospects of using the Lyman-alpha forest to constrain w, using semi-analytical techniques to model the intergalactic medium (IGM). A different value of w changes both the growth factor and the Hubble parameter as a function of time. The resulting change in the optical depth distribution affects the optical depth power spectrum, the number of regions of high transmission per unit redshift and the cross-correlation coefficient of spectra of quasar pairs. These can be detected in current data, provided we have independent estimates of the thermal state of the IGM, its ionization parameter and the baryon density.

We have performed hydrodynamical simulations to investigate the effects of galactic winds on the high-redshift (z=3) universe. Strong winds suppress the formation of low-mass galaxies significantly, and the metals carried by them produce C IV absorption lines with properties in reasonable agreement with observations. The winds have little effect on the statistics of the H I absorption lines, because the hot gas bubbles blown by the winds fill only a small fraction of the volume and because they tend to escape into the voids, thereby leaving the filaments that produce these lines intact.

The intergalactic medium at redshifts 2--6 can be studied observationally through the absorption features it produces in the spectra of background quasars. Most of the UV-absorption lines arise in mildly overdense regions, which can be simulated reliably with current hydrodynamical simulations. Comparison of observed and simulated spectra allows one to put contraints on the model's parameters.

Studies of absorption spectra of high-z QSOs have revealed that the intergalactic medium at z ~ 2-3 is enriched to ~ 0.1%-1% solar for gas densities more than a few times the mean cosmic density, but have not yet produced an accurate metallicity estimate, nor constrained variations in the metallicity with density, redshift, or spatial location. This paper discusses the ``pixel optical depth'' (POD) method of QSO spectrum analysis, using realistic simulated spectra from cosmological simulations. In this method, absorption in Ly-alpha is compared to corresponding metal absorption on a pixel-by-pixel basis, yielding for each analyzed spectrum a single statistical correlation encoding metal enrichment information. Our simulations allow testing and optimization of each step of the technique's implementation. Tests show that previous studies have probably been limited by C IV self-contamination and O VI contamination by HI lines; we have developed and tested an effective method of correcting for both contaminants. We summarize these and other findings, and provide a useful recipe for the POD technique's application to observed spectra. Our tests reveal that the POD technique applied to spectra of presently available quality is effective in recovering useful metallicity information even in underdense gas. We present an extension of the POD technique to directly recover the intergalactic metallicity as a function of gas density. For a given ionizing background, both the oxygen and carbon abundance can be measured with errors of at most a factor of a few over at least an order of magnitude in density, using a single high-quality spectrum.

The expansion of the Universe leads to a rapid drop in the mean hydrogen Ly-alpha opacity tau of the intergalactic medium (IGM), tau (1+z)^{3.8}, between redshifts 4 and 3. Measurements of the temperature evolution of the IGM and of the He II opacity both suggest that He II reionizes in this redshift range. We use hydrodynamical simulations to show that the temperature increase associated with He II reionization leads to a relatively sudden decrease in tau around the reionization epoch of ~ 10 per cent. We find clear evidence for such a feature in the evolution of tau determined from a sample of ~ 1100 quasars obtained from the SDSS. He II reionization starts at redshift ~ 3.4, and lasts for (Delta z) ~ 0.4. The increase in the IGM temperature also explains the widths of hydrogen absorption lines as measured in high-resolution spectra.

A variety of approximate schemes for modelling the low-density Intergalactic Medium (IGM) in the high-redshift Universe is compared to the results of a large high-resolution hydro-dynamical simulation. These schemes use either an analytical description of the dark matter distribution and the IGM or numerical simulations of the DM distributions combined with different approximate relations between dark matter field and the gas distribution. Schemes based on a filtering of the dark matter distribution with a global Jeans scale result in a rather poor description of the gas distribution. An adaptive filtering which takes into account the density/temperature dependence of the Jeans scale is required. A reasonable description of the gas distribution can be achieved using a fit of the mean relation between the dark matter and gas densities in the hydro-dynamical simulation to relate dark matter and gas distribution. In the hydro-dynamical simulations deviations from this mean relation are correlated with gradients in the dark matter peculiar velocity field indicative of shocks in the gas component. A scheme which takes into account this correlation results in a further improved gas distribution. Such adaptive filtering schemes applied to dark matter simulations will be very well suited for studies of statistical properties of the Lyalpha forest which investigate the IGM and the underlying dark matter distribution and require a large dynamic range and/or an extensive parameter study.

The temperature of the diffuse, photo-heated intergalactic medium (IGM) depends on its reionization history because the thermal time scales are long. The widths of the hydrogen Lyman-alpha absorption lines seen in the spectra of distant quasars that arise in the IGM can be used to determine its temperature. We use a wavelet analysis of the Lyman-alpha forest region of quasar spectra to demonstrate that there is a relatively sudden increase in the line widths between redshifts z=3.5 and 3.0, which we associate with entropy injection resulting from the reionization of HeII. The subsequent fall-off in temperature after z~3.5, is consistent with a thermal evolution dominated by adiabatic expansion. If, as expected, the temperature also drops rapidly after hydrogen reionization, then the high temperatures inferred from the line widths before HeII reionization imply that hydrogen reionization occurred below redshift z=9.

The temperature of the low-density intergalactic medium is set by the balance between adiabatic cooling resulting from the expansion of the universe, and photo-heating by the UV-background. We have analysed the Lyman-alpha forest of eleven high-resolution quasar spectra using wavelets, and find strong evidence of a marked jump in the temperature at the mean density, T_0, of 60 per cent around a redshift z=3.3, which we attribute to reionization of HeliumII. The jump can be seen in all three of our spectra that straddle redshift 3.3, at a significance of 99 per cent. Below z=3.1, our results are consistent with a smooth cooling down of the universe, as expected when adiabatic expansion dominates over photo-heating by a UV-background from QSOs and galaxies. We find no evidence of thermal fluctuations on scales greater than 5000 km/s larger than 50 per cent, which could be detected by our method, suggesting that the IGM follows a reasonably well-defined temperature-density relation. We demonstrate that the mean wavelet amplitude scales inversely with T_0 and calibrate the relation with hydrodynamical simulations. We find T_0= 12000K at z>3.6. Such high temperature suggest that Hydrogen reionization occured relatively recent.

PINOCCHIO (PINpointing Orbit-Crossing Collapsed Hierarchical Objects) is a new algorithm for identifying dark matter halos in a given numerical realisation of the linear density field in a hierarchical universe (Monaco et al. 2001). It is shown that Lagrangian perturbation theory, and in particular its ellipsoidal truncation, is able to predict accurately the collapse, in the orbit-crossing sense, of generic mass elements. Some points that have undergone orbit crossing are assigned to the network of filaments and sheets that connects the halos; it is demonstrated that this network resembles closely that found in N-body simulations. The code generates a catalogue of dark matter halos with known mass, position, velocity, merging history and angular momentum. It is shown that the predictions of the code are very accurate when compared with the results of large N-body simulations that cover a range of cosmological models, box sizes and numerical resolutions. The mass function is recovered with an accuracy of better than 10 per cent in number density for halos with at least 30-50 particles. A similar accuracy is reached in the estimate of the correlation length r_0. The good agreement is still valid on the object-by-object level, with 70-100 per cent of the objects with more than 50 particles in the simulations also identified by our algorithm. For these objects the masses are recovered with an error of 20-40 per cent, and positions and velocities with a root mean square error of ~1-2 Mpc (0.5-2 grid lengths) and ~100 km/s, respectively. The recovery of the angular momentum of halos is considerably noisier and accuracy at the statistical level is achieved only by introducing free parameters. The algorithm requires negligible computer time as compared with performing a numerical $N$-body simulation.

The temperature of the low-density intergalactic medium is set by the balance between adiabatic cooling resulting from the expansion of the universe, and photo-heating by the UV-background. A sudden injection of entropy from the reionization will increase the temperature of the gas, leading to a broadening of the hydrogen Lyman-alpha absorption lines produced in the IGM, and observed in the spectra of background quasars. We present a method based on wavelets to characterise objectively the line widths of such absorption lines. We use high resolution hydrodynamical simulations to demonstrate that the algorithm can detect changes in temperature of order of 50 per cent on scales > 5000 km/s. We apply the method to a UVES/VLT spectrum of quasar 0055--269 (emission redshift=3.7) and detect at the 99 per cent confidence level a sudden increase in temperature below redshift z=3.3, which we interpret as evidence for HeliumII reionization.

(abridged) We present an effective implementation of analytical calculations of the Lyalpha opacity distribution of the Intergalactic Medium (IGM) along multiple lines of sight (LOS) to distant quasars in a cosmological setting. This method fully accounts for the expected correlations between LOS and the cosmic variance in the large-scale modes of the dark matter distribution. Strong correlations extending up to ~300 kpc (proper) and more are found at redshift z~2->3, in agreement with observations. These correlations are investigated using the cross-correlation coefficient and the cross-power spectrum of the flux distribution along different LOS and by identifying coincident absorption features as fitted with a Voigt profile fitting routine. The cross-correlation coefficient between the LOS can be used to constrain the shape-parameter Gamma of the power spectrum if the temperature and the temperature density relation of the IGM can be determined indepedently. We also propose a new technique to recover the 3D linear dark matter power spectrum by integrating over 1D flux cross-spectra which is complementary to the usual `differentiation' of 1D auto-spectra. The cross-power spectrum suffers much less from errors uncorrelated in different LOS, like those introduced by continuum fitting. Investigations of the flux correlations in adjacent LOS should thus allow to extend studies of the dark matter power spectrum with the Lyalpha forest to significantly larger scales than is possible with flux auto-power spectra. 30 pairs with separation of 1-2 arcmin should be sufficient to determine the 1D cross-spectrum at scales of 60 Mpc/h with an accuracy of about 30% if the error is dominated by cosmic variance.

We study the ability of PINOCCHIO (PINpointing Orbit-Crossing Collapsed HIerarchical Objects) to predict the merging histories of dark matter (DM) haloes, comparing the PINOCCHIO predictions with the results of two large N-body simulations run from the same set of initial conditions. We focus our attention on quantities most relevant to galaxy formation and large-scale structure studies. PINOCCHIO is able to predict the statistics of merger trees with a typical accuracy of 20 per cent. Its validity extends to higher-order moments of the distribution of progenitors. The agreement is valid also at the object-by-object level, with 70-90 per cent of the progenitors cleanly recognised when the parent halo is cleanly recognised itself. Predictions are presented also for quantities that are usually not reproduced by semi-analytic codes, such as the two-point correlation function of the progenitors of massive haloes and the distribution of initial orbital parameters of merging haloes. For the accuracy of the prediction and for the facility with which merger histories are produced, PINOCCHIO provides a means to generate catalogues of DM haloes which is extremely competitive to large-scale N-body simulations, making it a suitable tool for galaxy formation and large-scale structure studies.

We present a new algorithm (PINOCCHIO, PINpointing Orbit-Crossing Collapsed HIerarchical objects) to predict accurately the formation and evolution of individual dark matter haloes in a given realization of an initial linear density field. Compared with the halo population formed in a large (360^3 particles) collisionless simulation of a CDM universe, our method is able to predict to better than 10 per cent statistical quantities such as the mass function, two-point correlation function and progenitor mass function of the haloes. Masses of individual haloes are estimated accurately as well, with errors typically of order 30 per cent in the mass range well resolved by the numerical simulation. These results show that the hierarchical formation of dark matter haloes can be accurately predicted using local approximations to the dynamics when the correlations in the initial density field are properly taken into account. The approach allows one to automatically generate a large ensemble of accurate merging histories of haloes with complete knowledge of their spatial distribution. The construction of the full merger tree for a 256^3 realisation requires a few hours of CPU-time on a personal computer, orders of magnitude faster than the corresponding N-body simulation would take, and without needing any extensive post-processing. The technique can be efficiently used, for instance, for generating the input for galaxy formation modeling.

Models for the formation of galaxies and clusters of galaxies require strong feedback in order to explain the observed properties of these systems. We investigate whether such feedback has observational consequences for the intergalactic medium, as probed in absorption towards background quasars. A typical quasar sight-line intersects one proto-cluster per unit redshift, and significant feedback from forming galaxies or AGN, heating the proto-cluster gas, will result in a large clearing of reduced absorption in the Lyman alpha-forest. Such a gap could be detected at redshift > 3 when the mean opacity is high. Feedback from Lyman-break galaxies in proto-clusters can be probed by the absorption lines produced in their winds. Strong feedback from galaxies has a major impact on the number and properties of absorption lines with neutral hydrogen column densities of 10^16/ cm^2. This feedback can be probed with multiple sight-lines and by studying the unsaturated higher-order lines of the Lyman series. Galactic winds from dwarf galaxies should break-up into clouds, in order not to over produce the number of absorption lines. These clouds can then coast to large distances.

Alignments in the angular momentum vectors of galaxies can induce large scale correlations in their projected orientations. Such alignments arise from the tidal torques exerted on neighboring proto-galaxies by the smoothly varying shear field. Weak gravitational lensing can also induce ellipticity correlations since the images of neighboring galaxies will be distorted coherently by the intervening mass distribution. Comparing these two sources of shape correlations, it is found that for current weak lensing surveys with a median redshift of 1, the intrinsic signal is a contaminant on the order of 1-10% of the measured signal. However, for shallower surveys, the intrinsic correlations dominate over the lensing signal. The distortions induced by lensing are curl-free, whereas those resulting from intrinsic alignments are not. This difference can be used to disentangle these two sources of ellipticity correlations. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors.

The distortion field defined by the ellipticities of galaxy shapes projected on the sky can be uniquely decomposed into a gradient and a curl component. If the observed ellipticities are induced by weak gravitational lensing, then the distortion field is curl free. Here we show that, in contrast, the distortion field resulting from intrinsic spin alignments is not curl free. This provides a powerful discriminant between lensing and intrinsic contributions to observed ellipticity correlations. We also show how these contributions can be disentangled statistically from the ellipticity correlations or computed locally from circular integrals of the ellipticity field. This allows for an unambiguous detection of intrinsic galaxy alignments in the data. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors. These techniques can be applied equally well to the polarization of the microwave background, where it can be used to separate curl-free scalar perturbations from those produced by gravity waves or defects.

Large scale correlations in the orientations of galaxies can result from alignments in their angular momentum vectors. These alignments arise from the tidal torques exerted on neighboring proto-galaxies by the smoothly varying shear field. We compute the predicted amplitude of such ellipticity correlations using the Zel'dovich approximation for a realistic distribution of galaxy shapes. Weak gravitational lensing can also induce ellipticity correlations since the images of neighboring galaxies will be distorted coherently. On comparing these two effects that induce shape correlations, we find that for current weak lensing surveys with a median redshift of z_m = 1, the intrinsic signal is of order 1 - 10 percent of the measured signal. However, for shallower surveys with z_m < 0.3, the intrinsic correlations dominate over the lensing signal. The distortions induced by lensing are curl-free, whereas those resulting from intrinsic alignments are not. This difference can be used to disentangle these two sources of ellipticity correlations.

We investigate the dependence of QSO Lyman-alpha absorption features on the temperature of the absorbing gas and on the amplitude of the underlying dark matter fluctuations. We use high-resolution hydrodynamic simulations in cold dark matter dominated cosmological models. In models with a hotter intergalactic medium (IGM), the increased temperature enhances the pressure gradients between low and high density regions and this changes the spatial distribution and the velocity field of the gas. Combined with more thermal broadening, this leads to significantly wider absorption features in hotter models. Cosmological models with little small scale power also have broader absorption features, because fluctuations on the scale of the Jeans length are still in the linear regime . Consequently, both the amplitude of dark matter fluctuations on small scales and thermal smoothing affect the flux decrement distribution in a similar way. However, the b-parameter distribution of Voigt profile fits, obtained by deblending the absorption features into a sum of thermally broadened lines, is largely independent of the amount of small scale power, but does depend strongly on the IGM temperature. The same is true for the two-point function of the flux and for the flux power spectrum on small scales. These three flux statistics are thus sensitive probes of the temperature of the IGM. We compare them computed for our models and obtained from a HIRES spectrum of the quasar Q1422+231 and conclude that the IGM temperature at z=3.25 is fairly high, T_0>= 15000K.

At redshifts z >~ 2, most of the baryons reside in the smooth intergalactic medium which is responsible for the low column density Lyman-alpha forest. This photoheated gas follows a tight temperature-density relation which introduces a cut-off in the distribution of widths of the Lyman-alpha absorption lines (b-parameters) as a function of column density. We have measured this cut-off in a sample of nine high resolution, high signal-to-noise quasar spectra, and determined the thermal evolution of the intergalactic medium in the redshift range 2.0-4.5. At redshift z ~ 3, the temperature at the mean density shows a peak and the gas becomes nearly isothermal. We interpret this as evidence for the reionization of HeII.

The temperature of the intergalactic medium (IGM) is an important factor in determining the line-widths of the absorption lines in the Lyman-alpha forest. We present a method to characterise the line-widths distribution using a decomposition of a Lyman-alpha spectrum in terms of discrete wavelets. Such wavelets form an orthogonal basis so the decomposition is unique. We demonstrate using hydrodynamic simulations that the mean and dispersion of the wavelet amplitudes is strongly correlated with both the temperature of the absorbing gas and its dependence on the gas density. Since wavelets are also localised in space, we are able to analyse the temperature distribution as a function of position along the spectrum. We illustrate how this method could be used to identify fluctuations in the IGM temperature which might result from late reionization or local effects.

The last few years have seen a dramatic improvement in our understanding of the origin of Lyman alpha absorption systems. Hydrodynamic numerical simulations of cold dark matter dominated universes have shown that the many properties of the Lyman alpha absorption systems can be explained by a photoionized, space-filling, intergalactic medium. Lyman alpha lines offer promising probes of the photoionizing background, the amplitude of the mass fluctuations at high redshift and the evolution of the equation of state of the intergalactic medium.

The thermal evolution of the intergalactic medium (IGM) depends on the reionization history of the universe. Numerical simulations indicate that the low density IGM, which is responsible for the low column density Ly-alpha forest, follows a well defined temperature-density relation. This results in a cut-off in the distribution of line widths as a function of column density. We use hydrodynamic simulations to calibrate the relation between the cut-off and the temperature-density relation and apply this relation to Keck spectra spanning a redshift range z=2-4.5. We find that the temperature peaks at z~3 and interpret this as evidence for reheating due to the reionization of helium.

Recent numerical simulations have lead to a paradigm shift in our understanding of the intergalactic medium, and the loss of a physical justification for Voigt profile fitting of the Lyman-alpha forest. Many individual lines seen in simulated spectra have significant departures from the Voigt profile, yet could be well fitted by a blend of two or more such lines. We discuss the expected effect on the line profiles due to ongoing gravitational structure formation and Hubble expansion. We develop a method to detect departures from Voigt profiles of the absorption lines in a statistical way and apply this method to simulated Lyman-alpha forest spectra, confirming that the profiles seen do statistically differ from Voigt profiles.

(Abridged) Numerical simulations indicate that the intergalactic medium (IGM) responsible for the low column density Lyman-alpha forest follows a well defined temperature-density relation. We demonstrate that such an equation of state results in a cutoff in the distribution of line widths (b-parameters) as a function of column density (N) for the low column density absorption lines. This explains the existence of the lower envelope which is clearly seen in scatter plots of the b(N)-distribution in observed QSO spectra. We show that the parameters of the cutoff in the b(N)-distribution are strongly correlated with the parameters of the underlying equation of state. We use simulations to determine these relations, which can then be applied to the observed cutoff in the b(N)-distribution to measure the equation of state of the IGM. Using Monte Carlo simulations of Keck spectra at z=3, we show that determining the slope of the equation of state will be difficult, but that the amplitude can be determined to within ten per cent, even from a single QSO spectrum. Measuring the evolution of the equation of state with redshift will allow us to put tight constraints on the reionization history of the universe.

Numerical simulations indicate that the smooth, photoionized intergalactic medium (IGM) responsible for the low column density Lyman-alpha forest follows a well defined temperature-density relation. We show that such an equation of state results in a cutoff in the distribution of line widths (b-parameters) as a function of column density (N) for the low column density absorption lines. This explains the existence of the lower envelope which is clearly seen in scatter plots of the b(N)-distribution in observed QSO spectra. The intercept and slope of this cutoff can be used to measure the equation of state of the IGM. Measuring the evolution of the equation of state with redshift will allow us to put tight constraints on the reionization history of the universe.

We have performed high-resolution hydrodynamic simulations of the Lyman-alpha forest in a variety of popular cold dark matter dominated cosmologies, including a low density and a vacuum-dominated model. The fluctuation amplitude of these models is chosen to match the observed abundance of galaxy clusters at low redshift. We assume that the intergalactic medium is photoionized and photoheated by a uniform UV-background with the required amplitude to give the observed mean hydrogen absorption. We produce simulated spectra, analyze them by fitting Voigt profiles and compare line statistics with those obtained from high resolution observations. All models give column density distributions in good agreement with observations. However, the distribution of line widths (the b-parameter distribution) reflects differences in the temperature of the intergalactic medium between the models, with colder models producing more narrow lines. All the models with a low baryon density, Omega_b h^2=0.0125, are too cold to produce a b-parameter distribution in agreement with observations. Models with a higher baryon density, Omega_b h^2=0.025, are hotter and provide better fits. Peculiar velocities contribute significantly to the line widths in models with low matter density, and this improves the agreement with observations further. We briefly discuss alternative mechanisms for reconciling the simulations with the observed b-parameter distributions.

(Abridged) We investigate the importance of several numerical artifacts such as lack of resolution on spectral properties of the Lyman alpha forest as computed from cosmological hydrodynamic simulations in a standard cold dark matter universe. We assume an ionising background produced by quasars as computed by Haardt & Madau. We use a new simulation code based on P3M and SPH, which we compare in detail with a modified version of HYDRA (Couchman et al.) and published results of TREESPH (Hernquist et al.). The agreement is very good between all three codes. We then use our new code to investigate several numerical effects such as resolution on spectral statistics deduced from Voigt profile fitting. Our highest resolution simulation has a mass resolution of 2.1x10^5 solar masses. The column density distribution is converged but the b-parameter distribution is only marginally converged. The simulation reproduces both the HI column density and b-parameter distribution when we assume a high baryon density, Omega_B h^2 > 0.028. In addition we need to impose a higher IGM temperature than predicted within our basic set of assumptions. The simulated HI optical depth is in good agreement with observations but the HeII optical depth is lower than observed. Fitting the latter requires a larger jump between the photon flux at the H and He edge than is present in the Haardt & Madau spectrum.

The low-redshift evolution of the intergalactic medium is investigated using hydrodynamic cosmological simulations. The assumed cosmological model is a critical density cold dark matter universe. The imposed uniform background of ionizing radiation has the amplitude, shape and redshift evolution as computed from the observed quasar luminosity function by Haardt & Madau. We have analysed simulated Lyman-alpha spectra using Voigt-profile fitting, mimicking the procedure with which quasar spectra are analysed. Our simulations reproduce the observed evolution of the number of Lyman-alpha absorption lines over the whole observed interval of z=0.5 to z=4. In particular, our simulations show that the decrease in the rate of evolution of Lyman-alpha absorption lines at z< 2, as observed by Hubble Space Telescope, can be explained by the steep decline in the photo-ionizing background resulting from the rapid decline in quasar numbers at low redshift.

The formalism that describes the non-linear growth of the angular momentum L of protostructures from tidal torques in a Friedmann Universe, as developed in a previous paper, is extended to include non-Gaussian initial conditions. We restrict our analysis here to a particular class of non-Gaussian primordial distributions, namely multiplicative models. In such models, strongly correlated phases are produced by obtaining the gravitational potential via a nonlinear local transformation of an underlying Gaussian random field. The dynamical evolution of the system is followed by describing the trajectories of fluid particles using second-order Lagrangian perturbation theory. In the Einstein-de Sitter universe, the lowest-order perturbative correction to the variance of the linear angular momentum of collapsing structures grows as t^8/3 for generic non-Gaussian statistics, which contrasts with the t^10/3 growth rate characteristic of Gaussian statistics. This is a consequence of the fact that the lowest-order perturbative spin contribution in the non-Gaussian case arises from the third moment of the gravitational potential, which is identically zero for a Gaussian field. Evaluating these corrections at the maximum expansion time of the collapsing structure, we find that these non-Gaussian and non-linear terms can be as high as the linear estimate, without the degree of non-Gaussianity as quantified by skewness and kurtosis of the density field being unacceptably large. The results suggest that higher-order terms in the perturbative expansion may contribute significantly to galactic spin which contrasts with the straightforward Gaussian case.

We have identified ten candidate intergalactic planetary nebulae in the Fornax galaxy cluster. These objects were found during observations in 1992 and 1993 in three fields chosen well away from any Fornax galaxy at 15 arcmin, 30 arcmin, and 45 arcmin from the centre of Fornax. We used the usual method of blinking images taken in a narrow OIII filter, with images taken in an adjacent broad filter. The measured fluxes in the narrow, broad, and I bands are consistent with these unresolved objects being planetary nebulae immersed in an intergalactic population of stars. Such a population is expected to arise as a consequence of tidal encounters between galaxies, and our findings strengthen the case for the existence of such tidal debris. The confirmation of some or all of these ten candidates as planetary nebulae would imply that intergalactic stars constitute a substantial fraction of all the stars in Fornax, up to an estimated 40 per cent. Intergalactic planetary nebulae could prove useful in probing the underlying cluster potential, since they would be far more abundant than galaxies. We discuss possible contamination of the sample by emission-line galaxies, but conclude that planetary nebulae is the most likely identification for the detected objects. keywords: intergalactic medium-- galaxies: interactions -- planetary nebulae: general -- clusters: individual: Fornax

The growth of the angular momentum L of protogalaxies induced by tidal torques is reconsidered within the Zel'dovich approximation. We obtain a general expression for the ensemble expectation value of the square of L in terms of the first and second invariant of the inertia tensor of the Lagrangian volume enclosing the protoobject's collapsing mass. We then specialize the formalism to the particular case in which this volume is centered on a peak of the smoothed Gaussian density field and approximated by an isodensity ellipsoid. The result is the appropriate analytical estimate for the rms angular momentum of peaks to be compared against simulations that make use of the Hoffman-Ribak algorithm to set up a constrained density field that contains a peak with given shape. Extending the work of Heavens & Peacock, we calculate the joint probability distribution function for several spin parameters and peak mass M using the distribution of peak shapes, for different initial power spectra. The values of observed specific angular momentum versus mass are well fitted by our theoretical isoprobability contours. In contrast, the observed lower values for the specific angular momentum for ellipticals of the same mass cannot be accounted for within our linear regime investigation, highlighting the importance of strongly non-linear phenomena to explain the spin of such objects.

We discuss the non-linear evolution of the angular momentum L acquired by protostructures, like protogalaxies and protoclusters, due to tidal interactions with the surrounding matter inhomogeneities. The primordial density distribution is assumed to be Gaussian and the non-linear dynamics of the collisionless mass fluid is followed using Lagrangian perturbation theory. For a Cold Dark Matter spectrum, the inclusion of the leading-order Lagrangian correction terms results in a value of the rms ensemble average <L^2>^{1/2} which is only a factor of 1.3 higher than the corresponding linear estimate, irrespective of the scale. Consequently, the predictions of linear theory are rather accurate in quantifying the evolution of the angular momentum of protostructures before collapse sets in. In the Einstein-de Sitter universe, the initial torque is a good estimate for the tidal torque over the whole period during which the object is spun up.

Smoothed particle hydrodynamics (SPH) is used to estimate accretion rates of mass, linear and angular momentum in a binary system where one component undergoes mass loss through a wind. Physical parameters are chosen such as to model the alleged binary precursors of barium stars, whose chemical peculiarities are believed to result from the accretion of the wind from a companion formerly on the asymptotic giant branch (AGB). The binary system modelled consists of a 3 solar masses AGB star on the main sequence, in a 3AU circular orbit. Three-dimensional simulations are performed for gases with polytropic indices gamma=1, 1.1 and 1.5, to bracket more realistic situations that would include radiative cooling. Mass accretion rates are found to depend on resolution and we estimate typical values of 1-2% for the gamma=1.5 case and 8% for the other models. The highest resolution obtained (with 400k particles) corresponds to an accretor of linear size 16 solar radii. Despite being (in the gamma = 1.5 case) about ten times smaller than theoretical estimates based on the Bondi-Hoyle prescription, the SPH accretion rates remain large enough to explain the pollution of barium stars. Uncertainties in the current SPH rates remain however, due to the simplified treatment of the wind acceleration mechanism, as well as to the absence of any cooling prescription and to the limited numerical resolution. Angular momentum transfer leads to significant spin up of the accretor and can account for the rapid rotation of HD165141, a barium star with a young white dwarf companion and a rotation rate unusually large among K giants.

The energy diffusion coefficients D_n(E) (n=1,2) for a system of equal mass particles moving self-consistently in an N-body realisation of a King model are computed from the probability per unit time, P(E, Delta E), that a star with initial energy E will undergo an energy change Delta E. In turn, P is computed from the number of times during the simulation that a particle in a state of given energy undergoes a transition to another state. These particle states are defined directly from the time evolution of E by identifying them with the event occuring between two local maxima in the E(t) curve. If one assumes next that energy changes are uncorrelated between different states, one can use diffusion theory to compute D_n(E). The simulations employ N=512, 2048,... , 32768 particles and are performed using an implementation of Aarseth's direct integrator N-body1 on a massively parallel computer. The more than seven million transitions measured in the largest N simulation provide excellent statistics. The numerically determined D(E)'s are compared against their theoretical counterparts which are computed from phase-space averaged rates of energy change due to independent binary encounters. The overall agreement between them is impressive over most of the energy range, notwithstanding the very different type of approximations involved, giving considerable support to the valid usage of these theoretical expressions to simulate dynamical evolution in Fokker-Planck type calculations.

A P3M (Particle-Particle, Particle-Mesh) algorithm to compute the gravitational force on a set of particles is described. The gravitational force is computed using Fast Fourier Transforms. This leads to an incorrect force when the distance between two particles is of the order of a grid cell. This incorrect force is subtracted exactly from all particles in parallel using convolution with the appropriate Green's function in real space in a time of order $N_T$, irrespective of the degree of clustering of particles. Next, the correct $1/r^2$ force is added for all neighbouring particles in parallel, leading to an accurate algorithm which runs efficiently on a highly parallel computer. A full force calculation for $128k$ particles on a $128^3$ grid in a mildly clustered situation requires approximately 196 seconds on a $8k$ Connection Machine 2 with 8MHz clock. This decreases to an estimated 9.8 seconds on a full-sized $64k$ CM200.