Showing 1–50 of 62 results for author: Tornatore, L

Dynamical friction and evolution of black holes in cosmological simulations: a new implementation in OpenGadget3

Authors: Alice Damiano, Milena Valentini, Stefano Borgani, Luca Tornatore, Giuseppe Murante, Antonio Ragagnin, Cinthia Ragone-Figueroa, Klaus Dolag

Abstract: We implement a sub-resolution prescription for the unresolved dynamical friction onto black holes (BHs) in the OpenGadget3 code. We carry out cosmological simulations of a volume of 16 cMpc3 and zoom-ins of a galaxy group and of a galaxy cluster. The advantages of our new technique are assessed in comparison to commonly adopted methods to hamper spurious BH displacements, i.e. repositioning onto a… ▽ More We implement a sub-resolution prescription for the unresolved dynamical friction onto black holes (BHs) in the OpenGadget3 code. We carry out cosmological simulations of a volume of 16 cMpc3 and zoom-ins of a galaxy group and of a galaxy cluster. The advantages of our new technique are assessed in comparison to commonly adopted methods to hamper spurious BH displacements, i.e. repositioning onto a local minimum of the gravitational potential and ad-hoc boosting of the BH particle dynamical mass. The newly-introduced dynamical friction correction provides centering of BHs on host halos which is at least comparable with the other techniques. It predicts half as many merger events with respect to the repositioning prescription, with the advantage of being less prone to leave sub-structures without any central BH. Simulations featuring our dynamical friction prescription produce a smaller (by up to 50% with respect to repositioning) population of wandering BHs and final BH masses in good agreement with observations. As for individual BH-BH interactions, our dynamical friction model captures the gradual inspiraling of orbits before the merger occurs. By contrast, the repositioning scheme, in its most classical renditions considered, describes extremely fast mergers, while the dynamical mass misrepresents the BHs' dynamics, introducing numerical scattering between the orbiting BHs. Given its performances in describing the centering of BHs within host galaxies and the orbiting of BH pair before their merging, our dynamical friction correction opens interesting applications for an accurate description of the evolution of BH demography within cosmological simulations of galaxy formation at different cosmic epochs and within different environments. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 19 pages, 10 figures, submitted to A&A

GrGadget: an N-body TreePM relativistic code for cosmological simulations

Authors: Eduardo Quintana-Miranda, Pierluigi Monaco, Luca Tornatore

Abstract: We present the merging of the Particle-Mesh (PM) relativistic Gevolution code with the TreePM Gadget-4 code, with the aim of studying general relativity effects in cosmology. Our code, called GrGadget, is able to track the evolution of metric perturbations in the weak field limit by using Gevolution's implementation of a relativistic PM in the Poisson gauge. To achieve this, starting from Gevoluti… ▽ More We present the merging of the Particle-Mesh (PM) relativistic Gevolution code with the TreePM Gadget-4 code, with the aim of studying general relativity effects in cosmology. Our code, called GrGadget, is able to track the evolution of metric perturbations in the weak field limit by using Gevolution's implementation of a relativistic PM in the Poisson gauge. To achieve this, starting from Gevolution we have written a C++ library called libgevolution, that allows a code to access and use the same abstractions and resources that Gevolution uses for its PM-only N-body simulations. The code works under the assumption that particle interactions at short distances can be approximated as Newtonian, so that we can combine the forces computed with a Newtonian Tree with those computed with a relativistic PM. The result is a TreePM simulation code that represents metric perturbations at the scales where they are relevant, while resolving non-linear structures. We validate our code by closely matching Gadget-4 forces, computed with the Tree switched off, with those computed with libgevolution in the Newtonian limit. With GrGadget we obtain a matter power spectrum that is compatible with Newtonian Gadget at small scales and contains GR features at large scales that are consistent with results obtained with Gevolution. We demonstrate that, due to the better resolution of the highly non-linear regime, the representation of the relativistic fields sampled on the mesh improves with respect to the PM-only simulations. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 15 pages, 14 figures

Impact of H$_{\rm 2}$-driven star formation and stellar feedback from low-enrichment environments on the formation of spiral galaxies

Authors: Milena Valentini, Klaus Dolag, Stefano Borgani, Giuseppe Murante, Umberto Maio, Luca Tornatore, Gian Luigi Granato, Cinthia Ragone-Figueroa, Andreas Burkert, Antonio Ragagnin, Elena Rasia

