Halos and galaxies in the standard cosmological model: results from the Bolshoi simulation
TL;DR: The first results from the new Bolshoi N-body cosmological LCDM simulation were presented in this paper, where the authors presented accurate approximations for statistics such as the halo mass function, the concentrations for distinct halos and sub-halos, abundance of halos as function of their circular velocity, the abundance and the spatial distribution of subhalos.
Abstract: We present the first results from the new Bolshoi N-body cosmological LCDM simulation that uses cosmological parameters favored by current observations. The Bolshoi simulation was done in a volume 250Mpc on a side using 8billion particles with mass and force resolution adequate to follow subhalos down to a completeness limit of Vcirc=50km/ s circular velocity. Using excellent statistics of halos and subhalos (10M at every moment and 50M over the whole history) we present accurate approximations for statistics such as the halo mass function, the concentrations for distinct halos and subhalos, abundance of halos as function of their circular velocity, the abundance and the spatial distribution of subhalos. We find that at high redshifts the concentration falls to a minimum of about 3.8 and then rises slightly for higher values of halo mass. We find that while the Sheth-Tormen approximation for the mass function of halos found by spherical overdensity is accurate at low redshifts, it over-predicts the abundance of halos by nearly an order of magnitude by z=10. We find that the number of subhalos scales with the circular velocity of the host halo as Vhost**0.5, and that subhalos have nearly the same radial distribution as dark matter particles at radii 0.3-2 times the host halo virial radius. The subhalo velocity function n(>V) behaves as V**(-3). We give normalization of this relation for different masses and redshifts. Finally, we use an abundance-matching procedure to assign r-band luminosities to dark matter halos as a function of halo Vcirc, and find that the luminosity-velocity relation is in remarkably good agreement with the observed Tully-Fisher relation for galaxies in the range 50-200km/s.
Citations
More filters
TL;DR: In this article, the authors review the range of complementary techniques and theoretical tools that allow astronomers to map the cosmic history of star formation, heavy element production, and reionization of the Universe from the cosmic "dark ages" to the present epoch.
Abstract: Over the past two decades, an avalanche of data from multiwavelength imaging and spectroscopic surveys has revolutionized our view of galaxy formation and evolution. Here we review the range of complementary techniques and theoretical tools that allow astronomers to map the cosmic history of star formation, heavy element production, and reionization of the Universe from the cosmic "dark ages" to the present epoch. A consistent picture is emerging, whereby the star-formation rate density peaked approximately 3.5 Gyr after the Big Bang, at z~1.9, and declined exponentially at later times, with an e-folding timescale of 3.9 Gyr. Half of the stellar mass observed today was formed before a redshift z = 1.3. About 25% formed before the peak of the cosmic star-formation rate density, and another 25% formed after z = 0.7. Less than ~1% of today's stars formed during the epoch of reionization. Under the assumption of a universal initial mass function, the global stellar mass density inferred at any epoch matches reasonably well the time integral of all the preceding star-formation activity. The comoving rates of star formation and central black hole accretion follow a similar rise and fall, offering evidence for co-evolution of black holes and their host galaxies. The rise of the mean metallicity of the Universe to about 0.001 solar by z = 6, one Gyr after the Big Bang, appears to have been accompanied by the production of fewer than ten hydrogen Lyman-continuum photons per baryon, a rather tight budget for cosmological reionization.
3,104 citations
Space Telescope Science Institute1, University of California, Santa Cruz2, Johns Hopkins University3, Rutgers University4, Durham University5, University of Nottingham6, Harvard University7, University of Innsbruck8, University of Michigan9, DSM10, University of Edinburgh11, University of Massachusetts Amherst12, California Institute of Technology13, UK Astronomy Technology Centre14, University of California, Irvine15, Swinburne University of Technology16, University of Arizona17, The Catholic University of America18, Goddard Space Flight Center19, Hebrew University of Jerusalem20, University of Victoria21, University of California, Berkeley22, Texas A&M University23, University of Notre Dame24, Carnegie Institution for Science25, Smithsonian Institution26, Yale University27, University of Missouri–Kansas City28, University of California, Riverside29, Max Planck Society30, University of Pittsburgh31, Inter-University Centre for Astronomy and Astrophysics32, University of Barcelona33, European Southern Observatory34, University of Minnesota35, National Research Council36, Western Kentucky University37, Stanford University38, Atacama Large Millimeter Submillimeter Array39, University of Missouri40
TL;DR: The Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) as discussed by the authors was designed to document the first third of galactic evolution, from z approx. 8 - 1.5 to test their accuracy as standard candles for cosmology.
Abstract: The Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, from z approx. 8 - 1.5. It will image > 250,000 distant galaxies using three separate cameras on the Hubble Space Tele8cope, from the mid-UV to near-IR, and will find and measure Type Ia supernovae beyond z > 1.5 to test their accuracy as standard candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10(exp 9) solar mass to z approx. 2, reaching the knee of the UV luminosity function of galaxies to z approx. 8. The survey covers approximately 800 square arc minutes and is divided into two parts. The CANDELS/Deep survey (5(sigma) point-source limit H =27.7mag) covers approx. 125 square arcminutes within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 50(sigma) point-source limit of H ? or approx. = 27.0 mag. Together with the Hubble Ultradeep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are non-proprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design.
2,088 citations
TL;DR: In this article, a robust method to constrain average galaxy star formation rates, star formation histories (SFHs), and the intracluster light (ICL) as a function of halo mass is presented.
Abstract: We present a robust method to constrain average galaxy star formation rates (SFRs), star formation histories (SFHs), and the intracluster light (ICL) as a function of halo mass. Our results are consistent with observed galaxy stellar mass functions, specific star formation rates (SSFRs), and cosmic star formation rates (CSFRs) from z = 0 to z = 8. We consider the effects of a wide range of uncertainties on our results, including those affecting stellar masses, SFRs, and the halo mass function at the heart of our analysis. As they are relevant to our method, we also present new calibrations of the dark matter halo mass function, halo mass accretion histories, and halo-subhalo merger rates out to z = 8. We also provide new compilations of CSFRs and SSFRs; more recent measurements are now consistent with the buildup of the cosmic stellar mass density at all redshifts. Implications of our work include: halos near 1012 M ☉ are the most efficient at forming stars at all redshifts, the baryon conversion efficiency of massive halos drops markedly after z ~ 2.5 (consistent with theories of cold-mode accretion), the ICL for massive galaxies is expected to be significant out to at least z ~ 1-1.5, and dwarf galaxies at low redshifts have higher stellar mass to halo mass ratios than previous expectations and form later than in most theoretical models. Finally, we provide new fitting formulae for SFHs that are more accurate than the standard declining tau model. Our approach places a wide variety of observations relating to the SFH of galaxies into a self-consistent framework based on the modern understanding of structure formation in ΛCDM. Constraints on the stellar mass-halo mass relationship and SFRs are available for download online.
2,085 citations
TL;DR: The Illustris Project as mentioned in this paper is a series of large-scale hydrodynamical simulations of galaxy formation, which includes primordial and metal-line cooling with self-shielding corrections, stellar evolution, stellar feedback, gas recycling, chemical enrichment, supermassive black hole growth, and feedback from active galactic nuclei.
Abstract: We introduce the Illustris Project, a series of large-scale hydrodynamical simulations of galaxy formation. The highest resolution simulation, Illustris-1, covers a volume of (106.5 Mpc)^3, has a dark mass resolution of 6.26 × 10^6 M_⊙, and an initial baryonic matter mass resolution of 1.26 × 10^6 M_⊙. At z = 0 gravitational forces are softened on scales of 710 pc, and the smallest hydrodynamical gas cells have an extent of 48 pc. We follow the dynamical evolution of 2 × 1820^3 resolution elements and in addition passively evolve 1820^3 Monte Carlo tracer particles reaching a total particle count of more than 18 billion. The galaxy formation model includes: primordial and metal-line cooling with self-shielding corrections, stellar evolution, stellar feedback, gas recycling, chemical enrichment, supermassive black hole growth, and feedback from active galactic nuclei. Here we describe the simulation suite, and contrast basic predictions of our model for the present-day galaxy population with observations of the local universe. At z = 0 our simulation volume contains about 40 000 well-resolved galaxies covering a diverse range of morphologies and colours including early-type, late-type and irregular galaxies. The simulation reproduces reasonably well the cosmic star formation rate density, the galaxy luminosity function, and baryon conversion efficiency at z = 0. It also qualitatively captures the impact of galaxy environment on the red fractions of galaxies. The internal velocity structure of selected well-resolved disc galaxies obeys the stellar and baryonic Tully–Fisher relation together with flat circular velocity curves. In the well-resolved regime, the simulation reproduces the observed mix of early-type and late-type galaxies. Our model predicts a halo mass dependent impact of baryonic effects on the halo mass function and the masses of haloes caused by feedback from supernova and active galactic nuclei.
2,012 citations
TL;DR: In this paper, the authors review the range of complementary techniques and theoretical tools that allow astronomers to map the cosmic history of star formation, heavy element production, and reionization of the Universe from the cosmic "dark ages" to the present epoch.
Abstract: Over the past two decades, an avalanche of data from multiwavelength imaging and spectroscopic surveys has revolutionized our view of galaxy formation and evolution. Here we review the range of complementary techniques and theoretical tools that allow astronomers to map the cosmic history of star formation, heavy element production, and reionization of the Universe from the cosmic "dark ages" to the present epoch. A consistent picture is emerging, whereby the star-formation rate density peaked approximately 3.5 Gyr after the Big Bang, at z~1.9, and declined exponentially at later times, with an e-folding timescale of 3.9 Gyr. Half of the stellar mass observed today was formed before a redshift z = 1.3. About 25% formed before the peak of the cosmic star-formation rate density, and another 25% formed after z = 0.7. Less than ~1% of today's stars formed during the epoch of reionization. Under the assumption of a universal initial mass function, the global stellar mass density inferred at any epoch matches reasonably well the time integral of all the preceding star-formation activity. The comoving rates of star formation and central black hole accretion follow a similar rise and fall, offering evidence for co-evolution of black holes and their host galaxies. The rise of the mean metallicity of the Universe to about 0.001 solar by z = 6, one Gyr after the Big Bang, appears to have been accompanied by the production of fewer than ten hydrogen Lyman-continuum photons per baryon, a rather tight budget for cosmological reionization.
1,626 citations
References
More filters
TL;DR: In this paper, the authors used a modern compilation of observational data along with the best available large-scale cosmological simulation of dark matter to fit all basic statistics of galaxies with circular velocities Vcirc > 80 km/s.
Abstract: It has long been regarded as difficult 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. We find that the standard cosmological model, used in conjunction with halo abundance matching (HAM) and simple dynamical corrections, fits all basic statistics of galaxies with circular velocities Vcirc > 80 km/s. Our observational constraint is the luminosity-velocity relation which allows all types of galaxies to be included. We have compiled data for a variety of galaxies ranging from dwarf irregulars to giant ellipticals. The data present a clear monotonic luminosity-velocity relation from 50 km/s to 500 km/s, with a bend below 80 km/s and a systematic offset between late- and early-type galaxies. For comparison to theory, we employ our LCDM "Bolshoi" simulation of dark matter, which has unprecedented mass and force resolution. We use halo abundance matching to assign rank-ordered galaxy luminosities to the dark matter halos. The resulting predictions for the luminosity-velocity relation are in excellent agreement with the available data on both early-type and late-type galaxies for the luminosity range from Mr = -14-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 is in agreement with the galaxy velocity function for 80-400 km/s. However, we find that the dark matter halos with Vcirc < 80 km/s are much more abundant than observed galaxies with the same Vcirc . We find that the two-point correlation function of galaxies in our model matches very well the results from the SDSS.
242 citations
TL;DR: In this article, the authors presented a new full sky temperature and polarization maps based on seven years of data from WMAP, which are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures.
Abstract: (Abridged) New full sky temperature and polarization maps based on seven years of data from WMAP are presented. The new results are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures. The improvements are described in detail. The seven year data set is well fit by a minimal six-parameter flat Lambda-CDM model. The parameters for this model, using the WMAP data in conjunction with baryon acoustic oscillation data from the Sloan Digital Sky Survey and priors on H_0 from Hubble Space Telescope observations, are: Omega_bh^2 = 0.02260 +-0.00053, Omega_ch^2 = 0.1123 +-0.0035, Omega_Lambda = 0.728 +0.015 -0.016, n_s = 0.963 +-0.012, tau = 0.087 +-0.014 and sigma_8 = 0.809 +-0.024 (68 % CL uncertainties). The temperature power spectrum signal-to-noise ratio per multipole is greater that unity for multipoles 3 sigma. These new WMAP measurements provide important tests of Big Bang cosmology.
223 citations
TL;DR: In this paper, the authors present an analysis of the distribution of structural properties for Milky Way-mass halos in the Millennium-II Simulation (MS-II) and compare their results to simulations of six halos at a much higher resolution.
Abstract: We present an analysis of the distribution of structural properties for Milky Way-mass halos in the Millennium-II Simulation (MS-II). This simulation of structure formation within the standard LCDM cosmology contains thousands of Milky Way-mass halos and has sufficient resolution to properly resolve many subhalos per host. It thus provides a major improvement in the statistical power available to explore the distribution of internal structure for halos of this mass. In addition, the MS-II contains lower resolution versions of the Aquarius Project halos, allowing us to compare our results to simulations of six halos at a much higher resolution. We study the distributions of mass assembly histories, of subhalo mass functions and accretion times, and of merger and stripping histories for subhalos capable of impacting disks at the centers of halos. We show that subhalo abundances are not well-described by Poisson statistics at low mass, but rather are dominated by intrinsic scatter. Using the masses of subhalos at infall and the abundance-matching assumption, there is less than a 10% chance that a Milky Way halo with M_vir =10^12 M_sun will host two galaxies as bright as the Magellanic Clouds. This probability rises to ~25% for a halo with M_vir=2.5 x 10^12 M_sun. The statistics relevant for disk heating are very sensitive to the mass range that is considered relevant. Mergers with infall mass : redshift zero virial mass greater than 1:30 could well impact a central galactic disk and are a near inevitability since z=2, whereas only half of all halos have had a merger with infall mass : redshift zero virial mass greater than 1:10 over this same period.
118 citations
Posted Content•
TL;DR: This work releases for public use a cosmological PM N-body code, a complete package of routines needed to set initial conditions, to run the code, and to analyze the results, and provides results of test runs.
Abstract: Particle-Mesh (PM) codes are still very useful tools for testing predictions of cosmological models in cases when extra high resolution is not very important. We release for public use a cosmological PM N-body code. We provide a complete package of routines needed to set initial conditions, to run the code, and to analyze the results. The package allows you to simulate models with numerous combinations of parameters: open/flat/closed background, with or without the cosmological constant, different values of the Hubble constant, with or without hot neutrinos, tilted or non-tilted initial spectra, different amount of baryons. Routines are included to measure the power spectrum and the density distribution function in your simulations, and a bound-density-maxima code for halo finding. We also provide results of test runs. A simulation with 256^3 mesh and 128^3 particles can be done in a couple of days on a typical workstation (70Mb of RAM are needed). To run simulations with 800^3 mesh and 256^3 particles one needs a computer with 1Gb memory and 1Gb disk space. The code has been successfully tested on an HP workstation and on a Sun workstation running Solaris. Most of the files (not tests) can be obtained from this ftp URL The package can be downloaded from this http URL We provide this tool as a service to the astronomical community, but we cannot guarantee results or publications.
91 citations
TL;DR: The basic ideas of MPI parallelization of the N-body Adaptive Refinement Tree (ART) code, which uses self-adaptive domain decomposition, is described, which requires little communications between processors and is very efficient for large cosmological simulations.
Abstract: We describe the basic ideas of MPI parallelization of the N-body Adaptive Refinement Tree (ART) code. The code uses self-adaptive domain decomposition where boundaries of the domains (parallelepipeds) constantly move—with many degrees of freedom—in the search of the minimum of CPU time. The actual CPU time spent by each MPI task on previous time-step is used to adjust boundaries for the next time-step. For a typical decomposition of 53 domains, the number of possible changes in boundaries is 384≈1040. We describe two algorithms of finding minimum of CPU time for configurations with a large number of domains. Each MPI task in our code solves the N-body problem where the large-scale distribution of matter outside of the boundaries of a domain is represented by relatively few temporary large particles created by other domains. At the beginning of a zero-level time-step, domains create and exchange large particles. Then each domain advances all its particles for many small time-steps. At the end of the large step, the domains decide where to place new boundaries and re-distribute particles. The scheme requires little communications between processors and is very efficient for large cosmological simulations.
38 citations