It has been claimed that the standard model of cosmology (LCDM) cannot easily account for a number of observations on relatively small scales, motivating extensions to the standard model. Here we introduce a new suite of cosmological simulations that systematically explores three plausible extensions: warm dark matter, self-interacting dark matter, and a running of the scalar spectral index of density fluctuations. Current observational constraints are used to specify the additional parameters that come with these extensions. We examine a large range of observable metrics on small scales, including the halo mass function, density and circular velocity profiles, the abundance of satellite subhaloes, and halo concentrations. For any given metric, significant degeneracies can be present between the extensions. In detail, however, the different extensions have quantitatively distinct mass and radial dependencies, suggesting that a multi-probe approach over a range of scales can be used to break the degeneracies. We also demonstrate that the relative effects on the radial density profiles in the different extensions (compared to the standard model) are converged down to significantly smaller radii than are the absolute profiles. We compare the derived cosmological trends with the impact of baryonic physics using the EAGLE and ARTEMIS simulations. Significant degeneracies are also present between baryonic physics and cosmological variations (with both having similar magnitude effects on some observables). Given the inherent uncertainties both in the modelling of galaxy formation physics and extensions to LCDM, a systematic and simultaneous exploration of both is strongly warranted.

Recent analyses of the cosmic microwave background (CMB) and the Lyman-alpha forest indicate a mild preference for a deviation from a power law primordial matter power spectrum (a so-called negative `running'). We use an extension to the BAHAMAS suite of cosmological hydrodynamic simulations to explore the effects that a running scalar spectral index has on large-scale structure (LSS), using Planck CMB constraints to initialize the simulations. We focus on 5 key statistics: i) the non-linear matter power spectrum ii) the halo mass function; iii) the halo two-point auto correlation function; iv) total mass halo density profiles; and v) the halo concentration-mass relation. In terms of the matter power spectrum, we find that a running scalar spectral index affects all k-scales examined in this study, with a negative (positive) running leading to an amplification (suppression) of power. These effects should be easily detectable with upcoming surveys such as LSST and Euclid. In the mass range sampled, a positive running leads to an increase in the mass of galaxy groups and clusters, with the favoured negative running leading to a decrease in mass of lower-mass (M <~ 10^13 M_solar) halos, but an increase for the most massive (M >~ 10^13 M_solar) halos. Changes in the mass are generally confined to 5-10% which, while not insignificant, cannot by itself reconcile the claimed tension between the primary CMB and cluster number counts. We find that running does not significantly affect the shapes of density profiles of matched halos, changing only their amplitude. Finally, we demonstrate that the observed effects on LSS due to a running scalar spectral index are separable from those of baryonic effects to typically a few percent precision.

On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes; so-called baryonic feedback. A lack of knowledge about the details of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We develop a cosmology-dependent model for the matter distribution that simultaneously accounts for the clustering of dark matter, gas and stars. We inform our model by comparing it to power spectra measured from the BAHAMAS suite of hydrodynamical simulations. As well as considering matter power spectra, we also consider spectra involving the electron-pressure field, which directly relates to the thermal Sunyaev-Zel'dovich (tSZ) effect. We fit parameters in our model so that it can simultaneously model both matter and pressure data and such that the distribution of gas as inferred from tSZ has influence on the matter spectrum predicted by our model. We present two variants; one that matches the feedback-induced suppression seen in the matter-matter power spectrum at the per-cent level and a second that matches the matter-matter data slightly less well (~2 per cent), but that is able to simultaneously model the matter-electron pressure spectrum at the ~15 per-cent level. We envisage our models being used to simultaneously learn about cosmological parameters and the strength of baryonic feedback using a combination of tSZ and lensing auto- and cross-correlation data.

Dark matter (DM) direct detection experiments aim to place constraints on the DM--nucleon scattering cross-section and the DM particle mass. These constraints depend sensitively on the assumed local DM density and velocity distribution function. While astrophysical observations can inform the former (in a model-dependent way), the latter is not directly accessible with observations. Here we use the high-resolution ARTEMIS cosmological hydrodynamical simulation suite of 42 Milky Way-mass halos to explore the spatial and kinematical distributions of the DM in the solar neighbourhood, and we examine how these quantities are influenced by substructures, baryons, the presence of dark discs, as well as general halo-to-halo scatter (cosmic variance). We also explore the accuracy of the standard Maxwellian approach for modelling the velocity distribution function. We find significant halo-to-halo scatter in the density and velocity functions which, if propagated through the standard halo model for predicting the DM detection limits, implies a significant scatter about the typically quoted limit. We also show that, in general, the Maxwellian approximation works relatively well for simulations that include the important gravitational effects of baryons, but is less accurate for collisionless (DM-only) simulations. Given the significant halo-to-halo scatter in quantities relevant for DM direct detection, we advocate propagating this source of uncertainty through in order to derive conservative DM detection limits.

The standard model of cosmology, the LCDM model, robustly predicts the existence of a multitude of dark matter 'subhaloes' around galaxies like the Milky Way. A wide variety of observations have been proposed to look for the gravitational effects such subhaloes would induce in observable matter. Most of these approaches pertain to the stellar or cool gaseous phases of matter. Here we propose a new approach, which is to search for the perturbations that such dark subhaloes would source in the warm/hot circumgalactic medium (CGM) around normal galaxies. With a combination of analytic theory, carefully-controlled high-resolution idealised simulations, and full cosmological hydrodynamical simulations, we calculate the expected signal and how it depends on important physical parameters (subhalo mass, CGM temperature, and relative velocity). We find that dark subhaloes enhance the local CGM pressure, density, and temperature, in order of decreasing magnitude of the effect. For the pressure, the fluctuations can vary in magnitude from tens of percent (for subhaloes with M_sub=10^10 Msun) to a few percent (for subhaloes with M_sub=10^8 Msun), although this depends strongly on the CGM temperature. The subhaloes also induce fluctuations in the velocity field ranging in magnitude from a few km/s up to 25 km/s. We propose that X-ray, Sunyaev-Zel'dovich effect, radio dispersion measure, and quasar absorption line observations can be used to measure these fluctuations and place constraints on the abundance and distribution of dark subhaloes, thereby placing constraints on the nature of dark matter.

A large body of work based on collisionless cosmological N-body simulations going back over two decades has advanced the idea that collapsed dark matter haloes have simple and approximately universal forms for their mass density and pseudo-phase space density (PPSD) distributions. However, a general consensus on the physical origin of these results has not yet been reached. In the present study, we explore to what extent the apparent universality of these forms holds when we vary the initial conditions (i.e., the primordial power spectrum of density fluctuations) away from the standard CMB-normalised case, but still within the context of LCDM with a fixed expansion history. Using simulations that vary the initial amplitude and shape, we show that the structure of dark matter haloes retains a clear memory of the initial conditions. Specifically, increasing (lowering) the amplitude of fluctuations increases (decreases) the concentration of haloes and, if pushed far enough, the density profiles deviate strongly from the NFW form that is a good approximation for the CMB-normalised case. Although, an Einasto form works well. Rather than being universal, the slope of the PPSD (or pseudo-entropy) profile steepens (flattens) with increasing (decreasing) power spectrum amplitude and can exhibit a strong halo mass dependence. Our results therefore indicate that the previously identified universality of the structure of dark matter haloes is mostly a consequence of adopting a narrow range of (CMB-normalised) initial conditions for the simulations. Our new suite provides a useful test-bench against which physical models for the origin of halo structure can be validated.

The observable properties of galaxies depend on both internal processes and the external environment. In terms of the environmental role, we still do not have a clear picture of the processes driving the transformation of galaxies. The use of proxies for environment (e.g., host halo mass, distance to the N^th nearest neighbour, etc.), as opposed to the real physical conditions (e.g., hot gas density) may bear some responsibility for this. Here we propose a new method that directly links galaxies to their local environments, by using spatial cross-correlations of galaxy catalogues with maps from large-scale structure surveys (e.g., thermal Sunyaev-Zel'dovich [tSZ] effect, diffuse X-ray emission, weak lensing of galaxies or the CMB). We focus here on the quenching of galaxies and its link to local hot gas properties. Maps of galaxy overdensity and quenched fraction excess are constructed from volume-limited SDSS catalogs, which are cross-correlated with tSZ effect and X-ray maps from Planck and ROSAT, respectively. Strong signals out to Mpc scales are detected for most cross-correlations and are compared to predictions from the EAGLE and BAHAMAS cosmological hydrodynamical simulations. The simulations successfully reproduce many, but not all, of the observed power spectra, with an indication that environmental quenching may be too efficient in the simulations. We demonstrate that the cross-correlations are sensitive to both the internal (e.g., AGN and stellar feedback) and external processes (e.g., ram pressure stripping, harassment, strangulation, etc.) responsible for quenching. The methods outlined in this paper can be adapted to other observables and, with upcoming surveys, will provide a stringent test of physical models for environmental transformation.

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

We evaluate the ability of Convolutional Neural Networks (CNNs) to predict galaxy cluster masses in the BAHAMAS hydrodynamical simulations. We train four separate single-channel networks using: stellar mass, soft X-ray flux, bolometric X-ray flux, and the Compton $y$ parameter as observational tracers, respectively. Our training set consists of $\sim$6400 synthetic cluster images generated from the simulation, while an additional $\sim$1600 images form a test set. We also train a `multi-channel' CNN by combining the four observational tracers. We utilize $\texttt{Keras}$ with a $\texttt{Tensorflow}$ backend to train the network, and all four converge within 2000 epochs. The cluster masses predicted from these networks are evaluated using the average fractional difference between predicted cluster mass and true cluster mass. The resulting predictions are especially precise for halo masses in the range $10^{13.25}M_{\odot}<M<10^{14.5}M_{\odot}$, where all five networks produce mean mass biases of order $\approx 1\%$ with a scatter on the mean bias of $\approx 0.5\%$. The network trained with Compton $y$ parameter maps yields the most precise predictions. We interpret the network's behaviour using two diagnostic tests to determine which features are used to predict cluster mass. The CNN trained with stellar mass images detect galaxies (not surprisingly), while CNNs trained with gas-based tracers utilise the shape of the signal to estimate cluster mass.

In this work we consider the impact of spatially-uniform but time-varying dark energy (or `dynamical dark energy', DDE) on large-scale structure in a spatially flat universe, using large cosmological hydrodynamical simulations that form part of the BAHAMAS project. As DDE changes the expansion history of the universe, it impacts the growth of structure. We explore variations in DDE that are constrained to be consistent with the cosmic microwave background. We find that DDE can affect the clustering of matter and haloes at the ~10% level (suppressing it for so-called `freezing' models, while enhancing it for `thawing' models), which should be distinguishable with upcoming large-scale structure surveys. DDE cosmologies can also enhance or suppress the halo mass function (with respect to LCDM) over a wide range of halo masses. The internal properties of haloes are minimally affected by changes in DDE, however. Finally, we show that the impact of baryons and associated feedback processes is largely independent of the change in cosmology and that these processes can be modelled separately to typically better than a few percent accuracy

We introduce the ARTEMIS simulations, a new set of 42 zoomed-in, high-resolution (baryon particle mass of ~ 2x10^4 Msun/h), hydrodynamical simulations of galaxies residing in haloes of Milky Way mass, simulated with the EAGLE galaxy formation code with re-calibrated stellar feedback. In this study, we analyse the structure of stellar haloes, specifically the mass density, surface brightness, metallicity, colour and age radial profiles, finding generally very good agreement with recent observations of local galaxies. The stellar density profiles are well fitted by broken power laws, with inner slopes of ~ -3, outer slopes of ~ -4 and break radii that are typically ~ 20-40 kpc. The break radii generally mark the transition between in situ formation and accretion-driven formation of the halo. The metallicity, colour and age profiles show mild large-scale gradients, particularly when spherically-averaged or viewed along the major axes. Along the minor axes, however, the profiles are nearly flat, in agreement with observations. Overall, the structural properties can be understood by two factors: that in situ stars dominate the inner regions and that they reside in a spatially-flattened distribution that is aligned with the disc. Observations targeting both the major and minor axes of galaxies are thus required to obtain a complete picture of stellar haloes.

We present results on the environmental dependence of the star-forming galaxy main sequence in 11 galaxy cluster fields at $1.0 < z < 1.5$ from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN) survey. We use a homogeneously selected sample of field and cluster galaxies whose membership is derived from dynamical analysis. Using [OII]-derived star formation rates (SFRs), we find that cluster galaxies have suppressed SFRs at fixed stellar mass in comparison to their field counterparts by a factor of 1.4 $\pm$ 0.1 ($\sim3.3\sigma$) across the stellar mass range: $9.0 < \log(M_{*} /M_{\odot}) < 11.2$. We also find that this modest suppression in the cluster galaxy star-forming main sequence is mass and redshift dependent: the difference between cluster and field increases towards lower stellar masses and lower redshift. When comparing the distribution of cluster and field galaxy SFRs to the star-forming main sequence, we find an overall shift towards lower SFRs in the cluster population, and note the absence of a tail of high SFR galaxies as seen in the field. Given this observed suppression in the cluster galaxy star-forming main sequence, we explore the implications for several scenarios such as formation time differences between cluster and field galaxies, and environmentally-induced star formation quenching and associated timescales.

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

The location of a galaxy cluster's centroid is typically derived from observations of the galactic and/or gas component of the cluster, but these typically deviate from the true centre. This can produce bias when observations are combined to study average cluster properties. Using data from the BAHAMAS cosmological hydrodynamic simulations we study this bias in both two and three dimensions for 2000 clusters over the $10^{13} - 10^{15} ~\mathrm{M_{\odot}}$ mass range. We quantify and model the offset distributions between observationally-motivated centres and the `true' centre of the cluster, which is taken to be the most gravitationally bound particle measured in the simulation. We fit the cumulative distribution function of offsets with an exponential distribution and a Gamma distribution fit well with most of the centroid definitions. The galaxy-based centres can be seen to be divided into a mis-centred group and a well-centred group, with the well-centred group making up about $60\%$ of all the clusters. Gas-based centres are overall less scattered than galaxy-based centres. We also find a cluster-mass dependence of the offset distribution of gas-based centres, with generally larger offsets for smaller mass clusters. We then measure cluster density profiles centred at each choice of the centres and fit them with empirical models. Stacked, mis-centred density profiles fit to the Navarro-Frenk-White dark-matter profile and Komatsu-Seljak gas profile show that recovered shape and size parameters can significantly deviate from the true values. For the galaxy-based centres, this can lead to cluster masses being underestimated by up to $10\%$.

