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

We study the characteristics of galaxy protoclusters using the latest L-galaxies semi-analytic model. Searching for protoclusters on a scale of $\sim 10 \, \mathrm{cMpc}$ gives an excellent compromise between the completeness and purity of their galaxy populations, leads to high distinction from the field in overdensity space, and allows accurate determination of the descendant cluster mass. This scale is valid over a range of redshifts and selection criteria. We present a procedure for estimating, given a measured galaxy overdensity, the protocluster probability and its descendant cluster mass for a range of modelling assumptions, particularly taking into account the shape of the measurement aperture. This procedure produces lower protocluster probabilities compared to previous estimates using fixed size apertures. The relationship between AGN and protoclusters is also investigated, and shows significant evolution with redshift; at $z \sim 2$ the fraction of protoclusters traced by AGN is high, but the fraction of all AGN in protoclusters is low, whereas at $z \geqslant 5$ the fraction of protoclusters containing AGN is low, but most AGN are in protoclusters. We also find indirect evidence for the emergence of a passive sequence in protoclusters at $z \sim 2$, and note that a significant fraction of all galaxies reside in protoclusters at $z \geqslant 2$, particularly the most massive.

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.

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

We 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 study the quenching of star formation as a function of redshift, environment and stellar mass in the galaxy formation simulations of Henriques et al. (2015), which implement an updated version of the Munich semi-analytic model (L-GALAXIES) on the two Millennium Simulations after scaling to a Planck cosmology. In this model massive galaxies are quenched by AGN feedback depending on both black hole and hot gas mass, and hence indirectly on stellar mass. In addition, satellite galaxies of any mass can be quenched by ram-pressure or tidal stripping of gas and through the suppression of gaseous infall. This combination of processes produces quenching efficiencies which depend on stellar mass, host halo mass, environment density, distance to group centre and group central galaxy properties in ways which agree qualitatively with observation. Some discrepancies remain in dense regions and close to group centres, where quenching still seems too efficient. In addition, although the mean stellar age of massive galaxies agrees with observation, the assumed AGN feedback model allows too much ongoing star formation at late times. The fact that both AGN feedback and environmental effects are stronger in higher density environments leads to a correlation between the quenching of central and satellite galaxies which roughly reproduces observed conformity trends inside haloes.

We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied to data stored in the Millennium data repository. We contrast model SFHs from the Millennium Simulation with observed ones from the VESPA algorithm as applied to the SDSS-7 catalogue. The overall agreement is good, with both simulated and SDSS galaxies showing a steeper SFH with increased stellar mass. The SFHs of blue and red galaxies, however, show poor agreement between data and simulations, which may indicate that the termination of star formation is too abrupt in the models. The mean star-formation rate (SFR) of model galaxies is well-defined and is accurately modelled by a double power law at all redshifts: SFR proportional to $1/(x^{-1.39}+x^{1.33})$, where $x=(t_a-t)/3.0\,$Gyr, $t$ is the age of the stars and $t_a$ is the loopback time to the onset of galaxy formation; above a redshift of unity, this is well approximated by a gamma function: SFR proportional to $x^{1.5}e^{-x}$, where $x=(t_a-t)/2.0\,$Gyr. Individual galaxies, however, show a wide dispersion about this mean. When split by mass, the SFR peaks earlier for high-mass galaxies than for lower-mass ones, and we interpret this downsizing as a mass-dependence in the evolution of the quenched fraction: the SFHs of star-forming galaxies show only a weak mass dependence.

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.

Using the L-Galaxies semi-analytic model we simultaneously fit the HI mass function, stellar mass function and galaxy colours. We find good fits to all three observations at z = 0 and to the stellar mass function and galaxy colours at z = 2. Using Markov Chain Monte Carlo (MCMC) techniques we adjust the L-Galaxies parameters to best fit the constraining data. In order to fit the HI mass function we must greatly reduce the gas surface density threshold for star formation, thus lowering the number of low HI mass galaxies. A simultaneous reduction in the star formation efficiency prevents the over production of stellar content. A simplified model in which the surface density threshold is eliminated altogether also provides a good fit to the data. Unfortunately, these changes weaken the fit to the Kennicutt-Schmidt relation and raise the star-formation rate density at recent times, suggesting that a change to the model is required to prevent accumulation of gas onto dwarf galaxies in the local universe.

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

