We explored a suite of high-resolution cosmological simulations from the First Billion Years Project (FiBY) at $z \geq 6$. All substructures within the simulations have been identified with the SUBFIND algorithm. From our analysis, two distinct groups of objects emerge. We hypothesise that the substructures in the first group, which appear to have a high baryon fraction ($f_{\rm b} \geq 0.95$), are possible infant GC candidates. Objects belonging to the second group have a high stellar fraction ($f_{\rm star} \geq 0.95$) and show a potential resemblance to infant ultra-faint dwarf galaxies. The high baryon fraction objects identified in this study are characterised by a stellar content similar to the one observed in present-day GCs, but they still contain a high gas fraction ($f_{\rm gas} \sim 0.95$) and a relatively low amount of dark matter. They are compact, dense systems. Their sizes are consistent with recent estimates based on the first observations of possible proto-GCs at high redshifts. These types of infant GC candidates appear to be more massive and more abundant in massive host galaxies, indicating that the assembly of galaxies via mergers may play an important role in building several GC-host scaling relations. Specifically, we express the relation between the mass of the most massive infant GC and its host stellar mass as $\log(M_{\rm cl}) = (0.31\pm0.15)\log(M_{\rm *,gal} + (4.17\pm1.06)$. We also report a new relation between the most massive infant GC and the parent specific star formation rate of the form $\log(M_{\rm cl}) = (0.85\pm0.30)\log(sSFR) + \alpha$ that describes the data at both low and high redshift. Finally, we assess the present-day GC mass (GC number) -- halo mass relation offers a satisfactory description of the behaviour of our infant GC candidates at high redshift, suggesting that such a relation may be set at formation.

We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, ${\gamma}_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of ${\gamma}_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of ${\gamma}_{\rm gNFW}$.

We investigate higher order symplectic integration strategies within Bayesian cosmic density field reconstruction methods. In particular, we study the fourth-order discretisation of Hamiltonian equations of motion (EoM). This is achieved by recursively applying the basic second-order leap-frog scheme (considering the single evaluation of the EoM) in a combination of even numbers of forward time integration steps with a single intermediate backward step. This largely reduces the number of evaluations and random gradient computations, as required in the usual second-order case. We restrict this study to the lognormal-Poisson model, applied to a full volume halo catalogue in real space on a cubical mesh of 1250 $h^{-1}$ Mpc side and 256$^3$ cells. Hence, we neglect selection effects, redshift space distortions, and displacements. We note that those observational and cosmic evolution effects can be accounted for in subsequent Gibbs-sampling steps within the COSMIC BIRTH algorithm. We find that going from the usual second to fourth-order in the leap-frog scheme shortens the burn-in phase by a factor of at least $\sim30$. This implies that 80-100 independent samples are obtained while the fastest second-order method converges. After convergence, the correlation lengths indicate an improvement factor of $\sim2.4$ fewer gradient computations. This gain in computational efficiency can help to go towards a full Bayesian analysis of the cosmological large-scale structure for upcoming galaxy surveys.

We have developed a new tool to analyse galaxies in the EAGLE simulations as close as possible to observations.We investigated the evolution of their kinematic properties by means of the angular momentum proxy parameter,$ \lambda_{Re} $for galaxies with $M_{*} \ge 5 \times 10^{9} M_{\odot}$ in the RefL0100N1504 simulation up to redshift two (z = 2). Galaxies in the simulation show a wide variety of kinematic features, similiar to those found in integral-field spectroscopic studies. At z=0 the distribution of galaxies in the ${\lambda}_{Re}-{\epsilon}$ plane is also in good agreement with results from observations. Scaling relations at z = 0 indicate that there is critical mass, $M_{crit} = 10^{10.3} M_{\odot}$, that divides two different regimes when we include the ${\lambda}_{Re}$ parameter. The simulation shows that the distribution of galaxies in the ${\lambda}_{Re}-{\epsilon}$ plane evolves with time until z = 2 when galaxies are equally distributed both in ${\lambda}_{Re}$ and ${\epsilon}$. We studied the evolution of ${\lambda}_{Re}$ with time and found that there is no connection between the angular momentum at z = 2 and z = 0. All systems reach their maximum ${\lambda}_{Re}$ at z = 1 and then steadily lose angular momentum regardless of their merger history, except for the high star-forming systems that sustain that maximum value over time. The evolution of the Re in galaxies that have not experienced any merger in the last 10 Gyr can be explained by their level of gas accretion.

By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-density contours using the procedure of Montes \& Trujillo (2019), and find that the shape of the stellar mass distribution follows that of the total matter even more closely than observed, although their radial profiles differ substantially. The ratio between stellar and total matter density profiles in circular apertures, shows a slope close to -1, with a small dependence on the cluster's total mass. We propose an indirect method to calculate the halo mass and mass density profile from the radial profile of the intra-cluster stellar mass density.

Our location in the Milky Way provides an exceptional opportunity to gain insight on the galactic evolution processes, and complement the information inferred from observations of external galaxies. Since the Milky Way is a barred galaxy, the study of motions of individual stars in the bulge and disc is useful to understand the role of the bar. The Gaia mission enables such study by providing the most precise parallaxes and proper motions to date. In this theoretical work, we explore the effects of the bar on the distribution of higher-order moments --the skewness and kurtosis-- of the proper motions by confronting two simulated galaxies, one with a bar and one nearly axisymmetric, with observations from the latest Gaia data release (GaiaDR2). We introduce the code ASGAIA to account for observational errors of Gaia in the kinematical structures predicted by the numerical models. As a result, we find clear imprints of the bar in the skewness distribution of the longitudinal proper motion $\mu_\ell$ in GaiaDR2, as well as other features predicted for the next Gaia data releases.

Early-type galaxies are considered to be the end-products of massive galaxy formation. Optical spectroscopic studies reveal that massive early-type galaxies formed the bulk of their stars over short timescales (1<Gyr) and at high redshift (z>2), followed by passive evolution to the present. However, their optical spectra are insensitive to constrain small episodes of recent star formation, since they are dominated by old stars. Fortunately, this problem can be tackled in the ultraviolet range. While recent studies that make use of ultraviolet absorption lines have suggested the presence of young stars in a few early-type galaxies, the age and mass fractions of young stars and their dependence on galaxy mass, is unknown. Here we report a detailed study of these young stellar populations, from high-quality stacked spectra of 28,663 galaxies from the BOSS survey, analysing optical and ultraviolet absorption lines simultaneously. We find that residual star formation is ubiquitous in massive early-type galaxies, measuring average mass fractions of ~0.5% in young stars in the last 2 Gyr of their evolution. This fraction shows a decreasing trend with galaxy stellar mass, consistent with a down-sizing scenario. We also find that synthetic galaxies from state-of-the-art cosmological numerical simulations significantly overproduce both intermediate and young stellar populations. Therefore, our results pose stringent constraints on numerical simulations of galaxy formation.

In this paper we demonstrate that the information encoded in \emph{one} single (sufficiently large) $N$-body simulation can be used to reproduce arbitrary numbers of halo catalogues, using approximated realisations of dark matter density fields with different initial conditions. To this end we use as a reference one realisation (from an ensemble of $300$) of the Minerva $N$-body simulations and the recently published Bias Assignment Method to extract the local and non-local bias linking the halo to the dark matter distribution. We use an approximate (and fast) gravity solver to generate $300$ dark matter density fields from the down-sampled initial conditions of the reference simulation and sample each of these fields using the halo-bias and a kernel, both calibrated from the arbitrarily chosen realisation of the reference simulation. We show that the power spectrum, its variance and the three-point statistics are reproduced within $\sim 2\%$ (up to $k\sim1.0\,h\,{\rm Mpc}^{-1}$), $\sim 5-10\%$ and $\sim 10\%$, respectively. Using a model for the real space power spectrum (with three free bias parameters), we show that the covariance matrices obtained from our procedure lead to parameter uncertainties that are compatible within $\sim 10\%$ with respect to those derived from the reference covariance matrix, and motivate approaches that can help to reduce these differences to $\sim 1\%$. Our method has the potential to learn from one simulation with moderate volumes and high-mass resolution and extrapolate the information of the bias and the kernel to larger volumes, making it ideal for the construction of mock catalogues for present and forthcoming observational campaigns such as Euclid or DESI.

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.

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 dark ages of the Universe end with the formation of the first generation of stars residing in primeval galaxies. These objects were the first to produce ultraviolet ionizing photons in a period when the cosmic gas changed from a neutral state to an ionized one, known as Epoch of Reionization (EoR). A pivotal aspect to comprehend the EoR is to probe the intertwined relationship between the fraction of ionizing photons capable to escape dark haloes, also known as the escape fraction ($f_{esc}$), and the physical properties of the galaxy. This work develops a sound statistical model suitable to account for such non-linear relationships and the non-Gaussian nature of $f_{esc}$. This model simultaneously estimates the probability that a given primordial galaxy starts the ionizing photon production and estimates the mean level of the $f_{esc}$ once it is triggered. The model was employed in the First Billion Years simulation suite, from which we show that the baryonic fraction and the rate of ionizing photons appear to have a larger impact on $f_{esc}$ than previously thought. A naive univariate analysis of the same problem would suggest smaller effects for these properties and a much larger impact for the specific star formation rate, which is lessened after accounting for other galaxy properties and non-linearities in the statistical model.

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.

Forthcoming large-scale spectroscopic surveys will soon provide data on thousands of galaxy clusters. It is important that the systematics of the various mass estimation techniques are well understood and calibrated. We compare three different dynamical mass estimators using the C-EAGLE galaxy clusters, a set of high resolution simulations with resolved galaxies a median total mass, $M_{200c} = 10^{14.7} \, \mathrm{M_\odot}$. We quantify the bias and scatter of the Jeans, virial, and caustic mass estimators using all galaxies with a stellar mass $M_*> 10^9 \, \mathrm{M_\odot}$, both in the ideal 3D case and in the more realistic projected case. On average we find our mass estimates are unbiased, though relative to the true mass within $r_{200c}$ the scatter is large with a range of $0.09$ - $0.15$ dex. We see a slight increase in the scatter when projecting the clusters. Selecting galaxies using the same criteria, we find no significant difference in the mass bias or scatter when comparing results from hydrodynamical and dark matter only simulations. However, selecting galaxies by stellar mass reduces the bias compared to selecting by total mass. Comparing X-ray derived hydrostatic and dynamical masses, the former are ${\sim} 30$ per cent lower. We find a slight dependence between substructure, measured using two different metrics, and mass bias. In conclusion, we find that dynamical mass estimators, when averaged together, are unbiased with a scatter of $0.11 \pm 0.02$ dex when including interloper galaxies and with no prior knowledge of $r_{200c}$.

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 compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.

Bars are common galactic structures in the local universe that play an important role in the secular evolution of galaxies, including the Milky Way. In particular, the velocity distribution of individual stars in our galaxy is useful to shed light on stellar dynamics, and provides information complementary to that inferred from the integrated light of external galaxies. However, since a wide variety of models reproduce the distribution of velocity and the velocity dispersion observed in the Milky Way, we look for signatures of the bar on higher-order moments of the line-of-sight velocity ($V_{los}$) distribution. We make use of two different numerical simulations --one that has developed a bar and one that remains nearly axisymmetric-- to compare them with observations in the latest APOGEE data release (SDSS DR14). This comparison reveals three interesting structures that support the notion that the Milky Way is a barred galaxy. A high skewness region found at positive longitudes constrains the orientation angle of the bar, and is incompatible with the orientation of the bar at $\ell=0^\circ$ proposed in previous studies. We also analyse the $V_{los}$ distributions in three regions, and introduce the Hellinger distance to quantify the differences among them. Our results show a strong non-Gaussian distribution both in the data and in the barred model, confirming the qualitative conclusions drawn from the velocity maps. In contrast to earlier work, we conclude it is possible to infer the presence of the bar from the kurtosis distribution.

For the same stellar mass, physically smaller star-forming galaxies are also metal richer (Ellison et al. 2008). What causes the relation remains unclear. The central star-forming galaxies in the EAGLE cosmological numerical simulation reproduce the observed trend. We use them to explore the origin of the relation assuming that the physical mechanism responsible for the anti-correlation between size and gas-phase metallicity is the same in the simulated and the observed galaxies. We consider the three most likely causes: (1) metal-poor gas inflows feeding the star-formation process, (2) metal-rich gas outflows particularly efficient in shallow gravitational potentials, and (3) enhanced efficiency of the star-formation process in compact galaxies. Outflows (2) and enhanced star-formation efficiency (3) can be discarded. Metal-poor gas inflows (1) cause the correlation in the simulated galaxies. Galaxies grow in size with time, so those that receive gas later are both metal poorer and larger, giving rise to the observed anti-correlation. As expected within this explanation, larger galaxies have younger stellar populations. We explore the variation with redshift of the relation, which is maintained up to, at least, redshift 8.

We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions.

Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of volumes up to $(100\,\rm Mpc)^3$, while Hydrangea consists of $24$ cosmological simulations of galaxy clusters and their environment. Together they provide a statistically significant sample in the stellar mass range $10^{9.5}\,\rm M_{\odot}-10^{12.3}\,\rm M_{\odot}$, of $16,358$ galaxies. We construct IFS-like cubes and measure stellar spin parameters, $\lambda_{\rm R}$, and ellipticities, allowing us to classify galaxies into slow/fast rotators as in observations. The simulations display a primary dependence of $F_{\rm SR}$ on stellar mass, with a weak dependence on environment. At fixed stellar mass, satellite galaxies are more likely to be SRs than centrals. $F_{\rm SR}$ shows a dependence on halo mass at fixed stellar mass for central galaxies, while no such trend is seen for satellites. We find that $\approx 70$% of SRs at $z=0$ have experienced at least one merger with mass ratio $\ge 0.1$, with dry mergers being at least twice more common than wet mergers. Individual dry mergers tend to decrease $\lambda_{\rm R}$, while wet mergers mostly increase it. However, $30$% of SRs at $z=0$ have not experienced mergers, and those inhabit halos with median spins twice smaller than the halos hosting the rest of the SRs. Thus, although the formation paths of SRs can be varied, dry mergers and/or halos with small spins dominate.

We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2<z<3$ in the \emph{Hubble} \emph{Ultra} \emph{Deep} \emph{Field} (HUDF), as well as empirical $A_V - M_{\star}$ relations predicted by a Calzetti-like law. Our results, combined with other literature data, suggest that the $A_V - M_{\star}$ relation does not evolve over the redshift range $0<z<5$, at least for galaxies with log$(M_{\star}/M_{\odot}) \gtrsim 9.5$. Finally, we present tentative evidence which suggests that the attenuation curve may become steeper at log$(M_{\star}/M_{\odot}) \lesssim 9.0$.

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.

The diffuse soft X-ray emissivity from galactic winds is computed during the Epoch of Reionization (EoR). We consider two analytic models, a pressure-driven wind and a superbubble model, and a 3D cosmological simulation including gas dynamics from the First Billion Years (FiBY) project. The analytic models are normalized to match the diffuse X-ray emissivity of star-forming galaxies in the nearby Universe. The cosmological simulation uses physically motivated star formation and wind prescriptions, and includes radiative transfer corrections. The models and the simulation all are found to produce sufficient heating of the Intergalactic Medium to be detectable by current and planned radio facilities through 21 cm measurements during the EoR. While the analytic models predict a 21 cm emission signal relative to the Cosmic Microwave Background sets in by $z_{\rm trans} \simeq 8 - 10$, the predicted signal in the FiBY simulation remains in absorption until reionization completes. The 21 cm absorption differential brightness temperature reaches a minimum of $\Delta T \simeq -130$ to $-200$ mK, depending on model. Allowing for additional heat from high mass X-ray binaries pushes the transition to emission to $z_{\rm trans} \simeq 10 - 12$, with shallower absorption signatures having a minimum of $\Delta T \simeq -110$ to $-140$ mK. The 21 cm signal may be a means of distinguishing between the wind models, with the superbubble model favouring earlier reheating. While an early transition to emission may indicate X-ray binaries dominate the reheating, a transition to emission as early as $z_{\rm trans} > 12$ would suggest the presence of additional heat sources.

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 present the results of a study investigating the dust attenuation law at $z\simeq 5$, based on synthetic spectral energy distributions (SEDs) calculated for a sample of N=498 galaxies drawn from the First Billion Years (FiBY) simulation project. The simulated galaxies at $z\simeq 5$, which have M$_{1500} \leq -18.0$ and $7.5 \leq \rm{log(M/M}_{\odot}\rm{)} \leq 10.2$, display a mass-dependent $\alpha$-enhancement, with a median value of $[\alpha/\rm{Fe}]_{z=5}~\simeq~4~\times~[\alpha/\rm{Fe}]_{Z_{\odot}}$. The median Fe/H ratio of the simulated galaxies is $0.14\pm0.05$ which, even including the effects of nebular continuum, produces steep intrinsic UV continuum slopes; $\langle \beta_{i} \rangle = -2.4 \pm 0.05$. Using a set of simple dust attenuation models, in which the wavelength-dependent attenuation is assumed to be of the form $A(\lambda) \propto \lambda^{n}$, we explore the parameter values which best reproduce the observed $z=5$ luminosity function (LF) and colour-magnitude relation (CMR). We find that a simple model in which the absolute UV attenuation is a linearly increasing function of log stellar mass, and the dust attenuation slope ($n$) is within the range $-0.7 \leq n \leq-0.3$, can successfully reproduce the LF and CMR over a wide range of stellar population synthesis model (SPS) assumptions. This range of attenuation curves is consistent with a power-law fit to the Calzetti attenuation law in the UV ($n=-0.55$), and other similarly `grey' star-forming galaxy attenuation curves recently derived at $z\simeq2$. In contrast, attenuation curves as steep as the Small Magellanic Cloud (SMC) extinction curve ($n=-1.24$) are formally ruled out. Finally, we show that our models are consistent with recent 1.3mm ALMA observations of the Hubble Ultra Deep Field (HUDF), and predict the form of the $z\simeq5$ IRX$-\beta$ relation.

We investigate the environment in which direct-collapse black holes may form by analysing a cosmological, hydrodynamical simulation that is part of the First Billion Years project. This simulation includes the most relevant physical processes leading to direct collapse of haloes, most importantly, molecular hydrogen depletion by dissociation of $H_2$ and $H^-$ from the evolving Lyman-Werner radiation field. We selected a sample of pristine atomic cooling haloes that have never formed stars in their past, have not been polluted with heavy elements and are cooling predominantly via atomic hydrogen lines. Amongst them we identified six haloes that could potentially harbour massive seed black holes formed via direct collapse (with masses in the range of $10^{4-6} M_{sun}$). These potential hosts of direct-collapse black holes form as satellites and are found within 15 physical kpc of proto-galaxies, with stellar masses in the range $10^{5-7} M_{sun}$ and maximal star formation rates of 0.1 Msun/yr over the past 5 Myr, and are exposed to the highest flux of Lyman-Werner radiation emitted from the neighbouring galaxies. It is the proximity to these proto-galaxies that differentiates these haloes from rest of the sample.

We present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution towards observable galaxies at $z \sim 6$. We focus on three different galaxies which end up in halos with masses $M_{h} = 2.4 \times10^{10}~h^{-1}\; M_{\odot}$ (Halo-10), $1.6 \times10^{11}~h^{-1}\; M_{\odot}$ (Halo-11) and $0.7 \times10^{12}~h^{-1} M_{\odot}$ (Halo-12) at z=6. Our simulations also probe impacts of different sub-grid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that star formation occurs intermittently due to supernova (SN) feedback at z > 10, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no-feedback or lower SF efficiency models can sustain galactic disk for long periods > 10 Myr. The expulsion of gas at the galactic center also affects the inner dark matter density profile. However, SN feedback does not seem to keep the shallow profile of dark matter for a long period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the star formation rates (SFR) and stellar masses of observed Lyman-$\alpha$ emitters (LAEs) at z = 7-8 fairly well given observational uncertainties. In addition, we investigate the effect of UV background radiation on star formation as an external feedback source, and find that earlier reionization extends the quenching time of star formation due to photo-ionization heating, but does not affect the stellar mass at z=6.

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.

Galaxies in clusters are strongly affected by their environment. They evolve according to several physical mechanisms that are active in clusters. Their efficiency can strongly depend on the orbital configuration of the galaxies. Our aim is to analyse the orbits of the galaxies in the cluster Abell 85, based on the study of the galaxy velocity anisotropy parameter. We have solved the Jeans equation under the assumption that the galaxies in A85 are collisionless objects, within the spherically symmetric gravitational potential of the virialized cluster. The mass of the cluster was estimated with X-ray and caustic analyses. We find that the anisotropy profile of the full galaxy population in A85 is an increasing monotonic function of the distance from the cluster centre: on average, galaxies in the central region (r/r200 < 0.3) are on isotropic orbits, while galaxies in the outer regions are on radial orbits. We also find that the orbital properties of the galaxies strongly depend on their stellar colour. In particular, blue galaxies are on less radial orbits than red galaxies. The different families of cluster galaxies considered here have the pseudo phase-space density profiles Q(r) and Qr(r) consistent with the profiles expected in virialized dark matter halos in $N$-body simulations. This result suggests that the galaxies in A85 have reached dynamical equilibrium within the cluster potential. Our results indicate that the origin of the blue and red colour of the different galaxy populations is the different orbital shape rather than the accretion time.

We compare the predicted conditional mass function (CMF) of dark matter halos from two theoretical prescriptions against numerical N-body simulations, both in overdense and underdense regions and at different Eulerian scales ranging from $5$ to $30\,h^{-1}\,$Mpc. In particular, we consider in detail a locally-implemented rescaling of the unconditional mass function (UMF) already discussed in the literature, and also a generalization of the standard rescaling method described in the extended Press-Schechter formalism. First, we test the consistency of these two rescalings by verifying the normalization of the CMF at different scales, and showing that none of the proposed cases provides a normalized CMF. In order to satisfy the normalization condition, we include a modification in the rescaling procedure. After this modification, the resulting CMF generally provides a better description of numerical results. We finally present an analytical fit to the ratio between the CMF and the UMF (also known as the matter-to-halo bias function) in underdense regions, which could be of special interest to speed-up the computation of the halo abundance when studying void statistics. In this case, the CMF prescription based on the locally-implemented rescaling provides a slightly better description of the numerical results when compared to the standard rescaling.

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {\AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

Local extremely metal-poor (XMP) galaxies are of particular astrophysical interest since they allow us to look into physical processes characteristic of the early Universe, from the assembly of galaxy disks to the formation of stars in conditions of low metallicity. Given the luminosity-metallicity relationship, all galaxies fainter than Mr < -13 are expected to be XMPs. Therefore, XMPs should be common in galaxy surveys. However, they are not, because several observational biases hamper their detection. This work compares the number of faint XMPs in the SDSS-DR7 spectroscopic survey with the expected number, given the known biases and the observed galaxy luminosity function. The faint end of the luminosity function is poorly constrained observationally, but it determines the expected number of XMPs. Surprisingly, the number of observed faint XMPs (around 10) is over-predicted by our calculation, unless the upturn in the faint end of the luminosity function is not present in the model. The lack of an upturn can be naturally understood if most XMPs are central galaxies in their low-mass dark matter halos, which are highly depleted in baryons due to interaction with the cosmic ultraviolet background and to other physical processes. Our result also suggests that the upturn towards low luminosity of the observed galaxy luminosity function is due to satellite galaxies.

We extract cosmological information from the anisotropic power spectrum measurements from the recently completed Baryon Oscillation Spectroscopic Survey (BOSS), extending the concept of clustering wedges to Fourier space. Making use of new FFT-based estimators, we measure the power spectrum clustering wedges of the BOSS sample by filtering out the information of Legendre multipoles l > 4. Our modelling of these measurements is based on novel approaches to describe non-linear evolution, bias, and redshift-space distortions, which we test using synthetic catalogues based on large-volume N-body simulations. We are able to include smaller scales than in previous analyses, resulting in tighter cosmological constraints. Using three overlapping redshift bins, we measure the angular diameter distance, the Hubble parameter, and the cosmic growth rate, and explore the cosmological implications of our full shape clustering measurements in combination with CMB and SN Ia data. Assuming a {\Lambda}CDM cosmology, we constrain the matter density to {\Omega}_m = 0.311 -0.010 +0.009 and the Hubble parameter to H_0 = 67.6 -0.6 +0.7 km s^-1 Mpc^-1, at a confidence level (CL) of 68 per cent. We also allow for non-standard dark energy models and modifications of the growth rate, finding good agreement with the {\Lambda}CDM paradigm. For example, we constrain the equation-of-state parameter to w = -1.019 -0.039 +0.048. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. 2016 to produce the final cosmological constraints from BOSS.

We introduce a new suite of radiation-hydrodynamical simulations of galaxy formation and reionization called Aurora. The Aurora simulations make use of a spatially adaptive radiative transfer technique that lets us accurately capture the small-scale structure in the gas at the resolution of the hydrodynamics, in cosmological volumes. In addition to ionizing radiation, Aurora includes galactic winds driven by star formation and the enrichment of the universe with metals synthesized in the stars. Our reference simulation uses 2x512^3 dark matter and gas particles in a box of size 25 comoving Mpc/h with a force softening scale of at most 0.28 kpc/h. It is accompanied by simulations in larger and smaller boxes and at higher and lower resolution, employing up to 2x1024^3 particles, to investigate numerical convergence. All simulations are calibrated to yield simulated star formation rate (SFR) functions in close agreement with observational constraints at redshift z = 7 and to achieve reionization at z = 8.3, which is consistent with the observed optical depth to reionization. We focus on the design and calibration of the simulations and present some first results. The median stellar metallicities of low-mass galaxies at z = 6 are consistent with the metallicities of dwarf galaxies in the Local Group, which are believed to have formed most of their stars at high redshifts. After reionization, the mean photoionization rate decreases systematically with increasing resolution. This coincides with a systematic increase in the abundance of neutral hydrogen absorbers in the IGM.

For several decades it has been known that stellar bars in disc galaxies can be triggered by interactions, or by internal processes such as dynamical instabilities. In this work, we explore the differences between these two mechanisms using numerical simulations. We perform two groups of simulations based on isolated galaxies, one group in which a bar develops naturally, and another group in which the bar could not develop in isolation. The rest of the simulations recreate 1:1 coplanar fly-by interactions computed with the impulse approximation. The orbits we use for the interactions represent the fly-bys in groups or clusters of different masses accordingly to the velocity of the encounter. In the analysis we focus on bars' amplitude, size, pattern speed and their rotation parameter, ${\cal R}=R_{CR}/R_{bar}$. The latter is used to define fast (${\cal R}<1.4$) and slow rotation (${\cal R}>1.4$). Compared with equivalent isolated galaxies we find that bars affected or triggered by interactions: (i) remain in the slow regime for longer; (ii) are more boxy in face-on views; (iii) they host kinematically hotter discs. Within this set of simulations we do not see strong differences between retrograde or prograde fly-bys. We also show that slow interactions can trigger bar formation.

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.

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 present a new deep spectroscopic catalogue for Abell 85, within 3.0 $\times$ 2.6 Mpc$^2$ and down to $M_{r} \sim M_{r}^* +6$. Using the Visible Multi-Object Spectrograph at the Very Large Telescope (VIMOS@VLT) and the AutoFiber 2 at the William Herschel Telescope (AF2@WHT), we obtained almost 1,430 new redshifts for galaxies with $m_r \leq 21$ mag and $\langle \mu_{e,r} \rangle \leq 24$ mag arcsec$^{-2}$. These redshifts, together with SDSS-DR6 and NED spectroscopic information, result in 460 confirmed cluster members. This dataset allows the study of the luminosity function (LF) of the cluster galaxies covering three orders of magnitudes in luminosities. The total and radial LFs are best modelled by a double Schechter function. The normalized LFs show that their bright ($M_{r} \leq -21.5$) and faint ($M_{r}\geq -18.0$) ends are independent of clustercentric distance and similar to the field LFs unlike the intermediate luminosity range ($-21.5 \leq M_{r} \leq -18.0$). Similar results are found for the LFs of the dominant types of galaxies: red, passive, virialized and early-infall members. On the contrary, the LFs of blue, star forming, non-virialized and recent-infall galaxies are well described by a single Schechter function. These populations contribute to a small fraction of the galaxy density in the innermost cluster region. However, in the outskirts of the cluster, they have similar densities to red, passive, virialized and early-infall members at the LF faint end. These results confirm a clear dependence of the colour and star formation of Abell 85 members in the cluster centric distance.

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

Massive high-redshift quiescent compact galaxies (nicknamed red nuggets) have been traditionally connected to present-day elliptical galaxies, often overlooking the relationships that they may have with other galaxy types. We use large bulge-disk decomposition catalogues based on the Sloan Digital Sky Survey (SDSS) to check the hypothesis that red nuggets have survived as compact cores embedded inside the haloes or disks of present-day massive galaxies. In this study, we designate a "compact core" as the bulge component that satisfies a prescribed compactness criterion. Photometric and dynamic mass-size and mass-density relations are used to show that, in the inner regions of galaxies at z ~ 0.1, there are "abundant" compact cores matching the peculiar properties of the red nuggets, an abundance comparable to that of red nuggets at z ~ 1.5. Furthermore, the morphology distribution of the present-day galaxies hosting compact cores is used to demonstrate that, in addition to the standard channel connecting red nuggets with elliptical galaxies, a comparable fraction of red nuggets might have ended up embedded in disks. This result generalises the inside-out formation scenario; present-day massive galaxies can begin as dense spheroidal cores (red nuggets), around which either a spheroidal halo or a disk are formed later.

Measurements of the redshift-space galaxy clustering have been a prolific source of cosmological information in recent years. Accurate covariance estimates are an essential step for the validation of galaxy clustering models of the redshift-space two-point statistics. Usually, only a limited set of accurate N-body simulations is available. Thus, assessing the data covariance is not possible or only leads to a noisy estimate. Further, relying on simulated realisations of the survey data means that tests of the cosmology dependence of the covariance are expensive. With these points in mind, this work presents a simple theoretical model for the linear covariance of anisotropic galaxy clustering observations with synthetic catalogues. Considering the Legendre moments (`multipoles') of the two-point statistics and projections into wide bins of the line-of-sight parameter (`clustering wedges'), we describe the modelling of the covariance for these anisotropic clustering measurements for galaxy samples with a trivial geometry in the case of a Gaussian approximation of the clustering likelihood. As main result of this paper, we give the explicit formulae for Fourier and configuration space covariance matrices. To validate our model, we create synthetic HOD galaxy catalogues by populating the haloes of an ensemble of large-volume N-body simulations. Using linear and non-linear input power spectra, we find very good agreement between the model predictions and the measurements on the synthetic catalogues in the quasi-linear regime.

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

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.

Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionizing photons needed to reionize the Universe, due to their high escape fractions of ionizing photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form, by computing escape fractions in more than 75000 haloes between redshifts 27 and 6 that were extracted from the First Billion Years project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the gas column density in a radius of 10 pc around the stellar populations, causing a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 1e8 solar mass results in low column densities that can be penetrated by radiation from young stars (age < 5 Myr). In haloes with virial mass above 1e8 solar mass supernova feedback is important, but only 30% of the haloes in this mass range have an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. This makes it very hard to predict the escape fraction on the basis of halo properties and results in a highly anisotropic escape fraction. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionization and is something current models of cosmic reionization should strive to take into account.

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.

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.

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 present a suite of cosmological radiation-hydrodynamical simulations of the assembly of galaxies driving the reionization of the intergalactic medium (IGM) at z >~ 6. The simulations account for the hydrodynamical feedback from photoionization heating and the explosion of massive stars as supernovae (SNe). Our reference simulation, which was carried out in a box of size 25 comoving Mpc/h using 2 x 512^3 particles, produces a reasonable reionization history and matches the observed UV luminosity function of galaxies. Simulations with different box sizes and resolutions are used to investigate numerical convergence, and simulations in which either SNe or photoionization heating or both are turned off, are used to investigate the role of feedback from star formation. Ionizing radiation is treated using accurate radiative transfer at the high spatially adaptive resolution at which the hydrodynamics is carried out. SN feedback strongly reduces the star formation rates (SFRs) over nearly the full mass range of simulated galaxies and is required to yield SFRs in agreement with observations. Photoheating helps to suppress star formation in low-mass galaxies, but its impact on the cosmic SFR is small. Because the effect of photoheating is masked by the strong SN feedback, it does not imprint a signature on the UV galaxy luminosity function, although we note that our resolution is insufficient to model star-forming minihaloes cooling through molecular hydrogen transitions. Photoheating does provide a strong positive feedback on reionization because it smooths density fluctuations in the IGM, which lowers the IGM recombination rate substantially. Our simulations demonstrate a tight non-linear coupling of galaxy formation and reionization, motivating the need for the accurate and simultaneous inclusion of photoheating and SN feedback in models of the early Universe.

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.