Abstract: The reservoir of molecular gas (H$_{\rm 2}$) represents the fuel for the star formation (SF) of a galaxy. Connecting the star formation rate (SFR) to the available H$_{\rm 2}$ is key to accurately model SF in cosmological simulations of galaxy formation. We investigate how modifying the underlying modelling of H$_{\rm 2}$ and the description of stellar feedback in low-metallicity environments (LMF… ▽ More The reservoir of molecular gas (H$_{\rm 2}$) represents the fuel for the star formation (SF) of a galaxy. Connecting the star formation rate (SFR) to the available H$_{\rm 2}$ is key to accurately model SF in cosmological simulations of galaxy formation. We investigate how modifying the underlying modelling of H$_{\rm 2}$ and the description of stellar feedback in low-metallicity environments (LMF, i.e. low-metallicity stellar feedback) in cosmological, zoomed-in simulations of a Milky Way-size halo influences the formation history of the forming, spiral galaxy and its final properties. We exploit two different models to compute the molecular fraction of cold gas (f$_{\rm H_{\rm 2}}$): $i)$ the theoretical model by Krumholz et al. (2009b) and $ii)$ the phenomenological prescription by Blitz & Rosolowsky (2006). We find that the model adopted to estimate f$_{\rm H_{\rm 2}}$ plays a key role in determining final properties and in shaping the morphology of the galaxy. The clumpier interstellar medium (ISM) and the more complex H$_{\rm 2}$ distribution that the Krumholz et al. (2009b) model predicts result in better agreement with observations of nearby disc galaxies. This shows how crucial it is to link the SFR to the physical properties of the star-forming, molecular ISM. The additional source of energy that LMF supplies in a metal-poor ISM is key in controlling SF at high redshift and in regulating the reservoir of SF across cosmic time. Not only is LMF able to regulate cooling properties of the ISM, but it also reduces the stellar mass of the galaxy bulge. These findings can foster the improvement of the numerical modelling of SF in cosmological simulations. △ Less

Submitted 28 October, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: accepted for publication in MNRAS

Galaxies in the central regions of simulated galaxy clusters

Authors: Antonio Ragagnin, Massimo Meneghetti, Luigi Bassini, Cinthia Ragone-Figueroa, Gian Luigi Granato, Giulia Despali, Carlo Giocoli, Giovanni Granata, Lauro Moscardini, Pietro Bergamini, Elena Rasia, Milena Valentini, Stefano Borgani, Francesco Calura, Klaus Dolag, Claudio Grillo, Amata Mercurio, Giuseppe Murante, Priyamvada Natarajan, Piero Rosati, Giuliano Taffoni, Luca Tornatore, Luca Tortorelli

Abstract: In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedbac… ▽ More In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedback scheme. We study the impact of these different setups on the subhalo abundances, their radial distribution, their density and mass profiles and the relation between the maximum circular velocity, which is a proxy for subhalo compactness. Regardless of the adopted numerical resolution and feedback model, subhaloes with masses Msub < 1e11Msun/h, the most relevant mass-range for galaxy-galaxy strong lensing, have maximum circular velocities ~30% smaller than those measured from strong lensing observations of Bergamini et al. (2019). We also find that simulations with less effective AGN energy feedback produce massive subhaloes (Msub> 1e11 Msun/h ) with higher maximum circular velocity and that their Vmax - Msub relation approaches the observed one. However the stellar-mass number count of these objects exceeds the one found in observations and we find that the compactness of these simulated subhaloes is the result of an extremely over-efficient star formation in their cores, also leading to larger-than-observed subhalo stellar mass. We conclude that simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged from the analysis of Meneghetti et al. (2020) persists. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the LCDM paradigm. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 11 pages, 10 figures, abstract is redacted to fit arXiv character count limit

Journal ref: A&A 665, A16 (2022)

The probability of galaxy-galaxy strong lensing events in hydrodynamical simulations of galaxy clusters

Authors: Massimo Meneghetti, Antonio Ragagnin, Stefano Borgani, Francesco Calura, Giulia Despali, Carlo Giocoli, Gian Luigi Granato, Claudio Grillo, Lauro Moscardini, Elena Rasia, Piero Rosati, Giuseppe Angora, Luigi Bassini, Pietro Bergamini, Gabriel B. Caminha, Giovanni Granata, Amata Mercurio, Robert Benton Metcalf, Priyamvada Natarajan, Mario Nonino, Giada Venusta Pignataro, Cinthia Ragone-Figueroa, Eros Vanzella, Ana Acebron, Klaus Dolag , et al. (5 additional authors not shown)

Abstract: Meneghetti et al. (2020) recently reported an excess of galaxy-galaxy strong lensing (GGSL) in galaxy clusters compared to expectations from the LCDM cosmological model. Theoretical estimates of the GGSL probability are based on the analysis of numerical hydrodynamical simulations in the LCDM cosmology. We quantify the impact of the numerical resolution and AGN feedback scheme adopted in cosmologi… ▽ More Meneghetti et al. (2020) recently reported an excess of galaxy-galaxy strong lensing (GGSL) in galaxy clusters compared to expectations from the LCDM cosmological model. Theoretical estimates of the GGSL probability are based on the analysis of numerical hydrodynamical simulations in the LCDM cosmology. We quantify the impact of the numerical resolution and AGN feedback scheme adopted in cosmological simulations on the predicted GGSL probability and determine if varying these simulation properties can alleviate the gap with observations. We repeat the analysis of Meneghetti et al. (2020) on cluster-size halos simulated with different mass and force resolutions and implementing several independent AGN feedback schemes. We find that improving the mass resolution by a factor of ten and twenty-five, while using the same galaxy formation model that includes AGN feedback, does not affect the GGSL probability. We find similar results regarding the choice of gravitational softening. On the contrary, adopting an AGN feedback scheme that is less efficient at suppressing gas cooling and star formation leads to an increase in the GGSL probability by a factor between three and six. However, we notice that such simulations form overly massive subhalos whose contribution to the lensing cross-section would be significant while their Einstein radii are too large to be consistent with the observations. The primary contributors to the observed GGSL cross-sections are subhalos with smaller masses, that are compact enough to become critical for lensing. The population with these required characteristics appears to be absent in simulations. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 13 pages, 11 figures. Submitted for publication on Astronomy and Astrophysics

