Search

The one-point probability distribution function (PDF) of the matter density field in the universe is a fundamental property that plays an essential role in cosmology for estimates such as gravitational weak lensing, non-linear clustering, massive production of mock galaxy catalogues, and testing predictions of cosmological models. Here we make a comprehensive analysis of the dark matter PDF, using a suite of ∼7000 N -body simulations that covers a wide range of numerical and cosmological parameters. We find that the PDF has a simple shape: it declines with density as a power-law P ∝ ρ−2, which is exponentially suppressed on both small and large densities. The proposed double-exponential approximation provides an accurate fit to all our N -body results for small filtering scales R < 5 h−1Mpc with rms density fluctuations σ > 1. In combination with the spherical infall model that works well for small fluctuations σ < 1, the PDF is now approximated with just few per cent errors over the range of 12 orders of magnitude – a remarkable example of precision cosmology. We find that at |${\sim } 5{-}10{{\ \rm per\ cent}}$| level the PDF explicitly depends on redshift (at fixed σ) and on cosmological density parameter Ωm. We test different existing analytical approximations and find that the often-used lognormal approximation is always 3–5 times less accurate than either the double-exponential approximation or the spherical infall model. [ABSTRACT FROM AUTHOR]

Dwarf spheroidal galaxies (dSphs) are extremely gas-poor, dark matter-dominated galaxies, which make them ideal to test the predictions of the cold dark matter (CDM) model. We argue that the removal of the baryonic component from gas-rich dwarf irregular galaxies, the progenitors of dSphs, can substantially reduce their central density. Thus, it may play an important role in alleviating one of the problems of the CDM model related to the structure of relatively massive satellite galaxies of the Milky Way (MW). Traditionally, collisionless cosmological N-body simulations are used when confronting theoretical predictions with observations. However, these simulations assume that the baryon fraction everywhere in the Universe is equal to the cosmic mean, an assumption which is incorrect for dSphs. We find that the combination of (i) the lower baryon fraction in dSphs compared to the cosmic mean and (ii) the concentration of baryons in the inner part of the MW halo can go a long way towards explaining the observed circular velocity profiles of dSphs. We perform controlled numerical simulations that mimic the effects of baryons. From these we find that the blowing away of baryons by ram pressure, when the dwarfs fall into larger galaxies, decreases the circular velocity profile of the satellite and reduces the density in the central ∼200–500 pc by a factor of $(1- f_{\rm b})^4\approx 0.5$, where $f_{\rm b}$ is the cosmological fraction of baryons. Additionally, the enhanced baryonic mass in the central regions of the parent galaxy generates tidal forces, which are larger than those experienced by subhaloes in traditional N-body simulations. Increased tidal forces substantially alter circular velocity profiles for satellites with pericentres less than 50 kpc. We show that these two effects are strong enough to bring the predictions of subhaloes from CDM simulations into agreement with the observed structure of MW dSphs, regardless of the details of the baryonic processes. [ABSTRACT FROM AUTHOR]

ABSTRACT We study the concentration of dark matter haloes and its evolution in N-body simulations of the standard Λ cold dark matter (ΛCDM) cosmology. The results presented in this paper are based on four large N-body simulations with ∼10 billion particles each: the Millennium-I and -II, Bolshoi and MultiDark simulations. The MultiDark (or BigBolshoi) simulation is introduced in this paper. This suite of simulations with high mass resolution over a large volume allows us to compute with unprecedented accuracy the concentration over a large range of scales (about six orders of magnitude in mass), which constitutes the state of the art of our current knowledge on this basic property of dark matter haloes in the ΛCDM cosmology. We find that there is consistency among the different simulation data sets, despite the different codes, numerical algorithms and halo/subhalo finders used in our analysis. We confirm a novel feature for halo concentrations at high redshifts: a flattening and upturn with increasing mass. The concentration c( M, z) as a function of mass and the redshift and for different cosmological parameters shows a remarkably complex pattern. However, when expressed in terms of the linear rms fluctuation of the density field σ( M, z), the halo concentration c(σ) shows a nearly universal simple U-shaped behaviour with a minimum at a well-defined scale at σ∼ 0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (σ < 1), the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-σ( M, z) relation can be accurately parametrized and provides an analytical model for the dependence of concentration on halo mass. When applied to galaxy clusters, our estimates of concentrations are substantially larger - by a factor up to 1.5 - than previous results from smaller simulations, and are in much better agreement with results of observations. [ABSTRACT FROM AUTHOR]

