Adrian Jenkins

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-

Stable clustering, the halo model and non-linear cosmological power spectra

We present the results of a large library of cosmological N-body simulations, using power-law initial spectra.

The Aquarius project: the subhaloes of galactic haloes

We have performed the largest ever particle simulation of a Milky Way sized dark matter halo, and present the most comprehensive convergence study for an individual dark matter halo carried out thus far. We have also simulated a sample of six ultrahighly resolved Milky Way sized haloes, allowing us to estimate the halo-to-halo scatter in substructure statistics. In our largest simulation, we resolve nearly 300 000 gravitationally bound subhaloes within the virialized region of the halo. Simulations of the same object differing in mass resolution by factors of up to 1800 accurately reproduce the largest subhaloes with the same mass, maximum circular velocity and position, and yield good convergence for the abundance and internal properties of dark matter substructures. We detect up to four generations of subhaloes within subhaloes, but contrary to recent claims, we find less substructure in subhaloes than in the main halo when regions of equal mean overdensity are compared. The overall substructure mass fraction is much lower in subhaloes than in the main halo. Extrapolating the main halos subhalo mass spectrum down to an Earth mass, we predict the mass fraction in substructure to be well below 3 per cent within 100 kpc, and to be below 0.1 per cent within the solar circle. The inner density profiles of subhaloes show no sign of converging to a fixed asymptotic slope and are well fitted by gently curving profiles of Einasto form. The mean concentrations of isolated haloes are accurately described by the fitting formula of Neto et al. down to maximum circular velocities of 1.5 km s(-1), an extrapolation over some five orders of magnitude in mass. However, at equal maximum circular velocity, subhaloes are more concentrated than field haloes, with a characteristic density that is typically similar to 2.6 times larger and increases with decreasing distance from halo centre.

The 2dF Galaxy Redshift Survey: power-spectrum analysis of the final data set and cosmological implications

We present a power-spectrum analysis of the final 2dF Galaxy Redshift Survey (2dFGRS), employing a direct Fourier method. The sample used comprises 221 414 galaxies with measured redshifts. We investigate in detail the modelling of the sample selection, improving on previous treatments in a number of respects. A new angular mask is derived, based on revisions to the photometric calibration. The redshift selection function is determined by dividing the survey according to rest-frame colour, and deducing a self-consistent treatment of k-corrections and evolution for each population. The covariance matrix for the power-spectrum estimates is determined using two different approaches to the construction of mock surveys, which are used to demonstrate that the input cosmological model can be correctly recovered. We discuss in detail the possible differences between the galaxy and mass power spectra, and treat these using simulations, analytic models and a hybrid empirical approach. Based on these investigations, we are confident that the 2dFGRS power spectrum can be used to infer the matter content of the universe. On large scales, our estimated power spectrum shows evidence for the ‘baryon oscillations’ that are predicted in cold dark matter (CDM) models. Fitting to a CDM model, assuming a primordial n s = 1 spectrum, h = 0.72 and negligible neutrino mass, the preferred

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.

The inner structure of ΛCDM haloes – III. Universality and asymptotic slopes

We investigate the mass profile of cold dark matter (ACDM) haloes using a suite of numerical simulations spanning five decades in halo mass, from dwarf galaxies to rich galaxy clusters. These haloes typically have a few million particles within the virial radius (r 200 ), allowing robust mass profile estimates down to radii <1 per cent of r 200 . Our analysis confirms the proposal of Navarro, Frenk & White (NFW) that the shape of the ACDM halo mass profiles differs strongly from a power law and depends little on mass. The logarithmic slope of the spherically averaged density profile, as measured by β = -d In p/d In r, decreases monotonically towards the centre and becomes shallower than isothermal (β < 2) inside a characteristic radius, r -2 . The fitting formula proposed by NFW provides a reasonably good approximation to the density and circular velocity profiles of individual haloes; circular velocities typically deviate from NFW best fits by <10 per cent over the radial range that is numerically well resolved. Alternatively, systematic deviations from the NFW best fits are also noticeable. Inside r -2 , the profile of simulated haloes becomes shallower with radius more gradually than predicted and, as a result, NFW fits tend to underestimate the dark matter density in these regions. This discrepancy has been interpreted as indicating a steeply divergent cusp with asymptotic inner slope, β 0 ≡ β(r = 0) ∼ 1.5. Our results suggest a different interpretation. We use the density and enclosed mass at our innermost resolved radii to place strong constraints on β 0 : density cusps as steep as r -1.5 are inconsistent with most of our simulations, although β 0 = 1 is still consistent with our data. Our density profiles show no sign of converging to a well-defined asymptotic inner power law. We propose a simple formula that reproduces the radial dependence of the slope better than the NFW profile, and so may minimize errors when extrapolating our results inward to radii not yet reliably probed by numerical simulations.

The inner structure of ΛCDM haloes - I. A numerical convergence study

