Thermodynamical properties of the ICM from hydrodynamical simulations

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.Comment: 26 pages, 12 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 13; work done by an international team at the International Space Science Institute (ISSI), Bern, organised by J.S. Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke

The baryon fraction in 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.Comment: 10 pages, to appear in MNRA

Joint deprojection of Sunyaev-Zeldovich and X-ray images of galaxy clusters

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)Comment: 17 pages, 15 figures, accepted by MNRA

Constraining cosmological models with cluster power spectra

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 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.Comment: LaTeX, 27 pages, 8 figures (EPS). Revised version (title changed, CHDM model added, discussion expanded). Accepted by New

Pointing to the minimum scatter: the generalized scaling relations for galaxy clusters

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.Comment: 6 pages, 3 figures. MNRAS in pres

Radiative feedback and cosmic molecular gas: the role of different radiative sources

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.Comment: 19 pages on MNRA

Reconstructing mass profiles of simulated galaxy clusters by combining Sunyaev-Zeldovich and X-ray images

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.Comment: 13 pages, 11 figures, submitted to MNRA