The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investigate the redshift evolution of metallicity, with particular focus on the cluster outskirts. The early enrichment model, in which the majority of metals are produced in the core of cluster progenitors at high redshift, suggests that metals in cluster outskirts have not significantly evolved since $z=2$. With the C-EAGLE sample, we find reasonable agreement with the early enrichment model as there is very little scatter in the metallicity abundance at large radius across the whole sample, out to at least $z=2$. The exception is Fe for which the radial dependence of metallicity was found to evolve at low redshift as a result of being mainly produced by Type Ia supernovae, which are more likely to be formed at later times than core-collapse supernovae. We also found considerable redshift evolution of metal abundances in the cores of the C-EAGLE clusters which has not been seen in other simulations or observation based metallicity studies. Since we find this evolution to be driven by accretion of low metallicity gas, it suggests that the interaction between outflowing, AGN heated material and the surrounding gas is important for determining the core abundances in clusters.

Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alignment of the halo or galaxy shape with its own mass. We report that backsplash galaxies, i.e. objects that have already passed through the cluster radius but are now located in the outskirts, show a stronger radial alignment than infalling objects. We further find that there exists a population of well radially aligned objects passing very close to the central cluster's centre which were found to be on highly radial orbit.

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

Upcoming wide-field surveys are well-suited to studying the growth of galaxy clusters by tracing galaxy and gas accretion along cosmic filaments. We use hydrodynamic simulations of volumes surrounding 324 clusters from \textsc{The ThreeHundred} project to develop a framework for identifying and characterising these filamentary structures, and associating galaxies with them. We define 3-dimensional reference filament networks reaching $5R_{200}$ based on the underlying gas distribution and quantify their recovery using mock galaxy samples mimicking observations such as those of the WEAVE Wide-Field Cluster Survey. Since massive galaxies trace filaments, they are best recovered by mass-weighting galaxies or imposing a bright limit (e.g. $>L^*$) on their selection. We measure the transverse gas density profile of filaments, derive a characteristic filament radius of $\simeq0.7$--$1~h^{-1}\rm{Mpc}$, and use this to assign galaxies to filaments. For different filament extraction methods we find that at $R>R_{200}$, $\sim15$--$20%$ of galaxies with $M_*>3 \times 10^9 M_{\odot}$ are in filaments, increasing to $\sim60%$ for galaxies more massive than the Milky-Way. The fraction of galaxies in filaments is independent of cluster mass and dynamical state, and is a function of cluster-centric distance, increasing from $\sim13$% at $5R_{200}$ to $\sim21$% at $1.5R_{200}$. As a bridge to the design of observational studies, we measure the purity and completeness of different filament galaxy selection strategies. Encouragingly, the overall 3-dimensional filament networks and $\sim67$% of the galaxies associated with them are recovered from 2-dimensional galaxy positions.

Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we compare simulations to observations in the correct way. Here we introduce the \textsc{Mock-X} analysis framework, a multiwavelength tool that generates synthetic images from cosmological simulations and derives halo properties via observational methods. We detail our methods for generating optical, Compton-$y$ and X-ray images. Outlining our synthetic X-ray image analysis method, we demonstrate the capabilities of the framework by exploring hydrostatic mass bias for the IllustrisTNG, BAHAMAS and MACSIS simulations. Using simulation derived profiles we find an approximately constant bias $b\approx0.13$ with cluster mass, independent of hydrodynamical method or subgrid physics. However, the hydrostatic bias derived from synthetic observations is mass-dependent, increasing to $b=0.3$ for the most massive clusters. This result is driven by a single temperature fit to a spectrum produced by gas with a wide temperature distribution in quasi-pressure equilibrium. The spectroscopic temperature and mass estimate are biased low by cooler gas dominating the emission, due to its quadratic density dependence. The bias and the scatter in estimated mass remain independent of the numerical method and subgrid physics. Our results are consistent with current observations and future surveys will contain sufficient samples of massive clusters to confirm the mass dependence of the hydrostatic bias.

In the outer regions of a galaxy cluster, galaxies may be either falling into the cluster for the first time, or have already passed through the cluster centre at some point in their past. To investigate these two distinct populations, we utilise TheThreeHundred project, a suite of 324 hydrodynamical resimulations of galaxy clusters. In particular, we study the 'backsplash population' of galaxies; those that have passed within $R_{200}$ of the cluster centre at some time in their history, but are now outside of this radius. We find that, on average, over half of all galaxies between $R_{200}$ and $2R_{200}$ from their host at $z=0$ are backsplash galaxies, but that this fraction is dependent on the dynamical state of a cluster, as dynamically relaxed clusters have a greater backsplash fraction. We also find that this population is mostly developed at recent times ($z\leq0.4$), and is dependent on the recent history of a cluster. Finally, we show that the dynamical state of a given cluster, and thus the fraction of backsplash galaxies in its outskirts, can be predicted based on observational properties of the cluster.

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history.

We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but still producing clusters which broadly match observations. At the same mass resolution as BAHAMAS, for the most massive clusters ($M_{500} > 10^{15}$ M$_{\odot}$), we find changes in the SPH method produce the greatest differences in the final halo, while the subgrid models dominate at lower mass. By calculating the mass of all of the clusters using different permutations of the pressure, temperature and density profiles, created with either the true simulated data or mock spectroscopic data, we find that the spectroscopic temperature causes a bias in the hydrostatic mass estimates which increases with the mass of the cluster, regardless of the SPH flavour used. For the most massive clusters, the estimated mass of the cluster using spectroscopic density and temperature profiles is found to be as low as 50 per cent of the true mass compared to $\sim$ 90 per cent for low mass clusters. When including a correction for non-thermal pressure, the spectroscopic hydrostatic mass estimates are less biased on average and the mass dependence of the bias is reduced, although the scatter in the measurements does increase.

We use TheThreeHundred project, a suite of 324 resimulated massive galaxy clusters embedded in a broad range of environments, to investigate (i) how the gas content of surrounding haloes correlates with phase-space position at $z=0$, and (ii) to investigate the role that ram pressure plays in this correlation. By stacking all 324 normalised phase-space planes containing 169287 haloes and subhaloes, we show that the halo gas content is tightly correlated with phase-space position. At $\sim\,1.5-2\,\text{R}_{\text{200}}$ of the cluster dark matter halo, we find an extremely steep decline in the halo gas content of infalling haloes and subhaloes irrespective of cluster mass, possibly indicating the presence of an accretion shock. We also find that subhaloes are particularly gas-poor, even in the cluster outskirts, which could indicate active regions of ongoing pre-processing. By modelling the instantaneous ram pressure experienced by each halo and subhalo at $z=0$, we show that the ram pressure intensity is also well correlated with phase-space position, which is again irrespective of cluster mass. In fact, we show that regions in the phase-space plane with high differential velocity between a halo or subhalo and its local gas environment, are almost mutually exclusive with high halo gas content regions. This suggests a causal link between the gas content of objects and the instantaneous ram pressure they experience, where the dominant factor is the differential velocity.