To provide a quantitative cosmological context to ongoing observational work on the formation histories and location of compact massive galaxies, we locate and study a sample of exceptionally compact systems in the Bolshoi simulation, using the dark matter structural parameters from a real, compact massive galaxy (NGC 1277) as a basis for our working criteria. We find that over 80% of objects in this nominal compact category are substructures of more massive groups or clusters, and that the probability of a given massive substructure being this compact increases significantly with the mass of the host structure; rising to ~30% for the most massive clusters in the simulation. Tracking the main progenitors of this subsample back to z=2, we find them all to be distinct structures with scale radii and densities representative of the population as a whole at this epoch. What does characterise their histories, in addition to mostly becoming substructures, is that they have almost all experienced below-average mass accretion since z=2; a third of them barely retaining, or even losing mass during the intervening 10 Gyr.

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 have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

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.

Long gamma-ray burst's (LGRB's) association to the death of massive stars suggest they could be used to probe the cosmic star formation history (CSFH) with high accuracy, due to their high luminosities. We utilise cosmological simulations from the First Billion Years project to investigate the biases between the CSFH and the LGRB rate at z>5, assuming various different models and constraints on the progenitors of LGRBs. We populate LGRBs using a selection based on environmental properties and demonstrate that the LGRB rate should trace the CSFH to high redshifts. The measured LGRB rate suggests that LGRBs have opening angles of theta_jet=0.1 deg, although the degeneracy with the progenitor model cannot rule out an underlying bias. We demonstrate that proxies that relate the LGRB rate with global LGRB host properties do not reflect the underlying LGRB environment, and are in fact a result of the host galaxy's spatial properties, such that LGRBs can exist in galaxies of solar metallicity. However, we find a class of host galaxies that have low stellar mass and are metal-rich, and that their metallicity dispersions would not allow low-metallicity environments. Detection of hosts with this set of properties would directly reflect the progenitor's environment. We predict that 10% of LGRBs per year are associated with this set of galaxies that would have forbidden line emission that could be detected by instruments on the James Webb Space Telescope. Such a discovery would place strong constraints on the collapsar model and suggest other avenues to be investigated.

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 study the effect of baryons on the inner dark matter profile of the first galaxies using the First Billion Years simulation between z=16-6 before secular evolution sets in. Using a large statistical sample from two simulations of the same volume and cosmological initial conditions, one with and one without baryons, we are able to directly compare haloes with their baryon-free counterparts, allowing a detailed study of the modifications to the dark matter density profile due to the presence of baryons during the first billion years of galaxy formation. For each of the ~ 5000 haloes in our sample we quantify the impact of the baryons using eta, defined as the ratio of dark matter mass enclosed in 100 pc in the baryonic run to its counterpart without baryons. During this epoch of rapid growth of galaxies, we find that many haloes of these first galaxies show an enhancement of dark matter in the halo centre compared to the baryon-free simulation, while many others show a deficit. We find that the mean value of eta is close to unity, but there is a large dispersion, with a standard deviation of 0.677. The enhancement is cyclical in time and tracks the star formation cycle of the galaxy; as gas falls to the centre and forms stars, the dark matter moves in as well. Supernova feedback then removes the gas, and the dark matter again responds to the changing potential. We study three physical models relating the motion of baryons to that of the dark matter: adiabatic contraction, dynamical friction, and rapid outflows. Abridged, see text for full abstract