We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $\gamma = -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $\gamma \simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$.

Running hydrodynamical simulations to produce mock data of large-scale structure and baryonic probes, such as the thermal Sunyaev-Zeldovich (tSZ) effect, at cosmological scales is computationally challenging. We propose to leverage the expressive power of deep generative models to find an effective description of the large-scale gas distribution and temperature. We train two deep generative models, a variational auto-encoder and a generative adversarial network, on pairs of matter density and pressure slices from the BAHAMAS hydrodynamical simulation. The trained models are able to successfully map matter density to the corresponding gas pressure. We then apply the trained models on 100 lines-of-sight from SLICS, a suite of N-body simulations optimised for weak lensing covariance estimation, to generate maps of the tSZ effect. The generated tSZ maps are found to be statistically consistent with those from BAHAMAS. We conclude by considering a specific observable, the angular cross-power spectrum between the weak lensing convergence and the tSZ effect and its variance, where we find excellent agreement between the predictions from BAHAMAS and SLICS, thus enabling the use of SLICS for tSZ covariance estimation.

We construct the mean thermal Sunyaev-Zel'dovich (tSZ) Comptonization y profile around Luminous Red Galaxies (LRGs) in the redshift range 0.16 < z < 0.47 from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) using the Planck y map. The mean central tSZ signal for the full sample is y ~ 1.8 * 10^(-7) and we detect tSZ emission out to ~30 arcmin, which is well beyond the 10 arcmin angular resolution of the y map and well beyond the virial radii of the LRGs. We compare the measured profile with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations. This comparison agrees well for models that include feedback from active galactic nuclei (AGN), but not with hydrodynamic models without this energetic feedback mechanism. This suggests that an additional heating mechanism is required over SNe feedback and star formation to explain the y data profile. We also compare our results with predictions based on the halo model with a universal pressure profile (UPP) giving the y signal. The predicted profile is consistent with the data, but only if we account for the clustering of haloes via a two-halo term and if halo masses are estimated using the mean stellar-to-halo mass (SHM) relation of Coupon et al. (2015) or Wang et al.(2016) estimated from gravitational lensing measurements. We also discuss the importance of scatter in the SHM relation on the model predictions.

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 search the Planck data for a thermal Sunyaev-Zel'dovich (tSZ) signal due to gas filaments between pairs of Luminous Red Galaxies (LRG's) taken from the Sloan Digital Sky Survey Data Release 12 (SDSS/DR12). We identify $\sim$260,000 LRG pairs in the DR12 catalog that lie within 6-10 $h^{-1} \mathrm{Mpc}$ of each other in tangential direction and within 6 $h^{-1} \mathrm{Mpc}$ in radial direction. We stack pairs by rotating and scaling the angular positions of each LRG so they lie on a common reference frame, then we subtract a circularly symmetric halo from each member of the pair to search for a residual signal between the pair members. We find a statistically significant (5.3$\sigma$) signal between LRG pairs in the stacked data with a magnitude $\Delta y = (1.31 \pm 0.25) \times 10^{-8}$. The uncertainty is estimated from two Monte Carlo null tests which also establish the reliability of our analysis. Assuming a simple, isothermal, cylindrical filament model of electron over-density with a radial density profile proportional to $r_c/r$ (as determined from simulations), where $r$ is the perpendicular distance from the cylinder axis and $r_c$ is the core radius of the density profile, we constrain the product of over-density and filament temperature to be $\delta_c \times (T_{\rm e}/10^7 \, {\rm K}) \times (r_c/0.5h^{-1} \, {\rm Mpc}) = 2.7 \pm 0.5$. To our knowledge, this is the first detection of filamentary gas at over-densities typical of cosmological large-scale structure. We compare our result to the BAHAMAS suite of cosmological hydrodynamic simulations (McCarthy et al. 2017) and find a slightly lower, but marginally consistent Comptonization excess, $\Delta y = (0.84 \pm 0.24) \times 10^{-8}$.

We present a novel statistic to extract cosmological information in weak lensing data: the lensing minima. We also investigate the effect of baryons on the cosmological constraints from peak and minimum counts. Using the \texttt{MassiveNuS} simulations, we find that lensing minima are sensitive to non-Gaussian cosmological information and are complementary to the lensing power spectrum and peak counts. For an LSST-like survey, we obtain $95\%$ credible intervals from a combination of lensing minima and peaks that are significantly stronger than from the power spectrum alone, by $44\%$, $11\%$, and $63\%$ for the neutrino mass sum $\sum m_\nu$, matter density $\Omega_m$, and amplitude of fluctuation $A_s$, respectively. We explore the effect of baryonic processes on lensing minima and peaks using the hydrodynamical simulations \texttt{BAHAMAS} and \texttt{Osato15}. We find that ignoring baryonic effects would lead to strong ($\approx 4 \sigma$) biases in inferences from peak counts, but negligible ($\approx 0.5 \sigma$) for minimum counts, suggesting lensing minima are a potentially more robust tool against baryonic effects. Finally, we demonstrate that the biases can in principle be mitigated without significantly degrading cosmological constraints when we model and marginalize the baryonic effects.

Upcoming weak lensing surveys require a detailed theoretical understanding of the matter power spectrum in order to derive accurate and precise cosmological parameter values. While galaxy formation is known to play an important role, its precise effects are currently unknown. We present a set of 92 matter power spectra from the OWLS, cosmo-OWLS and BAHAMAS simulation suites, including different $\Lambda$CDM cosmologies, neutrino masses, subgrid prescriptions and AGN feedback strengths. We conduct a detailed investigation of the dependence of the relative difference between the total matter power spectra in hydrodynamical and collisionless simulations on the effectiveness of stellar and AGN feedback, cosmology and redshift. The strength of AGN feedback can greatly affect the power on a range of scales, while a lack of stellar feedback can greatly increase the effectiveness of AGN feedback on large scales. We also examine differences in the initial conditions of hydrodynamic and N-body simulations that can lead to a ~1% discrepancy in the large-scale power, and furthermore show our results to be insensitive to cosmic variance. We present an empirical model capable of predicting the effect of galaxy formation on the matter power spectrum at z=0 to within 1% for k<1 h/Mpc, given only the mean baryon fraction in galaxy groups. Differences in group baryon fractions can also explain the quantitative disagreement between predictions from the literature. All total and dark matter only power spectra in this library will be made publicly available at powerlib.strw.leidenuniv.nl.

We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray-tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass $M_{\nu} =$ 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of $M_{\nu}$. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations.

Models of Cold Dark Matter predict that the distribution of dark matter in galaxy clusters should be cuspy, centrally concentrated. Constant density cores would be strong evidence for beyond-CDM physics, such as Self-Interacting Dark Matter (SIDM). An observable consequence would be oscillations of the Brightest Cluster Galaxy (BCG) in otherwise relaxed galaxy clusters. Offset BCGs have indeed been observed - but only interpreted via a simplified, analytic model of oscillations. We compare these observations to the BAHAMAS-SIDM suite of cosmological simulations, which include SIDM and a fully hydrodynamical treatment of star formation and feedback. We predict that the median offset of BCGs increases with the SIDM cross-section, cluster mass and the amount of stellar mass within 10kpc, while CDM exhibits no trend in mass. Interpolating between the simulated cross-sections, we find that the observations (of 10 clusters) is consistent with CDM at the ~1.5$\sigma$ level, and prefer cross-section $\sigma$/m < 0.12(0.39)cm$^2$/g at 68% (95%) confidence level. This is on the verge of ruling out velocity-independent dark matter self-interactions as the solution to discrepancies between the predicted and observed behaviour of dwarf galaxies, and will be improved by larger surveys by Euclid or SuperBIT.

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the expected statistical uncertainty of upcoming experiments at separations under 10 Mpc. Successful extraction of information in such a regime thus requires, at the very least, unbiased models for the impact of galaxy formation on the matter distribution, and can benefit from complementary observational priors. This work reviews the current state of the art in the modelling of baryons for cosmology, from numerical methods to approximate analytical prescriptions, and makes recommendations for studies in the next decade, including a discussion of potential probe combinations that can help constrain the role of baryons in cosmological studies. We focus, in particular, on the modelling of the matter power spectrum, $P(k,z)$, as a function of scale and redshift, and of the observables derived from this quantity. This work is the result of a workshop held at the University of Oxford in November of 2018.

We examine the origin of scatter in the relationship between the gas fraction and mass of dark matter haloes hosting present-day $\sim L^\star$ central galaxies in the EAGLE simulations. The scatter is uncorrelated with the accretion rate of the central galaxy's black hole (BH), but correlates strongly and negatively with the BH's mass, implicating differences in the expulsion of gas by active galactic nucleus feedback, throughout the assembly of the halo, as the main cause of scatter. Haloes whose central galaxies host undermassive BHs also tend to retain a higher gas fraction, and exhibit elevated star formation rates (SFRs). Diversity in the mass of central BHs stems primarily from diversity in the dark matter halo binding energy, as these quantities are strongly and positively correlated at fixed halo mass, such that $\sim L^\star$ galaxies hosted by haloes that are more (less) tightly-bound develop central BHs that are more (less) massive than is typical for their halo mass. Variations in the halo gas fraction at fixed halo mass are reflected in both the soft X-ray luminosity and thermal Sunyaev-Zel'dovich flux, suggesting that the prediction of a strong coupling between the properties of galaxies and their halo gas fractions can be tested with measurements of these diagnostics for galaxies with diverse SFRs but similar halo masses.

We investigate the connection between the morphology and internal kinematics of the stellar component of central galaxies with mass $M_\star > {10}^{9.5} {\rm M}_\odot$ in the EAGLE simulations. We compare several kinematic diagnostics commonly used to describe simulated galaxies, and find good consistency between them. We model the structure of galaxies as ellipsoids and quantify their morphology via the ratios of their principal axes. We show that the differentiation of blue star-forming and red quiescent galaxies using morphological diagnostics can be achieved with similar efficacy to the use of kinematical diagnostics, but only if one is able to measure both the flattening and the triaxiality of the galaxy. Flattened oblate galaxies exhibit greater rotational support than their spheroidal counterparts, but there is significant scatter in the relationship between morphological and kinematical diagnostics, such that kinematically-similar galaxies can exhibit a broad range of morphologies. The scatter in the relationship between the flattening and the ratio of the rotation and dispersion velocities ($v/\sigma$) correlates strongly with the anisotropy of the stellar velocity dispersion: at fixed $v/\sigma$, flatter galaxies exhibit greater dispersion in the plane defined by the intermediate and major axes than along the minor axis, indicating that the morphology of simulated galaxies is influenced significantly by the structure of their velocity dispersion. The simulations reveal that this anisotropy correlates with the intrinsic morphology of the galaxy's inner dark matter halo, i.e. the halo's morphology that emerges in the absence of dissipative baryonic physics. This implies the existence of a causal relationship between the morphologies of galaxies and that of their host dark matter haloes.

Weak gravitational lensing depends on the integrated mass along the line of sight. Baryons contribute to the mass distribution of galaxy clusters and the resulting mass estimates from lensing analysis. We use the cosmo-OWLS suite of hydrodynamic simulations to investigate the impact of baryonic processes on the bias and scatter of weak lensing mass estimates of clusters. These estimates are obtained by fitting NFW profiles to mock data using MCMC techniques. In particular, we examine the difference in estimates between dark matter-only runs and those including various prescriptions for baryonic physics. We find no significant difference in the mass bias when baryonic physics is included, though the overall mass estimates are suppressed when feedback from AGN is included. For lowest-mass systems for which a reliable mass can be obtained ($M_{200} \approx 2 \times 10^{14}$ $M_{\odot}$), we find a bias of $\approx -10$ per cent. The magnitude of the bias tends to decrease for higher mass clusters, consistent with no bias for the most massive clusters which have masses comparable to those found in the CLASH and HFF samples. For the lowest mass clusters, the mass bias is particularly sensitive to the fit radii and the limits placed on the concentration prior, rendering reliable mass estimates difficult. The scatter in mass estimates between the dark matter-only and the various baryonic runs is less than between different projections of individual clusters, highlighting the importance of triaxiality.

For more than two decades, the Navarro, Frenk, and White (NFW) model has stood the test of time; it has been used to describe the distribution of mass in galaxy clusters out to their outskirts. Stacked weak lensing measurements of clusters are now revealing the distribution of mass out to and beyond their virial radii, where the NFW model is no longer applicable. In this study we assess how well the parameterised Diemer & Kravstov (DK) density profile describes the characteristic mass distribution of galaxy clusters extracted from cosmological simulations. This is determined from stacked synthetic lensing measurements of the 50 most massive clusters extracted from the Cosmo-OWLS simulations, using the Dark Matter Only run and also the run that most closely matches observations. The characteristics of the data reflect the Weighing the Giants survey and data from the future Large Synoptic Survey Telescope (LSST). In comparison with the NFW model, the DK model favored by the stacked data, in particular for the future LSST data, where the number density of background galaxies is higher. The DK profile depends on the accretion history of clusters which is specified in the current study. Eventually however subsamples of galaxy clusters with qualities indicative of disparate accretion histories could be studied.

Recent studies have presented evidence for tension between the constraints on Omega_m and sigma_8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully-calibrated cosmological hydrodynamical simulations, which thus capture the back reaction of the baryons on the total matter distribution. We have extended the BAHAMAS simulations to include a treatment of massive neutrinos, which currently represents the best motivated extension to the standard model. We make synthetic thermal Sunyaev-Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: i) in general there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and iii) if one accounts for internal tensions in the Planck CMB dataset (by allowing the lensing amplitude, A_Lens, to vary), invoking a non-minimal neutrino mass, typically of 0.2-0.4 eV, can resolve the tension. This solution is fully consistent with separate constraints from the primary CMB and baryon acoustic oscillations.