Recent numerical studies of the dark matter density profiles of massive galaxy clusters ($M_{\rm halo} > 10^{15}$M$_{\odot}$) show that their median radial mass density profile remains unchanged up to $z > 1$, displaying a highly self-similar evolution. We verify this by using the data set of the THE THREE HUNDRED project, i.e. 324 cluster-sized haloes as found in full physics hydrodynamical simulations. We track the progenitors of the mass-complete sample of clusters at $z=0$, and find that their median shape is already in place by $z=2.5$. However, selecting a dynamically relaxed subsample ($\sim16$ per cent of the clusters), we observe a shift of the scale radius $r_s$ towards larger values at earlier times. Classifying the whole sample by formation time, this evolution is understood as a result of a two-phase halo mass accretion process. Early-forming clusters -- identified as relaxed today -- have already entered their slow accretion phase, hence their mass growth occurs mostly at the outskirts. Late-forming clusters -- which are still unrelaxed today -- are in their fast accretion phase, thus the central region of the clusters is still growing. We conclude that the density profile of galaxy clusters shows a profound self-similarity out to redshifts $z\sim2.5$. This result holds for both gas and total density profiles when including baryonic physics, as reported here for two rather distinct sub-grid models.

The relationship between galaxy properties and environment is a widely discussed topic within astrophysics. Here we use galaxy samples from hydrodynamical re-simulations to examine this relationship. We use the over-density ($\delta_1$) within a $1 h^{-1}{\rm Mpc}$ sphere around a galaxy to evaluate its environment. Then the relations between galaxy properties, such as specific star formation rate(sSFR), fraction of star forming galaxies, $g-r$ colour and $\delta_1$ are examined within three galactic samples formed from galaxies within large clusters, those in the vicinity of large clusters and those in the field. We find tight environmental correlations for these galaxy properties. In brief, galaxies in denser environments tend to be redder and are more likely to be quenched. This is consistent with observations. We find that although the sSFR decreases with $\delta_1$, this is mainly because that galaxies with higher stellar mass reside in environment with higher overdensity. At fixed over-density a galaxy's color is also independent of whether it lives within a cluster or within the field, but the relative fractions of the two samples varies dramatically with over-density and this drives an apparent evolution.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON.

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

Investigating the spin parameter distribution of subhaloes in two high resolution isolated halo simulations, re- cent work by Onions et al. suggested that typical subhalo spins are consistently lower than the spin distribution found for field haloes. To further examine this puzzle, we have analyzed simulations of a cosmological volume with sufficient resolution to resolve a significant subhalo population. We confirm the result of Onions et al. and show that the typical spin of a subhalo decreases with decreasing mass and increasing proximity to the host halo center. We interpret this as the growing influence of tidal stripping in removing the outer layers, and hence the higher angular momentum particles, of the subhaloes as they move within the host potential. Investigating the redshift dependence of this effect, we find that the typical subhalo spin is smaller with decreasing redshift. This indicates a temporal evolution as expected in the tidal stripping scenario.

A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

We study the systematic bias introduced when selecting the spectroscopic redshifts of brighter cluster galaxies to estimate the velocity dispersion of galaxy clusters from both simulated and observational galaxy catalogues. We select clusters with Ngal > 50 at five low redshift snapshots from a semi-analytic model galaxy catalogue, and from a catalogue of SDSS DR8 groups and clusters across the redshift range 0.021<z<0.098. We employ various selection techniques to explore whether the velocity dispersion bias is simply due to a lack of dynamical information or is the result of an underlying physical process occurring in the cluster, for example, dynamical friction. The velocity dispersions and stacked particle velocity distributions of the parent dark matter (DM) halos are compared to the corresponding cluster dispersions and galaxy velocity distribution. We find a clear bias between the halo and the semi-analytic galaxy cluster velocity dispersion on the order of sigma gal / sigma DM = 0.87-0.95 and a distinct difference in the stacked galaxy and DM particle velocity distribution. We identify a systematic underestimation of the velocity dispersions when imposing increasing absolute I-band magnitude limits. This underestimation is enhanced when using only the brighter cluster members for dynamical analysis on the order of 5-35%, indicating that dynamical friction is a serious source of bias when using galaxy velocities as tracers of the underlying gravitational potential. In contrast to the literature we find that the resulting bias is not only halo mass-dependent but that the nature of the dependence changes according to the galaxy selection strategy. We make a recommendation that, in the realistic case of limited availability of spectral observations, a strictly magnitude-limited sample should be avoided to ensure an unbiased estimate of the velocity dispersion.

There is strong evidence that supermassive black holes reside in all galaxies that contain a stellar spheroid and their mass is tightly correlated with properties such as stellar bulge mass and velocity dispersion. There are also strong theoretical arguments that feedback from supermassive black holes plays an important role in shaping the high mass end of the galaxy mass function, hence to accurately model galaxies we also need to model the black holes. We present a comparison of two black hole growth models implemented within a large-scale, cosmological SPH simulation including star formation and feedback. One model is a modified Bondi-Hoyle prescription that grows black holes based on the smooth density of local gas, while the other is the recently proposed Accretion Disc Particle (ADP) method. This model swallows baryonic particles that pass within an accretion radius of the black hole and adds them to a subgrid accretion disc. Black holes are then grown by material from this disc. We find that both models can reproduce local scaling relations, although the ADP model is offset from the observed relations at high black hole masses. The total black hole mass density agrees between models to within a factor of three, but both struggle to reproduce the black hole mass function. The simulated mass functions are too steep and underestimate the number of intermediate and high mass black holes. In addition, the ADP model swallows excessive amounts of material at the resolution of large-scale, cosmological simulations producing unrealistically large accretion discs. Future work needs to be performed to improve the black hole mass function within simulations. This should be done through the mass growth and feedback as they are strongly coupled and should not be treated as separate entities.

Galaxy environment is frequently discussed, but inconsistently defined. It is especially difficult to measure at high redshift where only photometric redshifts are available. With a focus on early forming proto-clusters, we use a semi-analytical model of galaxy formation to show how the environment measurement around high redshift galaxies is sensitive to both scale and metric, as well as to cluster viewing angle, evolutionary state, and the availability of either spectroscopic or photometric data. We use two types of environment metrics (nearest neighbour and fixed aperture) at a range of scales on simulated high-z clusters to see how "observed" overdensities compare to "real" overdensities. We also "observationally" identify z = 2 proto-cluster candidates in our model and track the growth histories of their parent halos through time, considering in particular their final state at z = 0. Although the measured environment of early forming clusters is critically dependent on all of the above effects (and in particular the viewing angle), we show that such clusters are very likely (< 90%) to remain overdense at z = 0, although many will no longer be among the most massive. Object to object comparisons using different methodologies and different data, however, require much more caution.

While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations (Knebe et al., 2011, Onions et al., 2012), this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three non-tracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite (Springel et al., 2008) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the {\em unbound tidal features} associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of $\sim20%$. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of $\sim40%$. Analysing the results of the substructure finders, we find that the general distribution of {\em substructures} differ significantly from the distribution of bound {\em subhaloes}. Most importantly, both bound and unbound {\em substructures} together constitute $\sim18%$ of the host halo mass, which is a factor of $\sim2$ higher than the fraction in self-bound {\em subhaloes}. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history.

We analyze environmental correlations using mark clustering statistics with the mock galaxy catalogue constructed by Muldrew et al. (Paper I). We find that mark correlation functions are able to detect even a small dependence of galaxy properties on the environment, quantified by the overdensity $1+\delta$, while such a small dependence would be difficult to detect by traditional methods. We then show that rank ordering the marks and using the rank as a weight is a simple way of comparing the correlation signals for different marks. With this we quantify to what extent fixed-aperture overdensities are sensitive to large-scale halo environments, nearest-neighbor overdensities are sensitive to small-scale environments within haloes, and colour is a better tracer of overdensity than is luminosity.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

