Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $\sigma_\mathrm{v}$ with the density $\rho$ can be used to define the pseudo-entropy $S(r)=\sigma_\mathrm{v}^2/\rho^{2/3}$, whose profile is well-described by a simple power-law $S\propto r^{\alpha}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profile as traced by substructures, while pseudo-entropy profiles as traced by DM particles are characterized by an almost universal behavior and can be used as a low-scatter ($\sim6-7\%$) mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius.

We present the first metal abundance profiles for a representative sample of massive clusters. Our measures extend to $R_{500}$ and are corrected for a systematic error plaguing previous outskirt estimates. Our profiles flatten out at large radii, admittedly not a new result, however the radial range and representative nature of our sample extends its import well beyond previous findings. We find no evidence of segregation between cool-core and non-cool-core systems beyond $\sim 0.3 R_{500}$, implying that, as was found for thermodynamic properties (Ghirardini et al, 2019), the physical state of the core does not affect global cluster properties. Our mean abundance within $R_{500}$ shows a very modest scatter, $< $15%, suggesting the enrichment process must be quite similar in all these massive systems. This is a new finding and has significant implications on feedback processes. Together with results from thermodynamic properties presented in a previous X-COP paper, it affords a coherent picture where feedback effects do not vary significantly from one system to another. By combing ICM with stellar measurements we have found the amount of Fe diffused in the ICM to be about ten times higher than that locked in stars. Although our estimates suggest, with some strength, that the measured iron mass in clusters is well in excess of the predicted one, systematic errors prevent us from making a definitive statement. Further advancements will only be possible when systematic uncertainties, principally those associated to stellar masses, both within and beyond $R_{500}$, can be reduced.

We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range $0<z<4$, comparing them to recent observational studies; we also investigated the effect of varying the parameters of the star formation model on galaxy properties such as SFR, star-formation efficiency, and gas fraction. We analyze a set of zoom-in cosmological hydrodynamical simulations centred on twelve clusters. The simulations are carried out with the GADGET-3 TreePM/SPH code which includes various subgrid models to treat unresolved baryonic physics, including AGN feedback. Simulations do not reproduce the high values of SFR observed within protoclusters cores, where the values of SFR are underpredicted by a factor $\gtrsim 4$ both at $z\sim2$ and $z\sim 4$. The difference arises as simulations are unable to reproduce the observed starburst population and is worsened at $z\sim 2$ because simulations underpredict the normalization of the main sequence of star forming galaxies (i.e., the correlation between stellar mass and SFR) by a factor of $\sim 3$. As the low normalization of the main sequence seems to be driven by an underestimated gas fraction, it remains unclear whether numerical simulations miss starburst galaxies due to a too low predicted gas fractions or too low star formation efficiencies. Our results are stable against varying several parameters of the star formation subgrid model and do not depend on the details of the AGN feedback.

Contradictory results have been reported on the time evolution of the alignment between clusters and their Brightest Cluster Galaxy (BCG). We study this topic by analyzing cosmological hydro-simulations of 24 massive clusters with $M_{200}|_{z=0} \gtrsim 10^{15}\, M_\odot$, plus 5 less massive with $1 \times 10^{14} \lesssim M_{200}|_{z=0} \lesssim 7 \times 10^{14}\, M_\odot$, which have already proven to produce realistic BCG masses. We compute the BCG alignment with both the distribution of cluster galaxies and the dark matter (DM) halo. At redshift $z=0$, the major axes of the simulated BCGs and their host cluster galaxy distributions are aligned on average within 20$^\circ$. The BCG alignment with the DM halo is even tighter. The alignment persists up to $z\lesssim2$ with no evident evolution. This result continues, although with a weaker signal, when considering the projected alignment. The cluster alignment with the surrounding distribution of matter ($3R_{200}$) is already in place at $z\sim4$ with a typical angle of $35^\circ$, before the BCG-Cluster alignment develops. The BCG turns out to be also aligned with the same matter distribution, albeit always to a lesser extent. These results taken together might imply that the BCG-Cluster alignment occurs in an outside-in fashion. Depending on their frequency and geometry, mergers can promote, destroy or weaken the alignments. Clusters that do not experience recent major mergers are typically more relaxed and aligned with their BCG. In turn, accretions closer to the cluster elongation axis tend to improve the alignment as opposed to accretions closer to the cluster minor axis.

Joint lensing and dynamical mass profile determinations of galaxy clusters are an excellent tool to constrain modification of gravity at cosmological scales. However, search for tiny departures from General Relativity calls for an accurate control of the systematics affecting the method. In this analysis we concentrate on the systematics in the reconstruction of mass profiles from the dynamics of cluster member galaxies, while assuming that lensing provides unbiased mass profile reconstructions. In particular, in the case study of linear $f(R)$ gravity, we aim at veryfying whether in realistic simulations of cluster formation a spurious detection of departure from GR can be detected due to violation of the main assumptions (e.g. dynamical equilibrium and spherical symmetry) on which the method is based. We aim at identifying and calibrating the impact of those systematics by analyzing a set of Dark Matter halos taken from $\Lambda$CDM N-body cosmological simulations performed with the GADGET-3 code. [...] If no selection criteria are applied, $\sim 60\%$ of clusters in a $\Lambda$CDM Universe (where GR is assumed) produce a spurious detection of modified gravity. We find that the probability of finding cluster in agreement with GR predictions $P_{GR}$ mainly depends on the properties of the halo's projected phase-space and on shape orientation of the cluster along the line-of-sight projection. We define two observational criteria which correlate with the probability to find clusters in agreement with GR predictions and which can be used to select [...] those objects that are more suitable for the application of the proposed method. In particular, we find that according to these criteria the percentage of spurious detection can be lowered down to $\sim 20\%$ in the best case. Our results are relevant in view of data that will be available with the next generation surveys.