We present a comprehensive set of convergence tests which explore the role of various numerical parameters on the equilibrium structure of a simulated dark matter halo. We report results obtained with two independent, state-of-the-art, multi-stepping, parallel N–body codes: PKDGRAV and GADGET. We find that convergent mass profiles can be obtained for suitable choices of the gravitational softening, timestep, force accuracy, initial redshift, and particle number. For softenings chosen so that particle discreteness effects are negligible, convergence in the circular velocity is obtained at radii where the following conditions are satisfied: (i) the timestep is much shorter than the local orbital timescale; (ii) accelerations do not exceed a characteristic acceleration imprinted by the gravitational softening; and (iii) enough particles are enclosed so that the collisional relaxation timescale is longer than the age of the universe. Convergence also requires sufficiently high initial redshift and accurate force computations. Poor spatial, time, or force resolution leads generally to systems with artificially low central density, but may also result in the formation of artificially dense central cusps. We have explored several adaptive time-stepping choices and obtained best results when individual timesteps are chosen according to the local acceleration and the gravitational softening (�ti / (ǫ/ai) 1/2 ), although further experimentation may yield better and more efficient criteria. The most stringent requirement for c is typically that imposed on the particle number by the collisional relaxation criterion, which implies that in order to estimate accurate circular velocities at radii where the density contrast may reach � 10 6 , the region must enclose of order 3000 particles (or more than a few

The subhalo populations of ΛCDM dark haloes

We investigate the subhalo populations of dark matter haloes in the concordance �CDM cosmology. We use a large cosmological simulation and a variety of high resolution resimulations of individual cluster and galaxy haloes to study the systematics of subhalo populations over ranges of 1000 in halo mass and 1000 in the ratio of subhalo to parent halo mass. The subhalo populations of different haloes are not scaled copies of each other, but vary systematically with halo properties. On average, the amount of substructure increases with halo mass. At fixed mass, it decreases with halo concentration and with halo formation redshift. These trends are comparable in size to the scatter in subhalo abundance between similar haloes. Averaged over all haloes of given mass, the abundance of low mass subhaloes per unit parent halo mass is independendent of parent mass. It is very similar to the abundance per unit mass of low mass haloes in the universe as a whole, once differing boundary definitions for subhaloes and haloes are accounted for. The radial distribution of subhaloes within their parent haloes is substantially less centrally concentrated than that of the dark matter. It varies at most weakly with the mass (or concentration) of the parent halo and not at all with subhalo mass. It does depend on the criteria used to define the subhalo population considered. About 90 per cent of present-day subhaloes were accreted after z = 1 and about 70 per cent after z = 0.5. Only about 8 per cent of the total mass of all haloes accreted at z = 1 survives as bound subhaloes at z = 0. For haloes accreted at z = 2, the survival mass fraction is just 2 per cent. Subhaloes seen near the centre of their parent typically were accreted earlier and retain less of their original mass than those seen near the edge. These strong systematics mean that comparison with galaxies in real clusters is only possible if the formation of the luminous component is modelled appropriately.

The diversity and similarity of simulated cold dark matter haloes

We study the mass, velocity dispersion and anisotropy profiles of Lambda cold dark matter (Lambda CDM) haloes using a suite of N-body simulations of unprecedented numerical resolution. The Aquarius Project follows the formation of six different galaxy-sized haloes simulated several times at varying numerical resolution, allowing numerical convergence to be assessed directly. The highest resolution simulation represents a single dark matter halo using 4.4 billion particles, of which 1.1 billion end up within the virial radius. Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The mass profile of Lambda CDM haloes deviates slightly but systematically from the form proposed by Navarro, Frenk & White. The spherically averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is gamma equivalent to - d ln /d ln r less than or similar to 1. Asymptotic inner slopes as steep as the recently claimed proportional to r-1.2 are clearly ruled out. The radial dependence of gamma is well approximated by a power law, gamma proportional to r alpha (the Einasto profile). The shape parameter, alpha, varies slightly but significantly from halo to halo, implying that the mass profiles of Lambda CDM haloes are not strictly universal: different haloes cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, /Sigma 3 proportional to r-1.875. The index here is identical to that of Bertschingers similarity solution for self-similar infall on to a point mass from an otherwise uniform Einstein-de Sitter universe. The origin of this striking behaviour is unclear, but its robustness suggests that it reflects a fundamental structural property of Lambda CDM haloes. Our conclusions are reliable down to radii below 0.4 per cent of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

The Santa Barbara Cluster Comparison Project: A Comparison of Cosmological Hydrodynamics Solutions

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including smoothed particle hydrodynamics (SPH) and grid methods with fixed, deformable, or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gasdynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be nonradiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying, although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of 2. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to within about 10%, and the ratio of specific dark matter kinetic to gas thermal energies agree to within about 5%. Various factors, including differences in the internal timing of the simulations, contribute to the spread in calculated cluster properties. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modeled.