We study the scaling relations between the baryonic content and total mass of groups of galaxies, as these systems provide a unique way to examine the role of non-gravitational processes in structure formation. Using Planck and ROSAT data, we conduct detailed comparisons of the stacked thermal Sunyaev-Zel'dovich (tSZ) effect and X-ray scaling relations of galaxy groups found in the the Galaxy And Mass Assembly (GAMA) survey and the BAHAMAS hydrodynamical simulation. We use weak gravitational lensing data from the Kilo Degree Survey (KiDS) to determine the average halo mass of the studied systems. We analyse the simulation in the same way, using realistic weak lensing, X-ray, and tSZ synthetic observations. Furthermore, to keep selection biases under control, we employ exactly the same galaxy selection and group identification procedures to the observations and simulation. Applying this careful comparison, we find that the simulations are in agreement with the observations, particularly with regards to the scaling relations of the lensing and tSZ results. This finding demonstrates that hydrodynamical simulation have reached the level of realism that is required to interpret observational survey data and study the baryon physics within dark matter haloes, where analytical modelling is challenging. Finally, using simulated data, we demonstrate that our observational processing of the X-ray and tSZ signals is free of significant biases. We find that our optical group selection procedure has, however, some room for improvement.

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.

[Abridged] Typical disc galaxies forming in a LambdaCDM cosmology encounter a violent environment, where they often experience mergers with massive satellites. The fact that disc galaxies are ubiquitous in the local Universe suggests that a quiescent history is not necessary for their formation. Modern cosmological simulations can now obtain relatively realistic populations of disc galaxies, but it still remains to be clarified how discs manage to survive massive mergers. Here we use a suite of high-resolution hydrodynamical simulations set in a LambdaCDM cosmology to elucidate the fate of discs encountering massive mergers. We extract a sample of approximately 100 disc galaxies and follow the changes in their post-merger morphologies, as tracked by their disc-to-total ratios (D/T). We also examine the relations between their present-day morphology, assembly history and gas fractions. We find that approximately half of present-day disc galaxies underwent at least one merger with a satellite of total mass exceeding the host system's stellar mass, a third had mergers with satellites of mass exceeding 3 times the host's stellar mass, and approximately one-sixth had mergers with satellites of mass exceeding 10 times of the host's stellar mass. These mergers lead to a sharp, but often temporary, decrease in the D/T of the hosts, implying that discs are usually disrupted but then quickly re-grow. To do so, high cold gas fractions are required post-merger, as well as a relatively quiescent recent history (over a few Gyrs before z=0). Our results show that discs can form via diverse merger pathways and that quiescent histories are not the dominant mode of disc formation.

Modern cosmological simulations rely heavily on feedback from active galactic nuclei (AGN) in order to stave off overcooling in massive galaxies and galaxy groups and clusters. An important independent test is whether or not the simulations capture the broad demographics of the observed AGN population. Here, we have used the cosmo-OWLS suite of cosmological hydrodynamical simulations to produce realistic synthetic catalogs of X-ray AGN out to $z$=3, with the aim of comparing the catalogs to the observed X-ray AGN population in the XXL survey and other recent surveys. We focused on the unabsorbed X-ray luminosity function (XLF), the Eddington ratio distribution, the black hole mass function, and the projected clustering of X-ray AGN. To compute the unabsorbed XLF of the simulated AGN, we used recent empirically-determined bolometric corrections. We show that the simulated AGN sample accurately reproduces the observed XLF over 3 orders of magnitude in X-ray luminosity in all redshift bins. To compare to the observed Eddington ratio distribution and the clustering of AGN, we produced detailed 'XMM-Newton-detected' catalogs of the simulated AGN. This requires the production of synthetic X-ray images extracted from light cones of the simulations that fold in the relevant instrumental effects of XMM-Newton. We apply a luminosity- and redshift-dependent obscuration function for the AGN and employ the same AGN detection algorithm as used for the real XXL survey. We demonstrate that the detected population of simulated AGN reproduces the observed Eddington ratio distribution and projected clustering from XXL quite well. We conclude that the simulations have a broadly realistic population of AGN and that our synthetic X-ray AGN catalogs should be useful for interpreting additional trends and as a helpful tool for quantifying AGN contamination in galaxy group and cluster X-ray surveys.

A striking signal of dark matter beyond the standard model is the existence of cores in the centre of galaxy clusters. Recent simulations predict that a Brightest Cluster Galaxy (BCG) inside a cored galaxy cluster will exhibit residual wobbling due to previous major mergers, long after the relaxation of the overall cluster. This phenomenon is absent with standard cold dark matter where a cuspy density profile keeps a BCG tightly bound at the centre. We test this hypothesis using cosmological simulations and deep observations of 10 galaxy clusters acting as strong gravitational lenses. Modelling the BCG wobble as a simple harmonic oscillator, we measure the wobble amplitude, $A_w$, in the BAHAMAS suite of cosmological hydrodynamical simulations, finding an upper limit for the CDM paradigm of $A_w < 2$kpc at the 95% confidence limit. We carry out the same test on the data finding a non-zero amplitude of $A_w=11.82^{+7.3}_{-3.0}$kpc, with the observations dis-favouring $A_w = 0$ at the 3$\sigma$ confidence level. This detection of BCG wobbling is evidence for a dark matter core at the heart of galaxy clusters. It also shows that strong lensing models of clusters cannot assume that the BCG is exactly coincident with the large scale halo. While our small sample of galaxy clusters already indicates a non-zero $A_w$, with larger surveys, e.g. Euclid, we will be able to not only to confirm the effect but also to use it to determine whether or not the wobbling finds its origin in new fundamental physics or astrophysical process.

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.

Observations with the Chandra X-ray Observatory are used to examine the hot gas properties within a sample of 10 galaxy groups selected from the Galaxy And Mass Assembly survey's optical Friends-of-Friends group catalogue. Our groups have been screened to eliminate spurious and unrelaxed systems, and the effectiveness of this procedure is demonstrated by the detection of intergalactic hot gas in 80 per cent of our sample. However, we find that 9 of the 10 are X-ray underluminous by a mean factor of $\sim$4 compared to typical X-ray-selected samples. Consistent with this, the majority of our groups have gas fractions that are lower and gas entropies somewhat higher than those seen in typical X-ray-selected samples. Two groups, which have high 2{\sigma} lower limits on their gas entropy, are candidates for the population of high-entropy groups predicted by some active galactic nucleus feedback models.

The Gaia astrometric mission may offer an unprecedented opportunity to discover new tidal streams in the Galactic halo. To test this, we apply nGC3, a great-circle-cell count method that combines position and proper motion data to identify streams, to ten mock Gaia catalogues of K giants and RR Lyrae stars constructed from cosmological simulations of Milky Way analogues. We analyse two sets of simulations, one using a combination of $N$-body and semi-analytical methods which has extremely high resolution, the other using hydro-dynamical methods, which captures the dynamics of baryons, including the formation of an in situ halo. These ten realisations of plausible Galactic merger histories allow us to assess the potential for the recovery of tidal streams in different Milky Way formation scenarios. We include the Gaia~selection function and observational errors in these mock catalogues. We find that the nGC3 method has a well-defined detection boundary in the space of stream width and projected overdensity, that can be predicted based on direct observables alone. We predict that about 4-13 dwarf galaxy streams can be detected in a typical Milky Way-mass halo with Gaia+nGC3, with an estimated efficiency of $>$80\% inside the detection boundary. The progenitors of these streams are in the mass range of the classical dwarf galaxies and may have been accreted as early as redshift $\sim3$. Finally, we analyse how different possible extensions of the Gaia mission will improve the detection of tidal streams.

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 use the cosmo-OWLS and BAHAMAS suites of cosmological hydrodynamical simulations to explore the separate and combined effects of baryon physics (particularly feedback from active galactic nuclei, AGN) and free-streaming of massive neutrinos on large-scale structure. We focus on five diagnostics: i) the halo mass function; ii) halo mass density profiles; iii) the halo mass-concentration relation; iv) the clustering of haloes; and v) the clustering of matter; and we explore the extent to which the effects of baryon physics and neutrino free-streaming can be treated independently. Consistent with previous studies, we find that both AGN feedback and neutrino free-streaming suppress the total matter power spectrum, although their scale and redshift dependencies differ significantly. The inclusion of AGN feedback can significantly reduce the masses of groups and clusters, and increase their scale radii. These effects lead to a decrease in the amplitude of the mass-concentration relation and an increase in the halo autocorrelation function at fixed mass. Neutrinos also lower the masses of groups and clusters while having no significant effect on the shape of their density profiles (thus also affecting the mass-concentration relation and halo clustering in a qualitatively similar way to feedback). We show that, with only a small number of exceptions, the combined effects of baryon physics and neutrino free-streaming on all five diagnostics can be estimated to typically better than a few percent accuracy by treating these processes independently (i.e., by multiplying their separate effects).

Mass is a fundamental property of galaxy groups and clusters. In theory weak gravitational lensing will enable an approximately unbiased measurement of mass, but parametric methods for extracting cluster masses from data require the additional knowledge of concentration. Measurements of both mass and concentration are limited by the degeneracy between the two parameters, particularly in low mass, high redshift systems where the signal-to-noise is low. In this paper we develop a hierarchical model of mass and concentration for mass inference we test our method on toy data and then apply it to a sample of galaxy groups and poor clusters down to masses of $\sim$ 1e13 M$_\odot$. Our fit and model gives a relationship among masses, concentrations and redshift that allow prediction of these parameters from incomplete and noisy future measurements. Additionally the underlying population can be used to infer an observationally based concentration-mass relation. Our method is equivalent to a quasi- stacking approach with the degree of stacking set by the data. We also demonstrate that mass and concentration derived from pure stacking can be offset from the population mean with differing values depending on the method of stacking.

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations to investigate the scatter and evolution of the global hot gas properties of large simulated populations of galaxy groups and clusters. Our aim is to compare the predictions of different physical models and to explore the extent to which commonly-adopted assumptions in observational analyses (e.g. self-similar evolution) are violated. We examine the relations between (true) halo mass and the X-ray temperature, X-ray luminosity, gas mass, Sunyaev-Zel'dovich (SZ) flux, the X-ray analogue of the SZ flux ($Y_X$) and the hydrostatic mass. For the most realistic models, which include AGN feedback, the slopes of the various mass-observable relations deviate substantially from the self-similar ones, particularly at late times and for low-mass clusters. The amplitude of the mass-temperature relation shows negative evolution with respect to the self-similar prediction (i.e. slower than the prediction) for all models, driven by an increase in non-thermal pressure support at higher redshifts. The AGN models predict strong positive evolution of the gas mass fractions at low halo masses. The SZ flux and $Y_X$ show positive evolution with respect to self-similarity at low mass but negative evolution at high mass. The scatter about the relations is well approximated by log-normal distributions, with widths that depend mildly on halo mass. The scatter decreases significantly with increasing redshift. The exception is the hydrostatic mass-halo mass relation, for which the scatter increases with redshift. Finally, we discuss the relative merits of various hot gas-based mass proxies.

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

The evolution of the large-scale distribution of matter is sensitive to a variety of fundamental parameters that characterise the dark matter, dark energy, and other aspects of our cosmological framework. Since the majority of the mass density is in the form of dark matter that cannot be directly observed, to do cosmology with large-scale structure one must use observable (baryonic) quantities that trace the underlying matter distribution in a (hopefully) predictable way. However, recent numerical studies have demonstrated that the mapping between observable and total mass, as well as the total mass itself, are sensitive to unresolved feedback processes associated with galaxy formation, motivating explicit calibration of the feedback efficiencies. Here we construct a new suite of large-volume cosmological hydrodynamical simulations (called BAHAMAS, for BAryons and HAloes of MAssive Systems) where subgrid models of stellar and Active Galactic Nucleus (AGN) feedback have been calibrated to reproduce the present-day galaxy stellar mass function and the hot gas mass fractions of groups and clusters in order to ensure the effects of feedback on the overall matter distribution are broadly correct. We show that the calibrated simulations reproduce an unprecedentedly wide range of properties of massive systems, including the various observed mappings between galaxies, hot gas, total mass, and black holes, and represent a significant advance in our ability to mitigate the primary systematic uncertainty in most present large-scale structure tests.

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

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.

Observations of galaxies in the local Universe have shown that both the ionized gas and the stars of satellites are more metal-rich than of equally massive centrals. To gain insight into the connection between this metallicity enhancement and other differences between centrals and satellites, such as their star formation rates, gas content, and growth history, we study the metallicities of >3600 galaxies with M_star > 10^10 M_sun in the cosmological hydrodynamical EAGLE 100 Mpc `Reference' simulation, including ~1500 in the vicinity of galaxy groups and clusters (M_200 >= 10^13 M_sun). The simulation predicts excess gas and stellar metallicities in satellites consistent with observations, except for stellar metallicities at M_star <~ 10^10.2 M_sun where the predicted excess is smaller than observed. The exact magnitude of the effect depends on galaxy selection, aperture, and on whether the metallicity is weighted by stellar mass or luminosity. The stellar metallicity excess in clusters is also sensitive to the efficiency scaling of star formation feedback. We identify stripping of low-metallicity gas from the galaxy outskirts, as well as suppression of metal-poor inflows towards the galaxy centre, as key drivers of the enhancement of gas metallicity. Stellar metallicities in satellites are higher than in the field as a direct consequence of the more metal-rich star forming gas, whereas stripping of stars and suppressed stellar mass growth, as well as differences in accreted vs. in-situ star formation between satellites and the field, are of secondary importance.

We present measurements of the spatial mapping between (hot) baryons and the total matter in the Universe, via the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) map from Planck and the weak gravitational lensing maps from the Red Sequence Cluster Survey (RCSLenS). The cross-correlations are performed on the map level where all the sources (including diffuse intergalactic gas) contribute to the signal. We consider two configuration-space correlation function estimators, $\xi^{ y-\kappa}$ and $\xi^ {y-\gamma_{t}}$, and a Fourier space estimator, $C_{\ell}^{y-\kappa}$, in our analysis. We detect a significant correlation out to three degrees of angular separation on the sky. Based on statistical noise only, we can report 13$\sigma$ and 17$\sigma$ detections of the cross-correlation using the configuration-space $y-\kappa$ and $y-\gamma_{t}$ estimators, respectively. Including a heuristic estimate of the sampling variance yields a detection significance of 6$\sigma$ and 8$\sigma$, respectively. A similar level of detection is obtained from the Fourier-space estimator, $C_{\ell}^{y-\kappa}$. As each estimator probes different dynamical ranges, their combination improves the significance of the detection. We compare our measurements with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations, where different galactic feedback models are implemented. We find that a model with considerable AGN feedback that removes large quantities of hot gas from galaxy groups and WMAP-7yr best-fit cosmological parameters provides the best match to the measurements. All baryonic models in the context of a Planck cosmology over-predict the observed signal. Similar cosmological conclusions are drawn when we employ a halo model with the observed `universal' pressure profile.

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

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.

