We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, ${\gamma}_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of ${\gamma}_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of ${\gamma}_{\rm gNFW}$.

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.

By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-density contours using the procedure of Montes \& Trujillo (2019), and find that the shape of the stellar mass distribution follows that of the total matter even more closely than observed, although their radial profiles differ substantially. The ratio between stellar and total matter density profiles in circular apertures, shows a slope close to -1, with a small dependence on the cluster's total mass. We propose an indirect method to calculate the halo mass and mass density profile from the radial profile of the intra-cluster stellar mass density.

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

The Sunyaev-Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain $>10,000$ simulated galaxy groups and clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (tSZ) effect. The third measure determines the lowest order relativistic correction to the tSZ signal, which is seeing increased observational relevance. Our procedure allows us to accurately model the relativistic SZ (rSZ) contribution and we show that a $\gtrsim 10\%-40\%$ underestimation of this rSZ cluster temperature is expected when applying standard X-ray relations. The correction also exhibits significant mass and redshift evolution, as we demonstrate here. We present the mass dependence of each temperature measure alongside their profiles and a short analysis of the temperature dispersion as derived from the aforementioned simulations. We also discuss a new relation connecting the temperature and Compton-$y$ parameter, which can be directly used for rSZ modelling. Simple fits to the obtained scaling relations and profiles are provided. These should be useful for future studies of the rSZ effect and its relevance to cluster cosmology.

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.

Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$\sigma$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters.

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.

Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function (LF) of cluster red sequence galaxies as a function of redshift. This is the first homogeneous optical/X-ray sample large enough to constrain the evolution of the luminosity function simultaneously in redshift ($0.1<z<1.05$) and cluster mass ($13.5 \le \rm{log_{10}}(M_{200crit}) \sim< 15.0$). We pay particular attention to completeness issues and the detection limit of the galaxy sample. We then apply a hierarchical Bayesian model to fit the cluster galaxy LFs via a Schecter function, including its characteristic break ($m^*$) to a faint end power-law slope ($\alpha$). Our method enables us to avoid known issues in similar analyses based on stacking or binning the clusters. We find weak and statistically insignificant ($\sim 1.9 \sigma$) evolution in the faint end slope $\alpha$ versus redshift. We also find no dependence in $\alpha$ or $m^*$ with the X-ray inferred cluster masses. However, the amplitude of the LF as a function of cluster mass is constrained to $\sim 20\%$ precision. As a by-product of our algorithm, we utilize the correlation between the LF and cluster mass to provide an improved estimate of the individual cluster masses as well as the scatter in true mass given the X-ray inferred masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, significantly improving the sensitivity of the analysis.