Accurate and precise measurements of masses of galaxy clusters are key to derive robust constraints on cosmological parameters. Rising evidence from observations, however, confirms that X-ray masses, obtained under the assumption of hydrostatic equilibrium, might be underestimated, as previously predicted by cosmological simulations. We analyse more than 300 simulated massive clusters, from `The Three Hundred Project', and investigate the connection between mass bias and several diagnostics extracted from synthetic X-ray images of these simulated clusters. We find that the azimuthal scatter measured in 12 sectors of the X-ray flux maps is a statistically significant indication of the presence of an intrinsic (i.e. 3D) clumpy gas distribution. We verify that a robust correction to the hydrostatic mass bias can be inferred when estimates of the gas inhomogeneity from X-ray maps (such as the azimuthal scatter or the gas ellipticity) are combined with the asymptotic external slope of the gas density or pressure profiles, which can be respectively derived from X-ray and millimetric (Sunyaev-Zeldovich effect) observations. We also obtain that mass measurements based on either gas density and temperature or gas density and pressure result in similar distributions of the mass bias. In both cases, we provide corrections that help reduce both the dispersion and skewness of the mass bias distribution. These are effective even when irregular clusters are included leading to interesting implications for the modelling and correction of hydrostatic mass bias in cosmological analyses of current and future X-ray and SZ cluster surveys.

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 present simulations of galaxy formation, based on the GADGET-3 code, in which a sub-resolution model for star formation and stellar feedback is interfaced with a new model for AGN feedback. Our sub-resolution model describes a multiphase ISM, accounting for hot and cold gas within the same resolution element: we exploit this feature to investigate the impact of coupling AGN feedback energy to the different phases of the ISM over cosmic time. Our fiducial model considers that AGN feedback energy coupling is driven by the covering factors of the hot and cold phases. We perform a suite of cosmological hydrodynamical simulations of disc galaxies ($M_{\rm halo, \, DM} \simeq 2 \cdot 10^{12}$ M$_{\odot}$, at $z=0$), to investigate: $(i)$ the effect of different ways of coupling AGN feedback energy to the multiphase ISM; $(ii)$ the impact of different prescriptions for gas accretion (i.e. only cold gas, both cold and hot gas, with the additional possibility of limiting gas accretion from cold gas with high angular momentum); $(iii)$ how different models of gas accretion and coupling of AGN feedback energy affect the coevolution of supermassive BHs and their host galaxy. We find that at least a share of the AGN feedback energy has to couple with the diffuse gas, in order to avoid an excessive growth of the BH mass. When the BH only accretes cold gas, it experiences a growth that is faster than in the case in which both cold and hot gas are accreted. If the accretion of cold gas with high angular momentum is reduced, the BH mass growth is delayed, the BH mass at $z=0$ is reduced by up to an order of magnitude, and the BH is prevented from accreting below $z \lesssim 2$, when the galaxy disc forms.

Recently, relations connecting the SMBH mass of central galaxies and global properties of the hosting cluster, such as temperature and mass, were observed. We investigate the correlation between SMBH mass and cluster mass and temperature, their establishment and evolution. We compare their scatter to that of the classical $M_{\rm BH}-M_{\rm BCG}$ relation. We study how gas accretion and BH-BH mergers contribute to SMBH growth across cosmic time. We employed 135 groups and clusters with a mass range $1.4\times 10^{13}M_{\odot}-2.5\times 10^{15} M_{\odot}$ extracted from a set of 29 zoom-in cosmological hydro-dynamical simulations where the baryonic physics is treated with various sub-grid models, including feedback by AGN. In our simulations we find that $M_{\rm BH}$ correlates well with $M_{500}$ and $T_{500}$, with the scatter around these relations compatible within $2\sigma$ with the scatter around $M_{\rm BH}-M_{\rm BCG}$ at $z=0$. The $M_{\rm BH}-M_{500}$ relation evolves with time, becoming shallower at lower redshift as a direct consequence of hierarchical structure formation. On average, in our simulations the contribution of gas accretion to the total SMBH mass dominates for the majority of the cosmic time ($z>0.4$), while in the last 2 Gyr the BH-BH mergers become a larger contributor. During this last process, substructures hosting SMBHs are disrupted in the merger process with the BCG and the unbound stars enrich the diffuse stellar component rather than increase BCG mass. From the results obtained in our simulations with simple sub-grid models we conclude that the scatter around the $M_{\rm BH}-T_{500}$ relation is comparable to the scatter around the $M_{\rm BH}-M_{\rm BCG}$ relation and that, given the observational difficulties related to the estimation of the BCG mass, clusters temperature and mass can be a useful proxy for the SMBHs mass, especially at high redshift.

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by $M_\bullet-L_{\rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_\bullet-\sigma_{\rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_\bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

Metals form an essential part of the Universe at all scales. Without metals we would not exist, and the Cosmos would look completely different. Metals are primarily born through nuclear processes in stars. They leave their cradles through winds or explosions, and then start their journey through space. This can lead them in and out of astronomical objects on all scales, ranging from comets, planets, stars, entire galaxies, groups and clusters of galaxies to the largest structures of the Universe. Their wanderings are fundamental in determining how these objects, and the entire universe, evolve. In addition, their bare presence can be used to trace what these structures look like. The scope of this paper is to highlight the most important open astrophysical problems that will be central in the next decades and for which a deep understanding of the Universe-wandering metals, their physical and kinematical states and their chemical composition represents the only viable solution. The majority of these studies can only be efficiently performed through High Resolution Spectroscopy in the soft X-ray band.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

High performance computing numerical simulations are today one of the more effective instruments to implement and study new theoretical models, and they are mandatory during the preparatory phase and operational phase of any scientific experiment. New challenges in Cosmology and Astrophysics will require a large number of new extremely computationally intensive simulations to investigate physical processes at different scales. Moreover, the size and complexity of the new generation of observational facilities also implies a new generation of high performance data reduction and analysis tools pushing toward the use of Exascale computing capabilities. Exascale supercomputers cannot be produced today. We discuss the major technological challenges in the design, development and use of such computing capabilities and we will report on the progresses that has been made in the last years in Europe, in particular in the framework of the ExaNeSt European funded project. We also discuss the impact of this new computing resources on the numerical codes in Astronomy and Astrophysics.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

We present cosmological hydrodynamical simulations including atomic and molecular non-equilibrium chemistry, multi-frequency radiative transfer (0.7-100 eV sampled over 150 frequency bins) and stellar population evolution to investigate the host candidates of the seeds of supermassive black holes coming from direct collapse of gas in primordial haloes (direct-collapse black holes, DCBHs). We consistently address the role played by atomic and molecular cooling, stellar radiation and metal spreading of C, N, O, Ne, Mg, Si, S, Ca, Fe, etc. from primordial sources, as well as their implications for nearby quiescent proto-galaxies under different assumptions for early source emissivity, initial mass function and metal yields. We find that putative DCBH host candidates need powerful primordial stellar generations, since common solar-like stars and hot OB-type stars are neither able to determine the conditions for direct collapse nor capable of building up a dissociating Lyman-Werner background radiation field. Thermal and molecular features of the identified DCBH host candidates in the scenario with very massive primordial stars seem favourable, with illuminating Lyman-Werner intensities featuring values of 1-50 J21. Nevertheless, additional non-linear processes, such as merger events, substructure formation, rotational motions and photo-evaporation, should inhibit pure DCBH formation in 2/3 of the cases. Local turbulence may delay gas direct collapse almost irrespectively from other environmental conditions. The impact of large Lyman-Werner fluxes at distances smaller than 5 kpc is severely limited by metal pollution.

We perform a suite of cosmological hydrodynamical simulations of disc galaxies, with zoomed-in initial conditions leading to the formation of a halo of mass $M_{\rm halo, \, DM} \simeq 2 \cdot 10^{12}$ M$_{\odot}$ at redshift $z=0$. These simulations aim at investigating the chemical evolution and the distribution of metals in a disc galaxy, and at quantifying the effect of $(i)$ the assumed IMF, $(ii)$ the adopted stellar yields, and $(iii)$ the impact of binary systems originating SNe Ia on the process of chemical enrichment. We consider either a Kroupa et al. (1993) or a more top-heavy Kroupa (2001) IMF, two sets of stellar yields and different values for the fraction of binary systems suitable to give rise to SNe Ia. We investigate stellar ages, SN rates, stellar and gas metallicity gradients, and stellar $\alpha$-enhancement in simulations, and compare predictions with observations. We find that a Kroupa et al. (1993) IMF has to be preferred when modelling late-type galaxies in the local universe. On the other hand, the comparison of stellar metallicity profiles and $\alpha$-enhancement trends with observations of Milky Way stars shows a better agreement when a Kroupa (2001) IMF is assumed. Comparing the predicted SN rates and stellar $\alpha$-enhancement with observations supports a value for the fraction of binary systems producing SNe Ia of $0.03$, at least for late-type galaxies and for the considered IMFs. Adopted stellar yields are crucial in regulating cooling and star formation, and in determining patterns of chemical enrichment for stars, especially for those located in the galaxy bulge.

Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics \texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations agrees well in trend and evolution with the observational results obtained from two data samples characterised by a wide range of masses and a large redshift coverage. We find that the iron abundance in the cluster core ($r<0.1R_{500}$) does not correlate with the temperature nor presents a significant evolution. The scale invariance is confirmed when the metallicity is related directly to the total mass. The slope of the best-fitting relations is shallow ($\beta\sim-0.1$) in the innermost regions ($r<0.5R_{500}$) and consistent with zero outside. We investigate the impact of the AGN feedback and find that it plays a key role in producing a constant value of the outskirts metallicity from groups to clusters. This finding additionally supports the picture of early enrichment.

Answers to the metal production of the Universe can be found in galaxy clusters, notably within their Intra-Cluster Medium (ICM). The X-ray Integral Field Unit (X-IFU) on board the next-generation European X-ray observatory Athena (2030s) will provide the necessary leap forward in spatially-resolved spectroscopy required to disentangle the intricate mechanisms responsible for this chemical enrichment. In this paper, we investigate the future capabilities of the X-IFU in probing the hot gas within galaxy clusters. From a test sample of four clusters extracted from cosmological hydrodynamical simulations, we present comprehensive synthetic observations of these clusters at different redshifts (up to z = 2) and within the scaled radius R500 performed using the instrument simulator SIXTE. Through 100 ks exposures, we demonstrate that the X-IFU will provide spatially-resolved mapping of the ICM physical properties with little to no biases (<5%) and well within statistical uncertainties. The detailed study of abundance profiles and abundance ratios within R500 also highlights the power of the X-IFU in providing constraints on the various enrichment models. From synthetic observations out to z = 2, we also quantify its ability to track the chemical elements across cosmic time with excellent accuracy, and thereby to investigate the evolution of metal production mechanisms as well as the link to the stellar initial mass-function. Our study demonstrates the unprecedented capabilities of the X-IFU in unveiling the properties of the ICM but also stresses the data analysis challenges faced by future high-resolution X-ray missions such as Athena.

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

State-of-the-art cosmological hydrodynamical simulations have star particles with typical mass between $\sim 10^8$ and $\sim 10^3$ M$_{\odot}$ according to resolution, and treat them as simple stellar populations. On the other hand, observations in nearby galaxies resolve individual stars and provide us with single star properties. An accurate and fair comparison between predictions from simulations and observations is a crucial task. We introduce a novel approach to consistently populate star particles with stars. We developed a technique to generate a theoretical catalogue of mock stars whose characteristics are derived from the properties of parent star particles from a cosmological simulation. Also, a library of stellar evolutionary tracks and synthetic spectra is used to mimic the photometric properties of mock stars. The aim of this tool is to produce a database of synthetic stars from the properties of parent star particles in simulations: such a database represents the observable stellar content of simulated galaxies and allows a comparison as accurate as possible with observations of resolved stellar populations. With this innovative approach we are able to provide a colour-magnitude diagram from a cosmological hydrodynamical simulation. This method is flexible and can be tailored to fit output of different codes used for cosmological simulations. Also, it is of paramount importance with ongoing survey data releases (e.g. GAIA and surveys of resolved stellar populations), and will be useful to predict properties of stars with peculiar chemical features and to compare predictions from hydrodynamical models with data of different tracers of stellar populations.

It has been known for decades that the observed number of baryons in the local universe falls about 30-40% short of the total number of baryons predicted by Big-Bang Nucleosynthesis, as inferred from density fluctuations of the Cosmic Microwave Background and seen during the first 2-3 billion years of the universe in the so called Lyman-alpha Forest. A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas - hydrogen - is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios. Indeed, despite the large observational efforts, only a few marginal claims of detection have been made so far. Here we report observations of two absorbers of highly ionized oxygen (OVII) in the high signal-to-noise-ratio X-ray spectrum of a quasar at redshift >0.4. These absorbers show no variability over a 2-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The OVII systems lie in regions characterized by large (x4 compared to average) galaxy over-densities and their number (down to the sensitivity threshold of our data), agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium (WHIM). We conclude that the missing baryons have been found.

We analyze the stellar growth of Brightest Cluster Galaxies (BCGs) produced by cosmological zoom-in hydrodynamical simulations of the formation of massive galaxy clusters. The evolution of the stellar mass content is studied considering different apertures, and tracking backwards either the main progenitor of the $z=0$ BCG or that of the cluster hosting the BCG at $z=0$. Both methods lead to similar results up to $z \simeq 1.5$. The simulated BCGs masses at $z=0$ are in agreement with recent observations. In the redshift interval from $z=1$ to $z=0$ we find growth factors 1.3, 1.6 and 3.6 for stellar masses within 30kpc, 50kpc and 10% of $R_{500}$ respectively. The first two factors, and in general the mass evolution in this redshift range, are in agreement with most recent observations. The last larger factor is similar to the growth factor obtained by a semi-analytical model (SAM). Half of the star particles that end up in the inner 50 kpc was typically formed by redshift $\sim$ 3.7, while the assembly of half of the BCGs stellar mass occurs on average at lower redshifts $\sim 1.5$. This assembly redshift correlates with the mass attained by the cluster at high $z \gtrsim 1.3$, due to the broader range of the progenitor clusters at high-$z$. The assembly redshift of BCGs decreases with increasing apertures. Our results are compatible with the {\it inside-out} scenario. Simulated BCGs could lack intense enough star formation (SF) at high redshift, while possibly exhibit an excess of residual SF at low redshift.

We implement a state-of-the-art treatment of the processes affecting the production and Interstellar Medium (ISM) evolution of carbonaceous and silicate dust grains within SPH simulations. We trace the dust grain size distribution by means of a two-size approximation. We test our method on zoom-in simulations of four massive ($M_{200} \geq 3 \times 10^{14} M_{\odot}$) galaxy clusters. We predict that during the early stages of assembly of the cluster at $z \gtrsim 3$, where the star formation activity is at its maximum in our simulations, the proto-cluster regions are rich of dusty gas. Compared to the case in which only dust production in stellar ejecta is active, if we include processes occurring in the cold ISM,the dust content is enhanced by a factor $2-3$. However, the dust properties in this stage turn out to be significantly different than those observationally derived for the {\it average} Milky Way dust, and commonly adopted in calculations of dust reprocessing. We show that these differences may have a strong impact on the predicted spectral energy distributions. At low redshift in star forming regions our model reproduces reasonably well the trend of dust abundances over metallicity as observed in local galaxies. However we under-produce by a factor of 2 to 3 the total dust content of clusters estimated observationally at low redshift, $z \lesssim 0.5$ using IRAS, Planck and Herschel satellites data. This discrepancy can be solved by decreasing the efficiency of sputtering which erodes dust grains in the hot Intracluster Medium (ICM).

Observational cosmology is passing through a unique moment of grandeur with the amount of quality data growing fast. However, in order to better take advantage of this moment, data analysis tools have to keep up the pace. Understanding the effect of baryonic matter on the large-scale structure is one of the challenges to be faced in cosmology. In this work, we have thoroughly studied the effect of baryonic physics on different lensing statistics. Making use of the Magneticum Pathfinder suite of simulations we show that the influence of luminous matter on the 1-point lensing statistics of point sources is significant, enhancing the probability of magnified objects with $\mu>3$ by a factor of $2$ and the occurrence of multiple-images by a factor $5-500$ depending on the source redshift and size. We also discuss the dependence of the lensing statistics on the angular resolution of sources. Our results and methodology were carefully tested in order to guarantee that our uncertainties are much smaller than the effects here presented.

The uniformity of the intra-cluster medium (ICM) enrichment level in the outskirts of nearby galaxy clusters suggests that chemical elements were deposited and widely spread into the intergalactic medium before the cluster formation. This observational evidence is supported by numerical findings from cosmological hydrodynamical simulations, as presented in Biffi et al. (2017), including the effect of thermal feedback from active galactic nuclei. Here, we further investigate this picture, by tracing back in time the spatial origin and metallicity evolution of the gas residing at z=0 in the outskirts of simulated galaxy clusters. In these regions, we find a large distribution of iron abundances, including a component of highly-enriched gas, already present at z=2. At z>1, the gas in the present-day outskirts was distributed over tens of virial radii from the the main cluster and had been already enriched within high-redshift haloes. At z=2, about 40% of the most Fe-rich gas at z=0 was not residing in any halo more massive than 1e11 Msun/h in the region and yet its average iron abundance was already 0.4, w.r.t. the solar value by Anders & Grevesse (1989). This confirms that the in situ enrichment of the ICM in the outskirts of present-day clusters does not play a significant role, and its uniform metal abundance is rather the consequence of the accretion of both low-metallicity and pre-enriched (at z>2) gas, from the diffuse component and through merging substructures. These findings do not depend on the mass of the cluster nor on its core properties.

We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and $Y_{X}$. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star formation and stellar feedback, and the addition of feedback by active galactic nuclei (AGN). We select clusters with $M_{500} > 10^{14} M_{\odot} E(z)^{-1}$, mimicking the typical selection of Sunyaev-Zeldovich samples. This permits to have a mass range large enough to enable robust fitting of the relations even at $z \sim 2$. The results of the analysis show a general agreement with observations. The values of the slope of the mass-gas mass and mass-temperature relations at $z=2$ are 10 per cent lower with respect to $z=0$ due to the applied mass selection, in the former case, and to the effect of early merger in the latter. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that cosmological studies through scaling relations should be limited to the redshift range $z=0-1$, where we find that the slope, the scatter, and the covariance matrix of the relations are stable. The scaling between mass and $Y_X$ is confirmed to be the most robust relation, being almost independent of the gas physics. At higher redshifts, the scaling relations are sensitive to the inclusion of AGNs which influences low-mass systems. The detailed study of these objects will be crucial to evaluate the AGN effect on the ICM.

We investigate the impact of galactic outflow modelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping initial conditions, simulation code and star formation (SF) model all fixed. We present simulations that are based on a version of the GADGET3 code where our sub-resolution model is coupled with an advanced implementation of Smoothed Particle Hydrodynamics that ensures a more accurate fluid sampling and an improved description of gas mixing and hydrodynamical instabilities. We quantify the strong interplay between the adopted hydrodynamic scheme and the sub-resolution model describing SF and feedback. We consider four different galactic outflow models, including the one introduced by Dalla Vecchia and Schaye (2012) and a scheme that is inspired by the Springel and Hernquist (2003) model. We find that the sub-resolution prescriptions adopted to generate galactic outflows are the main shaping factor of the stellar disc component at low redshift. The key requirement that a feedback model must have to be successful in producing a disc-dominated galaxy is the ability to regulate the high-redshift SF (responsible for the formation of the bulge component), the cosmological infall of gas from the large-scale environment, and gas fall-back within the galactic radius at low redshift, in order to avoid a too high SF rate at $z=0$.

Using a series of cosmological simulations that includes one dark-matter-only (DM-only) run, one gas cooling-star formation-supernovae feedback (CSF) run and one that additionally includes feedback from active galactic nuclei (AGNs), we classify the large-scale structures with both a velocity-shear-tensor code (Vweb) and a tidal-tensor code (Pweb). We find that the baryonic processes have almost no impact on large-scale structures -- at least not when classified using aforementioned techniques. More importantly, our results confirm that the gas component alone can be used to infer the filamentary structure of the Universe practically un-biased, which could be applied to cosmology constrains. In addition, the gas filaments are classified with its velocity (Vweb) and density (Pweb) fields, which can theoretically connect to the radio observations, such as HI surveys. This will help us to bias-freely link the radio observations with DM distributions at large scale.

Comparisons between observed and predicted strong lensing properties of galaxy clusters have been routinely used to claim either tension or consistency with $\Lambda$CDM cosmology. However, standard approaches to such cosmological tests are unable to quantify the preference for one cosmology over another. We advocate approximating the relevant Bayes factor using a marginal likelihood that is based on the following summary statistic: the posterior probability distribution function for the parameters of the scaling relation between Einstein radii and cluster mass, $\alpha$ and $\beta$. We demonstrate, for the first time, a method of estimating the marginal likelihood using the X-ray selected $z>0.5$ MACS clusters as a case in point and employing both N-body and hydrodynamic simulations of clusters. We investigate the uncertainty in this estimate and consequential ability to compare competing cosmologies, that arises from incomplete descriptions of baryonic processes, discrepancies in cluster selection criteria, redshift distribution, and dynamical state. The relation between triaxial cluster masses at various overdensities provide a promising alternative to the strong lensing test.

We perform a maximum likelihood kinematic analysis of the two dynamically relaxed galaxy clusters MACS J1206.2-0847 at $z=0.44$ and RXC J2248.7-4431 at $z=0.35$ to determine the total mass profile in modified gravity models, using a modified version of the MAMPOSSt code of Mamon, Biviano and Bou\'e. Our work is based on the kinematic and lensing mass profiles derived using the data from the Cluster Lensing And Supernova survey with Hubble (hereafter CLASH) and the spectroscopic follow-up with the Very Large Telescope (hereafter CLASH-VLT). We assume a spherical Navarro-Frenk-White (NFW hereafter) profile in order to obtain a constraint on the fifth force interaction range $\lambda$ for models in which the dependence of this parameter on the enviroment is negligible at the scale considered (i.e. $\lambda=const$) and fixing the fifth force strength to the value predicted in $f(R)$ gravity. We then use information from lensing analysis to put a prior on the other NFW free parameters. In the case of MACSJ 1206 the joint kinematic+lensing analysis leads to an upper limit on the effective interaction range $\lambda\le1.61\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$ on the marginalized distribution. For RXJ 2248 instead a possible tension with the $\Lambda$CDM model appears when adding lensing information, with a lower limit $\lambda\ge0.14\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$. This is consequence of the slight difference between the lensing and kinematic data, appearing in GR for this cluster, that could be explained in terms of modifications of gravity. We discuss the impact of systematics and the limits of our analysis as well as future improvements of the results obtained.[...]

We derive and compare the fractions of cool-core clusters in the {\em Planck} Early Sunyaev-Zel'dovich sample of 164 clusters with $z \leq 0.35$ and in a flux-limited X-ray sample of 100 clusters with $z \leq 0.30$, using {\em Chandra} observations. We use four metrics to identify cool-core clusters: 1) the concentration parameter: the ratio of the integrated emissivity profile within 0.15 $r_{500}$ to that within $r_{500}$, and 2) the ratio of the integrated emissivity profile within 40 kpc to that within 400 kpc, 3) the cuspiness of the gas density profile: the negative of the logarithmic derivative of the gas density with respect to the radius, measured at 0.04 $r_{500}$, and 4) the central gas density, measured at 0.01 $r_{500}$. We find that the sample of X-ray selected clusters, as characterized by each of these metrics, contains a significantly larger fraction of cool-core clusters compared to the sample of SZ selected clusters (44$\pm$7\% vs. 28$\pm$4\% using the concentration parameter in the 0.15--1.0 $r_{500}$ range, 61$\pm$8\% vs. 36$\pm$5\% using the concentration parameter in the 40--400 kpc range, 64$\pm$8\% vs. 38$\pm$5\% using the cuspiness, and 53$\pm$7\% vs. 39$\pm$5\% using the central gas density). Qualitatively, cool-core clusters are more X-ray luminous at fixed mass. Hence, our X-ray flux-limited sample, compared to the approximately mass-limited SZ sample, is over-represented with cool-core clusters. We describe a simple quantitative model that uses the excess luminosity of cool-core clusters compared to non-cool-core clusters at fixed mass to successfully predict the observed fraction of cool-core clusters in X-ray selected samples.

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at $z=0,0.5,1$, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods [abridged]

The forecasted accuracy of upcoming surveys of large-scale structure cannot be achieved without a proper quantification of the error induced by foreground removal (or other systematics like 0-point photometry offset). Because these errors are highly correlated on the sky, their influence is expected to be especially important at very large scales. In this work we quantify how the uncertainty in the visibility mask of a survey influences the measured power spectrum of a sample of tracers of the density field and its covariance matrix. We start from a very large set of 10,000 catalogs of dark matter halos in periodic cosmological boxes, produced with the PINOCCHIO approximate method. To make an analytic approach feasible, we assume luminosity-independent halo bias and an idealized geometry for the visibility mask. We find that the power spectrum of these biased tracers can be expressed as the sum of a cosmological term, a mask term and a term involving their convolution. The mask and convolution terms scale like $P\propto l^2\sigma_A^2$, where $\sigma_A^2$ is the variance of the uncertainty on the visibility mask. With $l=30-100$ Mpc$/h$ and $\sigma_A=5-20$\%, the mask term can be significant at $k\sim0.01-0.1\ h/$Mpc, and the convolution term can amount to $\sim 1-10$\% of the total. For the power spectrum covariance, the coupling of the convolution term with the other two gives rise to several mixed terms, that we quantify by difference using the mock catalogs. These are found to be of the same order of the mask covariance, and to introduce non-diagonal terms at large scales. Then, the power spectrum covariance matrix cannot be expressed as the sum of a cosmological and of a mask term. Our results lie down the theoretical bases to quantify the impact that uncertainties in the mask calibration have on the derivation of cosmological constraints from large spectroscopic surveys. [Abridged]

The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes included in the simulation and with the thermal properties of the core. These simulations, consisting of re-simulations of 29 Lagrangian regions performed with an upgraded version of the SPH GADGET-3 code, have been run including two different sets of baryonic physics: one accounts for radiative cooling, star formation, metal enrichment and supernova (SN) feedback, and the other one further includes the effects of feedback from active galactic nuclei (AGN). In agreement with observations, we find an anti-correlation between entropy and metallicity in cluster cores, and similar radial distributions of heavy-element abundances and abundance ratios out to large cluster-centric distances (~R180). In the outskirts, namely outside of ~0.2R180, we find a remarkably homogeneous metallicity distribution, with almost flat profiles of the elements produced by either SNIa or SNII. We investigated the origin of this phenomenon and discovered that it is due to the widespread displacement of metal-rich gas by early (z>2-3) AGN powerful bursts, acting on small high-redshift haloes. Our results also indicate that the intrinsic metallicity of the hot gas for this sample is on average consistent with no evolution between z=2 and z=0, across the entire radial range.

We analyze the radial pressure profiles, the ICM clumping factor and the Sunyaev-Zel'dovich (SZ) scaling relations of a sample of simulated galaxy clusters and groups identified in a set of hydrodynamical simulations based on an updated version of the TreePM-SPH GADGET-3 code. Three different sets of simulations are performed: the first assumes non-radiative physics, the others include, among other processes, AGN and/or stellar feedback. Our results are analyzed as a function of redshift, ICM physics, cluster mass and cluster cool-coreness or dynamical state. In general, the mean pressure profiles obtained for our sample of groups and clusters show a good agreement with X-ray and SZ observations. Simulated cool-core (CC) and non-cool-core (NCC) clusters also show a good match with real data. We obtain in all cases a small (if any) redshift evolution of the pressure profiles of massive clusters, at least back to z=1. We find that the clumpiness of gas density and pressure increases with the distance from the cluster center and with the dynamical activity. The inclusion of AGN feedback in our simulations generates values for the gas clumping ($\sqrt C_{\rho}\sim 1.2$ at $R_{200}$) in good agreement with recent observational estimates. The simulated $Y_{SZ}-M$ scaling relations are in good accordance with several observed samples, especially for massive clusters. As for the scatter of these relations, we obtain a clear dependence on the cluster dynamical state, whereas this distinction is not so evident when looking at the subsamples of CC and NCC clusters.

We present a method that extends the capabilities of the PINpointing Orbit-Crossing Collapsed HIerarchical Objects (PINOCCHIO) code, allowing it to generate accurate dark matter halo mock catalogues in cosmological models where the linear growth factor and the growth rate depend on scale. Such cosmologies comprise, among others, models with massive neutrinos and some classes of modified gravity theories. We validate the code by comparing the halo properties from PINOCCHIO against N-body simulations, focusing on cosmologies with massive neutrinos: $\nu\Lambda$CDM. We analyse the halo mass function, halo two-point correlation function, halo power spectrum and the moments of the halo density field, showing that PINOCCHIO reproduces the results from simulations with the same level of precision as the original code ($\sim5-10\%$). We demonstrate that the abundance of halos in cosmologies with massless and massive neutrinos from PINOCCHIO matches very well the outcome of simulations, and point out that PINOCCHIO can reproduce the $\Omega_\nu-\sigma_8$ degeneracy that affects the halo mass function. We show that the clustering properties of the halos from PINOCCHIO matches accurately those from simulations both in real and redshift-space, in the latter case up to $k=0.3~h~{\rm Mpc}^{-1}$. We finally point out that the first moments of the halo density field from simulations are precisely reproduced by PINOCCHIO. We emphasize that the computational time required by PINOCCHIO to generate mock halo catalogues is orders of magnitude lower than the one needed for N-body simulations. This makes this tool ideal for applications like covariance matrix studies within the standard $\Lambda$CDM model but also in cosmologies with massive neutrinos or some modified gravity theories.

Our study is meant to extend our knowledge of the galaxy color and luminosity segregation in velocity space (VCS and VLS, resp.), to clusters at intermediate and high redshift. Our sample is a collection of 41 clusters in the 0.4<~z<~1.5 redshift range, for a total of 4172 galaxies, 1674 member galaxies within 2R200 with photometric or spectroscopic information, as taken from the literature. We pay attention to perform homogeneous procedures to select cluster members, compute global cluster properties, in particular the LOS velocity dispersion sigmaV, and separate blue from red galaxies. We find evidence of VCS in clusters out to z~0.8 (at the 97%-99.99% c.l., depending on the test), in the sense that the blue galaxy population has a 10-20% larger sigmaV than the red galaxy population. Poor or no VCS is found in the High-z sample at z>=0.8. For the first time, we detect VLS in non-local clusters and confirm that VLS only affects the very luminous galaxies, with brighter galaxies having lower velocities. The threshold magnitude of VLS is ~m3+0.5, where m3 is the magnitude of the third brightest cluster galaxy, and current data suggest that the threshold value moves to fainter magnitudes at higher redshift. We also detect (marginal) evidence of VLS for blue galaxies. We conclude that the segregation effects, when their study is extended to distant clusters, can be important tracers of the galaxy evolution and cluster assembly and discuss the poor/no evidence of VCS at high redshift.

In this Letter we compare the abundance of member galaxies of a rich, nearby ($z=0.09$) galaxy cluster, Abell 2142, with that of halos of comparable virial mass extracted from sets of state-of-the-art numerical simulations, both collisionless at different resolutions and with the inclusion of baryonic physics in the form of cooling, star formation, and feedback by active galactic nuclei. We also use two semi-analytical models to account for the presence of orphan galaxies. The photometric and spectroscopic information, taken from the Sloan Digital Sky Survey Data Release 12 (SDSS DR12) database, allows us to estimate the stellar velocity dispersion of member galaxies of Abell 2142. This quantity is used as proxy for the total mass of secure cluster members and is properly compared with that of subhalos in simulations. We find that simulated halos have a statistically significant ($\gtrsim 7$ sigma confidence level) smaller amount of massive (circular velocity above $200\,{\rm km\, s^{-1}}$) subhalos, even before accounting for the possible incompleteness of observations. These results corroborate the findings from a recent strong lensing study of the Hubble Frontier Fields galaxy cluster MACS J0416 \citep{grillo2015} and suggest that the observed difference is already present at the level of dark matter (DM) subhalos and is not solved by introducing baryonic physics. A deeper understanding of this discrepancy between observations and simulations will provide valuable insights into the impact of the physical properties of DM particles and the effect of baryons on the formation and evolution of cosmological structures.

Using a suite of cosmology simulations of a sample of $> 120$ galaxy clusters with $\log(M_{DM, vir}) \le 14.5$. We compare clusters that form in purely dark matter run and their counterparts in hydro runs and investigate 4 independent parameters, that are normally used to classify dynamical state. We find that the virial ratio $\eta$ in hydro-dynamical runs is $\sim 10$ per cent lower than in the DM run, and there is no clear separation between the relaxed and unrelaxed clusters for any parameter. Further, using the velocity dispersion deviation parameter $\zeta$, which is defined as the ratio between cluster velocity dispersion $\sigma$ and the theoretical prediction $\sigma_t = \sqrt{G M_{total}/R}$, we find that there is a linear correlation between the virial ratio $\eta$ and this $\zeta$ parameter. We propose to use this $\zeta$ parameter, which can be easily derived from observed galaxy clusters, as a substitute of the $\eta$ parameter to quantify the cluster dynamical state.

Clusters of galaxies have the potential of providing powerful constraints on possible deviations from General Relativity. We use the catalogue of Sunyaev-Zel'dovich sources detected by Planck and consider a correction to the halo mass function for a f(R) class of modified gravity models, which has been recently found to reproduce well results from N-body simulations, to place constraints on the scalaron field amplitude at the present time, $f_{R}^0$. We find that applying this correction to different calibrations of the halo mass function produces upper bounds on $f_{R}^0$ tighter by more than an order of magnitude, ranging from $\log_{10}(-f_{R}^0) < -5.81$ to $\log_{10}(-f_{R}^0) < -4.40$ (95 % confidence level). This sensitivity is due to the different shape of the halo mass function, which is degenerate with the parameters used to calibrate the scaling relations between SZ observables and cluster masses. Any claim of constraints more stringent that the weaker limit above, based on cluster number counts, appear to be premature and must be supported by a careful calibration of the halo mass function and by a robust calibration of the mass scaling relations.

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.

We implement novel numerical models of AGN feedback in the SPH code GADGET-3, where the energy from a supermassive black hole (BH) is coupled to the surrounding gas in the kinetic form. Gas particles lying inside a bi-conical volume around the BH are imparted a one-time velocity (10,000 km/s) increment. We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h M_sun), which is initially evolved to form a dense cool core, having central T<10^6 K. A BH resides at the cluster center, and ejects energy. The feedback-driven fast wind undergoes shock with the slower-moving gas, which causes the imparted kinetic energy to be thermalized. Bipolar bubble-like outflows form propagating radially outward to a distance of a few 100 kpc. The radial profiles of median gas properties are influenced by BH feedback in the inner regions (r<20-50 kpc). BH kinetic feedback, with a large value of the feedback efficiency, depletes the inner cool gas and reduces the hot gas content, such that the initial cool core of the cluster is heated up within a time 1.9 Gyr, whereby the core median temperature rises to above 10^7 K, and the central entropy flattens. Our implementation of BH thermal feedback (using the same efficiency as kinetic), within the star-formation model, cannot do this heating, where the cool core remains. The inclusion of cold gas accretion in the simulations produces naturally a duty cycle of the AGN with a periodicity of 100 Myr.

In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved version of the code that allows for a better description of hydrodynamical instabilities and gas mixing processes. Evaluating the radial balance between the gravitational and hydrodynamical forces, via the gas accelerations generated, we effectively examine the level of HE in every object of the sample, its dependence on the radial distance from the center and on the classification of the cluster in terms of either cool-coreness or dynamical state. We find an average deviation of 10-20% out to the virial radius, with no evident distinction between cool-core and non-cool-core clusters. Instead, we observe a clear separation between regular and disturbed systems, with a more significant deviation from HE for the disturbed objects. The investigation of the bias between the hydrostatic estimate and the total gravitating mass indicates that, on average, this traces very well the deviation from HE, even though individual cases show a more complex picture. Typically, in the radial ranges where mass bias and deviation from HE are substantially different, the gas is characterized by a significant amount of random motions (>~30 per cent), relative to thermal ones. As a general result, the HE-deviation and mass bias, at given interesting distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We present results from multifrequency radiative hydrodynamical chemistry simulations addressing primordial star formation and related stellar feedback from various populations of stars, stellar energy distributions (SEDs) and initial mass functions. Spectra for massive stars, intermediate-mass stars and regular solar-like stars are adopted over a grid of 150 frequency bins and consistently coupled with hydrodynamics, heavy-element pollution and non-equilibrium species calculations. Powerful massive population III stars are found to be able to largely ionize H and, subsequently, He and He$^+$, causing an inversion of the equation of state and a boost of the Jeans masses in the early intergalactic medium. Radiative effects on star formation rates are between a factor of a few and 1 dex, depending on the SED. Radiative processes are responsible for gas heating and photoevaporation, although emission from soft SEDs has minor impacts. These findings have implications for cosmic gas preheating, primordial direct-collapse black holes, the build-up of "cosmic fossils" such as low-mass dwarf galaxies, the role of AGNi during reionization, the early formation of extended disks and angular-momentum catastrophe.

We present the latest version of Pinocchio, a code that generates catalogues of DM haloes in an approximate but fast way with respect to an N-body simulation. This code version extends the computation of particle and halo displacements up to 3rd-order Lagrangian Perturbation Theory (LPT), in contrast with previous versions that used Zeldovich approximation (ZA). We run Pinocchio on the same initial configuration of a reference N-body simulation, so that the comparison extends to the object-by-object level. We consider haloes at redshifts 0 and 1, using different LPT orders either for halo construction - where displacements are needed to decide particle accretion onto a halo or halo merging - or to compute halo final positions. We compare the clustering properties of Pinocchio haloes with those from the simulation by computing the power spectrum and 2-point correlation function (2PCF) in real and redshift space (monopole and quadrupole), the bispectrum and the phase difference of halo distributions. We find that 2LPT and 3LPT give noticeable improvement. 3LPT provides the best agreement with N-body when it is used to displace haloes, while 2LPT gives better results for constructing haloes. At the highest orders, linear bias is typically recovered at a few per cent level. In Fourier space and using 3LPT for halo displacements, the halo power spectrum is recovered to within 10 per cent up to $k_{max}\sim0.5\ h/$Mpc. The results presented in this paper have interesting implications for the generation of large ensemble of mock surveys aimed at accurately compute covariance matrices for clustering statistics.

We use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $\eta=\Psi/\Phi$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $\Phi$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $\Phi+\Psi$). We thus express $\eta$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $\eta(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $\eta$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analyzed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter halos monotonically increases with the halo mass and can be well described by a power-law of the form $M_{\rm HI}(M,z)\propto M^{3/4}$. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by $\sim50\%$, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within halos is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to halos not resolved in our simulations we derive astrophysical implications from the measurements of $\Omega_{\rm HI}(z)$: halos with circular velocities larger than $\sim25~{\rm km/s}$ are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of $\Omega_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations.

We investigate the galaxy population in simulated proto-cluster regions using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations. We select the most massive clusters at redshift $z=0$ from our set of simulations, and follow their main progenitors back in time. The analysis shows that proto-cluster regions are dominated by central galaxies and their number decreases with time as many become satellites, clustering around the central object. In agreement with observations, we find an increasing velocity dispersion with cosmic time, the increase being faster for satellites. The analysis shows that proto-clusters are very extended regions, $\gtrsim 20 \, Mpc$ at $z \gtrsim 1$. The fraction of galaxies in proto-cluster regions that are not progenitor of cluster galaxies varies with redshift, stellar mass and area considered. It is about 20-30 per cent for galaxies with stellar mass $\sim 10^9\,{\rm M}_{\sun}$, while negligible for the most massive galaxies considered. Nevertheless, these objects have properties similar to those of progenitors. We investigate the building-up of the passive-sequence in clusters, and find that their progenitors are on average always active at any redshift of interest of proto-clusters. The main mechanism which quenches their star formation is the removal of the hot gas reservoir at the time of accretion. The later galaxies are accreted (become satellite), and the more the cold gas available, the longer the time spent as active. Central galaxies are active over all redshift range considered, although a non-negligible fraction of them become passive at redshift $z<1$, due to strong feedback from Active Galactic Nuclei.

Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical sample of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data.

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 present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and in observations. Furthermore, we also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.

We study the characteristics of the galaxy cluster samples expected from the European Space Agency's Euclid satellite and forecast constraints on cosmological parameters describing a variety of cosmological models. The method used in this paper, based on the Fisher Matrix approach, is the same one used to provide the constraints presented in the Euclid Red Book (Laureijs et al.2011). We describe the analytical approach to compute the selection function of the photometric and spectroscopic cluster surveys. Based on the photometric selection function, we forecast the constraints on a number of cosmological parameter sets corresponding to different extensions of the standard LambdaCDM model. The dynamical evolution of dark energy will be constrained to Delta w_0=0.03 and Delta w_a=0.2 with free curvature Omega_k, resulting in a (w_0,w_a) Figure of Merit (FoM) of 291. Including the Planck CMB covariance matrix improves the constraints to Delta w_0=0.02, Delta w_a=0.07 and a FoM=802. The amplitude of primordial non-Gaussianity, parametrised by f_NL, will be constrained to \Delta f_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter \gamma, which signals deviations from GR, will be constrained to Delta \gamma=0.02, and the neutrino density parameter to Delta Omega_\nu=0.0013 (or Delta \sum m_\nu=0.01). We emphasise that knowledge of the observable--mass scaling relation will be crucial to constrain cosmological parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass calibration from weak lensing and the dynamics of cluster galaxies. This information will be further complemented by wide-area multi-wavelength external cluster surveys that will already be available when Euclid flies. [Abridged]

We compute and study the infrared and sub-mm properties of high redshift ($z \gtrsim 1$) simulated clusters and proto-clusters. The results of a large set of hydro-dynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry. Here, we have slightly generalized the code to adapt it to the present purpose. Then we have post-processed boxes of physical size 2 Mpc encompassing each of the 24 most massive clusters identified at z=0, at several redshifts between 0.5 and 3, producing IR and sub-mm mock images of these regions and SEDs of the radiation coming out from them. While this field is in its infancy from the observational point of view, rapid development is expected in the near future thanks to observations performed in the far IR and sub-mm bands. Notably, we find that in this spectral regime our prediction are little affected by the assumption required by this post-processing, and the emission is mostly powered by star formation rather than accretion onto super massive black hole (SMBH). The comparison with the little observational information currently available, highlights that the simulated cluster regions never attain the impressive star formation rates suggested by these observations. This problem becomes more intriguing taking into account that the brightest cluster galaxies (BCGs) in the same simulations turn out to be too massive. It seems that the interplay between the feedback schemes and the star formation model should be revised, possibly incorporating a positive feedback mode.

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

We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

Studying galaxy clusters through their Sunyaev-Zel'dovich (SZ) imprint on the Cosmic Microwave Background has many important advantages. The total SZ signal is an accurate and precise tracer of the total pressure in the intra-cluster medium and of cluster mass, the key observable for using clusters as cosmological probes. Band 5 observations with SKA-MID towards cluster surveys from the next generation of X-ray telescopes such as e-ROSITA and from Euclid will provide the robust mass estimates required to exploit these samples. This will be especially important for high redshift systems, arising from the SZ's unique independence to redshift. In addition, galaxy clusters are very interesting astrophysical systems in their own right, and the SKA's excellent surface brightness sensitivity down to small angular scales will allow us to explore the detailed gas physics of the intra-cluster medium.

We report the analysis of the Chandra observation of XDCP J0044.0-2033, a massive, distant (z=1.579) galaxy cluster discovered in the XDCP survey. The total exposure time of 380 ks with Chandra ACIS-S provides the deepest X-ray observation currently achieved on a massive, high redshift cluster. Extended emission from the Intra Cluster Medium (ICM) is detected at a very high significance level (S/N~20) on a circular region with a 44" radius, corresponding to $R_{ext}=375$ kpc at the cluster redshift. We perform an X-ray spectral fit of the ICM emission modeling the spectrum with a single-temperature thermal mekal model. Our analysis provides a global temperature $kT=6.7^{+1.3}_{-0.9}$ keV, and a iron abundance $Z_{Fe} = 0.41_{-0.26}^{+0.29}Z_{Fe_\odot}$ (error bars correspond to 1 $\sigma$). We fit the background-subtracted surface brightness profile with a single $\beta$-model out to 44", finding a rather flat profile with no hints of a cool core. We derive the deprojected electron density profile and compute the ICM mass within the extraction radius $R_{ext}=375$ kpc to be $M_{ICM}(r<R_{ext}) = (1.48 \pm 0.20) \times 10^{13} M_\odot$. Under the assumption of hydrostatic equilibrium and assuming isothermality within $R_{ext}$, the total mass is $M_{2500}= 1.23_{-0.27}^{+0.46} \times 10 ^{14} M_\odot$ for $R_{2500} = 240_{-20}^{+30}$ kpc. Extrapolating the profile at radii larger than the extraction radius $R_{ext}$ we find $M_{500} = 3.2_{-0.6}^{+0.9} \times 10 ^{14}M_\odot$ for $R_{500} = 562_{-37}^{+50}$ kpc. This analysis establishes the existence of virialized, massive galaxy clusters at redshift $z\sim 1.6$, paving the way to the investigation of the progenitors of the most massive clusters today. Given its mass and the XDCP survey volume, XDCP J0044.0-2033 does not create significant tension with the WMAP-7 $\Lambda$CDM cosmology.

We present results of cosmological simulations of disk galaxies carried out with the GADGET-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator (MUPPI) model. This description is based on simple multi-phase model of the interstellar medium at unresolved scales, where mass and energy flows among the components are explicitly followed by solving a system of ordinary differential equations. Thermal energy from SNe is injected into the local hot phase, so as to avoid that it is promptly radiated away. A kinetic feedback prescription generates the massive outflows needed to avoid the over-production of stars. We use two sets of zoomed-in initial conditions of isolated cosmological halos with masses (2-3) * 10^{12} Msun, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T \approx 0.2), extended stellar and gas disks, flat rotation curves and realistic values of stellar masses. Gas profiles are relatively flat, molecular gas is found to dominate at the centre of galaxies, with star formation rates following the observed Schmidt-Kennicutt relation. Stars kinematically belonging to the bulge form early, while disk stars show a clear inside-out formation pattern and mostly form after redshift z=2. However, the baryon conversion efficiencies in our simulations differ from the relation given by Moster et al. (2010) at a 3 sigma level, thus indicating that our stellar disks are still too massive for the Dark Matter halo in which they reside. Results are found to be remarkably stable against resolution. This further demonstrates the feasibility of carrying out simulations producing a realistic population of galaxies within representative cosmological volumes, at a relatively modest resolution.

(Abridged) We discuss the reliability of neutrino mass constraints, either active or sterile, from the combination of WMAP 9-year or Planck CMB data with BAO measurements from BOSS DR11, galaxy shear measurements from CFHTLenS, SDSS Ly-$\alpha$ forest constraints and galaxy cluster mass function from Chandra observations. To avoid model dependence of the constraints we perform a full likelihood analysis for all the datasets employed. As for the cluster data analysis we rely on to the most recent calibration of massive neutrino effects in the halo mass function and we explore the impact of the uncertainty in the mass bias and re-calibration of the halo mass function due to baryonic feedback processes on cosmological parameters. We find that none of the low redshift probes alone provide evidence for massive neutrinos in combination with CMB measurements, while a larger than $2\sigma$ detection of non zero neutrino mass, either active or sterile, is achieved combining cluster or shear data with CMB and BAO measurements. The preference for massive neutrino is larger in the sterile neutrino scenario, and for the combination of Planck, BAO, shear and cluster datasets we find that the vanilla $\Lambda$CDM model is rejected at more than $3\sigma$ and a sterile neutrino mass as motivated by accelerator anomaly is within the $2\sigma$ errors. Finally, results from the full data combination reflect the tension between the $\sigma_8$ constraints obtained from cluster and shear data and that inferred from Ly-$\alpha$ forest measurements; in the active neutrino scenario for both CMB datasets employed, the full data combination yields only an upper limits on $\sum m_\nu$, while assuming an extra sterile neutrino we still get preference for non-vanishing mass, $m_s^{\rm eff}=0.26^{+0.22}_{-0.24}$ eV, and dark contribution to the radiation content, $\Delta N_{\rm eff}=0.82\pm0.55$.

[Abridged.] We investigate baryon effects on the halo mass function (HMF), with emphasis on the role played by AGN feedback. Halos are identified with both Friends-of-Friends (FoF) and Spherical Overdensity (SO) algorithms. We embed the standard SO algorithm into a memory-controlled frame program and present the {\bf P}ython spher{\bf I}c{\bf A}l {\bf O}verdensity code --- {\small PIAO}. For both FoF and SO halos, the effect of AGN feedback is that of suppressing the HMFs to a level even below that of Dark Matter simulations. The ratio between the HMFs in the AGN and in the DM simulations is $\sim 0.8$ at overdensity $\Delta_c=500$, a difference that increases at higher overdensity $\Delta_c=2500$, with no significant redshift and mass dependence. A decrease of the halo masses ratio with respect to the DM case induces the decrease of the HMF in the AGN simulation. The shallower inner density profiles of halos in the AGN simulation witnesses that mass reduction is induced by the sudden displacement of gas induced by thermal AGN feedback. We provide fitting functions to describe halo mass variations at different overdensities, which can recover the HMFs with a residual random scatter $\lt 5$ per cent for halo masses larger than $10^{13} ~h^{-1}{\rm M_\odot}$.

Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10% to 30%. The largest bias originates by both violation of hydrostatic equilibrium and an additional temperature bias caused by inhomogeneities in the X-ray emitting intra-cluster medium (ICM). To elucidate on this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle-hydrodynamics (SPH) and adaptive-mesh-refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in no-radiative simulations, while it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with the observed one in the central regions of clusters. Around R_500 the temperature inhomogeneities of the SPH simulations can generate twice the typical hydrostatic-equilibrium mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.

(Abridged) We present a spectral analysis of a new, flux-limited sample of 72 X-ray selected clusters of galaxies identified with the X-ray Telescope (XRT) on board the Swift satellite down to a flux limit of ~10-14 erg/s/cm2 (SWXCS, Tundo et al. 2012). We carry out a detailed X-ray spectral analysis with the twofold aim of measuring redshifts and characterizing the properties of the Intra-Cluster Medium (ICM). Optical counterparts and spectroscopic or photometric redshifts are obtained with a cross-correlation with NED. Additional photometric redshifts are computed with a dedicated follow-up program with the TNG and a cross-correlation with the SDSS. We also detect the iron emission lines in 35% of the sample, and hence obtain a robust measure of the X-ray redshift zX. We use zX whenever the optical redshift is not available. Finally, for all the sources with measured redshift, background-subtracted spectra are fitted with a mekal model. We perform extensive spectral simulations to derive an empirical formula to account for fitting bias. The bias-corrected values are then used to investigate the scaling properties of the X-ray observables. Overall, we are able to characterize the ICM of 46 sources. The sample is mostly constituted by clusters with temperatures between 3 and 10 keV, plus 14 low-mass clusters and groups with temperatures below 3 keV. The redshift distribution peaks around z~0.25 and extends up to z~1, with 60% of the sample at 0.1<z<0.4. We derive the Luminosity-Temperature relation for these 46 sources, finding good agreement with previous studies. The quality of the SWXCS sample is comparable to other samples available in the literature and obtained with much larger X-ray telescopes. Our results have interesting implications for the design of future X-ray survey telescopes, characterised by good-quality PSF over the entire field of view and low background.

In this study, we present a detailed, statistical analysis of black hole growth and the evolution of active galactic nuclei (AGN) using cosmological hydrodynamic simulations run down to $z=0$. The simulations self-consistently follow radiative cooling, star formation, metal enrichment, black hole growth and associated feedback processes from both supernovae typeII/Ia and AGN. We consider two simulation runs, one with a large co-moving volume of $(500\ \mathrm{Mpc})^3$ and one with a smaller volume of $(68\ \mathrm{Mpc})^3$ but with a by a factor of almost 20 higher mass resolution. Consistently with previous results, our simulations can widely match observed black hole properties of the local Universe. Furthermore, our simulations can successfully reproduce the evolution of the bolometric AGN luminosity function for both the low-luminosity and the high-luminosity end up to $z=3.0$. In addition, the smaller but higher resolution run is able to match the observational data of the low bolometric luminosity end at higher redshifts $z=3-4$. We also perform a direct comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for a torus-level obscuration, and find a similarly good agreement. These results nicely demonstrate that the observed "anti-hierarchical" trend in the AGN number density evolution (i.e. the number densities of luminous AGN peak at higher redshifts than those of faint AGN) is self-consistently predicted by our simulations. Implications of this downsizing behaviour on active black holes, their masses and Eddington-ratios are discussed. Overall, the downsizing behaviour in the AGN number density as a function of redshift can be mainly attributed to the evolution of the gas density in the resolved vicinity of a (massive) black hole. (shortened)

We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy.

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2\rho)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

We study the role of feedback from supernovae and black holes in the evolution of the star formation rate function (SFRF) of z~4-7 galaxies. We use a new set of cosmological hydrodynamic simulations, ANGUS (AustraliaN GADGET-3 early Universe Simulations), run with a modified and improved version of the parallel TreePM-smoothed particle hydrodynamics code GADGET-3 called P-GADGET3(XXL), that includes a self-consistent implementation of stellar evolution and metal enrichment. In our simulations both Supernova (SN) driven galactic winds and Active Galactic Nuclei (AGN) act simultaneously in a complex interplay. The SFRF is insensitive to feedback prescription at z>5, meaning that it cannot be used to discriminate between feedback models during reionisation. However, the SFRF is sensitive to the details of feedback prescription at lower redshift. By exploring different SN driven wind velocities and regimes for the AGN feedback, we find that the key factor for reproducing the observed SFRFs is a combination of "strong" SN winds and early AGN feedback in low mass galaxies. Conversely, we show that the choice of initial mass function and inclusion of metal cooling have less impact on the evolution of the SFRF. When variable winds are considered, we find that a non-aggressive wind scaling is needed to reproduce the SFRFs at z>4. Otherwise, the amount of objects with low SFRs is greatly suppressed and at the same time winds are not effective enough in the most massive systems.

We use a suite of N-body simulations that incorporate massive neutrinos as an extra-set of particles to investigate their effect on the halo mass function. We show that for cosmologies with massive neutrinos the mass function of dark matter haloes selected using the spherical overdensity (SO) criterion is well reproduced by the fitting formula of Tinker et al. (2008) once the cold dark matter power spectrum is considered instead of the total matter power, as it is usually done. The differences between the two implementations, i.e. using $P_{\rm cdm}(k)$ instead of $P_{\rm m}(k)$, are more pronounced for large values of the neutrino masses and in the high end of the halo mass function: in particular, the number of massive haloes is higher when $P_{\rm cdm}(k)$ is considered rather than $P_{\rm m}(k)$. As a quantitative application of our findings we consider a Planck-like SZ-clusters survey and show that the differences in predicted number counts can be as large as $30\%$ for $\sum m_\nu = 0.4$ eV. Finally, we use the Planck-SZ clusters sample, with an approximate likelihood calculation, to derive Planck-like constraints on cosmological parameters. We find that, in a massive neutrino cosmology, our correction to the halo mass function produces a shift in the $\sigma_8(\Omega_{\rm m}/0.27)^\gamma$ relation which can be quantified as $\Delta \gamma \sim 0.05$ and $\Delta \gamma \sim 0.14$ assuming one ($N_\nu=1$) or three ($N_\nu=3$) degenerate massive neutrino, respectively. The shift results in a lower mean value of $\sigma_8$ with $\Delta \sigma_8 = 0.01$ for $N_\nu=1$ and $\Delta \sigma_8 = 0.02$ for $N_\nu=3$, respectively. Such difference, in a cosmology with massive neutrinos, would increase the tension between cluster abundance and Planck CMB measurements.

[Abridged] In this paper, we carry out a detailed analysis of the performance of two different methods to identify the diffuse stellar light in cosmological hydrodynamical simulations of galaxy clusters. One method is based on a dynamical analysis of the stellar component. The second method is closer to techniques commonly employed in observational studies. Both the dynamical method and the method based on the surface brightness limit criterion are applied to the same set of hydrodynamical simulations for a large sample about 80 galaxy clusters. We find significant differences between the ICL and DSC fractions computed with the two corresponding methods, which amounts to about a factor of two for the AGN simulations, and a factor of four for the CSF set. We also find that the inclusion of AGN feedback boosts the DSC and ICL fractions by a factor of 1.5-2, respectively, while leaving the BCG+ICL and BCG+DSC mass fraction almost unchanged. The sum of the BCG and DSC mass stellar mass fraction is found to decrease from ~80 per cent in galaxy groups to ~60 per cent in rich clusters, thus in excess of what found from observational analysis. We identify the average surface brightness limits that yields the ICL fraction from the SBL method close to the DSC fraction from the dynamical method. These surface brightness limits turn out to be brighter in the CSF than in the AGN simulations. This is consistent with the finding that AGN feedback makes BCGs to be less massive and with shallower density profiles than in the CSF simulations. The BCG stellar component, as identified by both methods, are slightly older and more metal-rich than the stars in the diffuse component.

We study the formation of the Intra-Cluster Light (ICL) using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations of groups and clusters. We assume that the ICL forms by (1) stellar stripping of satellite galaxies and (2) relaxation processes that take place during galaxy mergers. The fraction of ICL in groups and clusters predicted by our models ranges between 10 and 40 per cent, with a large halo-to-halo scatter and no halo mass dependence. We note, however, that our predicted ICL fractions depend on the resolution: for a set of simulations with particle mass one order of magnitude larger than that adopted in the high resolution runs used in our study, we find that the predicted ICL fractions are ~30-40 per cent larger than those found in the high resolution runs. On cluster scale, large part of the scatter is due to a range of dynamical histories, while on smaller scale it is driven by individual accretion events and stripping of very massive satellites, $M_{*} \gtrsim 10^{10.5} M_{\odot}$, that we find to be the major contributors to the ICL. The ICL in our models forms very late (below $z\sim 1$), and a fraction varying between 5 and 25 per cent of it has been accreted during the hierarchical growth of haloes. In agreement with recent observational measurements, we find the ICL to be made of stars covering a relatively large range of metallicity, with the bulk of them being sub-solar.

We perform a forecast analysis on how well a Euclid-like photometric galaxy cluster survey will constrain the total neutrino mass and effective number of neutrino species. We base our analysis on the Monte Carlo Markov Chains technique by combining information from cluster number counts and cluster power spectrum. We find that combining cluster data with CMB measurements from Planck improves by more than an order of magnitude the constraint on neutrino masses compared to each probe used independently. For the LCDM+m_nu model the 2 sigma upper limit on total neutrino mass shifts from M_nu < 0.35 eV using cluster data alone to M_nu < 0.031 eV when combined with CMB data. When a non-standard model with N_eff number of neutrino species is considered, we estimate N_eff<3.14 (95% CL), while the bounds on neutrino mass are relaxed to M_nu < 0.040 eV. This accuracy would be sufficient for a 2 sigma detection of neutrino mass even in the minimal normal hierarchy scenario. We also consider scenarios with a constant dark energy equation of state and a non-vanishing curvature. When these models are considered the error on M_nu is only slightly affected, while there is a larger impact of the order of ~ 15 % and ~ 20% respectively on the 2 sigma error bar of N_eff with respect to the standard case. We also treat the LCDM+m_nu+N_eff case with free nuisance parameters, which parameterize the uncertainties on the cluster mass determination. In this case, the upper bounds on M_nu are relaxed by a factor larger than two, M_nu < 0.083 eV (95% CL), hence compromising the possibility of detecting the total neutrino mass with good significance. We thus confirm the potential that a large optical/near-IR cluster survey, like that to be carried out by Euclid, could have in constraining neutrino properties [abridged].

We present an analysis of the properties of the ICM in an extended set of cosmological hydrodynamical simulations of galaxy clusters and groups performed with the TreePM+SPH GADGET-3 code. Besides a set of non-radiative simulations, we carried out two sets of simulations including radiative cooling, star formation, metal enrichment and feedback from supernovae, one of which also accounts for the effect of feedback from AGN resulting from gas accretion onto super-massive black holes. These simulations are analysed with the aim of studying the relative role played by SN and AGN feedback on the general properties of the diffuse hot baryons in galaxy clusters and groups: scaling relations, temperature, entropy and pressure radial profiles, and ICM chemical enrichment. We find that simulations including AGN feedback produce scaling relations that are in good agreement with X-ray observations at all mass scales. However, our simulations are not able to account for the observed diversity between CC and NCC clusters: unlike for observations, we find that temperature and entropy profiles of relaxed and unrelaxed clusters are quite similar and resemble more the observed behaviour of NCC clusters. As for the pattern of metal enrichment, we find that an enhanced level of iron abundance is produced by AGN feedback with respect to the case of purely SN feedback. As a result, while simulations including AGN produce values of iron abundance in groups in agreement with observations, they over-enrich the ICM in massive clusters. The efficiency of AGN feedback in displacing enriched gas from halos into the inter-galactic medium at high redshift also creates a widespread enrichment in the outskirts of clusters and produces profiles of iron abundance whose slope is in better agreement with observations.

This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.

(Abridged) We investigate two modes of coupling the feedback energy from a central AGN to the neighboring gas in galaxy simulations: kinetic - velocity boost, and thermal - heating. We formulate kinetic feedback models for energy-driven wind (EDW) and momentum-driven wind (MDW), using two free parameters: feedback efficiency epsilon_f, and AGN wind velocity v_w. A novel numerical algorithm is implemented in the SPH code GADGET-3, to prevent the expansion of a hole in the gas distribution around the BH. We perform simulations of isolated evolution and merger of disk galaxies, of Milky-Way mass as well as lower and higher masses. We find that in the isolated galaxy BH kinetic feedback generates intermittent bipolar jet-like gas outflows. We infer that current prescriptions for BH subgrid physics in galaxy simulations can grow the BH to observed values even in an isolated disk galaxy. The BH growth is enhanced in a galaxy merger. Comparing the [M_BH - sigma_star] relation obtained in our simulations with observational data, we conclude that it is possible to find parameter sets for a fit in all the models, except for the case with MDW feedback in a galaxy merger, in which the BH is always too massive. The BH thermal feedback implementation of Springel, Di Matteo & Hernquist (2005) within the multiphase star-formation model is found to have negligible impact on gas properties; and the effect claimed in all previous studies is attributed to gas depletion around the BH by the creation of an artificial hole. The BH mass accretion rate in our simulations exhibit heavy fluctuations. The star formation rate is quenched with feedback by removal of gas. The CGM gas at galactocentric distances (20 - 100)/h kpc are found to give the best metallicity observational diagnostic to distinguish between BH models.

We analyze the basic properties of Brightest Cluster Galaxies (BCGs) produced by state of the art cosmological zoom-in hydrodynamical simulations. These simulations have been run with different sub-grid physics included. Here we focus on the results obtained with and without the inclusion of the prescriptions for supermassive black hole (SMBH) growth and of the ensuing Active Galactic Nuclei (AGN) feedback. The latter process goes in the right direction of decreasing significantly the overall formation of stars. However, BCGs end up still containing too much stellar mass, a problem that increases with halo mass, and having an unsatisfactory structure. This is in the sense that their effective radii are too large, and that their density profiles feature a flattening on scales much larger than observed. We also find that our model of thermal AGN feedback has very little effect on the stellar velocity dispersions, which turn out to be very large. Taken together, these problems, which to some extent can be recognized also in other numerical studies typically dealing with smaller halo masses, indicate that on one hand present day sub-resolution models of AGN feedback are not effective enough in diminishing the global formation of stars in the most massive galaxies, but on the other hand they are relatively too effective in their centers. It is likely that a form of feedback generating large scale gas outflows from BCGs precursors, and a more widespread effect over the galaxy volume, can alleviate these difficulties.

[Abridged] In the past 15 years, the concentration-mass relation has been investigated diffusely in theoretical studies. On the other hand, only recently has this relation been derived from X-ray observations. When that happened, the results caused a certain level of concern: the X-ray normalizations and slopes were found significantly dissimilar from those predicted by theory. We analyzed 52 objects, simulated each time with different physical recipes for the baryonic component, as well as 60 synthetic X-ray images, to determine if these discrepancies are real or artificial. In particular, we investigate how the simulated concentration-mass relation depends (1) on the radial range used to derive the concentration, (2) on the presence of baryons in the simulations, and on the prescription used to reproduce the gas. Finally, we evaluate (3) how the results differ when adopting an X-ray approach for the analysis and (4) how the selection functions based on X-ray luminosity can impact the results. All effects studied go in the direction of alleviating the discrepancy between observations and simulations, although with different significance: while the fitting radial range and the baryonic component play only a minor role, the X-ray approach and selection function have profound repercussion on the resulting concentration-mass relation.

Major astrophysical questions related to the formation and evolution of structures, and more specifically of galaxy groups and clusters, will still be open in the coming decade and beyond: what is the interplay of galaxy, supermassive black hole, and intergalactic gas evolution in the most massive objects in the Universe - galaxy groups and clusters? What are the processes driving the evolution of chemical enrichment of the hot diffuse gas in large-scale structures? How and when did the first galaxy groups in the Universe, massive enough to bind more than 10^7 K gas, form? Focussing on the period when groups and clusters assembled (0.5<z<2.5), we show that, due to the continuum and line emission of this hot intergalactic gas at X-ray wavelengths, Athena+, combining high sensitivity with excellent spectral and spatial resolution, will deliver breakthrough observations in view of the aforementioned issues. Indeed, the physical and chemical properties of the hot intra-cluster gas, and their evolution across time, are a key to understand the co-evolution of galaxy and supermassive black hole within their environments.

As the nodes of the cosmic web, clusters of galaxies trace the large-scale distribution of matter in the Universe. They are thus privileged sites in which to investigate the complex physics of structure formation. However, the complete story of how these structures grow, and how they dissipate the gravitational and non-thermal components of their energy budget over cosmic time, is still beyond our grasp. Fundamental questions such as How do hot diffuse baryons accrete and dynamically evolve in dark matter potentials? How and when was the energy that we observe in the ICM generated and distributed? Where and when are heavy elements produced and how are they circulated? are still unanswered. Most of the cluster baryons exists in the form of a diffuse, hot, metal-enriched plasma that radiates primarily in the X-ray band (the intracluster medium, ICM), allowing the X-ray observations of the evolving cluster population to provide a unique opportunity to address these topics. Athena+ with its large collecting area and unprecedented combination of high spectral and angular resolution offers the only way to make major advances in answering these questions. Athena+ will show how the baryonic gas evolves in the dark matter potential wells by studying the motions and turbulence in the ICM. Athena+ will be able to resolve the accreting region both spatially and spectroscopically, probing the true nature and physical state of the X-ray emitting plasma. Athena+ has the capabilities to permit a definitive understanding of the formation and evolution of large-scale cosmic structure through the study of the cluster population.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

We present a new parallel implementation of the PINpointing Orbit Crossing-Collapsed HIerarchical Objects (PINOCCHIO) algorithm, a quick tool, based on Lagrangian Perturbation Theory, for the hierarchical build-up of Dark Matter halos in cosmological volumes. To assess its ability to predict halo correlations on large scales, we compare its results with those of an N-body simulation of a 3 Gpc/h box sampled with 2048^3 particles taken from the MICE suite, matching the same seeds for the initial conditions. Thanks to the FFTW libraries and to the relatively simple design, the code shows very good scaling properties. The CPU time required by PINOCCHIO is a tiny fraction (~1/2000) of that required by the MICE simulation. Varying some of PINOCCHIO numerical parameters allows one to produce a universal mass function that lies in the range allowed by published fits, although it underestimates the MICE mass function of FoF halos in the high mass tail. We compare the matter-halo and the halo-halo power spectra with those of the MICE simulation and find that these 2-point statistics are well recovered on large scales. In particular, when catalogs are matched in number density, agreement within ten per cent is achieved for the halo power spectrum. At scales k>0.1 h/Mpc, the inaccuracy of the Zel'dovich approximation in locating halo positions causes an underestimate of the power spectrum that can be modeled as a Gaussian factor with a damping scale of d=3 Mpc/h at z=0, decreasing at higher redshift. Finally, a remarkable match is obtained for the reduced halo bispectrum, showing a good description of nonlinear halo bias. Our results demonstrate the potential of PINOCCHIO as an accurate and flexible tool for generating large ensembles of mock galaxy surveys, with interesting applications for the analysis of large galaxy redshift surveys.

We compare three analytical prescriptions for merger times available from the literature to simulations of isolated mergers. We probe three different redshifts, and several halo concentrations, mass ratios, orbital circularities and orbital energies of the satellite. We find that prescriptions available in the literature significantly under-predict long timescales for mergers at high redshift. We argue that these results have not been highlighted previously either because the evolution of halo concentration of satellite galaxies has been neglected (in previous isolated merger simulations), or because long merger times and mergers with high initial orbital circularities are under-represented (for prescriptions based on cosmological simulations). Motivated by the evolution of halo concentration at fixed mass, an explicit dependence on redshift added as t_merger,modified(z) = (1+z)^0.44 t_merger to the prescription based on isolated mergers gives a significant improvement in the predicted merger times up to ~20 t_dyn in the redshift range 0<z<2. When this modified prescription is used to compute galaxy stellar mass functions, we find that it leads up to a 25 per cent increase in the number of low mass galaxies surviving at z=0, and a 10 per cent increase for more massive galaxies. This worsen the known over-prediction in the number of low mass galaxies by hierarchical models of galaxy formation.

The injection of turbulence in the circum-galactic medium at redshift z = 2 is investigated using the mesh-based hydrodynamic code Enzo and a subgrid-scale (SGS) model for unresolved turbulence. Radiative cooling and heating by a uniform Ultraviolet (UV) background are included in our runs and compared with the effect of turbulence modelling. Mechanisms of gas exchange between galaxies and the surrounding medium, as well as metal enrichment, are not taken into account, and turbulence is here driven solely by structure formation (mergers and shocks). We find that turbulence, both at resolved and SGS scales, impacts mostly the warm-hot intergalactic medium (WHIM), with temperature between 10^5 and 10^7 K, mainly located around collapsed and shock heated structures, and in filaments. Typical values of the ratio of turbulent to thermal pressure is 0.1 in the WHIM, corresponding to a volume-weighted average of the SGS turbulent to thermal Doppler broadening b_t / b_therm = 0.26, on length scales below the grid resolution of 25 kpc/h. In the diffuse intergalactic medium (IGM), defined in a range of baryon overdensity \delta\ between 1 and 50, the importance of turbulence is smaller, but grows as a function of gas density, and the Doppler broadening ratio is fitted by the function b_t / b_therm = 0.023 \delta^{0.58}.

Using a set of hydrodynamical simulations of 62 galaxy clusters and groups we study the ICM of inhomogeneities, focusing on the ones on the large scale that, unlike clumps, are the most difficult to identify. To this purpose we introduce the concept of residual clumpiness, C_R, that quantifies the large-scale inhomogeneity of the ICM. After showing that this quantity can be robustly defined for relaxed systems, we characterize how it varies with radius, mass and dynamical state of the halo. Most importantly, we observe that it introduces an overestimate in the determination of the density profile from the X-ray emission, which translates into a systematic overestimate of 6 (12)% in the measurement of M_gas at R_200 for our relaxed (perturbed) cluster sample. At the same time, the increase of C_R with radius introduces also a ~2% systematic underestimate in the measurement of the hydrostatic-equilibrium mass (M_he), which adds to the previous one generating a systematic ~8.5% overestimate in f_gas in our relaxed sample. Since the residual clumpiness of the ICM is not directly observable, we study its correlation with the azimuthal scatter in the X-ray surface brightness of the halo and in the y-parameter profiles. We find that their correlation is highly significant (r_S = 0.6-0.7), allowing to define the azimuthal scatter measured in the X-ray surface brightness profile and in the y-parameter as robust proxies of C_R. After providing a function that connects the two quantities, we obtain that correcting the observed gas density profiles using the azimuthal scatter eliminates the bias in the measurement of M_gas for relaxed objects, which becomes (0+/-2)% up to 2R_200, and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the systematics on the measurements of M_he and f_gas, although a significant halo to halo scatter remains. (abridged)

We present the first results from our pilot 500 ks Chandra-LETG Large Program observation of the soft X-ray brightest source in the z>=0.4 sky, the blazar 1ES 1553+113, aimed to secure the first uncontroversial detections of the missing baryons in the X-rays. We identify a total of 11 possible absorption lines, with single-line statistical significances between 2.2-4.1\sigma. Six of these lines are detected at high single-line statistical significance (3.6 <= sigma <= 4.1), while the remaining five are regarded as marginal detections in association with either other X-ray lines detected at higher significance and/or Far-Ultraviolet (FUV) signposts. In particular, five of these possible intervening absorption lines, are identified as CV and CVI K\alpha absorbers belonging to three WHIM systems at z_X = 0.312, z_X = 0.237 and <z_X> = 0.133, which also produce broad HI (and OVI for the z_X = 0.312 system) absorption in the FUV. For two of these systems (z_X = 0.312 and 0.237), the Chandra X-ray data led the a-posteriori discovery of physically consistent broad HI associations in the FUV, so confirming the power of the X-ray-FUV synergy for WHIM studies. The true statistical significances of these three X-ray absorption systems, after properly accounting for the number of redshift trials, are 5.8 sigma (z_X = 0.312; 6.3 sigma if the low-significance OV and CV K-beta associations are considered), 3.9 sigma (z_X = 0.237), and 3.8 sigma (\langle z_X \rangle = 0.133), respectively.

We carry out an analysis of a set of cosmological SPH hydrodynamical simulations of galaxy clusters and groups aimed at studying the total baryon budget in clusters, and how this budget is shared between the hot diffuse component and the stellar component. Using the TreePM+SPH GADGET-3 code, we carried out one set of non-radiative simulations, and two sets of simulations including radiative cooling, star formation and feedback from supernovae (SN), one of which also accounting for the effect of feedback from active galactic nuclei (AGN). The analysis is carried out with the twofold aim of studying the implication of stellar and hot gas content on the relative role played by SN and AGN feedback, and to calibrate the cluster baryon fraction and its evolution as a cosmological tool. We find that both radiative simulation sets predict a trend of stellar mass fraction with cluster mass that tends to be weaker than the observed one. However this tension depends on the particular set of observational data considered. Including the effect of AGN feedback alleviates this tension on the stellar mass and predicts values of the hot gas mass fraction and total baryon fraction to be in closer agreement with observational results. We further compute the ratio between the cluster baryon content and the cosmic baryon fraction, Y_b, as a function of cluster-centric radius and redshift. At R_500 we find for massive clusters with M_500>2\times10^{14} h^{-1} M_sun that Y_b is nearly independent of the physical processes included and characterized by a negligible redshift evolution: Y_{b,500}=0.85+/-0.03 with the error accounting for the intrinsic r.m.s. scatter within the set of simulated clusters. At smaller radii, R_2500, the typical value of Y_b slightly decreases, by an amount that depends on the physics included in the simulations, while its scatter increases by about a factor of two.

[Abridged] We present an analysis of the relation between the masses of cluster- and group-sized halos, extracted from $\Lambda$CDM cosmological N-body and hydrodynamic simulations, and their velocity dispersions, at different redshifts from $z=2$ to $z=0$. The main aim of this analysis is to understand how the implementation of baryonic physics in simulations affects such relation, i.e. to what extent the use of the velocity dispersion as a proxy for cluster mass determination is hampered by the imperfect knowledge of the baryonic physics. In our analysis we use several sets of simulations with different physics implemented. Velocity dispersions are determined using three different tracers, DM particles, subhalos, and galaxies. We confirm that DM particles trace a relation that is fully consistent with the theoretical expectations based on the virial theorem and with previous results presented in the literature. On the other hand, subhalos and galaxies trace steeper relations, and with larger values of the normalization. Such relations imply that galaxies and subhalos have a $\sim10$ per cent velocity bias relative to the DM particles, which can be either positive or negative, depending on halo mass, redshift and physics implemented in the simulation. We explain these differences as due to dynamical processes, namely dynamical friction and tidal disruption, acting on substructures and galaxies, but not on DM particles. These processes appear to be more or less effective, depending on the halo masses and the importance of baryon cooling, and may create a non-trivial dependence of the velocity bias and the $\soneD$--$\Mtwo$ relation on the tracer, the halo mass and its redshift. These results are relevant in view of the application of velocity dispersion as a proxy for cluster masses in ongoing and future large redshift surveys.

(Abridged) We explore new observationally-constrained sub-resolution models of galactic outflows and investigate their impact on the circumgalactic medium (CGM) over redshifts z = 2 - 4. We perform cosmological hydrodynamic simulations, including star formation, chemical enrichment, and four cases of SNe-driven outflows: no wind (NW), an energy-driven constant velocity wind (CW), a radially varying wind (RVWa) where the outflow velocity has a positive correlation with galactocentric distance (r), and a RVW with additional dependence on halo mass (RVWb). Overall, we find that the outflows expel metal-enriched gas away from galaxies, significantly quench star formation, and enrich the CGM. At z = 2, the radial profiles of gas properties around galaxy centers are most sensitive to the choice of the wind model for halo masses (10^9 - 10^11) M_sun. We infer that the RVWb model is similar to the NW case, except that it substantially enriches the CGM: the carbon metallicity (Z_C) is 10 times higher in RVWb than in NW at r > R_200; and the warm gas of 10^4 - 10^5 K and delta < 100 is enriched to 50 times higher than in NW. We also find that the impact of models CW and RVWa are similar, with the following differences. RVWa causes a greater suppression of star formation rate at z < 5, and has a higher fraction of low-density (delta < 10), warm-hot (10^4 - 10^6 K) gas than in CW. Outflows in CW produce a higher and earlier enrichment of some IGM phases than RVWa. By visual inspection, we note that the RVWa model shows more pronounced bipolar outflows and galactic disks. We present fitting formulae for [Z_C-delta] and [Z_C-r], also for the abundance of CIV as a function of r. We predict observational diagnostics to distinguish between different outflow scenarios: Z_C of the CGM gas at r = (30 - 300) kpc/h comoving, and CIV fraction of the inner gas at r < (4 - 5) kpc/h comoving.

(Abridged) We report the discovery of CXO J1415.2+3610, a distant (z~1.5) galaxy cluster serendipitously detected in a deep, high-resolution Chandra observation targeted to study the cluster WARP J1415.1+3612 at z=1.03. This is the highest-z cluster discovered with Chandra so far. Moreover, the total exposure time of 280 ks with ACIS-S provides the deepest X-ray observation currently achieved on a cluster at z>1.5. We perform an X-ray spectral fit of the extended emission of the intracluster medium (ICM) with XSPEC, and we detect at a 99.5% confidence level the rest frame 6.7-6.9 keV Iron K_\alpha line complex, from which we obtain z_X=1.46\pm0.025. The analysis of the z-3.6\mu m color-magnitude diagram shows a well defined sequence of red galaxies within 1' from the cluster X-ray emission peak with a color range [5 < z-3.6 \mu m < 6]. The photometric redshift obtained by spectral energy distribution (SED) fitting is z_phot=1.47\pm 0.25. After fixing the redshift to z=1.46, we perform the final spectral analysis and measure the average gas temperature with a 20% error, kT=5.8^{+1.2}_{-1.0} keV, and the Fe abundance Z_Fe = 1.3_{-0.5}^{+0.8}Z_\odot. We fit the background subtracted surface brightness with a single beta--model out to 35" and derive the deprojected electron density profile. The ICM mass is 1.09_{-0.2}^{+0.3}\times 10^{13} M_\odot within 300 kpc. The total mass is M_{2500}= 8.6_{-1.7}^{+2.1} \times 10 ^{13} M_\odot for R_{2500}=(220\pm 55) kpc. Extrapolating the profile at larger radii we find M_{500}= 2.1_{-0.5}^{+0.7} \times 10 ^{14} M_\odot for R_{500} = 510_{-50}^{+55}$ kpc. This analysis establishes CXOJ1415.2+3610 as one of the best characterized distant galaxy clusters based on X-ray data alone.

Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3 x R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from XMM-Newton in the region in which the profiles overlap (i.e., [0.1-1]R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is slightly flatter than most predictions from numerical simulations. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0, c500, gamma, alpha, beta] = [6.41, 1.81, 0.31, 1.33, 4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3 x R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.

The observed abundance of giant arcs produced by galaxy cluster lenses and the measured Einstein radii have presented a source of tension for LCDM. Previous cosmological tests for high-redshift clusters (z>0.5) have suffered from small number statistics in the simulated sample and the implementation of baryonic physics is likely to affect the outcome. We analyse zoomed-in simulations of a fairly large sample of cluster-sized objects, with Mvir > 3x10^14 Msun/h, identified at z=0.25 and z=0.5, for a concordance LCDM cosmology. We start with dark matter only simulations, and then add gas hydrodynamics, with different treatments of baryonic processes: non-radiative cooling, radiative cooling with star formation and galactic winds powered by supernova explosions, and finally including the effect of AGN feedback. We find that the addition of gas in non-radiative simulations does not change the strong lensing predictions significantly, but gas cooling and star formation together significantly increase the number of expected giant arcs and the Einstein radii, particularly for lower redshift clusters and lower source redshifts. Further inclusion of AGN feedback reduces the predicted strong lensing efficiencies such that the lensing probability distributions becomes closer to those obtained for simulations including only dark matter. Our results indicate that the inclusion of baryonic physics in simulations will not solve the arc-statistics problem at low redshifts, when the physical processes included provide a realistic description of cooling in the central regions of galaxy clusters. [Abridged]

We present the final results from the XMM-Newton validation follow-up of new Planck galaxy cluster candidates. We observed 15 new candidates, detected with signal-to-noise ratios between 4.0 and 6.1 in the 15.5-month nominal Planck survey. The candidates were selected using ancillary data flags derived from the ROSAT All Sky Survey (RASS) and Digitized Sky Survey all-sky maps, with the aim of pushing into the low SZ flux, high-z regime and testing RASS flags as indicators of candidate reliability. 14 new clusters were detected by XMM, including 2 double systems. Redshifts lie in the range 0.2 to 0.9, with 6 clusters at z>0.5. Estimated M500 range from 2.5 10^14 to 8 10^14 Msun. We discuss our results in the context of the full XMM validation programme, in which 51 new clusters have been detected. This includes 4 double and 2 triple systems, some of which are chance projections on the sky of clusters at different z. We find that association with a RASS-BSC source is a robust indicator of the reliability of a candidate, whereas association with a FSC source does not guarantee that the SZ candidate is a bona fide cluster. Nevertheless, most Planck clusters appear in RASS maps, with a significance greater than 2 sigma being a good indication that the candidate is a real cluster. The full sample gives a Planck sensitivity threshold of Y500 ~ 4 10^-4 arcmin^2, with indication for Malmquist bias in the YX-Y500 relation below this level. The corresponding mass threshold depends on z. Systems with M500 > 5 10^14 Msun at z > 0.5 are easily detectable with Planck. The newly-detected clusters follow the YX-Y500 relation derived from X-ray selected samples. Compared to X-ray selected clusters, the new SZ clusters have a lower X-ray luminosity on average for their mass. There is no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. (abridged)

In this Letter, we analyse the predicted physical properties of massive galaxies, in the framework of recent semi-analytic models of galaxy formation. All models considered account for winds driven by supernovae explosions and suppression of gas condensation at the centre of relatively massive haloes by active galactic nuclei (AGN). We show that, while these models successfully reproduce the old stellar populations observed for massive galaxies, they fail in reproducing their observed chemical abundances. This problem is alleviate but still present if AGN feedback is completely switched off. Moreover, in this case, model predictions fail in accounting for the old stellar ages of massive galaxies. We argue that the difficulty of semi-analytical models in simultaneously reproducing the observed ages and metallicities of massive galaxies, signals a fundamental problem with the schemes that are currently adopted to model star formation, feedback, and related recycling of gas and metals.

We present forecasts on the capability of future wide-area high-sensitivity X-ray surveys of galaxy clusters to yield constraints on the parameters defining the Dark Energy (DE) equation of state (EoS). Our analysis is carried out for future X-ray surveys which have enough sensitivity to provide accurate measurements of X-ray mass proxies and Fe-line based redshifts for about 2x10^4 clusters. We base our analysis on the Fisher Matrix formalism, by combining information on the cluster number counts and power spectrum, also including, for the first time in the analysis of the large scale cluster distribution, the effect of linear redshift-space distortions (RSDs). This study is performed with the main purpose of dissecting the cosmological information provided by geometrical and growth tests, which are both included in the analysis of number counts and clustering of galaxy clusters. We compare cosmological constraints obtained by assuming different levels of prior knowledge of the parameters which define the observable-mass X-ray relation. This comparison further demonstrates the fundamental importance of having a well calibrated observable-mass relation and, most importantly, its redshift evolution. Such a calibration can be achieved only by having at least $\sim 10^3$ net photon counts for each cluster included in the survey. We show that RSDs in the power spectrum analysis carry important cosmological information also when traced with galaxy clusters and the DE FoM increases by a factor of 8. Besides confirming the potential that large cluster surveys have in constraining the nature of DE, our analysis emphasizes that a full exploitation of the cosmological information carried by such surveys requires not only a large statistic but also a robust measurement of the mass proxies and redshifts for a significant fraction of the cluster sample, derived from the same X-ray survey data.

At present, 30-40 per cent of the baryons in the local Universe is still undetected. According to theoretical predictions, this gas should reside in filaments filling the large-scale structure (LSS) in the form of a Warm-Hot Intergalactic Medium (WHIM), at temperatures of 10^5 - 10^7 K, thus emitting in the soft X-ray energies via free-free interaction and line emission from heavy elements. In this work we characterize the properties of the X-ray emission of the WHIM, and the LSS in general, focusing on the influence of different physical mechanisms, namely galactic winds (GWs), black-hole feedback and star-formation, and providing estimates of possible observational constraints. To this purpose we use a set of cosmological hydrodynamical simulations that include a self-consistent treatment of star-formation and chemical enrichment of the intergalactic medium, that allows us to follow the evolution of different metal species. We construct a set of simulated light-cones to make predictions of the emission in the 0.3-10 keV energy range. We obtain that GWs increase by a factor of 2 the emission of both galaxy clusters and WHIM. The amount of oxygen at average temperature and, consequently, the amount of expected bright Ovii and Oviii lines is increased by a factor of 3 due to GWs and by 20 per cent when assuming a top-heavy IMF. We compare our results with current observational constraints and find that the emission from faint groups and WHIM should account from half to all of the unresolved X-ray background in the 1-2 keV band.

In this review, we describe our current understanding of cluster formation: from the general picture of collapse from initial density fluctuations in an expanding Universe to detailed simulations of cluster formation including the effects of galaxy formation. We outline both the areas in which highly accurate predictions of theoretical models can be obtained and areas where predictions are uncertain due to uncertain physics of galaxy formation and feedback. The former includes the description of the structural properties of the dark matter halos hosting cluster, their mass function and clustering properties. Their study provides a foundation for cosmological applications of clusters and for testing the fundamental assumptions of the standard model of structure formation. The latter includes the description of the total gas and stellar fractions, the thermodynamical and non-thermal processes in the intracluster plasma. Their study serves as a testing ground for galaxy formation models and plasma physics. In this context, we identify a suitable radial range where the observed thermal properties of the intra-cluster plasma exhibit the most regular behavior and thus can be used to define robust observational proxies for the total cluster mass. We put particular emphasis on examining assumptions and limitations of the widely used self-similar model of clusters. Finally, we discuss the formation of clusters in non-standard cosmological models, such as non-Gaussian models for the initial density field and models with modified gravity, along with prospects for testing these alternative scenarios with large cluster surveys in the near future.

Forthcoming experiments will enable us to determine tomographic shear spectra at a high precision level. Most predictions about them have until now been biased on algorithms yielding the expected linear and non-linear spectrum of density fluctuations. Even when simulations have been used, so-called Halofit (Smith et al 2003) predictions on fairly large scales have been needed. We wish to go beyond this limitation. We perform N-body and hydrodynamical simulations within a sufficiently large cosmological volume to allow a direct connection between simulations and linear spectra. While covering large length-scales, the simulation resolution is good enough to allow us to explore the high-l harmonics of the cosmic shear (up to l ~ 50000), well into the domain where baryon physics becomes important. We then compare shear spectra in the absence and in presence of various kinds of baryon physics, such as radiative cooling, star formation, and supernova feedback in the form of galactic winds. We distinguish several typical properties of matter fluctuation spectra in the different simulations and test their impact on shear spectra. We compare our outputs with those obtainable using approximate expressions for non--linear spectra, and identify substantial discrepancies even between our results and those of purely N-body results. Our simulations and the treatment of their outputs however enable us, for the first time, to obtain shear results taht are fully independent of any approximate expression, also in the high-l range, where we need to incorporate a non-linear power spectrum of density perturbations, and the effects of baryon physics. This will allow us to fully exploit the cosmological information contained in future high--sensitivity cosmic shear surveys, exploring the physics of cosmic shears via weak lensing measurements.

We present a detailed investigation of the effects that a direct interaction between Dark Energy (DE) and Cold Dark Matter (CDM) particles imprints on the Halo Mass Function (HMF) of groups and clusters of galaxies. Making use of the public halo catalogs of the {\small CoDECS} simulations, we derive the HMF for several different types of coupled DE scenarios both based on the FoF algorithm and on the SO halo identification for different values of the overdensity threshold $\Delta_{c}$. We compare the computed HMFs for coupled DE cosmologies with $\Lambda $CDM as well as with the predictions of the standard analytic fitting functions. Our results show that the standard fitting functions still reproduce reasonably well both the FoF and the SO HMFs of interacting DE cosmologies at intermediate masses and at low redshifts, once rescaled to the characteristic amplitude of linear density perturbations of each specific model as given by $\sigma_{8}$. However, we also find that such apparent degeneracy with $\sigma_{8}$ is broken both by the high-mass tail and by the redshift evolution of our HMFs, with deviations beyond $\sim 10%$ for most of the models under investigation. Furthermore, the discrepancy with respect to the predictions of standard fitting functions rescaled with the characteristic value of $\sigma_{8}$ shows -- for some models -- a strong dependence on the spherical overdensity threshold $\Delta_{c}$ used for the halo identification. We find that such effect is due to a significant increase of halo concentration at low redshifts in these models, that is however absent in the majority of the cosmological scenarios considered in this work. We can therefore conclude that the universality of the HMF is violated by cosmological models that feature a direct interaction between DE and CDM.

We present an analysis of the effects of baryon physics on the halo mass function. The analysis is based on simulations of a cosmological volume. Besides a Dark Matter (DM) only simulation, we also carry out two other hydrodynamical simulations. We identified halos using a spherical overdensity algorithm and their masses are computed at three different overdensities (with respect to the critical one), $\Delta_c=200$, 500 and 1500. We find the fractional difference between halo masses in the hydrodynamical and in the DM simulations to be almost constant, at least for halos more massive than $\log (M_{\Delta_c} / \hMsun)\geq 13.5$. In this range, mass increase in the hydrodynamical simulations is of about 4-5 per cent at $\Delta_c=500$ and $\sim 1$ - 2 per cent at $\Delta_c=200$. Quite interestingly, these differences are nearly the same for both radiative and non-radiative simulations. Such variations of halo masses induce corresponding variations of the halo mass function (HMF). At $z=0$, the HMFs for GH and CSF simulations are close to the DM one, with differences of $\mincir 3$ per cent at $\Delta_c = 200$, and $\simeq 7$ per cent at $\Delta_c=500$, with $\sim 10$ - 20 per cent differences reached at $\Delta_c = 1500$. At this higher overdensity, the increase of the HMF for the radiative case is larger by about a factor 2 with respect to the non--radiative case. Assuming a constant mass shift to rescale the HMF from the hydrodynamic to the DM simulations, brings the HMF difference with respect to the DM case to be consistent with zero. Our results have interesting implications to bracket uncertainties in the mass function calibration associated to the uncertain baryon physics, in view of cosmological applications of future large surveys of galaxy clusters. (Abridged)

We present results from high--resolution cosmological hydrodynamical simulations of a Milky--Way-sized halo, aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter-stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial radius. Gas particles can be accreted through three different channels, depending on the maximum temperature value, $T_{\rm max}$, reached during the particles' past evolution: a cold channel for $T_{\rm max}<2.5 \times 10^5$ K, a hot one for $T>10^6$K, and a warm one for intermediate values of $T_{\rm max}$. We find that the warm channel is at least as important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous findings that the cold mode dominates gas accretion at high redshift. We ascribe this difference to the different supernova feedback scheme implemented in our simulations. While results presented so far in the literature are based on uneffective SN thermal feedback schemes and/or the presence of a kinetic feedback, our simulations include only effective thermal feedback. We argue that observational detections of a warm accretion mode in the high--redshift circum-galactic medium would provide useful constraints on the nature of the feedback that regulates star formation in galaxies.

We present the results of a series of numerical simulations aimed to study the evolution of a disc galaxy within the global tidal field of a group environment. Both the disc galaxy and the group are modelled as multi-component, collision-less, N-body systems, composed by both dark matter and stars. In our simulations, the evolution of disc galaxies is followed as their orbits sink towards the group centre, under the effect of dynamical friction. We explore a broad parameter space, covering several aspects of the galaxy-group interaction that are potentially relevant to galaxy evolution. Namely, prograde and retrograde orbits, orbital eccentricities, disc inclination, role of a central bulge in discs, internal disc kinematics, and galaxy-to-group mass ratios. We find that significant disc transformations occur only after the mean density of the group, measured within the orbit of the galaxy, exceeds ~0.3-1 times the central mean density of the galaxy. The morphological evolution of discs is found to be strongly dependent on the initial inclination of the disc with respect to its orbital plane. That is, discs on face-on and retrograde orbits are shown to retain longer their disc structures and kinematics, in comparison to prograde discs. This suggests that after interacting with the global tidal field alone, a significant fraction of disc galaxies should be found in the central regions of groups. Prominent central bulges are not produced, and pre-existing bulges are not enhanced in discs after the interaction with the group. Assuming that most S0 are formed in group environments, this implies that prominent bulges should be formed mostly by young stars, created only after a galaxy has been accreted by a group. Finally, contrary to some current implementations of tidal stripping in semi-analytical models of galaxy evolution, we find that more massive galaxies suffer more tidal stripping.

We study several versions of the Schmidt-Kennicutt (SK) relation obtained for isolated spiral galaxies in TreeSPH simulations run with the GADGET3 code including the novel MUlti-Phase Particle Integrator (MUPPI) algorithm for star formation and stellar feedback. [...] The standard SK relation between surface densities of cold (neutral+molecular) gas and star formation rate of simulated galaxies shows a steepening at low gas surface densities, starting from a knee whose position depends on disc gas fraction: for more gas-rich discs the steepening takes place at higher surface densities. Because gas fraction and metallicity are typically related, this environmental dependence mimics the predictions of models where the formation of H2 is modulated by metallicity. The cold gas surface density at which HI and molecular gas surface densities equate can range from ~10 up to 34 Msun/pc^2. As expected, the SK relation obtained using molecular gas shows much smaller variations among simulations. We find that disc pressure is not well represented by the classical external pressure of a disc in vertical hydrostatic equilibrium. Instead is well fit by the expression P_fit = Sigma_cold sigma_cold kappa / 6, where the three quantities on the right-hand side are cold gas surface density, vertical velocity dispersion and epicyclic frequency. When the "dynamical" SK relation, i.e. the relation that uses gas surface density divided by orbital time, is considered, we find that all of our simulations stay on the same relation. We interpret this as a manifestation of the equilibrium between energy injection and dissipation in stationary galaxy discs, when energetic feedback is effective and pressure is represented by the expression given above. These findings further support the idea that a realistic model of the structure of galaxy discs should take into account energy injection by SNe. [Abridged]

[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.

We introduce a generalized scaling law, M_tot = 10^K A^a B^b, to look for the minimum scatter in reconstructing the total mass of hydrodynamically simulated X-ray galaxy clusters, given gas mass M_gas, luminosity L and temperature T. We find a locus in the plane of the logarithmic slopes $a$ and $b$ of the scaling relations where the scatter in mass is minimized. This locus corresponds to b_M = -3/2 a_M +3/2 and b_L = -2 a_L +3/2 for A=M_gas and L, respectively, and B=T. Along these axes, all the known scaling relations can be identified (at different levels of scatter), plus a new one defined as M_tot ~ (LT)^(1/2). Simple formula to evaluate the expected evolution with redshift in the self-similar scenario are provided. In this scenario, no evolution of the scaling relations is predicted for the cases (b_M=0, a_M=1) and (b_L=7/2, a_L=-1), respectively. Once the single quantities are normalized to the average values of the sample under considerations, the normalizations K corresponding to the region with minimum scatter are very close to zero. The combination of these relations allows to reduce the number of free parameters of the fitting function that relates X-ray observables to the total mass and includes the self-similar redshift evolution.

We study the amount and distribution of dark matter substructures within dark matter haloes, using a large set of high-resolution simulations ranging from group size to cluster size haloes, and carried our within a cosmological model consistent with WMAP 7-year data. In particular, we study how the measured properties of subhaloes vary as a function of the parent halo mass, the physical properties of the parent halo, and redshift. The fraction of halo mass in substructures increases with increasing mass. There is, however, a very large halo-to-halo scatter that can be explained only in part by a range of halo physical properties, e.g. concentration. At given halo mass, less concentrated haloes contain significantly larger fractions of mass in substructures because of the reduced strength of tidal disruption. Most of the substructure mass is located at the outskirts of the parent haloes, in relatively few massive subhaloes. This mass segregation appears to become stronger at increasing redshift, and should reflect into a more significant mass segregation of the galaxy population at different cosmic epochs. When haloes are accreted onto larger structures, their mass is significantly reduced by tidal stripping. Haloes that are more massive at the time of accretion (these should host more luminous galaxies) are brought closer to the centre on shorter time-scales by dynamical friction, and therefore suffer of a more significant stripping. The halo merger rate depends strongly on the environment with substructure in more massive haloes suffering more important mergers than their counterparts residing in less massive systems. This should translate into a different morphological mix for haloes of different mass.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We examine sources of scatter in scaling relations between galaxy cluster mass and thermal Sunyaev-Zeldovich (SZ) effect using cluster samples extracted from cosmological hydrodynamical simulations. Overall, the scatter of the mass-SZ scaling relation is well correlated with the scatter in the mass-halo concentration relation with more concentrated halos having stronger integrated SZ signals at fixed mass. Additional sources of intrinsic scatter are projection effects from correlated structures, which cause the distribution of scatter to deviate from log-normality and skew it towards higher inferred masses, and the dynamical state of clusters. We study the evolution of merging clusters based on simulations of 39 clusters and their cosmological environment with high time resolution. This sample enables us to study for the first time the detailed evolution of merging clusters around the scaling relation for a cosmologically representative distribution of merger parameters. Major mergers cause an asymmetric scatter such that the inferred mass of merging systems is biased low. We find mergers to be the dominant source of bias towards low inferred masses: over 50% of outliers on this side of the scaling relation underwent a major merger within the last Gigayear. As the fraction of dynamically disturbed clusters increases with redshift, our analysis indicates that mergers cause a redshift-dependent bias in scaling relations. Furthermore, we find the SZ morphology of massive clusters to be well correlated with the clusters' dynamical state, suggesting that morphology may be used to constrain merger fractions and identify merger-induced outliers of the scaling relation.

We present first results on PLCK G266.6-27.3, a galaxy cluster candidate detected at a signal-to-noise ratio of 5 in the Planck All Sky survey. An XMM-Newton validation observation has allowed us to confirm that the candidate is a bona fide galaxy cluster. With these X-ray data we measure an accurate redshift, z = 0.94 +/- 0.02, and estimate the cluster mass to be M_500 = (7.8 +/- 0.8)e+14 solar masses. PLCK G266.6-27.3 is an exceptional system: its luminosity of L_X(0.5-2.0 keV)=(1.4 +/- 0.05)e+45 erg/s, equals that of the two most luminous known clusters in the z > 0.5 universe, and it is one of the most massive clusters at z~1. Moreover, unlike the majority of high-redshift clusters, PLCK G266.6-27.3 appears to be highly relaxed. This observation confirms Planck's capability of detecting high-redshift, high-mass clusters, and opens the way to the systematic study of population evolution in the exponential tail of the mass function.

Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formation takes place. By necessity, these merger trees are constructed from a finite set of discrete "snapshots" of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of "snapshot" times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.

Magnetic fields in the intra-cluster medium of galaxy clusters have been studied in the past years through different methods. In the next years the up-coming generation of radio telescopes is going to provide new data that have the potential of setting constraints on the properties of magnetic fields in galaxy clusters. Here we present zoomed-in simulations for a set of massive galaxy clusters (M_v > 10^15 M_sun/h). This is an ideal sample to study the evolution of magnetic field during the process of structure formation in detail. Turbulent motions of the gas within the ICM will manifest themselves in a macroscopic magnetic resistivity eta_m, which has to be taken explicitly into account, especially at scales below the resolution limit. We have adapted the MHD GADGET code by Dolag & Stasyszyn (2009) to include the treatment of the magnetic resistivity and for the first time we have included non-ideal MHD equations to better follow the evolution of the magnetic field within galaxy clusters. We investigate which value of the magnetic resistivity eta_m is required to match the magnetic field profile derived from radio observations. We find that a value of eta_m ~ 6 10^27 cm^2/s is necessary to recover the shape of the magnetic field profile inferred from radio observations of the Coma cluster. This value agrees well with the expected level of turbulent motions within the ICM at our resolution limit. The magnetic field profiles can be fitted by a beta-model like profile (Cavaliere & Fusco-Femiano 1976), with small dispersion of the parameters. We find also that that the temperature, density and entropy profiles of the clusters depend on the magnetic resistivity constant, having flatter profiles in the inner regions when the magnetic resistivity increases.

It has been a long-standing question whether fossil groups are just sampling the tail of the distribution of ordinary groups, or whether they are a physically distinct class of objects, characterized by an unusual and special formation history. To study this question, we here investigate fossil groups identified in the hydrodynamical simulations of the GIMIC project, which consists of resimulations of five regions in the Millennium Simulation (MS) that are characterized by different large-scale densities, ranging from a deep void to a proto-cluster region. For comparison, we also consider semi-analytic models built on top of the MS, as well as a conditional luminosity function approach. We identify galaxies in the GIMIC simulations as groups of stars and use a spectral synthesis code to derive their optical properties. The X-ray luminosity of the groups is estimated in terms of the thermal bremsstrahlung emission of the gas in the host halos, neglecting metallicity effects. We focus on comparing the properties of fossil groups in the theoretical models and observational results, highlighting the differences between them, and trying to identify possible dependencies on environment for which our approach is particularly well set-up. We find that the optical fossil fraction in all of our theoretical models declines with increasing halo mass, and there is no clear environmental dependence. Combining the optical and X-ray selection criteria for fossil groups, the halo mass dependence of the fossil groups seen in optical vanishes. Over the GIMIC halo mass range we resolve best, 9.0\times1012 \sim 4.0\times1013 h-1 M, the central galaxies in the fossil groups show similar properties as those in ordinary groups, in terms of age, metallicity, color, concentration, and mass-to-light ratio. [abridged]

We present results based on an implementation of the Godunov Smoothed Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of the equations of moment and energy conservation, starting from the convolution of these equations with the interpolating kernel. The two most important aspects of the numerical implementation of these equations are (a) the appearance of fluid velocity and pressure obtained from the solution of the Riemann problem between each pair of particles, and (b the absence of an artificial viscosity term. We carry out three different controlled hydrodynamical three-dimensional tests, namely the Sod shock tube, the development of Kelvin-Helmholtz instabilities in a shear flow test, and the "blob" test describing the evolution of a cold cloud moving against a hot wind. The results of our tests confirm and extend in a number of aspects those recently obtained by Cha (2010): (i) GSPH provides a much improved description of contact discontinuities, with respect to SPH, thus avoiding the appearance of spurious pressure forces; (ii) GSPH is able to follow the development of gas-dynamical instabilities, such as the Kevin--Helmholtz and the Rayleigh-Taylor ones; (iii) as a result, GSPH describes the development of curl structures in the shear-flow test and the dissolution of the cold cloud in the "blob" test. We also discuss in detail the effect on the performances of GSPH of changing different aspects of its implementation. The results of our tests demonstrate that GSPH is in fact a highly promising hydrodynamic scheme, also to be coupled to an N-body solver, for astrophysical and cosmological applications. [abridged]

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

The ubiquitous presence of the Fe line complex in the X-ray spectra of galaxy clusters offers the possibility of measuring their redshift without resorting to spectroscopic follow-up observations. In this paper we assess the accuracy with which the redshift of galaxy clusters can be recovered from an X-ray spectral analysis of Chandra archival data. This study indicates a strategy to build large surveys of clusters whose identification and redshift measurement are both based on X-ray data alone. We apply a blind search for K--shell and L--shell Fe line complex in X-ray cluster spectra using Chandra archival observations of galaxy clusters. The Fe line in the ICM spectra can be detected by simply analyzing the C-statistics variation $\Delta C_{stat}$ as a function of the redshift parameter. We repeat the measurement under different conditions, and compare the X-ray derived redshift $z_X$ with the one obtained by means of optical spectroscopy $z_o$. We explore how a number of priors on metallicity and luminosity can be effectively used to reduce catastrophic errors. The $\Delta C_{stat}$ provides the most efficient means for discarding wrong redshift measures and to estimate the actual error on $z_X$. We identify a simple and efficient procedure for optimally measuring the redshifts from the X-ray spectral analysis of clusters of galaxies. When this procedure is applied to mock catalogs extracted from high sensitivity, wide-area cluster surveys, such as those proposed with Wide Field X-ray Telescope (WFXT) mission, it is possible to obtain a complete samples of X-ray clusters with reliable redshift measurements, thus avoiding time-consuming optical spectroscopic observations. This methodology will make it possible to trace cosmic growth by studying the evolution of the cluster mass function directly using X-ray data.

We present a detailed study of scaling relations between total cluster mass and three mass proxies based on X-ray observables: temperature of the intra-cluster medium, gas mass and the product of the two, Y_X. Our analysis is based on two sets of high-resolution hydrodynamical simulations performed with the TreePM-SPH GADGET code. The first set includes about 140 clusters with masses above 5x10^13 M_sun/h (30 having mass above 10^15 M_sun/h), that have been simulated with (i) non-radiative physics and including (ii) cooling, star formation, chemical enrichment and the effect of supernova feedback triggering galactic ejecta. This large statistics is used to quantify the robustness of the scaling relations, to determine their redshift evolution and to calibrate their intrinsic scatter and its distribution. We use a smaller set of clusters including 18 halos with masses above 5x10^13 M_sun/h to test the robustness of mass proxies against changing the physical processes included in simulations (thermal conduction, artificial viscosity, cooling and star formation, galactic winds and AGN feedback). We find the M-Y_X scaling relation to be the least sensitive one to variations of the ICM physics, with its slope and redshift evolution close to the self-similar model predictions. The distribution of the scatter around the best-fitting relations is close to a log-normal one. M_gas has the smallest scatter in mass, with values of sigma_lnM = 0.04-0.06, depending on the physics included in the simulation, and with a mild dependence on redshift. The M-T relation is the one with the largest scatter, with sigma_lnM > 0.1 at z=0, increasing to > 0.15 at z=1. The intrinsic scatter in the M-Y_X relation is slightly larger than in the M-M_gas relation. These results confirm that both Y_X and M_gas mass proxies are well suited for cosmological applications of future large X-ray surveys. [abridged]

The most accepted scenario for the origin of fossil groups (FGs) is that they are galaxy associations in which the merging rate was fast and efficient. These systems have assembled half of their mass at early epoch of the Universe, subsequently growing by minor mergers. They could contain a fossil record of the galaxy structure formation. We have started a project in order to characterize a large sample of FGs. In this paper we present the analysis of the fossil system RX J105453.3+552102. Optical deep images were used for studying the properties of the brightest group galaxy and for computing the photometric luminosity function of the group. We have also performed a detail dynamical analysis of the system based on redshift data for 116 galaxies. This galaxy system is located at z=0.47, and shows a quite large line-of-sight velocity dispersion \sigma_{v}~1000 km/s. Assuming the dynamical equilibrium, we estimated a virial mass of M ~ 10^{15} h_{70} M_{\odot}. No evidence of substructure was found within 1.4 Mpc radius. We found a statistically significant departure from Gaussianity of the group members velocities in the most external regions of the group. This could indicate the presence of galaxies in radial orbits in the external region of the group. We also found that the photometrical luminosity function is bimodal, showing a lack of M_{r} ~ -19.5 galaxies. The brightest group galaxy shows low Sersic parameter (n~2) and a small peculiar velocity. Indeed, our accurate photometry shows that the difference between the brightest and the second brightest galaxies is 1.9 mag in the r-band, while the classical definition of FGs is based on a magnitude gap of 2. We conclude that this fossil system does not follow the empirical definition of FGs. Nevertheless, it is a massive, old and undisturbed galaxy system with little infall of L^{*} galaxies since its initial collapse.

We list here the contents of the Proceedings of the "Wide Field X-ray Telescope" conference held in Bologna, Italy on 25-26 Nov 2009. The conference highlighted the scientific potential and discovery space provided by an X-ray mission concept characterized by a wide field-of-view (1 sq.deg.), large effective area (1 sq.mt.) and approximately constant PSF (~5 arcsec HEW) across the whole FOV. The index is in html form with clickable links to the individual contributions.

We investigate the properties of triply ionized Carbon (CIV) in the Intergalactic Medium using a set of high-resolution and large box-size cosmological hydrodynamic simulations of a $\Lambda$CDM model. We rely on a modification of the GADGET-2 code that self-consistently follows the metal enrichment mechanism by means of a detailed chemical evolution model. We focus on several numerical implementations of galactic feedback: galactic winds in the energy driven and momentum driven prescriptions and Active Galactic Nuclei (AGN) powered by gas accretion onto massive black holes. We extract mock IGM transmission spectra in neutral hydrogen (HI) and CIV and perform Voigt profile fitting. The results are then compared with high-resolution quasar (QSO) spectra obtained with the UVES spectrograph at the VLT and the HIRES spectrograph at Keck. We find that feedback has little impact on statistics related to the neutral hydrogen, while CIV is more affected by galactic winds and/or AGN feedback. When the same analysis is performed over observed and simulated CIV lines, we find reasonables good agreement between data and simulations over the column density range $N_{\rm CIV}=10^{12.5-15}$ cm$^{-2}$. Also the CIV line-widths distribution appears to be in agreement with the observed values, while the HI Doppler parameters, $b_{\rm HI}$, are in general too large showing that the diffuse cosmic web is heated more than what is inferred by observations. The simulation without feedback fails in reproducing the CIV systems at high column densities at all redshift, while the AGN feedback case agrees with observations only at $z<3$, when this form of feedback is particularly effective. We also present scatter plots in the $b-N$ and in the $N_{\rm CIV}-N_{\rm HI}$ planes, showing that there is rough agreement between observations and simulations only when feedback is taken into account.

The caustic technique uses galaxy redshifts alone to measure the escape velocity and mass profiles of galaxy clusters to clustrocentric distances well beyond the virial radius, where dynamical equilibrium does not necessarily hold. We provide a detailed description of this technique and analyse its possible systematic errors. We apply the caustic technique to clusters with mass M_200>=10^{14}h^{-1} M_sun extracted from a cosmological hydrodynamic simulation of a LambdaCDM universe. With a few tens of redshifts per squared comoving megaparsec within the cluster, the caustic technique, on average, recovers the profile of the escape velocity from the cluster with better than 10 percent accuracy up to r~4 r_200. The caustic technique also recovers the mass profile with better than 10 percent accuracy in the range (0.6-4) r_200, but it overestimates the mass up to 70 percent at smaller radii. This overestimate is a consequence of neglecting the radial dependence of the filling function F_beta(r). The 1-sigma uncertainty on individual escape velocity profiles increases from ~20 to ~50 percent when the radius increases from r~0.1 r_200 to ~4 r_200. Individual mass profiles have 1-sigma uncertainty between 40 and 80 percent within the radial range (0.6-4) r_200. We show that the amplitude of these uncertainties is completely due to the assumption of spherical symmetry, which is difficult to drop. Alternatively, we can apply the technique to synthetic clusters obtained by stacking individual clusters: in this case, the 1-sigma uncertainty on the escape velocity profile is smaller than 20 percent out to 4 r_200. The caustic technique thus provides reliable average profiles which extend to regions difficult or impossible to probe with other techniques.

We discuss the central role played by the X-ray study of hot baryons within galaxy clusters to reconstruct the assembly of cosmic structures and to trace the past history of star formation and accretion onto supermassive Black Holes (BHs). We shortly review the progress in this field contributed by the current generation of X-ray telescopes. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM: (a) When and how is entropy injected into the inter-galactic medium (IGM)? (b) What is the history of metal enrichment of the IGM? (c) What physical mechanisms determine the presence of cool cores in galaxy clusters? (d) How is the appearance of proto-clusters at z~2 related to the peak of star formation activity and BH accretion? (e) What do galaxy clusters tell us about the nature of primordial density perturbations and on the history of their growth? We show that the most efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. A similar survey can only be carried out with a Wide-Field X-ray Telescope (WFXT), which combines a high survey speed with a sharp PSF across the entire FoV. We emphasize the important synergies that WFXT will have with a number of future ground-based and space telescopes, covering from the radio to the X-ray bands. Finally, we discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

The Wide Field X-Ray Telescope (WFXT) is a medium-class mission designed to be 2-orders-of-magnitude more sensitive than any previous or planned X-ray mission for large area surveys and to match in sensitivity the next generation of wide-area optical, IR and radio surveys. Using an innovative wide-field X-ray optics design, WFXT provides a field of view of 1 square degree (10 times Chandra) with an angular resolution of 5" (Half Energy Width, HEW) nearly constant over the entire field of view, and a large collecting area (up to 1 m^2 at 1 keV, > 10x Chandra) over the 0.1-7 keV band. WFXTs low-Earth orbit also minimizes the particle background. In five years of operation, WFXT will carry out three extragalactic surveys at unprecedented depth and address outstanding questions in astrophysics, cosmology and fundamental physics. In this article, we illustrate the mission concept and the connection between science requirements and mission parameters.

We discuss some of the main open issues in the evolution of Active Galactic Nuclei which can be solved by the sensitive, wide area surveys to be performed by the proposed Wide Field X-ray Telescope mission.

We investigate the scientific impact of the Wide Field X-ray Telescope mission. We present simulated images and spectra of X-ray sources as observed from the three surveys planned for the nominal 5-year WFXT lifetime. The goal of these simulations is to provide WFXT images of the extragalactic sky in different energy bands based on accurate description of AGN populations, normal and star forming galaxies, groups and clusters of galaxies. The images are realized using a detailed PSF model, instrumental and physical backgrounds/foregrounds, accurate model of the effective area and the related vignetting effect. Thanks to this comprehensive modelization of the WFXT properties, the simulated images can be used to evaluate the flux limits for detection of point and extended sources, the effect of source confusion at very faint fluxes, and in general the efficiency of detection algorithms. We also simulate the spectra of the detected sources, in order to address specific science topics which are unique to WFXT. Among them, we focus on the characterization of the Intra Cluster Medium (ICM) of high-z clusters, and in particular on the measurement of the redshift from the ICM spectrum in order to build a cosmological sample of galaxy clusters. The end-to-end simulation procedure presented here, is a valuable tool in optimizing the mission design. Therefore, these simulations can be used to reliably characterize the WFXT discovery space and to verify the connection between mission requirements and scientific goals. Thanks to this effort, we can conclude on firm basis that an X-ray mission optimized for surveys like WFXT is necessary to bring X-ray astronomy at the level of the optical, IR, submm and radio wavebands as foreseen in the coming decade.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

We present forecasts for constraints on deviations from Gaussian distribution of primordial density perturbations from future X-ray surveys of galaxy clusters. Our analysis is based on computing the Fisher-Matrix for number counts and large-scale power spectrum of clusters. We consider a survey with high-sensitivity and wide-area to detect about 2.5 x 10^5 extended sources. Based on the self-calibration approach, and including Planck priors in our analysis, we constrain 9 cosmological parameters and 4 nuisance parameters, which define the relation between cluster mass and X-ray flux. Because of the scale dependence of large-scale bias induced by local-shape non-Gaussianity, we find that the power spectrum provides strong constraints on the non-Gaussianity f_nl parameter, which complement the stringent constraints on the power spectrum normalization, \sigma_8, from the number counts. To quantify the joint constraints on theese two parameters, that specify the timing of structure formation for a fixed background expansion, we define the figure-of-merit FoM_SFT = (det[Cov(sigma_8,f_nl)])^{-1/2}. We find that our surveys constrain deviations from Gaussianity with a precision of \Delta f_nl~10 at 1 \sigma confidence level, with FoM_SFT~39. We point out that constraints on f_nl are weakly sensitive to the uncertainties in the knowledge of nuisance parameters. As an application of non-Gaussian constraints from available data, we analyse the impact of positive skewness on the occurrence of XMMU-J2235, a massive distant cluster discovered at z~1.4. We confirm that in a WMAP-7 Gaussian Lambda CDM cosmology, within the survey volume, ~5x10^{-3} objects like this are expected to be found: to increase this probability by a factor of at least 10, one needs to evade either the available constraints on f_nl or sigma_8.

We present a new multi-phase sub-resolution model for star formation and feedback in SPH numerical simulations of galaxy formation. Our model, called MUPPI (MUlti-Phase Particle Integrator), describes each gas particle as a multi-phase system, with cold and hot gas phases, coexisting in pressure equilibrium, and a stellar component. Cooling of the hot tenuous gas phase feeds the cold gas phase. Stars are formed out of molecular gas with a given efficiency, which scales with the dynamical time of the cold phase. Our prescription for star formation is not based on imposing the Schmidt-Kennicutt relation, which is instead naturally produced by MUPPI. Energy from supernova explosions is deposited partly into the hot phase of the gas particles, and partly to that of neighboring particles. Mass and energy flows among the different phases of each particle are described by a set of ordinary differential equations which we explicitly integrate for each gas particle, instead of relying on equilibrium solutions. This system of equations also includes the response of the multi-phase structure to energy changes associated to the thermodynamics of the gas. We apply our model to two isolated disk galaxy simulations and two spherical cooling flows. MUPPI is able to reproduce the Schmidt-Kennicutt relation for disc galaxies. It also reproduces the basic properties of the inter-stellar medium in disc galaxies, the surface densities of cold and molecular gas, of stars and of star formation rate, the vertical velocity dispersion of cold clouds and the flows connected to the galactic fountains. Quite remarkably, MUPPI also provides efficient stellar feedback without the need to include a scheme of kinetic energy feedback. [abridged]

Member galaxies within galaxy clusters nowadays can be routinely identified in cosmological, hydrodynamical simulations using methods based on identifying self bound, locally over dense substructures. However, distinguishing the central galaxy from the stellar diffuse component within clusters is notoriously difficult, and in the center it is not even clear if two distinct stellar populations exist. Here, after subtracting all member galaxies, we use the velocity distribution of the remaining stars and detect two dynamically, well-distinct stellar components within simulated galaxy clusters. These differences in the dynamics can be used to apply an un-binding procedure which leads to a spatial separation of the two components into a cD and a diffuse stellar component (DSC). Applying our new algorithm to a cosmological, hydrodynamical simulation we find that -- in line with previous studies -- these two components have clearly distinguished spatial and velocity distributions as well as different star formation histories. We show that the DSC fraction -- which can broadly be associated with the observed intra cluster light -- does not depend on the virial mass of the galaxy cluster and is much more sensitive to the formation history of the cluster. We conclude that the separation of the cD and the DSC in simulations, based on our dynamical criteria, is more physically motivated than current methods which depend on implicit assumptions on a length scale associated with the cD galaxy and therefore represent a step forward in understanding the different stellar components within galaxy clusters. Our results also show the importance of analyzing the dynamics of the DSC to characterize its properties and understand its origin.

We present a detailed comparison between the galaxy populations within a massive cluster, as predicted by hydrodynamical SPH simulations and by a semi-analytic model (SAM) of galaxy formation. Both models include gas cooling and a simple prescription of star formation, which consists in transforming instantaneously any cold gas available into stars, while neglecting any source of energy feedback. We find that, in general, galaxy populations from SAMs and SPH have similar statistical properties, in agreement with previous studies. However, when comparing galaxies on an object-by-object basis, we find a number of interesting differences: a) the star formation histories of the brightest cluster galaxies (BCGs) from SAM and SPH models differ significantly, with the SPH BCG exhibiting a lower level of star formation activity at low redshift, and a more intense and shorter initial burst of star formation with respect to its SAM counterpart; b) while all stars associated with the BCG were formed in its progenitors in the semi-analytic model used here, this holds true only for half of the final BCG stellar mass in the SPH simulation, the remaining half being contributed by tidal stripping of stars from the diffuse stellar component associated with galaxies accreted on the cluster halo; c) SPH satellites can loose up to 90 per cent of their stellar mass at the time of accretion, due to tidal stripping, a process not included in the semi-analytic model used in this study; d) in the SPH simulation, significant cooling occurs on the most massive satellite galaxies and this lasts for up to 1 Gyr after accretion. This physical process is not included in the semi-analytic model used in our study, as well as in most of the models discussed in the recent literature.

We study the substructure statistics of a representative sample of galaxy clusters by means of two currently popular substructure characterisation methods, power ratios and centroid shifts. We use the 31 clusters from the REXCESS sample, compiled from the southern ROSAT All-Sky cluster survey REFLEX with a morphologically unbiased selection in X-ray luminosity and redshift, all of which have been reobserved with XMM-Newton. We investigate the uncertainties of the substructure parameters and examine the dependence of the results on projection effects, finding that the uncertainties of the parameters can be quite substantial. Thus while the quantification of the dynamical state of individual clusters with these parameters should be treated with extreme caution, these substructure measures provide powerful statistical tools to characterise trends of properties in large cluster samples. The centre shift parameter, w, is found to be more sensitive in general. For the REXCESS sample neither the occurence of substructure nor the presence of cool cores depends on cluster mass. There is a significant anti-correlation between the existence of substantial substructure and cool cores. The simulated clusters show on average larger substructure parameters than the observed clusters, a trend that is traced to the fact that cool regions are more pronounced in the simulated clusters, leading to stronger substructure measures in merging clusters and clusters with offset cores. Moreover, the frequency of cool regions is higher in the simulations than in the observations, implying that the description of the physical processes shaping cluster formation in the simulations requires further improvement.

(Abridged) We examine the radial entropy distribution and its scaling using 31 nearby galaxy clusters from the Representative XMM-Newton Cluster Structure Survey (REXCESS). The entropy profiles are robustly measured at least out to R_1000 in all systems and out to R_500 in 13 systems. Compared to theoretical expectations, the observed distributions show a radial and mass-dependent excess entropy that is greater and extends to larger radii in lower mass systems. At R_500, the mass dependence and entropy excess are both negligible within the uncertainties. Mirroring this behaviour, the scaling of gas entropy is shallower than self-similar in the inner regions, but steepens with radius, becoming consistent with self-similar at R_500. The dispersion in scaled entropy in the inner regions is linked to the presence of cool cores and dynamical activity; at larger radii the dispersion decreases by a factor of two and the dichotomy between subsamples disappears. Parameterising the profiles with a power law plus constant model, there are two peaks in central entropy K_0; however, we cannot distinguish between a bimodal or a left-skewed distribution. The outer slopes are correlated with system temperature; their distribution is unimodal with a median value of 0.98. Renormalising the dimensionless entropy profiles by the gas mass fraction profile f_gas(< R), leads to a remarkable reduction in the scatter, implying that gas mass fraction variations with radius and mass are the cause of the observed entropy properties. We discuss a tentative scenario to explain the behaviour of the entropy and gas mass fraction in the REXCESS sample, in which extra heating and merger mixing maintains an elevated central entropy level in the majority of the population, and a smaller fraction of systems develops a cool core.

We analyse the evolution of the properties of the low-redshift Intergalactic Medium (IGM) using high-resolution hydrodynamic simulations that include a detailed chemical evolution model. We focus on the effects that two different forms of energy feedback, strong galactic winds driven by supernova explosion and Active Galactic Nuclei (AGN) powered by gas accretion onto super-massive black holes (BHs), have on the thermo- and chemo-dynamical properties of of the low redshift IGM. We find that feedback associated to winds (W) and BHs leave distinct signatures in both the chemical and thermal history of the baryons, especially at redshift z<3 [..] We present results for the enrichment in terms of mass and metallicity distributions for the WHIM phase, both as a function of density and temperature. Finally, we compute the evolution of the relative abundances between different heavy elements, namely Oxygen, Carbon and Iron. While both C/O and O/Fe evolve differently at high redshifts for different feedback models, their values are similar at z=0 [..]. The sensitivity of WHIM properties on the implemented feedback scheme could be important both for discriminating between different feedback physics and for detecting the WHIM with future far-UV and X-ray telescopes.

Combining cosmological SPH simulations with a previously developed Lyman Alpha production/transmission model and the Early Reionization Model (ERM, reionization ends at redshift z~7), we obtain Lyman Alpha and UV Luminosity Functions (LFs) for Lyman Alpha Emitters (LAEs) for redshifts between 5.7 and 7.6. Matching model results to observations at z~5.7 requires escape fractions of Lyman Alpha, f_alpha=0.3, and UV (non-ionizing) continuum photons, f_c=0.22, corresponding to a color excess, E(B-V)=0.15. We find that (i) f_c increases towards higher redshifts, due the decreasing mean dust content of galaxies, (ii) the evolution of f_alpha/f_c hints at the dust content of the ISM becoming progressively inhomogeneous/clumped with decreasing redshift. The clustering photoionization boost is important during the initial reionization phases but has little effect on the Lyman Alpha LF for a highly ionized IGM. Halo (stellar) masses are in the range 10.0 < \log M_h < 11.8 (8.1 < \log M_* < 10.4) with M_h \propto M_*^{0.64}. The star formation rates are between 3-120 solar masses per year, mass-weighted mean ages are greater than 20 Myr at all redshifts, while the mean stellar metallicity increases from Z=0.12 to 0.22 solar metallicity from z~7.6 to z~5.7; both age and metallicity positively correlate with stellar mass. The brightest LAEs are all characterized by large star formation rates and intermediate ages (~200 Myr), while objects in the faint end of the Lyman Alpha LF show large age and star formation rate spreads. With no more free parameters, the Spectral Energy Distributions of three LAE at z~5.7 observed by Lai et al. (2007) are well reproduced by an intermediate age (182-220 Myr) stellar population and the above E(B-V) value.

We present a study of the effect of AGN feedback on metal enrichment and thermal properties of the intracluster medium (ICM) in hydrodynamical simulations. The cosmological simulations are performed for a set of clusters using a version of the TreePM-SPH Gadget code that follows chemo-dynamical evolution by accounting for metal enrichment by different stellar populations. Besides runs not including any efficient form of energy feedback, we carry out simulations including: (i) kinetic feedback in the form of galactic winds triggered by supernova explosions; (ii) AGN feedback from gas accretion onto super-massive black holes (BHs); (iii) AGN feedback in which a 'radio mode' is included. We find that AGN feedback is able to quench star formation in the brightest cluster galaxies at z<4 and provides correct temperature profiles in the central regions of galaxy groups. However, its effect is not sufficient to create cool cores in massive clusters. AGN feedback creates a widespread enrichment in the outskirts of clusters, thanks to its efficiency in displacing enriched gas from galactic halos to the inter-galactic medium at relatively high redshift. Iron abundance profiles are in better agreement with observations, with a more pristine enrichment of the ICM around and beyond the cluster virial regions. From the pattern of the relative abundances of Silicon and Iron, we conclude that a significant fraction of ICM enrichment in simulations is contributed by a diffuse population of intra-cluster stars. Our simulations also predict that profiles of Z_Si/Z_Fe abundance ratio do not increase at least out to 0.5 R_vir. Our results clearly show that different sources of energy feedback leave distinct imprints in the enrichment pattern of the ICM, that are more evident when looking at cluster external regions.

We review recent progress in the description of the formation and evolution of galaxy clusters in a cosmological context by using numerical simulations. We focus our presentation on the comparison between simulated and observed X-ray properties, while we will also discuss numerical predictions on properties of the galaxy population in clusters. Many of the salient observed properties of clusters, such as X-ray scaling relations, radial profiles of entropy and density of the intracluster gas, and radial distribution of galaxies are reproduced quite well. In particular, the outer regions of cluster at radii beyond about 10 per cent of the virial radius are quite regular and exhibit scaling with mass remarkably close to that expected in the simplest case in which only the action of gravity determines the evolution of the intra-cluster gas. However, simulations generally fail at reproducing the observed cool-core structure of clusters: simulated clusters generally exhibit a significant excess of gas cooling in their central regions, which causes an overestimate of the star formation and incorrect temperature and entropy profiles. The total baryon fraction in clusters is below the mean universal value, by an amount which depends on the cluster-centric distance and the physics included in the simulations, with interesting tensions between observed stellar and gas fractions in clusters and predictions of simulations. Besides their important implications for the cosmological application of clusters, these puzzles also point towards the important role played by additional physical processes, beyond those already included in the simulations. We review the role played by these processes, along with the difficulty for their implementation, and discuss the outlook for the future progress in numerical modeling of clusters.

The abundance and structure of dark matter subhalos has been analyzed extensively in recent studies of dark matter-only simulations, but comparatively little is known about the impact of baryonic physics on halo substructures. We here extend the SUBFIND algorithm for substructure identification such that it can be reliably applied to dissipative hydrodynamical simulations that include star formation. This allows, in particular, the identification of galaxies as substructures in simulations of clusters of galaxies, and a determination of their content of gravitationally bound stars, dark matter, and hot and cold gas. Using a large set of cosmological cluster simulations, we present a detailed analysis of halo substructures in hydrodynamical simulations of galaxy clusters, focusing in particular on the influence both of radiative and non-radiative gas physics, and of non-standard physics such as thermal conduction and feedback by galactic outflows. We also examine the impact of numerical nuisance parameters such as artificial viscosity parameterizations. We find that diffuse hot gas is efficiently stripped from subhalos when they enter the highly pressurized cluster atmosphere. This has the effect of decreasing the subhalo mass function relative to a corresponding dark matter-only simulation. These effects are mitigated in radiative runs, where baryons condense in the central subhalo regions and form compact stellar cores. However, in all cases, only a very small fraction, of the order of one percent, of subhalos within the cluster virial radii preserve a gravitationally bound hot gaseous atmosphere. (abridged)

(Abriged) We present the analysis of the baryonic content of 52 X-ray luminous galaxy clusters observed with Chandra in the redshift range 0.3-1.273. We use the deprojected X-ray surface brightness profiles and the measured values of the gas temperature to recover the gas and total mass profiles. By assuming that galaxy clusters are representative of the cosmic baryon budget, the distribution of the cluster baryon fraction in the hottest (T> 4 keV) systems as a function of redshift is used to constrain the cosmological parameters. We discuss how our constraints are affected by several systematics, namely the isothermality, the assumed baryon fraction in stars, the depletion parameter and the sample selection. By using only the cluster baryon fraction as a proxy for the cosmological parameters, we obtain that Omega is very well constrained at the value of 0.35 with a relative statistical uncertainty of 11% (1 sigma level; w=-1) and a further systematic error of about (-6,+7)%. On the other hand, constraints on Lambda (without the prior of flat geometry) and w (using the prior of flat geometry) are definitely weaker due to the presence of larger statistical and systematic uncertainties (of the order of 40 per cent on Lambda and larger than 50 per cent on w). If the WMAP 5-year best-fit results are assumed to fix the cosmological parameters, we limit the contributions expected from non-thermal pressure support and ICM clumpiness to be lower than about 10 per cent, leaving also room to accommodate baryons not accounted for either in the X-ray emitting plasma or in stars of the order of 18 per cent of the total cluster baryon budget.

We investigate the properties of Damped Lyman-alpha systems (DLAs) using high-resolution and large box-size cosmological hydrodynamical simulations of a LambdaCDM model. The numerical code used is a modification of GADGET-2 with a self consistent implementation of the metal enrichment mechanism (Tornatore et al. 2007). We explore the numerical convergence of some relevant physical quantities and we vary the parameters describing the properties of galactic winds; the initial stellar mass function; the linear dark matter power spectrum and the metal enrichment pattern of the IGM (Intergalactic Medium) around DLAs. We find that strong galactic winds with speed of about 600 km/s, in an energy-driven wind scenario, are needed in order to match the observed column density distribution function for DLAs and the evolution of the neutral hydrogen content with redshift. The momentum-driven implementation of the galactic wind model, that relates the speed and mass load in the wind to the properties of the dark matter haloes, shows a behaviour which is intermediate between the energy-driven galactic winds of small (100 km/s) and large (600 km/s) velocities. At z=3 the contribution of haloes of masses between 10^9 and 10^10 h^-1 M_sun, for DLAs below 10^20.8 cm^-2, to the column density distribution function, is significant. By interpolating physical quantities along line-of-sights through massive haloes we qualitatively show how different galactic wind models impact on the IGM around DLAs. Furthermore, we analyse statistics related to the velocity widths of SiII associated to DLAs: the metallicity in the wind seems to be rather clumpy and this produces an underestimation of the observed velocity widths. We outline possible solutions to this problem.

We discuss the central role played by X-ray studies to reconstruct the past history of formation and evolution of supermassive Black Holes (BHs), and the role they played in shaping the properties of their host galaxies. We shortly review the progress in this field contributed by the current X-ray and multiwavelength surveys. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM, as for X-ray observations, and the legacy of the SDSS, as for wide area surveys: 1) When and how did the first supermassive black holes form? 2) How does cosmic environment regulate nuclear activity (and star formation) across cosmic time? 3) What is the history of nuclear activity in a galaxy lifetime? We show that the most efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. Such a survey can only be carried out with a Wide-Field X-ray Telescope (WFXT) with a high survey speed, due to the combination of large field of view and large effective area, i.e., grasp, and sharp PSF. We emphasize the important synergies that WFXT will have with a number of future groundbased and space telescopes, covering from the radio to the X-ray bands and discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

Sensitive, wide-area X-ray surveys which would be possible with the WFXT will detect huge samples of virialized objects spanning the mass range from sub-groups to the most massive clusters, and extending in redshift to beyond z=2. These samples will be an excellent dataset for carrying out many traditional cosmological tests using the cluster mass function and power spectrum. Uniquely, WFXT will be able not only to detect clusters but also to make detailed X-ray measurements for a large number of clusters and groups right from the survey data. Very high quality measurements of the cluster mass function and spatial correlation over a very wide range of masses, spatial scales, and redshifts, will be useful for expanding the cosmological discovery space, and in particular, in searching for departures from the "concordant" Lambda-CDM cosmological model. Finding such departures would have far-reaching implications on our understanding of the fundamental physics which governs the Universe.

We discuss the central role played by the X-ray study of hot baryons within galaxy clusters to reconstruct the assembly of cosmic structures and to trace the past history of star formation and accretion onto supermassive Black Holes (BHs). We shortly review the progress in this field contributed by the current generation of X-ray telescopes. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM: (a) When and how is entropy injected into the inter-galactic medium (IGM)? (b) What is the history of metal enrichment of the IGM? (c) What physical mechanisms determine the presence of cool cores in galaxy clusters? (d) How is the appearance of proto-clusters at z~2 related to the peak of star formation activity and BH accretion? We show that a highly efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. A similar survey can only be carried out with a Wide-Field X-ray Telescope (WFXT), which combines a high survey speed with a sharp PSF across the entire FoV. We emphasize the important synergies that WFXT will have with a number of future ground-based and space telescopes, covering from the radio to the X-ray bands. Finally, we discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

While the exceptional sensitivity of Chandra and XMM-Newton has resulted in revolutionary studies of the Galactic neighborhood in the soft (<10 keV) X-ray band, there are many open questions. We discuss these issues and how they would be addressed by very wide-area (> 100 sq. deg.) X-ray surveys.

In the study of the process of cosmic structure formation numerical simulations are crucial tools to interface observational data to theoretical models and to investigate issues otherwise unexplored. Enormous advances have been achieved in the last years thanks to the availability of sophisticated codes, now allowing to tackle the problem of cosmic structure formation and subsequent evolution by covering larger and larger dynamical ranges. Moreover, computational cosmology is the ideal interpretative framework for the overwhelming amount of new data from extragalactic surveys and from large sample of individual objects. The Workshop Novicosmo 2008 "The Impact of Simulations in Cosmology and Galaxy Formation' held in SISSA was aimed at providing the state-of-the-art on the latest numerical simulations in Cosmology and in Galaxy Formation. Particular emphasis was given to the implementation of new physical processes in simulation codes, to the comparison between different codes and numerical schemes and how to use best supercomputing facilities of the next generation. Finally, the impact on our knowledge on the Physics of the Universe brought by this new channel of investigation has also been focused. The Workshop was divided in three sections corresponding (roughly) to three main areas of study: Reionization and Intergalactic medium; Dark and Luminous matter in galaxies; Clusters of galaxies and Large scale Structures. This paper will provide i) a short resume' of the scientific results of the Workshop ii) the complete list of the talks and the instructions on how to retrieve the .pdf of the related (powerpoint) presentations iii) a brief presentation of the associated Exhibition "Space Art"

We present a method to recover mass profiles of galaxy clusters by combining data on thermal Sunyaev-Zeldovich (tSZ) and X-ray imaging, thereby avoiding to use any information on X-ray spectroscopy. This method, which represents a development of the geometrical deprojection technique presented in Ameglio et al. (2007), implements the solution of the hydrostatic equilibrium equation. In order to quantify the efficiency of our mass reconstructions, we apply our technique to a set of hydrodynamical simulations of galaxy clusters. We propose two versions of our method of mass reconstruction. Method 1 is completely model-independent, while Method 2 assumes instead the analytic mass profile proposed by Navarro et al. (1997) (NFW). We find that the main source of bias in recovering the mass profiles is due to deviations from hydrostatic equilibrium, which cause an underestimate of the mass of about 10 per cent at r_500 and up to 20 per cent at the virial radius. Method 1 provides a reconstructed mass which is biased low by about 10 per cent, with a 20 per cent scatter, with respect to the true mass profiles. Method 2 proves to be more stable, reducing the scatter to 10 per cent, but with a larger bias of 20 per cent, mainly induced by the deviations from equilibrium in the outskirts. To better understand the results of Method 2, we check how well it allows to recover the relation between mass and concentration parameter. When analyzing the 3D mass profiles we find that including in the fit the inner 5 per cent of the virial radius biases high the halo concentration. Also, at a fixed mass, hotter clusters tend to have larger concentration. Our procedure recovers the concentration parameter essentially unbiased but with a scatter of about 50 per cent.

We analyse the evolution of the Intergalactic Medium (IGM) by means of an extended set of large box size hydrodynamical simulations which include pre-heating. We focus on the properties of the z~2 Lyman-alpha forest and on the population of clusters and groups of galaxies at z=0. We investigate the distribution of voids in the Lyman-alpha flux and the entropy-temperature relation of galaxy groups, comparing the simulation results to recent data from high-resolution quasar spectra and from X-ray observations. Pre-heating is included through a simple phenomenological prescription, in which at z=4 the entropy of all gas particles, whose overdensity exceeds a threshold value delta_h is increased to a minimum value K_fl. While the entropy level observed in the central regions of galaxy groups requires a fairly strong pre-heating, with K_fl>100 keV cm^2, the void statistics of the Lyman-alpha forest impose that this pre-heating should take place only in relatively high-density regions, in order not to destroy the cold filaments that give rise to the forest. We conclude that any injection of non-gravitational energy in the diffuse baryons should avoid low-density regions at high redshift and/or take place at relatively low redshift.

We present results from two high-resolution hydrodynamical simulations of proto-cluster regions at z~2.1. The simulations have been compared to observational results for the socalled Spiderweb galaxy system, the core of a putative proto-cluster region at z = 2.16, found around a radio galaxy. The simulated regions have been chosen so as to form a poor cluster with M200~10^14 h-1 Msun (C1) and a rich cluster with M200~2x10^15 h-1 Msun (C2) at z = 0. The simulated proto-clusters show evidence of ongoing assembly of a dominating central galaxy. The stellar mass of the brightest cluster galaxy (BCG) of the C2 system is in excess with respect to observational estimates for the Spiderweb galaxy, with a total star formation rate which is also larger than indicated by observations. We find that the projected velocities of galaxies in the C2 cluster are consistent with observations, while those measured for the poorer cluster C1 are too low compared to the observed velocities. We argue that the Spiderweb complex resemble the high-redshift progenitor of a rich galaxy cluster. Our results indicate that the included supernovae feedback is not enough to suppress star formation in these systems, supporting the need of introducing AGN feedback. According to our simulations, a diffuse atmosphere of hot gas in hydrostatic equilibrium should already be present at this redshift, and enriched at a level comparable to that of nearby galaxy clusters. The presence of this gas should be detectable with future deep X-ray observations.

We study the degree to which non-radiative gas dynamics affects the merger histories of haloes along with subsequent predictions from a semi-analytic model (SAM) of galaxy formation. To this aim, we use a sample of dark matter only and non-radiative SPH simulations of four massive clusters. The presence of gas-dynamical processes (e.g. ram-pressure from the hot ICM) makes haloes more fragile in the runs which include gas. This results in a 25 per cent decrease in the total number of subhaloes at z = 0. The impact on the galaxy population predicted by SAMs is complicated by the presence of "orphan" galaxies, i.e. galaxies whose parent substructures are reduced below the resolution limit of the simulation. In the model employed in our study, these galaxies survive for a residual merging time that is computed using a variation of the Chandrasekhar formula. Due to ram-pressure stripping, haloes in gas simulations tend to be less massive than their counterparts in the dark matter simulations. The resulting merging times for satellite galaxies are then longer in these simulations. On the other hand, the presence of gas influences the orbits of haloes making them on average more circular and therefore reducing the estimated merging times with respect to the dark matter only simulation. This effect is particularly significant for the most massive satellites and is (at least in part) responsible for the fact that brightest cluster galaxies in runs with gas have stellar masses which are about 25 per cent larger than those obtained from dark matter only simulations. Our results show that gas-dynamics has only a marginal impact on the statistical properties of the galaxy population, but that its impact on the orbits and merging times of haloes strongly influences the assembly of the most massive galaxies.

Here we discuss the X-ray properties of clusters of galaxies optically selected in the Red-Sequence Cluster Survey (RCS) observed with the Chandra satellite, at redshifts 0.6 < z < 1.2. We intended to assess the evolutionary stage of optically selected high-z clusters of galaxies, performing a spectral analysis of the diffuse emission from their ICM. We also investigated the distribution of AGN in their surroundings. The background subtracted spectra were analyzed and fitted with a single temperature model to measure average ICM temperature, X-ray bolometric luminosity and Fe abundance within typical radii between 200 and 350 kpc. We also analyzed the point source number density and spatial distribution in the RCS clusters fields as a function of the X-ray flux. We detected emission for the majority of the clusters, except for three, for which we have only marginal detection at ~3 Sigma. We find that the normalization of the L-T relation for RCS clusters is a factor of ~2 lower than the one for X-ray selected clusters. We confirm that the Fe abundance in the detected objects is consistent with that of X-ray selected clusters at the same redshift. We also found an excess of low-luminosity AGN towards the center of the clusters.

We present a comparison between simulation results and X-ray observational data on the evolution of the metallicity of the intra-cluster medium (ICM). The simulations of galaxy clusters were performed with the Tree-SPH Gadget2 code that includes a detailed model of chemical evolution, by assuming three different shapes for the stellar initial mass function (IMF), namely the Salpeter (1955), Kroupa (2001) and Arimoto-Yoshii (1987) IMF. Our simulations predict significant radial gradients of the Iron abundance, which extend over the whole cluster virialized region. At larger radii, we do not detect any flattening of the metallicity profiles. As for the evolution of the ICM metal (Iron) abundance out to z=1, we find that it is determined by the combined action of (i) the sinking of already enriched gas, (ii) the ongoing metal production in galaxies and (iii) the locking of ICM metals in newborn stars. As a result, rather than suppressing the metallicity evolution, stopping star formation at z=1 has the effect of producing an even too fast evolution of the emission-weighted ICM metallicity with too high values at low redshift. Finally, we compare simulations with the observed rate of type-Ia supernovae per unit B-band luminosity (SnU_B). We find that our simulated clusters do not reproduce the decreasing trend of SnU_B at low redshift, unless star formation is truncated at z=1.

We present the work of an international team at the International Space Science Institute (ISSI) in Bern that worked together to review the current observational and theoretical status of the non-virialised X-ray emission components in clusters of galaxies. The subject is important for the study of large-scale hierarchical structure formation and to shed light on the "missing baryon" problem. The topics of the team work include thermal emission and absorption from the warm-hot intergalactic medium, non-thermal X-ray emission in clusters of galaxies, physical processes and chemical enrichment of this medium and clusters of galaxies, and the relationship between all these processes. One of the main goals of the team is to write and discuss a series of review papers on this subject. These reviews are intended as introductory text and reference for scientists wishing to work actively in this field. The team consists of sixteen experts in observations, theory and numerical simulations.

The study of the metal enrichment of the intra-cluster and inter-galactic media (ICM and IGM) represents a direct means to reconstruct the past history of star formation, the role of feedback processes and the gas-dynamical processes which determine the evolution of the cosmic baryons. In this paper we review the approaches that have been followed so far to model the enrichment of the ICM in a cosmological context. While our presentation will be focused on the role played by hydrodynamical simulations, we will also discuss other approaches based on semi-analytical models of galaxy formation, also critically discussing pros and cons of the different methods. We will first review the concept of the model of chemical evolution to be implemented in any chemo-dynamical description. We will emphasise how the predictions of this model critically depend on the choice of the stellar initial mass function, on the stellar life-times and on the stellar yields. We will then overview the comparisons presented so far between X-ray observations of the ICM enrichment and model predictions. We will show how the most recent chemo-dynamical models are able to capture the basic features of the observed metal content of the ICM and its evolution. We will conclude by highlighting the open questions in this study and the direction of improvements for cosmological chemo-dynamical models of the next generation.

Modern hydrodynamical simulations offer nowadays a powerful means to trace the evolution of the X-ray properties of the intra-cluster medium (ICM) during the cosmological history of the hierarchical build up of galaxy clusters. In this paper we review the current status of these simulations and how their predictions fare in reproducing the most recent X-ray observations of clusters. After briefly discussing the shortcomings of the self-similar model, based on assuming that gravity only drives the evolution of the ICM, we discuss how the processes of gas cooling and non-gravitational heating are expected to bring model predictions into better agreement with observational data. We then present results from the hydrodynamical simulations, performed by different groups, and how they compare with observational data. As terms of comparison, we use X-ray scaling relations between mass, luminosity, temperature and pressure, as well as the profiles of temperature and entropy. The results of this comparison can be summarised as follows: (a) simulations, which include gas cooling, star formation and supernova feedback, are generally successful in reproducing the X-ray properties of the ICM outside the core regions; (b) simulations generally fail in reproducing the observed ``cool core'' structure, in that they have serious difficulties in regulating overcooling, thereby producing steep negative central temperature profiles. This discrepancy calls for the need of introducing other physical processes, such as energy feedback from active galactic nuclei, which should compensate the radiative losses of the gas with high density, low entropy and short cooling time, which is observed to reside in the innermost regions of galaxy clusters.

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.

(Abridged)The X-ray measurements of the ICM metallicity are becoming more frequent due to the availability of powerful X-ray telescope with excellent spatial and spectral resolutions. The information which can be extracted from the measurements of the alpha-elements, like Oxygen, Magnesium and Silicon with respect to the Iron abundance is extremely important to better understand the stellar formation and its evolutionary history. In this paper we investigate possible source of bias connected to the plasma physics when recovering metal abundances from X-ray spectra. To do this we analyze 6 simulated galaxy clusters processed through the new version of our X-ray MAp Simulator, which allows to create mock XMM-Newton EPIC MOS1 and MOS2 observations. By comparing the spectroscopic results to the input values we find that: i) Fe is recovered with high accuracy for both hot (T>3 keV) and cold (T<2 keV) systems; at intermediate temperatures, however, we find a systematic overestimate which depends on the number counts; ii) O is well recovered in cold clusters, while in hot systems its measure may overestimate by a factor up to 2-3; iii) Being a weak line, the measurement of Mg is always difficult; despite of this, for cold systems (T<2 keV) we do not find any systematic behavior, while for very hot systems (T>5 keV) the spectroscopic measurement may be strongly overestimated up to a factor of 4; iv) Si is well recovered for all the clusters in our sample. We investigate in detail the nature of the systematic effects and biases found. We conclude that they are mainly connected with the multi-temperature nature of the projected observed spectra and to the intrinsic limitation of the XMM-Newton EPIC spectral resolution that does not always allow to disentangle among the emission lines produced by different elements.

In order to study the process of cooling in dark-matter (DM) halos and assess how well simple models can represent it, we run a set of radiative SPH hydrodynamical simulations of isolated halos, with gas sitting initially in hydrostatic equilibrium within Navarro-Frenk-White (NFW) potential wells. [...] After having assessed the numerical stability of the simulations, we compare the resulting evolution of the cooled mass with the predictions of the classical cooling model of White & Frenk and of the cooling model proposed in the MORGANA code of galaxy formation. We find that the classical model predicts fractions of cooled mass which, after about two central cooling times, are about one order of magnitude smaller than those found in simulations. Although this difference decreases with time, after 8 central cooling times, when simulations are stopped, the difference still amounts to a factor of 2-3. We ascribe this difference to the lack of validity of the assumption that a mass shell takes one cooling time, as computed on the initial conditions, to cool to very low temperature. [...] The MORGANA model [...] better agrees with the cooled mass fraction found in the simulations, especially at early times, when the density profile of the cooling gas is shallow. With the addition of the simple assumption that the increase of the radius of the cooling region is counteracted by a shrinking at the sound speed, the MORGANA model is also able to reproduce for all simulations the evolution of the cooled mass fraction to within 20-50 per cent, thereby providing a substantial improvement with respect to the classical model. Finally, we provide a very simple fitting function which accurately reproduces the cooling flow for the first ~10 central cooling times. [Abridged]

We present two non-parametric deprojection methods aimed at recovering the three-dimensional density and temperature profiles of galaxy clusters from spatially resolved thermal Sunyaev-Zeldovich (tSZ) and X-ray surface brightness maps, thus avoiding the use of X-ray spectroscopic data. In both methods, clusters are assumed to be spherically symmetric and modeled with an onion-skin structure. The first method follows a direct geometrical approach. The second method is based on the maximization of a single joint (tSZ and X-ray) likelihood function, which allows one to fit simultaneously the two signals by following a Monte Carlo Markov Chain approach. These techniques are tested against a set of cosmological simulations of clusters, with and without instrumental noise. We project each cluster along the three orthogonal directions defined by the principal axes of the momentum of inertia tensor. This enables us to check any bias in the deprojection associated to the cluster elongation along the line of sight. After averaging over all the three projection directions, we find an overall good reconstruction, with a small (<~10 per cent) overestimate of the gas density profile. This turns into a comparable overestimate of the gas mass within the virial radius, which we ascribe to the presence of residual gas clumping. Apart from this small bias the reconstruction has an intrinsic scatter of about 5 per cent, which is dominated by gas clumpiness. Cluster elongation along the line of sight biases the deprojected temperature profile upwards at r<~0.2r_vir and downwards at larger radii. A comparable bias is also found in the deprojected temperature profile. Overall, this turns into a systematic underestimate of the gas mass, up to 10 percent. (Abridged)

Using the results of a cosmological hydrodynamical simulation of the concordance LambdaCDM model, we study the global properties of the Sunyaev-Zel'dovich (SZ) effects, both considering the thermal (tSZ) and the kinetic (kSZ) component. The simulation follows gravitation and gas dynamics and includes also several physical processes that affect the baryonic component, like a simple reionization scenario, radiative cooling, star formation and supernova feedback. Starting from the outputs of the simulation we create mock maps of the SZ signals due to the large structures of the Universe integrated in the range 0 < z < 6. We predict that the Compton y-parameter has an average value of (1.19 +/- 0.32) 10^-6 and is lognormally distributed in the sky; half of the whole signal comes from z < 1 and about 10 per cent from z > 2. The Doppler b-parameter shows approximately a normal distribution with vanishing mean value and a standard deviation of 1.6 10^-6, with a significant contribution from high-redshift (z > 3) gas. We find that the tSZ effect is expected to dominate the primary CMB anisotropies for l >~ 3000 in the Rayleigh-Jeans limit, while interestingly the kSZ effect dominates at all frequencies at very high multipoles (l >~ 7 10^4). We also analyse the cross-correlation between the two SZ effects and the soft (0.5-2 keV) X-ray emission from the intergalactic medium and we obtain a strong correlation between the three signals, especially between X-ray emission and tSZ effect (r_l ~ 0.8-0.9) at all angular scales.

We present the analysis of the temperature and metallicity profiles of 12 galaxy clusters in the redshift range 0.1--0.3 selected from the Chandra archive with at least ~20,000 net ACIS counts and kT>6 keV. We divide the sample between 7 Cooling-Core (CC) and 5 Non-Cooling-Core (NCC) clusters according to their central cooling time. We find that single power-laws can describe properly both the temperature and metallicity profiles at radii larger than 0.1 r_180 in both CC and NCC systems, showing the NCC objects steeper profiles outwards. A significant deviation is only present in the inner 0.1 r_180. We perform a comparison of our sample with the De Grandi & Molendi BeppoSAX sample of local CC and NCC clusters, finding a complete agreement in the CC cluster profile and a marginally higher value (at ~1sigma) in the inner regions of the NCC clusters. The slope of the power-law describing kT(r) within 0.1 r_180 correlates strongly with the ratio between the cooling time and the age of the Universe at the cluster redshift, being the slope >0 and tau_c/tau_age<=0.6 in CC systems.

We present cosmological hydrodynamical simulations of galaxy clusters aimed at studying the process of metal enrichment of the intra--cluster medium (ICM). These simulations have been performed by implementing a detailed model of chemical evolution in the Tree-SPH \gd code. This model allows us to follow the metal release from SNII, SNIa and AGB stars, by properly accounting for the lifetimes of stars of different mass, as well as to change the stellar initial mass function (IMF), the lifetime function and the stellar yields. As such, our implementation of chemical evolution represents a powerful instrument to follow the cosmic history of metal production. The simulations presented here have been performed with the twofold aim of checking numerical effects, as well as the impact of changing the model of chemical evolution and the efficiency of stellar feedback.

The largest uncertainty for cosmological studies using clusters of galaxies is introduced by our limited knowledge of the statistics of galaxy cluster structure, and of the scaling relations between observables and cluster mass. To improve on this situation we have started an XMM-Newton Large Programme for the in-depth study of a representative sample of 33 galaxy clusters, selected in the redshift range z=0.055 to 0.183 from the REFLEX Cluster Survey, having X-ray luminosities above 0.4 X 10^44 h_70^-2 erg s^-1 in the 0.1 - 2.4 keV band. This paper introduces the sample, compiles properties of the clusters, and provides detailed information on the sample selection function. We describe the selection of a nearby galaxy cluster sample that makes optimal use of the XMM-Newton field-of-view, and provides nearly homogeneous X-ray luminosity coverage for the full range from poor clusters to the most massive objects in the Universe. For the clusters in the sample, X-ray fluxes are derived and compared to the previously obtained fluxes from the ROSAT All-Sky Survey. We find that the fluxes and the flux errors have been reliably determined in the ROSAT All-Sky Survey analysis used for the REFLEX Survey. We use the sample selection function documented in detail in this paper to determine the X-ray luminosity function, and compare it with the luminosity function of the entire REFLEX sample. We also discuss morphological peculiarities of some of the sample members. The sample and some of the background data given in this introductory paper will be important for the application of these data in the detailed studies of cluster structure, to appear in forthcoming publications.

We present a Chandra analysis of the X-ray spectra of 56 clusters of galaxies at $z>0.3$, which cover a temperature range of $3> kT > 15$ keV. Our analysis is aimed at measuring the iron abundance in the ICM out to the highest redshift probed to date. We find that the emission-weighted iron abundance measured within $(0.15-0.3) R_{vir}$ in clusters below 5 keV is, on average, a factor of $\sim2$ higher than in hotter clusters, following $Z(T)\simeq 0.88 T^{-0.47} Z_\odot$, which confirms the trend seen in local samples. We made use of combined spectral analysis performed over five redshift bins at $0.3> z > 1.3$ to estimate the average emission weighted iron abundance. We find a constant average iron abundance $Z_{Fe}\simeq 0.25 Z_\odot$ as a function of redshift, but only for clusters at $z>0.5$. The emission-weighted iron abundance is significantly higher ($Z_{Fe}\simeq0.4 Z_\odot$) in the redshift range $z\simeq0.3-0.5$, approaching the value measured locally in the inner $0.15 R_{vir}$ radii for a mix of cool-core and non cool-core clusters in the redshift range $0.1<z<0.3$. The decrease in $Z_{Fe}$ with $z$ can be parametrized by a power law of the form $\sim(1+z)^{-1.25}$. The observed evolution implies that the average iron content of the ICM at the present epoch is a factor of $\sim2$ larger than at $z\simeq 1.2$. We confirm that the ICM is already significantly enriched ($Z_{Fe}\simeq0.25 Z_\odot$) at a look-back time of 9 Gyr. Our data provide significant constraints on the time scales and physical processes that drive the chemical enrichment of the ICM.

We study the origin of the diffuse stellar component (DSC) in 117 galaxy clusters extracted from a cosmological hydrodynamical simulation. We identify all galaxies present in the simulated clusters at 17 output redshifts, and then build the family trees for all the z=0 cluster galaxies. Then for each diffuse star particle identified at z=0, we look for the galaxy to which it once belonged at an earlier redshift, thus linking the presence of the diffuse stellar component to the galaxy formation history. The main results of our analysis are: (i) On average, half of the DSC star particles comes from galaxies associated with the family tree of the most massive galaxy (BCG), one quarter comes from the family trees of other massive galaxies, and the remaining quarter from dissolved galaxies. I.e., the formation of the DSC is parallel to the build-up of the BCG and other massive galaxies. (ii) Most DSC star particles become unbound during mergers in the formation history of the BGCs and of other massive galaxies, independent of cluster mass. Our results suggest that the tidal stripping mechanism is responsible only for a minor fraction of the DSC. (iii) At cluster radii larger than 250 kpc/h, the DSC fraction from the BCG is reduced and the largest contribution comes from the other massive galaxies; in the cluster outskirts, galaxies of all masses contribute to the DSC. (iv) The DSC does not have a preferred redshift of formation: however, most DSC stars are unbound at z<1. (v) The amount of DSC stars at z=0 does not corre late strongly with the global dynamical history of clusters, and increases weakly with cluster mass. [abridged]

We present a Chandra analysis of the X-ray spectra of 56 clusters of galaxies at z>0.3, which cover a temperature range 3<kT<15 keV. Our analysis is aimed at measuring the iron abundance in the intra-cluster medium (ICM) out to the highest redshift probed to date. We made use of combined spectral analysis performed over five redshift bins at 0.3<z<1.3 to estimate the average emission weighted iron abundance. We find that the emission weighted iron abundance measured within (0.15-0.3)R_vir in clusters below 5 keV is, on average, a factor of ~2 higher than in hotter clusters, following Z(T) = 0.88 T^{-0.47} Z_o, confirming the trend seen in local samples. We also find a constant average iron abundance Z_{Fe}~0.25Z_o as a function of redshift, but only for clusters at z>0.5. The emission weighted iron abundance is significantly higher (Z_{Fe}~0.4 Z_o) in the redshift range z~0.3-0.5, approaching the value measured locally in the inner 0.15 R_{vir} radii for a mix of cool-core and non cool-core clusters at 0.1<z<0.3. The decrease in metallicity with redshift can be parametrized by a power law of the form ~(1+z)^{-1.25}. The observed evolution implies that the average iron content of the ICM at present epoch is a factor of ~2 larger than at z~1.2. We confirm that the ICM is already significantly enriched (Z_{Fe} ~ 0.25 Z_o) at a look-back time of 9 Gyr. Our data provide significant constraints on the time scales and physical processes which drive the chemical enrichment of the ICM. (abridged)

Using a set of hydrodynamical simulations of 9 galaxy clusters with masses in the range 1.5 10^{14} M_sun < M_vir < 3.4 10^{15} M_sun, we have studied the density, temperature and X-ray surface brightness profiles of the intracluster medium in the regions around the virial radius. We have analyzed the profiles in the radial range well above the cluster core, the physics of which are still unclear and matter of tension between simulated and observed properties, and up to the virial radius and beyond, where present observations are unable to provide any constraints. We have modeled the radial profiles between 0.3 R_200 and 3 R_200 with power laws with one index, two indexes and a rolling index. The simulated temperature and [0.5-2] keV surface brightness profiles well reproduce the observed behaviours outside the core. The shape of all these profiles in the radial range considered depends mainly on the activity of the gravitational collapse, with no significant difference among models including extraphysics. The profiles steepen in the outskirts, with the slope of the power-law fit that changes from -2.5 to -3.4 in the gas density, from -0.5 to -1.8 in the gas temperature, and from -3.5 to -5.0 in the X-ray soft surface brightness. We predict that the gas density, temperature and [0.5-2] keV surface brightness values at R_200 are, on average, 0.05, 0.60, 0.008 times the measured values at 0.3 R_200. At 2 R_200, these values decrease by an order of magnitude in the gas density and surface brightness, by a factor of 2 in the temperature, putting stringent limits on the detectable properties of the intracluster-medium (ICM) in the virial regions.

We present a Chandra analysis of the X-ray spectra of 56 clusters of galaxies at z>0.3, which cover a temperature range of 3>kT>15 keV. Our analysis is aimed at measuring the iron abundance in the ICM out to the highest redshift probed to date. We find that the emission-weighted iron abundance measured within (0.15-0.3)R_vir in clusters below 5 keV is, on average, a factor of ~2 higher than in hotter clusters, following Z(T)~0.88T^-(0.47)Z_o, which confirms the trend seen in local samples. We made use of combined spectral analysis performed over five redshift bins at 0.3>z>1.3 to estimate the average emission weighted iron abundance. We find a constant average iron abundance Z_Fe~0.25Z_o as a function of redshift, but only for clusters at z>0.5. The emission-weighted iron abundance is significantly higher (Z_Fe~0.4Z_o) in the redshift range z~0.3-0.5, approaching the value measured locally in the inner 0.15R_vir radii for a mix of cool-core and non cool-core clusters in the redshift range 0.1<z<0.3. The decrease in Z_Fe with redshift can be parametrized by a power law of the form ~(1+z)^(-1.25). The observed evolution implies that the average iron content of the ICM at the present epoch is a factor of ~2 larger than at z=1.2. We confirm that the ICM is already significantly enriched (Z_Fe~0.25Z_o) at a look-back time of 9 Gyr. Our data provide significant constraints on the time scales and physical processes that drive the chemical enrichment of the ICM.

The high end of the stellar mass function of galaxies is observed to have little evolution since z~1. This represents a stringent constraint for merger--based models, aimed at explaining the evolution of the most massive galaxies in the concordance LambdaCDM cosmology. In this Letter we show that it is possible to remove the tension between the above observations and model predictions by allowing a fraction of stars to be scattered to the Diffuse Stellar Component (DSC) of galaxy clusters at each galaxy merger, as recently suggested by the analysis of N-body hydrodynamical simulations. To this purpose, we use the MORGANA model of galaxy formation in a minimal version, in which gas cooling and star formation are switched off after z=1. In this way, any predicted evolution of the galaxy stellar mass function is purely driven by mergers. We show that, even in this extreme case, the predicted degree of evolution of the high end of the stellar mass function is larger than that suggested by data. Assuming instead that a significant fraction, ~30 per cent, of stars are scattered in the DSC at each merger event, leads to a significant suppression of the predicted evolution, in better agreement with observational constraints, while providing a total amount of DSC in clusters which is consistent with recent observational determinations.

A study of the structural and scaling properties of the temperature distribution of the hot, X-ray emitting intra-cluster medium of galaxy clusters, and its dependence on dynamical state, can give insights into the physical processes governing the formation and evolution of structure. We analyse the X-ray temperature profiles from XMM-Newton observations of 15 nearby (z < 0.2) clusters, drawn from a statistically representative sample. The clusters cover a temperature range from 2.5 keV to 8.5 keV, and present a variety of X-ray morphologies. We derive accurate projected temperature profiles to ~ 0.5 R_200, and compare structural properties (outer slope, presence of cooling core) with a quantitative measure of the X-ray morphology as expressed by power ratios. We also compare the results to recent cosmological numerical simulations. Once the temperature profiles are scaled by an average cluster temperature (excluding the central region) and the estimated virial radius, the profiles generally decline in the region 0.1 R_200 < R < 0.5 R_200. The central regions show the largest scatter, attributable mostly to the presence of cool core clusters. There is good agreement with numerical simulations outside the core regions. We find no obvious correlations between power ratio and outer profile slope. There may however be a weak trend with the existence of a cool core, in the sense that clusters with a central temperature decrement appear to be slightly more regular. The present results lend further evidence to indicate that clusters are a regular population, at least outside the core region.

We present a study of the galaxy population predicted by hydrodynamical simulations for a set of 19 galaxy clusters based on the GADGET-2 Tree+SPH code. These simulations include gas cooling, star formation, a detailed treatment of stellar evolution and chemical enrichment, as well as SN energy feedback in the form of galactic winds. We compute the spectro-photometric properties of the simulated galaxies. All simulations have been performed for two choices of the stellar initial mass function: a standard Salpeter IMF, and a top-heavier IMF. Several of the observational properties of the galaxy population in nearby clusters are reproduced fairly well by simulations. A Salpeter IMF is successful in accounting for the slope and the normalization of the color-magnitude relation for the bulk of the galaxy population. Simulated clusters have a relation between mass and optical luminosity which generally agrees with observations, both in normalization and slope. We find that galaxies are generally bluer, younger and more star forming in the cluster outskirts, thus reproducing the observational trends. However, simulated clusters have a total number of galaxies which is significantly smaller than the observed one, falling short by about a factor 2-3. Finally, the brightest cluster galaxies are always predicted to be too massive and too blue, when compared to observations, due to gas overcooling in the core cluster regions, even in the presence of a rather efficient SN feedback.

In this Chapter I review the role that galaxy clusters play as tools to constrain cosmological parameters. I will concentrate mostly on the application of the mass function of galaxy clusters, while other methods, such as that based on the baryon fraction, are covered by other Chapters of the book. Since most of the cosmological applications of galaxy clusters rely on precise measurements of their masses, a substantial part of my Lectures concentrates on the different methods that have been applied so far to weight galaxy clusters. I provide in Section 2 a short introduction to the basics of cosmic structure formation. In Section 3 I describe the Press--Schechter (PS) formalism to derive the cosmological mass function, then discussing extensions of the PS approach and the most recent calibrations from N--body simulations. In Section 4 I review the methods to build samples of galaxy clusters at different wavelengths. Section 5 is devoted to the discussion of different methods to derive cluster masses. In Section 6 I describe the cosmological constraints, which have been obtained so far by tracing the cluster mass function with a variety of methods. Finally, I describe in Section 7 the future perspectives for cosmology with galaxy clusters and the challenges for clusters to keep playing an important role in the era of precision cosmology.

We study the efficiency and reliability of cluster mass estimators that are based on the projected phase-space distribution of galaxies in a cluster region. To this aim, we analyse a data-set of 62 clusters extracted from a concordance LCDM cosmological hydrodynamical simulation. Galaxies (or Dark Matter particles) are first selected in cylinders of given radius (from 0.5 to 1.5 Mpc/h) and ~200 Mpc/h length. Cluster members are then identified by applying a suitable interloper removal algorithm. Two cluster mass estimators are considered: the virial mass estimator (Mvir), and a mass estimator (Msigma) based entirely on the cluster velocity dispersion estimate. Mvir overestimates the true mass by ~10%, and Msigma underestimates the true mass by ~15%, on average, for sample sizes of > 60 cluster members. For smaller sample sizes, the bias of the virial mass estimator substantially increases, while the Msigma estimator becomes essentially unbiased. The dispersion of both mass estimates increases by a factor ~2 as the number of cluster members decreases from ~400 to ~20. The bias in the Mvir estimates is reduced in clusters without significant evidence for subclustering, and when only early-type galaxies are selected. Radially-dependent incompleteness can drastically affect Mvir estimates, but leaves the Msigma estimates almost unaffected. Other observational effects, like centering and velocity errors, and different observational apertures, have little effect on the mass estimates. (Abridged)

The angular-diameter distance D_A of a galaxy cluster can be measuread by combining its X-ray emission with the cosmic microwave background fluctution due to the Sunyaev-Zeldovich effect. The application of this distance indicator usually assumes that the cluster is spherically symmetric, the gas is distributed according to the isothermal beta-model, and the X-ray temperature is an unbiased measure of the electron temperature. We test these assumptions with galaxy clusters extracted from an extended set of cosmological N-body/hydrodynamical simulations of a LCDM concordance cosmology, which include the effect of radiative cooling, star formation and energy feedback from supernovae. We find that, due to the steep temperature gradients which are present in the central regions of simulated clusters, the assumption of isothermal gas leads to a significant underestimate of D_A. This bias is efficiently corrected by using the polytropic version of the beta-model to account for the presence of temperature gradients. In this case, once irregular clusters are removed, the correct value of D_A is recovered with a ~ 5 per cent accuracy on average, with a ~ 20 per cent intrinsic scatter due to cluster asphericity. This result is valid when using either the electron temperature or a spectroscopic-like temperature. When using instead the emission-weighted definition for the temperature of the simulated clusters, D_A is biased low by \~ 20 per cent. We discuss the implications of our results for an accurate determination of the Hubble constant H_0 and of the density parameter Omega_m. We find that H_0 can be potentially recovered with exquisite precision, while the resulting estimate of Omega_m, which is unbiased, has typical errors Delta(Omega_m) ~ 0.05.

We examine the systematics affecting the X-ray mass estimators applied to a set of five simulated galaxy clusters. They have been processed through the X-ray Map Simulator, X-MAS, to provide Chandra-like long exposures that are analyzed to reconstruct the gas temperature, density, and mass profiles used as input. We find that at R_2500 the mass profile obtained via a direct application of the hydrostatic equilibrium equation is consistent within 1 sigma with the actual mass; although we notice this estimator shows high statistical errors due to high level of Chandra background. Instead, the poorness of the beta-model in describing the gas density profile makes the evaluated masses to be underestimated by \sim 40 per cent with respect to the true mass, both with an isothermal and a polytropic temperature profile. We also test ways to recover the mass by adopting an analytic mass model, such as those proposed by Navarro et al. (1997) and Rasia et al. (2004), and fitting the temperature profile expected from the hydrostatic equilibrium equation to the observed one. These methods and the one of the hydrostatic equilibrium equation provide a more robust mass estimation than the ones based on the beta-model. In the present work the main limitation for a precise mass reconstruction is to ascribe to the relatively high level of the background chosen to reproduce the Chandra one. After artificially reducing it by a factor of 100, we find that the estimated mass significantly underestimates the true mass profiles. This is manly due (i) to the neglected contribution of the gas bulk motions to the total energy budget and (ii) to the bias towards lower values of the X-ray temperature measurements because of the complex thermal structure of the emitting plasma.

We present a detailed X-ray spectral analysis of the sources in the 1Ms catalog of the Chandra Deep Field South (CDFS) taking advantage of optical spectroscopy and photometric redshifts for 321 sources. As a default spectral model, we adopt a power law with slope Gamma with an intrinsic redshifted absorption N_H, a fixed Galactic absorption and an unresolved Fe emission line. For 82 X-ray bright sources, we perform the X-ray spectral analysis leaving both Gamma and N_H free. The weighted mean value is <Gamma>~ 1.75+-0.02, with an intrinsic dispersion of sigma~0.30. We do not find hints of a correlation between Gamma and the N_H. We detect the presence of a scattered component at soft energies in 8 sources, and a pure reflection spectrum, typical of Compton-thick AGN, in 14 sources (Compton-thick AGN candidates). The intrinsic N_H distribution shows a lognormal shape, peaking around log(N_H)~23.1 and with sigma~1.1. We find that the fraction of absorbed sources (with N_H>10^{22} cm^{-2}) in the sample is constant (at the level of about 75%) or moderately increasing with redshift. Finally, we compare the optical classification to the X-ray spectral properties, confirming that the correspondence of unabsorbed (absorbed) X-ray sources to optical Type I (Type II) AGN is accurate for at least 80% of the sources with spectral identification (1/3 of the total X-ray sample).

We study the properties of the diffuse X-ray background by using the results of a cosmological hydrodynamical simulation of the concordance LambdaCDM model. The simulation follows gravitational and gas dynamics and includes a treatment of physical processes like radiative cooling, star formation and supernova feedback. From the simulation outputs, we produce a set of two-dimensional maps of the intergalactic medium X-ray emission integrated over redshift. We find that the signal in the soft (0.5-2 keV) band is lognormally distributed with a mean intensity of about 4 10^-12 erg s^-1 cm^-2 deg^-2; approximately 40 per cent of the emission originates from warm-hot gas (defined as baryons with 10^5<T<10^7 K), and 90 per cent comes from structures at z<0.9. Since the spectrum is soft, being mostly provided by the intergalactic medium at low temperature, the total mean intensity in the hard (2-10 keV) X-ray band is smaller by a factor of about 4. In order to constrain the physical processes included in our simulation, we compare our results with the observed upper limit (1.2 +/- 0.3) 10^-12 erg s^-1 cm^-2 deg^-2 of the soft X-ray emission due to diffuse gas. To this purpose, we remove the contributions of observable extended objects (groups and clusters of galaxies) from the simulated maps by adopting different detectability criteria which are calibrated on the properties of systems at intermediate redshifts observed by Chandra. We show that the simulated diffuse soft X-ray emission is consistent with the present observed upper limit. However, if future measurements will decrease the level of the unresolved X-ray background by a factor of two, a more efficient feedback mechanism should be required to suppress the soft emission of the gas residing within filaments and group-size haloes.

We discuss an extended set of Tree+SPH simulations of galaxy clusters, with the goal of investigating the interplay between numerical resolution effects and star-formation/feedback processes. The simulated clusters span the mass range (0.1-2.3) 10^{15}Msun/h, with mass resolution varying by several decades. At the highest achieved resolution, we resolve the virial region of a Virgo-like cluster with more than 2 million gas particles and with at least as many dark-matter (DM) particles. Our resolution study confirms that, in the absence of an efficient feedback mechanism, runaway cooling leads to about 35 per cent of baryons in clusters to be locked up in long lived stars at our highest resolution, with no evidence of convergence. However, including feedback causes the fraction of cooled baryons to converge at about 15 per cent already at modest resolution. Feedback also stabilizes other gas-related quantities, such as radial profiles of entropy, gas density and temperature, against variations due to changes in resolution. We also investigate the influence of the gravitational force softening length, and that of numerical heating of the gas induced by two-body encounters between DM and lighter gas particles. We show that simulations where more DM than gas particles are used, show a significantly enhanced efficiency of star formation at z>3. Our results are important for establishing and delineating the regime of numerical reliability of the present generation of hydrodynamical simulations of galaxy clusters.

We study the baryon mass fraction in a set of hydrodynamical simulations of galaxy clusters performed using the Tree+SPH code GADGET-2. We investigate the dependence of the baryon fraction upon the radiative cooling, star formation, feedback through galactic winds, conduction and redshift. Both the cold stellar component and the hot X-ray emitting gas have narrow distributions that, at large cluster-centric distances r>R500, are nearly independent of the physics included in the simulations. Only the non-radiative runs reproduce the gas fraction inferred from observations of the inner regions (r ~ R2500) of massive clusters. When cooling is turned on, the excess star formation is mitigated by the action of galactic winds, but yet not by the amount required by observational data. The baryon fraction within a fixed overdensity increases slightly with redshift, independent of the physical processes involved in the accumulation of baryons in the cluster potential well. In runs with cooling and feedback, the increase in baryons is associated with a larger stellar mass fraction that arises at high redshift as a consequence of more efficient gas cooling. For the same reason, the gas fraction appears less concentrated at higher redshift. We discuss the possible cosmological implications of our results and find that two assumptions generally adopted, (1) mean value of Yb = fb / (Omega_b/Omega_m) not evolving with redshift, and (2) a fixed ratio between f_star and f_gas independent of radius and redshift, might not be valid. In the estimate of the cosmic matter density parameter, this implies some systematic effects of the order of Delta Omega_m/Omega_m < +0.15 for non-radiative runs and Delta Omega_m/Omega_m ~ +0.05 and < -0.05 for radiative simulations.

The amount and nature of dark energy (DE) can be tightly constrained by measuring the spatial correlation features and evolution of a sample of ~ 100,000 galaxy clusters over the redshift range 0<z < 1.5. Such an X-ray survey will discover all collapsed structures with mass above 3.5e14 solar masss at redshifts z<2 (i.e. the full range where such objects are expected) in the high Galactic latitude sky. Above this mass threshold the tight correlations between X-ray observables and mass allow direct interpretation of the data. We describe the constraints on Dark Energy that can be inferred from such a survey, using powerful self-calibration techniques to relate the X-ray observables (luminosity and temperature) to the underlying mass.

We present the results from the dynamical analysis of the cluster of galaxies RX J0152-1357, which shows a complex structure in its X-ray emission, with two major clumps in the central region, and a third clump in the Eastern region. Our analysis is based on redshift data for 187 galaxies. We compute the line-of-sight velocity dispersion of galaxies, sigma_V=1322^+74_-68 kms^-1, which is significantly larger than what is expected in the case of a relaxed cluster with an observed X-ray temperature of 5-6 keV. We find evidence that this cluster is far from dynamical equilibrium, as shown by the non Gaussianity of the velocity distribution, the presence of a velocity gradient and significant substructure. Our analysis shows that the high value of sigma_V is due to the complex structure of RX J0152-1357, i.e. to the presence of three galaxy clumps of different mean velocity. Using optical data we detect a low-velocity clump (with sigma_V=300-500 kms^-1) in the central SW region and a high--velocity clump (with sigma_V sim 700 kms^-1) in the E region, well corresponding to the SW and E peaks detected in the X-ray emission. The central NE X-ray peak is associated to the main galaxy structure with a velocity which is intermediate between those of the other two clumps and sigma_V sim 900 kms^-1. The mass of the whole system within 2 Mpc is estimated to lie in the range (1.2-2.2)X10^{15} M_sun. Analytic calculations based on the two-body model indicate that the system is most likely bound, currently undergoing merging. In particular, we suggest that the SW clump is not a small group, but rather the dense cluster-core of a massive cluster.

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

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. Here we show that, if the intra-cluster medium is thermally complex, T_ew fails at reproducing T_spec. We propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we also find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and beta-model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X--ray temperature function could bias low sigma_8 by 10-20 per cent. This may alleviate the tension between the value of sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

We present the results of an extensive spectroscopic survey of RX J0152.7-1357, one of the most massive distant clusters of galaxies known. Multi-object spectroscopy, carried out with FORS1 and FORS2 on the ESO Very Large Telescope (VLT), has allowed us to measure more than 200 redshifts in the cluster field and to confirm 102 galaxies as cluster members. The mean redshift of the cluster is $z=0.837 \pm 0.001$ and we estimate the velocity dispersion of the overall cluster galaxy distribution to be $\sim 1600 \mathrm{km \ s^{-1}}$. The distribution of cluster members is clearly irregular, with two main clumps that follow the X-ray cluster emission mapped by Chandra. A third clump of galaxies to the east of the central structure and at the cluster redshift has also been identified. The two main clumps have velocity dispersions of $\sim919$ and $\sim737 \mathrm{km s^{-1}}$ respectively, and the peculiar velocity of the two clumps suggests that they will merge into a single more massive cluster. A segregation in the star formation activity of the member galaxies is observed. All star forming galaxies are located outside the high-density peaks, which are populated only by passive galaxies. A population of red galaxies (belonging to the cluster red sequence) with clear post-starburst spectral features and [OII] ($\lambda$3727) emission lines is observed in the outskirts of the cluster. Two AGNs, which were previously confused with the diffuse X-ray emission from the intracluster medium in ROSAT and BeppoSAX observations, are found to be cluster members.

The fluctuations in the Cosmic Microwave Background (CMB) intensity due to the Sunyaev-Zeldovich (SZ) effect are the sum of a thermal and a kinetic contribution. Separating the two components to measure the peculiar velocity of galaxy clusters requires radio and microwave observations at three or more frequencies, and knowledge of the temperature T_e of the intracluster medium weighted by the electron number density. To quantify the systematics of this procedure, we extract a sample of 117 massive clusters at redshift z=0 from an N-body hydrodynamical simulation, with 2x480^3 particles, of a cosmological volume 192 Mpc/h on a side of a flat Cold Dark Matter model with Omega_0=0.3 and Lambda=0.7. Our simulation includes radiative cooling, star formation and the effect of feedback and galactic winds from supernovae. We find that (1) our simulated clusters reproduce the observed scaling relations between X-ray and SZ properties; (2) bulk flows internal to the intracluster medium affect the velocity estimate by less than 200 km/s in 93 per cent of the cases; (3) using the X-ray emission weighted temperature, as an estimate of T_e, can overestimate the peculiar velocity by 20-50 per cent, if the microwave observations do not spatially resolve the cluster. For spatially resolved clusters, the assumptions on the spatial distribution of the ICM, required to separate the two SZ components, still produce a velocity overestimate of 10-20 per cent, even with an unbiased measure of T_e. Thanks to the large size of our cluster samples, these results set a robust lower limit of 200 km/s to the systematic errors that will affect upcoming measures of cluster peculiar velocities with the SZ effect.

We investigate the spatial clustering of X-ray selected sources in the two deepest X-ray fields to date, namely the 2Msec Chandra Deep Field North (CDFN) and the 1Msec Chandra Deep Field South (CDFS). The projected correlation function w(r_p), measured on scales ~0.2-10 h^-1 Mpc for a sample of 240 sources with spectroscopic redshift in the CDFN and 124 sources in the CDFS at a median redshift of z~0.8, is used to constrain the amplitude and slope of the real space correlation function xi(r)=(r/r0)^-gamma. The clustering signal is detected at high confidence (>~ 7 sigma) in both fields. The amplitude of the correlation is found to be significantly different in the two fields, the correlation length r0 being 8.6 +- 1.2 h^-1 Mpc in the CDFS and 4.2 +- 0.4 h^-1 Mpc in the CDFN, while the correlation slope gamma is found to be flat in both fields: gamma=1.33 +- 0.11 in the CDFS and gamma=1.42 +- 0.07 in the CDFN (a flat Universe with Omega_m=0.3 and Omega_L=0.7 is assumed; 1 sigma Poisson error estimates are considered). The correlation function has been also measured separately for sources classified as AGN or galaxies. In both fields AGN have a median redshift of z~0.9 and a median 0.5-10 keV luminosity of L_x~10^43 erg s^-1, i.e. they are generally in the Seyfert luminosity regime. As in the case of the total samples, we found a significant difference in the AGN clustering amplitude between the two fields, the best fit correlation parameters being r0=10.3 +- 1.7 h^-1 Mpc, gamma=1.33 +- 0.14 in the CDFS, and r0=5.5 +- 0.6 h^-1 Mpc, gamma=1.50 +- 0.12 in the CDFN. Within each field no statistically significant difference is found between soft and hard X-ray selected sources or between type 1 and type 2 AGN. (abridged)

The detection of excess of soft X-ray or Extreme Ultraviolet (EUV) radiation, above the thermal contribution from the hot intracluster medium (ICM), has been a controversial subject ever since the initial discovery of this phenomenon. We use a large--scale hydrodynamical simulation of a concordance $\Lambda$CDM model, to investigate the possible thermal origin for such an excess in a set of 20 simulated clusters having temperatures in the range 1--7 keV. Simulated clusters are analysed by mimicking the observational procedure applied to ROSAT--PSPC data, which for the first time showed evidences for the soft X-ray excess. For cluster--centric distances $0.4< R/R_{\rm vir}< 0.7$ we detect a significant excess in most of the simulated clusters, whose relative amount changes from cluster to cluster and, for the same cluster, by changing the projection direction. In about 30 per cent of the cases, the soft X-ray flux is measured to be at least 50 per cent larger than predicted by the one--temperature plasma model. We find that this excess is generated in most cases within the cluster virialized regions. It is mainly contributed by low--entropy and high--density gas associated with merging sub--halos, rather than to diffuse warm gas. Only in a few cases the excess arises from fore/background groups observed in projection, while no evidence is found for a significant contribution from gas lying within large--scale filaments. We compute the distribution of the relative soft excess, as a function of the cluster--centric distance, and compare it with the observational result by Bonamente et al. (2003) for the Coma cluster. Similar to observations, we find that the relative excess increases with the distance from the cluster center, with no significant excess detected for $R<0.4R_{\rm vir}$. (abridged)

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. In a recent paper Mazzotta et al. show that, if the cluster is thermally complex, T_ew fails at reproducing T_spec, and propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T_sl relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and \beta--model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X-ray temperature function could bias low \sigma_8 by 10-20 per cent. This may alleviate the tension between the value of \sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

(ABRIDGED) We analyze the X-ray properties of a sample of local and high redshift galaxy clusters extracted from a large cosmological hydrodynamical simulation. This simulation has been realized using the Tree+SPH code GADGET-2 for a LambdaCDM model. In our analysis, we consider only objects with T_ew >2 keV and adopt an approach that mimics observations, associating with each measurement an error comparable with recent observations and providing best-fit results via robust techniques. Within the clusters, baryons are distributed among (i) a cold neutral phase, with a relative contribution that increases from less than 1 to 3 per cent at higher redshift, (ii) stars which contribute with about 20 per cent and (iii) the X-ray emitting plasma that contributes by 80 (76) per cent at z=0 (1) to the total baryonic budget. A depletion of the cosmic baryon fraction of ~7 (at z=0) and 5 (at z=1) per cent is measured at the virial radius, R_vir, in good agreement with adiabatic hydrodynamical simulations. We confirm that, also at redshift >0.5, power-law relations hold between gas temperature, T, bolometric luminosity, L, central entropy, S, gas mass, M_gas, and total gravitating mass, M_tot and that these relations are steeper than predicted by simple gravitational collapse. A significant, negative evolution in the L-T and L-M_tot relations and positive evolution in the S-T relation are detected at 0.5 < z < 1 in this set of simulated galaxy clusters. This is partially consistent with recent analyses of the observed properties of z>0.5 X-ray galaxy clusters. By fixing the slope to the values predicted by simple gravitational collapse, we measure at high redshift normalizations lower by 10-40 per cent in the L-T, M_tot-T, M_gas-T, f_gas-T and L-M_tot relations than the observed estimates.

We study the properties of the diffuse light in galaxy clusters forming in a large hydrodynamical cosmological simulation of the Lambda-CDM cosmology. The simulation includes a model for radiative cooling, star formation in dense cold gas, and feedback by SN-II explosions. We select clusters having mass M>10^(14) h^(-1) Msun and study the spatial distribution of their star particles. While most stellar light is concentrated in gravitationally bound galaxies orbiting in the cluster potential, we find evidence for a substantial diffuse component, which may account for the extended halos of light observed around central cD galaxies. We find that more massive simulated clusters have a larger fraction of stars in the diffuse light than the less massive ones. The intracluster light is more centrally concentrated than the galaxy light, and the stars in the diffuse component are on average older than the stars in cluster galaxies, supporting the view that the diffuse light is not a random sampling of the stellar population in the cluster galaxies. We thus expect that at least ~10% of the stars in a cluster may be distributed as intracluster light, largely hidden thus far due to its very low surface brightness.

We present results from Tree+SPH simulations of a galaxy cluster, aimed at studying the metal enrichment of the intra--cluster medium (ICM). The simulation code includes a fairly advanced treatment of star formation, as well as the release of energy feedback and detailed yields from both type-II and type-Ia supernovae, also accurately accounting for the lifetimes of different stellar populations. We perform simulations of a cluster with virial mass ~ 3.9x 10^14 Msun, to investigate the effect of varying the feedback strength and the stellar initial mass function (IMF). Although most of the models are able to produce acceptable amounts of Fe mass, we find that the profiles of the iron abundance are always steeper than observed. The [O/Fe] ratio is found to be sub--solar for a Salpeter IMF, with [O/Fe] -0.2 at R >~ 0.1R200, whereas increasing to super-solar values in central regions, as a result of recent star formation. Using a top--heavier IMF gives a larger [O/Fe] over the whole cluster, at variance with observations. On the other hand, the adoption of a variable IMF, which becomes top-heavier at z>2, provides a roughly solar [O/Fe] ratio. Our results indicate that our simulations still lack a feedback mechanism which should quench star formation at low redshift and transport metals away from the star forming regions.

We study the formation of clusters of galaxies using high-resolution hydrodynamic cosmological simulations that include the effect of thermal conduction with an effective isotropic conductivity of 1/3 the classical Spitzer value. We find that, both for a hot ($T_{\rm ew}\simeq 12$ keV) and several cold ($T_{\rm ew}\simeq 2$ keV) galaxy clusters, the baryonic fraction converted into stars does not change significantly when thermal conduction is included. However, the temperature profiles are modified, particularly in our simulated hot system, where an extended isothermal core is readily formed. As a consequence of heat flowing from the inner regions of the cluster both to its outer parts and into its innermost resolved regions, the entropy profile is altered as well. This effect is almost negligible for the cold cluster, as expected based on the strong temperature dependence of the conductivity. Our results demonstrate that while thermal conduction can have a significant influence on the properties of the intra--cluster medium of rich galaxy clusters, it appears unlikely to provide by itself a solution for the overcooling problem in clusters, or to explain the current discrepancies between the observed and simulated properties of the intra--cluster medium.

(Abriged) We present a study of the integrated physical properties of a sample of 28 X-ray galaxy clusters observed with Chandra at redshift between 0.4 and 1.3. In particular, we have twelve objects in the redshift range 0.4-0.6, five between 0.6 and 0.8, seven between 0.8 and 1 and four at z>1, compounding the largest sample available for such a study. We focus particularly on the properties and evolution of the X-ray scaling laws. We fit both a single and a double beta-model with the former which provides a good representation of the observed surface brightness profiles, indicating that these clusters do not show any significant excess in their central brightness. By using the best-fit parameters of the beta-model together with the measured emission-weighted temperature (in the range 3-11 keV), we recover gas luminosity, gas mass and total gravitating mass out to R500. We observe scaling relations steeper than expected from self-similar model by a significant (>3 sigma) amount in the L-T and Mgas-T relations. We detect hints of negative evolution in the L-T, Mgas-T and L-Mtot relations, thus suggesting that systems at higher redshift have lower X-ray luminosity and gas mass for fixed temperature. In particular, when the 16 clusters at z>0.6 are considered, the evolution becomes more evident and its power-law modelization is a statistically good description of the data.

We present results on the X-ray properties of clusters and groups of galaxies, extracted from a large hydrodynamical simulation. We used the GADGET code to simulate a LambdaCDM model within a box of 192 Mpc/h on a side, with 480^3 dark matter particles and as many gas particles. The simulation includes radiative cooling, star formation and supernova feedback. The simulated M-T relation is consistent with observations once we mimic the procedure for mass estimates applied to real clusters. Also, with the adopted choices of Omega_m=0.3 and sigma_8=0.8 the resulting XTF agrees with observational determinations. The L-T relation also agrees with observations for clusters with T>2 keV, with no change of slope at the scale of groups. The entropy in central cluster regions is higher than predicted by gravitational heating alone, the excess being almost the same for clusters and groups. The simulated clusters appear to have suffered some overcooling, with f*~0.2, thus about twice as large as the value observed. Interestingly, temperature profiles are found to steadily increase toward cluster centers. They decrease in the outer regions, much like observational data do at r>0.2r_vir, while not showing an isothermal regime followed by a smooth temperature decline in the innermost regions.

We present deep Chandra and XMM obervations of the galaxy cluster RDCS1252.9-2927, which was selected from the ROSAT Deep Cluster Survey (RDCS) and confirmed by extensive spectroscopy with the VLT at redshift z=1.237. With the Chandra data, the X-ray emission from the intra-cluster medium is well resolved and traced out to 500 kpc, thus allowing a measurement of the physical properties of the gas with unprecedented accuracy at this redshift. We detect a clear 6.7 keV Iron K line in the Chandra spectrum providing a redshift within 1% of the spectroscopic one. By augmenting our spectroscopic analysis with the XMM data (MOS detectors only), we significantly narrow down the 1 sigma error bar to 10% for the temperature and 30% for the metallicity, with best fit values kT = 6.0(+0.7,-0.5) keV, Z = 0.36(+0.12,-0.10) Z_sun. In the likely hypothesis of hydrostatic equilibrium, we measure a total mass of M_{500} = (1.9+-0.3)10^14 h_70^{-1}M_sun within R_{Delta=500} = 536 kpc. Overall, these observations imply that RDCS1252.9-2927 is the most X-ray luminous and likely the most massive bona-fide cluster discovered to date at z>1. When combined with current samples of distant clusters, these data lend further support to a mild evolution of the cluster scaling relations, as well the metallicity of the intra-cluster gas. Inspection of the cluster mass function in the current cosmological concordance model (h,Omega_m,Omega_Lambda)=(0.7,0.3,0.7) and sigma_8=0.7-0.8 shows that RDCS1252.9-2927 is an M* cluster at z=1.24, in keeping with number density expectations in the RDCS survey volume.

We present the analysis of the X-ray spectra of 18 distant clusters of galaxies with redshift 0.3 < z < 1.3. Most of them were observed with the Chandra satellite in long exposures ranging from 36 ks to 180 ks. For two of the z>1 clusters we also use deep XMM-Newton observations. Overall, these clusters probe the temperature range 3<kT<8 keV. Our analysis is aimed at deriving the iron abundance in the Intra Cluster Medium (ICM) out to the highest redshifts probed to date. Using a combined spectral fit of cluster subsamples in different redshift bins, we investigate the evolution of the mean ICM metallicity with cosmic epoch. We find that the mean Fe abundance at <z> =0.8 is Z =0.25 +0.04 -0.06 Zo, consistent with the local canonical metallicity value, Z ~0.3 Zo, within 1sigma confidence level. Medium and low temperature clusters (kT <5 keV) tend to have larger iron abundances than hot clusters. At redshift <z>~1.2 (4 clusters at z>1) we obtain a statistically significant detection of the Fe-K line only in one cluster (Z>0.10 Zo at the 90% c.l.). Combining all the current data set from Chandra and XMM at z>1, the average metallicity is measured to be <Z> = 0.21 +0.10 -0.05 Zo (1sigma error), thus suggesting no evolution of the mean iron abundance out to z~1.2.

We report the discovery of large scale structures of X-ray sources in the 1Msec observation of the Chandra Deep Field South. Two main structures appear as narrow (Delta_z < 0.02) spikes in the source redshift distribution at z=0.67 and z=0.73, respectively. Their angular distribution spans a region at least ~ 17 arcmin wide, corresponding to a physical size of 7.3 h_{70}^{-1} Mpc at a redshift of z ~ 0.7 (Omega_m=0.3, Omega_{Lambda}=0.7). These spikes are populated by 19 sources each, which are mainly identified as Active Galactic Nuclei (AGN). Two sources in each spike are extended in X-rays, corresponding to galaxy groups/clusters embedded in larger structures. The X-ray source redshift distribution shows other spikes, the most remarkable at z=1.04, 1.62 and 2.57. This is one of the first evidences for large scale structure traced by X-ray sources and for spatial clustering of X-ray selected AGN. The X-ray data have been complemented with the spectroscopic data from the K20 near infrared survey (Cimatti et al. 2002), which covers ~1/10 of the X-ray field. Also in this survey the source redshift distribution shows several spikes. Two narrow structures at z=0.67 and z=0.73 (again with Delta_z ~ 0.02) are the most significant features, containing 24 and 47 galaxies, respectively. While the K20 structure at z=0.73 is dominated by a standard galaxy cluster with a significant concentration around a central cD galaxy and morphological segregation, the galaxies at z=0.67 constitute a loose structure rather uniformly distributed along the K20 field.[abridged]

We discuss Tree+SPH simulations of galaxy clusters and groups, aimed at studying the effect of cooling and non-gravitational heating on observable properties of the ICM. We simulate at high resolution four halos,with masses in the range (0.2-4)10^{14}M_sol. We discuss the effects of using different SPH implementations and show that high resolution is mandatory to correctly follow the cooling pattern of the ICM. All of our heating schemes which correctly reproduce the X-ray scaling properties of clusters and groups do not succeed in reducing the fraction of collapsed gas below a level of 20 (30) per cent at the cluster (group) scale. Finally, gas compression in cooling cluster regions causes an increase of the temperature and a steepening of the temperature profiles, independent of the presence of non-gravitational heating processes. This is inconsistent with recent observational evidence for a decrease of gas temperature towards the center of relaxed clusters. Provided these discrepancies persist even for a more refined modeling of energy feedback, they may indicate that some basic physical process is still missing in hydrodynamical simulations.

Different determinations currently suggest scattered values for the power spectrum normalization on the scale of galaxy clusters, sigma_8. Here we concentrate on the constraints coming from the X-ray temperature and luminosity functions (XTF and XLF), and investigate several possible sources of discrepancies in the results. We conclude that the main source of error in both methods is the scaling relation involved, in particular the way its intrinsic scatter and systematic normalization are treated. For temperature derived constraints, we use a sample adapted from HIFLUGCS, and test for several sources of systematic error. We parameterize the mass-temperature relation with an overall factor T_ast, with 1.5 \le T_ast \le 1.9. After marginalising over the range of T_ast, we obtain a 68 per cent confidence range of sigma_8=0.77^{+0.05}_{-0.04} for a standard LambdaCDM model. For luminosity derived constraints we use the XLF from the REFLEX survey and explore how sensitive the final results are on the details of the mass-luminosity, M-L, conversion. Assuming a uniform systematic uncertainty of +/-20 per cent in the amplitude of the mass-luminosity relation by Reiprich and Bohringer, we derive sigma_8=0.79^{+0.06}_{-0.07} for the same standard LambdaCDM model. Although the XTF and XLF derived constraints agree very well with each other, we emphasize that such results can change by about 10--15 per cent, depending on how uncertainties in the L-T-M conversions are interpreted and included in the analysis. In order to achieve precision cosmology, with this probe, it is important to separate the uncertainty in the scaling relation into its intrinsic and overall normalization parts.

We apply the ZTRACE algorithm to the optical NOG and infra-red PSCz galaxy catalogues to reconstruct the pattern of primordial fluctuations that have generated our local Universe. We check that the density fields traced by the two catalogues are well correlated, and consistent with a linear relation (either in $\delta$ or in log(1+$\delta$)) with relative bias (of NOG with respect to PSCz) b_rel = 1.1 $\pm$ 0.1. The relative bias relation is used to fill the optical zone of avoidance at |b| < 20$^\circ$ using the PSCz galaxy density field. We perform extensive testing on simulated galaxy catalogues to optimize the reconstruction. The quality of the reconstruction is predicted to be good at large scales, up to a limiting wavenumber k_lim $\simeq$ 0.4 h/Mpc beyond which all information is lost. We find that the improvement due to the denser sampling of the optical catalogue is compensated by the uncertainties connected to the larger zone of avoidance. The initial conditions reconstructed from the NOG catalogue are found (analogously to those from the PSCz) to be consistent with Gaussian paradigm. We use the reconstructions to produce sets of initial conditions ready to be used for constrained simulations of our local Universe.

We carry out a comparison between observations and hydrodynamic simulations of entropy profiles of groups and clusters of galaxies. We use the Tree+SPH GADGET code to simulate four halos of sizes in the M_500 = 1.0 - 16.e13 h^-1 Msun range, corresponding to poor groups up to Virgo-like clusters. We concentrate on the effect of introducing radiative cooling, star formation, and a variety of non-gravitational heating schemes on the entropy structure and the stellar fraction. We show that all the simulations result in a correct entropy profile for the Virgo-like cluster. With the heating energy budget of ~0.7 keV/particle injected at z_h=3, we are also able to reproduce the entropy profiles of groups. We obtain the flat entropy cores as a combined effect of preheating and cooling, while we achieve the high entropy at outskirts by preheating. The resulting baryon fraction locked into stars is in the 25-30% range, compared to 35-40% in the case of no preheating. Heating at higher redshift, z_h=9, strongly delays the star-formation, but fails to produce a sufficiently high specific entropy.

Considerable progress has been made over the last decade in the study of the evolutionary trends of the population of galaxy clusters in the Universe. In this review we focus on observations in the X-ray band. X-ray surveys with the ROSAT satellite, supplemented by follow-up studies with ASCA and Beppo-SAX, have allowed an assessment of the evolution of the space density of clusters out to z~1, and the evolution of the physical properties of the intra-cluster medium out to z~0.5. With the advent of Chandra and Newton-XMM, and their unprecedented sensitivity and angular resolution, these studies have been extended beyond redshift unity and have revealed the complexity of the thermodynamical structure of clusters. The properties of the intra-cluster gas are significantly affected by non-gravitational processes including star formation and Active Galactic Nucleus (AGN) activity. Convincing evidence has emerged for modest evolution of both the bulk of the X-ray cluster population and their thermodynamical properties since redshift unity. Such an observational scenario is consistent with hierarchical models of structure formation in a flat low density universe with Omega_m=0.3 and sigma_8=0.7-0.8 for the normalization of the power spectrum. Basic methodologies for construction of X-ray-selected cluster samples are reviewed and implications of cluster evolution for cosmological models are discussed.

We present a set of gas-dynamical simulations of galaxy groups and clusters aimed at exploring the effect of non-gravitational heating. We use GASOLINE, a parallel Tree+SPH code, to simulate the formation of four cosmic halos with temperature 0.5<T<8 keV. Non-gravitational heating is implemented in two different ways: (1) by imposing a minimum entropy floor at a given redshift, 1<z<5; (2) by gradually heating gas, proportionally to the SN rate expected from semi-analytical modeling of galaxy formation. Our main results are the following. (a) An extra heating energy of about 1 keV per gas particle is required to reproduce the observed Lx-T relation, independent of whether it is provided so as to create an entropy floor of 50-100 keV cm^2, or is modulated in redshift; our SN feedback recipe provides only 1/3 keV/part. (b) The M-T relation is almost unaffected by non-gravitational heating and follows the M T^{3/2} scaling, with a normalization ~40% higher than observed, independent of the heating scheme. The inclusion of cooling in a run of a small group has the effects of increasing T_ew by ~30%, possibly reconciling simulated and observed M-T relations, and of decreasing Lx by ~40%. In spite of the inclusion of SN feedback energy, almost 40% of the gas becomes cold, in excess of current observational estimates. (abridged)

We compute the chemical and thermal history of the intra-cluster medium in rich and poor clusters under the assumption that supernovae (I, II) are the major responsible both for the chemical enrichment and the heating of the intra-cluster gas. We assume that only ellipticals and S0 galaxies contribute to the enrichment and heating of the intra-cluster gas through supernova driven winds and explore several prescriptions for describing the feed-back between supernovae and the interstellar medium in galaxies. We integrate then the chemical and energetical contributions from single cluster galaxies over the cluster luminosity function and derive the variations of these quantities as functions of the cosmic time. We reach the following conclusions: i) while type II supernovae dominates the chemical enrichment and energetics inside the galaxies, type Ia supernovae play a predominant role in the intra-cluster medium, ii) galaxy models, which reproduce the observed chemical abundances and abundance ratios in the intra-cluster medium, predict a maximum of 0.3-0.4 keV per particle of energy input, a result obtained by assuming that type Ia supernovae contribute 100% of their initial blast wave energy whereas type II supernovae contribute only by a few percents of their initial energy.

In this Paper we present the source catalog obtained from a 942 ks exposure of the Chandra Deep Field South (CDFS), using ACIS-I on the Chandra X-ray Observatory. Catalog generation proceeded simultaneously using two different methods; a method of our own design using a modified version of the SExtractor algorithm, and a wavelet transform technique developed specifically for Chandra observations. The detection threshold has been set in order to have less than 10 spurious sources, as assessed by extensive simulations. We subdivided the catalog into four sections. The primary list consists of objects common to the two detection methods. Two secondary lists contain sources which were detected by: 1) the SExtractor algorithm alone and 2) the wavelet technique alone. The fourth list consists of possible diffuse or extended sources. The flux limits at the aimpoint for the soft (0.5--2 keV) and hard (2--10 keV) bands are 5.5E-17 erg/s/cm^2 and 4.5E-16 erg/s/cm^2 respectively. The total number of sources is 346; out of them, 307 were detected in the 0.5--2 keV band, and 251 in the 2--10 keV band. We also present optical identifications for the catalogued sources. Our primary optical data is R band imaging to a depth of R~26.5 (Vega). We found that the R-band/Chandra offsets are small, ~1 arcsec. Coordinate cross-correlation finds 85% of the Chandra sources covered in R to have counterparts within the 3-sigma error box (>~1.5 arcsec depending on off-axis angle and signal-to-noise). The unidentified fraction of sources, approximately ~10--15 %, is close to the limit expected from the observed X-ray flux to R-band ratio distribution for the identified sample.

We present the main results from our 940 ksec observation of the Chandra Deep Field South (CDFS), using the source catalog described in an accompanying paper (Giacconi et al. 2001). We extend the measurement of source number counts to 5.5e-17 erg/cm^2/s in the soft 0.5-2 keV band and 4.5e-16 erg/cm^2/s in the hard 2-10 keV band. The hard band LogN-LogS shows a significant flattening (slope~=0.6) below ~1e-14 erg/cm^2/s, leaving at most 10-15% of the X-ray background (XRB) to be resolved, the main uncertainty lying in the measurement of the total flux of the XRB. On the other hand, the analysis in the very hard 5-10 keV band reveals a relatively steep LogN-LogS (slope ~=1.3) down to 1e-15 erg/cm^2/s. Together with the evidence of a progressive flattening of the average X-ray spectrum near the flux limit, this indicates that there is still a non negligible population of faint hard sources to be discovered at energies not well probed by Chandra, which possibly contribute to the 30 keV bump in the spectrum of the XRB. We use optical redshifts and identifications, obtained with the VLT, for one quarter of the sample to characterize the combined optical and X-ray properties of the CDFS sample. Different source types are well separated in a parameter space which includes X-ray luminosity, hardness ratio and R-K color. Type II objects, while redder on average than the field population, have colors which are consistent with being hosted by a range of galaxy types. Type II AGN are mostly found at z<~1, in contrast with predictions based on AGN population synthesis models, thus suggesting a revision of their evolutionary parameters.

We discuss our current progress in studying a sample of z>0.8 clusters of galaxies from the ROSAT Distant Cluster Survey. To date, we have Chandra observations for four of the ten clusters. We find that the morphology of two of these four are quite regular, with deviations from circular of less than 5%, while two are strikingly elliptical. When the temperatures and luminosities of our sample are grouped with six other high-redshift measurements, there is no measured evolution in the luminosity-temperature relation. We identify a number of X-ray emitting point sources that are potential cluster members. These could be sources of intracluster medium heating, adding the entropy necessary to explain the cluster luminosity-temperature relation.

In this paper we show preliminary results concerning the chemical and energetic enrichment of the ICM by means of supernova-driven wind models in elliptical galaxies. These models are obtained by taking into account new prescriptions about supernova remnant evolution in the interstellar medium. We find that models, which can reproduce the Fe abundance and the [$\alpha$/Fe] ratios observed in the ICM, predict that the SN energy input can provide about 0.3 keV per ICM particle. We have obtained this result by assuming that each SN explosion inject on the average into the ISM no more than 20% of its initial blast wave energy. The predicted energy per particle is not enough to break the cluster self-similarity but is more than predicted in previous models.

We present results from high-resolution Tree+SPH simulations of galaxy clusters and groups, aimed at studying the effect of non-gravitational heating on the entropy of the ICM. We simulate three systems, having emission-weighted temperature T=0.6,1 and 3 keV, with spatial resolution better than 1% of the virial radius. We consider the effect of different prescriptions for non-gravitational ICM heating, such as SN energy feedback, as predicted by semi-analytical models of galaxy formation, and two different minimum entropy floors, S_fl=50 and 100 keV cm^2, imposed at z=3. Simulations with only gravitational heating nicely reproduce predictions from self-similar ICM models, while extra heating is shown to break the self-similarity, by a degree which depends on total injected energy and on cluster mass. We use observational results on the excess entropy in central regions of galaxy systems, to constrain the amount of extra-heating required. We find that setting the entropy floor S_fl=50 keV cm^2, which corresponds to an extra heating energy of about 1 keV per particle, is able to reproduce the observed excess of ICM entropy.

We present results from 300 ks of X-ray observations of the Chandra Deep Field South. The field of the four combined exposures is now 0.1035 deg^2 and we reach a flux limit of 10^{-16} erg s^{-1} cm^{-2} in the 0.5-2 keV soft band and 10^{-15} erg s^{-1} cm^{-2} in the 2-10 keV hard band, thus a factor 2 fainter than the previous 120 ks exposure. The total catalogue is composed of 197 sources including 22 sources detected only in the hard band, 51 only in the soft band, and 124 detected in both bands. We have now the optical spectra for 86 optical counterparts. We compute the total contribution to the X-ray background in the 2-10 keV band, which now amounts to (1.45\pm 0.15)*10^{-11} erg cm^{-2} s^{-1} deg^{-2} (after the inclusion of the ASCA sources to account for the bright end) to a flux limit of 10^{-15} erg s^{-1} cm^{-2}. This corresponds to 60-90% of the unresolved hard X-ray background (XRB), given the uncertainties on its actual value.[ABRIDGED]

This Meeting featured the recent advancements in our understanding of galaxy clusters and the distant Universe, achieved by the past and new generation of X-ray satellites. I summarize here the main themes that have been discussed: (a) Clusters of galaxies as probes of cosmological models; (b) The physics of cosmic baryons trapped within the potential wells of galaxy clusters; (c) The origin of the cosmic X-ray background and the nature of the contributing sources.

We analyze the ROSAT Deep Cluster Survey (RDCS) to derive cosmological constraints from the evolution of the cluster X-ray luminosity distribution. The sample contains 103 galaxy clusters out to z=0.85 and flux-limit Flim=3 10^{-14} cgs (RDCS-3) in the [0.5-2.0] keV energy band, with a high-z extension containing four clusters at 0.90<z<1.26 and F>1 10^{-14} cgs (RDCS-1). Model predictions for the cluster mass function are converted into the X-ray luminosity function in two steps. First we convert mass into intra-cluster gas temperature by assuming hydrostatic equilibrium. Then temperature is converted into X-ray luminosity by using the most recent data on the Lx-T relation for nearby and distant clusters. These include the Chandra data for seven distant clusters at 0.57<z<1.27. From RDCS-3 we find \Omega_m=0.35+/-0.12 and \sigma_8=0.66+/-0.06 for a spatially flat Universe with cosmological constant, with no significant constraint on \Gamma . Even accounting for theoretical and observational uncertainties in the mass/X-ray luminosity conversion, an Einstein-de-Sitter model is always excluded at far more than the 3sigma level.

We present our first results from 130 ks of X-ray observations obtained with the Advanced CCD Imaging Spectrometer on the Chandra X-ray Observatory. We reach a flux of 2 X 10^(-16) erg s^(-1) cm^(-2) in the 0.5-2 keV soft band and 2 X 10^(-15) erg s^(-1) cm^(-2) in the 2-10 keV hard band. For the optical identification we carried out a survey in VRI with the FORS-1 imaging-spectrometer on the ANTU telescope (UT-1 at VLT) complete to R <= 26. This dataset was complemented with data from the ESO Imaging Survey (EIS) in the UBJK bands and the ESO Wide Field Imager Survey (WFI) in the B band. The positional accuracy of the X-ray detections is of order of 1'' in the central 6'. Optical identifications are found for ~90% of the sources. We obtain the cumulative spectra of the faint and bright X-ray sources in the sample. A power law fit in the range 2-10 keV using the galactic value of N_H ~ 8 x 10^(19) cm^(-2), yields a photon index of Gamma = 1.70+-0.06 and 1.35+-0.10 for the bright and faint sample respectively, showing a flattening of the spectrum at lower fluxes. We discuss the LogN-LogS relationship and the discrete source contribution to the integrated X-ray sky flux. In the soft band, the sources now detected at fluxes below 10^(-15) erg s^(-1) cm^(-2) contribute (4.0 +- 0.3) X 10^(-12) erg cm^(-2)s^(-1) deg^(-2) to the total XRB. The flux resolved in the hard band down to the flux limit of 2 X 10^(-15) erg s^(-1) cm^(-2) now contributes a total of 1.6 X 10^(-11) erg cm^(-2) s^(-1) deg^(-2) which amounts to a fraction of 60-80% of the total measured background. This result confirms that the XRB is due to the integrated contribution of discrete sources, but shows that there is still a relevant fraction (at least 20%) of the hard XRB to be resolved at fluxes below 10^(-15) erg s^(-1) cm^(-2). (ABRIDGED)

Clusters of galaxies outline the network of the distribution of visible matter in the Universe, marking the highest-mass knots where filamentary structures join together. If we observe the sky in X rays, clusters of galaxies stand out as cosmic lighthouses by virtue of a thin gas trapped and heated within their gravitational potential wells. This powerful emission is directly linked to the total gravitating mass they contain, such that they can be efficiently used as tracers of the cosmic mass distribution within a sizeable fraction of the observable Universe. Recent observational campaigns have for the first time used X-ray clusters to map cosmic structures on scales approaching 3 billion light years. The emerging picture is remarkably consistent with the expectation of a low-density Universe dominated by cold dark matter.

In this contribution the ongoing effort to build a statistical sample of clusters of galaxies over a wide range of redshifts to study the evolution of clusters and member galaxies is reviewed. The starting point for this project has been the list of candidate clusters identified from the I-band EIS-WIDE survey data. Since the completion of this survey, new optical/infrared observations have become available and have been used to confirm some of these candidates using the photometric data alone or in combination with the results of follow-up spectroscopic observations. Our preliminary results show that the yield of real physical associations from the original catalog is conservatively >60% and that a large sample of clusters in the southern hemisphere, extending to high-redshifts, is within reach.

The Chandra X-ray Observatory was used to obtain a 190 ks image of three high redshift galaxy clusters in one observation. The results of our analysis of these data are reported for the two z > 1 clusters in this Lynx field, including the most distant known X-ray selected cluster. Spatially-extended X-ray emission was detected from both these clusters, indicating the presence of hot gas in their intracluster media. A fit to the X-ray spectrum of RX J0849+4452, at z=1.26, yields a temperature of kT = 5.8^{+2.8}_{-1.7} keV. Using this temperature and the assumption of an isothermal sphere, the total mass of RX J0849+4452 is found to be 4.0^{+2.4}_{-1.9} X 10^{14} h_{65}^{-1} M_{\sun} within r = 1 h_{65}^{-1} Mpc. The T_x for RX J0849+4452 approximately agrees with the expectation based on its L_{bol} = 3.3^{+0.9}_{-0.5} X 10^{44}$ erg s$^{-1} according to the low redshift L_x - T_x relation. The very different distributions of X-ray emitting gas and of the red member galaxies in the two z > 1 clusters, in contrast to the similarity of the optical/IR colors of those galaxies, suggests that the early-type galaxies mostly formed before their host clusters.

We present an in-depth weak lensing analysis of the cluster MS1008 based on deep multicolor imaging obtained during the Science Verification of FORS1 at the VLT. The image quality (half arcsec seeing) and depth of the VLT images allow the shear signal to be mapped with high signal-to-noise and to be traced out to 1.2 h_50^-1 Mpc, near the edge of the 6'.8x6'.8 field of view. Using BVRI color information, as well as 81 redshifts in the field from the CNOC survey, background galaxies can be effectively separated from cluster and foreground objects. PSF distorsions are found to be moderate across the FORS images and thus easily removed. Due to the small statistical errors in the mass reconstruction, this dataset provides a testing ground where several systematic effects (e.g. mass-sheet degeneracy, redshift distribution of the background sources, cluster galaxy contamination) can be quantified. Several methods are used to remove the mass-sheet degeneracy which is found to dominate the systematic error budget. We measure a lower limit to the mass of 2.6 10^14 h_50^-1 M_sun within 1 h_50^-1 Mpc and a ``total'' mass of 5.3 10^14 h_50^-1 M_sun by fitting a softened isothermal sphere. The availability of the CNOC redshift data and X-ray observations on this cluster allow a comparison of different determinations of the mass radial profile. We find the lensing and X-ray measurements in excellent agreement, while the mass derived from the virial analysis is marginally (1-2 sigma) in agreement at radii where both methods are reliable.

We measure redshifts for 67 galaxies in the field of six cluster candidates from the ESO Imaging Survey (EIS). The cluster candidates are selected in the EIS patches C and D among those with estimated mean redshifts between 0.5 and 0.7. The observations were made with EFOSC2 at the 3.6m ESO telescope. In the six candidate cluster fields, we identify 19 possible sets of 2 to 7 galaxies in redshift space. In order to establish which of the 19 sets are likely to correspond to real dense systems we compare our counts with those expected from a uniform distribution of galaxies with given luminosity function. In order to take into account the effect of the Large Scale Structure, we modulate the probability computed from the luminosity function with random samplings of the Canada-France Redshift Survey. We find that four out of six candidate EIS clusters are likely to correspond to real systems in redshift space (> 95 % confidence level). Two of these systems have mean redshift in agreement with the redshift estimate given by the matched filter algorithm. The other two systems have significantly lower redshifts. We discuss the implications of our results in the context of our ongoing research projects aimed at defining high-redshift optically-selected cluster samples.

We present an analysis of the ENEAR sample of peculiar velocities of elliptical galaxies, obtained with D_n-\sigma distances. We use the velocity correlation function to analyze the statistics of the field-object's velocities, while the analysis of the cluster data is based on the estimate of their rms peculiar velocity, Vrms. The statistics of the model velocity field is parameterized by the amplitude, \eta_8=\sigma_8 \Omega_m^{0.6}, and by the shape parameter, \Gamma. From the velocity correlation statistics we obtain \eta_8=0.51{-0.09}{+0.24} for \Gamma=0.25 at the 2\sigma level. Even though less constraining, a consistent result is obtained by comparing the measured Vrms of clusters to linear theory predictions. For \Gamma=0.25 we find \eta_8=0.63{-0.19}{+0.22}$ at 1\sigma. Overall, our results point toward a statistical concordance of the cosmic flows traced by spirals and early-type galaxies, with galaxy distances estimated using TF and D_n-\sigma distance indicators, respectively.

We briefly review our current knowledge of the space density of distant X-ray clusters as measured by several ROSAT serendipitous surveys. We compare old and new determinations of the cluster X-ray Luminosity Function (XLF) at increasing redshifts, addressing the controversial issue of the evolution of its high end. We use complete subsamples, drawn from the ROSAT Deep Cluster Survey (RDCS), to quantify the statistical significance of the XLF evolution out to z ~1. A consistent observational picture emerges in which the bulk of the cluster population shows no significant evolution out to z ~1, whereas the most luminous systems (L_x >~ L* [0.5-2 keV] =~ 5x10^44 erg/s) were indeed rarer, at least at z >0.5, in keeping with the original findings of the EMSS. We also report on the recent spectroscopic identification of four clusters in the RDCS lying beyond z =1, the most distant X-ray clusters known to date, which set an interesting lower limit on the space density of clusters at z >1.

We use the ROSAT Deep Cluster Survey (RDCS) to trace the evolution of the cluster abundance out to $z\simeq 0.8$ and constrain cosmological models. We resort to a phenomenological prescription to convert masses into $X$-ray fluxes and apply a maximum-likelihood approach to the RDCS redshift- and luminosity-distribution. We find that, even changing the shape and the evolution on the $L_{bol}$-$T_X$ relation within the observational uncertainties, a critical density Universe is always excluded at more than $3\sigma$ level. By assuming a non-evolving $X$-ray luminosity-temperature relation with shape $L_{bol}\propto T_X^3$, it is $\Omega_m=0.35^{+0.35}_{-0.25}$ and $\sigma_8=0.76^{+0.38}_{-0.14}$ for flat models, with uncertainties corresponding to $3\sigma$ confidence levels.

In this paper we study the isocurvature mode contribution to the cosmic microwave background anisotropies and the large scale structure power spectrum, for a two-field model of inflation proposed by Linde and Mukhanov. We provide constraints on the parameters of the model by comparing its predictions with observations of the microwave background anisotropies, large scale structure data on the galaxy power spectrum, and the number density of nearby galaxy clusters. We find that such models are consistent with observations for a narrow range of parameters. As our main result, we find that only a very small isocurvature component is allowed, $\al \le 0.006$, for any underlying Friedmann model. Furthermore, we give the expected resolution with which the model parameters will be determined from future satellite missions like MAP and Planck, for a fiducial flat $\Lambda$CDM model. We find that Planck mission will be able to detect such small contributions, especially if polarization information is included. The isocurvature spectral index $\niso$ will also be determined with better than 8% precision.

Using EFOSC2 at the 3.6m ESO telescope, we obtained redshifts for 68 galaxies in the field of six cluster candidates from the ESO Imaging Survey (EIS). The cluster candidates were selected in the EIS patches C and D and have estimated mean redshifts between z=0.5 and z=0.7. In the six candidate cluster fields, we identify possible systems of galaxies in the redshift space. The likelihoods of these systems are established by comparison with random resamplings of the Canada-France Redshift Surveys, and using the redshift distribution expected from Postman et al.'s luminosity function, with Poggianti's K- and evolutionary-corrections. Four of the six candidate EIS clusters are found to correspond to a real system in the redshift space, with > 95 % probability. Two of them have a mean redshift in agreement with the estimate from the matched filter algorithm, while the other two have a significantly smaller redshift. The independent analysis of the V-I vs. I color-magnitude diagrams for five of our six cluster fields, supports our conclusions based on the spectroscopic data.

We present results of a statistical analysis of the SFI catalog of peculiar velocities, a recently completed survey of spiral field galaxies with I-band Tully-Fisher distances (Haynes et al. 1999). The velocity field statistic utilized is the velocity correlation function, $\psi_1(r)$ (Gorski et al. 1989). The analysis is performed in redshift space, so as to circumvent potential ambiguities connected with inhomogeneous Malmquist bias corrections. The results from the SFI sample are compared with linear-theory predictions. We generate a large set of mock samples, extracted from N-body simulations, which are used to assess the reliability of our analysis and to estimate the associated uncertainties. Defining $\eta_8=\sigma_8\Omega_0^{0.6}$, we find that the measured $\psi_1(r)$ implies a degenerate constraint in the $\eta_8-\Gamma$ plane, with $\eta_8=(0.3 +/- 0.1) (\Gamma/0.2)^{0.5}$, at the $2\sigma$ level, for the inverse Tully-Fisher (ITF) calibration presented in this paper. Model constraints are quite sensitive to the ITF calibration. The other ITF calibrations by Giovanelli et al. (1997) and da Costa et al. (1998) both give, for $\Gamma=0.2$, $\eta_8\simeq 0.6$ as the best-fitting value.

We present the results of the analysis of the internal velocity dispersions, \sigma_v, for the CNOC sample of distant galaxy clusters, based on an interlopers removal algorithm, which is different from that originally applied by Carlberg et al. (1996, C96). We find that the resulting \sigma_v values are consistent within <10% with the original C96 estimates. This result points in favor of a substantial robustness of currently applied methods for optical studies of the internal cluster dynamics. The resulting distribution of velocity dispersions is used to trace the redshift evolution of the cluster abundance with the aim of constraining the matter density parameter, \Omega_m. We find that constraints on \Omega_m are very sensitive to the adopted value of \sigma_8, as obtainable from the local cluster abundance: as \sigma_8 varies from 0.5 to 0.6 (for Omega_m=1), the best fitting Omega_m varies in the range 0.3-1.0.

We analyze a sample of 105 clusters having virial mass homogeneously estimated and for which galaxy magnitudes are available with a well defined high degree of completeness. In particular, we consider a subsample of 89 clusters with B_j band galaxy magnitudes taken from the COSMOS/UKST Southern Sky Object Catalogue. We compute cluster luminosities L_{B_j} within several clustercentric distances and within the virialization radius R_{vir}. In particular, we use the luminosity function and background counts estimated by Lumsden et al. (1997) on the Edinburgh/Durham Southern Galaxy Catalogue. We analyze the effect of several uncertainties connected to photometric data, fore/background removal, and extrapolation below the completeness limit of the photometry, in order to assess the robustness of our cluster luminosity estimates. We draw our results on the relations between luminosity and dynamical quantities from the COSMOS sample by considering mass and luminosities determined within the virialization radius. We find a very good correlation between cluster luminosity, L_{B_j}, and galaxy velocity dispersion, sigma_v, with L_{B_j} proportional to sigma_v^{2.1--2.3}. Our estimate of typical value for the mass-to-light ratio is M/L_{B_j} about 250 (in solar units). We do not find any correlation of M/L_{B_j} with cluster morphologies, i.e. Rood--Sastry and Bautz--Morgan types, and only a weak significant correlation with cluster richness. We find that mass has a slight, but significant, tendency to increase faster than the luminosity does, M proportional to L_{B_j}^{1.2--1.3}. We verify the robustness of this relation against a number of possible systematics.

The ROSAT Deep Cluster Survey (RDCS) has provided a new large deep sample of X-ray selected galaxy clusters. Observables such as the flux number counts n(S), the redshift distribution n(z) and the X-ray luminosity function (XLF) over a large redshift baseline (z\lesssim 0.8) are used here in order to constrain cosmological models. Our analysis is based on the Press-Schechter approach, whose reliability is tested against N-body simulations. Following a phenomenological approach, no assumption is made a priori on the relation between cluster masses and observed X-ray luminosities. As a first step, we use the local XLF from RDCS, along with the high-luminosity extension provided by the XLF from the BCS, in order to constrain the amplitude of the power spectrum, \sigma_8, and the shape of the local luminosity-temperature relation. We obtain \sigma_8=0.58 +/- 0.06 for Omega_0=1 for open models at 90% confidence level, almost independent of the L-T shape. The density parameter \Omega_0 and the evolution of the L-T relation are constrained by the RDCS XLF at z>0 and the EMSS XLF at z=0.33, and by the RDCS n(S) and n(z) distributions. By modelling the evolution for the amplitude of the L-T relation as (1+z)^A, an \Omega_0=1 model can be accommodated for the evolution of the XLF with 1<A<3 at 90% confidence level, while \Omega_0=0.4^{+0.3}_{-0.2} and \Omega_0<0.6 are implied by a non--evolving L-T for open and flat models, respectively.

We present an analysis of the 2-point correlation function, of the X-ray Brightest Abell-type Cluster sample (XBACs; Ebeling et al. 1998) and of the cosmological constraints that it provides. If \xi(r) is modelled as a power-law, we find r_0=26.0 +/- 4.5 Mpc/h and \gamma=2.0 +/-0.4, with errors corresponding to 2\sigma uncertainties. Only a marginal increase of the correlation amplitude is found as the flux limit is increased from 5x10^{-12} cgs to 12x10^{-12} cgs, thus indicating a weak dependence of the correlation amplitude on the cluster X-ray luminosity. Furthermore, we present a method to predict correlation functions for flux-limited X-ray cluster samples from cosmological models. The method is based on the analytical recipe by Mo & White (1996) and on an empirical approach to convert cluster fluxes into masses. We use a maximum-likelihood method to place constraints on the model parameter space. We find that the shape parameter of the power spectrum lies in the 2\sigma range 0.05<\Gamma<0.20. As for the amplitude of the power-spectrum, we find \sigma_8 in the range 0.4-0.8 for \Omega_0=1 and 0.8-2.0 for \Omega_0=0.3. This result is in agreement with, although less constraining than, results based on the local cluster abundance.

We have obtained new Tully-Fisher (TF) peculiar velocity measurements for 52 Abell galaxy clusters distributed throughout the sky between ~ 50 and 200 Mpc/h.The measurements are based on I band photometry and optical rotation curves for a sample of 522 spiral galaxies, from which an accurate TF template relation has been constructed. Individual cluster TF relations are referred to the template to compute cluster peculiar motions. The reflex motion of the Local Group of galaxies is measured with respect to the reference frame defined by our cluster sample and the distant portion of the Giovanelli et al. (1998) cluster set. We find the Local Group motion in this frame to be 565+/-113 km/s in the direction (l,b)=(267,26)+/-10 when peculiar velocities are weighted according to their errors. After optimizing the dipole calculation to sample equal volumes equally, the vector is 509+/-195 km/s towards (255,33)+/-22. Both solutions agree, to within 1-sigma or better, with the Local Group motion as inferred from the cosmic microwave background (CMB) dipole. Thus, the cluster sample as a whole moves slowly in the CMB reference frame, its bulk flow being at most 200 km/s.

Using extensive N-body simulations we estimate redshift space power spectra of clusters of galaxies for different cosmological models (SCDM, TCDM, CHDM, Lambda-CDM, OCDM, BSI, tau-CDM) and compare the results with observational data for Abell-ACO clusters. Our mock samples of galaxy clusters have the same geometry and selection functions as the observational sample which contains 417 clusters of galaxies in a double cone of galactic latitude |b| > 30 degrees up to a depth of 240 Mpc/h. The power spectrum has been estimated for wave numbers k in the range 0.03 < k < 0.2 h/Mpc. For k > k_max ~ 0.05 h/Mpc the power spectrum of the Abell-ACO clusters has a power-law shape, P(k)\propto k^n, with n ~ -1.9, while it changes sharply to a positive slope at k < k_max. By comparison with the mock catalogues SCDM, TCDM (n=0.9), and also OCDM with Omega_0 = 0.35 are rejected. Better agreement with observation can be found for the Lambda-CDM model with Omega_0 = 0.35 and h = 0.7 and the CHDM model with two degenerate neutrinos and Omega_HDM = 0.2 as well as for a CDM model with broken scale invariance (BSI) and the tau-CDM model. As for the peak in the Abell-ACO cluster power spectrum, we find that it does not represent a very unusual finding within the set of mock samples extracted from our simulations.

We present a new determination of the mass function of galaxy clusters, based on optical virial mass estimates for a large sample of 152 nearby (z\le0.15) Abell-ACO clusters, as provided by Girardi et al. (1998). This sample includes both data from the literature and the new ENACS data. The resulting mass function is reliably estimated for masses larger than M_{lim}\simeq 4\times 10^{14} \msun, while it is affected by sample incompleteness at smaller masses. We find N(>M_{lim})=(6.3\pm 1.2) 10^{-6} (h^{-1} Mpc)^{-3} for cluster masses estimated within a 1.5 h^{-1} radius. Our mass function is intermediate between the two previous estimates by Bahcall & Cen (1993) and by Biviano et al. (1993). Based on the Press-Schechter approach, we use this mass function to constrain the amplitude of the fluctuation power spectrum at the cluster scale. After suitably convolving the PS predictions with observational errors on cluster masses and COBE-normalizing the fluctuation power spectrum, we find sigma_8=(0.60\pm 0.04) Omega_0^{-0.46+0.09 Omega_0} for flat low-density models and sigma_8=(0.60\pm 0.04) Omega_0^{-0.48+0.17 Omega_0} for open models (at the 90% c.l.).

The combination of measurements of the local abundance of rich clusters of galaxies and its evolution to higher redshift offers the possibility of a direct measurement of $\Omega_0$ with little contribution from other cosmological parameters. We investigate the significance of recent claims that this evolution indicates that $\Omega_0$ must be small. The most recent cluster velocity dispersion function from a compilation including the ESO Northern Abell Cluster Survey (ENACS) results in a significantly higher normalization for models, corresponding to $\sigma_8\approx 0.6$ for $\Omega_0=1$, compared to the Eke, Cole, & Frenk result of $\sigma_8=0.52\pm 0.04$. Using the ENACS data for a $z=0$ calibration results in strong evolution in the abundance of clusters, and we find that the velocity dispersion function is consistent with $\Omega_0=1$. The results are dependent upon the choice and analysis of low-redshift and high-redshift data, so at present, the data is not good enough to determine $\Omega_0$ unambiguously.

We investigate the sampling and dipole convergence properties of flux-limited samples of mock X-ray clusters in relation to their underlying ``parent'' cluster distribution. To this purpose, we resort to numerical simulations of the cluster distribution and extract samples resembling the main observational features of X-ray selected cluster samples. The flux-limited samples, being quite sparse, underestimate the amplitude of the ``parent'' cluster dipole by approximately 15 per cent on average for Local Group-like observers. However, the general shape of their dipole amplitude profiles are in relatively good agreement. We also calculate the expected contribution of clusters, selected according to the relevant criteria, to the soft (0.1-2.4 keV) extragalactic X-ray background, using the ESO Key Project X-ray luminosity function, assuming a flat universe with vanishing cosmological constant. We obtain a value of about 10 per cent of the observed XRB flux.

We analyse the formation of cosmic structures in models where the dark matter is dominated by light gravitinos with mass of $ 100$ eV -- 1 keV, as predicted by gauge-mediated supersymmetry (SUSY) breaking models. After evaluating the number of degrees of freedom at the gravitinos decoupling ($g_*$), we compute the transfer function for matter fluctuations and show that gravitinos behave like warm dark matter (WDM) with free-streaming scale comparable to the galaxy mass scale. We consider different low-density variants of the WDM model, both with and without cosmological constant, and compare the predictions on the abundances of neutral hydrogen within high-redshift damped Ly--$\alpha$ systems and on the number density of local galaxy clusters with the corresponding observational constraints. We find that none of the models satisfies both constraints at the same time, unless a rather small $\Omega_0$ value ($\mincir 0.4$) and a rather large Hubble parameter ($\magcir 0.9$) is assumed. Furthermore, in a model with warm + hot dark matter, with hot component provided by massive neutrinos, the strong suppression of fluctuation on scales of $\sim 1\hm$ precludes the formation of high-redshift objects, when the low--$z$ cluster abundance is required. We conclude that all different variants of a light gravitino DM dominated model show strong difficulties for what concerns cosmic structure formation. This gives a severe cosmological constraint on the gauge-mediated SUSY breaking scheme.

We study the geometry and topology of the large-scale structure traced by galaxy clusters in numerical simulations of a box of side 320 $h^{-1}$ Mpc, and compare them with available data on real clusters. The simulations we use are generated by the Zel'dovich approximation, using the same methods as we have used in the first three papers in this series. We consider the following models to see if there are measurable differences in the topology and geometry of the superclustering they produce: (i) the standard CDM model (SCDM); (ii) a CDM model with $\Omega_0=0.2$ (OCDM); (iii) a CDM model with a `tilted' power spectrum having $n=0.7$ (TCDM); (iv) a CDM model with a very low Hubble constant, $h=0.3$ (LOWH); (v) a model with mixed CDM and HDM (CHDM); (vi) a flat low-density CDM model with $\Omega_0=0.2$ and a non-zero cosmological $\Lambda$ term ($\Lambda$CDM). We analyse these models using a variety of statistical tests based on the analysis of: (i) the Euler-Poincar\'{e} characteristic; (ii) percolation properties; (iii) the Minimal Spanning Tree construction. Taking all these tests together we find that the best fitting model is $\Lambda$CDM and, indeed, the others do not appear to be consistent with the data. Our results demonstrate that despite their biased and extremely sparse sampling of the cosmological density field, it is possible to use clusters to probe subtle statistical diagnostics of models which go far beyond the low-order correlation functions usually applied to study superclustering.

