Search

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., accepted for publication in MNRAS

INTRODUCTION Over the past two decades, ground- and space-based observations have unveiled thousands exoplanets and planetary systems around other stars in our Galaxy. About 5000 exoplanets are currently confirmed, in large part detected as transits by the Kepler and TESS missions. Launched in 2019, PLATO mission represents the first-step characterisation towards the understanding of the structural properties of these planets. Nevertheless, a significant boost for the detection of transiting Earth-analogues around bright stars is expected from the PLATO mission. However, it's only through remote atmospheric spectroscopy of potentially habitable rocky planets, that one of the main goals of exoplanetary science, the quest for life outside the Solar System, can be tackled. This observational challenge should be partly within reach of the recently launched JWST and the next ground-based astronomical observatory, E-ELT. To accomplish the demanding task of searching for and deciphering spectral signatures, a thorough and holistic observational and theoretical characterization of carefully selected rocky exoplanets is required.The selection, among the observationally reachable targets for high-resolution spectroscopy of thin atmospheres, requires habitability studies with climate models. These simulations will enable the identification of those exoplanets with the largest chance of potentially hosting a surface diffuse life, i.e. with the largest habitability, that must be evaluated over a wide range of mostly unknown conditions. A considerable effort of modelization that exploits all available observations will be needed in order to assess the global physical characterization of the selected exoplanets, and in particular precisely of their potential surface climate and habitability. THE MODEL Here we present EOS-ESTM, a flexible climate model aimed at simulating the surface and atmospheric conditions that characterize habitable planets. The model allows one to perform a fast exploration of the parameter space representative of planetary quantities, including those currently not measurable in rocky exoplanets. EOS-ESTM has been built up starting from ESTM (Vladilo et al. 2013, 2015), a seasonal-latitudinal EBM featuring an advanced treatment of surface and cloud components and a 2D (vertical and latitudinal) treatment of the energy transport. The main features of the model that we have implemented can be summarised as follows. Firstly, we have calculated the atmospheric radiative transfer using EOS (Simonetti et al. 2022), a procedure tailored for atmospheres of terrestrial-type planets, based on the opacity calculator HELIOS-K (Grimm & Heng 2015; Grimm et al. 2021) and the radiative transfer code HELIOS (Malik et al. 2017, 2019). Thanks to EOS, the ESTM radiative transfer can be now calculated for a variety of atmospheres with different bulk and greenhouse compositions, illuminated by stars with different SEDs. Then, we have upgraded the parameterizations that describe the clouds properties. New equations have been introduced for the albedo of the clouds and its dependence on the albedo of the underlying surface. The clouds coverage over ice is now a function of the global planetary ice coverage. A specific treatment for the transmittance and OLR forcing of clouds at very low temperature has been introduced. Lastly, we have introduced a generalized logistic function to estimate the ice coverage as a function of mean zonal surface temperature. Based on a detailed study of the ice distribution on Earth, the adopted algorithm discriminates between ice over lands and oceans. The albedo and thermal capacity of transitional ice is now estimated using the fractional ice coverage. RESULTS With the aim of providing a reference model for studies of habitable planets, we calibrated EOS-ESTM using a large set of Earth satellite and reanalysis data. The reference Earth model satisfies a variety of diagnostic tests, including mean latitudinal profiles of surface temperature (Figure 1), TOA albedo, OLR and ice coverage. Fig. 1. Mean annual latitude profile of surface temperature predicted by the reference Earth model The temperature profile is compared with ERA5 temperatures averaged in the period 2005-2015 (blue dots). To test the consistency of EOS-ESTM with previous studies of non-terrestrial climate conditions we performed a series of comparisons with a hierarchy of climate models, varying the levels of insolation (Figure 2), the stellar spectrum and planetary parameters (radius and rotation rate). Fig. 2. Comparison of global and annual mean surface temperature obtained from different climate Earth's models by increasing the solar constant. Red, solid line: EOS-ESTM (this work). Black, solid line: 3D model CAM4 (Wolf & Toon, 2015). Cyan, solid line: 3D model CAM3 (Wolf & Toon, 2014). Green, solid line: 3D model by Leconte et al. (2013). The application of EOS-ESTM to the case of a CO2-dominated atmosphere in maximum greenhouse conditions (Kasting et al. 1993) yields a detailed description of the transition to a snowball state that takes place when the insolation decreases in the proximity of the outer edge of the HZ. Thanks to the flexibility of our model we can explore how this transition develops in different planetary conditions (e.g. rotation rate, Figure 3). Fig. 3. Dependence on planetary rotation period rotation period, Prot, of the fractional ice coverage calculated at the outer edge of the HZ. The results were obtained for an Earth-like planet with a CO2-dominated, maximum greenhouse atmosphere, the remaining parameters being fixed to Earth values. REFERENCES Grimm S. L., Heng K., 2015, HELIOS-K: Opacity Calculator for Radiative Transfer (ascl:1503.004) Grimm S. L., et al., 2021, ApJS, 253, 30 Kasting J. F., Whitmire D. P., Reynolds R. T., 1993, Icarus, 101, 108 Leconte J., Forget F., Charnay B., Wordsworth R., Pottier A., 2013, Nature, 504, 268 Malik M., et al., 2017, AJ, 153, 56 Malik M., Kitzmann D., Mendonça J. M., Grimm S. L., Marleau G.-D., Linder E. F., Tsai S.-M., Heng K., 2019, AJ, 157, 170 Simonetti P., Vladilo G., Silva L., Maris M., Ivanovski S. L., Biasiotti L., Malik M., von Hardenberg J., 2022, ApJ, 925, 105 Vladilo G., Murante G., Silva L., Provenzale A., Ferri G., Ragazzini G., 2013, ApJ, 767, 65 Vladilo G., Silva L., Murante G., Filippi L., Provenzale A., 2015, ApJ, 804,50 Wolf E. T., Toon O. B., 2014, Geophys. Res. Lett., 41, 167 Wolf E. T., Toon O. B., 2015, Journal of Geophysical Research (Atmospheres), 120, 5775

Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalising the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume., Comment: 19 pages, 15 figures, accepted for publication in MNRAS

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., 11 pages, 10 figures, abstract is redacted to fit arXiv character count limit

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., Comment: 13 pages, 11 figures. Submitted for publication on Astronomy and Astrophysics