We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes.

Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations, designed to enable the study of high-mass galaxy formation and evolution in cluster environments. We find 69 $M_{\rm BH}(M_{\star})$ outliers at $z=0$, defined as those with $ \rm M_{BH} >10^{7} M_{\odot}$ and $\rm M_{BH}/\rm M_{\star}> 0.01$. This paper focusses on a sample of 5 extreme outliers, that have been selected based on their $\rm M_{BH}$ and $\rm M_{\star}$ values, which are comparable to the most recent estimates of observed positive outliers. This sample of 5 outliers, classified as `Black hole monster galaxies' (BMGs), was traced back in time to study their origin and evolution. In agreement with the results of previous simulations for lower-mass $\rm M_{BH}(\rm M_{\star})$ outliers, we find that these galaxies became outliers due to a combination of their early formation times and tidal stripping. For BMGs with $\rm M_{BH} > 10^9 M_{\odot}$, major mergers (with a stellar mass ratio of $\mu > 0.25$) at early times ($z>2$) precede the rapid growth of their supermassive BHs. Furthermore, the scatter in the relation between $\rm M_{BH}$ and stellar velocity dispersion, $\sigma$, correlates positively with the scatter in [Mg/Fe]($\sigma$). This indicates that the alpha enhancement of these galaxies, which is closely related to their star formation history, is related to the growth of their central BHs.

Forthcoming large-scale spectroscopic surveys will soon provide data on thousands of galaxy clusters. It is important that the systematics of the various mass estimation techniques are well understood and calibrated. We compare three different dynamical mass estimators using the C-EAGLE galaxy clusters, a set of high resolution simulations with resolved galaxies a median total mass, $M_{200c} = 10^{14.7} \, \mathrm{M_\odot}$. We quantify the bias and scatter of the Jeans, virial, and caustic mass estimators using all galaxies with a stellar mass $M_*> 10^9 \, \mathrm{M_\odot}$, both in the ideal 3D case and in the more realistic projected case. On average we find our mass estimates are unbiased, though relative to the true mass within $r_{200c}$ the scatter is large with a range of $0.09$ - $0.15$ dex. We see a slight increase in the scatter when projecting the clusters. Selecting galaxies using the same criteria, we find no significant difference in the mass bias or scatter when comparing results from hydrodynamical and dark matter only simulations. However, selecting galaxies by stellar mass reduces the bias compared to selecting by total mass. Comparing X-ray derived hydrostatic and dynamical masses, the former are ${\sim} 30$ per cent lower. We find a slight dependence between substructure, measured using two different metrics, and mass bias. In conclusion, we find that dynamical mass estimators, when averaged together, are unbiased with a scatter of $0.11 \pm 0.02$ dex when including interloper galaxies and with no prior knowledge of $r_{200c}$.

Machine learning (ML) techniques, in particular supervised regression algorithms, are a promising new way to use multiple observables to predict a cluster's mass or other key features. To investigate this approach we use the \textsc{MACSIS} sample of simulated hydrodynamical galaxy clusters to train a variety of ML models, mimicking different datasets. We find that compared to predicting the cluster mass from the $\sigma -M$ relation, the scatter in the predicted-to-true mass ratio is reduced by a factor of 4, from $0.130\pm0.004$ dex (${\simeq} 35$ per cent) to $0.031 \pm 0.001$ dex (${\simeq} 7$ per cent) when using the same, interloper contaminated, spectroscopic galaxy sample. Interestingly, omitting line-of-sight galaxy velocities from the training set has no effect on the scatter when the galaxies are taken from within $r_{200c}$. We also train ML models to reproduce estimated masses derived from mock X-ray and weak lensing analyses. While the weak lensing masses can be recovered with a similar scatter to that when training on the true mass, the hydrostatic mass suffers from significantly higher scatter of ${\simeq} 0.13$ dex (${\simeq} 35$ per cent). Training models using dark matter only simulations does not significantly increase the scatter in predicted cluster mass compared to training on simulated clusters with hydrodynamics. In summary, we find ML techniques to offer a powerful method to predict masses for large samples of clusters, a vital requirement for cosmological analysis with future surveys.

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter.

We present the Tiered Radio Extragalactic Continuum Simulation (T-RECS): a new simulation of the radio sky in continuum, over the 150 MHz-20 GHz range. T-RECS models two main populations of radio galaxies: Active Galactic Nuclei (AGNs) and Star-Forming Galaxies (SFGs), and corresponding sub-populations. Our model also includes polarized emission over the full frequency range, which has been characterised statistically for each population using the available information. We model the clustering properties in terms of probability distributions of hosting halo masses, and use lightcones extracted from a high-resolution cosmological simulation to determine the positions of haloes. This limits the sky area for the simulations including clustering to a 25deg2 field of view. We compare luminosity functions, number counts in total intensity and polarization, and clustering properties of our outputs to up-to-date compilations of data and find a very good agreement. We deliver a set of simulated catalogues, as well as the code to produce them, which can be used for simulating observations and predicting results from deep radio surveys with existing and forthcoming radio facilities, such as the Square Kilometre Array (SKA).

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.

We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions.

Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of volumes up to $(100\,\rm Mpc)^3$, while Hydrangea consists of $24$ cosmological simulations of galaxy clusters and their environment. Together they provide a statistically significant sample in the stellar mass range $10^{9.5}\,\rm M_{\odot}-10^{12.3}\,\rm M_{\odot}$, of $16,358$ galaxies. We construct IFS-like cubes and measure stellar spin parameters, $\lambda_{\rm R}$, and ellipticities, allowing us to classify galaxies into slow/fast rotators as in observations. The simulations display a primary dependence of $F_{\rm SR}$ on stellar mass, with a weak dependence on environment. At fixed stellar mass, satellite galaxies are more likely to be SRs than centrals. $F_{\rm SR}$ shows a dependence on halo mass at fixed stellar mass for central galaxies, while no such trend is seen for satellites. We find that $\approx 70$% of SRs at $z=0$ have experienced at least one merger with mass ratio $\ge 0.1$, with dry mergers being at least twice more common than wet mergers. Individual dry mergers tend to decrease $\lambda_{\rm R}$, while wet mergers mostly increase it. However, $30$% of SRs at $z=0$ have not experienced mergers, and those inhabit halos with median spins twice smaller than the halos hosting the rest of the SRs. Thus, although the formation paths of SRs can be varied, dry mergers and/or halos with small spins dominate.

We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range $14.0 \le \log_{10}(M_{\rm 200c}/\mathrm{M_\odot}) \le 15.4$, with their sophisticated sub-grid physics modelling and high numerical resolution (sub-kpc gravitational softening) making them ideal for this purpose. We find that selecting galaxies by their total mass results in a velocity dispersion that is 5-10 per cent higher than the dark matter particles. However, selecting galaxies by their stellar mass results in an almost unbiased ($<5$ per cent) estimator of the velocity dispersion. This result holds out to $z=1.5$ and is relatively insensitive to the choice of cluster aperture, varying by less than 5 per cent between $r_{\rm 500c}$ and $r_{\rm 200m}$. We show that the velocity bias is a function of the time spent by a galaxy inside the cluster environment. Selecting galaxies by their total mass results in a larger bias because a larger fraction of objects have only recently entered the cluster and these have a velocity bias above unity. Galaxies that entered more than $4 \, \mathrm{Gyr}$ ago become progressively colder with time, as expected from dynamical friction. We conclude that velocity bias should not be a major issue when estimating cluster masses from kinematic methods.

Using tens of thousands of halos realized in the BAHAMAS and MACSIS simulations produced with a consistent astrophysics treatment that includes AGN feedback, we validate a multi-property statistical model for the stellar and hot gas mass behavior in halos hosting groups and clusters of galaxies. The large sample size allows us to extract fine-scale mass--property relations (MPRs) by performing local linear regression (LLR) on individual halo stellar mass (${\rm M}_{\rm star}$) and hot gas mass (${\rm M}_{\rm gas}$) as a function of total halo mass (${\rm M}_{\rm halo}$). We find that: 1) both the local slope and variance of the MPRs run with mass (primarily) and redshift (secondarily); 2) the conditional likelihood, $p({\rm M}_{\rm star},\ {\rm M}_{\rm gas} | \ {\rm M}_{\rm halo}, z)$ is accurately described by a multivariate, log-normal distribution, and; 3) the covariance of ${\rm M}_{\rm star}$ and ${\rm M}_{\rm gas}$ at fixed ${\rm M}_{\rm halo}$ is generally negative, reflecting a partially closed baryon box model for high mass halos. We validate the analytical population model of Evrard et al. (2014), finding sub-percent accuracy in the log-mean halo mass selected at fixed property, $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm gas} \rangle$ or $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm star} \rangle$, when scale-dependent MPR parameters are employed. This work highlights the potential importance of allowing for running in the slope and scatter of MPRs when modeling cluster counts for cosmological studies. We tabulate LLR fit parameters as a function of halo mass at $z=0$, $0.5$ and 1 for two popular mass conventions.

We use a combination of full hydrodynamic and dark matter only simulations to investigate the effect that baryonic physics and selecting super-cluster regions have on the matter power spectrum, by re-simulating a sample of super-cluster sub-volumes. On large scales we find that the matter power spectrum measured from our super-cluster sample has at least twice as much power as that measured from our random sample. Our investigation of the effect of baryonic physics on the matter power spectrum is found to be in agreement with previous studies and is weaker than the selection effect over the majority of scales. In addition, we investigate the effect of targeting a cosmologically non-representative, super-cluster region of the sky on the weak lensing shear power spectrum. We do this by generating shear and convergence maps using a line of sight integration technique, which intercepts our random and super-cluster sub-volumes. We find the convergence power spectrum measured from our super-cluster sample has a larger amplitude than that measured from the random sample at all scales. We frame our results within the context of the Super-CLuster Assisted Shear Survey (Super-CLASS), which aims to measure the cosmic shear signal in the radio band by targeting a region of the sky that contains five Abell clusters. Assuming the Super-CLASS survey will have a source density of 1.5 galaxies/arcmin$^2$, we forecast a detection significance of $2.7^{+1.5}_{-1.2}$, which indicates that in the absence of systematics the Super-CLASS project could make a cosmic shear detection with radio data alone.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bah\'e et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of $1.8\times10^{6}\,\mathrm{M}_{\odot}$ and physical softening length of $0.7\,\mathrm{kpc}$. In this paper, we introduce the sample and present the low-redshift global properties of the clusters. We calculate the X-ray properties in a manner consistent with observational techniques, demonstrating the bias and scatter introduced by using estimated masses. We find the total stellar content and black hole masses of the clusters to be in good agreement with the observed relations. However, the clusters are too gas rich, suggesting that the AGN feedback model is not efficient enough at expelling gas from the high-redshift progenitors of the clusters. The X-ray properties, such as the spectroscopic temperature and the soft-band luminosity, and the Sunyaev-Zel'dovich properties are in reasonable agreement with the observed relations. However, the clusters have too high central temperatures and larger-than-observed entropy cores, which is likely driven by the AGN feedback after the cluster core has formed. The total metal content and its distribution throughout the ICM are a good match to the observations.

Cosmological parameter constraints from observations of time-delay lenses are becoming increasingly precise. However, there may be significant bias and scatter in these measurements due to, among other things, the so-called mass-sheet degeneracy. To estimate these uncertainties, we analyze strong lenses from the largest EAGLE hydrodynamical simulation. We apply a mass-sheet transformation to the radial density profiles of lenses, and by selecting lenses near isothermality, we find that the bias on H0 can be reduced to 5% with an intrinsic scatter of 10%, confirming previous results performed on a different simulation data set. We further investigate whether combining lensing observables with kinematic constraints helps to minimize this bias. We do not detect any significant dependence of the bias on lens model parameters or observational properties of the galaxy, but depending on the source--lens configuration, a bias may still exist. Cross lenses provide an accurate estimate of the Hubble constant, while fold (double) lenses tend to be biased low (high). With kinematic constraints, double lenses show bias and intrinsic scatter of 6% and 10%, respectively, while quad lenses show bias and intrinsic scatter of 0.5% and 10%, respectively. For lenses with a reduced $\chi^2 > 1$, a power-law dependence of the $\chi^2$ on the lens environment (number of nearby galaxies) is seen. Lastly, we model, in greater detail, the cases of two double lenses that are significantly biased. We are able to remove the bias, suggesting that the remaining biases could also be reduced by carefully taking into account additional sources of systematic uncertainty.

We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses Mstar > 10^10 M_Sun this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (Mstar < 10^10 M_Sun) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching. Conversely, the simulations predict an overabundance of massive galaxies in clusters compared to the field that persists to their far outskirts (> 5r_200c). This is caused by a significantly increased stellar mass fraction of (sub-)haloes in the cluster environment, by up to ~0.3 dex even well beyond r_200c. Haloes near clusters are also more concentrated than equally massive field haloes, but these two effects are largely uncorrelated.

We present the MAssive ClusterS and Intercluster Structures (MACSIS) project, a suite of 390 clusters simulated with baryonic physics that yields realistic massive galaxy clusters capable of matching a wide range of observed properties. MACSIS extends the recent BAHAMAS simulation to higher masses, enabling robust predictions for the redshift evolution of cluster properties and an assessment of the effect of selecting only the hottest systems. We study the observable-mass scaling relations and the X-ray luminosity-temperature relation over the complete observed cluster mass range. As expected, we find the slope of these scaling relations and the evolution of their normalization with redshift departs significantly from the self-similar predictions. However, for a sample of hot clusters with core-excised temperatures $k_{\rm{B}}T\geq5\,\rm{keV}$ the normalization and slope of the observable-mass relations and their evolution are significantly closer to self-similar. The exception is the temperature-mass relation, for which the increased importance of non-thermal pressure support and biased X-ray temperatures leads to a greater departure from self-similarity in the hottest systems. As a consequence, these also affect the slope and evolution of the normalization in the luminosity-temperature relation. The median hot gas profiles show good agreement with observational data at $z=0$ and $z=1$, with their evolution again departing significantly from the self-similar prediction. However, selecting a hot sample of clusters yields profiles that evolve significantly closer to the self-similar prediction. In conclusion, our results show that understanding the selection function is vital for robust calibration of cluster properties with mass and redshift.

We use the BAHAMAS and MACSIS hydrodynamic simulations to quantify the impact of baryons on the mass distribution and dynamics of massive galaxy clusters, as well as the bias in X-ray and weak lensing mass estimates. These simulations use the sub-grid physics models calibrated in the BAHAMAS project, which include feedback from both supernovae and active galactic nuclei. They form a cluster population covering almost two orders of magnitude in mass, with more than 3,500 clusters with masses greater than $10^{14}\,\mathrm{M}_\odot$ at $z=0$. We start by characterising the clusters in terms of their spin, shape and density profile, before considering the bias in both weak lensing and hydrostatic mass estimates. Whilst including baryonic effects leads to more spherical, centrally concentrated clusters, the median weak lensing mass bias is unaffected by the presence of baryons. In both the dark matter only and hydrodynamic simulations, the weak lensing measurements underestimate cluster masses by ${\approx}10\%$ for clusters with $M_{200}{\leq}10^{15}\mathrm{M}_\odot$ and this bias tends to zero at higher masses. We also consider the hydrostatic bias when using both the true density and temperature profiles, and those derived from X-ray spectroscopy. When using spectroscopic temperatures and densities, the hydrostatic bias decreases as a function of mass, leading to a bias of ${\approx}40\%$ for clusters with $M_{500}{\geq}10^{15}\,\mathrm{M}_\odot$. This is due to the presence of cooler gas in the cluster outskirts. Using mass weighted temperatures and the true density profile reduces this bias to $5{-}15\%$.

We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43 <z <0.7$) galaxies is inferred via their association with the dark-matter halos of 174 X-ray-selected galaxy clusters drawn from the XMM Cluster Survey (XCS). Halo masses are determined for each galaxy cluster based on X-ray temperature measurements, and range between ${\rm log_{10}} (M_{180}/M_{\odot}) = 13-15$. Our directly measured HODs are consistent with the HOD-model fits inferred via the galaxy-clustering analyses of Parejko et al. for the BOSS LOWZ sample and White et al. for the BOSS CMASS sample. Under the simplifying assumption that the other parameters that describe the HOD hold the values measured by these authors, we have determined a best-fit alpha-index of 0.91$\pm$0.08 and $1.27^{+0.03}_{-0.04}$ for the CMASS and LOWZ HOD, respectively. These alpha-index values are consistent with those measured by White et al. and Parejko et al. In summary, our study provides independent support for the HOD models assumed during the development of the BOSS mock-galaxy catalogues that have subsequently been used to derive BOSS cosmological constraints.

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.

We measure the evolution of the velocity dispersion--temperature ($\sigma_{\rm v}$--$T_{\rm X}$) relation up to $z = 1$ using a sample of 38 galaxy clusters drawn from the \textit{XMM} Cluster Survey. This work improves upon previous studies by the use of a homogeneous cluster sample and in terms of the number of high redshift clusters included. We present here new redshift and velocity dispersion measurements for 12 $z > 0.5$ clusters observed with the GMOS instruments on the Gemini telescopes. Using an orthogonal regression method, we find that the slope of the relation is steeper than that expected if clusters were self-similar, and that the evolution of the normalisation is slightly negative, but not significantly different from zero ($\sigma_{\rm v} \propto T^{0.86 \pm 0.14} E(z)^{-0.37 \pm 0.33}$). We verify our results by applying our methods to cosmological hydrodynamical simulations. The lack of evolution seen in our data is consistent with simulations that include both feedback and radiative cooling.

We describe updates to the \redmapper{} algorithm, a photometric red-sequence cluster finder specifically designed for large photometric surveys. The updated algorithm is applied to $150\,\mathrm{deg}^2$ of Science Verification (SV) data from the Dark Energy Survey (DES), and to the Sloan Digital Sky Survey (SDSS) DR8 photometric data set. The DES SV catalog is locally volume limited, and contains 786 clusters with richness $\lambda>20$ (roughly equivalent to $M_{\rm{500c}}\gtrsim10^{14}\,h_{70}^{-1}\,M_{\odot}$) and $0.2<z<0.9$. The DR8 catalog consists of 26311 clusters with $0.08<z<0.6$, with a sharply increasing richness threshold as a function of redshift for $z\gtrsim 0.35$. The photometric redshift performance of both catalogs is shown to be excellent, with photometric redshift uncertainties controlled at the $\sigma_z/(1+z)\sim 0.01$ level for $z\lesssim0.7$, rising to $\sim0.02$ at $z\sim0.9$ in DES SV. We make use of \emph{Chandra} and \emph{XMM} X-ray and South Pole Telescope Sunyaev-Zeldovich data to show that the centering performance and mass--richness scatter are consistent with expectations based on prior runs of \redmapper{} on SDSS data. We also show how the \redmapper{} \photoz{} and richness estimates are relatively insensitive to imperfect star/galaxy separation and small-scale star masks.

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.

Using the science verification data of the Dark Energy Survey (DES) for a new sample of 106 X-Ray selected clusters and groups, we study the stellar mass growth of Bright Central Galaxies (BCGs) since redshift 1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncertainties associated with cluster mass, redshift, and BCG stellar mass measurements into analysis of a redshift-dependent BCG-cluster mass relation, $m_{*}\propto(\frac{M_{200}}{1.5\times 10^{14}M_{\odot}})^{0.24\pm 0.08}(1+z)^{-0.19\pm0.34}$, and compare the observed relation to the model prediction. We estimate the average growth rate since $z = 1.0$ for BCGs hosted by clusters of $M_{200, z}=10^{13.8}M_{\odot}$, at $z=1.0$: $m_{*, BCG}$ appears to have grown by $0.13\pm0.11$ dex, in tension at $\sim 2.5 \sigma$ significance level with the $0.40$ dex growth rate expected from the semi-analytic model. We show that the buildup of extended intra-cluster light after $z=1.0$ may alleviate this tension in BCG growth rates.

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

We present results from a new set of 30 cosmological simulations of galaxy clusters, including the effects of radiative cooling, star formation, supernova feedback, black hole growth and AGN feedback. We first demonstrate that our AGN model is capable of reproducing the observed cluster pressure profile at redshift, z~0, once the AGN heating temperature of the targeted particles is made to scale with the final virial temperature of the halo. This allows the ejected gas to reach larger radii in higher-mass clusters than would be possible had a fixed heating temperature been used. Such a model also successfully reduces the star formation rate in brightest cluster galaxies and broadly reproduces a number of other observational properties at low redshift, including baryon, gas and star fractions; entropy profiles outside the core; and the X-ray luminosity-mass relation. Our results are consistent with the notion that the excess entropy is generated via selective removal of the densest material through radiative cooling; supernova and AGN feedback largely serve as regulation mechanisms, moving heated gas out of galaxies and away from cluster cores. However, our simulations fail to address a number of serious issues; for example, they are incapable of reproducing the shape and diversity of the observed entropy profiles within the core region. We also show that the stellar and black hole masses are sensitive to numerical resolution, particularly the gravitational softening length; a smaller value leads to more efficient black hole growth at early times and a smaller central galaxy.

We present novel statistical tools to cross-correlate frequency cleaned thermal Sunyaev-Zel'dovich (tSZ) maps and tomographic weak lensing (wl) convergence maps. Moving beyond the lowest order cross-correlation, we introduce a hierarchy of mixed higher-order statistics, the cumulants and cumulant correlators, to analyze non-Gaussianity in real space, as well as corresponding polyspectra in the harmonic domain. Using these moments, we derive analytical expressions for the joint two-point probability distribution function (2PDF) for smoothed tSZ (y_s) and convergence (\kappa_s) maps. The presence of tomographic information allows us to study the evolution of higher order {\em mixed} tSZ-weak lensing statistics with redshift. We express the joint PDFs p_{\kappa y}(\kappa_s,y_s) in terms of individual one-point PDFs (p_{\kappa}(\kappa_s), p_y(y_s)) and the relevant bias functions (b_{\kappa}(\kappa_s), b_y(y_s)). Analytical results for two different regimes are presented that correspond to the small and large angular smoothing scales. Results are also obtained for corresponding {\em hot spots} in the tSZ and convergence maps. In addition to results based on hierarchical techniques and perturbative methods, we present results of calculations based on the lognormal approximation. The analytical expressions derived here are generic and applicable to cross-correlation studies of arbitrary tracers of large scale structure including e.g. that of tSZ and soft X-ray background.

Recent results by the Planck collaboration have shown that cosmological parameters derived from the cosmic microwave background anisotropies and cluster number counts are in tension, with the latter preferring lower values of the matter density parameter, $\Omega_\mathrm{m}$, and power spectrum amplitude, $\sigma_8$. Motivated by this, we investigate the extent to which the tension may be ameliorated once the effect of baryonic depletion on the cluster mass function is taken into account. We use the large-volume Millennium Gas simulations in our study, including one where the gas is pre-heated at high redshift and one where the gas is heated by stars and active galactic nuclei (in the latter, the self-gravity of the baryons and radiative cooling are omitted). In both cases, the cluster baryon fractions are in reasonably good agreement with the data at low redshift, showing significant depletion of baryons with respect to the cosmic mean. As a result, it is found that the cluster abundance in these simulations is around 15 per cent lower than the commonly-adopted fit to dark matter simulations by Tinker et al (2008) for the mass range $10^{14}-10^{14.5}h^{-1} \mathrm{M}_\odot$. Ignoring this effect produces a significant artificial shift in cosmological parameters which can be expressed as $\Delta[\sigma_8(\Omega_\mathrm{m}/0.27)^{0.38}]\simeq -0.03$ at $z=0.17$ (the median redshift of the $\mathit{Planck}$ cluster sample) for the feedback model. While this shift is not sufficient to fully explain the $\mathit{Planck}$ discrepancy, it is clear that such an effect cannot be ignored in future precision measurements of cosmological parameters with clusters. Finally, we outline a simple, model-independent procedure that attempts to correct for the effect of baryonic depletion and show that it works if the baryon-dark matter back-reaction is negligible.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1.

We perform high resolution N-body+SPH simulations of isolated Milky-Way-like galaxies and major mergers between them, to investigate the effect of feedback from both an active galactic nucleus (AGN) and supernovae on the galaxy's evolution. Several AGN methods from the literature are used independently and in conjunction with supernova feedback to isolate the most important factors of these feedback processes. We find that in isolated galaxies, supernovae dominate the suppression of star formation but the star formation rate is unaffected by the presence of an AGN. In mergers the converse is true when models with strong AGN feedback are considered, shutting off star formation before a starburst can occur. AGN and supernovae simulated together suppress star formation only slightly more than if they acted independently. This low-level interaction between the feedback processes is due to AGN feedback maintaining the temperature of a hot halo of gas formed by supernovae. For each of the feedback processes the heating temperature is the dominant parameter rather than the overall energy budget or timing of heating events. Finally, we find that the black hole mass is highly resolution dependent, with more massive black holes found in lower resolution simulations.

We use numerical simulations to investigate how the statistical properties of dark matter (DM) haloes are affected by the baryonic processes associated with galaxy formation. We focus on how these processes influence the spin and shape of a large number of DM haloes covering a wide range of mass scales, from galaxies to clusters at redshifts zero and one, extending to dwarf galaxies at redshift two. The haloes are extracted from the OverWhelmingly Large Simulations, a suite of state-of-the-art high-resolution cosmological simulations run with a range of feedback prescriptions. We find that the median spin parameter in DM-only simulations is independent of mass, redshift and cosmology. At z=0 baryons increase the spin of the DM in the central region (<=0.25 r_200) by up to 30 per cent when feedback is weak or absent. This increase can be attributed to the transfer of angular momentum from baryons to the DM, but is no longer present at z=2. We also present fits to the mass dependence of the DM halo shape at both low and high redshift. At z=0 the sphericity (triaxiality) is negatively (positively) correlated with halo mass and both results are independent of cosmology. Interestingly, these mass-dependent trends are markedly weaker at z=2. While the cooling of baryons acts to make the overall DM halo more spherical, stronger feedback prescriptions tend to reduce the impact of baryons by reducing the central halo mass concentration. More generally, we demonstrate a strongly positive (negative) correlation between halo sphericity (triaxiality) and galaxy formation efficiency, with the latter measured using the central halo baryon fraction. In conclusion, our results suggest that the effects of baryons on the DM halo spin and shape are minor when the effects of cooling are mitigated, as required by realistic models of galaxy formation, although they remain significant for the inner halo.

At high angular frequencies, beyond the damping tail of the primary power spectrum, the dominant contribution to the power spectrum of cosmic microwave background (CMB) temperature fluctuations is the thermal Sunyaev-Zel'dovich (tSZ) effect. We investigate various important statistical properties of the Sunyaev-Zel'dovich maps, using well-motivated models for dark matter clustering to construct statistical descriptions of the tSZ effect to all orders enabling us to determine the entire probability distribution function (PDF). Any generic deterministic biasing scheme can be incorporated in our analysis and the effects of projection, biasing and the underlying density distribution can be analysed separately and transparently in this approach. We introduce the cumulant correlators as tools to analyse tSZ catalogs and relate them to corresponding statistical descriptors of the underlying density distribution. The statistics of hot spots in frequency-cleaned tSZ maps are also developed in a self-consistent way to an arbitrary order, to obtain results complementary to those found using the halo model. We also consider different beam sizes, to check the extent to which the PDF can be extracted from various observational configurations. The formalism is presented with two specific models for underlying matter clustering: (1) the hierarchical ansatz; and (2) the lognormal distribution. We find both models to be in very good agreement with the simulation results, though the lognormal model has an edge over the hierarchical model.

We measure the evolution of the X-ray luminosity-temperature (L_X-T) relation since z~1.5 using a sample of 211 serendipitously detected galaxy clusters with spectroscopic redshifts drawn from the XMM Cluster Survey first data release (XCS-DR1). This is the first study spanning this redshift range using a single, large, homogeneous cluster sample. Using an orthogonal regression technique, we find no evidence for evolution in the slope or intrinsic scatter of the relation since z~1.5, finding both to be consistent with previous measurements at z~0.1. However, the normalisation is seen to evolve negatively with respect to the self-similar expectation: we find E(z)^{-1} L_X = 10^{44.67 +/- 0.09} (T/5)^{3.04 +/- 0.16} (1+z)^{-1.5 +/- 0.5}, which is within 2 sigma of the zero evolution case. We see milder, but still negative, evolution with respect to self-similar when using a bisector regression technique. We compare our results to numerical simulations, where we fit simulated cluster samples using the same methods used on the XCS data. Our data favour models in which the majority of the excess entropy required to explain the slope of the L_X-T relation is injected at high redshift. Simulations in which AGN feedback is implemented using prescriptions from current semi-analytic galaxy formation models predict positive evolution of the normalisation, and differ from our data at more than 5 sigma. This suggests that more efficient feedback at high redshift may be needed in these models.

We explore the dynamical signatures imprinted by baryons on dark matter haloes during the formation process using the OverWhelmingly Large Simulations (OWLS), a set of state-of-the-art high-resolution cosmological hydrodynamical simulations. We present a detailed study of the effects of the implemented feedback prescriptions on the orbits of dark matter particles, stellar particles and subhaloes, analysing runs with no feedback, with stellar feedback and with feedback from supermassive black holes. We focus on the central regions (0.25 r_{200}) of haloes with virial masses ~ 6 x 10^{13} (~ 7 x 10^{11}) Msun/h at z = 0(2). We also investigate how the orbital content (relative fractions of the different orbital types) of these haloes depends on several key parameters such as their mass, redshift and dynamical state. The results of spectral analyses of the orbital content of these simulations are compared, and the change in fraction of box, tube and irregular orbits is quantified. Box orbits are found to dominate the orbital structure of dark matter haloes in cosmological simulations. There is a strong anticorrelation between the fraction of box orbits and the central baryon fraction. While radiative cooling acts to reduce the fraction of box orbits, strong feedback implementations result in a similar orbital distribution to that of the dark matter only case. The orbital content described by the stellar particles is found to be remarkably similar to that drawn from the orbits of dark matter particles, suggesting that either they have forgotten their dynamical history, or that subhaloes bringing in stars are not biased significantly with respect to the main distribution. The orbital content of the subhaloes is in broad agreement with that seen in the outer regions of the particle distributions.

Using a sample of 123 X-ray clusters and groups drawn from the XMM-Cluster Survey first data release, we investigate the interplay between the brightest cluster galaxy (BCG), its black hole, and the intra-cluster/group medium (ICM). It appears that for groups and clusters with a BCG likely to host significant AGN feedback, gas cooling dominates in those with Tx > 2 keV while AGN feedback dominates below. This may be understood through the sub-unity exponent found in the scaling relation we derive between the BCG mass and cluster mass over the halo mass range 10^13 < M500 < 10^15Msol and the lack of correlation between radio luminosity and cluster mass, such that BCG AGN in groups can have relatively more energetic influence on the ICM. The Lx - Tx relation for systems with the most massive BCGs, or those with BCGs co-located with the peak of the ICM emission, is steeper than that for those with the least massive and most offset, which instead follows self-similarity. This is evidence that a combination of central gas cooling and powerful, well fuelled AGN causes the departure of the ICM from pure gravitational heating, with the steepened relation crossing self-similarity at Tx = 2 keV. Importantly, regardless of their black hole mass, BCGs are more likely to host radio-loud AGN if they are in a massive cluster (Tx > 2 keV) and again co-located with an effective fuel supply of dense, cooling gas. This demonstrates that the most massive black holes appear to know more about their host cluster than they do about their host galaxy. The results lead us to propose a physically motivated, empirical definition of 'cluster' and 'group', delineated at 2 keV.

We present a list of 15 clusters of galaxies, serendipitously detected by the XMM Cluster Survey (XCS), that have a high probability of detection by the Planck satellite. Three of them already appear in the Planck Early Sunyaev-Zel'dovich (ESZ) catalogue. The estimation of the Planck detection probability assumes the flat Lambda cold dark matter (LambdaCDM) cosmology most compatible with 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) data. It takes into account the XCS selection function and Planck sensitivity, as well as the covariance of the cluster X-ray luminosity, temperature, and integrated comptonization parameter, as a function of cluster mass and redshift, determined by the Millennium Gas Simulations. We also characterize the properties of the galaxy clusters in the final data release of the XCS that we expect Planck will have detected by the end of its extended mission. Finally, we briefly discuss possible joint applications of the XCS and Planck data.

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.

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations.

The XMM Cluster Survey (XCS) is a serendipitous search for galaxy clusters using all publicly available data in the XMM-Newton Science Archive. Its main aims are to measure cosmological parameters and trace the evolution of X-ray scaling relations. In this paper we present the first data release from the XMM Cluster Survey (XCS-DR1). This consists of 503 optically confirmed, serendipitously detected, X-ray clusters. Of these clusters, 255 are new to the literature and 356 are new X-ray discoveries. We present 464 clusters with a redshift estimate (0.06 < z < 1.46), including 261 clusters with spectroscopic redshifts. In addition, we have measured X-ray temperatures (Tx) for 402 clusters (0.4 < Tx < 14.7 keV). We highlight seven interesting subsamples of XCS-DR1 clusters: (i) 10 clusters at high redshift (z > 1.0, including a new spectroscopically-confirmed cluster at z = 1.01); (ii) 67 clusters with high Tx (> 5 keV); (iii) 131 clusters/groups with low Tx (< 2 keV); (iv) 27 clusters with measured Tx values in the SDSS `Stripe 82' co-add region; (v) 78 clusters with measured Tx values in the Dark Energy Survey region; (vi) 40 clusters detected with sufficient counts to permit mass measurements (under the assumption of hydrostatic equilibrium); (vii) 105 clusters that can be used for applications such as the derivation of cosmological parameters and the measurement of cluster scaling relations. The X-ray analysis methodology used to construct and analyse the XCS-DR1 cluster sample has been presented in a companion paper, Lloyd-Davies et al. (2010).

