Gustavo Yepes

Toward a halo mass function for precision cosmology: The Limits of universality

We measure the mass function of dark matter halos in a large set of collisionless cosmological simulations of flat ΛCDM cosmology and investigate its evolution at -->z 2. Halos are identified as isolated density peaks, and their masses are measured within a series of radii enclosing specific overdensities. We argue that these spherical overdensity masses are more directly linked to cluster observables than masses measured using the friends-of-friends algorithm (FOF), and are therefore preferable for accurate forecasts of halo abundances. Our simulation set allows us to calibrate the mass function at -->z = 0 for virial masses in the range -->1011 h−1 M☉ ≤ M≤ 1015 h−1 M☉ to 5%, improving on previous results by a factor of 2-3. We derive fitting functions for the halo mass function in this mass range for a wide range of overdensities, both at -->z = 0 and earlier epochs. Earlier studies have sought to calibrate a universal mass function, in the sense that the same functional form and parameters can be used for different cosmologies and redshifts when expressed in appropriate variables. In addition to our fitting formulae, our main finding is that the mass function cannot be represented by a universal function at this level or accuracy. The amplitude of the universal function decreases monotonically by 20%-50%, depending on the mass definition, from -->z = 0 to 2.5. We also find evidence for redshift evolution in the overall shape of the mass function.

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.

The large-scale bias of dark matter halos: numerical calibration and model tests

We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flatCDM cosmology. Halos are identified using the spherical over density algorithm, which finds the mass around isolated peaks in the density field such that the m ean density istimes the background. We calibrate fitting functions for the large scale bias that are adaptable to any value ofwe examine. We find a � 6% scatter about our best fit bias relation. Our fitting functi ons couple to the halo mass functions of Tinker et. al. (2008) such that bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellip soidal collapse model of Sheth, Mo, & Tormen (2001), and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking l ength 0.2 are systematically lower than for halos with the canonical � = 200 overdensity by � 10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak- background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but � 20% residuals remain, both at high and low masses. Subject headings:cosmology:theory — methods:numerical — large scale structure of the universe

Haloes gone MAD: The Halo-Finder Comparison Project

We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends, spherical-overdensity and phase-space-based algorithms. We

MultiDark simulations: the story of dark matter halo concentrations and density profiles

Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations.

A Dynamical Classification of the Cosmic Web

In this paper, we propose a new dynamical classification of the cosmic web. Each point in space is classified in one of four possible web types: voids, sheets, filaments and knots. The classification is based on the evaluation of the deformation tensor (i.e. the Hessian of the gravitational potential) on a grid. The classification is based on counting the number of eigenvalues above a certain threshold, λth, at each grid point, where the case of zero, one, two or three such eigenvalues corresponds to void, sheet, filament or a knot grid point. The collection of neighbouring grid points, friends of friends, of the same web type constitutes voids, sheets, filaments and knots as extended web objects. A simple dynamical consideration of the emergence of the web suggests that the threshold should not be null, as in previous implementations of the algorithm. A detailed dynamical analysis would have found different threshold values for the collapse of sheets, filaments and knots. Short of such an analysis a phenomenological approach has been opted for, looking for a single threshold to be determined by analysing numerical simulations. Our cosmic web classification has been applied and tested against a suite of large (dark matter only) cosmological N-body simulations. In particular, the dependence of the volume and mass filling fractions on λth and on the resolution has been calculated for the four web types. We also study the percolation properties of voids and filaments. Our main findings are as follows. (i) Already at λth= 0.1 the resulting web classification reproduces the visual impression of the cosmic web. (ii) Between 0.2 ≲λth≲ 0.4, a system of percolated voids coexists with a net of interconnected filaments. This suggests a reasonable choice for λth as the parameter that defines the cosmic web. (iii) The dynamical nature of the suggested classification provides a robust framework for incorporating environmental information into galaxy formation models, and in particular to semi-analytical models.

Hydrodynamical simulations of galaxy formation: effects of supernova feedback

We numerically simulate some of the most critical physical processes in galaxy formation: The supernova feedback, in conjunction with gasdynamics and gravity, plays a crucial role in determining how galaxies arise within the context of a model for large-scale structure. Our treatment incorporates a multi-phase model of the interstellar medium and includes the effects of cooling, heating and metal enrichment by supernovae, and evaporation of cold clouds. The star formation happens inside the clouds of cold gas, which are produced via thermal instability. We simulate the galaxy formation in standard biased CDM model for a variety of parameters and for several resolutions in the range 2--20

THE VELOCITY FUNCTION IN THE LOCAL ENVIRONMENT FROM ΛCDM AND ΛWDM CONSTRAINED SIMULATIONS

h^{-1}

The halo mass function through the cosmic ages

kpc. In our picture, supernova feedback regulates the evolution of the gas components and star formation. The efficiency of cloud evaporation by supernova strongly influences star formation rates. This feedback results in a steady rate of star formation in large galaxies (mass larger than

A kinematic classification of the cosmic web

2-3x10^{11}\Msun