The evolution of galaxy cluster counts is a powerful probe of several fundamental cosmological parameters. A number of recent studies using this probe have claimed tension with the cosmology preferred by the analysis of the Planck primary CMB data, in the sense that there are fewer clusters observed than predicted based on the primary CMB cosmology. One possible resolution to this problem is systematic errors in the absolute halo mass calibration in cluster studies, which is required to convert the standard theoretical prediction (the halo mass function) into counts as a function of the observable (e.g., X-ray luminosity, Sunyaev-Zel'dovich flux, optical richness). Here we propose an alternative strategy, which is to directly compare predicted and observed cluster counts as a function of the one-dimensional velocity dispersion of the cluster galaxies. We argue that the velocity dispersion of groups/clusters can be theoretically predicted as robustly as mass but, unlike mass, it can also be directly observed, thus circumventing the main systematic bias in traditional cluster counts studies. With the aid of the BAHAMAS suite of cosmological hydrodynamical simulations, we demonstrate the potential of the velocity dispersion counts for discriminating even similar $\Lambda$CDM models. These predictions can be compared with the results from existing redshift surveys such as the highly-complete Galaxy And Mass Assembly (GAMA) survey, and upcoming wide-field spectroscopic surveys such as the Wide Area Vista Extragalactic Survey (WAVES) and the Dark Energy Survey Instrument (DESI).

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations, which includes different galactic feedback models, to predict the cross-correlation signal between weak gravitational lensing and the thermal Sunyaev-Zeldovich (tSZ) $y$-parameter. The predictions are compared to the recent detection reported by van Waerbeke and collaborators. The simulations reproduce the weak lensing-tSZ cross-correlation, $\xi_{y\kappa}(\theta)$, well. The uncertainty arising from different possible feedback models appears to be important on small scales only ($\theta \lesssim 10$ arcmin), while the amplitude of the correlation on all scales is sensitive to cosmological parameters that control the growth rate of structure (such as $\sigma_8$, $\Omega_m$ and $\Omega_b$). This study confirms our previous claim (in Ma et al.) that a significant proportion of the signal originates from the diffuse gas component in low-mass ($M_{\rm{halo}} \lesssim 10^{14} M_{\odot}$) clusters as well as from the region beyond the virial radius. We estimate that approximately 20$\%$ of the detected signal comes from low-mass clusters, which corresponds to about 30$\%$ of the baryon density of the Universe. The simulations also suggest that more than half of the baryons in the Universe are in the form of diffuse gas outside halos ($\gtrsim 5$ times the virial radius) which is not hot or dense enough to produce a significant tSZ signal or be observed by X-ray experiments. Finally, we show that future high-resolution tSZ-lensing cross-correlation observations will serve as a powerful tool for discriminating between different galactic feedback models.

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

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

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

The XXL survey is the largest survey carried out by XMM-Newton. Covering an area of 50deg$^2$, the survey contains $\sim450$ galaxy clusters out to a redshift $\sim$2 and to an X-ray flux limit of $\sim5\times10^{-15}erg\,s^{-1}cm^{-2}$. This paper is part of the first release of XXL results focussed on the bright cluster sample. We investigate the scaling relation between weak-lensing mass and X-ray temperature for the brightest clusters in XXL. The scaling relation is used to estimate the mass of all 100 clusters in XXL-100-GC. Based on a subsample of 38 objects that lie within the intersection of the northern XXL field and the publicly available CFHTLenS catalog, we derive the $M_{WL}$ of each system with careful considerations of the systematics. The clusters lie at $0.1<z<0.6$ and span a range of $ T\simeq1-5keV$. We combine our sample with 58 clusters from the literature, increasing the range out to 10keV. To date, this is the largest sample of clusters with $M_{WL}$ measurements that has been used to study the mass-temperature relation. The fit ($M\propto T^b$) to the XXL clusters returns a slope $b=1.78^{+0.37}_{-0.32}$ and intrinsic scatter $\sigma_{\ln M|T}\simeq0.53$; the scatter is dominated by disturbed clusters. The fit to the combined sample of 96 clusters is in tension with self-similarity, $b=1.67\pm0.12$ and $\sigma_{\ln M|T}\simeq0.41$. Overall our results demonstrate the feasibility of ground-based weak-lensing scaling relation studies down to cool systems of $\sim1keV$ temperature and highlight that the current data and samples are a limit to our statistical precision. As such we are unable to determine whether the validity of hydrostatic equilibrium is a function of halo mass. An enlarged sample of cool systems, deeper weak-lensing data, and robust modelling of the selection function will help to explore these issues further.

Traditionally, galaxy clusters have been expected to retain all the material accreted since their formation epoch. For this reason, their matter content should be representative of the Universe as a whole, and thus their baryon fraction should be close to the Universal baryon fraction. We make use of the sample of the 100 brightest galaxy clusters discovered in the XXL Survey to investigate the fraction of baryons in the form of hot gas and stars in the cluster population. We measure the gas masses of the detected halos and use a mass--temperature relation directly calibrated using weak-lensing measurements for a subset of XXL clusters to estimate the halo mass. We find that the weak-lensing calibrated gas fraction of XXL-100-GC clusters is substantially lower than was found in previous studies using hydrostatic masses. Our best-fit relation between gas fraction and mass reads $f_{\rm gas,500}=0.055_{-0.006}^{+0.007}\left(M_{\rm 500}/10^{14}M_\odot\right)^{0.21_{-0.10}^{+0.11}}$. The baryon budget of galaxy clusters therefore falls short of the Universal baryon fraction by about a factor of two at $r_{\rm 500}$. Our measurements require a hydrostatic bias $1-b=M_X/M_{\rm WL}=0.72_{-0.07}^{+0.08}$ to match the gas fraction obtained using lensing and hydrostatic equilibrium. Comparing our gas fraction measurements with the expectations from numerical simulations, our results favour an extreme feedback scheme in which a significant fraction of the baryons are expelled from the cores of halos. This model is, however, in contrast with the thermodynamical properties of observed halos, which might suggest that weak-lensing masses are overestimated. We note that a mass bias $1-b=0.58$ as required to reconcile Planck CMB and cluster counts should translate into an even lower baryon fraction, which poses a major challenge to our current understanding of galaxy clusters. [Abridged]

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 test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $\beta_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $\beta_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $\beta_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5\sigma$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts.

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7<log_{10}(M_{200}/[h^{-1}M_\odot])<15$) and a large range of separations ($-1<log_{10}(r/[h^{-1}Mpc])< 2$). We focus on two classes of alignments: the orientations of galaxies with respect to either the directions to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is a strongly decreasing function of the distance between galaxies. The orientation-direction alignment can remain significant up to ~100 Mpc, for galaxies hosted by the most massive haloes in our simulations. Galaxies hosted by more massive subhaloes show stronger alignment. At a fixed halo mass, more aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites is anisotropic and significantly aligned with the major axis of the main host halo. The major axis of satellite galaxies, when all stars are considered, are preferentially aligned towards the centre of the main host halo. The predicted projected direction-orientation alignment, $\epsilon_{g+}(r_{p})$, is in broad agreement with recent observations when only stars within the typical observable extent of a galaxy are used to define galaxy orientations. We find that the orientation-orientation alignment is weaker than the orientation-direction alignment on all scales. Overall, the strength of galaxy alignments depends strongly on the subset of stars that are used to measure the orientations of galaxies and it is always weaker than the alignment of the dark matter haloes. Thus, alignment models that use halo orientation as a direct proxy for galaxy orientation will overestimate the impact of intrinsic alignments on weak lensing analyses.

We present results from a subset of simulations from the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) suite in which the formulation of the hydrodynamics scheme is varied. We compare simulations that use the same subgrid models without re-calibration of the parameters but employing the standard GADGET flavour of smoothed particle hydrodynamics (SPH) instead of the more recent state-of-the-art ANARCHY formulation of SPH that was used in the fiducial EAGLE runs. We find that the properties of most galaxies, including their masses and sizes, are not significantly affected by the details of the hydrodynamics solver. However, the star formation rates of the most massive objects are affected by the lack of phase mixing due to spurious surface tension in the simulation using standard SPH. This affects the efficiency with which AGN activity can quench star formation in these galaxies and it also leads to differences in the intragroup medium that affect the X-ray emission from these objects. The differences that can be attributed to the hydrodynamics solver are, however, likely to be less important at lower resolution. We also find that the use of a time step limiter is important for achieving the feedback efficiency required to match observations of the low-mass end of the galaxy stellar mass function.

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with the local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations.