The XMM Cluster Survey (XCS) is a serendipitous search for galaxy clusters using all publicly available data in the XMM-Newton Science Archive. Its main aims are to measure cosmological parameters and trace the evolution of X-ray scaling relations. In this paper we describe the data processing methodology applied to the 5,776 XMM observations used to construct the current XCS source catalogue. A total of 3,675 > 4-sigma cluster candidates with > 50 background-subtracted X-ray counts are extracted from a total non-overlapping area suitable for cluster searching of 410 deg^2. Of these, 993 candidates are detected with > 300 background-subtracted X-ray photon counts, and we demonstrate that robust temperature measurements can be obtained down to this count limit. We describe in detail the automated pipelines used to perform the spectral and surface brightness fitting for these candidates, as well as to estimate redshifts from the X-ray data alone. A total of 587 (122) X-ray temperatures to a typical accuracy of < 40 (< 10) per cent have been measured to date. We also present the methodology adopted for determining the selection function of the survey, and show that the extended source detection algorithm is robust to a range of cluster morphologies by inserting mock clusters derived from hydrodynamical simulations into real XMM images. These tests show that the simple isothermal beta-profiles is sufficient to capture the essential details of the cluster population detected in the archival XMM observations. The redshift follow-up of the XCS cluster sample is presented in a companion paper, together with a first data release of 503 optically-confirmed clusters.