Journal ref: A&A 668, A188 (2022)

Dust evolution in zoom-in cosmological simulations of galaxy formation

Authors: Gian Luigi Granato, Cinthia Ragone-Figueroa, Antonela Taverna, Laura Silva, Milena Valentini, Stefano Borgani, Pierluigi Monaco, Giuseppe Murante, Luca Tornatore

Abstract: We present cosmological zoom-in hydro-dynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain size distribution, with the MUPPI star formation and feedback sub-resolution model. Our dust evolutio… ▽ More We present cosmological zoom-in hydro-dynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain size distribution, with the MUPPI star formation and feedback sub-resolution model. Our dust evolution model follows carbon and silicate dust separately. To distinguish differences induced by the chaotic behaviour of simulations from those genuinely due to different simulation set-up, we run each model six times, after introducing tiny perturbations in the initial conditions. With this method, we discuss the role of various dust-related physical processes and the effect of a few possible approximations adopted in the literature. Metal depletion and dust cooling affect the evolution of the system, causing substantial variations in its stellar, gas and dust content. We discuss possible effects on the Spectral Energy Distribution of the significant variations of the size distribution and chemical composition of grains, as predicted by our simulations during the evolution of the galaxy. We compare dust surface density, dust-to-gas ratio and small-to-big grain mass ratio as a function of galaxy radius and gas metallicity predicted by our fiducial run with recent observational estimates for three disc galaxies of different masses. The general agreement is good, in particular taking into account that we have not adjusted our model for this purpose. △ Less

Submitted 7 February, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 21 pages, 15 figures, MNRAS, accepted

Constraining the origin and models of chemical enrichment in galaxy clusters using the Athena X-IFU

Authors: F. Mernier, E. Cucchetti, L. Tornatore, V. Biffi, E. Pointecouteau, N. Clerc, P. Peille, E. Rasia, D. Barret, S. Borgani, E. Bulbul, T. Dauser, K. Dolag, S. Ettori, M. Gaspari, F. Pajot, M. Roncarelli, J. Wilms, C. Noûs

Abstract: 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… ▽ More 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. △ Less

Submitted 9 October, 2020; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 17 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 642, A90 (2020)

On the phase-space structure of galaxy clusters from cosmological simulations

Authors: I. Marini, A. Saro, S. Borgani, G. Murante, E. Rasia, K. Dolag, W. Lin, N. R. Napolitano, A. Ragagnin, L. Tornatore, Y. Wang

Abstract: 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 $σ_\mathrm{v}$ with the density $ρ$ can be used to define the pseudo-entropy $S(r)=σ_\mathrm{v}^2/ρ^{\,2/3}$, who… ▽ More 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 $σ_\mathrm{v}$ with the density $ρ$ can be used to define the pseudo-entropy $S(r)=σ_\mathrm{v}^2/ρ^{\,2/3}$, whose profile is well-described by a simple power-law $S\propto\,r^{\,α}$. 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. △ Less

Submitted 5 November, 2020; v1 submitted 10 July, 2020; originally announced July 2020.

Report number: Monthly Notices of the Royal Astronomical Society, Vol 507, Issue 4, (2021) 5780--5795

The DIANOGA simulations of galaxy clusters: characterizing star formation in proto-clusters

Authors: L. Bassini, E. Rasia, S. Borgani, G. L. Granato, C. Ragone-Figueroa, V. Biffi, A. Ragagnin, K. Dolag, W. Lin, G. Murante, N. R. Napolitano, G. Taffoni, L. Tornatore, Y. Wang

Abstract: 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 hydr… ▽ More 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. △ Less

Submitted 20 July, 2020; v1 submitted 24 June, 2020; originally announced June 2020.

Comments: 20 pages, 19 figures. Accepted for publication in A&A

Journal ref: A&A 642, A37 (2020)

Performance and energy footprint assessment of FPGAs and GPUs on HPC systems using Astrophysics application

Authors: David Goz, Georgios Ieronymakis, Vassilis Papaefstathiou, Nikolaos Dimou, Sara Bertocco, Francesco Simula, Antonio Ragagnin, Luca Tornatore, Igor Coretti, Giuliano Taffoni

