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

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 have developed a parallel Particle-Particle, Particle-Mesh (P3M) simulation code for the Cray T3E parallel supercomputer that is well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in cosmological contexts. The parallel code is based upon the public-domain serial Adaptive P3M code of Couchman et al. The algorithm resolves gravitational forces into a long-range component computed by discretizing the mass distribution and solving Poisson's equation on a grid using an FFT convolution method, and a short-range component computed by direct force summation for sufficiently close particle pairs. The code consists primarily of a particle-particle computation parallelized by domain decomposition over blocks of neighbour-cells, a more regular mesh calculation distributed in planes along one dimension, and several transformations between the two distributions. The load balancing of the P3M code is static, since this greatly aids the ongoing implementation of parallel adaptive refinements of the particle and mesh systems. Great care was taken throughout to make optimal use of the available memory, so that a version of the current implementation has been used to simulate systems of up to 10^9 particles with a 1024^3 mesh for the long-range force computation. These are the largest Cosmological N-body simulations of which we are aware. We discuss these memory optimizations as well as those motivated by computational performance. Performance results are very encouraging, and, even without refinements, the code has been used effectively for simulations in which the particle distribution becomes highly clustered as well as for other non-uniform systems of astrophysical interest.

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 have developed a parallel Particle-Particle, Particle-Mesh (P^3M) simulation code for the T3E well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in cosmological contexts. The parallel code is based upon the public-domain serial Adaptive P^3M code of Couchman et al(1). The algorithm resolves gravitational forces into a long range component computed by discretizing the mass distribution and solving Poisson's equation on a grid using an FFT convolution method, and a short range component computed by direct force summation for sufficiently close particle pairs. The code consists primarily of a particle-particle computation parallelized by domain decomposition over blocks of neighbor-cells, a more regular mesh calculation distributed in planes along one dimension, and several transformations between the two distributions. Great care was taken throughout to make optimal use of the available memory, so that the current implementation is capable of simulating systems approaching 10^9 particles using a 1024^3 mesh for the long range force computation. These are thus among the largest N-body simulations ever carried out. We discuss these memory optimizations as well as those motivated by computational performance. Results from production runs have been very encouraging, and even prior to the implimentation of the full adaptive scheme the code has been used effectively for simulations in which the particle distribution becomes highly clustered as well as for other non-uniform systems of astrophysical interest.