We use cosmological hydrodynamical simulations to investigate how the inclusion of physical processes relevant to galaxy formation (star formation, metal-line cooling, stellar winds, supernovae and feedback from Active Galactic Nuclei, AGN) change the properties of haloes, over four orders of magnitude in mass. We find that gas expulsion and the associated dark matter (DM) expansion induced by supernova-driven winds are important for haloes with masses M200 < 10^13 Msun, lowering their masses by up to 20% relative to a DM-only model. AGN feedback, which is required to prevent overcooling, has a significant impact on halo masses all the way up to cluster scales (M200 ~ 10^15 Msun). Baryonic physics changes the total mass profiles of haloes out to several times the virial radius, a modification that cannot be captured by a change in the halo concentration. The decrease in the total halo mass causes a decrease in the halo mass function of about 20%. This effect can have important consequences for abundance matching technique as well as for most semi-analytic models of galaxy formation. We provide analytic fitting formulae, derived from simulations that reproduce the observed baryon fractions, to correct halo masses and mass functions from DM-only simulations. The effect of baryonic physics (AGN feedback in particular) on cluster number counts is about as large as changing the cosmology from WMAP7 to Planck, even when a moderately high mass limit of M500 ~ 10^14 Msun is adopted. Thus, for precision cosmology the effects of baryons must be accounted for.

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.