We use a combination of three large N-body simulations to investigate the dependence of dark matter halo concentrations on halo mass and redshift in the WMAP year 5 cosmology. The median relation between concentration and mass is adequately described by a power-law for halo masses in the range 10^11 - 10^15 Msol/h and redshifts z < 2, regardless of whether the halo density profiles are fit using NFW or Einasto profiles. Compared with recent analyses of the Millennium Simulation, which uses a value of sigma_8 that is higher than allowed by WMAP5, z = 0 halo concentrations are reduced by factors ranging from 23 per cent at 10^11 Msol/h to 16 per cent at 10^14 Msol/h. The predicted concentrations are much lower than inferred from X-ray observations of groups and clusters.

We present a parameterized model of the intra-cluster medium that is suitable for jointly analysing pointed observations of the Sunyaev-Zel'dovich (SZ) effect and X-ray emission in galaxy clusters. The model is based on assumptions of hydrostatic equilibrium, the Navarro, Frenk and White (NFW) model for the dark matter, and a softened power law profile for the gas entropy. We test this entropy-based model against high and low signal-to-noise mock observations of a relaxed and recently-merged cluster from N-body/hydrodynamic simulations, using Bayesian hyper-parameters to optimise the relative statistical weighting of the mock SZ and X-ray data. We find that it accurately reproduces both the global values of the cluster temperature, total mass and gas mass fraction (fgas), as well as the radial dependencies of these quantities outside of the core (r > kpc). For reference we also provide a comparison with results from the single isothermal beta model. We confirm previous results that the single isothermal beta model can result in significant biases in derived cluster properties.