We have exploited the large-volume Millennium Gas cosmological N-body hydrodynamics simulations to study the SZ cluster population at low and high redshift, for three models with varying gas physics. We confirm previous results using smaller samples that the intrinsic (spherical) Y_{500}-M_{500} relation has very little scatter (sigma_{log_{10}Y}~0.04), is insensitive to cluster gas physics and evolves to redshift one in accord with self-similar expectations. Our pre-heating and feedback models predict scaling relations that are in excellent agreement with the recent analysis from combined Planck and XMM-Newton data by the Planck Collaboration. This agreement is largely preserved when r_{500} and M_{500} are derived using the hydrostatic mass proxy, Y_{X,500}, albeit with significantly reduced scatter (sigma_{log_{10}Y}~0.02), a result that is due to the tight correlation between Y_{500} and Y_{X,500}. Interestingly, this assumption also hides any bias in the relation due to dynamical activity. We also assess the importance of projection effects from large-scale structure along the line-of-sight, by extracting cluster Y_{500} values from fifty simulated 5x5 square degree sky maps. Once the (model-dependent) mean signal is subtracted from the maps we find that the integrated SZ signal is unbiased with respect to the underlying clusters, although the scatter in the (cylindrical) Y_{500}-M_{500} relation increases in the pre-heating case, where a significant amount of energy was injected into the intergalactic medium at high redshift. Finally, we study the hot gas pressure profiles to investigate the origin of the SZ signal and find that the largest contribution comes from radii close to r_{500} in all cases. The profiles themselves are well described by generalised Navarro, Frenk & White profiles but there is significant cluster-to-cluster scatter.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

We use numerical simulations to investigate, for the first time, the joint effect of feedback from supernovae (SNe) and active galactic nuclei (AGN) on the evolution of galaxy cluster X-ray scaling relations. Our simulations are drawn from the Millennium Gas Project and are some of the largest hydrodynamical N-body simulations ever carried out. Feedback is implemented using a hybrid scheme, where the energy input into intracluster gas by SNe and AGN is taken from a semi-analytic model of galaxy formation. This ensures that the source of feedback is a population of galaxies that closely resembles that found in the real universe. We show that our feedback model is capable of reproducing observed local X-ray scaling laws, at least for non-cool core clusters, but that almost identical results can be obtained with a simplistic preheating model. However, we demonstrate that the two models predict opposing evolutionary behaviour. We have examined whether the evolution predicted by our feedback model is compatible with observations of high-redshift clusters. Broadly speaking, we find that the data seems to favour the feedback model for z<0.5, and the preheating model at higher redshift. However, a statistically meaningful comparison with observations is impossible, because the large samples of high-redshift clusters currently available are prone to strong selection biases. As the observational picture becomes clearer in the near future, it should be possible to place tight constraints on the evolution of the scaling laws, providing us with an invaluable probe of the physical processes operating in galaxy clusters.

The Millennium Gas project aims to undertake smoothed-particle hydrodynamic resimulations of the Millennium Simulation, providing many hundred massive galaxy clusters for comparison with X-ray surveys (170 clusters with kTsl > 3 keV). This paper looks at the hot gas and stellar fractions of clusters in simulations with different physical heating mechanisms. These fail to reproduce cool-core systems but are successful in matching the hot gas profiles of non-cool-core clusters. Although there is immense scatter in the observational data, the simulated clusters broadly match the integrated gas fractions within r500 . In line with previous work, however, they fare much less well when compared to the stellar fractions, having a dependence on cluster mass that is much weaker than is observed. The evolution with redshift of the hot gas fraction is much larger in the simulation with early preheating than in one with continual feedback; observations favour the latter model. The strong dependence of hot gas fraction on cluster physics limits its use as a probe of cosmological parameters.

With the ever increasing resolution of N-body simulations, accurate subhalo detection is becoming essential in the study of the formation of structure, the production of merger trees and the seeding of semi-analytic models. To investigate the state of halo finders, we compare two different approaches to detecting subhaloes; the first based on overdensities in a halo and the second being adaptive mesh refinement. A set of stable mock NFW dark matter haloes were produced and a subhalo was placed at different radii within a larger halo. SUBFIND (a Friends-of-Friends based finder) and AHF (an adaptive mesh based finder) were employed to recover the subhalo. As expected, we found that the mass of the subhalo recovered by SUBFIND has a strong dependence on the radial position and that neither halo finder can accurately recover the subhalo when it is very near the centre of the halo. This radial dependence is shown to be related to the subhalo being truncated by the background density of the halo and originates due to the subhalo being defined as an overdensity. If the subhalo size is instead determined using the peak of the circular velocity profile, a much more stable value is recovered. The downside to this is that the maximum circular velocity is a poor measure of stripping and is affected by resolution. For future halo finders to recover all the particles in a subhalo, a search of phase space will need to be introduced.

We investigate the assembly of groups and clusters of galaxies using the Millennium dark matter simulation and the associated gas simulations and semi-analytic catalogues of galaxies. In particular, in order to find an observable quantity that could be used to identify early-formed groups, we study the development of the difference in magnitude between their brightest galaxies to assess the use of magnitude gaps as possible indicators. We select galaxy groups and clusters at redshift z=1 with dark matter halo mass M(R200) > 1E13/h Msun, and trace their properties until the present time (z=0). We consider only the systems with X-ray luminosity L_X> 0.25E42/h^2 erg/s at z=0. While it is true that a large magnitude gap between the two brightest galaxies of a particular group often indicates that a large fraction of its mass was assembled at an early epoch, it is not a necessary condition. More than 90% of fossil groups defined on the basis of their magnitude gaps (at any epoch between 0<z<1) cease to be fossils within 4 Gyr, mostly because other massive galaxies are assembled within their cores, even though most of the mass in their haloes might have been assembled at early times. We show that, compared to the conventional definition of fossil galaxy groups based on the magnitude gap Delta m(12)> 2 (in the R-band, within 0.5R200 of the centre of the group), an alternative criterion Delta m(14)>2.5 (within the same radius) finds 50% more early-formed systems, and those that on average retain their fossil phase longer. However, the conventional criterion performs marginally better at finding early-formed groups at the high-mass end of groups. Nevertheless, both criteria fail to identify a majority of the early-formed systems.

The joint likelihood of observable cluster signals reflects the astrophysical evolution of the coupled baryonic and dark matter components in massive halos, and its knowledge will enhance cosmological parameter constraints in the coming era of large, multi-wavelength cluster surveys. We present a computational study of intrinsic covariance in cluster properties using halo populations derived from Millennium Gas Simulations (MGS). The MGS are re-simulations of the original 500 Mpc/h Millennium Simulation performed with gas dynamics under two different physical treatments: shock heating driven by gravity only (GO) and a second treatment with cooling and preheating (PH). We examine relationships among structural properties and observable X-ray and Sunyaev-Zel'dovich (SZ) signals for samples of thousands of halos with M_200 > 5 \times 10^{13} Msun/h and z < 2. While the X-ray scaling behavior of PH model halos at low-redshift offers a good match to local clusters, the model exhibits non-standard features testable with larger surveys, including weakly running slopes in hot gas observable--mass relations and ~10% departures from self-similar redshift evolution for 10^14 Msun/h halos at redshift z ~ 1. We find that the form of the joint likelihood of signal pairs is generally well-described by a multivariate, log-normal distribution, especially in the PH case which exhibits less halo substructure than the GO model. At fixed mass and epoch, joint deviations of signal pairs display mainly positive correlations, especially the thermal SZ effect paired with either hot gas fraction (r=0.88/0.69 for PH/GO at z=0) or X-ray temperature (r=0.62/0.83). We discuss halo mass selection by signal pairs, and find a minimum mass scatter of 4% in the \PH model by combining thermal SZ and gas fraction measurements.