The observed clustering of galaxies and the cross-correlation of galaxies and mass provide important constraints on both cosmology and models of galaxy formation. Even though the dissipation and feedback processes associated with galaxy formation are thought to affect the distribution of matter, essentially all models used to predict clustering data are based on collisionless simulations. Here, we use large hydrodynamical simulations to investigate how galaxy formation affects the autocorrelation functions of galaxies and subhaloes, as well as their cross-correlation with matter. We show that the changes due to the inclusion of baryons are not limited to small scales and are even present in samples selected by subhalo mass. Samples selected by subhalo mass cluster ~10% more strongly in a baryonic run on scales r > 1Mpc/h, and this difference increases for smaller separations. While the inclusion of baryons boosts the clustering at fixed subhalo mass on all scales, the sign of the effect on the cross-correlation of subhaloes with matter can vary with radius. We show that the large-scale effects are due to the change in subhalo mass caused by the strong feedback associated with galaxy formation and may therefore not affect samples selected by number density. However, on scales r < r_vir significant differences remain after accounting for the change in subhalo mass. We conclude that predictions for galaxy-galaxy and galaxy-mass clustering from models based on collisionless simulations will have errors greater than 10% on sub-Mpc scales, unless the simulation results are modified to correctly account for the effects of baryons on the distributions of mass and satellites.

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 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 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 measure the dark matter halo masses of z~2.36 UV color-selected star-forming galaxies by matching the observed median HI Ly{\alpha} absorption around them, as observed in the spectra of background QSOs, to the absorption around haloes above a given mass in cosmological simulations. Focusing on transverse separations 0-2 pMpc and line of sight separations 154-616 km/s, we find a minimum halo mass of log10(Mmin/Msun)=11.6(+0.2)(-0.2), which is in good agreement with published halo mass estimates from clustering analyses. We verified that the measured halo mass is insensitive to a change in the cosmological parameters (WMAP1 vs. WMAP3) and to the inclusion of strong AGN feedback. One unique strength of this method is that it can be used in narrow field galaxy-QSO surveys, i.e. ~30 x 30 arcseconds. In addition, we find that the observed anisotropy in the 2-D HI Ly{\alpha} absorption distribution on scales of 1.5-2 pMpc is consistent with being a consequence of large-scale gas infall into the potential wells occupied by galaxies.