It has long been regarded as difficult if not impossible for a cosmological model to account simultaneously for the galaxy luminosity, mass, and velocity distributions. We revisit this issue using a modern compilation of observational data along with the best available large-scale cosmological simulation of dark matter (DM). We find that the standard cosmological model, used in conjunction with halo abundance matching (HAM) and simple dynamical corrections, fits--at least on average--all basic statistics of galaxies with circular velocities Vcirc > 80km s-1 calculated at a radius of ~10 kpc. Our primary observational constraint is the luminosity-velocity (LV) relation--which generalizes the Tully-Fisher and Faber-Jackson relations in allowing all types of galaxies to be included, and provides a fundamental benchmark to be reproduced by any theory of galaxy formation. We have compiled data for a variety of galaxies ranging from dwarf irregulars to giant ellipticals. The data present a clear monotonic LV relation from ~50 km s-1 to ~500 km s-1 with a bend below ~80 km s-1 and a systematic offset between late-and early-type galaxies. For comparison to theory, we employ our new ACDM "Bolshoi" simulation of DM, which has unprecedented mass and force resolution over a large cosmological volume, while using an up-to-date set of cosmological parameters. We use HAM to assign rank-ordered galaxy luminosities to the DM halos, a procedure that automatically fits the empirical luminosity function and provides a predicted LV relation that can be checked against observations. The adiabatic contraction of DM halos in response to the infall of the baryons is included as an optional model ingredient. The resulting predictions for the LV relation are in excellent agreement with the available data on both early-type and late-type galaxies for the luminosity range from Mr = -14 to Mr = -22. We also compare our predictions for the "cold" baryon mass (i.e., stars and cold gas) of galaxies as a function of circular velocity with the available observations, again finding a very good agreement. The predicted circular velocity function (VF) is also in agreement with the galaxy VF from 80 to 400 km s-1 using the HIPASS survey for late-type galaxies and Sloan Digital Sky Survey (SDSS) for early-type galaxies. However, in accord with other recent results, we find that the DM halos with Vcirc < 80 km s-1 are much more abundant than observed galaxies with the same Vcirc. Finally, we find that the two-point correlation function of bright galaxies in our model matches very well the results from the final data release of the SDSS, especially when a small amount of scatter is included in the HAM prescription. [ABSTRACT FROM AUTHOR]

We study the luminosity function, the peculiar velocities and the sizes of voids in the Local Volume (LV) in observational samples of galaxies which contain galaxies down to and to within the distance 4–8 Mpc. When we compare the results with the predictions of the standard cosmological λ cold dark matter (ΛCDM) model, we find that the theory faces a sever problem: it predicts a factor of 10 more dwarf haloes as compared with the observed number of dwarf galaxies. In the LV, we identify voids with sizes ranging from 1 to 4.5 Mpc and compare the observational distribution of void sizes with the voids in very high-resolution simulations of the ΛCDM model with Wilkinson Microwave Anisotropy Probe 1( WMAP1) and WMAP3 parameters. The theoretical void function matches the observations remarkably well only if we use haloes with circular velocities Vc larger than 40–45 km s−1[ ] for models with and [ ] for . We exclude the possibility that in the ΛCDM model haloes with circular velocities can host galaxies as bright as : there are too many small haloes in the ΛCDM model resulting in voids being too small as compared with the observations. The problem is that many of the observed dwarf galaxies have H i rotational velocities below 25 km s−1 that strictly contradicts the ΛCDM predictions. Thus, the ΛCDM model faces the same overabundance problem, which it had with the number of satellites in the Local Group (LG). We also estimate the rms deviations from the Hubble flow σH for galaxies at different distances from the LG and find that in most of our model LV candidates the rms peculiar velocities are consistent with observational values: for distances less than 3 Mpc and for distances less than 8 Mpc. At distances 4 (8) Mpc, the observed overdensities of galaxies are 3.5–5.5 (1.3–1.6) – significantly larger than typically assumed. [ABSTRACT FROM AUTHOR]