Abstract: New challenges in Astronomy and Astrophysics (AA) are urging the need for a large number of exceptionally computationally intensive simulations. "Exascale" (and beyond) computational facilities are mandatory to address the size of theoretical problems and data coming from the new generation of observational facilities in AA. Currently, the High Performance Computing (HPC) sector is undergoing a pr… ▽ More New challenges in Astronomy and Astrophysics (AA) are urging the need for a large number of exceptionally computationally intensive simulations. "Exascale" (and beyond) computational facilities are mandatory to address the size of theoretical problems and data coming from the new generation of observational facilities in AA. Currently, the High Performance Computing (HPC) sector is undergoing a profound phase of innovation, in which the primary challenge to the achievement of the "Exascale" is the power-consumption. The goal of this work is to give some insights about performance and energy footprint of contemporary architectures for a real astrophysical application in an HPC context. We use a state-of-the-art N-body application that we re-engineered and optimized to exploit the heterogeneous underlying hardware fully. We quantitatively evaluate the impact of computation on energy consumption when running on four different platforms. Two of them represent the current HPC systems (Intel-based and equipped with NVIDIA GPUs), one is a micro-cluster based on ARM-MPSoC, and one is a "prototype towards Exascale" equipped with ARM-MPSoCs tightly coupled with FPGAs. We investigate the behavior of the different devices where the high-end GPUs excel in terms of time-to-solution while MPSoC-FPGA systems outperform GPUs in power consumption. Our experience reveals that considering FPGAs for computationally intensive application seems very promising, as their performance is improving to meet the requirements of scientific applications. This work can be a reference for future platforms development for astrophysics applications where computationally intensive calculations are required. △ Less

Submitted 10 April, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Comments: 15 pages, 4 figures, 3 tables; Preprint (V2) submitted to MDPI (Special Issue: Energy-Efficient Computing on Parallel Architectures)

INAF Trieste Astronomical Observatory Information Technology Framework

Authors: S. Bertocco, D. Goz, L. Tornatore, A. Ragagnin, G. Maggio, F. Gasparo, C. Vuerli, G. Taffoni, M. Molinaro

Abstract: 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… ▽ More 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. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Impact of AGN feedback on galaxies and their multiphase ISM across cosmic time

Authors: Milena Valentini, Giuseppe Murante, Stefano Borgani, Gian Luigi Granato, Pierluigi Monaco, Fabrizio Brighenti, Luca Tornatore, Alessandro Bressan, Andrea Lapi

Abstract: 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 t… ▽ More 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. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Comments: Accepted for publication by MNRAS. Comments welcome

Direct N-body application on low-power and energy-efficient parallel architectures

Authors: D. Goz, G. Ieronymakis, V. Papaefstathiou, N. Dimou, S. Bertocco, A. Ragagnin, L. Tornatore, G. Taffoni, I. Coretti

Abstract: 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 astroph… ▽ More 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. △ Less

Submitted 31 October, 2019; originally announced October 2019.

Comments: 10 pages, 5 figure, 2 tables; The final publication will be available at IOS Press

Shall numerical astrophysics step into the era of Exascale computing?

Authors: Giuliano Taffoni, Giuseppe Murante, Luca Tornatore, David Goz, Stefano Borgani, Manolis Katevenis, Nikolaos Chrysos, Manolis Marazakis

Abstract: 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… ▽ More 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. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 3 figures, invited talk for proceedings of ADASS XXVI, accepted by ASP Conference Series

Black hole mass of central galaxies and cluster mass correlation in cosmological hydro-dynamical simulations

Authors: Luigi Bassini, Elena Rasia, Stefano Borgani, Cinthia Ragone-Figueroa, Veronica Biffi, Klaus Dolag, Massimo Gaspari, Gian Luigi Granato, Giuseppe Murante, Giuliano Taffoni, Luca Tornatore

Abstract: 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 merg… ▽ More 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σ$ 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. △ Less

Submitted 27 October, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 16 pages, 17 figures, 3 tables

Journal ref: A&A 630, A144 (2019)

Chemical evolution of disc galaxies from cosmological simulations

Authors: Milena Valentini, Stefano Borgani, Alessandro Bressan, Giuseppe Murante, Luca Tornatore, Pierluigi Monaco

Abstract: 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,… ▽ More 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 $α$-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 $α$-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 $α$-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. △ Less

Submitted 15 February, 2019; originally announced February 2019.

Comments: Accepted for publications in MNRAS. Comments welcome

Direct $N$-body code on low-power embedded ARM GPUs

Authors: David Goz, Sara Bertocco, Luca Tornatore, Giuliano Taffoni

Abstract: This work arises on the environment of the ExaNeSt project aiming at design and development of an exascale ready supercomputer with low energy consumption profile but able to support the most demanding scientific and technical applications. The ExaNeSt compute unit consists of densely-packed low-power 64-bit ARM processors, embedded within Xilinx FPGA SoCs. SoC boards are heterogeneous architectur… ▽ More This work arises on the environment of the ExaNeSt project aiming at design and development of an exascale ready supercomputer with low energy consumption profile but able to support the most demanding scientific and technical applications. The ExaNeSt compute unit consists of densely-packed low-power 64-bit ARM processors, embedded within Xilinx FPGA SoCs. SoC boards are heterogeneous architecture where computing power is supplied both by CPUs and GPUs, and are emerging as a possible low-power and low-cost alternative to clusters based on traditional CPUs. A state-of-the-art direct $N$-body code suitable for astrophysical simulations has been re-engineered in order to exploit SoC heterogeneous platforms based on ARM CPUs and embedded GPUs. Performance tests show that embedded GPUs can be effectively used to accelerate real-life scientific calculations, and that are promising also because of their energy efficiency, which is a crucial design in future exascale platforms. △ Less

