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

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 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.

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

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

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.

[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 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

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]

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

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

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disk galaxies that are too compact with respect to the observations. High-resolution N-body/SPH simulations in a cosmological context are presented including cold gas and dark matter. A halo with quiet merging activity since z~3.8 and with a high spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disk. We show that the gas and the dark matter have similar specific angular momenta until a merger event occurs at z~2 with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and falls quickly into the center. Dynamical friction plays a minor role,in contrast to previous claims. In fact, after this event a new extended disk begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disk galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered star burst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected as a result of these processes, late-type galaxies could form with a dominant extended disk component, resulting from late infall, a small bulge-to-disk ratio and a low baryon fraction, in agreement with observations., Comment: 7 pages, 5 figures, submitted to MNRAS. Request for high resolution figures to the authors

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

X-Ray observations of groups and clusters of galaxies show that the Intra-Cluster Medium (ICM) in their cores is hotter than expected from cosmological numerical simulations of cluster formation which include star formation, radiative cooling and SN feedback. We investigate the effect of the injection of supersonic AGN jets into the ICM using axisymmetric hydrodynamical numerical simulations. A simple model for the ICM, describing the radial properties of gas and the gravitational potential in cosmological N-Body+SPH simulations of one cluster and three groups of galaxies at redshift z=0, is obtained and used to set the environment in which the jets are injected. We varied the kinetic power of the jet and the emission-weighted X-Ray temperature of the ICM. The jets transfer their energy to the ICM mainly by the effects of their terminal shocks. A high fraction of the injected energy can be deposited through irreversible processes in the cluster gas, up to 75% in our simulations. We show how one single, powerful jet can reconcile the predicted X-Ray properties of small groups, e.g. the Lx-Tx relation, with observations. We argue that the interaction between AGN jets and galaxy groups and cluster atmospheres is a viable feedback mechanism., 18 pages, 14 figures, to appear on Astronomy and Astrophysics

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

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

We evaluate the dependence of the cluster correlation length r_c on the mean intercluster separation D_c, for three models with critical matter density, vanishing vacuum energy (Lambda = 0) and COBE normalized: a tilted CDM (tCDM) model (n=0.8) and two blue mixed models with two light massive neutrinos yielding Omega_h = 0.26 and 0.14 (MDM1 and MDM2, respectively). All models approach the observational value of sigma_8 (and, henceforth, the observed cluster abundance) and are consistent with the observed abundance of Damped Lyman_alpha systems. Mixed models have a motivation in recent results of neutrino physics; they also agree with the observed value of the ratio sigma_8/sigma_25, yielding the spectral slope parameter Gamma, and nicely fit LCRS reconstructed spectra. We use parallel AP3M simulations, performed in a wide box (side 360/h Mpc) and with high mass and distance resolution, enabling us to build artificial samples of clusters, whose total number and mass range allow to cover the same D_c interval inspected through APM and Abell cluster clustering data. We find that the tCDM model performs substantially better than n=1 critical density CDM models. Our main finding, however, is that mixed models provide a surprisingly good fit of cluster clustering data., Comment: 22 pages + 10 Postscript figures. Accepted for publication in ApJ

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

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

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

The bar formation is still an open problem in modern astrophysics. In this paper we present numerical simulation performed with the aim of analyzing the growth of the bar instability inside stellar-gaseous disks, where the star formation is triggered, and a central black hole is present. The aim of this paper is to point out the impact of such a central massive black hole on the growth of the bar. We use N-body-SPH simulations of the same isolated disk-to-halo mass systems harboring black holes with different initial masses and different energy feedback on the surrounding gas. We compare the results of these simulations with the one of the same disk without black hole in its center. We make the same comparison (disk with and without black hole) for a stellar disk in a fully cosmological scenario. A stellar bar, lasting 10 Gyrs, is present in all our simulations. The central black hole mass has in general a mild effect on the ellipticity of the bar but it is never able to destroy it. The black holes grow in different way according their initial mass and their feedback efficiency, the final values of the velocity dispersions and of the black hole masses are near to the phenomenological constraints., Comment: 10 pages, 8 figures, accepted for pubblication in "Astrophysics and Space Science"