We present deep J and Ks band photometry of 20 high redshift galaxy clusters between z=0.8-1.5, 19 of which are observed with the MOIRCS instrument on the Subaru Telescope. By using near-infrared light as a proxy for stellar mass we find the surprising result that the average stellar mass of Brightest Cluster Galaxies (BCGs) has remained constant at ~9e11MSol since z~1.5. We investigate the effect on this result of differing star formation histories generated by three well known and independent stellar population codes and find it to be robust for reasonable, physically motivated choices of age and metallicity. By performing Monte Carlo simulations we find that the result is unaffected by any correlation between BCG mass and cluster mass in either the observed or model clusters. The large stellar masses imply that the assemblage of these galaxies took place at the same time as the initial burst of star formation. This result leads us to conclude that dry merging has had little effect on the average stellar mass of BCGs over the last 9-10 Gyr in stark contrast to the predictions of semi-analytic models, based on the hierarchical merging of dark matter haloes, which predict a more protracted mass build up over a Hubble time. We discuss however that there is potential for reconciliation between observation and theory if there is a significant growth of material in the intracluster light over the same period.

We use Chandra X-ray and Spitzer infrared observations to explore the AGN and starburst populations of XMMXCS J2215.9-1738 at z=1.46, one of the most distant spectroscopically confirmed galaxy clusters known. The high resolution X-ray imaging reveals that the cluster emission is contaminated by point sources that were not resolved in XMM observations of the system, and have the effect of hardening the spectrum, leading to the previously reported temperature for this system being overestimated. From a joint spectroscopic analysis of the Chandra and XMM data, the cluster is found to have temperature T=4.1_-0.9^+0.6 keV and luminosity L_X=(2.92_-0.35^+0.24)x10^44 erg/s extrapolated to a radius of 2 Mpc. As a result of this revised analysis, the cluster is found to lie on the sigma_v-T relation, but the cluster remains less luminous than would be expected from self-similar evolution of the local L_X-T relation. Two of the newly discovered X-ray AGN are cluster members, while a third object, which is also a prominent 24 micron source, is found to have properties consistent with it being a high redshift, highly obscured object in the background. We find a total of eight >5 sigma 24 micron sources associated with cluster members (four spectroscopically confirmed, and four selected using photometric redshifts), and one additional 24 micron source with two possible optical/near-IR counterparts that may be associated with the cluster. Examining the IRAC colors of these sources, we find one object is likely to be an AGN. Assuming that the other 24 micron sources are powered by star formation, their infrared luminosities imply star formation rates ~100 M_sun/yr. We find that three of these sources are located at projected distances of <250 kpc from the cluster center, suggesting that a large amount of star formation may be taking place in the cluster core, in contrast to clusters at low redshift.

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

In this paper, we undertake a study to determine what insight can be reliably gleaned from the comparison of the X-ray and the weak lensing mass maps of galaxy clusters. We do this by investigating the 2D substructure within three high-resolution cosmological simulations of galaxy clusters. Our main results focus on non-radiative gas dynamics, but we also consider the effects of radiative cooling at high redshift. For our analysis, we use a novel approach, based on unsharp-masking, to identify substructures in 2D surface mass density and X-ray surface brightness maps. At full resolution (~ 15 h^-1 kpc), this technique is capable of identifying almost all self-bound dark matter subhaloes with M>10^12 h^-1 M_sun. We also report a correlation between the mass of a subhalo and the area of its corresponding 2D detection; such a correlation, once calibrated, could provide a useful estimator for substructure mass. Comparing our 2D mass and X-ray substructures, we find a surprising number of cases where the matching fails: around one third of galaxy-sized substructures have no X-ray counterpart. Some interesting cases are also found at larger masses, in particular the cores of merging clusters where the situation can be complex. Finally, we degrade our mass maps to what is currently achievable with weak-lensing observations (~100 h^-1 kpc at z=0.2). While the completeness mass limit increases by around an order of magnitude, a mass-area correlation remains. Our paper clearly demonstrates that the next generation of lensing surveys should start to reveal a wealth of information on cluster substructure. (Abridged)

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

Using large numbers of simulations of the microwave sky, incorporating the Cosmic Microwave Background (CMB) and the Sunyaev-Zel'dovich (SZ) effect due to clusters, we investigate the statistics of the power spectrum at microwave frequencies between spherical multipoles of 1000 and 10000. From these virtual sky maps, we find that the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than would be expected from purely Gaussian realizations, and has a distribution that is significantly skewed towards higher values, especially when small map sizes are used. The standard deviation is also increased by around 10 percent compared to the trispectrum calculation due to the clustering of galaxy clusters. We also consider the effects of including residual point sources and uncertainties in the gas physics. This has implications for the excess power measured in the CMB power spectrum by the Cosmic Background Imager and BIMA experiments. Our results indicate that the observed excess could be explained using a lower value of $\sigma_8$ than previously suggested, however the effect is not enough to match $\sigma_8=0.825$. The uncertainties in the gas physics could also play a substantial role. We have made our maps of the SZ effect available online.