We use the "Evolution and assembly of galaxies and their environments" (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $\rho_*(r) \propto r^{-3}$. Stars dominate the mass density for $r<10~\rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($r\gtrsim3.0~\rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy.

The evolution of the metal content of galaxies and its relations to other global properties [such as total stellar mass (M*), circular velocity, star formation rate (SFR), halo mass, etc.] provides important constraints on models of galaxy formation. Here we examine the evolution of metallicity scaling relations of simulated galaxies in the Galaxies-Intergalactic Medium Interaction Calculation suite of cosmological simulations. We make comparisons to observations of the correlation of gas-phase abundances with M* (the mass-metallicity relation, MZR), as well as with both M* and SFR or gas mass fraction (the so-called 3D fundamental metallicity relations, FMRs). The simulated galaxies follow the observed local MZR and FMRs over an order of magnitude in M*, but overpredict the metallicity of massive galaxies (log M* > 10.5), plausibly due to inefficient feedback in this regime. We discuss the origin of the MZR and FMRs in the context of galactic outflows and gas accretion. We examine the evolution of mass-metallicity relations defined using different elements that probe the three enrichment channels (SNII, SNIa, and AGB stars). Relations based on elements produced mainly by SNII evolve weakly, whereas those based on elements produced preferentially in SNIa/AGB exhibit stronger evolution, due to the longer timescales associated with these channels. Finally, we compare the relations of central and satellite galaxies, finding systematically higher metallicities for satellites, as observed. We show this is due to the removal of the metal poor gas reservoir that normally surrounds galaxies and acts to dilute their gas-phase metallicity (via cooling/accretion onto the disk), but is lost due to ram pressure stripping for satellites.

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $\Lambda$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass.

The thermal Sunyaev-Zel'dovich (tSZ) effect offers a means of probing the hot gas in and around massive galaxies and galaxy groups and clusters, which is thought to constitute a large fraction of the baryon content of the Universe. The Planck collaboration recently performed a stacking analysis of a large sample of `locally brightest galaxies' (LBGs) and, surprisingly, inferred an approximately self-similar relation between the tSZ flux and halo mass. At face value, this implies that the hot gas mass fraction is independent of halo mass, a result which is in apparent conflict with resolved X-ray observations. We test the robustness of the inferred trend using synthetic tSZ maps generated from cosmological hydrodynamical simulations and using the same tools and assumptions applied in the Planck study. We show that, while the detection and the estimate of the `total' flux (within $5 r_{500}$) is reasonably robust, the inferred flux originating from within $r_{500}$ (i.e. the limiting radius to which X-ray observations typically probe) is highly sensitive to the assumed pressure distribution of the gas. Using our most realistic simulations with AGN feedback, that reproduce a wide variety of X-ray and optical properties of groups and clusters, we estimate that the derived tSZ flux within $r_{500}$ is biased high by up to to an order of magnitude for haloes with masses $M_{500} \sim 10^{13}$ M$_{\odot}$. Moreover, we show that the AGN simulations are consistent with the total tSZ flux-mass relation observed with Planck, whereas a self-similar model is ruled out.

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation.

We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EAGLE simulation suite using the GALAXEV population synthesis models. We take into account obscuration by dust in birth clouds and diffuse ISM using a two-component screen model, following the prescription of Charlot and Fall. We compare models in which the dust optical depth is constant to models where it depends on gas metallicity, gas fraction and orientation. The colours of EAGLE galaxies for the more sophisticated models are in broad agreement with those of observed galaxies. In particular, EAGLE produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, with approximately the correct fraction of galaxies in each population and with g-r colours within 0.1 magnitudes of those observed. Luminosity functions from UV to NIR wavelengths differ from observations at a level comparable to systematic shifts resulting from a choice between Petrosian and Kron photometric apertures. Despite the generally good agreement there are clear discrepancies with observations. The blue cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, consistent with the modest underestimate of the passive fraction in massive EAGLE galaxies. There is also a moderate excess of bright blue galaxies compared to observations. The overall level of agreement with the observed colour distribution suggests that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that are comparable to those of observed galaxies.

We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

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.

Star formation is observed to be suppressed in group and cluster galaxies compared to the field. To gain insight into the quenching process, we have analysed ~2000 galaxies formed in the GIMIC suite of cosmological hydrodynamical simulations. The time of quenching varies from ~2 Gyr before accretion (first crossing of r200,c) to >4 Gyr after, depending on satellite and host mass. Once begun, quenching is rapid (>~ 500 Myr) in low-mass galaxies (M* < 10^10 M_Sun), but significantly more protracted for more massive satellites. The simulations predict a substantial role of outflows driven by ram pressure -- but not tidal forces -- in removing the star-forming interstellar matter (ISM) from satellite galaxies, especially dwarfs (M* ~ 10^9 M_Sun) where they account for nearly two thirds of ISM loss in both groups and clusters. Immediately before quenching is complete, this fraction rises to ~80% even for Milky Way analogues (M* ~ 10^10.5 M_Sun) in groups (M_host ~ 10^13.5 M_Sun). We show that (i) ISM stripping was significantly more effective at early times than at z = 0; (ii) approximately half the gas is stripped from `galactic fountains' and half directly from the star forming disk; (iii) galaxies undergoing stripping experience ram pressure up to ~100 times the average at a given group/cluster-centric radius, because they are preferentially located in overdense ICM regions. Remarkably, stripping causes at most half the loss of the extended gas haloes surrounding our simulated satellites. These results contrast sharply with the current picture of strangulation -- removal of the ISM through star formation after stripping of the hot halo -- being the dominant mechanism quenching group and cluster satellites.

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 use cosmological hydrodynamical simulations to investigate how the inclusion of physical processes relevant to galaxy formation (star formation, metal-line cooling, stellar winds, supernovae and feedback from Active Galactic Nuclei, AGN) change the properties of haloes, over four orders of magnitude in mass. We find that gas expulsion and the associated dark matter (DM) expansion induced by supernova-driven winds are important for haloes with masses M200 < 10^13 Msun, lowering their masses by up to 20% relative to a DM-only model. AGN feedback, which is required to prevent overcooling, has a significant impact on halo masses all the way up to cluster scales (M200 ~ 10^15 Msun). Baryonic physics changes the total mass profiles of haloes out to several times the virial radius, a modification that cannot be captured by a change in the halo concentration. The decrease in the total halo mass causes a decrease in the halo mass function of about 20%. This effect can have important consequences for abundance matching technique as well as for most semi-analytic models of galaxy formation. We provide analytic fitting formulae, derived from simulations that reproduce the observed baryon fractions, to correct halo masses and mass functions from DM-only simulations. The effect of baryonic physics (AGN feedback in particular) on cluster number counts is about as large as changing the cosmology from WMAP7 to Planck, even when a moderately high mass limit of M500 ~ 10^14 Msun is adopted. Thus, for precision cosmology the effects of baryons must be accounted for.

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

We present a new suite of large-volume cosmological hydrodynamical simulations called cosmo-OWLS. They form an extension to the OverWhelmingly Large Simulations (OWLS) project, and have been designed to help improve our understanding of cluster astrophysics and non-linear structure formation, which are now the limiting systematic errors when using clusters as cosmological probes. Starting from identical initial conditions in either the Planck or WMAP7 cosmologies, we systematically vary the most important `sub-grid' physics, including feedback from supernovae and active galactic nuclei (AGN). We compare the properties of the simulated galaxy groups and clusters to a wide range of observational data, such as X-ray luminosity and temperature, gas mass fractions, entropy and density profiles, Sunyaev-Zel'dovich flux, I-band mass-to-light ratio, dominance of the brightest cluster galaxy, and central massive black hole (BH) masses, by producing synthetic observations and mimicking observational analysis techniques. These comparisons demonstrate that some AGN feedback models can produce a realistic population of galaxy groups and clusters, broadly reproducing both the median trend and, for the first time, the scatter in physical properties over approximately two decades in mass ($10^{13} \lesssim M_{500} \lesssim 10^{15}~\textrm{M}_{\odot}$) and 1.5 decades in radius ($0.05 \lesssim r/r_{500} \lesssim 1.5$). However, in other models, the AGN feedback is too violent (even though they reproduce the observed BH scaling relations), implying calibration of the models is required. The production of realistic populations of simulated groups and clusters, as well as models that bracket the observations, opens the door to the creation of synthetic surveys for assisting the astrophysical and cosmological interpretation of cluster surveys, as well as quantifying the impact of selection effects.

(Abridged) The amplitude of the thermal Sunyaev Zel'dovich effect (tSZ) power spectrum is extremely sensitive to the abundance of galaxy clusters and therefore to fundamental cosmological parameters that control their growth, such as sigma_8 and Omega_m. Here we explore the sensitivity of the tSZ power spectrum to important non-gravitational ('sub-grid') physics by employing the cosmo-OWLS suite of large-volume cosmological hydrodynamical simulations, run in both the Planck and WMAP7 best-fit cosmologies. On intermediate and small angular scales (ell > ~1000, or theta < ~10 arcmin), accessible with the South Pole Telescope and the Atacama Cosmology Telescope, the predicted tSZ power spectrum is highly model dependent, with AGN feedback having a particularly large effect. However, at large scales, observable with the Planck telescope, the effects of sub-grid physics are minor. Comparing the simulations with observations, we find a significant amplitude offset on all measured angular scales (including large scales), if the Planck best-fit cosmology is assumed by the simulations. This is shown to be a generic result for all current tSZ models. By contrast, if the WMAP7 cosmology is adopted, there is full consistency with the Planck power spectrum measurements on large scales and agreement at the 2 sigma level with the SPT/ACT measurements at intermediate scales for our fiducial AGN model, which Le Brun et al. (2014) have shown reproduces the 'resolved' properties of the local cluster population remarkably well. These findings strongly suggest that there are significantly fewer massive galaxy clusters than expected for the Planck best-fit cosmology, which is consistent with recent measurements of the tSZ number counts. Our findings therefore pose a significant challenge to the cosmological parameter values preferred (and/or the model adopted) by the Planck primary CMB analyses.

The observed clustering of galaxies and the cross-correlation of galaxies and mass provide important constraints on both cosmology and models of galaxy formation. Even though the dissipation and feedback processes associated with galaxy formation are thought to affect the distribution of matter, essentially all models used to predict clustering data are based on collisionless simulations. Here, we use large hydrodynamical simulations to investigate how galaxy formation affects the autocorrelation functions of galaxies and subhaloes, as well as their cross-correlation with matter. We show that the changes due to the inclusion of baryons are not limited to small scales and are even present in samples selected by subhalo mass. Samples selected by subhalo mass cluster ~10% more strongly in a baryonic run on scales r > 1Mpc/h, and this difference increases for smaller separations. While the inclusion of baryons boosts the clustering at fixed subhalo mass on all scales, the sign of the effect on the cross-correlation of subhaloes with matter can vary with radius. We show that the large-scale effects are due to the change in subhalo mass caused by the strong feedback associated with galaxy formation and may therefore not affect samples selected by number density. However, on scales r < r_vir significant differences remain after accounting for the change in subhalo mass. We conclude that predictions for galaxy-galaxy and galaxy-mass clustering from models based on collisionless simulations will have errors greater than 10% on sub-Mpc scales, unless the simulation results are modified to correctly account for the effects of baryons on the distributions of mass and satellites.

We employ a large sample of 20171 optically-selected groups and clusters at 0.15 < z < 0.4 in the SDSS to investigate how the stacked stellar mass fraction varies across a wide range of total mass, $M_{500}$. Our study improves upon previous observational studies in a number of important ways, including having a much larger sample size, an explicit inclusion of the intracluster light (ICL) component, and a thorough examination of the accuracy of our total mass estimates via comparisons to simulations and weak lensing observations. We find that the stellar mass fraction depends only weakly on total mass and that the contribution of ICL to the total stellar mass fraction is significant (typically 20-40 per cent). Both of these findings are in excellent accordance with the predictions of cosmological simulations. Under the assumption of a Chabrier (Salpeter) IMF, the derived star formation efficiency ($f_{star}$/$f_{b}$, where $f_b=\Omega_b/\Omega_m$) is relatively low at 8 per cent (14 per cent) and is consistent with the global star formation efficiency of semi-analytic models that reproduce the galaxy stellar mass function. When our measured stellar mass fractions are combined with the observed relation between hot gas mass fraction and total mass from X-ray observations, our results imply that galaxy groups have significantly lower baryon fractions than massive clusters. Ejection of gas due to energetic AGN feedback (most likely at high redshift) provides a plausible mechanism for explaining the trends we observe.

We investigate the relationship between the halo mass, M_200, and concentration, c, for a sample of 26 group- and cluster-scale strong gravitational lenses. In contrast with previous results, we find that these systems are only ~ 0.1 dex more over-concentrated than similar-mass halos from dark matter simulations; the concentration of a halo with M_200 = 10^14 M_sun is log c = 0.78\pm0.05, while simulations of halos with this mass at similar redshifts (z ~ 0.4) predict log c ~ 0.56 - 0.71. We also find that we are unable to make informative inference on the slope of the M_200-c relation in spite of our large sample size; we note that the steep slopes found in previous studies tend to follow the slope in the covariance between M_200 and c, indicating that these results may be measuring the scatter in the data rather than the intrinsic signal. Furthermore, we conclude that our inability to constrain the M_200-c slope is due to a limited range of halo masses, as determined by explicitly modelling our halo mass distribution, and we suggest that other studies may be producing biased results by using an incorrect distribution for their halo masses.

Mechanical feedback via Active Galactic Nuclei (AGN) jets in the centres of galaxy groups and clusters is a crucial ingredient in current models of galaxy formation and cluster evolution. Jet feedback is believed to regulate gas cooling and thus star formation in the most massive galaxies, but a robust physical understanding of this feedback mode is currently lacking. The large collecting area, excellent spectral resolution and high spatial resolution of Athena+ will provide the breakthrough diagnostic ability necessary to develop this understanding, via: (1) the first kinematic measurements on relevant spatial scales of the hot gas in galaxy, group and cluster haloes as it absorbs the impact of AGN jets, and (2) vastly improved ability to map thermodynamic conditions on scales well-matched to the jets, lobes and gas disturbances produced by them. Athena+ will therefore determine for the first time how jet energy is dissipated and distributed in group and cluster gas, and how a feedback loop operates in group/cluster cores to regulate gas cooling and AGN fuelling. Athena+ will also establish firmly the cumulative impact of powerful radio galaxies on the evolution of baryons from the epoch of group/cluster formation to the present day.

Models of galaxy formation in a CDM universe predict that massive galaxies are surrounded by a hot, quasi-hydrostatic circumgalactic corona of slowly cooling gas, predominantly accreted from the IGM. This prediction is borne out by the cosmological hydrodynamical simulations of Crain et al., which reproduce scaling relations between the X-ray and optical properties of nearby disc galaxies. Such coronae are metal poor, but observations of the X-ray emitting circumgalactic medium (CGM) of local galaxies typically indicate enrichment to near-solar iron abundance, potentially signalling a shortcoming in galaxy formation models. We show here that, while the hot CGM of galaxies formed in the simulations is metal poor in a mass-weighted sense, its X-ray luminosity-weighted metallicity is often close to solar. This bias arises because the soft X-ray emissivity of a typical 0.1 keV corona is dominated by collisionally-excited metal ions that are synthesised in stars and recycled into the hot CGM. We find that these metals are ejected primarily by stars that form in-situ to the main progenitor of the galaxy, rather than in satellites or external galaxies. The enrichment of the hot CGM therefore proceeds in an inside-out fashion throughout the assembly of the galaxy: metals are transported from the central galaxy by SNe-driven winds and convection over several gigayears, establishing a strong negative radial metallicity gradient. Whilst metal ions synthesised by stars are necessary to produce the X-ray emissivity that enables the hot CGM of isolated galaxies to be detected, the electrons that collisionally excite them are equally important. Since our simulations indicate that the electron density of hot coronae is dominated by the metal-poor gas accreted from the IGM, we infer that the hot CGM observed via X-ray emission is the outcome of both hierarchical accretion and stellar recycling.

In the local Universe, galaxies in groups and clusters contain less gas and are less likely to be forming stars than their field counterparts. This effect is not limited to the central group/cluster regions, but is shown by recent observations to persist out to several virial radii. To gain insight into the extent and cause of this large-scale environmental influence, we use a suite of high-resolution cosmological hydrodynamic simulations to analyse galaxies around simulated groups and clusters of a wide range of mass (log M/M_sun = [13.0, 15.2]). In qualitative agreement with the observations, we find a systematic depletion of both hot and cold gas and a decline in the star forming fraction of galaxies as far out as ~ 5 r200 from the host centre. While a substantial fraction of these galaxies are on highly elliptical orbits and are not infalling for the first time (~ 50 per cent at 2 r200, independent of host mass) or are affected by `pre-processing' (less than 10 per cent of galaxies around groups, increasing to ~ 50 per cent around a massive cluster), even a combination of these indirect mechanisms does not fully account for the environmental influence, particularly in the case of the hot gas content. Direct ram pressure interaction with an extended gas `halo' surrounding groups and clusters is shown to be sufficiently strong to strip the hot gas atmospheres of infalling galaxies out to ~ 5 r200. We show that this influence is highly anisotropic, with ram pressure along filaments enhanced by up to a factor of 100 despite significant co-flow of gas and galaxies.

We present adaptive optics assisted, spatially resolved spectroscopy of a sample of nine H-alpha-selected galaxies at z=0.84--2.23 drawn from the HiZELS narrow-band survey. These galaxies have star-formation rates of 1-27Mo/yr and are therefore representative of the typical high-redshift star-forming population. Our ~kpc-scale resolution observations show that approximately half of the sample have dynamics suggesting that the ionised gas is in large, rotating disks. We model their velocity fields to infer the inclination-corrected, asymptotic rotational velocities. We use the absolute B-band magnitudes and stellar masses to investigate the evolution of the B-band and stellar mass Tully-Fisher relationships. By combining our sample with a number of similar measurements from the literature, we show that, at fixed circular velocity, the stellar mass of star-forming galaxies has increased by a factor 2.5 between z=2 and z=0, whilst the rest-frame B-band luminosity has decreased by a factor ~6 over the same period. Together, these demonstrate a change in mass-to-light ratio in the B-band of Delta(M/L_B)/(M/L_B)_(z=0) \sim 3.5 between z=1.5 and z=0, with most of the evolution occurring below z=1. We also use the spatial variation of [NII]/Halpha to show that the metallicity of the ionised gas in these galaxies declines monotonically with galacto-centric radius, with an average Delta(log O/H)/DeltaR=-0.027+/-0.005dex/kpc. This gradient is consistent with predictions for high-redshift disk galaxies from cosmologically based hydrodynamic simulations.

We examine the rotation rates, sizes, and star formation (SF) efficiencies of a representative population of simulated disc galaxies extracted from the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamic simulations. These simulations include efficient, but energetically feasible supernova feedback, but have not been tuned in any way to produce 'realistic' disc galaxies. Yet, they generate a large number of discs, without requiring extremely high resolution. Over the wide galaxy stellar mass range, 9.0 < log10[Mstar (Msun)] < 10.5, the simulations reproduce the observed Tully-Fisher relation, the rotation curves of disc galaxies in bins of stellar mass, the mass-size relation of disc galaxies, the optical rotation to virial circular velocity ratio (Vopt/Vvir), and the SF efficiencies of disc galaxies as inferred from stacked weak lensing and stacked satellite kinematics observations. They also reproduce the specific star formation rates of ~L* galaxies but predict too low levels of star formation for low-mass galaxies, which is plausibly due to the finite resolution of the simulations. At higher stellar masses, log10[Mstar (\Msun)] > 10.6, the simulated galaxies are too concentrated and have too high SF efficiencies. We conjecture that this shortcoming reflects the neglect of feedback from accreting supermassive black holes in these simulations. We conclude that it is possible to generate a representative population of disc galaxies that reproduces many of the observed trends of local disc galaxies using standard numerical hydrodynamic techniques and a plausible implementation of the "subgrid" astrophysical processes thought to be relevant to galaxy formation.

We use matter power spectra from cosmological hydrodynamic simulations to quantify the effect of baryon physics on the weak gravitational lensing shear signal. The simulations consider a number of processes, such as radiative cooling, star formation, supernovae and feedback from active galactic nuclei (AGN). Van Daalen et al. (2011) used the same simulations to show that baryon physics, in particular the strong feedback that is required to solve the overcooling problem, modifies the matter power spectrum on scales relevant for cosmological weak lensing studies. As a result, the use of power spectra from dark matter simulations can lead to significant biases in the inferred cosmological parameters. We show that the typical biases are much larger than the precision with which future missions aim to constrain the dark energy equation of state, w_0. For instance, the simulation with AGN feedback, which reproduces X-ray and optical properties of groups of galaxies, gives rise to a ~40% bias in w_0. We demonstrate that the modification of the power spectrum is dominated by groups and clusters of galaxies, the effect of which can be modelled. We consider an approach based on the popular halo model and show that simple modifications can capture the main features of baryonic feedback. Despite its simplicity, we find that our model, when calibrated on the simulations, is able to reduce the bias in w_0 to a level comparable to the size of the statistical uncertainties for a Euclid-like mission. While observations of the gas and stellar fractions as a function of halo mass can be used to calibrate the model, hydrodynamic simulations will likely still be needed to extend the observed scaling relations down to halo masses of 10 ^12 M_sun/h.

Ram pressure stripping of the hot gas that surrounds normal galaxies as they fall into groups and clusters (also referred to as `strangulation' or `starvation') is generally thought to shut down star formation on a time scale of a few Gyr. However, it has recently been suggested, on the basis of X-ray-optical scaling relations of galaxies in the field and the group/cluster environment, that confinement pressure by the intra-cluster medium can actually lead to an increase in the mass of hot gas surrounding these galaxies. We investigate the competition between pressure confinement and ram pressure stripping for satellite galaxies in orbit about galaxy groups and clusters using simple analytic models and detailed cosmological hydrodynamic simulations. It is found that, independent of host mass, ram pressure is generally dominant over confinement pressure --- only ~16 % of galaxies find themselves in the reverse situation. Furthermore, these galaxies have, on average, less hot gas than ram-pressure dominated ones, contrary to simple expectations. This is explained by the fact that the small number of galaxies which are confinement dominated are typically at first or second apocentre and have therefore already been maximally affected by ram pressure stripping around first pericentre. Our results are shown to be insensitive to host halo mass; we argue that the same is true for uncertain sub-grid processes, such as feedback.

(Abridged) We quantify the bias and scatter in galaxy cluster masses and concentrations derived from an idealised mock weak gravitational lensing (WL) survey, and their effect on the cluster mass-concentration relation. For this, we simulate WL distortions on a population of background galaxies due to a large (~3000) sample of galaxy cluster haloes extracted from the Millennium Simulation at z~0.2. This study takes into account the influence of shape noise, cluster substructure and asphericity as well as correlated large-scale structure, but not uncorrelated large-scale structure along the line of sight and observational effects. We find a small, but non-negligble, negative median bias in both mass and concentration at a level of ~5%, the exact value depending both on cluster mass and radial survey range. Both the mass and concentration derived from WL show considerable scatter about their true values. This scatter has, even for the highest mass clusters of M200 > 10^14.8 M_sun, a level of ~30% and ~20% for concentration and mass respectively and increases strongly with decreasing cluster mass. For a typical survey analysing 30 galaxies per arcmin^2 over a radial range from 30" to 15' from the cluster centre, the derived M200-c relation has a slope and normalisation too low compared to the underlying true (3D) relation by ~40% and ~15% respectively. The scatter and bias in mass are shown to reflect a departure at large radii of the true WL shear/matter distribution of the simulated clusters from the NFW profile adopted in modelling the mock observations. Orientation of the triaxial cluster haloes dominates the concentration scatter (except at low masses, where galaxy shape noise becomes dominant), while the bias in c is mostly due to substructure within the virial radius.

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 new catalogue of 55,121 groups and clusters centred on Luminous Red Galaxies from SDSS DR7 in the redshift range 0.15<z<0.4. We provide halo mass estimates for each of these groups derived from a calibration between the optical richness of bright galaxies (M_r<-20.5) within 1 Mpc, and X-ray-derived mass for a small subset of 129 groups and clusters with X-ray measurements. We derive the mean (stacked) surface number density profiles of galaxies as a function of total halo mass in different mass bins. We find that derived profiles can be well-described by a projected NFW profile with a concentration parameter (<c>~2.6) that is approximately a factor of two lower than that of the dark matter (as predicted by N-body cosmological simulations) and nearly independent of halo mass. Interestingly, in spite of the difference in shape between the galaxy and dark matter radial distributions, both exhibit a high degree of self-similarity. A self-consistent comparison to several recent semi-analytic models of galaxy formation indicates that: (1) beyond ~0.3 r_500 current models are able to reproduce both the shape and normalisation of the satellite profiles; and (2) within ~0.3 r_500 the predicted profiles are sensitive to the details of the satellite-BCG merger timescale calculation. The former is a direct result of the models being tuned to match the global galaxy luminosity function combined with the assumption that the satellite galaxies do not suffer significant tidal stripping, even though their surrounding dark matter haloes can be removed through this process. Combining our results with measurements of the intracluster light should provide a way to inform theoretical models on the efficacy of the tidal stripping and merging processes.

We estimate the masses of elliptical galaxies out to five effective radii using planetary nebulae and globular clusters as tracers. A sample of 15 elliptical galaxies with a broad variation in mass is compiled from the literature. A distribution function-maximum likelihood analysis is used to estimate the overall potential slope, normalisation and velocity anisotropy of the tracers. We assume power-law profiles for the potential and tracer density and a constant velocity anisotropy. The derived potential power-law indices lie in between the isothermal and Keplerian regime and vary with mass: there is tentative evidence that the less massive galaxies have steeper potential profiles than the more massive galaxies. We use stellar mass-to-light ratios appropriate for either a Chabrier/KTG (Kroupa, Tout & Gilmore) or Salpeter initial mass function to disentangle the stellar and dark matter components. The fraction of dark matter within five effective radii increases with mass, in agreement with several other studies. We employ simple models to show that a combination of star formation efficiency and baryon extent are able to account for this trend. These models are in good agreement with both our measurements out to five effective radii and recent SLACS measurements within one effective radii when a universal Chabrier/KTG initial mass function is adopted.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r < ~30 kpc. In the present study we show that this component was primarily formed in a proto-disc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, that is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axis ratio of ~0.6 at radii of < ~30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ vs. accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a 'dual stellar halo'.

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the formation of stellar spheroids of Milky Way-mass disc galaxies. The simulations contain accurate treatments of metal-dependent radiative cooling, star formation, supernova feedback, and chemodynamics, and the large volumes that have been simulated yield an unprecedentedly large sample of ~400 simulated L_* disc galaxies. The simulated galaxies are surrounded by low-mass, low-surface brightness stellar haloes that extend out to ~100 kpc and beyond. The diffuse stellar distributions bear a remarkable resemblance to those observed around the Milky Way, M31 and other nearby galaxies, in terms of mass density, surface brightness, and metallicity profiles. We show that in situ star formation typically dominates the stellar spheroids by mass at radii of r < 30 kpc, whereas accretion of stars dominates at larger radii and this change in origin induces a change in slope of the surface brightness and metallicity profiles, which is also present in the observational data. The system-to-system scatter in the in situ mass fractions of the spheroid, however, is large and spans over a factor of 4. Consequently, there is a large degree of scatter in the shape and normalisation of the spheroid density profile within r < 30 kpc (e.g., when fit by a spherical powerlaw profile the indices range from -2.6 to -3.4). We show that the in situ mass fraction of the spheroid is linked to the formation epoch of the system. Dynamically older systems have, on average, larger contributions from in situ star formation, although there is significant system-to-system scatter in this relationship. Thus, in situ star formation likely represents the solution to the longstanding failure of pure accretion-based models to reproduce the observed properties of the inner spheroid.

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.

Galaxy groups are not scaled down versions of massive galaxy clusters - the hot gas in groups (known as the intragroup medium, IGrM for short) is, on average, less dense than the intracluster medium, implying that one or more non-gravitational processes (e.g., radiative cooling, star formation, and/or feedback) has had a relatively larger effect on groups. In the present study, we compare a number of cosmological hydrodynamic simulations that form part of the OverWhelmingly Large Simulations project to isolate and quantify the effects of cooling and feedback from supernovae (SNe) and active galactic nuclei (AGN) on the gas. This is achieved by comparing Lagrangian thermal histories of the gas in the different runs, which were all started from identical initial conditions. While radiative cooling, star formation, and SN feedback are all necessary ingredients, only runs that also include AGN feedback are able to successfully reproduce the optical and X-ray properties of groups and low-mass clusters. We isolate how, when, and exactly what gas is heated by AGN. Interestingly, we find that the gas that constitutes the present-day IGrM is that which was not strongly heated by AGN. Instead, the low median density/high median entropy of the gas in present-day groups is achieved by the ejection of lower entropy gas from low-mass progenitor galaxies at high redshift (primarily 2 < z < 4). This corresponds to the epoch when supermassive black holes accreted most of their mass, typically at a rate that is close to the Eddington limit (i.e., when the black holes are in a `quasar mode').

The existence of hot, X-ray luminous gaseous coronae surrounding present day L* galaxies is a generic prediction of galaxy formation theory in the cold dark matter cosmogony. While extended X-ray emission has been known to exist around elliptical galaxies for a long time, diffuse extra-planar emission has only recently been detected around disc galaxies. We compile samples of elliptical and disc galaxies that have Chandra and XMM-Newton measurements, and compare the scaling of the coronal X-ray luminosity (L_X) with both the K-band luminosity (L_K) and the coronal X-ray temperature (T_X). The X-ray flux measurements are corrected for non-thermal point source contamination by spatial excision and spectral subtraction for resolved and unresolved sources respectively. We find that the properties of the extended X-ray emission from galaxies of different morphological types are similar: for both elliptical and disc galaxies, the L_X - L_K and L_X - T_X relations have similar slope, normalisation and scatter. The observed universality of coronal X-ray properties suggests that the bulk of this emission originates from gas that has been accreted, shock-heated and compressed during the assembly of the galaxy and that outflows triggered by stellar processes make only a minor contribution to the X-ray emission. This reservoir of cooling gas is a potential source of morphological transformation; it provides a fresh supply of material for discs to grow around galaxies of all morphological types.

The existence of X-ray luminous gaseous coronae around massive disc galaxies is a long-standing prediction of galaxy formation theory in the cold dark matter cosmogony. This prediction has garnered little observational support, with non-detections commonplace and detections for only a relatively small number of galaxies which are much less luminous than expected. We investigate the coronal properties of a large sample of bright, disc-dominated galaxies extracted from the GIMIC suite of cosmological hydrodynamic simulations recently presented by Crain et al. Remarkably, the simulations reproduce the observed scalings of X-ray luminosity with K-band luminosity and star formation rate and, when account is taken of the density structure of the halo, with disc rotation velocity as well. Most of the star formation in the simulated galaxies (which have realistic stellar mass fractions) is fuelled by gas cooling from a quasi-hydrostatic hot corona. However, these coronae are more diffuse, and of a lower luminosity, than predicted by the analytic models of White & Frenk because of a substantial increase in entropy at z ~ 1-3. Both the removal of low entropy gas by star formation and energy injection from supernovae contribute to this increase in entropy, but the latter is dominant for halo masses M_200 <~ 10^(12.5) Msun. Only a small fraction of the mass of the hot gas is outflowing as a wind but, because of its high density and metallicity, it contributes disproportionally to the X-ray emission. The bulk of the X-ray emission, however, comes from the diffuse quasi-hydrostatic corona which supplies the fuel for ongoing star formation in discs today. Future deep X-ray observations with high spectral resolution (e.g. with NeXT/ASTRO-H or IXO) should be able to map the velocity structure of the hot gas and test this fundamental prediction of current galaxy formation theory.

[Abridged] The relatively recent insight that energy input from supermassive black holes (BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies motivates us to examine its effects on the scale of galaxy groups. At present, groups contain most of the galaxies and a significant fraction of the overall baryon content of the universe. To explore the effects of BH feedback on groups, we analyse two high resolution cosmological hydro simulations from the OverWhelmingly Large Simulations project. While both include galactic winds driven by supernovae, only one includes feedback from BHs. We compare the properties of the simulated groups to a wide range of observational data, including hot gas radial profiles and gas mass fractions (fgas), luminosity-mass-temperature (L-M-T) scaling relations, K-band luminosity of the group and its central brightest galaxy (CBG), SFRs and ages of the CBG, and gas/stellar metallicities. Both runs yield entropy profiles similar to the data, while the run without AGN feedback yields highly peaked temperature profiles, in discord with the observations. Energy input from BHs significantly reduces fgas for groups with masses less than ~10^14 Msun, yielding fgas-T and L-T relations that are in agreement with the data. The run without AGN feedback suffers from the well known overcooling problem; the resulting K-band luminosities are much larger than observed. By contrast, the run that includes BH feedback yields K-band luminosities and CBG SFRs and ages in agreement with current estimates. Both runs yield very similar gas-phase metallicities that match X-ray data, but they predict very different stellar metallicities. Based on the above, galaxy groups provide a compelling case that BH feedback is a crucial ingredient in the formation of massive galaxies.

We investigate the impact of AGN feedback on the gravitational lensing properties of a sample of galaxy clusters with masses in the range 10^14 - 10^15 Msol, using state-of-the-art simulations. Adopting a ray-tracing algorithm, we compute the cross-section of giant arcs from clusters simulated with dark-matter only physics (DM), dark matter plus gas with cooling and star formation (CSF), and dark matter plus gas with cooling, star formation and AGN feedback (CSFBH). Once AGN feedback is included, baryonic physics boosts the strong lensing cross-section by much less than previously estimated using clusters simulated with only cooling and star formation. For a cluster with a virial mass of 7.4 X 10^14 Msol, inclusion of baryonic physics without feedback can boost the cross-section by as much as a factor of 3, in agreement with previous studies, whereas once AGN feedback is included this maximal figure falls to a factor of 2 at most. Typically, clusters simulated with DM and CSFBH physics have similar cross-sections for the production of giant arcs. We also investigate how baryonic physics affects the weak lensing properties of the simulated clusters by fitting NFW profiles to synthetic weak lensing data sets using a Markov Chain Monte Carlo approach, and by performing non-parametric mass reconstructions. Without the inclusion of AGN feedback, measured concentration parameters can be much larger than those obtained with AGN feedback, which are similar to the dark-matter only case.

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.

The problem of detecting dark matter filaments in the cosmic web is considered. Weak lensing is an ideal probe of dark matter, and therefore forms the basis of particularly promising detection methods. We consider and develop a number of weak lensing techniques that could be used to detect filaments in individual or stacked cluster fields, and apply them to synthetic lensing data sets in the fields of clusters from the Millennium Simulation. These techniques are multipole moments of the shear and convergence, mass reconstruction, and parameterized fits to filament mass profiles using a Markov Chain Monte Carlo approach. In particular, two new filament detection techniques are explored (multipole shear filters and Markov Chain Monte Carlo mass profile fits), and we outline the quality of data required to be able to identify and quantify filament profiles. We also consider the effects of large scale structure on filament detection. We conclude that using these techniques, there will be realistic prospects of detecting filaments in data from future space-based missions. The methods presented in this paper will be of great use in the identification of dark matter filaments in future surveys.

We use the galaxy stellar mass and halo merger tree information from the semi-analytic model galaxy catalogue of Font et al. (2009) to examine the accretion of galaxies into a large sample of groups and clusters, covering a wide range in halo mass (10E12.9 to 10E15.3 Msun/h), and selected from each of four redshift epochs (z=0, 0.5, 1.0 and 1.5). We find that clusters at all examined redshifts have accreted a significant fraction of their final galaxy populations through galaxy groups. A 10E14.5 Msun/h mass cluster at z=0 has, on average, accreted ~ 40% of its galaxies (Mstellar > 10E9 Msun/h) from halos with masses greater than 10E13 Msun/h. Further, the galaxies which are accreted through groups are more massive, on average, than galaxies accreted through smaller halos or from the field population. We find that at a given epoch, the fraction of galaxies accreted from isolated environments is independent of the final cluster or group mass. In contrast, we find that observing a cluster of the same halo mass at each redshift epoch implies different accretion rates of isolated galaxies, from 5-6 % per Gyr at z=0 to 15% per Gyr at z=1.5. We find that combining the existence of a Butcher Oemler effect at z=0.5 and the observations that galaxies within groups display significant environmental effects with galaxy accretion histories justifies striking conclusions. Namely, that the dominant environmental process must begin to occur in halos of 10E12 -- 10E13 Msun/h, and act over timescales of > 2 Gyrs. This argues in favor of a mechanism like "strangulation", in which the hot halo of a galaxy is stripped upon infalling into a more massive halo . This simple model predicts that by z=1.5 galaxy groups and clusters will display little to no environmental effects.

(Abridged) The thermal state of the intracluster medium results from a competition between gas cooling and heating. The heating comes from two distinct sources: gravitational heating from the collapse of the dark matter halo and thermal input from galaxy/black hole formation. However, a long standing problem has been that cosmological simulations based on smoothed particle hydrodynamics (SPH) and Eulerian mesh codes predict different results even when cooling and galaxy/black hole heating are switched off. Clusters formed in SPH simulations show near powerlaw entropy profiles, while those formed in mesh simulations develop a core and do not allow gas to reach such low entropies. Since the cooling rate is closely connected to the minimum entropy of the gas, the differences are of potentially key importance. In this paper, we investigate the origin of this discrepancy. By comparing simulations run using the GADGET-2 SPH code and the FLASH adaptive Eulerian mesh code, we show that the discrepancy arises during the idealised merger of two clusters. The difference is not sensitive to the resolution of our simulations, nor is it is due differences in the gravity solvers, Galilean non-invariance of the mesh code, or an effect of unsuitable artificial viscosity in the SPH code. Instead, we find that the difference is inherent to the treatment of eddies and fluid instabilities. These are suppressed in the SPH simulations, while the cluster mergers generate strong vortices in the mesh simulations that efficiently mix the fluid and erase the low entropy gas. Consequently, particles in the SPH simulations retain a close connection to their initial entropy, while this connection is much weaker in the mesh simulations. We discuss the potentially profound implications of these results.

Only ~10% of baryons in the universe are in the form of stars, yet most models of luminous structure formation have concentrated on the properties of the luminous stellar matter. In this paper we focus on the "flip side" of galaxy formation and investigate the properties of the material that is not presently locked up in galaxies. This "by-product" of galaxy formation can be observed as an X-ray emitting plasma (the intracluster medium, hereafter ICM) in groups and clusters, and we present a version of the Durham semi-analytic galaxy formation model GALFORM that allows us to investigate the properties of the ICM. As we would expect on the basis of gravitational scaling arguments, the previous model (presented in Bower et al. 2006) fails to reproduce even the most basic observed properties of the ICM; however, we present a simple modification to the model to allow for heat input into the ICM from the AGN "radio mode" feedback. This heating acts to expel gas from the X-ray luminous central regions of the host halo. With this modification, the model reproduces the observed gas mass fractions and luminosity-temperature relation of groups and clusters. Introducing the heating process into the model requires changes to a number of model parameters in order to retain a good match to the observed galaxy properties. With the revised parameters, the best fitting luminosity function is comparable to that presented in Bower et al. (2006). The new model makes a fundamental step forward, providing a unified model of galaxy and cluster ICM formation. However, the detailed comparison with the data is not completely satisfactory, and we highlight key areas for improvement.

Current models of galaxy formation predict satellite galaxies in groups and clusters that are redder than observed. We investigate the effect on the colours of satellite galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as the satellites orbit within their parent halo. We incorporate a model of the stripping process based on detailed hydrodynamic simulations within the Durham semi-analytic model of galaxy formation. The simulations show that the environment in groups and clusters is less aggressive than previously assumed. The main uncertainty in the model is the treatment of gas expelled by supernovae. With reasonable assumptions for the stripping of this material, we find that satellite galaxies are able to retain a significant fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs of cold, star forming gas and remaining blue for a relatively long period of time. A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is established and the colours of the satellite galaxies are in good agreement with the data. In addition, our model naturally accounts for the observed dependence of satellite colours on environment, from small groups to high mass clusters.

(Abridged) We use the simulations presented in Poole et al. 2006 to examine the effects of mergers on compact cool cores in X-ray clusters. We propose a scheme for classifying the morphology of clusters based on their surface brightness and entropy profiles. Three dominant morphologies emerge: two hosting compact cores and central temperatures which are cool (CCC systems) or warm (CWC systems) and one hosting extended cores which are warm (EWC systems). We find that CCC states are disrupted only after direct collisions between cluster cores in head-on collisions or during second pericentric passage in off-axis mergers. By the time they relax, our remnant cores have generally been heated to warm core (CWC or EWC) states but subsequently recover CCC states. The only case resulting in a long-lived EWC state is a slightly off-axis 3:1 merger for which the majority of shock heating occurs during the accretion of a low-entropy stream formed from the disruption of the secondary's core. Compression prevents core temperatures from falling until after relaxation thus explaining the observed population of relaxed CWC systems with no need to invoke AGN feedback. The morphological segregation observed in the L_x-T_x and beta-r_c scaling relations is reflected in our simulations as well. However, none of the cases we have studied produce sufficiently high remnant central entropies to account for the most under-luminous EWC systems observed. Lastly, systems which initially host central metallicity gradients do not yield merger remnants with flat metallicity profiles. Taken together, these results suggest that once formed, compact core systems are remarkably stable against disruption from mergers. It remains to be demonstrated exactly how the sizable observed population of extended core systems was formed.

(Abridged) X-ray clusters are conventionally divided into two classes: "cool core" (CC) clusters and "non-cool core" (NCC) clusters. Yet relatively little attention has been given to the origins of this dichotomy and, in particular, to the energetics and thermal histories of the two classes. We develop a model for the entropy profiles of clusters starting from the configuration established by gravitational shock heating and radiative cooling. At large radii, gravitational heating accounts for the observed profiles and their scalings well. However, at small and intermediate radii, radiative cooling and gravitational heating cannot be combined to explain the observed profiles of either type of cluster. The inferred entropy profiles of NCC clusters require that material is preheated prior to cluster collapse in order to explain the absence of low entropy (cool) material in these systems. We show that a similar modification is also required in CC clusters in order to match their properties at intermediate radii. In CC clusters, this modification is unstable, and an additional process is required to prevent cooling below a temperature of a few keV. We show that this can be achieved by adding a self-consistent AGN feedback loop in which the lowest-entropy, most rapidly cooling material is heated so that it rises buoyantly to mix with material at larger radii. The resulting model does not require fine tuning and is in excellent agreement with a wide variety of observational data. Some of the other implications of this model are briefly discussed.

(Abridged) We demonstrate that the tenet of hierarchical structure growth leads directly to a robust, falsifiable prediction for the correlation between stellar fraction (fstar) and total system mass (M500) of galaxy groups and clusters. This prediction is relatively insensitive to the details of baryonic physics or cosmological parameters. In particular, if the fstar-M500 relation is fixed and does not evolve with redshift, CDM models predict the logarithmic slope of this relation to be b>-0.3. This constraint can be weakened if the fstar-M500 relation evolves strongly, but this implies more stars must be formed in situ in groups at low redshift. Conservatively requiring that at least half the stars in groups were formed by z=1, the constraint from evolution models is b>-0.35. Since the most massive clusters (M500=1E15 Msun) are observed to have fstar=0.01, this means that groups with M500=5E13 Msun must have fstar<0.03. Recent observations by Gonzalez et al. (2007) indicate a much steeper relation, with fstar>0.04 in groups, leading to b=-0.64. If confirmed, this would rule out hierarchical structure formation models: today's clusters could not have been built from today's groups, or even from the higher-redshift progenitors of those groups. We perform a careful analysis of these and other data to identify the most important systematic uncertainties in their measurements. Although correlated uncertainties on stellar and total masses might explain the steep observed relation, the data are only consistent with theory if the observed group masses are systematically underestimated.

We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous halos of galaxies as they fall into massive groups and clusters. The sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ~30% of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott (1972), except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disk) and takes into account that stripping is not instantaneous but occurs on a characteristic timescale. The model reproduces the results of the simulations to within approximately 10% at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10% of the group or cluster into which it is falling (in which case the model under-predicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically-simulated galaxies and as a way to improve current semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.

The universal baryonic mass fraction (Omega_b/Omega_m) can be sensitively constrained using X-ray observations of galaxy clusters. In this paper, we compare the baryonic mass fraction inferred from measurements of the cosmic microwave background with the gas mass fractions (f_gas) of a large sample of clusters taken from the recent literature. In systems cooler than 4 keV, f_gas declines as the system temperature decreases. However, in higher temperature systems, f_gas(r500) converges to approx. (0.12 +/- 0.02)(h/0.72)^{-1.5}, where the uncertainty reflects the systematic variations between clusters at r500. This is significantly lower than the maximum-likelihood value of the baryon fraction from the recently released WMAP 3-year results. We investigate possible reasons for this discrepancy, including the effects of radiative cooling and non-gravitational heating, and conclude that the most likely solution is that Omega_m is higher than the best-fit WMAP value (we find Omega_m = 0.36^{+0.11}_{-0.08}), but consistent at the 2-sigma level. Degeneracies within the WMAP data require that sigma_8 must also be greater than the maximum likelihood value for consistency between the data sets.

The cosmic background due to the Sunyaev-Zeldovich (SZ) effect is expected to be the largest signal at mm and cm wavelengths at a resolution of a few arcminutes. We investigate some simple statistics of SZ maps and their scaling with the normalization of the matter power spectrum, sigma_8, as well as the effects of the unknown physics of the intracluster medium on these statistics. We show that the SZ background provides a significant background for SZ cluster searches, with the onset of confusion occurring around 10^{14} h^{-1} solar masses in a cosmology-dependent way, where confusion is defined as typical errors in recovered flux larger than 20%. The confusion limit, corresponds to the mass at which there are roughly ten clusters per square degree, with this number nearly independent of cosmology and cluster gas physics. Typical errors grow quickly as lower mass objects are included in the catalog. We also point out that there is nothing in particular about the rms of the filtered map that makes it especially well-suited for capturing aspects of the SZ effect, and other indicators of the one-point SZ probability distribution function are at least as well suited for the task. For example, the full width at half maximum of the one point probability distribution has a field-to-field scatter that is about 60% that of the rms. The simplest statistics of SZ maps are largely unaffected by cluster physics such aspreheating, although the impact of preheating is clear by eye in the maps.Studies aimed at learning about the physics of the intracluster medium will apparently require more specialized statistical indicators.

We use the suite of simulations presented in Poole et al (2006) to examine global X-ray and Sunyaev-Zel'dovich (SZ) observables for systems of merging relaxed X-ray clusters. In all of the scaling relations generated from the properties we have studied, we observe a generic evolution (in good qualitative agreement with previous authors): a rapid transient roughly along the mass scaling relations, a subsequent slow drift across the scatter until virialization, followed by a slow evolution along and up the mass scaling relations as cooling recovers in the cluster cores. However, this drift is not sufficient to account for the observed scatter in the scaling relations. We also study the effects of mergers on several theoretical temperature measures of the intracluster medium: emission weighted measures (T_ew), the spectroscopic-like measure proposed by Mazzotta et al (2004; T_sl) and plasma model fits to the integrated spectrum of the system (T_spec). We find that T_sl tracks T_spec for the entire duration of our mergers, illustrating that it remains a good tool for observational comparison even for highly disturbed systems. Furthermore, the transient temperature increases produced during first and second pericentric passage are 15-40% larger for T_ew than for T_sl or T_spec. This suggests that the effects of transient temperature increases on sigma_8 and Omega_M derived by Randall et al (2002) are over estimated. Lastly, we examine the X-ray SZ proxy proposed by Kravtsov et al (2006) and find that the tight mass scaling relation they predict remains secure through the entire duration of a merger event, independent of projection effects.

In McCarthy et al. (2007, astro-ph/0701335) we present an analysis of entropy generation during cluster mergers, using SPH simulations. This conference proceeding contribution summarizes some of the most important findings of that paper, without all the technical details. Our main result is that during cluster mergers, entropy is generated in two major shock episodes, that heat the gas from the inside-out, in a way that distributes most of the energy within the more massive clump. This is different from the "onion-skin" model adopted in many simple analytic models of cluster formation.

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

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

Using a set of high-resolution N-body/SPH cosmological simulations with identical initial conditions but run with different numerical setups, we investigate the influence of baryonic matter on the mass distribution of dark halos when radiative cooling is NOT included. We compare the concentration parameters of about 400 massive halos with virial mass from $10^{13}$ \Msun to $7.1 \times 10^{14}$ \Msun. We find that the concentration parameters for the total mass and dark matter distributions in non radiative simulations are on average larger by ~3% and 10% than those in a pure dark matter simulation. Our results indicate that the total mass density profile is little affected by a hot gas component in the simulations. After carefully excluding the effects of resolutions and spurious two-body heating between dark matter and gas particles, we conclude that the increase of the dark matter concentration parameters is due to interactions between baryons and dark matter. We demonstrate this with the aid of idealized simulations of two-body mergers. The results of individual halos simulated with different mass resolutions show that the gas profiles of densities, temperature and entropy are subjects of mass resolution of SPH particles. In particular, we find that in the inner parts of halos, as the SPH resolution increases the gas density becomes higher but both the entropy and temperature decrease.

We report on the analysis of a suite of SPH simulations (incorporating cooling and star formation) of mergers involving idealised X-ray clusters whose initial conditions resemble relaxed clusters with cool compact cores observed by Chandra and XMM. The simulations sample the most interesting, theoretically plausible, range of impact parameters and progenitor mass ratios. We find that all mergers evolve via a common progression. We illustrate this progression in the projected gas density, X-ray surface brightness, SZ, temperature, and gas entropy maps. Several different classes of transient ``cold front''-like features can arise over the course of a merger. We find that all of these classes are present in Chandra and XMM observations of merging systems and propose a naming scheme for these features: ``comet-like'' tails, bridges, plumes, streams and edges. In none of the cases considered do the initial cool compact cores of the primary and the secondary get destroyed during the course of the mergers. We quantify the evolving morphology of our mergers using centroid variance, power ratios and offset between the X-ray and the projected mass maps. We find that the centroid variance best captures the dynamical state of the cluster. Placing the system at z=0.1, we find that all easily identified observable traces of the secondary disappear from a simulated 50 ks Chandra image following the second pericentric passage. The system, however, takes approximately 2 additional Gyrs to relax and virialize. Temperature fluctuations at the level of 20% can persist in the final systems well past the point of virialization, suggesting that that the existence of temperature fluctuations, in and of themselves, do not necessarily indicate a disturbed or unrelaxed system.

(Abridged) We revisit the scaling relationships between the dark matter mass and observed X-ray luminosity and temperature of galaxy clusters and groups in the local Universe. Specifically, we compare recent observations with analytic models of the intracluster medium in which the gas entropy distribution has been shifted by a variable amount, K0, to investigate the origin of the scatter in these scaling relations, and its influence on the luminosity and temperature functions. We find that variations in halo concentration or formation epoch (which might determine the time available for low entropy gas to cool out) are insufficient to explain the amount of scatter in the mass-luminosity relation. Instead, a range of entropy floors at a fixed halo mass, spanning approximately 50 keV sq.cm to 700 keV sq.cm, is required to match the data. This range is likely related to the variance in heating and/or cooling efficiency from halo to halo. We demonstrate that these models are consistent with the observed temperature and luminosity functions of clusters, if scatter in the mass-luminosity relation is accounted for. Finally, we present predictions for the redshift evolution of these scaling relations and luminosity/temperature functions. Comparison with recent data at z<0.7 shows reasonable agreement with a model that assumes a median entropy floor of K0=200 keV sq.cm. When observations are extended to group scales (kT < 1 keV), this evolution will have the potential to discriminate between an entropy floor that is independent of redshift (for example, in a preheating scenario) and one that depends on the cooling time of the halo.

High resolution X-ray observations indicate that the entropy profiles in the central regions of some massive cooling flow clusters are well approximated by powerlaws. McCarthy and coworkers recently accounted for this trend with an analytic model that includes the detailed effects of radiative cooling. Interestingly, these authors found that cooling (and subsequent inflow of the gas) naturally establishes approximate steady-state powerlaw entropy profiles in the cores of clusters. However, the origin of this behavior and its dependence on initial conditions have yet to be elucidated. In the present study, we explain this trend in the context of the self-similar cooling wave model developed previously by Bertschinger (1989). It is shown that the logarithmic slope of the entropy profile in the cores of relaxed cooling flow clusters is given by a simple analytic function that depends only on the logarithmic slopes of the local gravitational potential and the cooling function. We discuss a number of potentially interesting uses of the above result, including: (1) a way of measuring the shapes of gravitational potentials of cooling flow clusters (which may, for example, be compared with the results of the standard hydrostatic equilibrium method); (2) a simple method for constructing realistic analytic cluster models that include the effects of radiative cooling; and (3) a test of the reliability of cooling routines implemented in analytic models and hydrodynamic simulations.

(Abridged) Theoretical models that include only gravitationally-driven processes fail to match the observed mean X-ray properties of clusters. As a result, there has recently been increased interest in models in which either radiative cooling or entropy injection play a central role in mediating the properties of the intracluster medium. Both sets of models give reasonable fits to the mean properties of clusters, but cooling only models result in fractions of cold baryons in excess of observationally established limits and the simplest entropy injection models do not treat the "cooling core" structure present in many clusters and cannot account for entropy profiles revealed by recent X-ray observations. We consider models that marry radiative cooling with entropy injection, and confront model predictions for the global and structural properties of massive clusters with the latest X-ray data. The models successfully and simultaneously reproduce the observed L-T and L-M relations, yield detailed entropy, surface brightness, and temperature profiles in excellent agreement with observations, and predict a cooled gas fraction that is consistent with observational constraints. The model also provides a possible explanation for the significant intrinsic scatter present in the L-T and L-M relations and provides a natural way of distinguishing between clusters classically identified as "cooling flow" clusters and dynamically relaxed "non-cooling flow" clusters. The former correspond to systems that had only mild levels (< 300 keV cm^2) of entropy injection, while the latter are identified as systems that had much higher entropy injection. This is borne out by the entropy profiles derived from Chandra and XMM-Newton.

Photoevaporation may provide an explanation for the short lifetimes of disks around young stars. With the exception of neutral oxygen lines, the observed low-velocity forbidden line emission from T Tauri stars can be reproduced by photoevaporating models. The natural formation of a gap in the disk at several AU due to photoevaporation and viscous spreading provides a possible halting mechanism for migrating planets and an explanation for the abundance of observed planets at these radii.

We examine the ability of photoevaporative disk winds to explain the low-velocity components observed in the forbidden line spectra of low-mass T Tauri stars. Using the analytic model of Shu, Johnstone, & Hollenbach (1993) and Hollenbach et al. (1994) as a basis, we examine the characteristics of photoevaporative outflows with hydrodynamic simulations. General results from the simulations agree well with the analytic predictions, although some small differences are present. Most importantly, the flow of material from the disk surface develops at smaller radii than in the analytic approximations and the flow-velocity from the disk surface is only one-third the sound speed. A detailed presentation of observational consequences of the model is given, including predicted line widths, blue-shifts, and integrated luminosities of observable sulfur and nitrogen emission lines. We demonstrate that these predictions are in agreement with current observational data on the low-velocity forbidden line emission of ionized species from T Tauri stars. This is in contrast with magnetic wind models, which systematically under-predict these forbidden line luminosities. However, the present model cannot easily account for the luminosities of neutral oxygen lines in T Tauri stars.

Recent X-ray observations have highlighted clusters that lack entropy cores. At first glance, these results appear to invalidate the preheated ICM models. We show that a self-consistent preheating model, which factors in the effects of radiative cooling, is in excellent agreement with the observations. Moreover, the model naturally explains the intrinsic scatter in the L-T relation, with ``cooling flow'' and ``non-cooling flow'' systems corresponding to mildly and strongly preheated systems, respectively. We discuss why preheating ought to be favoured over merging as a mechanism for the origin of ``non-cooling flow'' clusters.

Studies of cluster X-ray scaling relations have led to suggestions that non-gravitational processes, e.g., radiative cooling and/or "preheating", have significantly modified the entropy of the intracluster medium (ICM). For the first time, we test this hypothesis through a comparison of predicted thermal Sunyaev-Zeldovich (SZ) effect scaling relations with available data from the literature. One of the relations that we explore, in principle, depends solely on SZ effect observations, thus offering an X-ray independent probe of the ICM. A detailed comparison of the theoretical relations with the largest compilation of high redshift SZ effect data to date indicates that the presence of an entropy floor is favored by the data. Furthermore, the inferred level of that floor, 300 keV cm^2 or higher, is comparable to that found in studies of X-ray scaling relations of nearby massive clusters. Thus, we find no evidence for significant evolution of the entropy floor out to z = 0.7. We further demonstrate that the high quality data to be obtained from the upcoming Sunyaev-Zeldovich Array (SZA) and the (soon-to-be) upgraded Owens Valley Radio Observatory (OVRO) array will open powerful new windows into the properties of the ICM. Specifically, the new measurements will allow for accurate measurements of the ICM entropy for even the most distant galaxy clusters.

A common feature of the X-ray bubbles observed in Chandra images of some "cooling flow" clusters is that they appear to be surrounded by bright, cool shells. Temperature maps of a few nearby luminous clusters reveal that the shells consist of the coolest gas in the clusters -- much cooler than the surrounding medium. Using simple models, we study the effects of this cool emission on the inferred cooling flow properties of clusters. We find that the introduction of bubbles into model clusters that do NOT have cooling flows results in temperature and surface brightness profiles that resemble those seen in nearby "cooling flow" clusters. They also approximately reproduce the recent XMM-Newton and Chandra observations of a high minimum temperature of 1-3 keV. Hence, bubbles, if present, must be taken into account when inferring the physical properties of the ICM. In the case of some clusters, bubbles may account entirely for these observed features, calling into question their designation as clusters with cooling flows. However, since not all nearby "cooling flow" clusters show bubble-like features, we suggest that there may be a diverse range of physical phenomena that give rise to the same observed features.

X-ray observations of an "entropy floor" in nearby groups and clusters of galaxies offer evidence that important non-gravitational processes, such as radiative cooling and/or "preheating", have strongly influenced the evolution of the intracluster medium (ICM). We examine how the presence of an entropy floor modifies the thermal Sunyaev-Zeldovich (SZ) effect. A detailed analysis of scaling relations between X-ray and SZ effect observables and also between the two primary SZ effect observables is presented. We find that relationships between the central Compton parameter and the temperature or mass of a cluster are extremely sensitive to the presence of an entropy floor. The same is true for correlations between the integrated Compton parameter and the X-ray luminosity or the central Compton parameter. In fact, if the entropy floor is as high as inferred in recent analyses of X-ray data, a comparison of these correlations with both current and future SZ effect observations should show a clear signature of this excess entropy. Moreover, because the SZ effect is redshift-independent, the relations can potentially be used to track the evolution of the cluster gas and possibly discriminate between the possible sources of the excess entropy. To facilitate comparisons with observations, we provide analytic fits to these scaling relations.

Recent efforts to account for the observed X-ray luminosity - temperature relation of galaxy clusters has led to suggestions that the ICM has an apparent ``entropy floor'' at or above the level of 300 keV cm^2. Here, we propose new tests based on the thermal Sunyaev-Zeldovich effect and on the cluster gas mass - temperature trend (from X-ray data) to probe the level of excess entropy in the ICM. We show that these new tests lend further support to the case for a high entropy floor in massive clusters.

Recent X-ray observations have been used to demonstrate that the cluster gas mass-temperature relation is steeper than theoretical self-similar predictions drawn from numerical simulations that consider the evolution of the cluster gas through the effects of gravity and shock heating alone. One possible explanation for this is that the gas mass fraction is not constant across clusters of different temperature, as is usually assumed. Observationally, however, there is no compelling evidence for gas mass fraction variation, especially in the case of hot clusters. Seeking an alternative physical explanation for the observed trends, we investigate the role in the cluster gas mass - temperature relation of the preheating of the intracluster medium by some arbitrary source for clusters with emission-weighted mean temperatures of greater than about 3 keV. Making use of the physically-motivated, analytic models developed in 2002 by Babul and coworkers, we find that preheating does, indeed, lead to a steeper relation. This is in agreement with previous theoretical studies on the relation. However, in apparent conflict with these studies, we argue that a ``high'' level of entropy injection is required to match observations. In particular, an entropy floor of greater than about 300 keV cm^2 is required. We also present a new test, namely, the study of the relation within different fixed radii. This allows one to indirectly probe the density profiles of clusters, since it samples different fractions of the virial radius for clusters of different temperature. This test also confirms that a high level of preheating is required to match observations.

We report the results of experiments aimed at reducing the major problem with cooling flow models of rich cluster X-ray sources: the fact that most of the cooled gas or its products have not been found. Here we show that much of the X-ray emission usually attributed to cooling flows can, in fact, be modeled by a power-law component which is indicative of a source(s) other than thermal bremsstrahlung from the intracluster medium. We find that adequate simultaneous fits to ROSAT PSPCB and ASCA GIS/SIS spectra of the central regions of ten clusters are obtained for two-component models that include a thermal plasma component that is attributable to hot intracluster gas and a power-law component that is likely generated by compact sources and/or extended non-thermal emission. For five of the clusters that purportedly have massive cooling flows, the best-fit models have power-law components that contribute $\sim$ 30 % of the total flux (0.14 - 10.0 keV) within the central 3 arcminutes. Because cooling flow mass deposition rates are inferred from X-ray fluxes, our finding opens the possibility of significantly reducing cooling rates.