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

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

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 present results of including heating and cooling in cosmological simulations of the $\Lambda$CDM cosmology and demonstrate their effects on scaling laws of galaxy clusters. The scaling relations when radiative cooling is included are in good agreement with observations but the fraction of cooled gas is on the upper limit allowed by observations. On the contrary, the preheating model has a more realistic cooled fraction but the scaling relations are less well produced.

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