Search

The chemical enrichment of the Universe at all scales is related to stellar winds and explosive supernovae phenomena. Metals produced by stars and later spread at the mega-parsec scale through the intra-cluster medium (ICM) become a fossil record of the chemical enrichment of the Universe and of the dynamical and feedback mechanisms determining their circulation. As demonstrated by the results of the soft X-ray spectrometer onboard Hitomi, high resolution X-ray spectroscopy is the path to to differentiate among the models that consider different metal production mechanisms, predict the outcoming yields, and are function of the nature, mass, and/or initial metallicity of their stellar progenitor. Transformational results shall be achieved through improvements in the energy resolution and effective area of X-ray observatories to detect rare metals (e.g. Na, Al) and constrain yet uncertain abundances (e.g. C, Ne, Ca, Ni). The X-ray Integral Field Unit (X-IFU) instrument onboard the next-generation European X-ray observatory Athena is expected to deliver such breakthroughs. Starting from 100 ks of synthetic observations of 12 abundance ratios in the ICM of four simulated clusters, we demonstrate that the X-IFU will be capable of recovering the input chemical enrichment models at both low ($z = 0.1$) and high ($z = 1$) redshifts, while statistically excluding more than 99.5% of all the other tested combinations of models. By fixing the enrichment models which provide the best fit to the simulated data, we also show that the X-IFU will constrain the slope of the stellar initial mass function within $\sim$12%. These constraints will be key ingredients in our understanding of the chemical enrichment of the Universe and its evolution., Comment: 17 pages, 6 figures, accepted for publication in A&A

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 profiles as traced by substructures, while DM particles profiles are characterized by an almost universal behavior, thus suggesting that the level of pseudo-entropy could represent a potential low-scatter 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 studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range $0

INAF Trieste Astronomical Observatory (OATs) has a long tradition in information technology applied to Astronomical and Astrophysical use cases, particularly for what regards computing for data reduction, analysis and simulations; data and archives management; space missions data processing; design and software development for ground-based instruments. The ensemble of these activities, in the last years, pushed the need to acquire new computing resources and technologies and to deep competences in theirs management. In this paper we describe INAF-OATs computing centre technological stuff, our involvement in different EU Projects both in the path of building of EOSC, the European Open Science Cloud; in the design and prototyping of new Exascale supercomputers in Europe and the main research activities carried on using our computing centre.

The aim of this work is to quantitatively evaluate the impact of computation on the energy consumption on ARM MPSoC platforms, exploiting CPUs, embedded GPUs and FPGAs. One of them possibly represents the future of High Performance Computing systems: a prototype of an Exascale supercomputer. Performance and energy measurements are made using a state-of-the-art direct $N$-body code from the astrophysical domain. We provide a comparison of the time-to-solution and energy delay product metrics, for different software configurations. We have shown that FPGA technologies can be used for application kernel acceleration and are emerging as a promising alternative to "traditional" technologies for HPC, which purely focus on peak-performance than on power-efficiency., Comment: 10 pages, 5 figure, 2 tables; The final publication will be available at IOS Press

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., Comment: 22 pages, 15 figures

The architecture of Exascale computing facilities, which involves millions of heterogeneous processing units, will deeply impact on scientific applications. Future astrophysical HPC applications must be designed to make such computing systems exploitable. The ExaNeSt H2020 EU-funded project aims to design and develop an exascale ready prototype based on low-energy-consumption ARM64 cores and FPGA accelerators. We participate to the design of the platform and to the validation of the prototype with cosmological N-body and hydrodynamical codes suited to perform large-scale, high-resolution numerical simulations of cosmic structures formation and evolution. We discuss our activities on astrophysical applications to take advantage of the underlying architecture., Comment: submitted to ASP

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., Comment: 18 pages; submitted to MNRAS

[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., Comment: 17 pages, 7 figures, 7 tables, 2 appendices, submitted to New Journal of Physics

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], Comment: 17 pages, 11 figures, 4 tables, to appear in MNRAS

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., Comment: 22 pages, 20 figures, minor revisions, accepted for publication in MNRAS