[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of "downsizing" (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.

We test four commonly used astrophysical simulation codes; Enzo, Flash, Gadget and Hydra, using a suite of numerical problems with analytic initial and final states. Situations similar to the conditions of these tests, a Sod shock, a Sedov blast and both a static and translating King sphere occur commonly in astrophysics, where the accurate treatment of shocks, sound waves, supernovae explosions and collapsed haloes is a key condition for obtaining reliable validated simulations. We demonstrate that comparable results can be obtained for Lagrangian and Eulerian codes by requiring that approximately one particle exists per grid cell in the region of interest. We conclude that adaptive Eulerian codes, with their ability to place refinements in regions of rapidly changing density, are well suited to problems where physical processes are related to such changes. Lagrangian methods, on the other hand, are well suited to problems where large density contrasts occur and the physics is related to the local density itself rather than the local density gradient.

The evolution of present-day fossil galaxy groups is studied in the Millennium Simulation. Using the corresponding Millennium gas simulation and semi-analytic galaxy catalogues, we select fossil groups at redshift zero according to the conventional observational criteria, and trace the haloes corresponding to these groups backwards in time, extracting the associated dark matter, gas and galaxy properties. The space density of the fossils from this study is remarkably close to the observed estimates and various possibilities for the remaining discrepancy are discussed. The fraction of X-ray bright systems which are fossils appears to be in reasonable agreement with observation, and the simulations predict that fossil systems will be found in significant numbers (3-4% of the population) even in quite rich clusters. We find that fossils assemble a higher fraction of their mass at high redshift, compared to non-fossil groups, with the ratio of the currently assembled halo mass to final mass, at any epoch, being about 10 to 20% higher for fossils. This supports the paradigm whereby fossils represent undisturbed, early-forming systems in which large galaxies have merged to form a single dominant elliptical.

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.

The physical processes that define the spine of the galaxy cluster X-ray luminosity -- temperature (L-T) relation are investigated using a large hydrodynamical simulation of the Universe. This simulation models the same volume and phases as the Millennium Simulation and has a linear extent of 500 h^{-1} Mpc. We demonstrate that mergers typically boost a cluster along but also slightly below the L-T relation. Due to this boost we expect that all of the very brightest clusters will be near the peak of a merger. Objects from near the top of the L-T relation tend to have assembled much of their mass earlier than an average halo of similar final mass. Conversely, objects from the bottom of the relation are often experiencing an ongoing or recent merger.

(Abridged) The thermal history of the intracluster medium (ICM) is complex. Heat input from cluster mergers, AGN, and galaxy winds offsets and may even halt the cooling of the ICM. Consequently, the processes that set the properties of the ICM play a key role in determining how galaxies form. In this paper we focus on the shock heating of the ICM during cluster mergers, with the eventual aim of incorporating this mechanism into semi-analytic models of galaxy formation. We use a suite of hydrodynamic simulations to track the evolution of the ICM in idealised two-body mergers. We find the heating of the ICM can be understood relatively simply by considering the evolution of the gas entropy during the mergers. We examine the processes that generate the entropy in order to understand why previous analytic shock heating models failed. We find that: (1) The energy that is thermalised in the collision greatly exceeds the kinetic energy available when the systems first touch. The smaller system penetrates deep into the potential well before it is disrupted. (2) For unequal mass mergers, most of the energy is thermalised in the more massive component. The heating of the smaller system is minor and its gas sinks to the centre of the final system. (3) The bulk of the entropy generation occurs in two distinct episodes. The first episode occurs following core collision, when a shock wave is generated that propagates outwards from the centre. This causes the combined system to expand rapidly and overshoot hydrostatic equilibrium. The second episode occurs as this material is shock heated as it re-collapses. This revised model for entropy generation significantly improves our physical understanding of cosmological gas simulations.

We follow the growth of baryonic structure in the presence of a magnetic field within an approximate cosmological magneto-hydrodynamic simulation, produced by adding an (isotropic) magnetic pressure related to the local gas pressure. We perform an ensemble of these simulations to follow the amplification of the field with time. By using a variety of initial field strengths and changing the slope of the power law that governs the way the field grows with increasing density we span the range of current observations and demonstrate the size of the effect realistic magnetic fields could have on the central density of groups and clusters. A strong magnetic field significantly reduces the central gas density which, in turn, reduces observable quantities such as the X-ray luminosity.

The comparison between observations of galaxy clusters thermo-dynamical properties and theoretical predictions suggests that non-gravitational heating needs to be added into the models. We implement an internally self-consistent heating scheme into GADGET-2 for the third (and fourth) run of the Millennium gas project (Pearce et al. in preparation), a set of four hydrodynamical cosmological simulations with N=2(5x10^8) particles and with the same volume (L=500 h-1 Mpc) and structures as the the N-body Millennium Simulation (Springel et al. 2005). Our aim is to reproduce the observed thermo-dynamical properties of galaxy clusters.

Using the Millennium N-body simulation we explore how the shape and angular momentum of galaxy dark matter haloes surrounding the largest cosmological voids are oriented. We find that the major and intermediate axes of the haloes tend to lie parallel to the surface of the voids, whereas the minor axis points preferentially in the radial direction. We have quantified the strength of these alignments at different radial distances from the void centres. The effect of these orientations is still detected at distances as large as 2.2 R_void from the void centre. Taking a subsample of haloes expected to contain disc-dominated galaxies at their centres we detect, at the 99.9% confidence level, a signal that the angular momentum of those haloes tends to lie parallel to the surface of the voids. Contrary to the alignments of the inertia axes, this signal is only detected in shells at the void surface (1<R<1.07 R_void) and disappears at larger distances. This signal, together with the similar alignment observed using real spiral galaxies (Trujillo, Carretero & Patiri 2006), strongly supports the prediction of the Tidal Torque theory that both dark matter haloes and baryonic matter have acquired, conjointly, their angular momentum before the moment of turnaround.

We investigate the baryon fraction in dark matter haloes formed in non-radiative gas-dynamical simulations of the LambdaCDM cosmogony. By combining a realisation of the Millennium Simulation (Springel et al.) with a simulation of a smaller volume focussing on dwarf haloes, our study spans five decades in halo mass, from 10^10 Msun/h to 10^15 Msun/h. We find that the baryon fraction within the halo virial radius is typically 90% of the cosmic mean, with an rms scatter of 6%, independently of redshift and of halo mass down to the smallest resolved haloes. Our results show that, contrary to the proposal of Mo et al. (2005), pre-virialisation gravitational heating is unable to prevent the collapse of gas within galactic and proto-galactic haloes, and confirm the need for non-gravitational feedback in order to reduce the efficiency of gas cooling and star formation in dwarf galaxy haloes. Simulations including a simple photoheating model (where a gas temperature floor of T_{floor} = 2x10^4 K is imposed from z=11) confirm earlier suggestions that photoheating can only prevent the collapse of baryons in systems with virial temperatures T_{200} < ~2.2 T_{floor} ~ 4.4x10^4 K (corresponding to a virial mass of M_{200} ~ 10^10 Msun/h and a circular velocity of V_{200} ~ 35 km/s). Photoheating may thus help regulate the formation of dwarf spheroidals and other galaxies at the extreme faint-end of the luminosity function, but it cannot, on its own, reconcile the abundance of sub-L* galaxies with the vast number of dwarf haloes expected in the LambdaCDM cosmogony. The lack of evolution or mass dependence seen in the baryon fraction augurs well for X-ray cluster studies that assume a universal and non-evolving baryon fraction to place constraints on cosmological parameters.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

We present deep, large area B and r' imaging for a sample of 49 brightest cluster galaxies (BCGs). The clusters were selected by their x-ray luminosity and redshift to form two volume limited samples, one with mean redshift ~ 0.07 and one at a mean redshift ~ 0.17. For each cluster the data cover 41' by 41'. We discuss our data reduction techniques in detail, and show that we can reliably measure the surface brightness at the levels of mu_B ~ 29 and mu_r' ~ 28. For each galaxy we present the B and r' images together with the surface brightness profile, B-r' colour, eccentricity and position angle as a function of radius. We investigate the distribution of positional offsets between the optical centroid of the BCG and the centre of the X-ray emission, and conclude that the mass profiles are cuspy, and do not have extended cores. We also introduce a method to objectively identify the transition from BCG to extended envelope of intra-cluster light, using the Petrosian index as a function of radius.

Galaxy groups have been under-studied relative to their richer counterparts -- clusters. The Group Evolution Multiwavelength Study (GEMS) aims to redress some the balance. Here we describe the GEMS sample selection and resulting sample of 60 nearby (distance < 130 Mpc) galaxy groups and our multiwavelength dataset of X-ray, optical and HI imaging. ROSAT X-ray images of each group are presented. GEMS also utilizes near-infrared imaging from the 2MASS survey and optical spectra from the 6dFGS. These observational data are complemented by mock group catalogues generated from the latest LCDM simulations with gas physics included. Existing GEMS publications are briefly highlighted as are future publication plans.

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 present predictions for the angular correlation function and redshift distribution for SHADES, the SCUBA HAlf-Degree Extragalactic Survey, which will yield a sample of around 300 sub-mm sources in the 850 micron waveband in two separate fields. Complete and unbiased photometric redshift information on these sub-mm sources will be derived by combining the SCUBA data with i) deep radio imaging already obtained with the VLA, ii) guaranteed-time Spitzer data at mid-infrared wavelengths, and iii) far-infrared maps to be produced by BLAST, the Balloon-borne Large-Aperture Sub-millimeter Telescope. Predictions for the redshift distribution and clustering properties of the final anticipated SHADES sample have been computed for a wide variety of models, each constrained to fit the observed number counts. Since we are dealing with around 150 sources per field, we use the sky-averaged angular correlation function to produce a more robust fit of a power-law shape w(theta)=(theta/A)^{-delta} to the model data. Comparing the predicted distributions of redshift and of the clustering amplitude A and slope delta, we find that models can be constrained from the combined SHADES data with the expected photometric redshift information.

