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.

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.

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.

We investigate the coupled formation and evolution of galaxies and their embedded supermassive black holes using state-of-the-art hydrodynamic simulations of cosmological structure formation. For the first time, we self-consistently follow the dark matter dynamics, radiative gas cooling, star formation, as well as black hole growth and associated feedback processes, starting directly from initial conditions appropriate for the LambdaCDM cosmology. Our modeling of the black hole physics is based on an approach we have developed in simulations of isolated galaxy mergers. Here we examine: (i) the predicted global history of black hole mass assembly (ii) the evolution of the local black hole-host mass correlations and (iii) the conditions that allow rapid growth of the first quasars, and the properties of their hosts and descendants today. We find a total black hole mass density in good agreement with observational estimates. The black hole accretion rate density peaks at lower redshift and evolves more strongly at high redshift than the star formation rate density, but the ratio of black hole to stellar mass densities shows only a moderate evolution at low redshifts. We find strong correlations between black hole masses and properties of the stellar systems, agreeing well with the measured local M_BH-sigma and M_BH -M_* relationships, but also suggesting (dependent on the mass range) a weak evolution with redshift in the normalization and the slope. Our simulations also produce massive black holes at high redshift, due to extended periods of exponential growth in regions that collapse early and exhibit strong gas inflows. These first supermassive BH systems however are not necessarily the most massive ones today, since they are often overtaken in growth by quasars that form later. (abridged)

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.

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 $\Lambda$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.

The cold dark matter model has become the leading theoretical paradigm for the formation of structure in the Universe. Together with the theory of cosmic inflation, this model makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a novel framework for the quantitative physical interpretation of such surveys. This combines the largest simulation of the growth of dark matter structure ever carried out with new techniques for following the formation and evolution of the visible components. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with next generation surveys.

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular $\Lambda$ 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 $\rho(<r)/\rho(r_{\rm eff}) \propto \exp[(r/r_{\rm eff})^\alpha]$ where $r_{\rm eff}$ is the effective radius of a void and $\alpha\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)

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 $\Lambda$ 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)

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.