We present high resolution FLAMES/VLT spectroscopy of intracluster planetary nebula (ICPN) candidates, targeting three new fields in the Virgo cluster core with surface brightness down to mu_B = 28.5. Based on the projected phase space information we separate the old and 12 newly-confirmed PNs into galaxy and intracluster components. The M87 PNs are confined to the extended stellar envelope of M87, within a projected radius of ~ 160 kpc, while the ICPNs are scattered across the whole surveyed region between M87 and M86. The velocity dispersions determined from the M87 PNs at projected radii of 60 kpc and 144 kpc show that the galaxy's velocity dispersion profile decreases in the outer halo, down to 78 +/- 25 km/s. A Jeans model for the M87 halo stars in the gravitational potential traced by the X-ray emission fits the observed velocity dispersion profile only if the stellar orbits are strongly radially anisotropic (beta ~= 0.4 at r ~= 10 kpc increasing to 0.8 at the outer edge), and if additionally the stellar halo is truncated at ~= 150 kpc average elliptical radius. From the spatial and velocity distribution of the ICPNs we infer that M87 and M86 are falling towards each other and that we may be observing them just before the first close pass. The inferred luminosity-specific PN numbers for the M87 halo and the ICL are in the range of values observed for old (> 10 Gyr) stellar populations (abridged)., Accepted for publication in Astronomy and Astrophysics. 16 pages, 14 figures and 4 tables

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

We aim at finding a cosmologically motivated infall law to understand if the LambdaCDM cosmology can reproduce the main chemical characteristics of a Milky Way-like spiral galaxy. In this work we test several different gas infall laws, starting from that suggested in the two-infall model for the chemical evolution of the Milky Way by Chiappini et al., but focusing on laws derived from cosmological simulations which follows a concordance LambdaCDM cosmology. By means of a detailed chemical evolution model for the solar vicinity, we study the effects of the different gas infall laws on the abundance patterns and the G-dwarf metallicity distribution. The cosmological gas infall law predicts two main gas accretion episodes. By means of this cosmologically motivated infall law, we study the star formation rate, the SNIa and SNII rate, the total amount of gas and stars in the solar neighbourhood and the behaviour of several chemical abundances. We find that the results of the two-infall model are fully compatible with the evolution of the Milky Way with cosmological accretion laws. A gas assembly history derived from a DM halo, compatible with the formation of a late-type galaxy from the morphological point of view, can produce chemical properties in agreement with the available observations., Comment: This paper has 26 pages, 19 figures and 5 tables

Aims: We want to investigate the growth of bar instability in stellar disks embedded in a suitable dark matter halo evolving in a fully consistent cosmological framework. Methods: We perform seven cosmological simulations to emphasise the role of both the disk-to-halo mass ratio and of the Toomre parameter, Q, on the evolution of the disk.We also compare our fully cosmological cases with corresponding isolated simulations where the same halo, is extracted from the cosmological scenario and evolved in physical coordinates. Results: A long living bar, lasting about 10 Gyr, appears in all our simulations. In particular, disks expected to be stable according to classical criteria, form indeed weak bars. We argue that such a result is due to the dynamical properties of our cosmological halo which is far from stability and isotropy, typical of the classical halos used in literature; it is dynamically active, endowed of substructures and infall. Conclusions: At least for mild self-gravitating disks, the study of the bar instability using isolated isotropic halos, in gravitational equilibrium, can lead to misleading results. Furthermore, the cosmological framework is needed for quantitatively investigating such an instability., Comment: Astronomy & Astrophysics, accepted, 19 pages, 21 figures

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

We study, by numerical simulations, the propagation of an axisymmetric supersonic jet in an isothermal King atmosphere and we analyse the evolution of the resulting X-ray properties and their dependence on the jet physical parameters. We show the existence of two distinct regimes of interaction, with strong and weak shocks. In the first case shells of enhanced X-ray emission are to be expected, whereas in the second case we expect deficit of X-ray emission coincident with the cocoon. Analysing the results of our simulations we find that the jet kinetic power is the main parameter controlling the transition between the two regimes. We also discuss, in the same scheme, the ICM heating induced by the jet propagation, considering its effects on the observed relations between the cluster X-ray luminosity and temperature and between cluster entropy and temperature.

