-
Supermassive Black Holes and Their Environments
Authors:
Joerg M. Colberg,
Tiziana Di Matteo
Abstract:
We make use of the first high--resolution hydrodynamic simulations of structure formation which self-consistently follows the build up of supermassive black holes introduced in Di Matteo et al. (2007) to investigate the relation between black holes (BH), host halo and large--scale environment. There are well--defined relations between halo and black hole masses and between the activities of gala…
▽ More
We make use of the first high--resolution hydrodynamic simulations of structure formation which self-consistently follows the build up of supermassive black holes introduced in Di Matteo et al. (2007) to investigate the relation between black holes (BH), host halo and large--scale environment. There are well--defined relations between halo and black hole masses and between the activities of galactic nuclei and halo masses at low redshifts. A large fraction of black holes forms anti--hierarchically, with a higher ratio of black hole to halo mass at high than at low redshifts. At $z=1$, we predict group environments (regions of enhanced local density) to contain the highest mass and most active (albeit with a large scatter) BHs while the rest of the BH population to be spread over all densities from groups to filaments and voids. Density dependencies are more pronounced at high rather than low redshift. These results are consistent with the idea that gas rich mergers are likely the main regulator of quasar activity. We find star formation to be a somewhat stronger and tighter function of local density than BH activity, indicating some difference in the triggering of the latter versus the former. There exists a large number of low--mass black holes, growing slowly predominantly through accretion, which extends all the way into the most underdense regions, i.e. in voids.
△ Less
Submitted 10 April, 2008;
originally announced April 2008.
-
The Aspen--Amsterdam Void Finder Comparison Project
Authors:
Joerg M. Colberg,
Frazer Pearce,
Caroline Foster,
Erwin Platen,
Riccardo Brunino,
Mark Neyrinck,
Spyros Basilakos,
Anthony Fairall,
Hume Feldman,
Stefan Gottloeber,
Oliver Hahn,
Fiona Hoyle,
Volker Mueller,
Lorne Nelson,
Manolis Plionis,
Cristiano Porciani,
Sergei Shandarin,
Michael S. Vogeley,
Rien van de Weygaert
Abstract:
Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recen…
▽ More
Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.
△ Less
Submitted 10 April, 2008; v1 submitted 6 March, 2008;
originally announced March 2008.
-
On the importance of high redshift intergalactic voids
Authors:
Matteo Viel,
Joerg M. Colberg,
Tae-Sun Kim
Abstract:
We investigate the properties of one--dimensional flux ``voids'' (connected regions in the flux distribution above the mean flux level) by comparing hydrodynamical simulations of large cosmological volumes with a set of observed high--resolution spectra at z ~ 2. After addressing the effects of box size and resolution, we study how the void distribution changes when the most significant cosmolog…
▽ More
We investigate the properties of one--dimensional flux ``voids'' (connected regions in the flux distribution above the mean flux level) by comparing hydrodynamical simulations of large cosmological volumes with a set of observed high--resolution spectra at z ~ 2. After addressing the effects of box size and resolution, we study how the void distribution changes when the most significant cosmological and astrophysical parameters are varied. We find that the void distribution in the flux is in excellent agreement with predictions of the standard LCDM cosmology, which also fits other flux statistics remarkably well. We then model the relation between flux voids and the corresponding one--dimensional gas density field along the line--of--sight and make a preliminary attempt to connect the one--dimensional properties of the gas density field to the three--dimensional dark matter distribution at the same redshift. This provides a framework that allows statistical interpretations of the void population at high redshift using observed quasar spectra, and eventually it will enable linking the void properties of the high--redshift universe with those at lower redshifts, which are better known.
△ Less
Submitted 20 February, 2008; v1 submitted 8 January, 2008;
originally announced January 2008.
-
Quantifying Cosmic Superstructures
Authors:
J. M. Colberg
Abstract:
The Large Scale Structure (LSS) found in galaxy redshift surveys and in computer simulations of cosmic structure formation shows a very complex network of galaxy clusters, filaments, and sheets around large voids. Here, we introduce a new algorithm, based on a Minimal Spanning Tree, to find basic structural elements of this network and their properties. We demonstrate how the algorithm works usi…
▽ More
The Large Scale Structure (LSS) found in galaxy redshift surveys and in computer simulations of cosmic structure formation shows a very complex network of galaxy clusters, filaments, and sheets around large voids. Here, we introduce a new algorithm, based on a Minimal Spanning Tree, to find basic structural elements of this network and their properties. We demonstrate how the algorithm works using simple test cases and then apply it to haloes from the Millennium Run simulation (Springel et al. 2005). We show that about 70% of the total halo mass is contained in a structure composed of more than 74,000 individual elements, the vast majority of which are filamentary, with lengths of up to 15 Mpc/h preferred. Spatially more extended structures do exist, as do examples of what appear to be sheet-like configurations of matter. What is more, LSS appears to be composed of a fixed set of basic building blocks. The LSS formed by mass selected subsamples of haloes shows a clear correlation between the threshold mass and the mean extent of major branches, with cluster-size haloes forming structures whose branches can extend to almost 200 Mpc/h - the backbone of LSS to which smaller branches consisting of smaller haloes are attached.
△ Less
Submitted 20 November, 2006;
originally announced November 2006.
-
The Influence of Large-Scale Structure on Halo Shapes and Alignments
Authors:
G. Altay,
J. M. Colberg,
R. A. C. Croft
Abstract:
Alignments of galaxy clusters (the Binggeli effect), as well as of galaxies themselves have long been studied both observationally and theoretically. Here we test the influence of large-scales structures and tidal fields on the shapes and alignments of cluster-size and galaxy-size dark matter halos. We use a high-resolution N-body simulation of a $Λ$CDM universe, together with the results of Col…
▽ More
Alignments of galaxy clusters (the Binggeli effect), as well as of galaxies themselves have long been studied both observationally and theoretically. Here we test the influence of large-scales structures and tidal fields on the shapes and alignments of cluster-size and galaxy-size dark matter halos. We use a high-resolution N-body simulation of a $Λ$CDM universe, together with the results of Colberg et al. (2005), who identified filaments connecting pairs of clusters. We find that cluster pairs connected by a filament are strongly aligned with the cluster-cluster axis, whereas unconnected ones are not. For smaller, galaxy-size halos, there also is an alignment signal, but its strength is independent of whether the halo is part of an obvious large-scale structure. Additionally, we find no measureable dependence of galaxy halo shape on membership of a filament. We also quantify the influence of tidal fields and find that these do correlate strongly with alignments of halos. The alignments of most halos are thus caused by tidal fields, with cluster-size halos being strongly aligned through the added mechanism of infall of matter from filaments.
△ Less
Submitted 11 May, 2006;
originally announced May 2006.
-
Voids in a $Λ$CDM Universe
Authors:
Joerg M. Colberg,
Ravi K. Sheth,
Antonaldo Diaferio,
Liang Gao,
Naoki Yoshida
Abstract:
We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular $Λ$ Cold Dark Matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of -0.8 or less. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement. Voids very c…
▽ More
We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular $Λ$ Cold Dark Matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of -0.8 or less. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement. Voids very clearly correspond to minima in the smoothed initial density field. We find a universal void mass profile of the form $ρ(<r)/ρ(r_{\rm eff}) \propto \exp[(r/r_{\rm eff})^α]$ where $r_{\rm eff}$ is the effective radius of a void and $α\sim 2$. The mass function of haloes in voids is steeper than that of haloes that populate denser regions. The abundances of void haloes seem to evolve somewhat more strongly between redshifts $\sim 1$ and 0 than the global abundances of haloes. (abridged)
△ Less
Submitted 5 April, 2005; v1 submitted 7 September, 2004;
originally announced September 2004.
-
Inter-cluster filaments in a $Λ$CDM Universe
Authors:
Joerg M. Colberg,
K. Simon Krughoff,
Andrew J. Connolly
Abstract:
The large--scale structure (LSS) in the Universe comprises a complicated filamentary network of matter. We study this network using a high--resolution simulation of structure formation of a $Λ$ Cold Dark Matter cosmology. We investigate the distribution of matter between neighbouring large haloes whose masses are comparable to massive clusters of galaxies. We identify a total of 228 filaments be…
▽ More
The large--scale structure (LSS) in the Universe comprises a complicated filamentary network of matter. We study this network using a high--resolution simulation of structure formation of a $Λ$ Cold Dark Matter cosmology. We investigate the distribution of matter between neighbouring large haloes whose masses are comparable to massive clusters of galaxies. We identify a total of 228 filaments between neighbouring clusters. Roughly half of the filaments are either warped or lie off the cluster--cluster axis. We find that straight filaments on the average are shorter than warped ones. More massive clusters are connected to more filaments than less massive ones on average. This finding indicates that the most massive clusters form at the intersections of the filamentary backbone of LSS. For straight filaments, we compute mass profiles. Radial profiles show a fairly well--defined radius, $r_s$, beyond which the profiles follow an $r^{-2}$ power law fairly closely. For the majority of filaments, $r_s$ lies between 1.5 $h^{-1}$ Mpc and 2.0 $h^{-1}$ Mpc. The enclosed overdensity inside $r_s$ varies between a few times up to 25 times mean density, independent of the length of the filaments. Along the filaments' axes, material is not distributed uniformly. Towards the clusters, the density rises, indicating the presence of the cluster infall regions. In addition, we also find some sheet--like connections between clusters. In roughly a fifth of all cluster--cluster connections where we could not identify a filament or sheet, projection effects lead to filamentary structures in the projected mass distribution. (abridged)
△ Less
Submitted 4 February, 2005; v1 submitted 29 June, 2004;
originally announced June 2004.
-
Galaxy Clusters in Hubble Volume Simulations: Cosmological Constraints from Sky Survey Populations
Authors:
A. E. Evrard,
T. J. MacFarland,
H. M. P. Couchman,
J. M. Colberg,
N. Yoshida,
S. D. M. White,
A. R. Jenkins,
C. S. Frenk,
F. R. Pearce,
J. A. Peacock,
P. A. Thomas
Abstract:
We use giga-particle N-body simulations to study galaxy cluster populations in Hubble Volumes of LCDM (Omega_m=0.3, Omega_Lambda=0.7) and tCDM (Omega_m=1) world models. Mapping past light-cones of locations in the computational space, we create mock sky surveys of dark matter structure to z~1.4 over 10,000 sq deg and to z~0.5 over two full spheres. Calibrating the Jenkins mass function at z=0 wi…
▽ More
We use giga-particle N-body simulations to study galaxy cluster populations in Hubble Volumes of LCDM (Omega_m=0.3, Omega_Lambda=0.7) and tCDM (Omega_m=1) world models. Mapping past light-cones of locations in the computational space, we create mock sky surveys of dark matter structure to z~1.4 over 10,000 sq deg and to z~0.5 over two full spheres. Calibrating the Jenkins mass function at z=0 with samples of ~1.5 million clusters, we show that the fit describes the sky survey counts to <~20% acccuracy over all redshifts for systems larger than poor groups (M>5e13 Msun/h). Fitting the observed local temperature function determines the ratio beta of specific thermal energies in dark matter and intracluster gas. We derive a scaling with power spectrum normalization beta \propto sigma8^{5/3}, and measure a 4% error on sigma8 arising from cosmic variance in temperature-limited cluster samples. Considering distant clusters, the LCDM model matches EMSS and RDCS X-ray-selected survey observations under economical assumptions for intracluster gas evolution. Using transformations of mass-limited cluster samples that mimic sigma8 variation, we explore SZ search expectations for a 10 sq deg survey complete above 10^{14} Msun/h. Cluster counts are shown to be extremely sensitive to sigma8 uncertainty while redshift statistics, such as the sample median, are much more stable. For LCDM, the characteristic temperature at fixed sky surface density is a weak function of redshift, implying an abundance of hot clusters at z>1. Assuming constant beta, four kT>8 keV clusters lie at z>2 and 40 kT>5 keV clusters lie at z>3 on the whole sky. Detection of Coma-sized clusters at z>1 violate LCDM at 95% confidence if their surface density exceeds 0.003 per sq deg, or 120 on the whole sky.
△ Less
Submitted 11 March, 2002; v1 submitted 10 October, 2001;
originally announced October 2001.
-
Public Release of N-body simulation and related data by the Virgo consortium
Authors:
C. S. Frenk,
J. M. Colberg,
H. M. P. Couchman,
G. Efstathiou,
A. E. Evrard,
A. Jenkins,
T. J. MacFarland,
B. Moore,
J. A. Peacock,
F. R. Pearce,
P. A. Thomas,
S. D. M. White,
N. Yoshida
Abstract:
We are making available on the WWW a selection of the archived data from N-body simulations carried out by the Virgo consortium and related groups. This currently includes: (i) time-slice, lightcone and cluster data from the two $10^9$-particle Hubble volume simulations described by Evrard 1998; (ii) time-slice data from simulations of 4 different cold dark matter cosmological models with…
▽ More
We are making available on the WWW a selection of the archived data from N-body simulations carried out by the Virgo consortium and related groups. This currently includes: (i) time-slice, lightcone and cluster data from the two $10^9$-particle Hubble volume simulations described by Evrard 1998; (ii) time-slice data from simulations of 4 different cold dark matter cosmological models with $256^3$ particles analysed by Jenkins et al 1998; (iii) Dark halo catalogs, merger trees and galaxy catalogs from the GIF project described by Kauffmann et al 1999. Basic software is supplied to read the data. The data can be accessed from: http://www.mpa-garching.mpg.de/Virgo/data_download.html
△ Less
Submitted 25 July, 2000;
originally announced July 2000.
-
Mass function of dark matter halos
Authors:
A. Jenkins,
C. S. Frenk,
S. D. M. White,
J. M. Colberg,
S. Cole,
A. E. Evrard,
H. M. P. Couchman,
N. Yoshida
Abstract:
We combine data from a number of N-body simulations to predict the abundance of dark halos in Cold Dark Matter universes over more than 4 orders of magnitude in mass. A comparison of different simulations suggests that the dominant uncertainty in our results is systematic and is smaller than 10--30% at all masses, depending on the halo definition used. In particular, our ``Hubble Volume'' simula…
▽ More
We combine data from a number of N-body simulations to predict the abundance of dark halos in Cold Dark Matter universes over more than 4 orders of magnitude in mass. A comparison of different simulations suggests that the dominant uncertainty in our results is systematic and is smaller than 10--30% at all masses, depending on the halo definition used. In particular, our ``Hubble Volume'' simulations of \tcdm and \lcdm cosmologies allow the abundance of massive clusters to be predicted with uncertainties well below those expected in all currently planned observational surveys. We show that for a range of CDM cosmologies and for a suitable halo definition, the simulated mass function is almost independent of epoch, of cosmological parameters, and of initial power spectrum when expressed in appropriate variables. This universality is of exactly the kind predicted by the familiar Press-Schechter model, although this model predicts a mass function shape which differs from our numerical results, overestimating the abundance of ``typical'' halos and underestimating that of massive systems.
△ Less
Submitted 25 September, 2000; v1 submitted 12 May, 2000;
originally announced May 2000.
-
Clustering of Galaxy Clusters in CDM Universes
Authors:
J. M. Colberg,
S. D. M. White,
N. Yoshida,
T. J. MacFarland,
A. Jenkins,
C. S. Frenk,
F. R. Pearce,
A. E. Evrard,
H. M. P. Couchman,
G. Efstathiou,
J. A. Peacock,
P. A. Thomas
Abstract:
We use very large cosmological N--body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ($τ$CDM) and a flat low density model with a cosmological constant ($Λ$CDM). Our simulations each use $10^9$ particles to follow…
▽ More
We use very large cosmological N--body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ($τ$CDM) and a flat low density model with a cosmological constant ($Λ$CDM). Our simulations each use $10^9$ particles to follow the mass distribution within cubes of side $2h^{-1}$Gpc ($τ$CDM) and $3h^{-1}$Gpc ($Λ$CDM) with a force resolution better than $10^{-4}$ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from $r_0=12$ -- 13$h^{-1}$Mpc for samples with mean separation $d_{\rm c}=30h^{-1}$Mpc to $r_0=22$-- 27$h^{-1}$Mpc for samples with $d_{\rm c}=100h^{-1}$Mpc. The lower value here corresponds to $τ$CDM and the upper to $Λ$CDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10% using the simple model of Sheth, Mo & Tormen (2000). This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the APM cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for $Λ$CDM. The observed power spectrum (for a cluster sample of mean separation $d_{\rm c}=31h^{-1}$Mpc) lies significantly above the predictions of both models.
△ Less
Submitted 11 July, 2000; v1 submitted 11 May, 2000;
originally announced May 2000.
-
Galaxy clustering determined from numerical cosmological simulations
Authors:
Adrian Jenkins,
C. S. Frenk,
F. R. Pearce,
P. A. Thomas,
J. M. Colberg,
S. D. M. White,
H. M. P. Couchman,
J. A. Peacock,
G. Efstathiou,
A. H. Nelson
Abstract:
We have simulated the growth of structure in two 100 Mpc boxes for LCDM and SCDM universes. These N-body/SPH simulations include a gaseous component which is able to cool radiatively. A fraction of the gas cools into cold dense objects which we identify as galaxies. In this article we give a preliminary analysis of the clustering behaviour of these galaxies concentrating on the LCDM model. We fi…
▽ More
We have simulated the growth of structure in two 100 Mpc boxes for LCDM and SCDM universes. These N-body/SPH simulations include a gaseous component which is able to cool radiatively. A fraction of the gas cools into cold dense objects which we identify as galaxies. In this article we give a preliminary analysis of the clustering behaviour of these galaxies concentrating on the LCDM model. We find a galaxy correlation function which is very close to a power law and which evolves relatively little with redshift. The pairwise dispersions of the galaxies are significantly lower than the dark matter. The LCDM model gives a surprisingly good match to the observational determinations of the galaxy correlation function and pairwise dispersions.
△ Less
Submitted 2 June, 1999;
originally announced June 1999.
-
Cosmological galaxy formation
Authors:
F. R. Pearce,
C. S. Frenk,
A. Jenkins,
J. M. Colberg,
P. A. Thomas,
H. M. P. Couchman,
S. D. M. White,
G. Efstathiou,
J. A. Peacock,
A. H. Nelson
Abstract:
We discuss the first results from two successful simulations of galaxy formation within a cosmological volume. With over 2000 large objects forming in each we have sufficient numbers to reliably produce both galaxy correlation and luminosity functions. We find that the observed galaxy counts are well fitted by these models and that the galaxies display an almost un-evolving correlation function…
▽ More
We discuss the first results from two successful simulations of galaxy formation within a cosmological volume. With over 2000 large objects forming in each we have sufficient numbers to reliably produce both galaxy correlation and luminosity functions. We find that the observed galaxy counts are well fitted by these models and that the galaxies display an almost un-evolving correlation function back to a redshift of 3 which closely resembles the featureless observed form and amplitude.
△ Less
Submitted 2 June, 1999;
originally announced June 1999.
-
A simulation of galaxy formation and clustering
Authors:
F. R. Pearce,
A. Jenkins,
C. S. Frenk,
J. M. Colberg,
S. D. M. White,
P. A. Thomas,
H. M. P. Couchman,
J. A. Peacock,
G. Efstathiou
Abstract:
We discuss early results from the first large N-body/hydrodynamical simulation to resolve the formation of galaxies in a cold dark matter universe. The simulation follows the formation of galaxies by gas cooling within dark halos of mass a few times $10^{11}\Msun$ and above, in a flat universe with a positive cosmological constant. Over 2200 galaxies form in our simulated volume of…
▽ More
We discuss early results from the first large N-body/hydrodynamical simulation to resolve the formation of galaxies in a cold dark matter universe. The simulation follows the formation of galaxies by gas cooling within dark halos of mass a few times $10^{11}\Msun$ and above, in a flat universe with a positive cosmological constant. Over 2200 galaxies form in our simulated volume of $(100 \Mpc)^3$. Assigning luminosities to the model galaxies using a spectral population synthesis model results in a K-band luminosity function in excellent agreement with observations. The two-point correlation function of galaxies in the simulation evolves very little since $z=3$ and has a shape close to a power-law over four orders of magnitude in amplitude. At the present day, the galaxy correlation function in the simulation is antibiased relative to the mass on small scales and unbiased on large scales. It provides a reasonable match to observations.
△ Less
Submitted 13 May, 1999;
originally announced May 1999.
-
Clustering of Galaxies in a Hierarchical Universe: III. Mock Redshift Surveys
Authors:
Antonaldo Diaferio,
Guinevere Kauffmann,
Joerg M. Colberg,
Simon D. M. White
Abstract:
This is the third paper in a series which combines N-body simulations and semi-analytic modelling to provide a fully spatially resolved simulation of the galaxy formation and clustering processes. Here we extract mock redshift surveys from our simulations: a Cold Dark Matter model with either Omega_0=1 (tauCDM) or Omega_0=0.3 and Lambda=0.7 (LambdaCDM). We compare the mock catalogues with the no…
▽ More
This is the third paper in a series which combines N-body simulations and semi-analytic modelling to provide a fully spatially resolved simulation of the galaxy formation and clustering processes. Here we extract mock redshift surveys from our simulations: a Cold Dark Matter model with either Omega_0=1 (tauCDM) or Omega_0=0.3 and Lambda=0.7 (LambdaCDM). We compare the mock catalogues with the northern region (CfA2N) of the Center for Astrophysics (CfA) Redshift Surveys. We study the properties of galaxy groups and clusters identified using standard observational techniques and we study the relation of these groups to real virialised systems. Most features of CfA2N groups are reproduced quite well by both models with no obvious dependence on Omega_0. Redshift space correlations and pairwise velocities are also similar in the two cosmologies. The luminosity functions predicted by our galaxy formation models depend sensitively on the treatment of star formation and feedback. For the particular choices of Paper I they agree poorly with the CfA survey. To isolate the effect of this discrepancy on our mock redshift surveys, we modify galaxy luminosities in our simulations to reproduce the CfA luminosity function exactly. This adjustment improves agreement with the observed abundance of groups, which depends primarily on the galaxy luminosity density, but other statistics, connected more closely with the underlying mass distribution, remain unaffected. Regardless of the luminosity function adopted, modest differences with observation remain. These can be attributed to the presence of the ``Great Wall'' in the CfA2N. It is unclear whether the greater coherence of the real structure is a result of cosmic variance, given the relatively small region studied, or reflects a physical deficiency of the models.
△ Less
Submitted 29 March, 1999; v1 submitted 1 December, 1998;
originally announced December 1998.
-
Clustering of Galaxies in a Hierarchical Universe: II. Evolution to High Redshift
Authors:
Guinevere Kauffmann,
Joerg M. Colberg,
Antonaldo Diaferio,
Simon D. M. White
Abstract:
In hierarchical cosmologies the evolution of galaxy clustering depends both on cosmological quantities such as Omega and Lambda, which determine how dark matter halos form and evolve, and on the physical processes - cooling, star formation and feedback - which drive the formation of galaxies within these merging halos. In this paper, we combine dissipationless cosmological N-body simulations and…
▽ More
In hierarchical cosmologies the evolution of galaxy clustering depends both on cosmological quantities such as Omega and Lambda, which determine how dark matter halos form and evolve, and on the physical processes - cooling, star formation and feedback - which drive the formation of galaxies within these merging halos. In this paper, we combine dissipationless cosmological N-body simulations and semi-analytic models of galaxy formation in order to study how these two aspects interact. We focus on the differences in clustering predicted for galaxies of differing luminosity, colour, morphology and star formation rate and on what these differences can teach us about the galaxy formation process. We show that a "dip" in the amplitude of galaxy correlations between z=0 and z=1 can be an important diagnostic. Such a dip occurs in low-density CDM models because structure forms early and dark matter halos of 10**12 solar masses, containing galaxies with luminosities around L*, are unbiased tracers of the dark matter over this redshift range; their clustering amplitude then evolves similarly to that of the dark matter. At higher redshifts bright galaxies become strongly biased and the clustering amplitude increases again. In high density models, structure forms late and bias evolves more rapidly. As a result, the clustering amplitude of L* galaxies remains constant from z=0 to 1. The strength of these effects is sensitive to sample selection. The dip is weaker for galaxies with lower star formation rates, redder colours, higher luminosities and earlier morphological types. We explain why this is the case and how it is related to the variation with redshift of the abundance and and environment of the observed galaxies. Studies of clustering evolution as a function of galaxy properties should place strong constraints on models.
△ Less
Submitted 18 September, 1998; v1 submitted 14 September, 1998;
originally announced September 1998.
-
Galaxy Clusters in the Hubble Volume Simulations
Authors:
Joerg M. Colberg,
Simon D. M. White,
Thomas J. MacFarland,
Adrian Jenkins,
Carlos S. Frenk,
F. R. Pearce,
August E. Evrard,
H. M. P. Couchman,
George Efstathiou,
John A. Peacock,
Peter A. Thomas
Abstract:
We report on analyses of cluster samples obtained from the Hubble Volume Simulations. These simulations, an $Ω=1$ model named $τ$CDM and a flat low $Ω$ model with a cosmological constant ($Λ$CDM), comprise the largest computational efforts to date in numerical cosmology. We investigate the presence of massive galaxy clusters at $z\approx 0.8$. The $τ$CDM model fails to form clusters at such a re…
▽ More
We report on analyses of cluster samples obtained from the Hubble Volume Simulations. These simulations, an $Ω=1$ model named $τ$CDM and a flat low $Ω$ model with a cosmological constant ($Λ$CDM), comprise the largest computational efforts to date in numerical cosmology. We investigate the presence of massive galaxy clusters at $z\approx 0.8$. The $τ$CDM model fails to form clusters at such a redshift. However, due to the small number of observed clusters around $z\approx 0.8$ and the uncertainties in the determinations of their masses, this conclusion still is somewhat preliminary. We produce cluster catalogs at $z=0$ for both cosmologies and investigate their two--point correlation function $ξ$. We show that the relationship between the mean density of subsamples of clusters, expressed via their mean separation $d_{\rm c}$, and the correlation length $r_0$, defined through $ξ(r_0) = 1$, is not linear but turns over gently for large $d_{\rm c}$. An analytic prediction by Mo & White (1996) overpredicts $r_0$. The results from the analysis of the APM cluster data by Croft et al. (1997) are nicely matched by the $Λ$CDM model.
△ Less
Submitted 24 August, 1998;
originally announced August 1998.
-
Clustering of Galaxies in a Hierarchical Universe: I. Methods and Results at z=0
Authors:
Guinevere Kauffmann,
Joerg M. Colberg,
Antonaldo Diaferio,
Simon D. M. White
Abstract:
We introduce a new technique for following the formation and evolution of galaxies in cosmological N-body simulations. Dissipationless simulations are used to track the formation and merging of dark matter halos as a function of redshift. Simple prescriptions, taken directly from semi-analytic models of galaxy formation, are adopted for cooling, star formation, supernova feedback and the merging…
▽ More
We introduce a new technique for following the formation and evolution of galaxies in cosmological N-body simulations. Dissipationless simulations are used to track the formation and merging of dark matter halos as a function of redshift. Simple prescriptions, taken directly from semi-analytic models of galaxy formation, are adopted for cooling, star formation, supernova feedback and the merging of galaxies within the halos. This scheme enables us to study the clustering properties of galaxies and to investigate how selection by type, colour or luminosity influences the results. In this paper, we study properties of the galaxy distribution at z=0. These include luminosity functions, colours, correlation functions, pairwise peculiar velocities, cluster M/L ratios and star formation rates. We focus on two variants of a CDM cosmology: a high- density model with Gamma=0.21 (TCDM) and a low-density model with Omega=0.3 and Lambda=0.7 (LCDM). Both are normalized to reproduce the I-band Tully-Fisher relation near a circular velocity of 220 km/s. Our results depend strongly both on this normalization and on the prescriptions for star formation and feedback. Very different assumptions are required to obtain an acceptable model in the two cases. For TCDM, efficient feedback is required to suppress the growth of galaxies low-mass field halos. Without it, there are too many galaxies and the correlation function turns over below 1 Mpc. For LCDM, feedback must be weak, otherwise too few L* galaxies are produced and the correlation function is too steep. Given the uncertainties in modelling some of the key physical processes, we conclude that it is not yet possible to draw conclusions about the values of cosmological parameters from studies of this kind. Further work on global star formation and feedback effects is required to narrow the range of possibilities
△ Less
Submitted 21 May, 1998;
originally announced May 1998.
-
Peculiar Velocities of Galaxy Clusters
Authors:
J. M. Colberg,
S. D. M. White,
T. J. MacFarland,
A. Jenkins,
F. R. Pearce,
C. S. Frenk,
P. A. Thomas,
H. M. P. Couchman
Abstract:
We investigate the peculiar velocities predicted for galaxy clusters by theories in the cold dark matter family. A widely used hypothesis identifies rich clusters with high peaks of a suitably smoothed version of the linear density fluctuation field. Their peculiar velocities are then obtained by extrapolating the similarly smoothed linear peculiar velocities at the positions of these peaks. We…
▽ More
We investigate the peculiar velocities predicted for galaxy clusters by theories in the cold dark matter family. A widely used hypothesis identifies rich clusters with high peaks of a suitably smoothed version of the linear density fluctuation field. Their peculiar velocities are then obtained by extrapolating the similarly smoothed linear peculiar velocities at the positions of these peaks. We test these ideas using large high resolution N-body simulations carried out within the Virgo supercomputing consortium. We find that at early times the barycentre of the material which ends up in a rich cluster is generally very close to a high peak of the initial density field. Furthermore the mean peculiar velocity of this material agrees well with the linear value at the peak. The late-time growth of peculiar velocities is, however, systematically underestimated by linear theory. At the time clusters are identified we find their rms peculiar velocity to be about 40% larger than predicted. Nonlinear effects are particularly important in superclusters. These systematics must be borne in mind when using cluster peculiar velocities to estimate the parameter combination $σ_8Ω^{0.6}$.
△ Less
Submitted 6 May, 1998;
originally announced May 1998.
-
Linking Cluster Formation to Large Scale Structure
Authors:
J. M. Colberg,
S. D. M. White,
A. Jenkins,
F. R. Pearce
Abstract:
We use two high resolution CDM simulations to show that (i) when clusters of galaxies form the infall pattern of matter is not random but shows clear features which are correlated in time; (ii) in addition, the infall patterns are correlated with the cluster's surrounding Large Scale Structure; (iii) Large Scale Structure shows a mix of both filaments and sheets; (iv) the amount of mass in filam…
▽ More
We use two high resolution CDM simulations to show that (i) when clusters of galaxies form the infall pattern of matter is not random but shows clear features which are correlated in time; (ii) in addition, the infall patterns are correlated with the cluster's surrounding Large Scale Structure; (iii) Large Scale Structure shows a mix of both filaments and sheets; (iv) the amount of mass in filaments is slightly larger for a low $Ω$ model.
△ Less
Submitted 5 November, 1997;
originally announced November 1997.
-
Genus Statistics of the Virgo N-body simulations and the 1.2-Jy Redshift Survey
Authors:
V. Springel,
S. D. M. White,
J. M. Colberg,
H. M. P. Couchman,
G. P. Efstathiou,
C. S. Frenk,
A. R. Jenkins,
F. R. Pearce,
A. H. Nelson,
J. A. Peacock,
P. A. Thomas
Abstract:
We study the topology of the Virgo N-body simulations and compare it to the 1.2-Jy redshift survey of IRAS galaxies by means of the genus statistic. Four high-resolution simulations of variants of the CDM cosmology are considered: a flat standard model (SCDM), a variant of it with more large-scale power (tCDM), and two low density universes, one open (OCDM) and one flat (LCDM). The fully sampled…
▽ More
We study the topology of the Virgo N-body simulations and compare it to the 1.2-Jy redshift survey of IRAS galaxies by means of the genus statistic. Four high-resolution simulations of variants of the CDM cosmology are considered: a flat standard model (SCDM), a variant of it with more large-scale power (tCDM), and two low density universes, one open (OCDM) and one flat (LCDM). The fully sampled N-body simulations are examined down to strongly nonlinear scales, both with spatially fixed smoothing, and with an adaptive smoothing technique. While the tCDM, LCDM, and OCDM simulations have very similar genus statistics in the regime accessible to fixed smoothing, they can be separated with adaptive smoothing at small mass scales. In order to compare the N-body models with the 1.2-Jy survey, we extract large ensembles of mock catalogues from the simulations. These mock surveys are used to test for systematic effects in the genus analysis and to establish the distribution of errors of the genus curve. We find that a simple multivariate analysis of the genus measurements is compromised both by non-Gaussian distributed errors and by noise that dominates the covariance matrix. We therefore introduce a principal components analysis of the genus curve. With a likelihood ratio test we find that the 1.2-Jy data favours the LCDM, tCDM and OCDM models compared to SCDM. When genus measurements for different smoothing scales are combined, the SCDM model can be excluded at a 99% confidence level, while the other three models fit the 1.2-Jy data well. (abridged)
△ Less
Submitted 31 October, 1997;
originally announced October 1997.
-
On Arcs and Omega
Authors:
Matthias Bartelmann,
Andreas Huss,
Joerg M. Colberg,
Adrian Jenkins,
Frazer R. Pearce
Abstract:
The gravitational lens effect of galaxy clusters can produce large arcs from source galaxies in their background. Typical source redshifts of ~ 1 require clusters at z ~ 0.3 for arcs to form efficiently. Given the cluster abundance at the present epoch, the fewer clusters exist at z ~ 0.3 the higher Omega_0 is, because the formation epoch of galaxy clusters strongly depends on Omega_0. In additi…
▽ More
The gravitational lens effect of galaxy clusters can produce large arcs from source galaxies in their background. Typical source redshifts of ~ 1 require clusters at z ~ 0.3 for arcs to form efficiently. Given the cluster abundance at the present epoch, the fewer clusters exist at z ~ 0.3 the higher Omega_0 is, because the formation epoch of galaxy clusters strongly depends on Omega_0. In addition, at fixed Omega_0, clusters are less concentrated, and hence less efficient lenses, when the cosmological constant is positive, Omega_Lambda > 0. Numerical cluster simulations show that the expected number of arcs on the sky is indeed a sensitive function of Omega_0 and Omega_Lambda. The numerical results are compatible with the statistics of observed arcs only in a universe with low matter density, Omega_0 ~ 0.3, and zero cosmological constant. Other models fail by one or two orders of magnitude, rendering arc statistics a sensitive probe for cosmological parameters.
△ Less
Submitted 24 September, 1997;
originally announced September 1997.
-
Evolution of structure in cold dark matter universes
Authors:
The Virgo Consortium,
:,
A. Jenkins,
C. S. Frenk,
F. R. Pearce,
P. A. Thomas,
J. M. Colberg,
S. D. M. White,
H. M. P. Couchman,
J. A. Peacock,
G. Efstathiou,
A. H. Nelson
Abstract:
We present an analysis of the clustering evolution of dark matter in four cold dark matter (CDM) cosmologies. We use a suite of high resolution, 17-million particle, N-body simulations which sample volumes large enough to give clustering statistics with unprecedented accuracy. We investigate both a flat and an open model with Omega_0=0.3, and two models with Omega=1, one with the standard CDM po…
▽ More
We present an analysis of the clustering evolution of dark matter in four cold dark matter (CDM) cosmologies. We use a suite of high resolution, 17-million particle, N-body simulations which sample volumes large enough to give clustering statistics with unprecedented accuracy. We investigate both a flat and an open model with Omega_0=0.3, and two models with Omega=1, one with the standard CDM power spectrum and the other with the same power spectrum as the Omega_0=0.3 models. The amplitude of primordial fluctuations is set so that the models reproduce the observed abundance of rich galaxy clusters by the present day. The mass 2-point correlation function and power spectrum of all the simulations differ significantly from those of the observed galaxy distribution, in both shape and amplitude. Thus, for any of these models to provide an acceptable representation of reality, the distribution of galaxies must be biased relative to the mass in a non-trivial, scale-dependent, fashion. In the Omega=1 models the required bias is always greater than unity, but in the Omega_0=0.3 models an "antibias" is required on scales smaller than \sim 5\hmpc. The mass correlation functions in the simulations are well fit by recently published analytic models. The velocity fields are remarkably similar in all the models, whether they be characterised as bulk flows, 1-particle or pairwise velocity dispersions. This similarity is a direct consequence of our adopted normalisation. The small-scale pairwise velocity dispersion of the dark matter is somewhat larger than recent determinations from galaxy redshift surveys, but the bulk-flows predicted by our models are broadly in agreement with most available data.
△ Less
Submitted 9 January, 1998; v1 submitted 2 September, 1997;
originally announced September 1997.
-
Arc statistics with realistic cluster potentials. IV. Clusters in different cosmologies
Authors:
Matthias Bartelmann,
Andreas Huss,
Joerg M. Colberg,
Adrian Jenkins,
Frazer R. Pearce
Abstract:
We use numerical simulations of galaxy clusters in different cosmologies to study their ability to form large arcs. The cosmological models are: Standard CDM (SCDM; Omega_0=1, Omega_Lambda=0); tauCDM with reduced small-scale power (parameters as SCDM, but with a smaller shape parameter of the power spectrum); open CDM (OCDM; Omega_0=0.3, Omega_Lambda=0); and spatially flat, low-density CDM (Lamb…
▽ More
We use numerical simulations of galaxy clusters in different cosmologies to study their ability to form large arcs. The cosmological models are: Standard CDM (SCDM; Omega_0=1, Omega_Lambda=0); tauCDM with reduced small-scale power (parameters as SCDM, but with a smaller shape parameter of the power spectrum); open CDM (OCDM; Omega_0=0.3, Omega_Lambda=0); and spatially flat, low-density CDM (LambdaCDM; Omega_0=0.3, Omega_Lambda=0.7). All models are normalised to the local number density of rich clusters. Simulating gravitational lensing by these clusters, we compute optical depths for the formation of large arcs. For large arcs with length-to-width ratio >= 10, the optical depth is largest for OCDM. Relative to OCDM, the optical depth is lower by about an order of magnitude for LambdaCDM, and by about two orders of magnitude for S/tauCDM. These differences originate from the different epochs of cluster formation across the cosmological models, and from the non-linearity of the strong lensing effect. We conclude that only the OCDM model can reproduce the observed arc abundance well, while the other models fail to do so by orders of magnitude.
△ Less
Submitted 15 July, 1997;
originally announced July 1997.
-
The structure of galaxy clusters in different cosmologies
Authors:
Virgo Consortium,
:,
Peter A. Thomas,
J. M. Colberg,
H. M. P. Couchman,
G. P. Efstathiou,
C. S. Frenk,
A. R. Jenkins,
A. H. Nelson,
R. M. Hutchings,
J. A. Peacock,
F. R. Pearce,
S. D. M. White
Abstract:
We investigate the internal structure of clusters of galaxies in high-resolution N-body simulations of 4 different cosmologies. There is a higher proportion of disordered clusters in critical-density than in low-density universes, although the structure of relaxed clusters is very similar in each. Crude measures of substructure, such as the shift in the position of the centre-of-mass as the dens…
▽ More
We investigate the internal structure of clusters of galaxies in high-resolution N-body simulations of 4 different cosmologies. There is a higher proportion of disordered clusters in critical-density than in low-density universes, although the structure of relaxed clusters is very similar in each. Crude measures of substructure, such as the shift in the position of the centre-of-mass as the density threshold is varied, can distinguish the two in a sample of just 20 or so clusters; it is harder to differentiate between clusters in open and flat models with the same density parameter. Most clusters are in a quasi-steady state within the virial radius and are well-described by the density profile of Navarro, Frenk & White (1995).
△ Less
Submitted 2 July, 1997;
originally announced July 1997.
-
The $Ω$ Dependence in the Equations of Motion
Authors:
Adi Nusser,
Jörg M. Colberg
Abstract:
We show that the equations of motion governing the evolution of a collisionless gravitating system of particles in an expanding universe can be cast in a form which is almost independent of the cosmological density parameter, $Ω$, and the cosmological constant, $Λ$. The new equations are expressed in terms of a time variable $τ\equiv \ln D$, where $D$ is the linear rate of growth of density fluc…
▽ More
We show that the equations of motion governing the evolution of a collisionless gravitating system of particles in an expanding universe can be cast in a form which is almost independent of the cosmological density parameter, $Ω$, and the cosmological constant, $Λ$. The new equations are expressed in terms of a time variable $τ\equiv \ln D$, where $D$ is the linear rate of growth of density fluctuations. The weak dependence on the density parameter is proportional to $ε=Ω^{-0.2}-1$ times the difference between the peculiar velocity (with respect to $τ$) of particles and the gravity field. In the general case, the effect of this weak $Ω$ dependence is to enhance the rate of evolution of density perturbations in dense regions. In a flat universe with $Λ\ne 0$, this enhancement is less pronounced than in an open universe with $Λ=0$ and the same $Ω$. Using the spherical collapse model, we find that the increase of the $rms$ density fluctuations in a low $Ω$ universe relative to that in a flat universe with the same linear normalization is $\sim 0.01 ε(Ω) < δ^3 >$, where $δ$ is the density field in the flat universe. The equations predict that the smooth average velocity field scales like $Ω^{0.6}$ while the local velocity dispersion (rms value) scales, approximately, like $Ω^{0.5}$. High resolution N-body simulations confirm these results and show that density fields, when smoothed on scales slightly larger than clusters, are insensitive to the cosmological model. Halos in an open model simulation are more concentrated than halos of the same $M/Ω$ in a flat model simulation.
△ Less
Submitted 16 May, 1997;
originally announced May 1997.
-
The Virgo Consortium: The evolution and formation of galaxy clusters
Authors:
J. M. Colberg,
S. D. M. White,
A. Jenkins,
F. R. Pearce,
C. S. Frenk,
P. A. Thomas,
R. Hutchings,
H. M. P. Couchman,
J. A. Peacock,
G. P. Efstathiou,
A. H. Nelson
Abstract:
We report on work done by the Virgo consortium, an international collaboration set up in order to study the formation and evolution of Large Scale Structure using N-body simulations on the latest generation of parallel supercomputers. We show results of 256^3 particle simulations of the formation of clusters in four Dark Matter models with different cosmological parameters. Normalizing the model…
▽ More
We report on work done by the Virgo consortium, an international collaboration set up in order to study the formation and evolution of Large Scale Structure using N-body simulations on the latest generation of parallel supercomputers. We show results of 256^3 particle simulations of the formation of clusters in four Dark Matter models with different cosmological parameters. Normalizing the models such that one obtains the correct abundance of rich clusters yields an interesting result: The peculiar velocities of the clusters are almost independent of Ω, and depend only weakly on Γ, the shape parameter of the power spectrum. Thus, it is nearly impossible to distinguish between high and low Ωmodels on the basis of the peculiar velocities.
△ Less
Submitted 10 February, 1997;
originally announced February 1997.
-
The Virgo consortium: simulations of dark matter and galaxy clustering
Authors:
The Virgo Consortium,
:,
A. Jenkins,
C. S. Frenk,
F. R. Pearce,
P. A. Thomas,
R. Hutchings,
J. M. Colberg,
S. D. M. White,
H. M. P. Couchman,
J. A. Peacock,
G. P. Efstathiou,
A. H. Nelson
Abstract:
We report on work in progress by the Virgo consortium, a collaboration set up to carry out large simulations of the formation of galaxies and large-scale structure exploiting the latest generation of parallel supercomputers. We show results of $256^3$ particle N-body simulations of the clustering evolution of dark matter in four cold dark matter models with different cosmological parameters. The…
▽ More
We report on work in progress by the Virgo consortium, a collaboration set up to carry out large simulations of the formation of galaxies and large-scale structure exploiting the latest generation of parallel supercomputers. We show results of $256^3$ particle N-body simulations of the clustering evolution of dark matter in four cold dark matter models with different cosmological parameters. The high resolution and large volume of these simulations allows us to determine reliably the mass autocorrelation function for pair separations in the range $40\hkpc$ to $20\hmpc$. Comparison of these with the observed galaxy correlation function shows that for any of these models to be viable, the distribution of galaxies must be biased relative to the distribution of mass in a non-trivial, scale-dependent fashion. In particular, low $Ω_0$ models require the galaxies to be more ``weakly'' clustered than the mass at small and intermediate pair separations. Simulations which include the evolution of gas show that cold gas knots form with approximately the abundance expected on theoretical grounds, although a few excessively massive objects grow near the centres of rich clusters. The locations where these cold gas knots form are, in general, biased relative to the distribution of mass in a scale-dependent way. Some of these biases have the required sign but they are, for the most part, weaker than is necessary for agreement with observations. The antibias present in our low $Ω_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.
△ Less
Submitted 25 October, 1996;
originally announced October 1996.