We use numerical simulations to investigate how the statistical properties of dark matter (DM) haloes are affected by the baryonic processes associated with galaxy formation. We focus on how these processes influence the spin and shape of a large number of DM haloes covering a wide range of mass scales, from galaxies to clusters at redshifts zero and one, extending to dwarf galaxies at redshift two. The haloes are extracted from the OverWhelmingly Large Simulations, a suite of state-of-the-art high-resolution cosmological simulations run with a range of feedback prescriptions. We find that the median spin parameter in DM-only simulations is independent of mass, redshift and cosmology. At z=0 baryons increase the spin of the DM in the central region (<=0.25 r_200) by up to 30 per cent when feedback is weak or absent. This increase can be attributed to the transfer of angular momentum from baryons to the DM, but is no longer present at z=2. We also present fits to the mass dependence of the DM halo shape at both low and high redshift. At z=0 the sphericity (triaxiality) is negatively (positively) correlated with halo mass and both results are independent of cosmology. Interestingly, these mass-dependent trends are markedly weaker at z=2. While the cooling of baryons acts to make the overall DM halo more spherical, stronger feedback prescriptions tend to reduce the impact of baryons by reducing the central halo mass concentration. More generally, we demonstrate a strongly positive (negative) correlation between halo sphericity (triaxiality) and galaxy formation efficiency, with the latter measured using the central halo baryon fraction. In conclusion, our results suggest that the effects of baryons on the DM halo spin and shape are minor when the effects of cooling are mitigated, as required by realistic models of galaxy formation, although they remain significant for the inner halo.