THE ''LIfe in Space'' project was funded in the wake of the Italian Space Agency's proposal for the development of a network of institutions and laboratories conceived to implement Italian participation in space astrobiology experiments. Of primary concern for this project is the study of the origin of life in the Universe, a focus that will promote investigation into prebiotic chemistry in various possible scenarios, whether in polar or nonpolar solvents (e.g., Titan's environment). Such results will link with study of the effects of simulated space conditions on possible chemical bio- signatures. The limits of life as we know it will be investigated in ground-based experiments with microorganisms that have already demonstrated their resistance to extreme environments on Earth and to real or simulated space conditions. The potential survival of microorganisms will also be examined with up-to-date molecular methods. The ability of some microorganisms to produce atmospheric and surface biosignatures when exposed to simulated conditions will be tested and compared with the possible existence of bio- signatures on potentially habitable exoplanets. Furthermore, the search for potentially habitable exoplanets, with space- based observational methods, will be optimized by way of dedicated climate models with the capacity to predict the detectability of atmospheric biosignatures for a broad range of planetary conditions. The project embraces the four most important topics in astrobiological research, as listed below, along with relevant contributions from the participating Italian institutions. Origins and evolution of organic compounds of biological significance in space (comets, asteroids, rocky planets, and moons); Prebiotic syntheses, origin of life, and early life; The limits of life and biological habitability: origin, evolution and adaptation of life in extreme environments on Earth and in space; Biomarkers for life detection in the Solar System and on exoplanets.