Recently, peculiar velocity measurements became available for a new sample of galaxy clusters. From an accurately calibrated Tully-Fisher relation for spiral galaxies, we compute the rms cluster peculiar velocity and compare it to the linear theory predictions of COBE-normalized low-density and open CDM models (LambdaCDM and OCDM, respectively). Confidence levels for model rejection are estimated using a Monte Carlo procedure to generate for each model a large ensemble of artificial data sets. Following Zaroubi et al. (1997), we express our results in terms of constraints on the (Omega_0,n_pr,h) parameter space. Such constraints turn into sigma_8 Omega_0^{0.6}=0.50^{+0.25}_{-0.14} at the 90% c.l., thus in agreement with results from cluster abundance. We show that our constraints are also consistent with those implied by the shape of the galaxy power spectrum within a rather wide range for the values of the model parameters. Finally, we point out that our findings disagree at about the 3sigma level with respect to those by Zaroubi et al. (1997), based on the Mark III catalogue, which tend to prefer larger Omega_0 values within the CDM class of models.

We compare the distribution of internal velocity dispersions of galaxy clusters for an observational sample to those obtained from a set of N-body simulations of seven COBE-normalised cosmological scenarios: the standard CDM (SCDM) and a tilted (n=0.85) CDM (TCDM) model, a CHDM model with 25% of massive neutrinos, two low-density LCDM models with Omega_0=0.3 and 0.5, two open OCDM models with Omega_0=0.4 and 0.6. Simulated clusters are observed in projection so as to reproduce the main observational biases and are analysed by applying the same algorithm for interlopers removal and velocity dispersion estimate as for the reference observational sample. Velocity dispersions for individual clusters can be largely affected by observational biases in a model-dependent way: models in which clusters had less time to virialize show larger discrepancies between 3D and projected velocity dispersions. From the comparison with real clusters we find that both SCDM and TCDM largely overproduce clusters. The CHDM model marginally overproduces clusters and requires a somewhat larger sigma_8 than a purely CDM model in order to produce the same cluster abundance. The LCDM model with Omega_0=0.3 agrees with data, while the open model with Omega_0=0.4 and 0.6 underproduces and marginally overproduces clusters, respectively.