Addresses the emptiness of giant diameter voids found in the distribution of bright galaxies. Prediction of small dark matter haloes by cosmological models; Abundance of dwarf galaxies inside voids; Similarity of the structural elements in the large-scale structure of the galactic distribution.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 h−1 Gpc)3 and mass resolution of 1.5 × 109 h−1 M⊙) by applying the semi-analytic models galacticus, sag, and sage to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities – noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite-mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models – in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by ongoing and future photometric and redshift surveys. All model data consisting of a range of galaxy properties – including broad-band SDSS magnitudes – are publicly available. [ABSTRACT FROM AUTHOR]

N-body cosmological simulations are an essential tool to understand the observed distribution of galaxies. We use the MultiDark simulation suite, run with the Planck cosmological parameters, to revisit the mass and velocity functions. At redshift z = 0, the simulations cover four orders of magnitude in halo mass from ~1011M⊙ with 8783 874 distinct haloes and 532 533 subhaloes. The total volume used is ~515 Gpc³, more than eight times larger than in previous studies. We measure and model the halo mass function, its covariance matrix w.r.t halo mass and the large-scale halo bias. With the formalism of the excursion-set mass function, we explicit the tight interconnection between the covariance matrix, bias and halo mass function. We obtain a very accurate (<2 per cent level) model of the distinct halo mass function. We also model the subhalo mass function and its relation to the distinct halo mass function. The set of models obtained provides a complete and precise framework for the description of haloes in the concordance Planck cosmology. Finally, we provide precise analytical fits of the Vmax maximum velocity function up to redshift z < 2.3 to push for the development of halo occupation distribution using Vmax. The data and the analysis code are made publicly available in the Skies and Universes data base. [ABSTRACT FROM AUTHOR]

We present results on the clustering of 282 068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4 < z < 0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ∼500 to ∼90 h−1 Mpc. We compare these estimates with the expectations of the flat Λ cold dark matter (ΛCDM) standard cosmological model with parameters compatible with Wilkinson Microwave Anisotropy Probe 7 data. We use the MultiDark cosmological simulation, one of the largest N-body runs presently available, together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the ΛCDM model gives a reasonable description to the observed correlation functions at z ≈ 0.5, which is remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a ≳10 per cent deviation in the correlation functions for scales ≲ 1 and ∼10–40 h−1 Mpc. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12 per cent of the ‘galaxies’ in the abundance-matched sample are satellites inhabiting central haloes with mass M ≳ 1014 h−1 M⊙. Using the MultiDark simulation, we also study the real-space halo bias b of the matched catalogue finding that b = 2.00 ± 0.07 at large scales, consistent with the one obtained using the measured BOSS-projected correlation function. Furthermore, the linear large-scale bias, defined using the extrapolated linear matter power spectrum, depends on the number density n of the abundance-matched sample as b = −0.048 − (0.594 ± 0.02)log10(n/ h3 Mpc−3). Extrapolating these results to baryon acoustic oscillation scales, we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ∼2–4 per cent dips at ≳ 3σ level for wavenumbers k ≳ 0.1 h Mpc−1 in the linear large-scale bias. [ABSTRACT FROM AUTHOR]

ABSTRACT We make detailed theoretical predictions for the assembly properties of the Local Group (LG) in the standard Λ cold dark matter cosmological model. We use three cosmological N-body dark matter simulations from the Constrained Local Universe Simulations project, which are designed to reproduce the main dynamical features of the matter distribution down to the scale of a few Mpc around the LG. Additionally, we use the results of an unconstrained simulation with a 60 times larger volume to calibrate the influence of cosmic variance. We characterize the mass aggregation history (MAH) for each halo by three characteristic times: the formation, assembly and last major merger times. A major merger is defined by a minimal mass ratio of 10: 1. We find that the three LGs share a similar MAH with formation and last major merger epochs placed on average ≈10-12 Gyr ago. Between 12 and 17 per cent of the haloes in the mass range 5 × 1011 < Mh < 5 × 1012 h−1 M⊙ have a similar MAH. In a set of pairs of haloes within the same mass range, a fraction of 1-3 per cent share similar formation properties as both haloes in the simulated LG. An unsolved question posed by our results is the dynamical origin of the MAH of the LGs. The isolation criteria commonly used to define LG-like haloes in unconstrained simulations do not narrow down the halo population into a set with quiet MAHs, nor does a further constraint to reside in a low-density environment. The quiet MAH of the LGs provides a favourable environment for the formation of disc galaxies like the Milky Way and M31. The timing for the beginning of the last major merger in the Milky Way dark matter halo matches with the gas-rich merger origin for the thick component in the galactic disc. Our results support the view that the specific large- and mid-scale environments around the LG play a critical role in shaping its MAH and hence its baryonic structure at present. [ABSTRACT FROM AUTHOR]