Using the DIANOGA hydrodynamical zoom-in simulation set of galaxy clusters, we analyze the dynamics traced by stars belonging to the Brightest Cluster Galaxies (BCGs) and their surrounding diffuse component, forming the intracluster light (ICL), and compare it to the dynamics traced by dark matter and galaxies identified in the simulations. We compute scaling relations between the BCG and cluster velocity dispersions and their corresponding masses (i.e. $M_\mathrm{BCG}^{\star}$- $\sigma_\mathrm{BCG}^{\star}$, $M_{200}$- $\sigma_{200}$, $M_\mathrm{BCG}^{\star}$- $M_{200}$, $\sigma_\mathrm{BCG}^{\star}$- $\sigma_{200}$), we find in general a good agreement with observational results. Our simulations also predict $\sigma_\mathrm{BCG}^{\star}$- $\sigma_{200}$ relation to not change significantly up to redshift $z=1$, in line with a relatively slow accretion of the BCG stellar mass at late times. We analyze the main features of the velocity dispersion profiles, as traced by stars, dark matter, and galaxies. As a result, we discuss that observed stellar velocity dispersion profiles in the inner cluster regions are in excellent agreement with simulations. We also report that the slopes of the BCG velocity dispersion profile from simulations agree with what is measured in observations, confirming the existence of a robust correlation between the stellar velocity dispersion slope and the cluster velocity dispersion (thus, cluster mass) when the former is computed within $0.1 R_{500}$. Our results demonstrate that simulations can correctly describe the dynamics of BCGs and their surrounding stellar envelope, as determined by the past star-formation and assembly histories of the most massive galaxies of the Universe., Comment: 16 pages, 14 figures