We present results from the first cosmological simulations which study the onset of primordial, metal-free (population III), cosmic star formation and the transition to the present-day, metal-rich star formation (population II-I), including molecular (H$_2$, HD, etc.) evolution, tracing the injection of metals by supernov{\ae} into the surrounding intergalactic medium and following the change in the initial stellar mass function (IMF) according to the metallicity of the corresponding stellar population. Our investigation addresses the role of a wide variety of parameters (critical metallicity for the transition, IMF slope and range, SN/pair-instability SN metal yields, star formation threshold, resolution, etc.) on the metal-enrichment history and the associated transition in the star formation mode. All simulations present common trends. Metal enrichment is very patchy, with rare, unpolluted regions surviving at all redshifts, inducing the simultaneous presence of metal-free and metal-rich star formation regimes. As a result of the rapid pollution within high-density regions due to the first SN/pair-instability SN, local metallicity is quickly boosted above the critical metallicity for the transition. The population III regime lasts for a very short period during the first stages of star formation ($\sim 10^7\,\rm yr$), and its average contribution to the total star formation rate density drops rapidly below $\sim 10^{-3}-10^{-2}$.

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., Comment: 22 pages, 24 figures, 2 tables, to appear in MNRAS

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., Comment: 13 pages, 9 figures, accepted to MNRAS

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., Comment: 6 pages, 4 figures, accepted for publication in MNRAS (Letters)

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., Comment: 9 pages, 7 figures, to appear in MNRAS

We present a study on the origin of the metallicity evolution of the intra-cluster medium (ICM) by applying a semi-analytic model of galaxy formation to N-body/SPH (smoothed particle hydrodynamic) non-radiative numerical simulations of clusters of galaxies. The semi-analytic model includes gas cooling, star formation, supernovae feedback and metal enrichment, and is linked to the diffuse gas of the underlying simulations so that the chemical properties of gas particles are dynamically and consistently generated from stars in the galaxies. This hybrid model let us have information on the spatial distribution of metals in the ICM. The results obtained for a set of clusters with virial masses of ~1.5*10^15 h^{-1} M_sun contribute to the theoretical interpretation of recent observational X-ray data, which indicate a decrease of the average iron content of the intra-cluster gas with increasing redshift. We find that this evolution arises mainly as a result of a progressive increase of the iron abundance within ~0.15 R_vir. The clusters have been considerably enriched by z~1 with very low contribution from recent star formation. Low entropy gas that has been enriched at high redshift sinks to the cluster centre contributing to the evolution of the metallicity profiles., Comment: 11 pages, 6 figures, accepted for publication in MNRAS

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., Comment: 25 pages, 11 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 18; 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. Bleeker

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], Comment: 15 pages, accepted by MNRAS

(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., Comment: (e.g.: 17 pages, 8 figures, accepted for publication in the Astrophysical Journal, updated discussion to match published version-new section:6.3)

A PopIII/Pop II transition from massive to normal stars is predicted to occur when the metallicity of the star forming gas crosses the critical range Z_cr = 10^(-5 +/- 1) Z_sun. To investigate the cosmic implications of such process we use numerical simulations which follow the evolution, metal enrichment and energy deposition of both Pop III and Pop II stars. We find that: (i) due to inefficient heavy element transport by outflows and slow "genetic" transmission during hierarchical growth, large fluctuations around the average metallicity arise; as a result Pop III star formation continues down to z=2.5, but at a low peak rate of 10^-5 M_sun yr^-1 Mpc^-3 occurring at z~6 (about 10^-4 of the PopII one); (ii) Pop III star formation proceeds in a "inside-out" mode in which formation sites are progressively confined at the periphery of collapsed structures, where the low gas density and correspondingly long free-fall timescales result in a very inefficient astration. These conclusions strongly encourage deep searches for pristine star formation sites at moderate (2

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., Comment: to appear on MNRAS

This submission has been withdrawn by arXiv administrators because it is a duplicate of 0704.2182.

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., Comment: 15 pages, 17 figures, to appear in MNRAS

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., Comment: 15 pages, 10 figures, to appear in MNRAS; version with higher resolution figures available at http://www.daut.univ.trieste.it/borgani/ClusterRuns/paper_rv1.ps.gz

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), Comment: 10 pages, to appear in A&A

(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., Comment: 13 pages, MNRAS in press

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., Comment: 14 pages, 12 figures, MNRAS, in press. Figures 3 and 4 now contain more recent observational data. Other minor revisions according to referee's comments

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., Comment: 4 pages, 3 figures

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., Comment: to appear in MNRAS Letters

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., Comment: 20 pages, 14 figures, to appear in MNRAS

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., Comment: To appear in MNRAS: 18 pages, 13 figures. Paper with high-resolution figures available at http://www.daut.univ.trieste.it/borgani

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., Comment: 5 page, A&A in press