We create mock pencil-beam redshift surveys from very large cosmological $N$-body simulations of two Cold Dark Matter cosmogonies, an Einstein-de Sitter model ($\tau$CDM) and a flat model with $\Omega_0 =0.3$ and a cosmological constant ($\Lambda$CDM). We use these to assess the significance of the apparent periodicity discovered by Broadhurst et al. (1990). Simulation particles are tagged as `galaxies' so as to reproduce observed present-day correlations. They are then identified along the past light-cones of hypothetical observers to create mock catalogues with the geometry and the distance distribution of the Broadhurst et al. data. We produce 1936 (2625) quasi-independent catalogues from our $\tau$CDM ($\Lambda$CDM) simulation. A couple of large clumps in a catalogue can produce a high peak at low wavenumbers in the corresponding one-dimensional power spectrum, without any apparent large-scale periodicity in the original redshift histogram. Although the simulated redshift histograms frequently display regularly spaced clumps, the spacing of these clumps varies between catalogues and there is no `preferred' period over our many realisations. We find only a 0.72 (0.49) per cent chance that the highest peak in the power spectrum of a $\tau$CDM ($\Lambda$CDM) catalogue has a peak-to-noise ratio higher than that in the Broadhurst et al. data. None of the simulated catalogues with such high peaks shows coherently spaced clumps with a significance as high as that of the real data. We conclude that in CDM universes, the kind of regularity observed by Broadhurst et al. has a priori probability well below $10^{-3}$.

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.

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

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 ($\tau$CDM) and a flat low density model with a cosmological constant ($\Lambda$CDM). Our simulations each use $10^9$ particles to follow the mass distribution within cubes of side $2h^{-1}$Gpc ($\tau$CDM) and $3h^{-1}$Gpc ($\Lambda$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 $\tau$CDM and the upper to $\Lambda$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 $\Lambda$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.

Colombi et al. 1999 (paper I) investigated the counts-in-cells statistics and their respective errors in the $\tau$CDM Virgo Hubble Volume simulation. This extremely large $N$-body experiment also allows a numerical investigation of the {\em cosmic distribution function}, $\Upsilon(\tA)$ itself for the first time. For a statistic $A$, $\Upsilon(\tA)$ is the probability density of measuring the value $\tA$ in a finite galaxy catalog. $\Upsilon$ was evaluated for the distribution of counts-in-cells, $P_N$, the factorial moments, $F_k$, and the cumulants, $\xiav$ and $S_N$'s, using the same subsamples as paper I. While paper I concentrated on the first two moments of $\Upsilon$, i.e. the mean, the cosmic error and the cross-correlations, here the function $\Upsilon$ is studied in its full generality, including a preliminary analysis of joint distributions $\Upsilon(\tA,\tB)$. The most significant, and reassuring result for the analyses of future galaxy data is that the cosmic distribution function is nearly Gaussian provided its variance is small. A good practical criterion for the relative cosmic error is that $\Delta A/A \la 0.2$. This means that for accurate measurements, the theory of the cosmic errors, presented by Szapudi & Colombi (1996) and Szapudi, Colombi & Bernardeau (1999), and confirmed empirically by paper I, is sufficient for a full statistical description and thus for a maximum likelihood rating of models. As the cosmic error increases, the cosmic distribution function $\Upsilon$ becomes increasingly skewed and is well described by a generalization of the lognormal distribution. The cosmic skewness is introduced as an additional free parameter. (...more in paper...)

Counts-in-cells are measured in the $\tau$CDM Virgo Hubble Volume simulation. This large N-body experiment has 10^9 particles in a cubic box of size 2000 h^{-1} Mpc. The unprecedented combination of size and resolution allows for the first time a realistic numerical analysis of the cosmic errors and cosmic correlations of statistics related to counts-in-cells measurements, such as the probability distribution function P_N itself, its factorial moments F_k and the related cumulants $\xiav$ and S_N's. These statistics are extracted from the whole simulation cube, as well as from 4096 sub-cubes of size 125 h^{-1}Mpc, each representing a virtual random realization of the local universe. The measurements and their scatter over the sub-volumes are compared to the theoretical predictions of Colombi, Bouchet & Schaeffer (1995) for P_0, and of Szapudi & Colombi (1996, SC) and Szapudi, Colombi & Bernardeau (1999a, SCB) for the factorial moments and the cumulants. The general behavior of experimental variance and cross-correlations as functions of scale and order is well described by theoretical predictions, with a few percent accuracy in the weakly non-linear regime for the cosmic error on factorial moments. (... more in paper >...)

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.

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.

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.

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.

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.

We report on analyses of cluster samples obtained from the Hubble Volume Simulations. These simulations, an $\Omega=1$ model named $\tau$CDM and a flat low $\Omega$ model with a cosmological constant ($\Lambda$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 $\tau$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 $\xi$. 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 $\xi(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 $\Lambda$CDM model.

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

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 $\sigma_8\Omega^{0.6}$.

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.

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 $\Omega$ model.

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)

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.

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.

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).

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, $\Omega$, and the cosmological constant, $\Lambda$. The new equations are expressed in terms of a time variable $\tau\equiv \ln D$, where $D$ is the linear rate of growth of density fluctuations. The weak dependence on the density parameter is proportional to $\epsilon=\Omega^{-0.2}-1$ times the difference between the peculiar velocity (with respect to $\tau$) of particles and the gravity field. In the general case, the effect of this weak $\Omega$ dependence is to enhance the rate of evolution of density perturbations in dense regions. In a flat universe with $\Lambda\ne 0$, this enhancement is less pronounced than in an open universe with $\Lambda=0$ and the same $\Omega$. Using the spherical collapse model, we find that the increase of the $rms$ density fluctuations in a low $\Omega$ universe relative to that in a flat universe with the same linear normalization is $\sim 0.01 \epsilon(\Omega) < \delta^3 >$, where $\delta$ is the density field in the flat universe. The equations predict that the smooth average velocity field scales like $\Omega^{0.6}$ while the local velocity dispersion (rms value) scales, approximately, like $\Omega^{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/\Omega$ in a flat model simulation.

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 \Omega, and depend only weakly on \Gamma, the shape parameter of the power spectrum. Thus, it is nearly impossible to distinguish between high and low \Omega models on the basis of the peculiar velocities.

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 $\Omega_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 $\Omega_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.