We investigate the properties of the intracluster medium (ICM) that forms within N-body/hydrodynamical simulations of galaxy clusters in a \LambdaCDM cosmology. When radiative cooling and a simple model for galactic feedback are included, our clusters have X-ray luminosities and temperatures in good agreement with observed systems, demonstrating the required excess entropy in their cores. More generally, cooling and feedback increases the entropy of the ICM everywhere, albeit without significantly affecting the slope of the profile (S prop. to r) at large radii. The temperature of the ICM is only modestly increased by these processes, with projected temperature profiles being in reasonable agreement with the observations. Star/galaxy formation is still too efficient in our simulations, however, and so our gas mass fractions are around 60 per cent of the observed value at r_{2500}. Finally, we examine the reliability of using the hydrostatic equilibrium equation to estimate cluster masses and find that it underpredicts the true mass of our clusters by up to 20 per cent, due to incomplete thermalisation of the gas. Feedback reduces this discrepancy, however, with estimates being accurate to within 10 per cent out to r_{500}.

By utilising the inbuilt citation counts from NASA's astrophysics data system (ADS) I derive how many citations refereed articles receive as a function of time since publication. After five years, one paper in a hundred has accumulated 91 or more citations, a figure which rises to 145 citations after ten years. By adding up the number of citations active researchers have received over the past five years I have estimated their relative impact upon the field both for raw citations and citations weighted by the number of authors per paper.

We explore a scenario where metal-poor globular clusters (GCs) are enriched by the first supernovae in the Universe. If the first stars in a 10^7 Msun dark halo were very massive (>180 Msun), then a pair instability supernova from a single massive star can produce sufficient iron to enrich 10^6 Msun of pristine, primordial gas to [Fe/H] ~ -2. In such a scenario, where a single massive star acts as a seed for halo GCs, the accurate abundance analysis of GC stars would allow a direct measurement of the Population III initial mass. Using the latest theoretical yields for zero metallicity stars in the mass range 140-260 Msun, we find that the metals expelled from a ~230 Msun star are consistent with [Si/Fe] and [Ca/Fe] observed in GC stars. However, no single star in this mass range can simultaneously explain all halo GC heavy-element abundance ratios, such as [V/Fe], [Ti/Fe] and [Ni/Fe]. These require a combination masses for the Population III stellar progenitors. The various observational consequences of this scenario are discussed.

We present the results of a large library of cosmological N-body simulations, using power-law initial spectra. The nonlinear evolution of the matter power spectra is compared with the predictions of existing analytic scaling formulae based on the work of Hamilton et al. The scaling approach has assumed that highly nonlinear structures obey `stable clustering' and are frozen in proper coordinates. Our results show that, when transformed under the self-similarity scaling, the scale-free spectra define a nonlinear locus that is clearly shallower than would be required under stable clustering. Furthermore, the small-scale nonlinear power increases as both the power-spectrum index n and the density parameter Omega decrease, and this evolution is not well accounted for by the previous scaling formulae. This breakdown of stable clustering can be understood as resulting from the modification of dark-matter haloes by continuing mergers. These effects are naturally included in the analytic `halo model' for nonlinear structure; using this approach we are able to fit both our scale-free results and also our previous CDM data. This approach is more accurate than the commonly-used Peacock--Dodds formula and should be applicable to more general power spectra. Code to evaluate nonlinear power spectra using this method is available from http://as1.chem.nottingham.ac.uk/~res/software.html Following publication, we will make the power-law simulation data available through the Virgo website http://www.mpa-garching.mpg.de/Virgo