We forecast the constraints on the values of sigma_8, Omega_m, and cluster scaling relation parameters which we expect to obtain from the XMM Cluster Survey (XCS). We assume a flat Lambda-CDM Universe and perform a Monte Carlo Markov Chain analysis of the evolution of the number density of galaxy clusters that takes into account a detailed simulated selection function. Comparing our current observed number of clusters shows good agreement with predictions. We determine the expected degradation of the constraints as a result of self-calibrating the luminosity-temperature relation (with scatter), including temperature measurement errors, and relying on photometric methods for the estimation of galaxy cluster redshifts. We examine the effects of systematic errors in scaling relation and measurement error assumptions. Using only (T,z) self-calibration, we expect to measure Omega_m to +-0.03 (and Omega_Lambda to the same accuracy assuming flatness), and sigma_8 to +-0.05, also constraining the normalization and slope of the luminosity-temperature relation to +-6 and +-13 per cent (at 1sigma) respectively in the process. Self-calibration fails to jointly constrain the scatter and redshift evolution of the luminosity-temperature relation significantly. Additional archival and/or follow-up data will improve on this. We do not expect measurement errors or imperfect knowledge of their distribution to degrade constraints significantly. Scaling-relation systematics can easily lead to cosmological constraints 2sigma or more away from the fiducial model. Our treatment is the first exact treatment to this level of detail, and introduces a new `smoothed ML' estimate of expected constraints.

The current consensus is that galaxies begin as small density fluctuations in the early Universe and grow by in situ star formation and hierarchical merging. Stars begin to form relatively quickly in sub-galactic sized building blocks called haloes which are subsequently assembled into galaxies. However, exactly when this assembly takes place is a matter of some debate. Here we report that the stellar masses of brightest cluster galaxies, which are the most luminous objects emitting stellar light, some 9 billion years ago are not significantly different from their stellar masses today. Brightest cluster galaxies are almost fully assembled 4-5 Gyrs after the Big Bang, having grown to more than 90% of their final stellar mass by this time. Our data conflict with the most recent galaxy formation models based on the largest simulations of dark matter halo development. These models predict protracted formation of brightest cluster galaxies over a Hubble time, with only 22% of the stellar mass assembled at the epoch probed by our sample. Our findings suggest a new picture in which brightest cluster galaxies experience an early period of rapid growth rather than prolonged hierarchical assembly.

We present a study of the morphological fractions and color-magnitude relation in the most distant X-ray selected galaxy cluster currently known, XMMXCS J2215.9-1738 at z=1.46, using a combination of optical imaging data obtained with the Hubble Space Telescope Advanced Camera for Surveys, and infrared data from the Multi-Object Infrared Camera and Spectrograph, mounted on the 8.2m Subaru telescope. We find that the morphological mix of the cluster galaxy population is similar to clusters at z~1: approximately ~62% of the galaxies identified as likely cluster members are ellipticals or S0s; and ~38% are spirals or irregulars. We measure the color-magnitude relations for the early type galaxies, finding that the slope in the z_850-J relation is consistent with that measured in the Coma cluster, some ~9 Gyr earlier, although the uncertainty is large. In contrast, the measured intrinsic scatter about the color-magnitude relation is more than three times the value measured in Coma, after conversion to rest frame U-V. From comparison with stellar population synthesis models, the intrinsic scatter measurements imply mean luminosity weighted ages for the early type galaxies in J2215.9-1738 of ~3 Gyr, corresponding to the major epoch of star formation coming to an end at z_f = 3-5. We find that the cluster exhibits evidence of the `downsizing' phenomenon: the fraction of faint cluster members on the red sequence expressed using the Dwarf-to-Giant Ratio (DGR) is 0.32+/-0.18. This is consistent with extrapolation of the redshift evolution of the DGR seen in cluster samples at z < 1. In contrast to observations of some other z > 1 clusters, we find a lack of very bright galaxies within the cluster.

We demonstrate that the large scatter in the ultra-violet (UV) colours of intermediate-mass early-type galaxies in the local Universe and the inferred low-level recent star formation in these objects can be reproduced by minor mergers in the standard LCDM cosmology. Numerical simulations of mergers with mass ratios less than or equal to 1:4, with reasonable assumptions for the ages, metallicities and dust properties of the merger progenitors, produce good agreement to the observed UV colours of the early-type population, if the infalling satellites are assumed to have (cold) gas fractions greater than 20%. Early-types that satisfy (NUV-r) < 3.8 are likely to have experienced mergers with mass ratios between 1:4 and 1:6 within the last ~1.5 Gyrs, while those that satisfy 3.8<(NUV-r)<5.5 are consistent with either recent mergers with mass ratios < 1:6 or mergers with higher mass ratios that occurred more than ~1.5 Gyrs in the past. We demonstrate that the early-type colour-magnitude relations and colour distributions in both the UV and optical spectral ranges are consistent with the expected frequency of minor merging activity in the standard LCDM cosmology at low redshift. We present a strong plausibility argument for minor mergers to be the principal mechanism behind the large UV scatter and associated low-level recent star formation observed in early-type galaxies in the nearby Universe.

Substructures, expected in cold dark matter haloes, have been proposed to explain the anomalous flux ratios in gravitational lenses. About 25% of lenses in the Cosmic Lens All-Sky Survey (CLASS) appear to have luminous satellites within ~ 5 kpc/h of the main lensing galaxies, which are usually at redshift z ~ 0.2-1. In this work we use the Millennium Simulation combined with galaxy catalogues from semi-analytical techniques to study the predicted frequency of such satellites in simulated haloes. The fraction of haloes that host bright satellites within the (projected) central regions is similar for red and blue hosts and is found to increase as a function of host halo mass and redshift. Specifically, at z = 1, about 11% of galaxy-sized haloes (with masses between 10^{12} M_sun/h and 10^{13} M_sun/h) host bright satellite galaxies within a projected radius of 5 kpc/h. This fraction increases to about 17% (25%) if we consider bright (all) satellites of only group-sized haloes (with masses between 10^{13} M_sun/h and 10^{14} M_sun/h). These results are roughly consistent with the fraction (~ 25%) of CLASS lensing galaxies observed to host luminous satellites. At z = 0, only ~ 3% of galaxy-sized haloes host bright satellite galaxies. The fraction rises to ~ 6%, (10%) if we consider bright (all) satellites of only group-sized haloes at z = 0. However, most of the satellites found in the inner regions are `orphan' galaxies where the dark matter haloes have been completely stripped. Thus the agreement crucially depends on the true survival rate of these `orphan' galaxies. We also discuss the effects of numerical resolution and cosmologies on our results.

The internal dynamics of a dark matter structure may have the remarkable property that the local temperature in the structure depends on direction. This is parametrized by the velocity anisotropy beta which must be zero for relaxed collisional structures, but has been shown to be non-zero in numerical simulations of dark matter structures. Here we present a method to infer the radial profile of the velocity anisotropy of the dark matter halo in a galaxy cluster from X-ray observables of the intracluster gas. This non-parametric method is based on a universal relation between the dark matter temperature and the gas temperature which is confirmed through numerical simulations. We apply this method to observational data and we find that beta is significantly different from zero at intermediate radii. Thus we find a strong indication that dark matter is effectively collisionless on the dynamical time-scale of clusters, which implies an upper limit on the self-interaction cross-section per unit mass sigma/m < 1 cm2/g. Our results may provide an independent way to determine the stellar mass density in the central regions of a relaxed cluster, as well as a test of whether a cluster is in fact relaxed.

The relativistic Sunyaev-Zel'dovich (SZ) effect offers a method, independent of X-ray, for measuring the temperature of the intracluster medium (ICM) in the hottest systems. Here, using N-body/hydrodynamic simulations of three galaxy clusters, we compare the two quantities for a non-radiative ICM, and for one that is subject both to radiative cooling and strong energy feedback from galaxies. Our study has yielded two interesting results. Firstly, in all cases, the SZ temperature is hotter than the X-ray temperature and is within ten per cent of the virial temperature of the cluster. Secondly, the mean SZ temperature is less affected by cooling and feedback than the X-ray temperature. Both these results can be explained by the SZ temperature being less sensitive to the distribution of cool gas associated with cluster substructure. A comparison of the SZ and X-ray temperatures (measured for a sample of hot clusters) would therefore yield interesting constraints on the thermodynamic structure of the intracluster gas.

(Context) Central active galactic nuclei (AGN) are supposed to play a key role in the evolution of their host galaxies. In particular, the dynamical and physical properties of the gas core must be affected by the injected energy. (Aims) Our aim is to study the effects of an AGN on the dark matter profile and on the central stellar light distribution in massive early type galaxies. (Methods) By performing self-consistent N-body simulations, we assume in our analysis that periodic bipolar outbursts from a central AGN can induce harmonic oscillatory motions on both sides of the gas core. (Results) Using realistic AGN properties, we find that the motions of the gas core, driven by such feedback processes, can flatten the dark matter and/or stellar profiles after 4-5 Gyr. These results are consistent with recent observational studies that suggest that most giant elliptical galaxies have cores or are "missing light" in their inner part. Since stars behave as a "collisionless" fluid similar to dark matter, the density profile both of stars and dark matter should be affected in a similar way, leading to an effective reduction in the central brightness.

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.

We present new spectroscopic observations of the most distant X-ray selected galaxy cluster currently known, XMMXCS J2215.9-1738 at z=1.457, obtained with the DEIMOS instrument at the W. M. Keck Observatory, and the FORS2 instrument on the ESO Very Large Telescope. Within the cluster virial radius, as estimated from the cluster X-ray properties, we increase the number of known spectroscopic cluster members to 17 objects, and calculate the line of sight velocity dispersion of the cluster to be 580+/-140 km/s. We find mild evidence that the velocity distribution of galaxies within the virial radius deviates from a single Gaussian. We show that the properties of J2215.9-1738 are inconsistent with self-similar evolution of local X-ray scaling relations, finding that the cluster is underluminous given its X-ray temperature, and that the intracluster medium contains ~2-3 times the kinetic energy per unit mass of the cluster galaxies. These results can perhaps be explained if the cluster is observed in the aftermath of an off-axis merger. Alternatively, heating of the intracluster medium through supernovae and/or Active Galactic Nuclei activity, as is required to explain the observed slope of the local X-ray luminosity-temperature relation, may be responsible.