Submitted 24 January, 2019; originally announced January 2019.

Comments: 16 pages, 7 figures, 1 table, accepted for publication in the Computing Conference 2019 proceedings

Journal ref: Intelligent Computing. CompCom 2019. Advances in Intelligent Systems and Computing, vol 997. Springer, Cham

Astrophysical code migration into Exascale Era

Authors: David Goz, Sara Bertocco, Luca Tornatore, Giuliano Taffoni

Abstract: The ExaNeSt and EuroExa H2020 EU-funded projects aim to design and develop an exascale ready computing platform prototype based on low-energy-consumption ARM64 cores and FPGA accelerators. We participate in the application-driven design of the hardware solutions and prototype validation. To carry on this work we are using, among others, Hy-Nbody, a state-of-the-art direct N-body code. Core algorit… ▽ More The ExaNeSt and EuroExa H2020 EU-funded projects aim to design and develop an exascale ready computing platform prototype based on low-energy-consumption ARM64 cores and FPGA accelerators. We participate in the application-driven design of the hardware solutions and prototype validation. To carry on this work we are using, among others, Hy-Nbody, a state-of-the-art direct N-body code. Core algorithms of Hy-Nbody have been improved in such a way to increasingly fit them to the exascale target platform. Waiting for the ExaNest prototype release, we are performing tests and code tuning operations on an ARM64 SoC facility: a SLURM managed HPC cluster based on 64-bit ARMv8 Cortex-A72/Cortex-A53 core design and powered by a Mali-T864 embedded GPU. In parallel, we are porting a kernel of Hy-Nbody on FPGA aiming to test and compare the performance-per-watt of our algorithms on different platforms. In this paper we describe how we re-engineered the application and we show first results on ARM SoC. △ Less

Submitted 2 December, 2018; originally announced December 2018.

Comments: 4 pages, 1 figure, 1 table; proceedings of ADASS XXVIII, accepted by ASP Conference Series

Athena X-IFU synthetic observations of galaxy clusters to probe the chemical enrichment of the Universe

Authors: E. Cucchetti, E. Pointecouteau, P. Peille, N. Clerc, E. Rasia, V. Biffi, S. Borgani, L. Tornatore, K. Dolag, M. Roncarelli, M. Gaspari, S. Ettori, E. Bulbul, T. Dauser, J. Wilms, F. Pajot, D. Barret

Abstract: 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 enrichm… ▽ More 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. △ Less

Submitted 24 September, 2018; originally announced September 2018.

Comments: 22 pages, 15 figures

Journal ref: A&A 620, A173 (2018)

Colour-magnitude diagram in simulations of galaxy formation

Authors: Milena Valentini, Alessandro Bressan, Stefano Borgani, Giuseppe Murante, Léo Girardi, Luca Tornatore

Abstract: 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 simulatio… ▽ More 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. △ Less

Submitted 28 July, 2018; v1 submitted 30 April, 2018; originally announced May 2018.

Comments: Accepted for publication in MNRAS. Comments welcome

Dust Evolution in Galaxy Cluster Simulations

Authors: Eda Gjergo, Gian Luigi Granato, Giuseppe Murante, Cinthia Ragone-Figueroa, Luca Tornatore, Stefano Borgani

Abstract: 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 t… ▽ More 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). △ Less

Submitted 12 June, 2018; v1 submitted 18 April, 2018; originally announced April 2018.

Comments: 19 pages, 15 figures, accepted on June 8th, 2018, by MNRAS

Cosmological Simulations in Exascale Era

Authors: D. Goz, L. Tornatore, G. Taffoni, G. Murante

Abstract: 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… ▽ More 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. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Comments: submitted to ASP

The history of chemical enrichment in the intracluster medium from cosmological simulations

Authors: V. Biffi, S. Planelles, S. Borgani, D. Fabjan, E. Rasia, G. Murante, L. Tornatore, K. Dolag, G. L. Granato, M. Gaspari, A. M. Beck

Abstract: 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 incl… ▽ More 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. △ Less

Submitted 27 January, 2017; originally announced January 2017.

Comments: 18 pages; submitted to MNRAS

Simulating realistic disk galaxies with a novel sub-resolution ISM model

Authors: Giuseppe Murante, Pierluigi Monaco, Stefano Borgani, Luca Tornatore, Klaus Dolag, David Goz

Abstract: 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… ▽ More 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. △ Less

Submitted 14 November, 2014; v1 submitted 13 November, 2014; originally announced November 2014.

Comments: 19 pages, 21 figures, MNRAS accepted

Elemental Abundances in Milky Way-like Galaxies from a Hierarchical Galaxy Formation Model

