John C. Helly

Simulations of the formation, evolution and clustering of galaxies and quasars

The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,1603 particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.Numerical simulations are a primary theoretical tool to study the nonlinear gravitational growth of structure in the Universe, and to link the initial conditions of cold dark matter (CDM) cosmogonies to observations of galaxies at the present day. Without direct numerical simulation, the hierarchical build-up of structure with its threedimensional dynamics would be largely inaccessible. Since the dominant mass component, the dark matter, is assumed to consist of weakly interacting elementary particles that interact only gravitationally, such simulations use a set of discrete point particles to represent the collisionless dark matter fluid. This representation as an N-body system is obviously only a coarse approximation, and im-

Breaking the hierarchy of galaxy formation

Recent observations of the distant Universe suggest that much of the stellar mass of bright galaxies was already in place at z > 1. This presents a challenge for models of galaxy formation because massive halos are assembled late in the hierarchical clustering process intrinsic to the cold dark matter (CDM) cosmology. In this paper, we discuss a new implementation of the Durham semi-analytic model of galaxy formation in which feedback due to active galactic nuclei (AGN) is assumed to quench cooling flows in massive halos. This mechanism naturally creates a break in the local galaxy luminosity function at bright magnitudes. The model is implemented within the Millennium N-body simulation. The accurate dark matter merger trees and large number of realisations of the galaxy formation process enabled by this simulation result in highly accurate statistics. After adjusting the values of the physical parameters in the model by reference to the properties of the local galaxy population, we investigate the evolution of the K-band luminosity and galaxy stellar mass functions. We calculate the volume-averaged star formation rate density of the Universe as a function of redshift and the way in which this is apportioned amongst galaxies of different mass. The model robustly predicts a substantial population of massive galaxies out to redshift z � 5 and a star formation rate density which rises at least out to z � 2 in objects of all masses. Although observational data on these properties have been cited as evidence for “anti-hierarchical” galaxy formation, we find that when AGN feedback is taken into account, the fundamentally hierarchical CDM model provides a very good match to these observations.

The EAGLE project: simulating the evolution and assembly of galaxies and their environments

We introduce the Virgo Consortiums EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite 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 ≲0.2 dex over the full mass range, 108<M∗/M⊙≲1011, 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∗≳109M⊙, 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.

Galactic stellar haloes in the CDM model

We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high-resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function (LF) and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1 < z < 7 from less than five significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals (10^(7)–10^(8) M_⊙). In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses (~10^(8)–10^(9) M_⊙) and density/surface-brightness profiles of the stellar haloes (from 10 to 100 kpc) are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. The central regions (<10 kpc) of our haloes are highly prolate (c/a ~ 0.3), although we find one example of a massive accreted thick disc. Metallicity gradients in our haloes are typically significant only where the halo is built from a small number of satellites. We contrast the ages and metallicities of halo stars with surviving satellites, finding broad agreement with recent observations.

A marked correlation function analysis of halo formation times in the Millennium Simulation

We study the environmental dependence of the formation epoch of dark matter haloes in the Millennium Simulation: a ten billion particle N-body simulation of standard Lambda cold dark matter cosmology. A sensitive test of this dependence – the marked correlation function – reveals highly significant evidence that haloes of a given mass form earlier in denser regions. We define a marked cross-correlation function, which helps quantify how this effect depends upon the choice of the halo population used to define the environment. The mean halo formation redshift as a function of the local overdensity in dark matter is also well determined, and we see an especially clear dependence for galaxy-sized haloes. This contradicts one of the basic predictions of the excursion set model of structure formation, even though we see that this theory predicts other features of the distribution of halo formation epochs rather well. It also invalidates an assumption usually employed in the popular halo, or halo occupation distribution, models of galaxy clustering, namely that the distribution of halo properties is a function of halo mass but not of halo environment.

Generating dark matter halo merger trees

We present a new Monte Carlo algorithm to generate merger trees describing the formation history of dark matter haloes. The algorithm is a modification of the algorithm of Cole et al. used in the galform semi-analytic galaxy formation model. As such, it is based on the Extended Press–Schechter theory and so should be applicable to hierarchical models with a wide range of power spectra and cosmological models. It is tuned to be in accurate agreement with the conditional mass functions found in the analysis of merger trees extracted from the Λ cold dark matter Millennium N-body simulation. We present a comparison of its predictions not only with these conditional mass functions, but also with additional statistics of the Millennium Simulation halo merger histories. In all cases, we find it to be in good agreement with the Millennium Simulation and thus it should prove to be a very useful tool for semi-analytic models of galaxy formation and for modelling hierarchical structure formation in general. We have made our merger tree generation code and code to navigate the trees available at http://star-www.dur.ac.uk/~cole/merger_trees.

The distribution of satellite galaxies: the great pancake

The 11 known satellite galaxies within 250 kpc of the Milky Way lie close to a great circle on the sky. We use high-resolution N-body simulations of galactic dark matter haloes to test if this remarkable property can be understood within the context of the cold dark matter (CDM) cosmology. We construct halo merger trees from the simulations and use a semi-analytic model to follow the formation of satellite galaxies. We find that in all six of our simulations, the 11 brightest satellites are indeed distributed along thin, disc-like structures analogous to that traced by the satellites of the Milky Way. This is in sharp contrast to the overall distributions of dark matter in the halo and of subhaloes within it, which, although triaxial, are not highly aspherical. We find that the spatial distribution of satellites is significantly different from that of the most massive subhaloes but is similar to that of the subset of subhaloes that had the most massive progenitors at earlier times. The elongated disc-like structure delineated by the satellites has its long axis aligned with the major axis of the dark matter halo. We interpret our results as reflecting the preferential infall of satellites along the spines of a few filaments of the cosmic web.

Evolution of galaxy stellar masses and star formation rates in the EAGLE simulations

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

A unified multiwavelength model of galaxy formation

\Lambda

Lightcone mock catalogues from semi-analytic models of galaxy formation – I. Construction and application to the BzK colour selection

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.