We review the prospects for dark matter (DM) in two classes of supersymmetric (SUSY) models which are characterized by the different mechanism of SUSY breaking, namely the more common supergravity models with very large scale of SUSY breaking and the recent schemes where SUSY is broken at a relatively low scale and gravitinos are likely to be the lightest SUSY particle. We point out that the former scheme is in general associated with the cold dark matter (CDM) scenario, while the latter predicts a warm dark matter (WDM) dominated Universe.

(Shortened Abstract). We analyze a redshift sample of Abell/ACO clusters and compare them with numerical simulations based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, on which we estimate cosmic variance effects. We analyse correlation statistics, the probability density function, and supercluster properties from percolation analysis. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 < \Gamma < 0.25 and an effective spectral index at 20Mpc/h in the range [-1.1,-0.9] are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press--Schechter approach, whose reliability is explicitly tested against N--body simulations. We conclude that, of the cosmological models considered here, the only viable models are either Cold+Hot DM ones with \Omega_\nu = [0.2-0.3], better if shared between two massive neutrinos, and flat low-density CDM models with \Omega_0 = [0.3-0.5].

We apply several statistical estimators to high-resolution N-body simulations of two currently viable cosmological models: a mixed dark matter model, having $\Omega_\nu=0.2$ contributed by two massive neutrinos (C+2\nuDM), and a Cold Dark Matter model with Cosmological Constant (\LambdaCDM) with $\Omega_0=0.3$ and h=0.7. Our aim is to compare simulated galaxy samples with the Perseus-Pisces redshift survey (PPS). We consider the n-point correlation functions (n=2-4), the N-count probability functions P_N, including the void probability function P_0, and the underdensity probability function U_\epsilon (where \epsilon fixes the underdensity threshold in percentage of the average). We find that P_0 (for which PPS and CfA2 data agree) and P_1 distinguish efficiently between the models, while U_\epsilon is only marginally discriminatory. On the contrary, the reduced skewness and kurtosis are, respectively, S_3\simeq 2.2 and S_4\simeq 6-7 in all cases, quite independent of the scale, in agreement with hierarchical scaling predictions and estimates based on redshift surveys. Among our results, we emphasize the remarkable agreement between PPS data and C+2\nuDM in all the tests performed. In contrast, the above \LambdaCDM model has serious difficulties in reproducing observational data if galaxies and matter overdensities are related in a simple way.