Authors: Gabriella De Lucia, Luca Tornatore, Carlos S. Frenk, Amina Helmi, Julio F. Navarro, Simon D. M. White

Abstract: We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this ap… ▽ More We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope. △ Less

Submitted 29 September, 2014; v1 submitted 29 July, 2014; originally announced July 2014.

Comments: 20 pages, 13 figures, minor changes to match version in press in MNRAS

Simulated star formation rate functions at z~4-7, and the role of feedback in high-z galaxies

Authors: Edoardo Tescari, Antonios Katsianis, Stuart Wyithe, Klaus Dolag, Luca Tornatore, Paramita Barai, Matteo Viel, Stefano Borgani

Abstract: 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 includ… ▽ More 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. △ Less

Submitted 6 January, 2014; v1 submitted 17 December, 2013; originally announced December 2013.

Comments: 18 pages, 3 tables, 8 figures, accepted for publication in MNRAS

Galactic Winds in Cosmological Simulations of the Circumgalactic Medium

Authors: Paramita Barai, Matteo Viel, Stefano Borgani, Edoardo Tescari, Luca Tornatore, Klaus Dolag, Madhura Killedar, Pierluigi Monaco, Valentina D'Odorico, Stefano Cristiani

Abstract: (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 varyin… ▽ More (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. △ Less

Submitted 21 January, 2013; v1 submitted 12 October, 2012; originally announced October 2012.

Comments: Version accepted for publication in MNRAS. Manuscript file with higher resolution Fig. 2 can be found here: http://adlibitum.oats.inaf.it/barai/AllPages/Images-Movies/RadialVarWind.pdf

Baryon Census in Hydrodynamical Simulations of Galaxy Clusters

Authors: Susana Planelles, Stefano Borgani, Klaus Dolag, Stefano Ettori, Dunja Fabjan, Giuseppe Murante, Luca Tornatore

Abstract: 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 coolin… ▽ More 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. △ Less

Submitted 11 February, 2013; v1 submitted 23 September, 2012; originally announced September 2012.

Comments: 17 pages, 8 figures, 4 tables. Additional figure, minor modifications and references added to match the accepted version

How Baryonic Processes affect Strong Lensing properties of Simulated Galaxy Clusters

Authors: Madhura Killedar, Stefano Borgani, Massimo Meneghetti, Klaus Dolag, Dunja Fabjan, Luca Tornatore

Abstract: 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 l… ▽ More 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] △ Less

Submitted 28 August, 2012; originally announced August 2012.

Comments: 19 pages, 18 figures, 1 table, Accepted for publication in MNRAS

Tomographic weak lensing shear spectra from large N-body and hydrodynamical simulations

Authors: Luciano Casarini, Silvio A. Bonometto, Stefano Borgani, Klaus Dolag, Giuseppe Murante, Marino Mezzetti, Luca Tornatore, Giuseppe La Vacca

Abstract: 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… ▽ More 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. △ Less

Submitted 11 May, 2012; v1 submitted 23 March, 2012; originally announced March 2012.

Comments: 13 pages, 19 figures, A&A in press

Lensing and X-ray mass estimates of clusters (SIMULATION)

Authors: E. Rasia, M. Meneghetti, R. Martino, S. Borgani, A. Bonafede, K. Dolag, S. Ettori, D. Fabjan, C. Giocoli, P. Mazzotta, J. Merten, M. Radovich, L. Tornatore

Abstract: [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 resimulat… ▽ More [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. △ Less

Submitted 7 January, 2012; originally announced January 2012.

Comments: 17 pages, 7 figures, 7 tables, 2 appendices, submitted to New Journal of Physics

The effects of baryons on the halo mass function

Authors: Weiguang Cui, Stefano Borgani, Klaus Dolag, Giuseppe Murante, Luca Tornatore

Abstract: 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 o… ▽ More 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), $Δ_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_{Δ_c} / \hMsun)\geq 13.5$. In this range, mass increase in the hydrodynamical simulations is of about 4-5 per cent at $Δ_c=500$ and $\sim 1$ - 2 per cent at $Δ_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 $Δ_c = 200$, and $\simeq 7$ per cent at $Δ_c=500$, with $\sim 10$ - 20 per cent differences reached at $Δ_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) △ Less

Submitted 4 April, 2012; v1 submitted 13 November, 2011; originally announced November 2011.

Comments: 11 pages, 7 figures, minor changes, version accepted by MNRAS

X-ray mass proxies from hydrodynamic simulations of galaxy clusters (paper I)

Authors: D. Fabjan, S. Borgani, E. Rasia, A. Bonafede, K. Dolag, G. Murante, L. Tornatore

Abstract: 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… ▽ More 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] △ Less

Submitted 16 February, 2011; v1 submitted 14 February, 2011; originally announced February 2011.

Comments: 17 pages, 11 figures, 4 tables, to appear in MNRAS

Cosmic evolution of the CIV in high-resolution hydrodynamic simulations

Authors: E. Tescari, M. Viel, V. D'Odorico, S. Cristiani, F. Calura, S. Borgani, L. Tornatore