The contribution of stars in galaxies to cosmic reionisation depends on the star formation history in the Universe, the abundance of galaxies during reionisation, the escape fraction of ionising photons and the clumping factor of the inter-galactic medium (IGM). We compute the star formation rate and clumping factor during reionisation in a cosmological volume using a high-resolution hydrodynamical simulation. We post-process the output with detailed radiative transfer simulations to compute the escape fraction of ionising photons. Together, this gives us the opportunity to assess the contribution of galaxies to reionisation self-consistently. The strong mass and redshift dependence of the escape fraction indicates that reionisation occurred between z=15 and z=10 and was mainly driven by proto-galaxies forming in dark-matter haloes with masses between 1e7 and 1e8 solar mass. More massive galaxies that are rare at these redshifts and have significantly lower escape fractions contribute less photons to the reionisation process than the more-abundant low-mass galaxies. Star formation in the low-mass haloes is suppressed by radiative feedback from reionisation, therefore these proto-galaxies only contribute when the part of the Universe they live in is still neutral. After z~10, massive galaxies become more abundant and provide most of the ionising photons. In addition, we find that Population (Pop) III stars are too short-lived and not frequent enough to have a major contribution to reionisation. Although the stellar component of the proto-galaxies that produce the bulk of ionising photons during reionisation is too faint to be detected by the James Webb Space Telescope (JWST), these sources are brightest in the H-alpha and Ly-alpha recombination lines, which will likely be detected by JWST in deep surveys.

We study for the first time the environment of massive black hole (BH) seeds (~10^4-5 Msun) formed via the direct collapse of pristine gas clouds in massive haloes (>10^7 Msun) at z>6. Our model is based on the evolution of dark matter haloes within a cosmological N-body simulation, combined with prescriptions for the formation of BH along with both Pop III and Pop II stars. We calculate the spatially-varying intensity of Lyman Werner (LW) radiation from stars and identify the massive pristine haloes in which it is high enough to shut down molecular hydrogen cooling. In contrast to previous BH seeding models with a spatially constant LW background, we find that the intensity of LW radiation due to local sources, J_local, can be up to 10^6 times the spatially averaged background in the simulated volume and exceeds the critical value, J_crit, for the complete suppression of molecular cooling, in some cases by 4 orders of magnitude. Even after accounting for possible metal pollution in a halo from previous episodes of star formation, we find a steady rise in the formation rate of direct collapse (DC) BHs with decreasing redshift from 10^{-3}/Mpc^3/z at z=12 to 10^{-2}/Mpc^3/z at z=6. The onset of Pop II star formation at z~16 simultaneously marks the onset of the epoch of DCBH formation, as the increased level of LW radiation from Pop II stars is able to elevate the local levels of the LW intensity to J_local > J_crit while Pop III stars fail to do so at any time. The number density of DCBHs is sensitive to the number of LW photons and can vary by an order of magnitude at z=6 after accounting for reionisation feedback. Haloes hosting DCBHs are more clustered than similar massive counterparts that do not host DCBHs, especially at redshifts z>10. We also show that planned surveys with JWST should be able to detect the supermassive stellar precursors of DCBHs.

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.

Cosmological simulations make use of sub-grid recipes for the implementation of galactic winds driven by massive stars because direct injection of supernova energy in thermal form leads to strong radiative losses, rendering the feedback inefficient. We argue that the main cause of the catastrophic cooling is a mismatch between the mass of the gas in which the energy is injected and the mass of the parent stellar population. Because too much mass is heated, the temperatures are too low and the cooling times too short. We use analytic arguments to estimate, as a function of the gas density and the numerical resolution, the minimum heating temperature that is required for the injected thermal energy to be efficiently converted into kinetic energy. We then propose and test a stochastic implementation of thermal feedback that uses this minimum temperature increase as an input parameter and that can be employed in both particle- and grid-based codes. We use smoothed particle hydrodynamics simulations to test the method on models of isolated disc galaxies in dark matter haloes with total mass 10^10 and 10^12 h^-1 solar masses. The thermal feedback strongly suppresses the star formation rate and can drive massive, large-scale outflows without the need to turn off radiative cooling temporarily. In accord with expectations derived from analytic arguments, for sufficiently high resolution the results become insensitive to the imposed temperature jump and also agree with high-resolution simulations employing kinetic feedback.

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.

With the first metal enrichment by Population (Pop) III supernovae (SNe), the formation of the first metal-enriched, Pop II stars becomes possible. In turn, Pop III star formation and early metal enrichment are slowed by the high energy radiation emitted by Pop II stars. Thus, through the SNe and radiation they produce, Populations II and III coevolve in the early Universe, one regulated by the other. We present large (4 Mpc)^3, high resolution cosmological simulations in which we self-consistently model early metal enrichment and the stellar radiation responsible for the destruction of the coolants (H2 and HD) required for Pop III star formation. We find that the molecule-dissociating stellar radiation produced both locally and over cosmological distances reduces the Pop III star formation rate at z > 10 by up to an order of magnitude compared to the case in which this radiation is not included. However, we find that the effect of LW feedback is to enhance the amount of Pop II star formation. We attribute this to the reduced rate at which gas is blown out of dark matter haloes by SNe in the simulation with LW feedback, which results in larger reservoirs for metal-enriched star formation. Even accounting for metal enrichment, molecule-dissociating radiation and the strong suppression of low-mass galaxy formation due to reionization at z < 10, we find that Pop III stars are still formed at a rate of ~ 10^-5 M_sun yr^-1 Mpc^-3 down to z ~ 6. This suggests that the majority of primordial pair-instability SNe that may be uncovered in future surveys will be found at z < 10. We also find that the molecule-dissociating radiation emitted from Pop II stars may destroy H2 molecules at a high enough rate to suppress gas cooling and allow for the formation of supermassive primordial stars which collapse to form ~ 100,000 solar mass black holes.

One of the most promising explanations for the origin of the billion solar mass black holes (BHs) inferred to power quasars at redshifts z > 6 is that supermassive stars (SMSs) with masses > 10,000 solar masses collapse to form the seed BHs from which they grow. Here we review recent theoretical advances which provide support for this scenario. Firstly, given sufficiently high accretion rates of gas into the cores of primordial protogalaxies, it appears that neither the high energy radiation emitted from the stellar surface nor the limited lifetime of SMSs can prevent their growth to masses of up to > 100,000 solar masses. Secondly, recent cosmological simulations suggest that the high fluxes of molecule-dissociating radiation which may be required in order to achieve such high accretion rates may be more common in the early universe than previously thought. We conclude that the majority of supermassive BHs may originate from SMSs at high redshifts.

We explore the dynamical signatures imprinted by baryons on dark matter haloes during the formation process using the OverWhelmingly Large Simulations (OWLS), a set of state-of-the-art high-resolution cosmological hydrodynamical simulations. We present a detailed study of the effects of the implemented feedback prescriptions on the orbits of dark matter particles, stellar particles and subhaloes, analysing runs with no feedback, with stellar feedback and with feedback from supermassive black holes. We focus on the central regions (0.25 r_{200}) of haloes with virial masses ~ 6 x 10^{13} (~ 7 x 10^{11}) Msun/h at z = 0(2). We also investigate how the orbital content (relative fractions of the different orbital types) of these haloes depends on several key parameters such as their mass, redshift and dynamical state. The results of spectral analyses of the orbital content of these simulations are compared, and the change in fraction of box, tube and irregular orbits is quantified. Box orbits are found to dominate the orbital structure of dark matter haloes in cosmological simulations. There is a strong anticorrelation between the fraction of box orbits and the central baryon fraction. While radiative cooling acts to reduce the fraction of box orbits, strong feedback implementations result in a similar orbital distribution to that of the dark matter only case. The orbital content described by the stellar particles is found to be remarkably similar to that drawn from the orbits of dark matter particles, suggesting that either they have forgotten their dynamical history, or that subhaloes bringing in stars are not biased significantly with respect to the main distribution. The orbital content of the subhaloes is in broad agreement with that seen in the outer regions of the particle distributions.

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.