ABSTRACT 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 further introduce a robust (and publicly available) suite of test scenarios that allow halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo-finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo, and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Through a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase-space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high-resolution cosmological volume, we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity and peak of the rotation curve). We further suggest to utilize the peak of the rotation curve, vmax, as a proxy for mass, given the arbitrariness in defining a proper halo edge. [ABSTRACT FROM AUTHOR]

We investigate the origin and evolution of fossil groups in a concordance ΛCDM cosmological simulation. We consider haloes with masses between and , and study the physical mechanisms that lead to the formation of the large gap in magnitude between the brightest and the second most bright group member, which is typical for these fossil systems. Fossil groups are found to have high dark matter concentrations, which we can relate to their early formation time. The large magnitude gaps arise after the groups have built up half of their final mass, due to merging of massive group members. We show that the existence of fossil systems is primarily driven by the relatively early infall of massive satellites, and that we do not find a strong environmental dependence for these systems. In addition, we find tentative evidence for fossil group satellites falling in on orbits with typically lower angular momentum, which might lead to a more efficient merger on to the host. We find a population of groups at higher redshifts that go through a ‘fossil phase’: a stage where they show a large magnitude gap, which is terminated by renewed infall from their environment. [ABSTRACT FROM AUTHOR]

We present a comparison of the properties of galaxies in the most underdense regions of the Universe, where the galaxy number density is less than 10 per cent of the mean density, with galaxies from more typical regions. We have compiled a sample of galaxies in 46 large nearby voids that were identified using the Sloan Digital Sky Survey DR4, which provides the largest coverage of the sky. We study the u− r colour distribution, morphology, specific star formation rate (SFR) and radial number density profiles for a total of 495 galaxies fainter than located inside the voids and compare these properties with a control sample of field galaxies. We show that there is an excess of blue galaxies inside the voids. However, inspecting the properties of blue and red galaxies separately, we find that galaxy properties such as colour distribution, bulge-to-total ratios and concentrations are remarkably similar between the void and overall sample. The void galaxies also show the same specific SFR at fixed colour as the control galaxies. We compare our results with the predictions of cosmological simulations of galaxy formation using the Millennium Run semi-analytic galaxy catalogue. We show that the properties of the simulated galaxies in large voids are in reasonably good agreement with those found in similar environments in the real Universe. To summarize, in spite of the fact that galaxies in voids live in the least dense large-scale environment, this environment makes very little impact on the properties of galaxies. [ABSTRACT FROM AUTHOR]

We present a statistical analysis of voids in the two-degree Field Galaxy Redshift Survey (2dFGRS). In order to detect the voids, we have developed two robust algorithms. We define voids as non-overlapping maximal spheres empty of haloes or galaxies with mass or luminosity above a given value. We search for voids in cosmological N-Body simulations to test the performance of our void finders. We obtain and analyse the void statistics for several volume-limited samples for the North Galactic Pole (NGP) and the South Galactic Pole (SGP) constructed from the 2dFGRS full data release. We find that the results obtained from the NGP and the SGP are statistically compatible. From the results of several statistical tests we conclude that voids are essentially uncorrelated, with at most a mild anticorrelation and that at the 99.5 per cent confidence level there is a dependence of the void number density on redshift. We develop a technique to correct the distortion caused by the fact that we use the redshift as the radial coordinate. We calibrate this technique with mock catalogues and find that the correction might be of some relevance to carry out accurate inferences from void statistics. We study the statistics of the galaxies inside nine nearby voids. We find that galaxies in voids are not randomly distributed: they form structures like filaments. We also obtain the galaxy number density profile in voids. This profile follow a similar but steeper trend to that followed by haloes in voids. [ABSTRACT FROM AUTHOR]