Various differences in galaxy cluster properties derived from X-ray and weak lensing observations have been highlighted in the literature. One such difference is the observation of mass concentrations in lensing maps which have no X-ray counterparts (e.g. Jee, White, Ford et al. 2005). We investigate this issue by identifying substructures in maps of projected total mass (analogous to weak lensing mass reconstructions) and maps of projected X-ray surface brightness for three simulated clusters. We then compare the 2D mass substructures with both 3D subhalo data and the 2D X-ray substructures. Here we present preliminary results from the first comparison, where we have assessed the impact of projecting the data on subhalo identification.

We present results from a study of the X-ray cluster population that forms within the CLEF cosmological hydrodynamics simulation, a large N-body/SPH simulation of the Lambda CDM cosmology with radiative cooling, star formation and feedback. The scaled projected temperature and entropy profiles at z=0 are in good agreement with recent high-quality observations of cool core clusters, suggesting that the simulation grossly follows the processes that structure the intracluster medium (ICM) in these objects. Cool cores are a ubiquitous phenomenon in the simulation at low and high redshift, regardless of a cluster's dynamical state. This is at odds with the observations and so suggests there is still a heating mechanism missing from the simulation. Using a simple, observable measure of the concentration of the ICM, which correlates with the apparent mass deposition rate in the cluster core, we find a large dispersion within regular clusters at low redshift, but this diminishes at higher redshift, where strong "cooling-flow" systems are absent in our simulation. Consequently, our results predict that the normalisation and scatter of the luminosity-temperature relation should decrease with redshift; if such behaviour turns out to be a correct representation of X-ray cluster evolution, it will have significant consequences for the number of clusters found at high redshift in X-ray flux-limited surveys.

We report the discovery of XMMXCS J2215.9-1738, a massive galaxy cluster at z =1.45, which was found in the XMM Cluster Survey. The cluster candidate was initially identified as an extended X-ray source in archival XMM data. Optical spectroscopy shows that 6 galaxies within a 60 arcsec diameter region lie at z = 1.45 +/- 0.01. Model fits to the X-ray spectra of the extended emission yield kT = 7.4 (+2.7,-1.8) keV (90 % confidence); if there is an undetected central X-ray point source then kT = 6.5 (+2.6,-1.8) keV. The bolometric X-ray luminosity is Lx = 4.4 (+0.8,-0.6) x 10^44 ergs/s over a 2 Mpc radial region. The measured Tx, which is the highest known for a cluster at z > 1, suggests that this cluster is relatively massive for such a high redshift. The redshift of XMMXCS J2215.9-1738 is the highest currently known for a spectroscopically-confirmed cluster of galaxies.

We investigate the redshift dependence of X-ray cluster scaling relations drawn from three hydrodynamic simulations of the LCDM cosmology: a Radiative model that incorporates radiative cooling of the gas, a Preheating model that additionally heats the gas uniformly at high redshift, and a Feedback model that self-consistently heats cold gas in proportion to its local star-formation rate. While all three models are capable of reproducing the observed local Lx-Tx relation, they predict substantially different results at high redshift (to z=1.5), with the Radiative, Preheating and Feedback models predicting strongly positive, mildly positive and mildly negative evolution, respectively. The physical explanation for these differences lies in the structure of the intracluster medium. All three models predict significant temperature fluctuations at any given radius due to the presence of cool subclumps and, in the case of the Feedback simulation, reheated gas. The mean gas temperature lies above the dynamical temperature of the halo for all models at z=0, but differs between models at higher redshift with the Radiative model having the lowest mean gas temperature at z=1.5. We have not attempted to model the scaling relations in a manner that mimics the observational selection effects, nor has a consistent observational picture yet emerged. Nevertheless, evolution of the scaling relations promises to be a powerful probe of the physics of entropy generation in clusters. First indications are that early, widespread heating is favored over an extended period of heating that is associated with galaxy formation.

The radial entropy profile of the hot gas in clusters of galaxies tends to follow a power law in radius outside of the cluster core. Here we present a simple formula giving both the normalization and slope for the power-law entropy profiles of clusters that form in the absence of non-gravitational processes such as radiative cooling and subsequent feedback. It is based on seventy-one clusters drawn from four separate cosmological simulations, two using smoothed-particle hydrodynamics (SPH) and two using adaptive-mesh refinement (AMR), and can be used as a baseline for assessing the impact of non-gravitational processes on the intracluster medium outside of cluster cores. All the simulations produce clusters with self-similar structure in which the normalization of the entropy profile scales linearly with cluster temperature, and these profiles are in excellent agreement outside of 0.2 r_200. Because the observed entropy profiles of clusters do not scale linearly with temperature, our models confirm that non-gravitational processes are necessary to break the self-similarity seen in the simulations. However, the core entropy levels found by the two codes used here significantly differ, with the AMR code producing nearly twice as much entropy at the centre of a cluster.

We use hydrodynamical simulations of galaxy clusters and groups to study the effect of pre-heating on the entropy structure of the ICM. Our simulations account for non-gravitational heating of the gas either by imposing a minimum entropy floor at redshift z=3, or by considering feedback by galactic winds powered by supernova (SN) energy. In the adiabatic simulations we find that the entropy is increased out to the external regions as a consequence of the transition from clumpy to smooth accretion induced by extra heating. This result is in line with the predictions of the semi-analytical model by Voit et al. However, the introduction of radiative cooling substantially reduces this entropy amplification effect. While galactic winds of increasing strength are effective in regulating star formation, they have a negligible effect on the entropy profile of cluster-sized halos. Only in models where the action of the winds is complemented with diffuse heating corresponding to a pre-collapse entropy do we find a sizable entropy amplification out to the virial radius of the groups. Observational evidence for entropy amplification in the outskirts of galaxy clusters and groups therefore favours a scenario for feedback that distributes heating energy in a more diffuse way than predicted by the model for galactic winds from SN explosions explored here.

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.

Using the statistics of pixel optical depths, we compare HI, CIV and CIII absorption in a set of six high quality z ~ 3-4 quasar absorption spectra to that in spectra drawn from two different state-of-the-art cosmological simulations that include galactic outflows. We find that the simulations predict far too little CIV absorption unless the UVB is extremely soft, and always predict far too small CIII/CIV ratios. We note, however, that much of the enriched gas is in a phase (T ~ 10^5-10^7 K, overdensity ~ 0.1-10, Z >~0.1 Z_sol) that should cool by metal line emission -- which was not included in our simulations. When the effect of cooling is modeled, the predicted CIV absorption increases substantially, but the CIII/CIV ratios are still far too small because the density of the enriched gas is too low. Finally, we find that the predicted metal distribution is much too inhomogeneous to reproduce the observed probability distribution of CIV absorption. These findings suggest that strong z <~ 6 winds cannot fully explain the observed enrichment, and that an additional (perhaps higher-z) contribution is required.

Preliminary results are presented from the CLEF hydrodynamics simulation, a large (N=2(428)^3 particles within a 200 Mpc/h comoving box) simulation of the LCDM cosmology that includes both radiative cooling and a simple model for galactic feedback. Specifically, we focus on the X-ray properties of the simulated clusters at z=0 and demonstrate a reasonable level of agreement between simulated and observed cluster scaling relations.

The CLEF-SSH simulation project is an international collaboration, involving the IAS, LATT and the University of Sussex, to produce large hydrodynamic simulations of large-scale structure that implement realistic models of radiative gas cooling and energy feedback. The objective is to use these simulations to study the physics of galaxy clusters and to construct maps of large angular size of the Sunyaev-Zel'dovich (SZ) effect, for the preparation of future experiments. Here we present results from a first run, the CLEF hydrodynamics simulation, which features 2(428)^3 particles of gas and dark matter inside a comoving box with 200 Mpc/h on a side.

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

[ABRIDGED] The XMM Cluster Survey (XCS) is predicted to detect thousands of clusters observed serendipitously in XMM-Newton pointings. We investigate automating optical follow-up of cluster candidates using the SDSS public archive, concentrating on 42 XMM observations that overlap the SDSS First Data Release. Restricting to the inner 11 arcminutes of the XMM field-of-view gives 637 unique X-ray sources across a 3.09 deg^2 region with SDSS coverage. The log N-log S relation indicates survey completeness to a flux limit of around 1x10^{-14} erg s^-1 cm^-2. We have cross-correlated XMM point sources with SDSS quasars, finding 103 unique matches from which we determine a matching radius of 3.8 arcseconds. Using this we make immediate identifications of roughly half of all XMM point sources as quasars (159) or stars (55). We have estimated the typical error on SDSS-determined photometric cluster redshifts to be 5% for relaxed systems, and 11% for disturbed systems, by comparing them to spectroscopic redshifts for eight previously-known clusters. The error for disturbed systems may be problematic for cluster surveys relying on photometric redshifts alone. We use the False Discovery Rate to select 41 XMM sources (25 point-like, 16 extended) statistically associated with SDSS galaxy overdensities. Of the extended sources, 5 are new cluster candidates and 11 are previously-known clusters (0.044<z<0.782). There are 83 extended X-ray sources not associated with SDSS galaxy overdensities, which are strong candidates for new high-redshift clusters. We highlight two previously-known z>1 clusters rediscovered by our wavelet-based detection software. These results show that the SDSS can provide useful automated follow-up to X-ray cluster surveys, but cannot provide all the optical follow-up necessary for X-ray cluster surveys.