Abstract: 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 $Λ$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… ▽ More 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 $Λ$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. △ Less

Submitted 22 November, 2010; v1 submitted 9 July, 2010; originally announced July 2010.

Comments: 22 pages, 20 figures, minor revisions, accepted for publication in MNRAS

A decline and fall in the future of Italian Astronomy?

Authors: Angelo Antonelli, Vincenzo Antonuccio-Delogu, Andrea Baruffolo, Stefano Benetti, Simone Bianchi, Andrea Biviano, Annalisa Bonafede, Marco Bondi, Stefano Borgani, Angela Bragaglia, Massimo Brescia, John Robert Brucato, Gianfranco Brunetti, Riccardo Brunino, Michele Cantiello, Viviana Casasola, Rossella Cassano, Alberto Cellino, Gabriele Cescutti, Andrea Cimatti, Andrea Comastri, Edvige Corbelli, Giovanni Cresci, Serena Criscuoli, Stefano Cristiani , et al. (67 additional authors not shown)

Abstract: 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.… ▽ More 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. △ Less

Submitted 8 July, 2010; originally announced July 2010.

Comments: Also available at http://adoptitaastronom.altervista.org/index.html

The transition from population III to population II-I star formation

Authors: U. Maio, B. Ciardi, K. Dolag, L. Tornatore, S. Khochfar

Abstract: 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æ into the surrounding intergalactic medium and following the change in the i… ▽ More 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æ 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}$. △ Less

Submitted 25 March, 2010; originally announced March 2010.

The impact of feedback on the low redshift Intergalactic Medium

Authors: Luca Tornatore, Stefano Borgani, Matteo Viel, Volker Springel

Abstract: 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 (… ▽ More 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. △ Less

Submitted 13 November, 2009; v1 submitted 3 November, 2009; originally announced November 2009.

Comments: accepted by MNRAS; replaced to match the accepted version

Journal ref: Mon.Not.Roy.Astron.Soc.402:1911,2010

Simulating the effect of AGN feedback on the metal enrichment of galaxy clusters

Authors: D. Fabjan, S. Borgani, L. Tornatore, A. Saro, G. Murante, K. Dolag

Abstract: 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 includ… ▽ More 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. △ Less

Submitted 3 September, 2009; originally announced September 2009.

Comments: 22 pages, 24 figures, 2 tables, to appear in MNRAS

Lyman Alpha Emitter Evolution in the Reionization Epoch

Authors: P. Dayal, A. Ferrara, A. Saro, R. Salvaterra, S. Borgani, L. Tornatore

Abstract: 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… ▽ More 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. △ Less

Submitted 2 November, 2009; v1 submitted 2 July, 2009; originally announced July 2009.

Comments: 13 pages, 9 figures, accepted to MNRAS

The onset of star formation in primordial haloes

Authors: Umberto Maio, Benedetta Ciardi, Naoki Yoshida, Klaus Dolag, Luca Tornatore

Abstract: Star formation remains an unsolved problem in astrophysics. Numerical studies of large-scale structure simulations cannot resolve the whole process and their approach usually assumes that only gas denser than a typical threshold can host and form stars. We investigate the onset of cosmological star formation and compare several very-high-resolution, three-dimensional, N-body/SPH simulations that… ▽ More Star formation remains an unsolved problem in astrophysics. Numerical studies of large-scale structure simulations cannot resolve the whole process and their approach usually assumes that only gas denser than a typical threshold can host and form stars. We investigate the onset of cosmological star formation and compare several very-high-resolution, three-dimensional, N-body/SPH simulations that include non-equilibrium, atomic and molecular chemistry, star formation prescriptions, and feedback effects. We study how primordial star formation depends on gas density thresholds, cosmological parameters and initial set-ups. For mean-density initial conditions, we find that standard low-density star-formation threshold (0.2 h^2/cm3) models predict the onset of star formation at z~25-31, depending on the adopted cosmology. In these models stars are formed regardless of the time between the moment when the threshold is reached and the effective runaway collapse. At high redshift, this time interval represents a significant fraction of the Hubble time and thus this assumption can induce large artificial off-sets to the onset of star formation. Choosing higher density thresholds (135 h^2/cm3) allows the entire cooling process to be followed, and the onset of star formation is then estimated to be at redshift z~12-16. When isolated, rare, high-density peaks are considered, the chemical evolution is much faster and the first star formation episodes occur at z > 40, almost regardless of the choice for the density threshold. These results could have implications for the formation redshift of the first cosmological objects, as inferred from direct numerical simulations of mean-density environments, and on the studies of the reionization history of the universe. △ Less

Submitted 24 June, 2009; v1 submitted 2 June, 2009; originally announced June 2009.

Comments: 10 pages, 1 table, 5 figures; in press. Minor changes done

Damped Lyman-alpha systems in high-resolution hydrodynamical simulations

Authors: Edoardo Tescari, Matteo Viel, Luca Tornatore, Stefano Borgani

Abstract: 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… ▽ More 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. △ Less

Submitted 22 April, 2009; originally announced April 2009.

Comments: 21 pages, 18 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 397 (2009) 411

Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations

Authors: Robert P. C. Wiersma, Joop Schaye, Tom Theuns, Claudio Dalla Vecchia, Luca Tornatore

Abstract: (Abridged) We present an implementation of stellar evolution and chemical feedback for smoothed particle hydrodynamics (SPH) simulations. We consider the timed release of individual elements by both massive (Type II supernovae and stellar winds) and intermediate mass stars (Type Ia supernovae and asymptotic giant branch stars). We illustrate the results of our method using a suite of cosmologica… ▽ More (Abridged) We present an implementation of stellar evolution and chemical feedback for smoothed particle hydrodynamics (SPH) simulations. We consider the timed release of individual elements by both massive (Type II supernovae and stellar winds) and intermediate mass stars (Type Ia supernovae and asymptotic giant branch stars). We illustrate the results of our method using a suite of cosmological simulations that include new prescriptions for radiative cooling, star formation, and galactic winds. Radiative cooling is implemented element-by-element, in the presence of an ionizing radiation background, and we track all 11 elements that contribute significantly to the radiative cooling. We contrast two reasonable definitions of the metallicity of a resolution element and find that while they agree for high metallicities, there are large differences at low metallicities. We argue the discrepancy is indicative of the lack of metal mixing caused by the fact that metals are stuck to particles. We argue that since this is a (numerical) sampling problem, solving it using a poorly constrained physical process such as diffusion could have undesired consequences. We demonstrate that the two metallicity definitions result in redshift z = 0 stellar masses that can differ by up to a factor of two, because of the sensitivity of the cooling rates to the elemental abundances. We find that by z = 0 most of the metals are locked up in stars. The gaseous metals are distributed over a very wide range of gas densities and temperatures. The shock-heated warm-hot intergalactic medium has a relatively high metallicity of ~ 10^-1 Z_sun that evolves only weakly and is therefore an important reservoir of metals. △ Less

Submitted 6 July, 2009; v1 submitted 10 February, 2009; originally announced February 2009.

Comments: 29 pages (17 pages main text, 12 pages appendices) 25 figures. Accepted for publication in MNRAS, minor changes to submitted version

Journal ref: MNRAS 399 (2009) 574-600

Simulating the formation of a proto-cluster at z~2

Authors: A. Saro, S. Borgani, L. Tornatore, G. De Lucia, K. Dolag, G. Murante

Abstract: 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… ▽ More 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. △ Less

Submitted 15 October, 2008; originally announced October 2008.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS (Letters)

Evolution of the metal content of the intra-cluster medium with hydrodynamical simulations

Authors: D. Fabjan, L. Tornatore, S. Borgani, A. Saro, K. Dolag

Abstract: 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), Krou… ▽ More 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. △ Less