We investigate the bar instability in stellar disks embedded inside a dark matter halo evolving in a fully cosmological framework. We present cosmological simulations in which the role of the disk-halo mass ratio and the Q parameter are emphasized. Disks expected to be stable according to classical criteria form indeed weak bars. This result is due to dynamical properties of tha cosmological halo, which is far from stability and isotropy.

We present and analyse a linear, non‐resonant phenomenon of wave generation by vortices in smooth shear flows. The phenomenon is closely related to non‐normality of linear dynamics of perturbations in such flows and can be well interpreted by the use of the non‐modal approach, i.e. by following in time the linear dynamics of spatial Fourier harmonics of the vortex mode perturbations. We will consider both a planar two‐dimensional flow with linear velocity profile and a keplerian disk flow. We will also present numerical simulations in both configurations. These results may be important for the dynamics of accretion disks.

We show that dark matter haloes, in n--body simulations, have a boundary layer (BL) with precise features. In particular, it encloses all dynamically stable mass while, outside it, dynamical stability is lost soon. Particles can pass through such BL, which however acts as a confinement barrier for dynamical properties. BL is set by evaluating kinetic and potential energies (T(r) and W(r)) and calculating R=-2T/W. Then, on BL, R has a minimum which closely approaches a maximum of w= -dlog W/dlog r. Such $Rw$ ``requirement'' is consistent with virial equilibrium, but implies further regularities. We test the presence of a BL around haloes in spatially flat CDM simulations, with or without cosmological constant. We find that the mass M_c, enclosed within the radius r_c, where the $Rw$ requirement is fulfilled, closely approaches the mass M_{dyn}, evaluated from the velocities of all particles within r_c, according to the virial theorem. Using r_c we can then determine an individual density contrast Delta_c for each virialized halo, which can be compared with the "virial" density contrast $\Delta_v ~178 \Omega_m^{0.45}$ (Omega_m: matter density parameter) obtained assuming a spherically symmetric and unperturbed fluctuation growth. The spread in Delta_c is wide, and cannot be neglected when global physical quantities related to the clusters are calculated, while the average Delta_c is ~25 % smaller than the corresponding Delta_v; moreover if $M_{dyn}$ is defined from the radius linked to Delta_v, we have a much worse fit with particle mass then starting from {\it Rw} requirement., Comment: 4 pages, 5 figures, contribution to the XXXVIIth Rencontres de Moriond, The Cosmological Model, Les Arc March 16-23 2002, to appear in the proceedings

As a contribution to the study of the habitability of extrasolar planets, we implemented a 1-D Energy Balance Model (EBM), the simplest seasonal model of planetary climate, with new prescriptions for most physical quantities. Here we apply our EBM to investigate the surface habitability of planets with an Earth-like atmospheric composition but different levels of surface pressure. The habitability, defined as the mean fraction of the planet's surface on which liquid water could exist, is estimated from the pressure-dependent liquid water temperature range, taking into account seasonal and latitudinal variations of surface temperature. By running several thousands of EBM simulations we generated a map of the habitable zone (HZ) in the plane of the orbital semi-major axis, a, and surface pressure, p, for planets in circular orbits around a Sun-like star. As pressure increases, the HZ becomes broader, with an increase of 0.25 AU in its radial extent from p=1/3 bar to p=3 bar. At low pressure, the habitability is low and varies with a; at high pressure, the habitability is high and relatively constant inside the HZ. We interpret these results in terms of the pressure dependence of the greenhouse effect, the effciency of horizontal heat transport, and the extent of the liquid water temperature range. Within the limits discussed in the paper, the results can be extended to planets in eccentric orbits around non-solar type stars. The main characteristics of the pressure-dependent HZ are modestly affected by variations of planetary properties, particularly at high pressure., 58 pages, 12 figures, 3 tables, accepted for publication in ApJ. More figures can be found at http://adlibitum.oats.inaf.it/astrobiology/planhab/