We determine the Minkowski functionals for a sample of Abell/ACO clusters, 401 with measured and 16 with estimated redshifts. The four Minkowski functionals (including the void probability function and the mean genus) deliver a global description of the spatial distribution of clusters on scales from $10$ to $60\hMpc$ with a clear geometric interpretation. Comparisons with mock catalogues of N--body simulations using different variants of the CDM model demonstrate the discriminative power of the description. The standard CDM model and the model with tilted perturbation spectrum cannot generate the Minkowski functionals of the cluster data, while a model with a cosmological constant and a model with breaking of the scale invariance of perturbations (BSI) yield compatible results.

In theories with a gauge-mediated mechanism of supersymmetry breaking the gravitino is likely to be the lightest superparticle and, hence, a candidate for dark matter. We show that the decay of the next-to-lightest superparticle into a gravitino can yield a non-thermal population of gravitinos which behave as a hot dark matter component. Together with the warm component, which is provided by the population of gravitinos of thermal origin, they can give rise to viable schemes of mixed dark matter. This realization has some specific and testable features both in particle physics and astrophysics. We outline under which conditions the mechanism remains viable even when R parity is broken.

We investigate the properties of large--scale structure predicted in a class of mixed dark matter models in which the volatile component (made of particles with high rms velocity) derives from the decay of a heavier particle. Such models based on cold+volatile dark matter (CVDM) differ from the standard mixture of CDM and massive neutrinos, usually known as CHDM, in that they involve a component which has a non--thermal phase space distribution function. As a consequence, and differently from CHDM models, the value of the redshift at which volatile particles become non relativistic, $z_{nr}$, can be varied almost independently of the volatile fraction, $\Omega_X$. We compute transfer functions for a selection of such models, having $0.1\le \Omega_X\le 0.5$ and different values of $z_{nr}$. Using linear theory and assuming a scale--free primordial spectrum, we compare such models with observational constraints on large--scale galaxy clustering and bulk flows, as well as on the abundance of galaxy clusters and high--redshift damped Ly$\alpha$ systems. We find that these constraints enable us to discriminate between different $\Omega_X$ and $z_{nr}$; within the range of the models inspected, those which can be most easily accommodated by the data correspond to the parameter choice $\Omega_X\simeq 0.2$ and $z_{nr} \simeq 2\times 10^4\Omega_X$.