We compare the results of two techniques used to calculate the evolution of cooling gas during galaxy formation: Smooth Particle Hydrodynamics (SPH) simulations and semi-analytic modelling. We improve upon the earlier statistical comparison of Benson et al. by taking halo merger histories from the dark matter component of the SPH simulation, which allows us to compare the evolution of galaxies on an object-by-object basis in the two treatments. We use a ``stripped-down'' version of the semi-analytic model described by Helly et al. which includes only shock heating and radiative cooling of gas and which is adjusted to mimic the resolution and other parameters of a comparison SPH simulation as closely as possible. We compare the total mass of gas that cools in halos of different mass as a function of redshift as well as the masses and spatial distribution of individual ``galaxies.'' At redshift z=0, the cooled gas mass in well-resolved halos agrees remarkably well (to better than ~20%) in the SPH simulation and stripped-down semi-analytic model. At high redshift, resolution effects in the simulation become increasingly important and, as a result, more gas tends to cool in low mass halos in the SPH simulation than in the semi-analytic model. The cold gas mass function of individual galaxies in the two treatments at z=0 also agrees very well and, when the effects of mergers are accounted for, the masses of individual galaxies and their 2-point correlation functions are also in excellent agreement in the two treatments. Thus, our comparison confirms and extends the earlier conclusion of Benson et al. that SPH simulations and semi-analytic models give consistent results for the evolution of cooling galactic gas.

We calculate X-ray properties of present-day galaxy clusters from hydrodynamical cosmological simulations of the LCDM cosmology and compare these with recent X-ray observations. Results from three simulations are presented, each of which uses the same initial conditions: a standard adiabatic, Non-radiative model, a Radiative model that includes radiative cooling of the gas, and a Preheating model that also includes cooling but in addition impulsively heats the gas prior to cluster formation. At the end of the simulations, the global cooled baryon fractions in the latter two runs are 15 per cent and 0.4 per cent respectively which bracket the recent result from the K-band luminosity function. We construct cluster catalogues which consist of over 500 clusters and are complete in mass down to 1.18*10^{13} Msun/h. While clusters in the Non-radiative model behave in accord with the self-similar picture, those of the other two models reproduce key aspects of the observed X-ray properties: the core entropy, temperature-mass and luminosity-temperature relations are all in good agreement with recent observations. This agreement stems primarily from an increase in entropy with respect to the Non-radiative clusters. Although the physics affecting the intra-cluster medium is very different in the two models, the resulting cluster entropy profiles are very similar.

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 investigate phenomenological models of star formation and supernova feedback in N-body/SPH simulations of galaxy formation. First, we compare different prescriptions in the literature for turning cold gas into stars neglecting feedback effects. We find that most prescriptions give broadly similar results: the ratio of cold gas to stars in the final galaxies is primarily controlled by the range of gas densities where star formation is allowed to proceed efficiently. In the absence of feedback, the fraction of gas that cools is much too high resulting, for example, in a K-band luminosity function that is much brighter than observed. This problem is ameliorated by including a feedback model which either imparts radial kinetic perturbations to galactic gas or directly reheats such material and prevents it from cooling for a certain period of time. In both these models, a significant fraction of cold gas is heated and expelled from haloes with an efficiency that varies inversely with with halo circular velocity. Increasing the resolution of a simulation allows a wider dynamic range in mass to be followed, but the average properties of the resolved galaxy population remain largely unaffected. However, as the resolution is increased, more and more gas is reheated by small galaxies; our results suggest that convergence requires the full mass range of galaxies to be resolved.

We use hydrodynamical simulations to assess the impact of radiative cooling and `pre-heating' on predictions for the Sunyaev--Zel'dovich (SZ) effect. Cooling significantly reduces both the mean SZ signal and its angular power spectrum, while pre-heating can give a higher mean distortion while leaving the angular power spectrum below that found in a simulation without heating or cooling. We study the relative contribution from high and low density gas, and find that in the cooling model about 60 per cent of the mean thermal distortion arises from low overdensity gas. We find that haloes dominate the thermal SZ power spectrum in all models, while in the cooling simulation the kinetic SZ power spectrum originates predominantly in lower overdensity gas.

Recent measurements of the K-band luminosity function now provide us with strong, reliable constraints on the fraction of baryons which have cooled. Globally, this fraction is only about 5%, and there is no strong evidence that it is significantly higher in clusters. Without an effective sub-grid feedback prescription, the cooled gas fraction in any numerical simulation exceeds these observational constraints, and increases with increasing resolution. This compromises any discussion of galaxy and cluster properties based on results of simulations which include cooling but do not implement an effective feedback mechanism.

We have performed cosmological simulations in a $\Lambda$CDM cosmology with and without radiative cooling, in order to study the effect of cooling on the cluster scaling relations. Our simulations consist of 4.1 million particles each of gas and dark matter within a box-size of 100 $h^{-1}$ Mpc and the run with cooling is the largest of its kind to have been evolved to $z=0$. Our cluster catalogues both consist of over 400 objects and are complete in mass down to $\sim 10^{13} h^{-1} {\rm M_{\odot}}$. We contrast the emission-weighted temperature-mass ($T_{\rm ew}-M$) and bolometric luminosity-temperature ($L_{\rm bol}-T_{\rm ew}$) relations for the simulations at $z=0$. We find that radiative cooling \emph{increases} the temperature of intracluster gas and \emph{decreases} its total luminosity, in agreement with the results of Pearce et al. Furthermore, the temperature dependence of these effects flattens the slope of the $T_{\rm ew}-M$ relation and steepens the slope of the $L_{\rm bol}-T_{\rm ew}$ relation. Inclusion of radiative cooling in the simulations is sufficient to reproduce the observed X-ray scaling relations without requiring excessive non-gravitational energy injection.

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 present results of large N-body-hydrodynamic simulations of galaxy formation. Our simulations follow the formation of galaxies in cubic volumes of side 100Mpc, in two versions of the cold dark matter (CDM) cosmogony: the standard, Omega=1 SCDM model and the flat, Omega=0.3 LCDM model. Over 2000 galaxies form in each of these simulations. We examine the rate at which gas cools and condenses into dark matter halos. This roughly tracks the cosmic star formation rate inferred from observations at various redshifts. Galaxies in the simulations form gradually over time in the hierarchical fashion characteristic of the CDM cosmogony. In the LCDM model, substantial galaxies first appear at z=5 and the population builds up rapidly until z=1 after which the rate of galaxy formation declines as cold gas is consumed and the cooling time of hot gas increases. In both cosmologies, the present-day K-band luminosity function of the simulated galaxies resembles observations. The galaxy autocorrelation functions differ significantly from those of the dark matter. At the present epoch there is little bias in either model between galaxies and dark matter on large scales, but a significant anti-bias on scales of 1Mpc/h and a positive bias on scales of 100kpc/h. The galaxy correlation function evolves little with redshift in the range z=0-3, and depends on the luminosity of the galaxy sample. The projected pairwise velocity dispersion of the galaxies is much lower than that of the dark matter on scales less than 2Mpc/h. Applying a virial mass estimator to the largest galaxy clusters recovers the cluster virial masses in an unbiased way. (Abridged)

We compare the statistical properties of galaxies found in two different models of hierarchical galaxy formation: the semi-analytic model of Cole et al. and the smoothed particle hydrodynamics (SPH) simulations of Pearce et al. Using a `stripped-down' version of the semi-analytic model which mimics the resolution of the SPH simulations and excludes physical processes not included in them, we find that the two models produce an ensemble of galaxies with remarkably similar properties, although there are some differences in the gas cooling rates and in the number of galaxies that populate halos of different mass. The full semi-analytic model, which has effectively no resolution limit and includes a treatment of star formation and supernovae feedback, produces somewhat different (but readily understandable) results. Agreement is particularly good for the present-day global fractions of hot gas, cold dense (i.e. galactic) gas and uncollapsed gas, for which the SPH and stripped-down semi-analytic calculations differ by at most 25%. In the most massive halos, the stripped-down semi-analytic model predicts, on the whole, up to 50% less gas in galaxies than is seen in the SPH simulations. The two techniques apportion this cold gas somewhat differently amongst galaxies in a given halo. This difference can be tracked down to the greater cooling rate in massive halos in the SPH simulation compared to the semi-analytic model. (abridged)

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 have extracted over 400 clusters, covering more than 2 decades in mass, from three simulations of the TCDM cosmology. This represents the largest, uniform catalogue of simulated clusters ever produced. The clusters exhibit a wide variety of density-profiles. Only a minority are well-fit in their outer regions by the widely used density profile of Navarro, Frenk & White (1977). Others have steeper outer density profiles, show sharp breaks in their density profiles, or have significant substructure. If we force a fit to the NFW profile, then the best-fit concentrations decline with increasing mass, but this is driven primarily by an increase in substructure as one moves to higher masses. The measured temperature-mass relations for properties measured within a sphere enclosing a fixed overdensity all follow the self-similar form, T\propto M^{2/3}, however the normalisation is lower than in observed clusters. The temperature-mass relations for properties measured within a fixed physical radius are significantly steeper then this. Both can be accurately predicted using the NFW model.

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.