In this work we develop a new method to turn a state-of-the-art hydrodynamical cosmological simulation of galaxy formation (HYD) into a simple semi-analytic model (SAM). This is achieved by summarizing the efficiencies of accretion, cooling, star formation, and feedback given by the HYD, as functions of the halo mass and redshift. Surprisingly, by turning the HYD into a SAM, we conserve the mass of individual galaxies, with deviations at the level of 0.1 dex, on an object-by-object basis. This is true for all redshifts, and for the mass of stars and gas components, although the agreement reaches 0.2 dex for satellite galaxies at low redshift. We show that the same level of accuracy is obtained even in case the SAM uses only one phase of gas within each galaxy. Moreover, we demonstrate that the formation history of one massive galaxy provides sufficient information for the SAM to reproduce the population of galaxies within the entire cosmological box. The reasons for the small scatter between the HYD and SAM galaxies are: a) The efficiencies are matched as functions of the halo mass and redshift, meaning that the evolution within merger-trees agrees on average. b) For a given galaxy, efficiencies fluctuate around the mean value on time scales of 0.2-2 Gyr. c) The various mass components of galaxies are obtained by integrating the efficiencies over time, averaging out these fluctuations. We compare the efficiencies found here to standard SAM recipes and find that they often deviate significantly. For example, here the HYD shows smooth accretion that is less effective for low mass haloes, and is always composed of hot or dilute gas; cooling is less effective at high redshift; and star formation changes only mildly with cosmic time. The method developed here can be applied in general to any HYD, and can thus serve as a common language for both HYDs and SAM (Abridged).

Observations of quasar absorption line systems reveal that the z=3 intergalactic medium (IGM) is polluted by heavy elements down to HI optical depths tau_HI<<10. What is not yet clear, however, is what fraction of the volume needs to be enriched by metals and whether it suffices to enrich only regions close to galaxies in order to reproduce the observations. We use gas density fields derived from large cosmological simulations, together with synthetic quasar spectra and imposed, model metal distributions to investigate what enrichment patterns can reproduce the observed median optical depth of CIV as a function of tau_HI. Our models can only satisfy the observational constraints if the z=3 IGM was primarily enriched by galaxies that reside in low-mass (m_tot<10^10 M_sun) haloes that can eject metals out to distances >10^2 kpc. Galaxies in more massive haloes cannot possibly account for the observations as they are too rare for their outflows to cover a sufficiently large fraction of the volume. Galaxies need to enrich gas out to distances that are much greater than the virial radii of their host haloes. Assuming the metals to be well mixed on small scales, our modeling requires that the fractions of the simulated volume and baryonic mass that are polluted with metals are, respectively, >10% and >50% in order to match observations.

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations.

We perform simulations of feedback from supernovae and black holes with smoothed particle hydrodynamics (SPH). Such strong perturbations are inaccurately handled with standard time integration schemes, leading to poor energy conservation, a problem that is commonly overlooked. We show for the first time that, in the absence of radiative cooling, concordance of thermal and kinetic feedback are achieved when using an accurate time integration. In order to preserve the concordance of feedback methods when using a more efficient time integration scheme - as for instance the hierarchical time-step scheme - we implement a modified version of the time-step limiter proposed by Saitoh & Makino (2009). We apply the limiter to general test cases, and first show that this scheme violates energy conservation up to almost four orders of magnitude when energy is injected at random times. To tackle this issue, we find necessary, not only to ensure a fast information propagation, but also to enforce a prompt response of the system to the energy perturbation. The method proposed here to handle strong feedback events enables us to achieve energy conservation at percent level in all tests, even if all the available energy is injected into only one particle. We argue that concordance of feedback methods can be achieved in numerical simulations only if the time integration scheme preserve a high energy conservation level.

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.

We investigate the correlation between nine different dark matter halo properties using a rank correlation analysis and a Principal Component Analysis for a sample of haloes spanning five orders of magnitude in mass. We consider mass and dimensionless measures of concentration, age, relaxedness, sphericity, triaxiality, substructure, spin, and environment, where the latter is defined in a way that makes it insensitive to mass. We find that concentration is the most fundamental property. Except for environment, all parameters are strongly correlated with concentration. Concentration, age, substructure, mass, sphericity and relaxedness can be considered a single family of parameters, albeit with substantial scatter. In contrast, spin, environment, and triaxiality are more independent, although spin does correlate strongly with substructure and both spin and triaxiality correlate substantially with concentration. Although mass sets the scale of a halo, all other properties are more sensitive to concentration.

Upcoming weak lensing surveys, such as LSST, EUCLID, and WFIRST, aim to measure the matter power spectrum with unprecedented accuracy. In order to fully exploit these observations, models are needed that, given a set of cosmological parameters, can predict the non-linear matter power spectrum at the level of 1% or better for scales corresponding to comoving wave numbers 0.1<k<10 h/Mpc. We have employed the large suite of simulations from the OWLS project to investigate the effects of various baryonic processes on the matter power spectrum. In addition, we have examined the distribution of power over different mass components, the back-reaction of the baryons on the CDM, and the evolution of the dominant effects on the matter power spectrum. We find that single baryonic processes are capable of changing the power spectrum by up to several tens of per cent. Our simulation that includes AGN feedback, which we consider to be our most realistic simulation as, unlike those used in previous studies, it has been shown to solve the overcooling problem and to reproduce optical and X-ray observations of groups of galaxies, predicts a decrease in power relative to a dark matter only simulation ranging, at z=0, from 1% at k~0.3 h/Mpc to 10% at k~1 h/Mpc and to 30% at k~10 h/Mpc. This contradicts the naive view that baryons raise the power through cooling, which is the dominant effect only for k>70 h/Mpc. Therefore, baryons, and particularly AGN feedback, cannot be ignored in theoretical power spectra for k>0.3 h/Mpc. It will thus be necessary to improve our understanding of feedback processes in galaxy formation, or at least to constrain them through auxiliary observations, before we can fulfil the goals of upcoming weak lensing surveys.

We use a combination of three large N-body simulations to investigate the dependence of dark matter halo concentrations on halo mass and redshift in the WMAP year 5 cosmology. The median relation between concentration and mass is adequately described by a power-law for halo masses in the range 10^11 - 10^15 Msol/h and redshifts z < 2, regardless of whether the halo density profiles are fit using NFW or Einasto profiles. Compared with recent analyses of the Millennium Simulation, which uses a value of sigma_8 that is higher than allowed by WMAP5, z = 0 halo concentrations are reduced by factors ranging from 23 per cent at 10^11 Msol/h to 16 per cent at 10^14 Msol/h. The predicted concentrations are much lower than inferred from X-ray observations of groups and clusters.

The cosmic star formation rate is observed to drop sharply after redshift z=2. We use a large, cosmological, smoothed particle hydrodynamics simulation to investigate how this decline is related to the evolution of gas accretion and to outflows driven by active galactic nuclei (AGN). We find that the drop in the star formation rate follows a corresponding decline in the global cold-mode accretion rate density onto haloes, but with a delay of order the gas consumption time scale in the interstellar medium. Here we define cold-mode (hot-mode) accretion as gas that is accreted and whose temperature has never exceeded (did exceed) 10^5.5 K. In contrast to cold-mode accretion, which peaks at z~3, the hot mode continues to increase to z~1 and remains roughly constant thereafter. By the present time, the hot mode strongly dominates the global accretion rate onto haloes. Star formation does not track hot-mode halo accretion because most of the hot halo gas never accretes onto galaxies. AGN feedback plays a crucial role by preferentially preventing gas that entered haloes in the hot mode from accreting onto their central galaxies. Consequently, in the absence of AGN feedback, gas accreted in the hot mode would become the dominant source of fuel for star formation and the drop off in the cosmic star formation rate would be much less steep.

(Abridged) We study the rate at which gas accretes onto galaxies and haloes and investigate whether the accreted gas was shocked to high temperatures before reaching a galaxy. For this purpose we use a suite of large cosmological, hydrodynamical simulations from the OWLS project. We improve on previous work by considering a wider range of halo masses and redshifts, by distinguishing accretion onto haloes and galaxies, by including important feedback processes, and by comparing simulations with different physics. The specific rate of gas accretion onto haloes is, like that for dark matter, only weakly dependent on halo mass. For halo masses Mhalo>>10^11 Msun it is relatively insensitive to feedback processes. In contrast, accretion rates onto galaxies are determined by radiative cooling and by outflows driven by supernovae and active galactic nuclei. Galactic winds increase the halo mass at which the central galaxies grow the fastest by about two orders of magnitude to Mhalo~10^12 Msun. Gas accretion is bimodal, with maximum past temperatures either of order the virial temperature or <~10^5 K. The fraction of gas accreted on to haloes in the hot mode is insensitive to feedback and metal-line cooling. It increases with decreasing redshift, but is mostly determined by halo mass, increasing gradually from less than 10% for ~10^11 Msun to greater than 90% at 10^13 Msun. In contrast, for accretion onto galaxies the cold mode is always significant and the relative contributions of the two accretion modes are more sensitive to feedback and metal-line cooling. The majority of stars present in any mass halo at any redshift were formed from gas accreted in the cold mode, although the hot mode contributes typically over 10% for Mhalo>~10^11 Msun. Galaxies, but not necessarily their gaseous haloes, are predominantly fed by gas that did not experience an accretion shock when it entered the host halo.

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.