We investigate the growth over time of 20 massive (> 3 keV) clusters in a hydrodynamical simulation of the LambdaCDM cosmology with radiative cooling. The clusters show a variety of formation histories: some accrete most of their mass in major mergers; others more gradually. During major mergers the long-term (temporally-smoothed) luminosity increases such that the cluster moves approximately along the Lx-Tx relation; between times it slowly decreases, tracking the drift of the Lx-Tx relation. We identify several different kinds of short-term luminosity and temperature fluctuations associated with major mergers including double-peaked mergers in which the global intracluster medium merges first (Lx and Tx increase together) and then the cluster cores merge (Lx increases and Tx decreases). At both luminosity peaks, clusters tend to appear spherical and relaxed, which may lead to biases in high-redshift, flux-limited samples. There is no simple relationship between scatter in the Lx-Tx relation and either recent or overall merger activity or cluster formation time. The scatter in the Lx-M and Tx-M relations is reduced if properties are measured within R_500 rather than R_vir.

The Sunyaev-Zel'dovich effect is a powerful new tool for finding and studying clusters at high redshift, particularly in combination with their X-ray properties. In this paper we quantify the expected scaling relations between these properties using numerical simulations with various models for heating and cooling of the cluster gas. For a Non-radiative model, we find scaling relations in good agreement with self-similar predictions: $Y\propto T_X^{5/2}$ and $Y \propto L_X^{5/4}$. Our main results focus on predictions from Cooling and Preheating simulations, shown by Muanwong et al. (2002) to provide a good match to the X-ray scaling relations at z=0. For these runs we find slopes of approximately $Y \proptoT_X^3$ and $Y \propto L_X$, steeper and flatter than the self-similar scalings respectively. We also study the redshift evolution of the scaling relations and find the slopes show no evidence of evolution out to redshifts well beyond one, while the normalizations of relations between the SZ signal and X-ray properties do show evolution relative to that expected from self-similarity, particularly at z<1.

The entropy of the intracluster medium at large radii has been shown recently to deviate from the self-similar scaling with temperature. Using N-body/hydrodynamic simulations of the LCDM cosmology, we demonstrate that this deviation is evidence that feedback processes are important in generating excess entropy in clusters. While radiative cooling increases the entropy of intracluster gas, resulting in a good match to the data in the centres of clusters, it produces an entropy-temperature relation closer to the self-similar scaling at larger radii. A model that includes feedback from galaxies, however, not only stabilises the cooling rate in the simulation, but is capable of reproducing the observed scaling behaviour both in cluster cores and at large radii. Feedback modifies the entropy distribution in clusters due to its increasing ability at expelling gas from haloes with decreasing mass. The strength of the feedback required, as suggested from our simulations, is consistent with supernova energetics, providing a large fraction of the energy reaches low-density regions and is originally contained within a small mass of gas.

The amplitude of density perturbations, for the currently-favoured LambdaCDM cosmology, is constrained using the observed properties of galaxy clusters. The catalogue used is that of Ikebe et al. (2002). The cluster temperature to mass relation is obtained via N-body/hydrodynamical simulations including radiative cooling and preheating of cluster gas, which we have previously shown to reproduce well the observed temperature--mass relation in the innermost parts of clusters (Thomas et al. 2002). We generate and compare mock catalogues via a Monte Carlo method, which allows us to constrain the relation between X-ray temperature and luminosity, including its scatter, simultaneously with cosmological parameters. We find a luminosity-temperature relation in good agreement with the results of Ikebe et al. (2002), while for the matter power spectrum normalization, we find $\sigma_8 = 0.78_{-0.06}^{+0.30}$ at 95 per cent confidence for $\Omega_0 = 0.35$. Scaling to WMAP's central value of $\Omega_0 = 0.27$ would give a best-fit value of $\sigma_8 \simeq 0.9$.

Using hydrodynamical simulations of the LCDM cosmology, that include both radiative cooling and a phenomenological model for star formation and supernovae feedback, we investigate the impact of galaxy formation on the X-ray properties of groups at redshift zero. Motivated by the observed ``break'' in the L_x-T_x relation at kT~1-2keV, our feedback model is based on the assumption that supernovae imprint a temperature scale on the hot gas, with the star formation rate and corresponding reheated gas mass then depending only on the available energy budget. We demonstrate that a strong feedback model with a heating temperature comparable to this break (kT_=2keV), and an energy budget twice that available from supernovae (epsilon=2), raises the core entropy of groups sufficiently to produce an adequate match to their observed X-ray properties. A lower value of epsilon increases the star formation rate without significantly affecting the X-ray properties of groups, and a model with epsilon~0.1 reproduces the observed fraction of baryons in stars. However, a heating temperature which is lower than the virial temperatures of the groups leads to an excess of cooling gas that boosts their X-ray luminosities, due to the failure of the reheated material to escape from the gravitational potential. A limited study of numerical resolution effects reveals that the temperature of poorly-resolved objects is underestimated, therefore (in our case) a fully-resolved group population would lead to a steeper L_x-T_x relation, bringing our results into even better agreement with the observations.

Temperature maps are presented of the 9 largest clusters in the mock catalogues of Muanwong et al. for both the Preheating and Radiative models. The maps show that clusters are not smooth, featureless systems, but contain a variety of substructure which should be observable. The surface brightness contours are generally elliptical and features that are seen include cold clumps, hot spiral features, and cold fronts. Profiles of emission-weighted temperature, surface brightness and emission-weighted pressure across the surface brightness discontinuities seen in one of the bimodal clusters are consistent with the cold front in Abell 2142 observed by Markevitch et al.

We have performed hydrodynamical simulations to investigate the effects of galactic winds on the high-redshift (z=3) universe. Strong winds suppress the formation of low-mass galaxies significantly, and the metals carried by them produce C IV absorption lines with properties in reasonable agreement with observations. The winds have little effect on the statistics of the H I absorption lines, because the hot gas bubbles blown by the winds fill only a small fraction of the volume and because they tend to escape into the voids, thereby leaving the filaments that produce these lines intact.

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.

It has become increasingly apparent that traditional hydrodynamical simulations of galaxy clusters are unable to reproduce the observed properties of galaxy clusters, in particular overpredicting the mass corresponding to a given cluster temperature. Such overestimation may lead to systematic errors in results using galaxy clusters as cosmological probes, such as constraints on the density perturbation normalization sigma_8. In this paper we demonstrate that inclusion of additional gas physics, namely radiative cooling and a possible preheating of gas prior to cluster formation, is able to bring the temperature-mass relation in the innermost parts of clusters into good agreement with recent determinations by Allen, Schmidt & Fabian using Chandra data.

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 use the billion-particle Hubble Volume simulations to make statistical predictions for the distribution of galaxy clusters that will be observed by the Planck Surveyor satellite through their effect on the cosmic microwave background -- the Sunyaev-Zel'dovich effect. We utilize the lightcone datasets for both critical density (tauCDM) and flat low-density (LambdaCDM) cosmologies: a `full-sky' survey out to $z \sim 0.5$, two `octant' datasets out to beyond $z=1$ and a 100 square degree dataset extending to $z \sim 4$. Making simple, but robust, assumptions regarding both the thermodynamic state of the gas and the detection of objects against an unresolved background, we present the expected number of SZ sources as a function of redshift and angular size, and also by flux (for both the thermal and kinetic effects) for 3 of the relevant HFI frequency channels. We confirm the expectation that Planck will detect around $5\times 10^4$ clusters, though the exact number is sensitive to the choice of several parameters including the baryon fraction, and also to the cluster density profile, so that either cosmology may predict more clusters. We also find that the majority of detected sources should be at $z<1.5$, and we estimate that around one per cent of clusters will be spatially resolved by Planck, though this has a large uncertainty.

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.

In this paper, we have extended the entropy-driven model of cluster evolution developed by Bower (1997) in order to be able to predict the evolution of galaxy clusters for a range of cosmological scenarios. We have applied this model to recent measurements of the evolution of the L_x-T normalisation and X-ray luminosity function in order to place constraints on cosmological parameters. We find that these measurements alone do not select a particular cosmological frame-work. An additional constraint is required on the effective slope of the power spectrum to break the degeneracy that exists between this and the background cosmology. We have therefore included a theoretical calculation of the \Omega_0 dependence on the power spectrum, based on the cold dark matter paradigm, which infers \Omega_0<0.55 (0.1<\Omega_0<0.7 for \Omega_0+\Lambda_0=1), at the 95% confidence level. Alternatively, an independent measurement of the slope of the power spectrum from galaxy clustering requires \Omega_0<0.6 (\Omega_0<0.65 for \Omega_0+\Lambda_0=1), again to 95% confidence. The rate of entropy evolution is insensitive to the values of \Omega_0 considered, although is sensitive to changes in the distribution of the intracluster medium.