In this paper, we investigate the effect of cooling on the X-ray properties of galaxy clusters. We have performed N-body, hydrodynamical simulations both with and without the effects of radiative cooling, but neglecting the effects of star formation and feedback. We show that radiative cooling produces an inflow of high-entropy gas from the outer parts of the cluster, thus \emph{raising} the cluster temperature and \emph{decreasing} the X-ray luminosity. With radiative cooling clusters are on average three to five times less luminous in X-rays than the same cluster simulated without cooling. However, we do not produce a large constant-density core in either the gas or the dark matter distributions. Our results contradict previous work in which cooling raises the X-ray luminosity and deposits an unreasonably large amount of mass in the central cluster galaxy. We achieve this by selecting our numerical resolution in such a way that a reasonable fraction of the baryonic material cools and by decoupling the hot and cold gas in our simulations, a first step towards modelling multiphase gas. We emphasise that globally cooling a sensible amount of material is vital and the presence or absence of massive central concentrations of cold baryonic material has a dramatic effect upon the resultant X-ray properties of the clusters.

Numerical simulations of galaxy formation require a number of parameters. Some of these are intrinsic to the numerical integration scheme (eg the timestep), while others describe the physical model (eg the gas metallicity). In this paper, we present results of a systematic exploration of the effects of varying a subset of these parameters on simulations of galaxy formation. We use N-body and ``Smoothed Particle Hydrodynamics'' techniques to follow the evolution of cold dark matter and gas in a small volume. We compare a fiducial model to 24 different simulations, in which one parameter at a time is varied, focussing on properties such as the relative fraction of hot and cold gas, and the abundance and masses of galaxies. We find that for reasonable choices of numerical values, many parameters have relatively little effect on the galaxies, with the notable exception of the parameters that control the resolution of the simulation and the efficiency with which gas cools.

We compare the structure and substructure of dark matter halos in model universes dominated by collisional, strongly self interacting dark matter (SIDM) and collisionless, weakly interacting dark matter (CDM). While SIDM virialised halos are more nearly spherical than CDM halos, they can be rotationally flattened by as much as 20% in their inner regions. Substructure halos suffer ram-pressure truncation and drag which are more rapid and severe than their gravitational counterparts tidal stripping and dynamical friction. Lensing constraints on the size of galactic halos in clusters are a factor of two smaller than predicted by gravitational stripping, and the recent detection of tidal streams of stars escaping from the satellite galaxy Carina suggests that its tidal radius is close to its optical radius of a few hundred parsecs --- an order of magnitude smaller than predicted by CDM models but consistent with SIDM. The orbits of SIDM satellites suffer significant velocity bias $\sigma_{SIDM}/\sigma_{CDM}=0.85$ and are more circular than CDM, $\beta_{SIDM}} \approx 0.5$, in agreement with the inferred orbits of the Galaxy's satellites. In the limit of a short mean free path, SIDM halos have singular isothermal density profiles, thus in its simplest incarnation SIDM is inconsistent with galactic rotation curves.

The Hydra $N$-body hydrodynamics code has been modified to model, from the end of the Sedov phase, the effects of supernovae on the surrounding medium. The motivation is to investigate the feedback of energy into the interstellar/intergalactic medium. We compare our results for supernova remnants (SNRs) in a uniform medium to previous detailed work on the late evolution of SNRs. The code is found to reproduce the bulk characteristics of SNRs well. Results on the effects of a single central SNR on Plummer clouds are presented. The feedback of kinetic energy and the percentage mass loss can be parameterised in terms of the cloud mass and characteristic radius in a simple way. The kinetic energy fraction returned to the ISM from a SNR is $<3$ per cent. The removal of gas from the cold, dense phase and the addition of energy due to the lowering of the potential energy of a cloud is at least as significant, if not much more so, than the kinetic energy leaving a cloud.

In this paper we use state-of-the-art N-body hydrodynamic simulations of a cosmological volume of side 100Mpc to produce many galaxy clusters simultaneously in both the standard cold dark matter (SCDM) cosmology and a cosmology with a positive cosmological constant (LCDM). We have performed simulations of the same volume both with and without the effects of radiative cooling, but in all cases neglect the effects of star formation and feedback. With radiative cooling clusters are on average five times less luminous in X-rays than the same cluster simulated without cooling. The importance of the mass of the central galaxy in determining the X-ray luminosity is stressed.

We present a study of hydrodynamic drag forces in smoothed particle simulations. In particular, the deceleration of a resolution-limited cold clump of gas moving through a hot medium is examined. It is found that the drag for subsonic velocities exceeds that predicted by simple physical approximations. This is shown to be a result of the hydrodynamical method which encourages the accretion of particles from the hot medium onto a shell around the cold clump, effectively increasing the radius of the clump. For sonic and supersonic velocities, the drag is shown to be dependent on the effective cross section of the clump. The consequences for cosmological simulations are discussed.

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

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

We analyse the performance of twelve different implementations of Smoothed Particle Hydrodynamics (SPH) using seven tests designed to isolate key hydrodynamic elements of cosmological simulations which are known to cause the SPH algorithm problems. In order, we consider a shock tube, spherical adiabatic collapse, cooling flow model, drag, a cosmological simulation, rotating cloud-collapse and disc stability. In the implementations special attention is given to the way in which force symmetry is enforced in the equations of motion. We study in detail how the hydrodynamics are affected by different implementations of the artificial viscosity including those with a shear-correction modification. We present an improved first-order smoothing-length update algorithm that is designed to remove instabilities that are present in the Hernquist and Katz (1989) algorithm. For all tests we find that the artificial viscosity is the most important factor distinguishing the results from the various implementations. The second most important factor is the way force symmetry is achieved in the equation of motion. Most results favour a kernel symmetrization approach. The exact method by which SPH pressure forces are included has comparatively little effect on the results. Combining the equation of motion presented in Thomas and Couchman (1992) with a modification of the Monaghan and Gingold (1983) artificial viscosity leads to an SPH scheme that is both fast and reliable.

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.

