Julio F. Navarro

A Universal density profile from hierarchical clustering

We use high-resolution N-body simulations to study the equilibrium density profiles of dark matter halos in hierarchically clustering universes. We find that all such profiles have the same shape, independent of the halo mass, the initial density fluctuation spectrum, and the values of the cosmological parameters. Spherically averaged equilibrium profiles are well fitted over two decades in radius by a simple formula originally proposed to describe the structure of galaxy clusters in a cold dark matter universe. In any particular cosmology, the two scale parameters of the fit, the halo mass and its characteristic density, are strongly correlated. Low-mass halos are significantly denser than more massive systems, a correlation that reflects the higher collapse redshift of small halos. The characteristic density of an equilibrium halo is proportional to the density of the universe at the time it was assembled. A suitable definition of this assembly time allows the same proportionality constant to be used for all the cosmologies that we have tested. We compare our results with previous work on halo density profiles and show that there is good agreement. We also provide a step-by-step analytic procedure, based on the Press-Schechter formalism, that allows accurate equilibrium profiles to be calculated as a function of mass in any hierarchical model.

The Structure of cold dark matter halos

High resolution N-body simulations show that the density profiles of dark matter halos formed in the standard CDM cosmogony can be fit accurately by scaling a simple “universal” profile. Regardless of their mass, halos are nearly isothermal over a large range in radius, but significantly shallower than r -2 near the center and steeper than r -2 in the outer regions. The characteristic overdensity of a halo correlates strongly with halo mass in a manner consistent with the mass dependence of the epoch of halo formation. Matching the shape of the rotation curves of disk galaxies with this halo structure requires (i) disk mass-to-light ratios to increase systematically with luminosity, (ii) halo circular velocities to be systematically lower than the disk rotation speed, and (iii) that the masses of halos surrounding bright galaxies depend only weakly on galaxy luminosity. This offers an attractive explanation for the puzzling lack of correlation between luminosity and dynamics in observed samples of binary galaxies and of satellite companions of bright spiral galaxies, suggesting that the structure of dark matter halos surrounding bright spirals is similar to that of cold dark matter halos.

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-

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.

Simulations of X-ray clusters

We present simulations of the formation and evolution of galaxy clusters in the Cold Dark Matter cosmogony. Clusters with a wide range of mass were selected from previous N-body models, and were resimulated at higher resolution using a combined N-body/Smooth Particle Hydrodynamics code. The effects of radiative cooling on the gas are neglected. While many present-day clusters are predicted to be undergoing mergers, the density profiles of those that are approximately in equilibrium are all very similar, both for the gas and for the dark matter. These profiles show no sign of a uniform density core and steepen gradually from the centre outwards. The standard

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

\beta

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

-model is a reasonable fit over most of the radius range observable in real clusters. However, the value obtained for the slope parameter

A Recipe for galaxy formation

\beta_f

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

increases with the outermost radius of the fit. Temperature profiles of different simulated clusters are also similar. Typically the temperature is almost uniform in the regions which emit most of the X-ray flux but drops at larger radii. The gas temperature and dark matter velocity dispersion in equilibrium clusters give values of

The Origin of Star Formation Gradients in Rich Galaxy Clusters

\beta_T\equiv \mu m_p\sigma_{DM}^2/kT