The next generation of astronomical facilities will be able to retrieve exoplanetary atmospheric spectra in increasing quantity and of increasing quality. Radiative transfer (RT) models of these atmospheres is essential both for interpreting observational data and for linking these data to the planetary physical state with the aid of dedicated climate models. So far, a large effort has been placed in modelling the atmospheres of giant planets, which are the most easily accessible to observations. Now times are ripe to extend these studies to treat the relatively thin atmospheres of terrestrial-type exoplanets, which are the most likely targets for the search of atmospheric biomarkers. Here we present a procedure to perform radiative transfer calculations for terrestrial-type exoplanets with temperate physical conditions (Simonetti et al. in preparation). The procedure is based on HELIOS and HELIOS-K, which are novel, flexible and publicly available codes developed by the University of Bern (Grimm & Heng, 2015; Malik et al., 2017, 2019; Grimm et al., 2021) as a part of the Exoclime Simulation Platform (ESP) repository. These codes make full use of the computing power of Graphics Processing Units (GPUs, colloquially known as graphics cards) being therefore much faster (up to one order of magnitude, see Grimm et al. 2021) than other similar codes and are integrated with a variety of molecular and atomic line repositories such as HITRAN (Gordon et al., 2017), HITEMP (Rothman et al., 2010) and Kurucz. Until now, HELIOS has been mostly applied to study Jupiter-like planets. The main features of the procedure that we have implemented for the treatment of rocky, habitable planets can be summarized as follows. First, we added the treatment of the continuum absorption features of a variety of gases, in particular H2O (Clough et al., 1989; Mlawer et al., 2012) and CO2 (Gruszka & Borysow, 1997; Baranov et al., 2004; Baranov, 2018). These continua strongly contribute to the overall opacity of Earth-like atmospheres and cannot be neglected. Second, we paid special attention to the sub-Lorentzian profile of CO2 absorption lines, testing the effects of different recipes (Perrin & Hartmann, 1989; Tonkov et al., 1996). Third, we considered different hypotheses regarding the convective lapse rate of the troposphere. On these basis we: (i) tested the robustness of HELIOS and HELIOS-K against changes in model variables and (ii) compared them with other codes already published and used in the same context (e.g. LBLRTM Clough et al., 2005), as done by Yang et al. (2016). One of the main goals of this work is to provide a new and fast radiative transfer treatment for the ESTM, an energy balance climate model with upgraded treatment of the vertical and horizontal energy transport Vladilo et al. (2015). The ESTM is extremely flexible and allows for a rapid exploration of the planetary and atmospheric parameter space, providing us the ability to map quantitative indices of habitability on these parameters (Silva et al., 2017). The flexibility of both HELIOS and ESTM will allow us to test a wide variety of atmospheric compositions, which have applications in the study both of exoplanets and of ancient Earth and Mars. Moreover, the HELIOS procedure adapted to terrestrial-type atmospheres can be used to generate synthetic TOA fluxes useful to link the conditions at the planet's surface with quantities that will become observable with future generations of instruments, such as secondary eclipse spectra and direct IR emission spectra from terrestrial-type exoplanets (see e.g. Quanz et al., 2021). Finally, the output of the same procedure can be applied to other codes in the ESP repository, such as the THOR GCM (Mendonca et al., 2016; Deitrick et al., 2020). Figure 1 shows the TOA albedo obtained for three different stellar spectral classes for different values of surface temperature and stellar zenith angle, for an atmosphere of 1 bar of CO2 and a relative H2O humidity of 100%, as obtained by HELIOS using the procedure presented in Simonetti et al. (in preparation). The surface albedo was set to 0.15. Figure 2 shows the relation between OLR and surface temperature for different radiative transfer models for an Earth-like atmosphere composed by N2, O2, 360 ppmv of CO2, 1.8 ppmv of CH4 and a temperature-dependent quantity of H2O (relative humidity equal to 100%). The thick red curve labelled "HELIOS" has been obtained with the novel procedure presented in Simonetti et al. (in preparation). The data relative to the other curves have been taken from Yang et al. (2016). REFERENCES: Baranov, Y. I. 2018, Journal of Molecular Spectroscopy, 345, 11 Baranov, Y. I., La erty, W. J., & Fraser, G. T. 2004, Journal of Molecular Spectroscopy, 228, 432 Clough, S. A., Kneizys, F. X., & Davies, R. W. 1989, Atmospheric Research, 23, 229 Clough, S. A., Shephard, M. W., Mlawer, E. J., et al. 2005, JQSRT, 91, 233 Deitrick, R., Mendonca, J. M., Schro enegger, U., et al. 2020, ApJS, 248, 30 Gordon, I. E., Rothman, L. S., Hill, C., et al. 2017, JQSRT, 203, 3 Grimm, S. L., & Heng, K. 2015, ApJ, 808, 182 Grimm, S. L., Malik, M., Kitzmann, D., et al. 2021, ApJS, 253, 30 Gruszka, M., & Borysow, A. 1997, Icarus, 129, 172 Malik, M., Kitzmann, D., Mendonca, J. M., et al. 2019, AJ, 157, 170 Malik, M., Grosheintz, L., Mendonca, J. M., et al. 2017, AJ, 153, 56 Mendonca, J. M., Grimm, S. L., Grosheintz, L., & Heng, K. 2016, ApJ, 829, 115 Mlawer, E. J., Payne, V. H., Moncet, J. L., et al. 2012, Philosophical Transactions of the Royal Society of London Series A, 370, 2520 Perrin, M. Y., & Hartmann, J. M. 1989, JQSRT, 42, 311 Quanz, S. P., Ottiger, M., Fontanet, E., et al. 2021, arXiv e-prints, arXiv:2101.07500 Rothman, L. S., Gordon, I. E., Barber, R. J., et al. 2010, JQSRT, 111, 2139 Silva, L., Vladilo, G., Murante, G., & Provenzale, A. 2017, MNRAS, 470, 2270 Tonkov, M. V., Filippov, N. N., Bertsev, V. V., et al. 1996, Applied Optics, 35, 4863 Vladilo, G., Silva, L., Murante, G., Filippi, L., & Provenzale, A. 2015, ApJ, 804, 50 Yang, J., Leconte, J., Wolf, E. T., et al. 2016, ApJ, 826, 222