Recent data on the high--redshift abundance of damped Ly$\alpha$ systems are compared with theoretical predictions for `blue' (i.e. $n>1$) Mixed Dark Matter models. The results show that decreasing the hot component fraction $\Omega_\nu$ and/or increasing the primordial spectral index $n$ implies an earlier epoch of cosmic structure formation. However, we also show that varying $\Omega_\nu$ and $n$ in these directions makes the models barely consistent with the observed abundance of galaxy clusters. Therefore, requiring at the same time observational constraints on damped Ly$\alpha$ systems and cluster abundance to be satisfied represents a challenge for the Mixed Dark Matter class of models.

We study the large-scale velocity fields traced by galaxy clusters in numerical simulations of a box of side 960 $h^{-1}$ Mpc, and compare them with available data on real clusters. In order to test the reliability of the simulations, which are based on an optimized version of the Zel'dovich approximation, we compare their cluster velocities with those of "exact" N-body simulations, and find a remarkable agreement between the two according to a variety of statistical tests. We analyse Cold Dark Matter models with density parameter in the range $0.2\le \Omega_0\le 1$, both with and without the cosmological constant term to provide a flat geometry. We also simulate a Cold+Hot Dark Matter model, with 30% provided by the hot component. Comparison with real data is performed by applying tests based on the cumulative velocity frequency distribution (CVFD) and bulk flow statistics. For the CVFD, we use observational velocity data from different authors. By merging all available data in a combined sample, the CVFD test is able to exclude only open models with $\Omega_0\le 0.4$ and flat models with $\Omega_0=0.2$. However, the analysis of individual observational samples gives contradictory results; in particular the recent Giovanelli data, which imply much lower velocities (with significantly smaller uncertainties) than the other data, are actually only consistent with the open $\Omega_0=0.4$ and the flat $\Omega_0=0.2$ models, i.e. the models which are excluded by the other data. Furthermore, CVFD and bulk flow analyses of the Branchini et al. reconstructed velocity data disfavour precisely those models accepted on the grounds of Giovanelli's sample. Finally, we confirm that the Lauer & Postman bulk flow determination is an extremely rare event in the cosmological models we have analysed.

We use numerical simulations to investigate the behaviour of the dipole moment of the spatial distribution of different kinds of mass tracers. We select density peaks of the simulated matter distribution with mean separations of 38 and 30 $h^{-1}$ Mpc to represent two samples of rich clusters of galaxies. We extract, from the same simulations, samples selected to mimic the full 3D galaxy distribution of IR galaxies, and the flux--limited IRAS and QDOT galaxy samples. We compare the dipole moments of these ``galaxy'' and ``cluster'' samples in order to assess the effects of sampling uncertainties and shot--noise on the relationship between the ``true'' underlying galaxy dipole and the dipoles obtained for clusters and for the flux--limited galaxy samples. The results of this analysis demonstrate that the dipoles of both the IRAS and QDOT--like catalogues should trace the full 3D dipole shape fairly accurately, with the loss however of about 15--20$\%$ of the total 3D dipole amplitude. Furthermore, using a simple argument based on linear perturbation theory, on the linear biasing assumption and on the amplitude of the cluster dipole relative to that of galaxies, we can estimate their relative biasing factors quite accurately and in agreement with results obtained by other methods.

These Lecture Notes are devoted to an introductory description of some of the most widely applied statistical methods for the analysis of the Large-Scale Structure (LSS) of the Universe. Rather than providing technical details about the implementation of such methods, I concentrate more on their significance and on the statistical information they provide in the comparison between data and LSS formation models. I introduce the concept of correlation functions and their connection with the count-in-cell statistics of a point distribution. Also some geometrical descriptions of the LSS, like the void probability function and the topological genus characteristics, are briefly reviewed.

We study the dipole for a redshift sample of Abell/ACO clusters. To elucidate the constraints it places on dark matter models, we use numerical simulations based on the Zel'dovich approximation. We run 20 realizations of each of six different dark matter models: four of these have a density parameter $\Omega_0=1$, while the other two have $\Omega_0=0.2$, one with and one without a cosmological constant term. For the Abell/ACO sample, we evaluate the parameter $\beta=\Omega_0^{0.6}/b_{ cl}$, which reaches its asymptotic value at $R=160 Mpc/h$. We find $\beta=0.15\pm 0.04$ and $0.25\pm 0.06$ in LG and CMB frames respectively. After identifying in the simulations LG-like observers, we construct mock cluster samples. We find that an alignment between the cluster dipole and observer velocity, such as that observed (about 20 degrees) should not be expected necessarily in the simulations. This, together with the large observer-to-observer variance estimates of $\beta$, makes it difficult to place any firm constraints on cosmological models. Our results demonstrate that the large amplitude and convergence depth of the cluster dipole cannot be taken as strong evidence either for or against a low--density Universe.

We present simulations of the cluster distribution in several dark matter models, using an optimized version of the truncated Zel'dovich approximation (TZA). We compare them with N-body cluster simulations and find that the TZA provides a very accurate description of the cluster distribution as long as fluctuations on the cluster mass scale are in the mildly non-linear regime. The simulated dark matter models are: Standard CDM (SCDM), Tilted CDM (TCDM) with n=0.7, Cold+Hot DM (CHDM) with 30% of hot component, low Hubble constant (h=0.3) CDM (LOWH) and a spatially flat low-density CDM model with \Omega_0=0.2. We compare the simulations with a redshift sample of Abell/ACO clusters, using the integral of the 2-point correlation function and the probability density function. We find that the best models at reproducing the data are CHDM and LCDM. All the other models are ruled out. The reduced skewness S_3 is fairly constant with S_3=1.9, independent of the DM model and consistent with observational data. The abundances of clusters predicted using the Press--Schechter theory provide strong constraints: only the CHDM, LOWH and LCDM models appear to produce the correct number-density of clusters.

We use the moments of counts of neighbors as given by the Generalized Correlation Integrals, to study the clustering properties of Dark Matter Halos (DH) in Cold Dark Matter (CDM) and Cold+Hot Dark Matter (CHDM) models. We compare the results with those found in the CfA and SSRS galaxy catalogs. We show that if we apply the analysis in redshift space, both models reproduce equally well the observed clustering of galaxies. Mass segregation is also found in the models: more massive DHs are more clustered compared with less massive ones. In redshift space, this mass segregation is reduced by a factor 2-3 due to the peculiar velocities. Observational catalogs give an indication of luminosity and size segregation, which is consistent with the predictions of the models. Because the mass segregation is smaller in redshift space, it is suggestive that the real luminosity or size segregation of galaxies could be significantly larger than what it is found in redshift catalogs.

We estimate the variance and the skewness of the cluster distribution in several dark matter (DM) models. The cluster simulations are based on the Zel'dovich approximation, the low computational cost of which allows us to run 50 random realizations of each model. We compare our results with those coming from a similar analysis of a redshift sample of Abell/ACO clusters. Within the list of the considered models, we find that only the model based on Cold+Hot DM (with $\Omega_{\rm hot}=0.3$) provides a good fit to the data. The standard CDM model and the low-density ($\Omega_{\circ}=0.2$) CDM models, both with and without a cosmological constant term ($\Omega_\Lambda =0.8$) are ruled out. The tilted CDM model with primordial spectral index $n=0.7$ and a low Hubble constant ($h=0.3$) CDM model are only marginally consistent with the data.

We present extended simulations of the distribution of galaxy clusters in different dark matter models, using an optimized version of the Zel'dovich approximation. Six different dark matter models are studied: Standard CDM, Open CDM, Tilted CDM, Cold + Hot DM, low-density CDM with cosmological constant, and low Hubble constant ($h=0.3$) CDM. We compare simulations with an Abell/ACO redshift sample. We find that the models that best reproduce the clustering of the real data sets are the Cold + Hot DM model and the CDM model with a cosmological constant. The probability density function of all models and data is always well approximated by a lognormal distribution. We also find that the linear biasing parameter for the simulated clusters is nearly constant over a large range of scales, but its value depends on the model. We also note that the abundances of clusters predicted by the Press-Schechter theory provide strong constraints on these models: only the Cold + Hot DM model and the low-H CDM model appear to produce the correct number-density of clusters. Taking this constraint together with the cluster clustering statistics leads one to conclude that the best of our models is the Cold + Hot DM scenario.

Damped Ly$\alpha$ systems provide possibly the most significant evidence for early structure formation, and thus a stringent constraint on the Cold + Hot Dark Matter (CHDM) cosmology. Using the numbers of halos in N-body simulations to normalize Press-Schechter (PS) estimates of the number densities of protogalaxies as a function of redshift, we find that CHDM with $\Omega_c/\Omega_\nu/\Omega_b = 0.6/0.3/0.1$ is compatible with the damped Ly$\alpha$ data only at $z<2.5$, but that it is probably incompatible with data at $z>3$. The predictions of CHDM are quite sensitive to the neutrino fraction. We find that $\Omega_c/\Omega_\nu/\Omega_b = 0.725/0.20/0.075$ is compatible with the $z>3$ data. With one massive neutrino species, this corresponds to lowering the neutrino mass from 7.0 to 4.7 eV, for $H_0=50\kmsMpc$ and $T=2.726$ K. By analysing our numerical simulations with different resolutions and box sizes as well as those of Ma \& Bertchinger (1994), we show that for the CHDM models with $\Omega_\nu$=0.2--0.3 the PS approximation should be used with Gaussian filter with $\delta_c=1.3-1.4$ if one tries to recover the total mass of a collapsed halo and to include nonlinear effects, due to waves both longer and shorter than those within the simulation box.

We use the void probability statistics to study the redshift-space galaxy distribution as described by a volume-limited subsample of the Perseus-Pisces survey. We compare the results with the same analysis realized on artificial samples, extracted from high-resolution N-body simulations by reproducing the observational biases of the real data set. Simulations are run for the Cold+HotDM model (CHDM) and for unbiased and biased (b=1.5) CDM models in a 50 Mpc/h box. We identify galaxies as residing in peaks of the evolved density field. We fragment overmerged structures into individual galaxies so as to reproduce both the correct luminosity function (after assuming M/ L values for the resulting galaxy groups) and the two-point correlation function. Our main result is that a void-probability function (VPF) from the standard CHDM model with fractions 60% cold, 30% hot, 10% barions, exceeds the observational VPF with a high confidence level. CDM models produce smaller VPF independent of the biasing parameter. We verify the robustness of this result against changing the observer position in the simulations and the galaxy identification in the evolved density field.

Although the Cold + Hot Dark Matter (CHDM) cosmology provides perhaps the best fit of any model to all the available data at the current epoch, CHDM produces structure at relatively low redshifts and thus could be ruled out if there were evidence for formation of massive objects at high redshifts. Damped Ly$\alpha$ systems are abundant in quasar absorption spectra and thus provide possibly the most significant evidence for early structure formation, and thus perhaps the most stringent constraint on CHDM. Using the numbers of halos in N-body simulations to normalize Press-Schechter estimates of the number densities of protogalaxies as a function of redshift, we find that CHDM with $\Omega_c/\Omega_\nu/\Omega_b = 0.6/0.3/0.1$ is compatible with the damped Ly$\alpha$ data at $\le 2.5$, but that it is probably incompatible with the limited $z>3$ damped Ly$\alpha$ data. The situation is uncertain because there is very little data for $z>3$, and also it is unclear whether all damped Ly$\alpha$ systems are associated with collapsed protogalaxies. The predictions of CHDM are quite sensitive to the hot (neutrino) fraction, and we find that $\Omega_c/\Omega_\nu/\Omega_b = 0.675/0.25/0.075$ is compatible even with the $z>3$ data. This corresponds to lowering the neutrino mass from 6.8 to 5.7 eV, for $H_0=50\kmsMpc$. In CHDM, the higher redshift damped Ly$\alpha$ systems are predicted to have lower masses, which can be checked by measuring the velocity widths of the associated metal line systems.

The aim of this review article is to give a comprehensive description of the scaling properties detected for the distribution of cosmic structures. Due to the great variety of statistical methods to describe the large-scale structure of the Universe, I will mainly concentrate on those methods which reveal remarkable regularities and scaling in the structure of the Universe. Although in most cases I prefer not to enter into the technical aspects of how implementing such methods, more details will be furnishes about the description of galaxy clustering in terms of fractal concepts. Despite recent claims for a Universe, which behaves like a fractal at arbitrarily large scales, I will show that the fractal language can be usefully employed to disprove this picture. The emerging scenario is that of a Universe, which behaves like a self-similar structure at small scales, where fractality is dynamically generated by non-linear gravitational clustering, while preserving large-scale homogeneity. Nevertheless, even at scales $\magcir 10\hm$, where gravity still acts linearly, the distribution of galaxy clusters shows remarkable scale-invariant features, which could give precise hints about the initial conditions for the evolution of the large-scale structure of the Universe.

We apply fractal analysis methods to investigate the scaling properties in the Abell and ACO catalogs of rich galaxy clusters. We also discuss different technical aspects of the method when applied to data sets with small number of points as the cluster catalogs. Results are compared with simulations based on the Zel'dovich approximation. We limit our analysis to scales less than 100 $\hm$. The cluster distribution show a scale invariant multifractal behavior in a limited scale range. For the Abell catalog this range is 15--60$\hm$, while for the ACO sample it extends to smaller scales. Despite this difference in the extension of the scale--range where scale--invariant clustering takes place, both samples are characterized by remarkably similar multifractal spectra in the corresponding scaling regime. In particular, the correlation dimension turns out to be $D_2\simeq 2.2$ for both Abell and ACO clusters. Although it is difficult to point out the scale at which homogeneity is reached with the present size of these redshift surveys, our results indicate that the cluster distribution shows a tendency to homogeneity at large scales, disproving the picture of a pure scale invariant fractal structure extending to arbitrarily large distances.

We re-estimate the angular 2-point galaxy correlation function from the Lick galaxy catalogue. We argue that the large-scale gradients observed in the Lick catalogue are dominated by real clustering and therefore they should not be subtracted prior to the estimation of the 2-p correlation function. We find that if no such correction is introduced the galaxy correlations are perfectly consistent with the those found in the APM survey. Thus, the long standing discrepancy between the Lick and APM angular correlations is lifted.

We show how to simulate the clustering of rich clusters of galaxies using a technique based on the Zel'dovich approximation. This method well reproduces the spatial distribution of clusters obtainable from full N-body simulations at a fraction of the computational cost. We use an ensemble of large--scale simulations to assess the level and statistical significance of cluster clustering in open, tilted and flat versions of the Cold Dark Matter (CDM) model, as well as a model comprising a mixture of Cold and Hot Dark Matter (CHDM). We find the open and flat CDM models are excluded by the data. The tilted CDM model with a slight tilt is in marginal agreement, while larger tilt produces the right amount of clustering; CHDM is the best of all our models at reproducing the observations of cluster clustering. We find that {\em all} our models display a systematically weaker relationship between clustering length and mean cluster separation than seems to be implied by observations. We also note that the cluster bias factor, is not constant in any of the models, showing that one needs to be very careful when relating cluster clustering statistics to primordial density fluctuations.

We use high resolution PM N-body simulations to follow the development of non-linear clustering in a flat Universe, dominated by Cold + Hot Dark Matter (CHDM) with 60% of CDM, 30% of HDM and 10% of baryons; a simulation box of 100 Mpc a side ($h=0.5$) is used. We analyze two CHDM simulations with $b =1.5$ (COBE normalization). We also compare them with CDM simulations with $b=1.5$ and $b=1$. We evaluate high-order correlation functions and the void-probability-function (VPF). Correlation functions are obtained both from counts in cells and counts of neighbors. The analysis is made for DM particles and for galaxies, identified as massive halos in the evolved density field. We also check the effects of dynamical evolution and redshift space distortions. We find that clustering of DM particles exhibits deviations from the hierarchical scaling, which decrease somewhat in redshift space. Galaxies follow hierarchical scaling far more closely, with coefficients S_3=2.5 and S_4=7.5, in general agreement with observational results, just marginally affected by redshift distortions and on the choices for the initial spectra. The hierarchical scaling of galaxy clustering is confirmed by the VPF analysis and in all the cases considered a good agreement with observational results is obtained.

We analyze the statistical properties of bubble models for the large-scale distribution of galaxies. To this aim, we realize static simulations, in which galaxies are mostly randomly arranged in the regions surrounding bubbles. As a first test, we realize simulations of the Lick map, by suitably projecting the three-dimensional simulations. In this way, we are able to safely compare the angular correlation function implied by a bubbly geometry to that of the APM sample. We find that several bubble models provide an adequate amount of large-scale correlation, which nicely fits that of APM galaxies. Further, we apply the statistics of the count-in-cell moments to the three-dimensional distribution and compare them with available observational data on variance, skewness and kurtosis. Based on our purely geometrical constructions, we find that a well defined hierarchical scaling of higher order moments up to scales $\sim 70\hm$. The overall emerging picture is that the bubbly geometry is well suited to reproduce several aspects of large-scale clustering.

We perform a detailed investigation of the statistical properties of the projected distribution of galaxy clusters obtained in Cold Dark Matter (CDM) models with both Gaussian and skewed primordial density fluctuations. We use N-body simulations to construct a set artificial Lick maps. An objective cluster--finding algorithm is used to identify clusters of different richness. For Gaussian models, the overall number of clusters is too small in the standard CDM case, but a model with higher normalisation fares much better; non--Gaussian models with negative skewness also fit faily well. We apply several statistical tests to compare real and simulated cluster samples, such as the 2-point correlation function, the minimal spanning tree construction, the multifractal analysis and the skewness of cell counts. The emerging picture is that Gaussian models, even with a higher normalization, are in trouble. Skew-positive models are also ruled out, while skew-negative models can reproduce the observed clustering of galaxy clusters in the CDM framework.

We generate artificial Lick maps using $N$--body simulations and compare the angular correlation function, $w(\vartheta)$, measured from the simulations with the APM correlation. For the Gaussian CDM model, neither the standard biassed model nor a more evolved model (as suggested by the COBE data), reproduce the correlations on large angular scales. We come to a similar conclusion about CDM models with positively skewed initial fluctuation distributions. In contrast, models with initially negatively skewed fluctuations produce a $w(\vartheta)$ that declines much more gently on large scales. Such models are therefore in principle capable of reconciling the lack of large--scale power of the CDM spectrum with the observed clustering of APM galaxies.