We present an analysis of the iron abundance in the hot gas surrounding galaxy groups and clusters. To do this, we first compile and homogenise a large dataset of 79 low-redshift (|z| = 0.03) systems (159 individual measurements) from the literature. Our analysis accounts for differences in aperture size, solar abundance, and cosmology, and scales all measurements using customised radial profiles for the temperature (T), gas density, and iron abundance (Z). We then compare this dataset to groups and clusters in the L-Galaxies galaxy evolution model. Our homogenised dataset reveals a tight T-Z relation for clusters, with a scatter in Z of only 0.10 dex and a slight negative gradient. After examining potential measurement biases, we conclude that at least some of this negative gradient has a physical origin. Our model suggests greater accretion of hydrogen in the hottest systems, via stripping of gas from infalling satellites, as a cause. At lower temperatures, L-Galaxies over-estimates Z in groups, indicating that metal-rich gas removal (via e.g. AGN feedback) is required. L-Galaxies provides a reasonable match to the observed Z in the intracluster medium (ICM) of the hottest clusters from at least z ~ 1.3 to 0.3. This is achieved without needing to modify any of the galactic chemical evolution (GCE) model parameters. However, the Z in intermediate-temperature clusters appears to be under-estimated in our model at z = 0. The merits and problems with modifying the GCE modelling to correct this are discussed.

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

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

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

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

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.

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

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

We update the treatment of chemical evolution in the Munich semi-analytic model, L-GALAXIES. Our new implementation includes delayed enrichment from stellar winds, supernovae type II (SNe-II) and supernovae type Ia (SNe-Ia), as well as metallicity-dependent yields and a reformulation of the associated supernova feedback. Two different sets of SN-II yields and three different SN-Ia delay-time distributions (DTDs) are considered, and eleven heavy elements (including O, Mg and Fe) are self-consistently tracked. We compare the results of this new implementation with data on a) local, star-forming galaxies, b) Milky Way disc G dwarfs, and c) local, elliptical galaxies. We find that the z=0 gas-phase mass-metallicity relation is very well reproduced for all forms of DTD considered, as is the [Fe/H] distribution in the Milky Way disc. The [O/Fe] distribution in the Milky Way disc is best reproduced when using a DTD with less than or equal to 50 per cent of SNe-Ia exploding within ~400 Myrs. Positive slopes in the mass-[alpha/Fe] relations of local ellipticals are also obtained when using a DTD with such a minor `prompt' component. Alternatively, metal-rich winds that drive light alpha elements directly out into the circumgalactic medium also produce positive slopes for all forms of DTD and SN-II yields considered. Overall, we find that the best model for matching the wide range of observational data considered here should include a power-law SN-Ia DTD, SN-II yields that take account of prior mass loss through stellar winds, and some direct ejection of light alpha elements out of galaxies.

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

We introduce a new physical recipe into the De Lucia and Blaizot version of the Munich semi-analytic model built upon the Millennium dark matter simulation: the tidal stripping of stellar material from satellite galaxies during mergers. To test the significance of the new physical process we apply a Monte Carlo Markov Chain parameter estimation technique constraining the model with the $K$-band luminosity function, $B-V$ colours and the black hole-bulge mass relation. The differences in parameter correlations, and in the allowed regions in likelihood space, reveal the impact of the new physics on the basic ingredients of the model, such as the star-formation laws, feedback recipes and the black hole growth model. With satellite disruption in place, we get a model likelihood four times higher than in the original model, indicating that the new process seems to be favoured by observations. This is achieved mainly due to a reduction in black hole growth that produces a better agreement between the properties of central black holes and host galaxies. Compared to the best-fit model without disruption, the new model removes the excess of dwarf galaxies in the original recipe with a more modest supernova heating. The new model is now consistent with the three observational data sets used to constrain it, while significantly improving the agreement with observations for the distribution of metals in stars. Moreover, the model now follows the build up of intra-cluster light.

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.

We present results of a new feedback scheme implemented in the Munich galaxy formation model. The new scheme includes a dynamical treatment of galactic winds powered by supernovae explosions and stellar winds in a cosmological context. We find that such a scheme is a good alternative to empirically-motivated recipes for feedback in galaxy formation. Model results are in good agreement with the observed luminosity functions and stellar mass function for galaxies in the local Universe. In particular, the new scheme predicts a number density of dwarfs that is lower than in previous models. This is a consequence of a new feature of the model, that allows an estimate of the amount of mass and metals that haloes can permanently deposit into the IGM. This loss of material leads to the suppression of star formation in small haloes and therefore to the decrease in the number density of dwarf galaxies. The model is able to reproduce the observed mass-stellar metallicity and luminosity-gas metallicity relationships. This demonstrates that our scheme provides a significant improvement in the treatment of the feedback in dwarf galaxies. Despite these successes, our model does not reproduce the observed bimodality in galaxy colours and predicts a larger number of bright galaxies than observed. Finally, we investigate the efficiency of metal injection in winds and in the intergalactic medium. We find that galaxies that reside in haloes with M_vir < 10^12 M_sun/h may deposit most of their metal mass into the intergalactic medium, while groups and clusters at z=0 have lost at most a few percent of their metals before the bulk of the halo mass was accreted.