-
Uchuu-$ν^2$GC galaxies and AGN: Cosmic variance forecasts of high-redshift AGN for JWST, Euclid, and LSST
Authors:
Taira Oogi,
Tomoaki Ishiyama,
Francisco Prada,
Manodeep Sinha,
Darren Croton,
Sofía A. Cora,
Eric Jullo,
Anatoly A. Klypin,
Masahiro Nagashima,
J. López Cacheiro,
José Ruedas,
Masakazu A. R. Kobayashi,
Ryu Makiya
Abstract:
Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-bo…
▽ More
Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$ν^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$ν^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST).
△ Less
Submitted 5 August, 2023; v1 submitted 29 July, 2022;
originally announced July 2022.
-
The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
Authors:
Tomoaki Ishiyama,
Francisco Prada,
Anatoly A. Klypin,
Manodeep Sinha,
R. Benton Metcalf,
Eric Jullo,
Bruno Altieri,
Sofía A. Cora,
Darren Croton,
Sylvain de la Torre,
David E. Millán-Calero,
Taira Oogi,
José Ruedas,
Cristian A. Vega-Martínez
Abstract:
We introduce the Uchuu suite of large high-resolution cosmological $N$-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion ($12800^3$) dark matter particles in a box of side-length 2.0 Gpc/h, with particle mass $3.27 \times 10^{8}$ Msun/h. The highest resolution simulation, Shin-Uchuu, consists of 262 billion ($6400^3$) particles in a box of side-length 140 Mpc/h, with…
▽ More
We introduce the Uchuu suite of large high-resolution cosmological $N$-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion ($12800^3$) dark matter particles in a box of side-length 2.0 Gpc/h, with particle mass $3.27 \times 10^{8}$ Msun/h. The highest resolution simulation, Shin-Uchuu, consists of 262 billion ($6400^3$) particles in a box of side-length 140 Mpc/h, with particle mass $8.97 \times 10^{5}$ Msun/h. Combining these simulations we can follow the evolution of dark matter halos and subhalos spanning those hosting dwarf galaxies to massive galaxy clusters across an unprecedented volume. In this first paper, we present basic statistics, dark matter power spectra, and the halo and subhalo mass functions, which demonstrate the wide dynamic range and superb statistics of the Uchuu suite. From an analysis of the evolution of the power spectra we conclude that our simulations remain accurate from the Baryon Acoustic Oscillation scale down to the very small. We also provide parameters of a mass-concentration model, which describes the evolution of halo concentration and reproduces our simulation data to within 5 per cent for halos with masses spanning nearly eight orders of magnitude at redshift 0<z<14. There is an upturn in the mass-concentration relation for the population of all halos and of relaxed halos at z>0.5, whereas no upturn is detected at z<0.5. We make publicly available various $N$-body products as part of Uchuu Data Release 1 on the Skies & Universes site. Future releases will include gravitational lensing maps and mock galaxy, X-ray cluster, and active galactic nuclei catalogues.
△ Less
Submitted 13 July, 2021; v1 submitted 29 July, 2020;
originally announced July 2020.
-
A semi-analytical perspective on massive galaxies at $z\sim0.55$
Authors:
D. Stoppacher,
F. Prada,
A. D. Montero-Dorta,
S. Rodríguez-Torres,
A. Knebe,
G. Favole,
W. Cui,
A. J. Benson,
C. Behrens,
A. A. Klypin
Abstract:
The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMAS…
▽ More
The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMASS photometric selection from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The comparison of the GALACTICUS CMASS-mock with BOSS-CMASS data allows us to explore different aspects of the massive galaxy population at $0.5<z<0.6$, including the galaxy-halo connection and the galaxy clustering. We find good agreement between our modelled galaxies and observations regarding the galaxy-halo connection, but our CMASS-mock over-estimates the clustering amplitude of the 2-point correlation function, due to a smaller number density compared to BOSS, a lack of blue objects, and a small intrinsic scatter in stellar mass at fixed halo mass of $<0.1$ dex. To alleviate this problem, we construct an alternative mock catalogue mimicking the CMASS colour-magnitude distribution by randomly down-sampling the SAM catalogue. This CMASS-mock reproduces the clustering of CMASS galaxies within 1$σ$ and shows some environmental dependency of star formation properties that could be connected to the quenching of star formation and the assembly bias.
△ Less
Submitted 14 February, 2019;
originally announced February 2019.
-
MultiDark-Galaxies: data release and first results
Authors:
Alexander Knebe,
Doris Stoppacher,
Francisco Prada,
Christoph Behrens,
Andrew Benson,
Sofia A. Cora,
Darren J. Croton,
Nelson D. Padilla,
Andrés N. Ruiz,
Manodeep Sinha,
Adam R. H. Stevens,
Cristian A. Vega-Martínez,
Peter Behroozi,
Violeta Gonzalez-Perez,
Stefan Gottlöber,
Anatoly A. Klypin,
Gustavo Yepes,
Harry Enke,
Noam I. Libeskind,
Kristin Riebe,
Matthias Steinmetz
Abstract:
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 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data…
▽ More
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 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) 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 on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.
△ Less
Submitted 23 October, 2017;
originally announced October 2017.
-
nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics
Authors:
Chia-Hsun Chuang,
Cheng Zhao,
Francisco Prada,
Emiliano Munari,
Santiago Avila,
Albert Izard,
Francisco-Shu Kitaura,
Marc Manera,
Pierluigi Monaco,
Steven Murray,
Alexander Knebe,
Claudia G. Scoccola,
Gustavo Yepes,
Juan Garcia-Bellido,
Felipe A. Marin,
Volker Muller,
Ramin Skibba,
Martin Crocce,
Pablo Fosalba,
Stefan Gottlober,
Anatoly A. Klypin,
Chris Power,
Charling Tao,
Victor Turchaninov
Abstract:
We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used w…
▽ More
We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.
△ Less
Submitted 16 June, 2015; v1 submitted 24 December, 2014;
originally announced December 2014.
-
Hunting down systematics in baryon acoustic oscillations after cosmic high noon
Authors:
Francisco Prada,
Claudia G. Scoccola,
Chia-Hsun Chuang,
Gustavo Yepes,
Anatoly A. Klypin,
Francisco-Shu Kitaura,
Stefan Gottlober,
Cheng Zhao
Abstract:
Future dark energy experiments will require better and more accurate theoretical predictions for the baryonic acoustic oscillations (BAO) signature in the spectrum of cosmological perturbations. Here, we use large N-body simulations of the \LambdaCDM Planck cosmology to study any possible systematic shifts and damping in BAO due to the impact of nonlinear gravitational growth of structure, scale d…
▽ More
Future dark energy experiments will require better and more accurate theoretical predictions for the baryonic acoustic oscillations (BAO) signature in the spectrum of cosmological perturbations. Here, we use large N-body simulations of the \LambdaCDM Planck cosmology to study any possible systematic shifts and damping in BAO due to the impact of nonlinear gravitational growth of structure, scale dependent and non-local bias, and redshift-space distortions. The effect of cosmic variance is largely reduced by dividing the tracer power spectrum by that from a BAO-free simulation starting with the same phases. This permits us to study with unprecedented accuracy (better than 0.02% for dark matter and 0.07% for low-bias halos) small shifts of the pristine BAO wavenumbers towards larger k, and non-linear damping of BAO wiggles in the power spectrum of dark matter and halo populations in the redshift range z=0-1. For dark matter, we provide an accurate parametrization of the evolution of αas a function of the linear growth factor D(z). For halo samples, with bias ranging from 1.2 to 2.8, we measure a typical BAO shift of ~0.25%, observed in real-space, which does not show an appreciable evolution with redshift within the uncertainties. Moreover, we report a constant shift as a function of halo bias. We find a different evolution of the damping of the acoustic feature in all halo samples as compared to dark matter with haloes suffering less damping, and also find some weak dependence on bias. A larger BAO shift and damping is measured in redshift-space which can be well explained by linear theory due to redshift-space distortions. A clear modulation in phase with the acoustic scale is observed in the scale-dependent halo bias due to the presence of the baryonic acoustic oscillations.
△ Less
Submitted 14 November, 2014; v1 submitted 17 October, 2014;
originally announced October 2014.
-
Gravitationally Consistent Halo Catalogs and Merger Trees for Precision Cosmology
Authors:
Peter S. Behroozi,
Risa H. Wechsler,
Hao-Yi Wu,
Michael T. Busha,
Anatoly A. Klypin,
Joel R. Primack
Abstract:
We present a new algorithm for generating merger trees and halo catalogs which explicitly ensures consistency of halo properties (mass, position, and velocity) across timesteps. Our algorithm has demonstrated the ability to improve both the completeness (through detecting and inserting otherwise missing halos) and purity (through detecting and removing spurious objects) of both merger trees and ha…
▽ More
We present a new algorithm for generating merger trees and halo catalogs which explicitly ensures consistency of halo properties (mass, position, and velocity) across timesteps. Our algorithm has demonstrated the ability to improve both the completeness (through detecting and inserting otherwise missing halos) and purity (through detecting and removing spurious objects) of both merger trees and halo catalogs. In addition, our method is able to robustly measure the self-consistency of halo finders; it is the first to directly measure the uncertainties in halo positions, halo velocities, and the halo mass function for a given halo finder based on consistency between snapshots in cosmological simulations. We use this algorithm to generate merger trees for two large simulations (Bolshoi and Consuelo) and evaluate two halo finders (ROCKSTAR and BDM). We find that both the ROCKSTAR and BDM halo finders track halos extremely well; in both, the number of halos which do not have physically consistent progenitors is at the 1-2% level across all halo masses. Our code is publicly available at http://code.google.com/p/consistent-trees . Our trees and catalogs are publicly available at http://hipacc.ucsc.edu/Bolshoi/ .
△ Less
Submitted 11 January, 2013; v1 submitted 19 October, 2011;
originally announced October 2011.
-
Halo Contraction Effect in Hydrodynamic Simulations of Galaxy Formation
Authors:
Oleg Y. Gnedin,
Daniel Ceverino,
Nickolay Y. Gnedin,
Anatoly A. Klypin,
Andrey V. Kravtsov,
Robyn Levine,
Daisuke Nagai,
Gustavo Yepes
Abstract:
The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the h…
▽ More
The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.
△ Less
Submitted 29 August, 2011;
originally announced August 2011.
-
Halo concentrations in the standard LCDM cosmology
Authors:
Francisco Prada,
Anatoly A. Klypin,
Antonio J. Cuesta,
Juan E. Betancort-Rijo,
Joel Primack
Abstract:
We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 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…
▽ More
We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 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 halos in the LCDM cosmology. We find that there is consistency among the different simulation data sets. 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 sigma(M,z), the halo concentration c(sigma) shows a nearly-universal simple U-shaped behaviour with a minimum at a well defined scale at sigma=0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (sigma < 1) the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-sigma(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. (abridged)
△ Less
Submitted 27 April, 2011;
originally announced April 2011.
-
Statistics of Satellite Galaxies Around Milky Way-Like Hosts
Authors:
Michael T. Busha,
Risa H. Wechsler,
Peter S. Behroozi,
Brian F. Gerke,
Anatoly A. Klypin,
Joel R. Primack
Abstract:
We calculate the probability that a Milky-Way-like halo in the standard cosmological model has the observed number of Magellanic Clouds (MCs). The statistics of the number of MCs in the LCDM model are in good agreement with observations of a large sample of SDSS galaxies. Under the sub-halo abundance matching assumption of a relationship with small scatter between galaxy r-band luminosities and ha…
▽ More
We calculate the probability that a Milky-Way-like halo in the standard cosmological model has the observed number of Magellanic Clouds (MCs). The statistics of the number of MCs in the LCDM model are in good agreement with observations of a large sample of SDSS galaxies. Under the sub-halo abundance matching assumption of a relationship with small scatter between galaxy r-band luminosities and halo internal velocities v_max, we make detailed comparisons to similar measurements using SDSS DR7 data by Liu et al. (2010). Models and observational data give very similar probabilities for having zero, one, and two MC-like satellites. In both cases, Milky Way-luminosity hosts have just a \sim 10% chance of hosting two satellites similar to the Magellanic Clouds. In addition, we present a prediction for the probability for a host galaxy to have Nsats satellite galaxies as a function of the magnitudes of both the host and satellite. This probability and its scaling with host properties is significantly different from that of mass-selected objects because of scatter in the mass- luminosity relation and because of variations in the star formation efficiency with halo mass.
△ Less
Submitted 22 November, 2011; v1 submitted 29 November, 2010;
originally announced November 2010.
-
The dark outside: the density profile of dark matter haloes beyond the virial radius
Authors:
Hugo Tavio,
Antonio J. Cuesta,
Francisco Prada,
Anatoly A. Klypin,
Miguel A. Sanchez-Conde
Abstract:
We present an approximation for the average density profile of dark matter haloes in the LambdaCDM cosmological model, which is accurate to within 10--15% even for large radial distances from 0.05R_{vir} up to 10R_{vir} for halo masses ranging from 10^{11.5} to 10^{15.0} h^{-1}M_{Sun}. We propose a modified form of the Navarro, Frenk & White (NFW) approximation: ρ(r)=ρ_{NFW}(r)+ A (r/R_{vir})^{-…
▽ More
We present an approximation for the average density profile of dark matter haloes in the LambdaCDM cosmological model, which is accurate to within 10--15% even for large radial distances from 0.05R_{vir} up to 10R_{vir} for halo masses ranging from 10^{11.5} to 10^{15.0} h^{-1}M_{Sun}. We propose a modified form of the Navarro, Frenk & White (NFW) approximation: ρ(r)=ρ_{NFW}(r)+ A (r/R_{vir})^{-1}+B (1+r/R_{vir})^{-1}. This generalized expression, which is applicable to the external regions of dark matter haloes, only very slightly affects the density in the inner regions of haloes. The strong correlation among the different parameters in the model allows us to describe the profile in terms of just one parameter: the virial mass. We integrate our density profile to derive the enclosed mass in a sphere of a given radius and compare it with the NFW results. We find that the NFW underestimates the enclosed mass by more than 50% at 10R_{vir}, whereas our model reproduces the results from numerical simulations to within 2% accuracy even at this distance. We also use this new approximation to study the weak gravitational lensing and to obtain an analytic expression for the tangential shear. This allows us to quantify the contribution to the shear from the outer regions of the density profile. For the first time we calculate the difference between the tangential shear calculated via the NFW profile and the corresponding result when the external regions of haloes in cosmological simulations are taken into account. We find a 4% difference for all the mass ranges under study.
△ Less
Submitted 18 July, 2008;
originally announced July 2008.
-
Connecting Galaxy Disk and Extended Halo Gas Kinematics
Authors:
G. G. Kacprzak,
C. W. Churchill,
C. C. Steidel,
D. Ceverino,
A. A. Klypin,
M. T. Murphy
Abstract:
We have explored the galaxy disk/extended halo gas kinematic relationship using rotation curves (Keck/ESI) of ten intermediate redshift galaxies which were selected by MgII halo gas absorption observed in quasar spectra. Previous results of six edge-on galaxies, probed along their major axis, suggest that observed halo gas velocities are consistent with extended disk-like halo rotation at galact…
▽ More
We have explored the galaxy disk/extended halo gas kinematic relationship using rotation curves (Keck/ESI) of ten intermediate redshift galaxies which were selected by MgII halo gas absorption observed in quasar spectra. Previous results of six edge-on galaxies, probed along their major axis, suggest that observed halo gas velocities are consistent with extended disk-like halo rotation at galactocentric distances of 25-72 kpc. Using our new sample, we demonstrate that the gas velocities are by and large not consistent with being directly coupled to the galaxy kinematics. Thus, mechanisms other than co-rotation dynamics (i.e., gas inflow, feedback, galaxy-galaxy interactions, etc.) must be invoked to account for the overall observed kinematics of the halo gas. In order to better understand the dynamic interaction of the galaxy/halo/cosmic web environment, we performed similar mock observations of galaxies and gaseous halos in Lambda-CDM cosmological simulations. We discuss an example case of a z=0.92 galaxy with various orientations probing halo gas at a range of positions. The gas dynamics inferred using simulated quasar absorption lines are consistent with observational data.
△ Less
Submitted 25 October, 2007;
originally announced October 2007.
-
Voids in the Local Volume: a limit on appearance of a galaxy in a DM halo
Authors:
Anton V. Tikhonov,
Anatoly A. Klypin
Abstract:
Current explanation of the overabundance of dark matter subhalos in the Local Group (LG) indicates that there maybe a limit on mass of a halo, which can host a galaxy. This idea can be tested using voids in the distribution of galaxies: at some level small voids should not contain any (even dwarf) galaxies. We use observational samples complete to M_B = -12 with distances less than 8 Mpc to cons…
▽ More
Current explanation of the overabundance of dark matter subhalos in the Local Group (LG) indicates that there maybe a limit on mass of a halo, which can host a galaxy. This idea can be tested using voids in the distribution of galaxies: at some level small voids should not contain any (even dwarf) galaxies. We use observational samples complete to M_B = -12 with distances less than 8 Mpc to construct the void function (VF): the distribution of sizes of voids empty of any galaxies. There are ~30 voids with sizes ranging from 1 to 5 Mpc. We then study the distribution of dark matter halos in very high resolution simulations of the LCDM model. The theoretical VF matches the observations remarkably well only if we use halos with circular velocities larger than 45 +/- 10 km/s. This agrees with the Local Group predictions. There are smaller halos in the voids, but they should not produce any luminous matter. Small voids look quite similar to their giant cousins: the density has a minimum at the center of a void and it increases as we get closer to the border. Small nonluminous halos inside the void form a web of tiny filaments. Thus, both the Local Group data and the nearby voids indicate that isolated halos below 45 +/- 10 km/s must not host galaxies and that small (few Mpc) voids are truly dark.
△ Less
Submitted 24 July, 2007;
originally announced July 2007.
-
A Large Dark Matter Core in the Fornax Dwarf Spheroidal Galaxy?
Authors:
Louis E. Strigari,
James S. Bullock,
Manoj Kaplinghat,
Andrey V. Kravtsov,
Oleg Y. Gnedin,
Kevork Abazajian,
Anatoly A. Klypin
Abstract:
We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits…
▽ More
We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.
△ Less
Submitted 28 March, 2006;
originally announced March 2006.
-
Constraining the Projected Radial Distribution of Galactic Satellites with the Sloan Digital Sky Survey
Authors:
Jacqueline Chen,
Andrey V. Kravtsov,
Francisco Prada,
Erin S. Sheldon,
Anatoly A. Klypin,
Michael R. Blanton,
Jonathan Brinkmann,
Aniruddha R. Thakar
Abstract:
We use the Sloan Digital Sky Survey (SDSS) spectroscopic sample to constrain the projected radial distribution of satellites around isolated ~ L* galaxies. We employ mock galaxy catalogs derived from high-resolution cosmological simulations to investigate the effects of interloper contamination and show that interlopers significantly bias the estimated slope of the projected radial distribution…
▽ More
We use the Sloan Digital Sky Survey (SDSS) spectroscopic sample to constrain the projected radial distribution of satellites around isolated ~ L* galaxies. We employ mock galaxy catalogs derived from high-resolution cosmological simulations to investigate the effects of interloper contamination and show that interlopers significantly bias the estimated slope of the projected radial distribution of satellites. We also show that the distribution of interlopers around galaxies is expected to be non-uniform in velocity space because galaxies are clustered and reside in crowded environments. Successful methods of interloper contamination correction should therefore take into account environments of the host galaxies. Two such new methods are presented and the most reliable of them is used to correct for interloper contamination in analyses of the SDSS galaxy sample. The best fit power-law slope of the interloper-corrected surface density distribution of satellites, Sigma(R) ~ R^alpha, in the volume-limited SDSS sample is alpha ~ -1.7 +/- 0.1, independent of the galaxy and satellite luminosities. Comparison with Lambda-CDM simulations shows that the radial distribution of the SDSS satellites is more concentrated than that of subhalos around galaxy-sized halos, especially at R < 100 kpc/h. The predicted dark matter radial distribution is somewhat more concentrated than the profile of the SDSS satellites, but the difference is not statistically significant for our sample.
△ Less
Submitted 22 March, 2006; v1 submitted 14 December, 2005;
originally announced December 2005.
-
How far do they go? The outer structure of dark matter halos
Authors:
F. Prada,
A. A. Klypin,
E. Simonneau,
J. Betancort-Rijo,
S. Patiri,
S. Gottloeber,
M. A. Sanchez-Conde
Abstract:
We study the density profiles of collapsed galaxy-size dark matter halos with masses 1e11-5e12 Msun focusing mostly on the halo outer regions from the formal virial radius Rvir up to 5-7Rvir. We find that isolated halos in this mass range extend well beyond Rvir exhibiting all properties of virialized objects up to 2-3Rvir: relatively smooth density profiles and no systematic infall velocities.…
▽ More
We study the density profiles of collapsed galaxy-size dark matter halos with masses 1e11-5e12 Msun focusing mostly on the halo outer regions from the formal virial radius Rvir up to 5-7Rvir. We find that isolated halos in this mass range extend well beyond Rvir exhibiting all properties of virialized objects up to 2-3Rvir: relatively smooth density profiles and no systematic infall velocities. The dark matter halos in this mass range do not grow as one naively may expect through a steady accretion of satellites, i.e., on average there is no mass infall. This is strikingly different from more massive halos, which have large infall velocities outside of the virial radius. We provide accurate fit for the density profile of these galaxy-size halos. For a wide range (0.01-2)Rvir of radii the halo density profiles are fit with the approximation rho=rho_s exp(-2n[x^{1/n}-1])+rho_m, where x=r/r_s, rho_m is the mean matter density of the Universe, and the index n is in the range n=6-7.5. These profiles do not show a sudden change of behavior beyond the virial radius. For larger radii we combine the statistics of the initial fluctuations with the spherical collapse model to obtain predictions for the mean and most probable density profiles for halos of several masses. The model give excellent results beyond 2-3 formal virial radii.
△ Less
Submitted 17 June, 2005;
originally announced June 2005.
-
The Anisotropic Distribution of Galactic Satellites
Authors:
Andrew R. Zentner,
Andrey V. Kravtsov,
Oleg Y. Gnedin,
Anatoly A. Klypin
Abstract:
We present a study of the spatial distribution of subhalos in galactic dark matter halos using dissipationless cosmological simulations of the concordance LCDM model. We find that subhalos are distributed anisotropically and are preferentially located along the major axes of the triaxial mass distributions of their hosts. The Kolmogorov-Smirnov probability for drawing our simulated subhalo sampl…
▽ More
We present a study of the spatial distribution of subhalos in galactic dark matter halos using dissipationless cosmological simulations of the concordance LCDM model. We find that subhalos are distributed anisotropically and are preferentially located along the major axes of the triaxial mass distributions of their hosts. The Kolmogorov-Smirnov probability for drawing our simulated subhalo sample from an isotropic distribution is P_KS \simeq 1.5 \times 10^{-4}. An isotropic distribution of subhalos is thus not the correct null hypothesis for testing the CDM paradigm. The nearly planar distribution of observed Milky Way (MW) satellites is marginally consistent (probability \simeq 0.02) with being drawn randomly from the subhalo distribution in our simulations. Furthermore, if we select the subhalos likely to be luminous, we find a distribution that is consistent with the observed MW satellites. In fact, we show that subsamples of the subhalo population with a centrally-concentrated radial distribution, similar to that of the MW dwarfs, typically exhibit a comparable degree of planarity. We explore the origin of the observed subhalo anisotropy and conclude that it is likely due to (1) preferential accretion of subhalos along filaments, often closely aligned with the major axis of the host halo, and (2) evolution of satellite orbits within the prolate, triaxial potentials typical of CDM halos. Agreement between predictions and observations requires the major axis of the outer dark matter halo of the Milky Way to be nearly perpendicular to the disk. We discuss possible observational tests of such disk-halo alignment with current large galaxy surveys.
△ Less
Submitted 23 April, 2005; v1 submitted 24 February, 2005;
originally announced February 2005.
-
Response of dark matter halos to condensation of baryons: cosmological simulations and improved adiabatic contraction model
Authors:
Oleg Y. Gnedin,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Daisuke Nagai
Abstract:
The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. In contrast, halos in the hierarchical structure formation scenarios grow via multiple violent me…
▽ More
The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. In contrast, halos in the hierarchical structure formation scenarios grow via multiple violent mergers and accretion along filaments, and particle orbits in the halos are highly eccentric. We study the effects of the cooling of gas in the inner regions of halos using high-resolution cosmological simulations which include gas dynamics, radiative cooling, and star formation. We find that the dissipation of gas indeed increases the density of dark matter and steepens its radial profile in the inner regions of halos compared to the case without cooling. For the first time, we test the adiabatic contraction model in cosmological simulations and find that the standard model systematically overpredicts the increase of dark matter density in the inner 5% of the virial radius. We show that the model can be improved by a simple modification of the assumed invariant from M(r)r to M(r_av)r, where r and r_av are the current and orbit-averaged particle positions. This modification approximately accounts for orbital eccentricities of particles and reproduces simulation profiles to within 10-20%. We present analytical fitting functions that accurately describe the transformation of the dark matter profile in the modified model and can be used for interpretation of observations.
△ Less
Submitted 15 December, 2004; v1 submitted 9 June, 2004;
originally announced June 2004.
-
The tumultuous lives of galactic dwarfs and the missing satellites problem
Authors:
Andrey V. Kravtsov,
Oleg Y. Gnedin,
Anatoly A. Klypin
Abstract:
Hierarchical Cold Dark Matter (CDM) models predict that Milky Way (MW) sized halos contain hundreds of dense low-mass dark satellites, an order of magnitude more than the number of observed satellites in the Local Group (LG). If the CDM paradigm is correct, we need to understand why most of these halos failed to form stars and become galaxies. We analyze the dynamical evolution of the satellite…
▽ More
Hierarchical Cold Dark Matter (CDM) models predict that Milky Way (MW) sized halos contain hundreds of dense low-mass dark satellites, an order of magnitude more than the number of observed satellites in the Local Group (LG). If the CDM paradigm is correct, we need to understand why most of these halos failed to form stars and become galaxies. We analyze the dynamical evolution of the satellite halos in a high-resolution cosmological simulation of MW sized halos in the LCDM cosmology. We find that about 10% of the substructure halos with the present masses <10^8-10^9 Msun (circular velocities Vmax<30 km/s) had considerably larger masses and circular velocities when they formed at redshifts z>2. After the initial period of mass accretion in isolation, these objects experience dramatic mass loss due to tidal stripping, in some cases even before they are accreted by their host halo. This can explain how the smallest dwarf spheroidal galaxies of the LG were able to build up a sizable stellar mass in their seemingly shallow potential wells. We propose a new model in which all of the luminous dwarf spheroidals in the Local Group are descendants of the relatively massive (>10^9 Msun) high-redshift systems, in which the gas could cool efficiently by atomic line emission and which were not significantly affected by the extragalactic ultraviolet radiation. We present a simple galaxy formation model based on the trajectories extracted from the simulation, which accounts for the bursts of star formation after strong tidal shocks and the inefficiency of gas cooling in halos with virial temperatures Tvir<~10^4 K. Our model reproduces the abundance, spatial distribution, and morphological segregation of the observed Galactic satellites. The results are insensitive to the redshift of reionization.
△ Less
Submitted 14 May, 2004; v1 submitted 7 January, 2004;
originally announced January 2004.
-
Density profiles of LCDM clusters
Authors:
Argyro Tasitsiomi,
Andrey V. Kravtsov,
Stefan Gottloeber,
Anatoly A. Klypin
Abstract:
We analyze the mass accretion histories (MAHs) and density profiles of cluster- size halos with virial masses of 0.6-2.5x10^14/h Msun in a flat LCDM cosmology. In agreement with previous studies,we find that the concentration of the density distribution is tightly correlated with the halo's MAH and its formation redshift.During the period of fast mass growth the concentration remains approximate…
▽ More
We analyze the mass accretion histories (MAHs) and density profiles of cluster- size halos with virial masses of 0.6-2.5x10^14/h Msun in a flat LCDM cosmology. In agreement with previous studies,we find that the concentration of the density distribution is tightly correlated with the halo's MAH and its formation redshift.During the period of fast mass growth the concentration remains approximately constant and low c_v ~ 3-4,while during the slow accretion stages it increases with decreasing redshift as c_v ~ (1+z)^-1.We consider fits of three widely discussed analytic density profiles to the simulated clusters focusing on the most relaxed inner regions.We find that there is no unique best fit analytic profile for all the systems.If,however,a cluster is best fit by a particular analytic profile at z=0,the same is usually true at earlier epochs out to z ~ 1-2.The local logarithmic slope of the density profiles at 3% of the virial radius ranges from -1.2 to -2.0,a remarkable diversity for the relatively narrow mass range of our cluster sample.For all the studied clusters the logarithmic slope becomes shallower with decreasing radius without reaching an asymptotic value down to the smallest resolved scale < ~1% of the virial radius.We do not find a clear correlation of the inner slope with the formation redshift or the shape of the halo's MAH.We do find, however,that during the period of rapid mass growth the density profiles can be well described by a single power law rho(r) ~ r^-gamma with gamma ~ 1.5-2.The relatively shallow power law slopes result in low concentrations at these stages of evolution.This indicates that the inner power law like density distribution of halos is built up during the periods of rapid mass accretion and active merging,while outer steeper profile is formed when the mass accretion is slow.
△ Less
Submitted 23 May, 2004; v1 submitted 4 November, 2003;
originally announced November 2003.
-
The Dark Side of the Halo Occupation Distribution
Authors:
Andrey V. Kravtsov,
Andreas A. Berlind,
Risa H. Wechsler,
Anatoly A. Klypin,
Stefan Gottloeber,
Brandon Allgood,
Joel R. Primack
Abstract:
We analyze the halo occupation distribution (HOD), the probability for a halo of mass M to host a number of subhalos N, and two-point correlation function of galaxy-size dark matter halos using high-resolution dissipationless simulations of the concordance flat LCDM model. The halo samples include both the host halos and the subhalos, distinct gravitationally-bound halos within the virialized re…
▽ More
We analyze the halo occupation distribution (HOD), the probability for a halo of mass M to host a number of subhalos N, and two-point correlation function of galaxy-size dark matter halos using high-resolution dissipationless simulations of the concordance flat LCDM model. The halo samples include both the host halos and the subhalos, distinct gravitationally-bound halos within the virialized regions of larger host systems. We find that the first moment of the HOD, <N>(M), has a complicated shape consisting of a step, a shoulder, and a power law high-mass tail. The HOD can be described by a Poisson statistics at high halo masses but becomes sub-Poisson for <N><4. We show that the HOD can be understood as a combination of the probability for a halo of mass M to host a central galaxy and the probability to host a given number Ns of satellite galaxies. The former can be approximated by a step-like function, while the latter can be well approximated by a Poisson distribution, fully specified by its first moment <Ns>(M). We find that <Ns>~M^b with b~1 for a wide range of number densities, redshifts, and different power spectrum normalizations. This formulation provides a simple but accurate model for the halo occupation distribution found in simulations. At z=0, the two-point correlation function (CF) of galactic halos can be well fit by a power law down to ~100/h kpc with an amplitude and slope similar to those of observed galaxies. At redshifts z>~1, we find significant departures from the power-law shape of the CF at small scales. If the deviations are as strong as indicated by our results, the assumption of the single power law often used in observational analyses of high-redshift clustering is likely to bias the estimates of the correlation length and slope of the correlation function.
△ Less
Submitted 14 May, 2004; v1 submitted 29 August, 2003;
originally announced August 2003.
-
Constrained Simulations of the Real Universe: II. Observational Signatures of Intergalactic Gas in the Local Supercluster Region
Authors:
Andrey V. Kravtsov,
Anatoly A. Klypin,
Yehuda Hoffman
Abstract:
We present results of gasdynamics+N-body constrained cosmological simulations of the Local Supercluster region (LSC; about 30/h Mpc around the Virgo cluster), which closely mimic the real Universe within 100 Mpc by imposing constraints from Mark III catalog of galaxy peculiar velocities. The simulations are used to study the properties and possible observational signatures of intergalactic mediu…
▽ More
We present results of gasdynamics+N-body constrained cosmological simulations of the Local Supercluster region (LSC; about 30/h Mpc around the Virgo cluster), which closely mimic the real Universe within 100 Mpc by imposing constraints from Mark III catalog of galaxy peculiar velocities. The simulations are used to study the properties and possible observational signatures of intergalactic medium in the LSC region. We find that, in agreement with previous unconstrained simulations, about 30% of the gas in this region is in the warm/hot phase at T~10^5-10^7 K, and about 40% in the diffuse phase at T<10^5 K in low-density regions. The best prospects for detection of the warm/hot intergalactic medium of the LSC located in filaments and in the vicinity of virialized regions of groups and clusters are through absorption in resonant lines of OVII and OVIII in soft X-rays and in the OVI doublet in UV. If intergalactic gas in filaments (rho/<rho>~ 1-10) is enriched to typical metallicities of >0.05, the column densities of OVI, OVII, and OVIII along a random line of sight near the North Galactic Pole, especially near the supergalactic plane, have a significant probability to be in the range detectable by the current (FUSE,XMM) and future (Constellation-X) instruments. (Abridged)
△ Less
Submitted 5 September, 2001;
originally announced September 2001.
-
The origin of angular momentum in dark matter halos
Authors:
M. Vitvitska,
A. A. Klypin,
A. V. Kravtsov,
R. H. Wechsler,
J. R. Primack,
J. S. Bullock
Abstract:
We propose a new explanation for the origin of angular momentum in galaxies and their dark halos, in which the halos obtain their spin through the cumulative acquisition of angular momentum from satellite accretion. In our model, the build-up of angular momentum is a random walk process associated with the mass assembly history of the halo's major progenitor. We assume no correlation between the…
▽ More
We propose a new explanation for the origin of angular momentum in galaxies and their dark halos, in which the halos obtain their spin through the cumulative acquisition of angular momentum from satellite accretion. In our model, the build-up of angular momentum is a random walk process associated with the mass assembly history of the halo's major progenitor. We assume no correlation between the angular momenta of accreted objects. Using the extended Press-Schechter approximation, we calculate the growth of mass, angular momentum, and spin parameter $λ$ for many halos. Our random walk model reproduces the key features of the angular momentum of halos found in N-body simulations: a lognormal distribution in $λ$ with an average of $<λ> \approx 0.04$, independent of mass and redshift. The evolution of the spin parameter in individual halos in this model is quite different from the steady increase with time of angular momentum in the tidal torque picture. We find both in N-body simulations and in our random walk model that the value of $λ$ changes significantly with time for a halo's major progenitor. It typically has a sharp increase due to major mergers, and a steady decline during periods of gradual accretion of small satellites. The model predicts that on average the $λ$ of halos which had major mergers after redshift $z=2$ should be substantially larger than the $λ$ of those which did not. Perhaps surprisingly, this suggests that halos that host late-forming elliptical galaxies should rotate faster than halos of spiral galaxies.
△ Less
Submitted 30 August, 2002; v1 submitted 20 May, 2001;
originally announced May 2001.
-
Spin Profile of Galactic Halos and Disk Formation
Authors:
A. Dekel,
J. S. Bullock,
C. Porciani,
A. V. Kravtsov,
T. S. Kolatt,
A. A. Klypin,
J. R. Primack
Abstract:
We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter λ. The angular momentum…
▽ More
We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter λ. The angular momentum tends to be aligned throughout the halo and of cylindrical symmetry. Even if angular momentum is conserved during baryonic infall, the resultant disk density profile is predicted to deviate from exponential, with a denser core and an extended tail. A slightly corrected version of the scaling relation due to linear tidal-torque theory is used to explain the origin of a typical power-law profile in shells, j(M) \propto M^s with s \gsim 1. While linear theory crudely predicts the amplitudes of halo spins, it is not a good predictor of their directions. Independently, mergers of halos are found to produce a similar profile due to j transfer from the orbit to the product halo via dynamical friction and tidal stripping. The halo spin is correlated with having a recent major merger, though this correlation is weakened by mass loss. These two effects are due to a correlation between the spins of neighboring halos and their orbit, leading to prograde mergers.
△ Less
Submitted 31 October, 2000;
originally announced November 2000.
-
A universal angular momentum profile for galactic halos
Authors:
James S. Bullock,
Avishai Dekel,
Tsafrir S. Kolatt,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Cristiano Porciani,
Joel R. Primack
Abstract:
[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global sp…
▽ More
[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global spin parameter λ. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of μ. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.
△ Less
Submitted 2 November, 2000; v1 submitted 31 October, 2000;
originally announced November 2000.
-
Origin and Destiny of Dark Matter Halos: Cosmological Matter Exchange and Metal Enrichment
Authors:
Tsafrir S. Kolatt,
James S. Bullock,
Avishai Dekel,
Joel R. Primack,
Yair Sigad,
Andrey V. Kravtsov,
Anatoly A. Klypin
Abstract:
[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virial…
▽ More
[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virialized halos. We follow the time evolution of the mass distribution and flows between these components and provide fitting functions for the exchange rates. We use our derived exchange rates to gauge the importance of metal redistribution in the universe due solely to gravity-induced interactions. The diffuse metallicity in clusters is predicted to be ~40% that in isolated galaxies (~55% of groups) at z=0, and should be lower only slightly by z=1, consistent with observations. The metallicity of the diffuse media in poor groups is expected to be lower by a factor of 5 by z~2, in agreement with the observed metallicity of damped Ly$α$ systems. The metallicity of the background IGM is predicted to be (1-3)x10^{-4} that of z=0 clusters, also consistent with observations. The agreement of predicted and observed trends indicates that gravitational interaction alone may play an important role in metal enrichment of the intra-cluster and intergalactic media.
△ Less
Submitted 11 October, 2000;
originally announced October 2000.
-
Interaction rates of dark-matter halos and subhalos
Authors:
Tsafrir S. Kolatt,
James S. Bullock,
Yair Sigad,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Avishai Dekel
Abstract:
[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo coll…
▽ More
[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo collisions dominate the total collision rate, and even at z=3 they are involved in more than 30% of all collisions. About 40% of the "major" mergers are between subhalos. Therefore subhalo interactions must be incorporated in models of structure formation. We find that the collision rate between halos in physical density units, is $\propto (1+z)^δ$, with $δ= 3-4$, in agreement with earlier simulations and most observational data. We test previous analytic estimates of the interaction rates of subhalos based on statistical models, which could be very inaccurate because of the small number of subhalos and the variation of conditions within small host halos. We find that, while such statistical estimates may severely overestimate the rate within hosts $< 10^{13} \hsolmass$, typical of high redshifts, they are valid for larger hosts regardless of the number of subhalos in them. We find the Makino & Hut (1997) estimate of the subhalo merger rate to be valid for hosts $ > 6x10^{11}\hsolmass$ at all redshifts. The collision rate between subhalos and the central object of their host halo is approximated relatively well using the timescale for dynamical friction in circular orbits. This approximation fails in ~40% of the cases, partly because of deviations from the assumption of circular orbits and partly because of the invalidity of the assumption that the subhalo mass is negligible.
△ Less
Submitted 11 October, 2000;
originally announced October 2000.
-
Velocity and Mass Functions of Galactic Halos: Evolution and Environmental Dependence
Authors:
Yair Sigad,
Tsafrir S. Kolatt,
James S. Bullock,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Avishai Dekel
Abstract:
We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $Λ$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular vel…
▽ More
We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $Λ$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular velocity function. The distribution functions in the velocity range 100--500 km/s are well fit by a power-law with two parameters, slope and amplitude. We present the parameter dependence on redshift and provide useful fitting formulae. The amplitudes of the mass functions decrease with z, but, contrary to naive expectation, the comoving density of halos of a fixed velocity ~200 km/s actually increases out to z=5. This is because high-z halos are denser, so a fixed velocity corresponds to a smaller mass. The slope of the velocity function at z=0 is as steep as ~ -4, and the mass and velocity functions of distinct halos steepen with increasing z, while the functions of subhalos do not steepen with $z$, and become even flatter at z>2. A simple observable prediction is that the slope of the velocity function of isolated galaxies is steeper than that of galaxies in groups by as much as unity, reflecting the density biasing of high-velocity halos. We confirm that the Press-Schechter approximation typically overestimates the halo mass function by a factor of ~2, while modified approximations provide improved predictions.
△ Less
Submitted 16 May, 2000;
originally announced May 2000.
-
Profiles of dark haloes: evolution, scatter, and environment
Authors:
James S. Bullock,
Tsafrir S. Kolatt,
Yair Sigad,
Rachel S. Somerville,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Avishai Dekel
Abstract:
We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the r…
▽ More
We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.
△ Less
Submitted 10 October, 2000; v1 submitted 15 August, 1999;
originally announced August 1999.
-
Quantifying the evolution of higher-order clustering
Authors:
Jens Schmalzing,
Stefan Gottloeber,
Anatoly A. Klypin,
Andrey V. Kravtsov
Abstract:
We use a high-resolution dissipationless simulation to study the evolution of the dark matter and halo distributions in a spatially flat cosmological model dominated by a cosmological constant $Λ$ and cold dark matter ($Λ$CDM). In order to quantify the evolution of structure, we calculate the Minkowski functionals of the halos and the dark matter component at various redshifts. A comparison of M…
▽ More
We use a high-resolution dissipationless simulation to study the evolution of the dark matter and halo distributions in a spatially flat cosmological model dominated by a cosmological constant $Λ$ and cold dark matter ($Λ$CDM). In order to quantify the evolution of structure, we calculate the Minkowski functionals of the halos and the dark matter component at various redshifts. A comparison of Minkowski functionals and the more standard correlation function analysis shows that the Minkowski functionals contain information about correlation functions of arbitrary order. While little evolution of the Minkowski functionals of halos between $z=4$ and $z=0$ is observed, we find that the Minkowski functionals of the dark matter evolve rapidly with time. The difference of the Minkowski functionals of halos and dark matter can be interpreted as a scale dependent bias. This implies that scale-dependent bias is a property of not only the two-point halo correlation function, but also of correlation functions of higher order.
△ Less
Submitted 29 June, 1999;
originally announced June 1999.
-
Young galaxies: what turns them on?
Authors:
Tsafrir S. Kolatt,
James S. Bullock,
Rachel S. Somerville,
Yair Sigad,
Patrik Jonsson,
Andrey V. Kravtsov,
Anatoly A. Klypin,
Joel R. Primack,
Sandra M. Faber,
Avishai Dekel,
.
Abstract:
Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to…
▽ More
Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $Λ$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos.
△ Less
Submitted 7 June, 1999;
originally announced June 1999.
-
Where are the missing galactic satellites?
Authors:
Anatoly A. Klypin,
Andrey V. Kravtsov,
Octavio Valenzuela,
Francisco Prada
Abstract:
Using published data, we have compiled the circular velocity (Vc) distribution function (VDF) of galaxy satellites in the Local Group. We find that within the volumes of radius of 570 kpc (400/h kpc for h=0.7) centered on the Milky Way and Andromeda, the average VDF is roughly approximated as n(>Vc)~ 45(Vc/10 km/s)^{-1} h^3 Mpc^{-3} for Vc in the range ~10-70 km/s. The observed VDF is compared w…
▽ More
Using published data, we have compiled the circular velocity (Vc) distribution function (VDF) of galaxy satellites in the Local Group. We find that within the volumes of radius of 570 kpc (400/h kpc for h=0.7) centered on the Milky Way and Andromeda, the average VDF is roughly approximated as n(>Vc)~ 45(Vc/10 km/s)^{-1} h^3 Mpc^{-3} for Vc in the range ~10-70 km/s. The observed VDF is compared with results of high-resolution cosmological simulations. We find that the VDF in models is very different from the observed one: n(>Vc)~1200(Vc/10 km/s)^{-2.75}h^3 Mpc^{-3}. Cosmological models thus predict that a halo of the size of our Galaxy should have about 50 dark matter satellites with circular velocity >20 km/s and mass >3x10^8/h Msun within a 570 kpc radius. This number is significantly higher than the approximate dozen satellites actually observed around our Galaxy. The observed and predicted VDFs cross at ~50 km/s, indicating that the predicted abundance of satellites with Vc> 50 km/s is in reasonably good agreement with observations. We conclude, therefore, that unless a large fraction of the Local Group satellites has been missed in observations, there is a dramatic discrepancy between observations and hierarchical models, regardless of the model parameters. We discuss several possible explanations for this discrepancy including identification of some satellites with the High Velocity Clouds observed in the Local Group, and the existence of dark satellites that failed to accrete gas and form stars due either to the expulsion of gas in the supernovae-driven winds or to gas heating by the intergalactic ionizing background. (Abridged)
△ Less
Submitted 20 January, 1999; v1 submitted 18 January, 1999;
originally announced January 1999.
-
Formation of Dark Matter Halos
Authors:
J. R. Primack,
J. S. Bullock,
A. A. Klypin,
A. V. Kravtsov
Abstract:
This article concerns the formation and structure of dark matter halos, including (1) their radial density profiles, (2) their abundance, and (3) their merger rates. The last topic may be relevant to the nature of the small, bright, high-redshift galaxies discovered by the Lyman break technique. (1) Study of a statistical sample of galaxy-mass dark halos in high-resolution Adaptive Refinement Tr…
▽ More
This article concerns the formation and structure of dark matter halos, including (1) their radial density profiles, (2) their abundance, and (3) their merger rates. The last topic may be relevant to the nature of the small, bright, high-redshift galaxies discovered by the Lyman break technique. (1) Study of a statistical sample of galaxy-mass dark halos in high-resolution Adaptive Refinement Tree simulations shows that they have a central density profile p(r) \propto 1/r^g with g \approx 0.2, in agreement with data on dark-matter-dominated disk galaxies. We present recent, higher resolution results on this. (2) Another important new result is that the Press-Schechter approximation predicts about twice as many galaxy-mass halos at z=0 as are present in large dissipationless N-body simulations; more generally, PS overpredicts the abundance of M ~< 0.1M_* halos at all redshifts. (3) Finally, we discuss the assembly of these halos, in particular the merger rate of (sub-)halos at high redshift and the distribution of the starbursts that these mergers are likely to trigger. If most of the Lyman-break galaxies are such starbursts, this perhaps resolves the apparent paradox that these galaxies appear to cluster like massive halos (~10^{12} M_\odot), while their relatively low linewidths and their spectral energy distributions suggest that they have relatively low mass (few X 10^{10} M_\odot) and young ages (few X 10^8 yr). It also predicts much more star formation at high redshift in CDM-type hierarchical models for structure formation than if only quiescent star formation is included.
△ Less
Submitted 11 December, 1998;
originally announced December 1998.
-
Halo evolution in a cosmological environment
Authors:
S. Gottloeber,
A. A. Klypin,
A. V. Kravtsov
Abstract:
We study the formation and evolution of the dark matter halos in a LCDM model. The dynamics of 16.8 million DM particles is followed numerically in a box of 60 Mpc/h with the dynamic range of 32,000 in spatial resolution. The high resolution of the simulation allows us to examine evolution of both isolated and satellite halos in a representative volume.
We present results on the evolution of v…
▽ More
We study the formation and evolution of the dark matter halos in a LCDM model. The dynamics of 16.8 million DM particles is followed numerically in a box of 60 Mpc/h with the dynamic range of 32,000 in spatial resolution. The high resolution of the simulation allows us to examine evolution of both isolated and satellite halos in a representative volume.
We present results on the evolution of velocity function of DM halos, compare it with the Press-Schechter function and discuss the evolution of power spectrum of matter and halo distributions. The mass evolution of isolated virialized objects determined from the simulation is in good agreement with predictions of the extended Press-Schechter models. However, satellite halos evolve very different: for some of them the mass decreases with time, which happens if the halos fall into clusters or groups. We discuss the dependence of the correlation function of halo populations on their environment and merging history.
△ Less
Submitted 27 October, 1998;
originally announced October 1998.
-
The Cores of Dark Matter Dominated Galaxies: theory vs. observations
Authors:
Andrey V. Kravtsov,
Anatoly A. Klypin,
James S. Bullock,
Joel R. Primack
Abstract:
We use the rotation curves of a sample of dark matter dominated dwarf and low-surface brightness (LSB) late-type galaxies to study their radial mass distributions. We find that the shape of the rotation curves is remarkably similar for all (both dwarf and LSB) galaxies in the sample, suggesting a self-similar density distribution of their dark matter (DM) halos. This shape can be reproduced well…
▽ More
We use the rotation curves of a sample of dark matter dominated dwarf and low-surface brightness (LSB) late-type galaxies to study their radial mass distributions. We find that the shape of the rotation curves is remarkably similar for all (both dwarf and LSB) galaxies in the sample, suggesting a self-similar density distribution of their dark matter (DM) halos. This shape can be reproduced well by a density profile with a shallow central cusp (rho(r)\propto 1/r^gamma, gamma\approx 0.2-0.4), corresponding to a steeply rising velocity curve (v(r)\propto r^g, g\approx 0.9-0.8). We further show that the observed shape of the rotation curves is well matched by the density profiles of DM halos formed in very high resolution simulations of the CDM, LCDM, and CHDM models of structure formation. This is surprising in light of several previous studies which suggested that the structure of simulated DM halos is inconsistent with the dynamics of dwarf galaxies. We discuss possible explanations for this discrepancy. We show that although the mass distribution in the hierarchically formed halos is on average consistent with the shape of rotation curves of dark matter dominated galaxies, the scatter of the individual profiles around the average is substantial and should not be neglected in comparisons with the data. Finally, we show that the DM halos in our hierarchical simulations and the real galaxies in our sample exhibit a similar decrease in their characteristic densities with increasing characteristic radial scales, and increase in their maximum rotation velocities with increasing radii at which their maximum velocity occurs. (Abridged)
△ Less
Submitted 25 May, 1998; v1 submitted 19 August, 1997;
originally announced August 1997.
-
Adaptive Refinement Tree - a new high-resolution N-body code for cosmological simulations
Authors:
Andrey V. Kravtsov,
Anatoly A. Klypin,
Alexei M. Khokhlov
Abstract:
We present a new high-resolution N-body algorithm for cosmological simulations. The algorithm employs a traditional particle-mesh technique on a cubic grid and successive multilevel relaxations on the finer meshes, introduced recursively in a fully adaptive manner in the regions where the density exceeds a predefined threshold. The mesh is generated to effectively match an arbitrary geometry of…
▽ More
We present a new high-resolution N-body algorithm for cosmological simulations. The algorithm employs a traditional particle-mesh technique on a cubic grid and successive multilevel relaxations on the finer meshes, introduced recursively in a fully adaptive manner in the regions where the density exceeds a predefined threshold. The mesh is generated to effectively match an arbitrary geometry of the underlying density field -- a property particularly important for cosmological simulations. In a simulation the mesh structure is not created at every time step but is properly adjusted to the evolving particle distribution. The algorithm is fast and effectively parallel. We present a detailed description of the methodology, implementation, and tests of the code.
We further use the code to study the structure of dark matter halos in high-resolution (2/h kpc) simulations of standard CDM (Omega=1, h=0.5, sigma_8=0.63) and LCDM (Omega_Lambda=1-Omega_0=0.7, h=0.7, sigma_8=1.0) models. We find that halo density profiles in both CDM and LCDM models are well fitted by the analytical model presented recently by Navarro et al., which predicts a singular [$ρ(r)\propto r^{-1}$] behavior of the halo density profiles at small radii. We therefore conclude that halos formed in the $Λ$CDM model have structure similar to CDM halos and thus cannot explain the dynamics of the central parts of dwarf spiral galaxies, as inferred from the galaxies' rotation curves.
△ Less
Submitted 28 July, 1997; v1 submitted 24 January, 1997;
originally announced January 1997.
-
TWO-POINT ANGULAR CORRELATION FUNCTION FOR THE GREEN BANK 4.85 GHZ SKY SURVEY
Authors:
Brian L. Kooiman,
Jack O. Burns,
Anatoly A. Klypin
Abstract:
This paper presents an angular correlation analysis of the Green Bank 4.85 GHz radio catalog (Gregory \& Condon 1991) of 54,579 sources (S \gsim 25 mJy). The Green Bank catalog is found to be complete to S $\geq$ 35 mJy over 20$^{\circ} \leq δ< 74^{\circ}$, 0$^h \leq α< 24^h$, and Galactic latitude $|b| \geq 10^{\circ}$. The 2-point angular correlation function shows evidence for the clustering of…
▽ More
This paper presents an angular correlation analysis of the Green Bank 4.85 GHz radio catalog (Gregory \& Condon 1991) of 54,579 sources (S \gsim 25 mJy). The Green Bank catalog is found to be complete to S $\geq$ 35 mJy over 20$^{\circ} \leq δ< 74^{\circ}$, 0$^h \leq α< 24^h$, and Galactic latitude $|b| \geq 10^{\circ}$. The 2-point angular correlation function shows evidence for the clustering of radio sources, with a power-law distribution consistent with a slope $γ= -0.8$. This may well provide the $first$ detection of an angular correlation in a large area, complete deep radio survey.
△ Less
Submitted 26 February, 1995;
originally announced February 1995.
-
Galaxy Group Analysis - a Robust Discriminator between Cosmological Models: Cold+HOT Dark Matter and CDM Confront Cfa1
Authors:
R. Nolthenius,
A. A. Klypin,
J. R. Primack
Abstract:
We present techniques for identifying and analyzing galaxy groups and show that they provide a powerful discriminator between cosmological models. We apply these to high-resolution PM N-body simulations of structure formation in CHDM (cold/hot/baryon fractions=.6/.3/.1, b=1.5 (COBE norm) and two CDM models; b=1.5 and b=1.0 (COBE norm). Groups are identified with the Nolthenius (1993) algorithm.…
▽ More
We present techniques for identifying and analyzing galaxy groups and show that they provide a powerful discriminator between cosmological models. We apply these to high-resolution PM N-body simulations of structure formation in CHDM (cold/hot/baryon fractions=.6/.3/.1, b=1.5 (COBE norm) and two CDM models; b=1.5 and b=1.0 (COBE norm). Groups are identified with the Nolthenius (1993) algorithm. We find: (1) Properties of groups are a powerful and robust discriminator between these Gaussian models whose spatial and velocity properties differ on Mpc scales. We test robustness against several methods for assigning luminosity to dark matter halos, for merging CfA1 galaxies, and for breaking up massive over-mergers (2) When allowance is made for the high large scale power P(k) present in the CfA1 data, CHDM at our $\Omnu=0.30$ produces slightly too many groups and groups too high a fraction of the galaxies, while the fraction grouped in CDM is far too low. Lowering $\Omnu$ would appear to produce close agreement with all measures, save one: for all simulations, median group sizes are x1.4-1.7 larger than equivalently selected CfA1 groups. (3) The standard group M/L method gives $Ω\simeq 0.1$ for CHDM and CfA1, and $Ω\simeq 0.35$ for CDM. A conspiracy of three effects produce these very low values, and we conclude that low observed $Ω$'s need not argue for a low $Ω$ universe. When overmergers are broken up, the median virial- to-DM mass of 3D selected groups is $\sim1$ for all simulations. $M_{DM}>10^{14}\Msun$ groups appear virialized in all simulations. Global velocity biases $b_v$ are similar to previous studies. Within 3D-selected groups, CHDM and CDM b=1.5 show a stronger bias of $b_v=0.7-0.8$, while CDM b=1.0 shows group $b_v$'s $\simeq 1$.
△ Less
Submitted 29 October, 1994;
originally announced October 1994.