(Abridged) We investigate the importance of several numerical artifacts such as lack of resolution on spectral properties of the Lyman alpha forest as computed from cosmological hydrodynamic simulations in a standard cold dark matter universe. We assume an ionising background produced by quasars as computed by Haardt & Madau. We use a new simulation code based on P3M and SPH, which we compare in detail with a modified version of HYDRA (Couchman et al.) and published results of TREESPH (Hernquist et al.). The agreement is very good between all three codes. We then use our new code to investigate several numerical effects such as resolution on spectral statistics deduced from Voigt profile fitting. Our highest resolution simulation has a mass resolution of 2.1x10^5 solar masses. The column density distribution is converged but the b-parameter distribution is only marginally converged. The simulation reproduces both the HI column density and b-parameter distribution when we assume a high baryon density, Omega_B h^2 > 0.028. In addition we need to impose a higher IGM temperature than predicted within our basic set of assumptions. The simulated HI optical depth is in good agreement with observations but the HeII optical depth is lower than observed. Fitting the latter requires a larger jump between the photon flux at the H and He edge than is present in the Haardt & Madau spectrum.

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

I investigate the process of converting gas into stars within the framework of a standard cosmological model. By examining the set of objects grown in a combined N-body plus smoothed particle hydrodynamics simulation with those obtained in similar models where some of the cold, dense gas was replaced by collisionless ``star'' particles I show that it is possible to make this substitution without affecting the subsequent gas cooling rate. With even the most basic star forming criteria the masses of isolated objects are nearly identical to the mass of cold, dense gas found within the same objects in a non-star forming run. No evidence is found to support the contention that converting gas into stars might affect the amount of cold gas obtained in a simulation by retarding the cooling rate within those objects where stars have already formed. In practice, because cold gas can be reheated by shocks but stars remain as such whatever happens the masses of the largest objects found in the star forming runs are generally higher than those in the standard run. Finally, I demonstrate that an excellent match to the observed star formation rate can be achieved with even a very basic star formation prescription.

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

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

A revised version of Hydra, an adaptive particle-particle, particle-mesh plus smoothed particle hydrodynamics code developed by; Hugh Couchman, University of Western Ontario, Frazer Pearce, University of Durham, Peter Thomas, University of Sussex, is now available as a tar file from; http://coho.astro.uwo.ca/pub/hydra/hydra.html or http://star-www.maps.susx.ac.uk/~pat/hydra/hydra.html .

We describe the first parallel implementation of an adaptive particle-particle, particle-mesh code with smoothed particle hydrodynamics. Parallelisation of the serial code, ``Hydra'', is achieved by using CRAFT, a Cray proprietary language which allows rapid implementation of a serial code on a parallel machine by allowing global addressing of distributed memory. The collisionless variant of the code has already completed several 16.8 million particle cosmological simulations on a 128 processor Cray T3D whilst the full hydrodynamic code has completed several 4.2 million particle combined gas and dark matter runs. The efficiency of the code now allows parameter-space explorations to be performed routinely using $64^3$ particles of each species. A complete run including gas cooling, from high redshift to the present epoch requires approximately 10 hours on 64 processors. In this paper we present implementation details and results of the performance and scalability of the CRAFT version of Hydra under varying degrees of particle clustering.

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.

We intend to make sets of cosmological simulations available, with a large number of different output times that may be placed side-by-side to produce a complete history of the universe stretching back to high redshift. Currently there is only one series of runs under preparation as a trial. Gauging by the response, we will consider further releases. The first series consists of N-body simulations with CDM-like spectra (Gamma=0.25) in a 100/h Mpc box, in three different cosmologies: Einstein-de Sitter, low-density flat and low-density open. The parameters of each run are: 1) Run 501: Omega=1.0, lambda=0.0, sigma_8=0.64. N=128**3. 38 time-slices back to z=5.8 (a further 11 back to z=23.4 on request). 2) Run 502: Omega=0.3, lambda=0.7, sigma_8=1.22. N=86**3. 49 time-slices back to z=3.6. 3) Run 503: Omega=0.3, lambda=0.0, sigma_8=1.06. N=86**3. 49 time-slices back to z=6.2 Check out our Web pages for more details: http://coho.astro.uwo.ca/pub/consort.html

We present an implementation of Smoothed Particle Hydrodynamics (SPH) in an adaptive-mesh PPPM algorithm. The code evolves a mixture of purely gravitational particles and gas particles. The code retains the desirable properties of previous PPPM--SPH implementations; speed under light clustering, naturally periodic boundary conditions and accurate pairwise forces. Under heavy clustering the cycle time of the new code is only 2--3 times slower than for a uniform particle distribution, overcoming the principal disadvantage of previous implementations\dash a dramatic loss of efficiency as clustering develops. A 1000 step simulation with 65,536 particles (half dark, half gas) runs in one day on a Sun Sparc10 workstation. The choice of time integration scheme is investigated in detail. A simple single-step Predictor--Corrector type integrator is most efficient. A method for generating an initial distribution of particles by allowing a a uniform temperature gas of SPH particles to relax within a periodic box is presented. The average SPH density that results varies by $\sim\pm1.3$\%. We present a modified form of the Layzer--Irvine equation which includes the thermal contribution of the gas together with radiative cooling. Tests of sound waves, shocks, spherical infall and collapse are presented. Appropriate timestep constraints sufficient to ensure both energy and entropy conservation are discussed. A cluster simulation, repeating Thomas and

We present a practical guide to Smoothed Particle Hydrodynamics (\SPH) and its application to astrophysical problems. Although remarkably robust, \SPH\ must be used with care if the results are to be meaningful since the accuracy of \SPH\ is sensitive to the arrangement of the particles and the form of the smoothing kernel. In particular, the initial conditions for any \SPH\ simulation must consist of particles in dynamic equilibrium. We describe some of the numerical difficulties that may be encountered when using \SPH, and how these may be overcome. Through our experience in using \SPH\ code to model convective stars, galaxy clusters and large scale structure problems we have developed many diagnostic tests. We give these here as an aid to rapid identification of errors, together with a list of basic prerequisites for the most efficient implementation of \SPH.

Several simple mergers between model galaxy clusters containing a mixture of gas and dark matter are examined, testing the coupling of the gas to the underlying collisionless material. The gas is shocked, irreversibly dissipating the energy fed into it by the collisionless component and forms a resolved constant-density core. For the dark matter, however, admixture of phase space vacuum is not very efficient and a constant-density core is not produced. In the final state the central gas has little residual kinetic energy, indicating that streaming motions do not help to support the gas.

Orbits in triaxial ellipsoidal potentials precess about either the major or minor axis of the ellipsoid. In standard perturbation theory it can be shown that a circular orbit will precess about the minor axis if its angular momentum vector lies in a region bounded by two great circles which pass through the intermediate axis and which are inclined with minimum separation $i_T$ from the minor axis. We test the accuracy of the standard formula for $i_T$ by performing orbit integrations to determine $i_S$, the simulated turnover angle corresponding to $i_T$. We reach two principal conclusions: (i) $i_S$ is usually greater than $i_T$, by as much as 12 degrees even for moderate triaxialities, $A/1.2<B<C/0.8$. This reduces the expected frequency of polar rings. (ii) $i_S$ is not a single, well-defined number but can vary by a few degrees depending upon the initial phase of the orbit. This means that there is a reasonable probability for capture of gas onto orbits which precess about both axes. Interactions can then lead to substantial loss of angular momentum and subsequent infall to the galactic centre.