Submitted 15 February, 2008; originally announced February 2008.

Comments: 9 pages, 7 figures, to appear in MNRAS

On the Dynamical Origin of the ICM Metallicity Evolution

Authors: S. A. Cora, L. Tornatore, P. Tozzi, K. Dolag

Abstract: 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… ▽ More 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. △ Less

Submitted 7 February, 2008; originally announced February 2008.

Comments: 11 pages, 6 figures, accepted for publication in MNRAS

The chemical enrichment of the ICM from hydrodynamical simulations

Authors: S. Borgani, D. Fabjan, L. Tornatore, S. Schindler, K. Dolag, A. Diaferio

Abstract: 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 c… ▽ More 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. △ Less

Submitted 7 January, 2008; originally announced January 2008.

Comments: 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

How does gas cool in DM halos?

Authors: M. Viola, P. Monaco, S. Borgani, G. Murante, L. Tornatore

Abstract: 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 ev… ▽ More 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] △ Less

Submitted 12 October, 2007; originally announced October 2007.

Comments: 15 pages, accepted by MNRAS

X-MAS2: Study Systematics on the ICM Metallicity Measurements

Authors: E. Rasia, P. Mazzotta, H. Bourdin, S. Borgani, L. Tornatore, S. Ettori, K. Dolag, L. Moscardini

Abstract: (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 form… ▽ More (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. △ Less

Submitted 1 November, 2007; v1 submitted 17 July, 2007; originally announced July 2007.

Comments: (e.g.: 17 pages, 8 figures, accepted for publication in the Astrophysical Journal, updated discussion to match published version-new section:6.3)

Population III stars: hidden or disappeared ?

Authors: L. Tornatore, A. Ferrara, R. Schneider

Abstract: 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 inefficien… ▽ More 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<z<5) redshifts where metal free stars are likely to be hidden. △ Less

Submitted 10 July, 2007; originally announced July 2007.

Comments: to appear in MNRAS

Chemical enrichment of galaxy clusters from hydrodynamical simulations

Authors: L. Tornatore, S. Borgani, K. Dolag, F. Matteucci

Abstract: 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… ▽ More 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. △ Less

Submitted 12 May, 2007; originally announced May 2007.

Comments: to appear on MNRAS