Using 324 numerically modelled galaxy clusters as provided by THE THREE HUNDRED project, we study the evolution of the kinematic properties of the stellar component of haloes on first infall. We select objects with M$_{\textrm{star}}>5\times10^{10} h^{-1}M_{\odot}$ within $3R_{200}$ of the main cluster halo at $z=0$ and follow their progenitors. We find that although haloes are stripped of their dark matter and gas after entering the main cluster halo, there is practically no change in their stellar kinematics. For the vast majority of our `galaxies' -- defined as the central stellar component found within the haloes that form our sample -- their kinematic properties, as described by the fraction of ordered rotation, and their position in the specific stellar angular momentum$-$stellar mass plane $j_{\rm star}$ -- M$_{\rm star}$, are mostly unchanged by the influence of the central host cluster. However, for a small number of infalling galaxies, stellar mergers and encounters with remnant stellar cores close to the centre of the main cluster, particularly during pericentre passage, are able to spin-up their stellar component by $z=0$., 8 pages, 5 figures; accepted for publication in A&A

The quest for inhabited worlds beyond the Solar System is focussed on rocky exoplanets, many of which are being discovered via transit and radial-velocity surveys. The detection of life in such planets via atmospheric biosignatures is challenging and ambient conditions that maximize the production and detectability of atmospheric biosignatures should be preferred in the selection of targets for spectroscopic observations. In this presentation we discuss how climate models that predict the temperature distribution on the planetary surface and the absorption properties of the planetary atmosphere can be used to narrow the search for exoplanets able to sustain multicellular organisms and hence, potentially, intelligent life.

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

We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range $0, Comment: 20 pages, 19 figures. Accepted for publication in A&A

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., Accepted for publication by MNRAS. Comments welcome

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

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)., 19 pages, 15 figures, accepted on June 8th, 2018, by MNRAS

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

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

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

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

In an effort to derive temperature based criteria of habitability for multicellular life, we investigated the thermal limits of terrestrial poikilotherms, i.e. organisms whose body temperature and the functioning of all vital processes is directly affected by the ambient temperature. Multicellular poikilotherms are the most common and evolutionarily ancient form of complex life on earth. The thermal limits for their active metabolism and reproduction are bracketed by the temperature interval 0C, Comment: 33 pages, 8 figures, abstract abridged, accepted for publication in the International Journal of Astrobiology

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

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., Accepted for publication in MNRAS. Comments welcome

Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics \texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations agrees well in trend and evolution with the observational results obtained from two data samples characterised by a wide range of masses and a large redshift coverage. We find that the iron abundance in the cluster core ($r, MNRAS, accepted

We measure and quantify properties of galactic outflows and diffuse gas at $z \geq 1$ in cosmological hydrodynamical simulations. Our novel sub-resolution model, MUPPI, implements supernova feedback using fully local gas properties, where the wind velocity and mass loading are not given as input. We find the following trends at $z = 2$ by analysing central galaxies having a stellar mass higher than $10^{9} M_{\odot}$. The outflow velocity and mass outflow rate ($\dot{M}_{\rm out}$) exhibit positive correlations with galaxy mass and with the star formation rate (SFR). However, most of the relations present a large scatter. The outflow mass loading factor ($\eta$) is between $0.2 - 10$. The comparison Effective model generates a constant outflow velocity, and a negative correlation of $\eta$ with halo mass. The number fraction of galaxies where outflow is detected decreases at lower redshifts, but remains more than $80 \%$ over $z = 1 - 5$. High SF activity at $z \sim 2 - 4$ drives strong outflows, causing the positive and steep correlations of velocity and $\dot{M}_{\rm out}$ with SFR. The outflow velocity correlation with SFR becomes flatter at $z = 1$, and $\eta$ displays a negative correlation with halo mass in massive galaxies. Our study demonstrates that both the MUPPI and Effective models produce significant outflows at $\sim 1 / 10$ of the virial radius; at the same time shows that the properties of outflows generated can be different from the input speed and mass loading in the Effective model. Our MUPPI model, using local properties of gas in the sub-resolution recipe, is able to develop galactic outflows whose properties correlate with global galaxy properties, and consistent with observations., Comment: 23 pages, 16 figures, accepted for publication in MNRAS

We present results of cosmological simulations of disc galaxies carried out with the GADGET-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator model. This description is based on a simple multiphase 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 supernovae 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 overproduction of stars. We use two sets of zoomed-in initial conditions of isolated cosmological haloes with masses (2-3) × 1012 M☉, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T ̃ 0.2), extended stellar and gas discs, 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 disc 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. at a 3σ level, thus indicating that our stellar discs 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.

Kepler-452b is currently the best example of an Earth-size planet in the habitable zone of a sun-like star, a type of planet whose number of detections is expected to increase in the future. Searching for biosignatures in the supposedly thin atmospheres of these planets is a challenging goal that requires a careful selection of the targets. Under the assumption of a rocky-dominated nature for Kepler-452b, we considered it as a test case to calculate a temperature-dependent habitability index, $h_{050}$, designed to maximize the potential presence of biosignature-producing activity (Silva et al.\ 2016). The surface temperature has been computed for a broad range of climate factors using a climate model designed for terrestrial-type exoplanets (Vladilo et al.\ 2015). After fixing the planetary data according to the experimental results (Jenkins et al.\ 2015), we changed the surface gravity, CO$_2$ abundance, surface pressure, orbital eccentricity, rotation period, axis obliquity and ocean fraction within the range of validity of our model. For most choices of parameters we find habitable solutions with $h_{050}>0.2$ only for CO$_2$ partial pressure $p_\mathrm{CO_2} \lesssim 0.04$\,bar. At this limiting value of CO$_2$ abundance the planet is still habitable if the total pressure is $p \lesssim 2$\,bar. In all cases the habitability drops for eccentricity $e \gtrsim 0.3$. Changes of rotation period and obliquity affect the habitability through their impact on the equator-pole temperature difference rather than on the mean global temperature. We calculated the variation of $h_{050}$ resulting from the luminosity evolution of the host star for a wide range of input parameters. Only a small combination of parameters yield habitability-weighted lifetimes $\gtrsim 2$\,Gyr, sufficiently long to develop atmospheric biosignatures still detectable at the present time., 15 pages, 6 figures, MNRAS accepted

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

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

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

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

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

We present predictions of Spectral Energy Distributions (SEDs), from the UV to the FIR, of simulated galaxies at $z=0$. These were obtained by post-processing the results of an N-body+hydro simulation of a small cosmological volume, that uses the Multi-Phase Particle Integrator (MUPPI) for star formation and stellar feedback, with the GRASIL-3D radiative transfer code, that includes reprocessing of UV light by dust. Physical properties of galaxies resemble observed ones, though with some tension at small and large stellar masses. Comparing predicted SEDs of simulated galaxies with different samples of local galaxies, we find that these resemble observed ones, when normalised at 3.6 $\mu$m. A comparison with the Herschel Reference Survey shows that, when binning galaxies in Star Formation Rate (SFR), average SEDs are reproduced to within a factor of $\sim2$ even in normalization, while binning in stellar mass highlights the same tension that is present in the stellar mass -- SFR plane. We use our sample to investigate the correlation of IR luminosity in Spitzer and Herschel bands with several galaxy properties. SFR is the quantity that best correlates with IR light up to $160\ \mu$m, while at longer wavelengths better correlations are found with molecular mass and, at $500\ \mu$m, with dust mass. However, using the position of the FIR peak as a proxy for cold dust temperature, we assess that heating of cold dust is mostly determined by SFR, with stellar mass giving only a minor contribution. We finally show how our sample of simulated galaxies can be used as a guide to understand the physical properties and selection biases of observed samples., Comment: 18 pages, 11 figures, Submitted to MNRAS

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

We implement novel numerical models of AGN feedback in the SPH code GADGET-3, where the energy from a supermassive black hole (BH) is coupled to the surrounding gas in the kinetic form. Gas particles lying inside a bi-conical volume around the BH are imparted a one-time velocity (10,000 km/s) increment. We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h M_sun), which is initially evolved to form a dense cool core, having central T, 22 pages, 11 figures, version accepted for publication in MNRAS, references and minor revisions added

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

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

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (, 24 pages, 13 figures, MNRAS submitted

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

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

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

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

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

We have investigated the chemo-dynamical evolution of a Milky Way-like disk galaxy, AqC4, produced by a cosmological simulation integrating a sub-resolution ISM model. We evidence a global inside-out and upside-down disk evolution, that is consistent with a scenario where the “thin disk” stars are formed from the accreted gas close to the galactic plane, while the older “thick disk” stars are originated in situ at higher heights. Also, the bar appears the most effective heating mechanism in the inner disk. Finally, no significant metallicity-rotation correlation has been observed, in spite of the presence of a negative [Fe/H] radial gradient.

The theory of the Milky Way formation, in the framework of the ΛCDM model, predicts galactic stellar halos to be built from multiple accretion events starting from the first structure to collapse in the Universe.Evidences in the past few decades have indicated that the Galactic halo consists of two overlapping structural components, an inner and an outer halo. We provide a set of numerical N-body simulations aimed to study the formation of the outer Milky Way (MW) stellar halo through accretion events between a (bulgeless) MW-like system and a satellite galaxy. After these minor mergers take place, in several orbital configurations, we analyze the signal left by satellite stars in the rotation velocity distribution. The aim is to explore the orbital conditions of the mergers where a signal of retrograde rotation in the outer part of the halo can be obtained, in order to give a possible explanation of the observed rotational properties of the MW stellar halo.Our results show that the dynamical friction has a fundamental role in assembling the final velocity distributions originated by different orbits and that retrograde satellites moving on low inclination orbits deposit more stars in the outer halo regions and therefore can produce the counter-rotating behavior observed in the outer MW halo.

We provide a set of numerical N-body simulations for studying the formation of the outer Milky Ways's stellar halo through accretion events. After simulating minor mergers of prograde and retrograde orbiting satellite halo with a Dark Matter main halo, we analyze the signal left by satellite stars in the rotation velocity distribution. The aim is to explore the orbital conditions where a retrograde signal in the outer part of the halo can be obtained, in order to give a possible explanation of the observed rotational properties of the Milky Way stellar halo. Our results show that, for satellites more massive than $\sim 1/40$ of the main halo, the dynamical friction has a fundamental role in assembling the final velocity distributions resulting from different orbits and that retrograde satellites moving on low inclination orbits deposit more stars in the outer halo regions end therefore can produce the counter-rotating behavior observed in the outer Milky Way halo., 5 pages, 3 figures, ApJL, accepted

Thick disks have been observed in many disk galaxies and our Galaxy, the Milky Way, also presents a thick disk whose main spatial, kinematic, and chemical features of this population are well established. However, the origin of this ancient component is still unclear in spite the many studies carried out and several formation scenarios proposed until now. For the first time to our knowledge, we found evidence of a kinematics-metallicity correlation, of about 40–50 km s-1 per dex, amongst thick disk stars at 1 kpc

This is the third of a series of papers presenting the first attempt to analyze the growth of the bar instability in a consistent cosmological scenario. In the previous two articles we explored the role of the cosmology on stellar disks, and the impact of the gaseous component on a disk embedded in a cosmological dark matter halo. The aim of this paper is to point out the impact of the star formation on the bar instability inside disks having different gas fractions. We perform cosmological simulations of the same disk-to-halo mass systems as in the previous works where the star formation was not triggered. We compare the results of the new simulations with the previous ones to investigate the effect of the star formation by analysing the morphology of the stellar components, the bar strength, the behaviour of the pattern speed. We follow the gas and the central mass concentration during the evolution and their impact on the bar strength. In all our cosmological simulations a stellar bar, lasting 10 Gyr, is still living at z=0. The central mass concentration of gas and of the new stars has a mild action on the ellipticity of the bar but is not able to destroy it; at z=0 the stellar bar strength is enhanced by the star formation. The bar pattern speed is decreasing with the disk evolution., 10 pages, 21 figures, A&A accepted

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]