Galaxy groups are not scaled down versions of massive galaxy clusters - the hot gas in groups (known as the intragroup medium, IGrM for short) is, on average, less dense than the intracluster medium, implying that one or more non-gravitational processes (e.g., radiative cooling, star formation, and/or feedback) has had a relatively larger effect on groups. In the present study, we compare a number of cosmological hydrodynamic simulations that form part of the OverWhelmingly Large Simulations project to isolate and quantify the effects of cooling and feedback from supernovae (SNe) and active galactic nuclei (AGN) on the gas. This is achieved by comparing Lagrangian thermal histories of the gas in the different runs, which were all started from identical initial conditions. While radiative cooling, star formation, and SN feedback are all necessary ingredients, only runs that also include AGN feedback are able to successfully reproduce the optical and X-ray properties of groups and low-mass clusters. We isolate how, when, and exactly what gas is heated by AGN. Interestingly, we find that the gas that constitutes the present-day IGrM is that which was not strongly heated by AGN. Instead, the low median density/high median entropy of the gas in present-day groups is achieved by the ejection of lower entropy gas from low-mass progenitor galaxies at high redshift (primarily 2 < z < 4). This corresponds to the epoch when supermassive black holes accreted most of their mass, typically at a rate that is close to the Eddington limit (i.e., when the black holes are in a `quasar mode').

We study the properties of simulated high-redshift galaxies using cosmological N-body/gasdynamical runs from the OverWhelmingly Large Simulations (OWLS) project. The runs contrast several feedback implementations of varying effectiveness: from no-feedback, to supernova-driven winds to powerful AGN-driven outflows. These different feedback models result in large variations in the abundance and structural properties of bright galaxies at z=2. We find that feedback affects the baryonic mass of a galaxy much more severely than its spin, which is on average roughly half that of its surrounding dark matter halo in our runs. Feedback induces strong correlations between angular momentum content and galaxy mass that leave their imprint on galaxy scaling relations and morphologies. Encouragingly, we find that galaxy disks are common in moderate-feedback runs, making up typically ~50% of all galaxies at the centers of haloes with virial mass exceeding 1e11 M_sun. The size, stellar masses, and circular speeds of simulated galaxies formed in such runs have properties that straddle those of large star-forming disks and of compact early-type galaxies at z=2. Once the detailed abundance and structural properties of these rare objects are well established it may be possible to use them to gauge the overall efficacy of feedback in the formation of high redshift galaxies.

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.

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)

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.

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

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

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.

Galaxy formation models typically assume that the size and rotation speed of galaxy disks are largely dictated by the mass, concentration, and spin of their surrounding dark matter haloes. Equally important, however, are the fraction of baryons in the halo that collect into the central galaxy, as well as the net angular momentum that they are able to retain during its assembly process. We explore the latter using a set of four large cosmological N-body/gasdynamical simulations drawn from the OWLS (OverWhelmingly Large Simulations) project. These runs differ only in their implementation of feedback from supernovae. We find that, when expressed as fractions of their virial values, galaxy mass and net angular momentum are tightly correlated. Galaxy mass fractions, m_d=M_gal/M_vir, depend strongly on feedback, but only weakly on halo mass or spin over the halo mass range explored here (M_vir>1e11 h^{-1}M_sun). The angular momentum of a galaxy, j_d=J_gal/J_vir, correlates with m_d in a manner that is insensitive to feedback and that deviates strongly from the simple j_d = m_d assumption often adopted in semi-analytic models of galaxy formation. The m_d-j_d correlation implies that, in a given halo, galaxy disk size is maximal when the central galaxy makes up a substantial fraction (~20%-30%) of all baryons within the virial radius. At z=2, such systems may host gaseous disks with radial scale lengths as large as those reported for star-forming disks by the SINS survey, even in moderately massive haloes of average spin. Extended disks at z=2 may thus signal the presence of systems where galaxy formation has been particularly efficient, rather than the existence of haloes with unusually high spin parameter.

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

We use the OWLS (OverWhelmingly Large Simulations) set of cosmological Nbody/gasdynamical simulations to study the properties of simulated galaxies at z=2. We focus on the effect of supernova feedback from evolving stars on the baryonic mass and angular momentum content of galaxies that assemble at the center of 1e11-3e12 h^{-1}M_\odot halos. Our main finding is that the mass and angular momentum of such galaxies are strongly coupled, in a way that is approximately independent of feedback: varying the feedback implementation leads, in a given halo, to large variations in galaxy mass but leaves the galaxy mass-angular momentum correlation largely unaltered. In particular, the ratio between the angular momentum of a galaxy and that of its surrounding halo (j_d=J_gal/J_vir) correlates closely with the galaxy mass (expressed in units of the virial mass of the halo; m_d=M_gal/M_vir). This correlation differs substantially from the m_d=j_d assumption commonly adopted in semianalytic models of galaxy formation. We use these results to infer the sizes of disk galaxies at z=2 expected in the LCDM scenario and to interpret recent observations of extended disks at z~2 by the SINS collaboration

Feedback from star formation is thought to play a key role in the formation and evolution of galaxies, but its implementation in cosmological simulations is currently hampered by a lack of numerical resolution. We present and test a sub-grid recipe to model feedback from massive stars in cosmological smoothed particle hydrodynamics simulations. The energy is distributed in kinetic form among the gas particles surrounding recently formed stars. The impact of the feedback is studied using a suite of high-resolution simulations of isolated disc galaxies embedded in dark halos with total mass 10^{10} and 10^{12} Msol/h. We focus in particular on the effect of pressure forces on wind particles within the disc, which we turn off temporarily in some of our runs to mimic a recipe that has been widely used in the literature. We find that this popular recipe gives dramatically different results because (ram) pressure forces on expanding superbubbles determine both the structure of the disc and the development of large-scale outflows. Pressure forces exerted by expanding superbubbles puff up the disc, giving the dwarf galaxy an irregular morphology and creating a galactic fountain in the massive galaxy. Hydrodynamic drag within the disc results in a strong increase of the effective mass loading of the wind for the dwarf galaxy, but quenches much of the outflow in the case of the high-mass galaxy.

When averaged over large scales, star formation in galaxies is observed to follow the empirical Kennicutt-Schmidt (KS) law for surface densities above a constant threshold. While the observed law involves surface densities, theoretical models and simulations generally work with volume density laws (i.e. Schmidt laws). We derive analytic relations between star formation laws expressed in terms of surface densities, volume densities, and pressures and we show how these relations depend on parameters such as the effective equation of state of the multiphase interstellar medium. Our analytic relations enable us to implement observed surface density laws into simulations. Because the parameters of our prescription for star formation are observables, we are not free to tune them to match the observations. We test our theoretical framework using high-resolution simulations of isolated disc galaxies that assume an effective equation of state for the multiphase interstellar medium. We are able to reproduce the star formation threshold and both the slope and the normalisation of arbitrary input KS laws without tuning any parameters and with very little scatter, even for unstable galaxies and even if we use poor numerical resolution. Moreover, we can do so for arbitrary effective equations of state. Our prescription therefore enables simulations of galaxies to bypass our current inability to simulate the formation of stars. On the other hand, the fact that we can reproduce arbitrary input thresholds and KS laws, rather than just the particular ones picked out by nature, indicates that simulations that lack the physics and/or resolution to simulate the multiphase interstellar medium can only provide limited insight into the origin of the observed star formation laws.

The intergalactic magnetic field within filaments should be less polluted by magnetised outflows from active galaxies than magnetic fields in clusters. Therefore, filaments may be a better laboratory to study magnetic field amplification by structure formation than galaxy clusters which typically host many more active galaxies. We present highly resolved cosmological AMR simulations of magnetic fields in the cosmos and make predictions about the evolution and structure of magnetic fields in filaments. Comparing our results to observational evidence for magnetic fields in filaments suggests that amplification of seed fields by gravitational collapse is not sufficient to produce IGM fields. Finally, implications for cosmic ray transport are discussed.

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.