Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $\sigma_\mathrm{v}$ with the density $\rho$ can be used to define the pseudo-entropy $S(r)=\sigma_\mathrm{v}^2/\rho^{2/3}$, whose profile is well-described by a simple power-law $S\propto r^{\alpha}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profile as traced by substructures, while pseudo-entropy profiles as traced by DM particles are characterized by an almost universal behavior and can be used as a low-scatter ($\sim6-7\%$) mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius.

We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range $0<z<4$, comparing them to recent observational studies; we also investigated the effect of varying the parameters of the star formation model on galaxy properties such as SFR, star-formation efficiency, and gas fraction. We analyze a set of zoom-in cosmological hydrodynamical simulations centred on twelve clusters. The simulations are carried out with the GADGET-3 TreePM/SPH code which includes various subgrid models to treat unresolved baryonic physics, including AGN feedback. Simulations do not reproduce the high values of SFR observed within protoclusters cores, where the values of SFR are underpredicted by a factor $\gtrsim 4$ both at $z\sim2$ and $z\sim 4$. The difference arises as simulations are unable to reproduce the observed starburst population and is worsened at $z\sim 2$ because simulations underpredict the normalization of the main sequence of star forming galaxies (i.e., the correlation between stellar mass and SFR) by a factor of $\sim 3$. As the low normalization of the main sequence seems to be driven by an underestimated gas fraction, it remains unclear whether numerical simulations miss starburst galaxies due to a too low predicted gas fractions or too low star formation efficiencies. Our results are stable against varying several parameters of the star formation subgrid model and do not depend on the details of the AGN feedback.

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.

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.

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.

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.

High performance computing numerical simulations are today one of the more effective instruments to implement and study new theoretical models, and they are mandatory during the preparatory phase and operational phase of any scientific experiment. New challenges in Cosmology and Astrophysics will require a large number of new extremely computationally intensive simulations to investigate physical processes at different scales. Moreover, the size and complexity of the new generation of observational facilities also implies a new generation of high performance data reduction and analysis tools pushing toward the use of Exascale computing capabilities. Exascale supercomputers cannot be produced today. We discuss the major technological challenges in the design, development and use of such computing capabilities and we will report on the progresses that has been made in the last years in Europe, in particular in the framework of the ExaNeSt European funded project. We also discuss the impact of this new computing resources on the numerical codes in Astronomy and Astrophysics.

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.

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<0.1R_{500}$) does not correlate with the temperature nor presents a significant evolution. The scale invariance is confirmed when the metallicity is related directly to the total mass. The slope of the best-fitting relations is shallow ($\beta\sim-0.1$) in the innermost regions ($r<0.5R_{500}$) and consistent with zero outside. We investigate the impact of the AGN feedback and find that it plays a key role in producing a constant value of the outskirts metallicity from groups to clusters. This finding additionally supports the picture of early enrichment.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

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.

We characterize the properties of the intergalactic medium (IGM) around a sample of galaxies extracted from state-of-the-art hydrodynamical simulations of structure formation in a cosmological volume of 25 Mpc comoving at $z\sim 2$. The simulations are based on two different subresolution schemes for star formation and supernova feedback: the MUlti-Phase Particle Integrator (MUPPI) scheme and the Effective Model. We develop a quantitative and probabilistic analysis based on the apparent optical depth method of the properties of the absorbers as a function of impact parameter from their nearby galaxies: in such a way we probe different environments from circumgalactic medium to low-density filaments. Absorbers' properties are then compared with a spectroscopic observational data set obtained from high-resolution quasar spectra. Our main focus is on the N$_{\rm{ CIV}}$-N$_{\rm { HI}}$ relation around simulated galaxies: the results obtained with MUPPI and the Effective model are remarkably similar, with small differences only confined to regions at impact parameters $b = [1-3] \times r_{\rm {vir}}$. Using $\mbox{C IV}$ as a tracer of the metallicity, we obtain evidence that the observed metal absorption systems have the highest probability to be confined in a region of 150-400 kpc around galaxies. Near-filament environments have instead metallicities too low to be probed by present-day telescopes, but could be probed by future spectroscopical studies. Finally we compute $\mbox{C IV}$ covering fractions which are in agreement with observational data.

We analyze the stellar growth of Brightest Cluster Galaxies (BCGs) produced by cosmological zoom-in hydrodynamical simulations of the formation of massive galaxy clusters. The evolution of the stellar mass content is studied considering different apertures, and tracking backwards either the main progenitor of the $z=0$ BCG or that of the cluster hosting the BCG at $z=0$. Both methods lead to similar results up to $z \simeq 1.5$. The simulated BCGs masses at $z=0$ are in agreement with recent observations. In the redshift interval from $z=1$ to $z=0$ we find growth factors 1.3, 1.6 and 3.6 for stellar masses within 30kpc, 50kpc and 10% of $R_{500}$ respectively. The first two factors, and in general the mass evolution in this redshift range, are in agreement with most recent observations. The last larger factor is similar to the growth factor obtained by a semi-analytical model (SAM). Half of the star particles that end up in the inner 50 kpc was typically formed by redshift $\sim$ 3.7, while the assembly of half of the BCGs stellar mass occurs on average at lower redshifts $\sim 1.5$. This assembly redshift correlates with the mass attained by the cluster at high $z \gtrsim 1.3$, due to the broader range of the progenitor clusters at high-$z$. The assembly redshift of BCGs decreases with increasing apertures. Our results are compatible with the {\it inside-out} scenario. Simulated BCGs could lack intense enough star formation (SF) at high redshift, while possibly exhibit an excess of residual SF at low redshift.

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

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.

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.

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

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

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.

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.

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.

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.

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

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.

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<10^6 K. A BH resides at the cluster center, and ejects energy. The feedback-driven fast wind undergoes shock with the slower-moving gas, which causes the imparted kinetic energy to be thermalized. Bipolar bubble-like outflows form propagating radially outward to a distance of a few 100 kpc. The radial profiles of median gas properties are influenced by BH feedback in the inner regions (r<20-50 kpc). BH kinetic feedback, with a large value of the feedback efficiency, depletes the inner cool gas and reduces the hot gas content, such that the initial cool core of the cluster is heated up within a time 1.9 Gyr, whereby the core median temperature rises to above 10^7 K, and the central entropy flattens. Our implementation of BH thermal feedback (using the same efficiency as kinetic), within the star-formation model, cannot do this heating, where the cool core remains. The inclusion of cold gas accretion in the simulations produces naturally a duty cycle of the AGN with a periodicity of 100 Myr.

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.

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<T<50C. The same interval applies to the photosynthetic production of oxygen, an essential ingredient of complex life, and for the generation of atmospheric biosignatures. Analysis of the main mechanisms responsible for the thermal thresholds of terrestrial life suggests that the same mechanisms would apply to other forms of chemical life. We propose a habitability index for complex life, h050, representing the mean orbital fraction of planetary surface that satisfies the temperature limits 0C<T<50C. With the aid of a climate model tailored for the calculation of the surface temperature of Earth-like planets, we calculated h050 as a function of planet insolation S, and atmospheric columnar mass Natm, for a few earth-like atmospheric compositions. By displaying h050 as a function of S and Natm, we built up an atmospheric mass habitable zone (AMHZ) for complex life. At variance with the classic habitable zone, the inner edge of the complex life HZ is not affected by the uncertainties inherent to the calculation of the runaway greenhouse limit. The complex life HZ is significantly narrower than the HZ of dry planets. Our calculations illustrate how changes in ambient conditions dependent on S and Natm, such as temperature excursions and surface dose of secondary particles of cosmic rays, may influence the type of life potentially present at different epochs of planetary evolution inside the AMHZ.

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 (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

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.

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.

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.

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

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.

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.

We introduce a novel Earth-like planet surface temperature model (ESTM) for habitability studies based on the spatial-temporal distribution of planetary surface temperatures. The ESTM adopts a surface Energy Balance Model complemented by: radiative-convective atmospheric column calculations, a set of physically-based parameterizations of meridional transport, and descriptions of surface and cloud properties more refined than in standard EBMs. The parameterization is valid for rotating terrestrial planets with shallow atmospheres and moderate values of axis obliquity (epsilon >= 45^o). Comparison with a 3D model of atmospheric dynamics from the literature shows that the equator-to-pole temperature differences predicted by the two models agree within ~5K when the rotation rate, insolation, surface pressure and planet radius are varied in the intervals 0.5 <= Omega/Omega_o <= 2, 0.75 <= S/S_o <= 1.25, 0.3 <= p/(1 bar) <= 10, and 0.5 <= R/R_o <= 2, respectively. The ESTM has an extremely low computational cost and can be used when the planetary parameters are scarcely known (as for most exoplanets) and/or whenever many runs for different parameter configurations are needed. Model simulations of a test-case exoplanet (Kepler-62e) indicate that an uncertainty in surface pressure within the range expected for terrestrial planets may impact the mean temperature by ~60 K. Within the limits of validity of the ESTM, the impact of surface pressure is larger than that predicted by uncertainties in rotation rate, axis obliquity, and ocean fractions. We discuss the possibility of performing a statistical ranking of planetary habitability taking advantage of the flexibility of the ESTM.

We compute and study the infrared and sub-mm properties of high redshift ($z \gtrsim 1$) simulated clusters and proto-clusters. The results of a large set of hydro-dynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry. Here, we have slightly generalized the code to adapt it to the present purpose. Then we have post-processed boxes of physical size 2 Mpc encompassing each of the 24 most massive clusters identified at z=0, at several redshifts between 0.5 and 3, producing IR and sub-mm mock images of these regions and SEDs of the radiation coming out from them. While this field is in its infancy from the observational point of view, rapid development is expected in the near future thanks to observations performed in the far IR and sub-mm bands. Notably, we find that in this spectral regime our prediction are little affected by the assumption required by this post-processing, and the emission is mostly powered by star formation rather than accretion onto super massive black hole (SMBH). The comparison with the little observational information currently available, highlights that the simulated cluster regions never attain the impressive star formation rates suggested by these observations. This problem becomes more intriguing taking into account that the brightest cluster galaxies (BCGs) in the same simulations turn out to be too massive. It seems that the interplay between the feedback schemes and the star formation model should be revised, possibly incorporating a positive feedback mode.

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

We present a quantification of the properties of bars in two N-body+SPH cosmological simulations of spiral galaxies, named GA and AqC. The initial conditions were obtained using the zoom-in technique and represent two dark matter (DM) halos of $2-3\times10^{12}\ {\rm M}_\odot$, available at two different resolutions. The resulting galaxies are presented in the companion paper of Murante et al. (2014). We find that the GA galaxy has a bar of length $8.8$ kpc, present at the two resolution levels even though with a slightly different strength. Classical bar signatures (e.g. pattern of streaming motions, high $m=2$ Fourier mode with roughly constant phase) are consistently found at both resolutions. Though a close encounter with a merging satellite at $z\sim0.6$ (mass ratio $1:50$) causes a strong, transient spiral pattern and some heating of the disk, we find that bar instability is due to secular process, caused by a low Toomre parameter $Q\lesssim1$ due to accumulation of mass in the disk. The AqC galaxy has a slightly different history: it suffers a similar tidal disturbance due to a merging satellite at $z\sim0.5$ but with a mass ratio of $1:32$, that triggers a bar in the high-resolution simulation, while at low resolution the merging is found to take place at a later time, so that both secular evolution and merging are plausible triggers for bar instability.

We present results of cosmological simulations of disk galaxies carried out with the GADGET-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator (MUPPI) model. This description is based on simple multi-phase model of the interstellar medium at unresolved scales, where mass and energy flows among the components are explicitly followed by solving a system of ordinary differential equations. Thermal energy from SNe is injected into the local hot phase, so as to avoid that it is promptly radiated away. A kinetic feedback prescription generates the massive outflows needed to avoid the over-production of stars. We use two sets of zoomed-in initial conditions of isolated cosmological halos with masses (2-3) * 10^{12} Msun, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T \approx 0.2), extended stellar and gas disks, flat rotation curves and realistic values of stellar masses. Gas profiles are relatively flat, molecular gas is found to dominate at the centre of galaxies, with star formation rates following the observed Schmidt-Kennicutt relation. Stars kinematically belonging to the bulge form early, while disk stars show a clear inside-out formation pattern and mostly form after redshift z=2. However, the baryon conversion efficiencies in our simulations differ from the relation given by Moster et al. (2010) at a 3 sigma level, thus indicating that our stellar disks are still too massive for the Dark Matter halo in which they reside. Results are found to be remarkably stable against resolution. This further demonstrates the feasibility of carrying out simulations producing a realistic population of galaxies within representative cosmological volumes, at a relatively modest resolution.

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.

[Abridged.] We investigate baryon effects on the halo mass function (HMF), with emphasis on the role played by AGN feedback. Halos are identified with both Friends-of-Friends (FoF) and Spherical Overdensity (SO) algorithms. We embed the standard SO algorithm into a memory-controlled frame program and present the {\bf P}ython spher{\bf I}c{\bf A}l {\bf O}verdensity code --- {\small PIAO}. For both FoF and SO halos, the effect of AGN feedback is that of suppressing the HMFs to a level even below that of Dark Matter simulations. The ratio between the HMFs in the AGN and in the DM simulations is $\sim 0.8$ at overdensity $\Delta_c=500$, a difference that increases at higher overdensity $\Delta_c=2500$, with no significant redshift and mass dependence. A decrease of the halo masses ratio with respect to the DM case induces the decrease of the HMF in the AGN simulation. The shallower inner density profiles of halos in the AGN simulation witnesses that mass reduction is induced by the sudden displacement of gas induced by thermal AGN feedback. We provide fitting functions to describe halo mass variations at different overdensities, which can recover the HMFs with a residual random scatter $\lt 5$ per cent for halo masses larger than $10^{13} ~h^{-1}{\rm M_\odot}$.

Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10% to 30%. The largest bias originates by both violation of hydrostatic equilibrium and an additional temperature bias caused by inhomogeneities in the X-ray emitting intra-cluster medium (ICM). To elucidate on this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle-hydrodynamics (SPH) and adaptive-mesh-refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in no-radiative simulations, while it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with the observed one in the central regions of clusters. Around R_500 the temperature inhomogeneities of the SPH simulations can generate twice the typical hydrostatic-equilibrium mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.

We study the evolution of disc galaxies in group environments under the effect of both the global tidal field and close-encounters between galaxies, using controlled N-body simulations of isolated mergers. We find that close-range encounters between galaxies are less frequent and less damaging to disc galaxies than originally expected, since they mostly occur when group members have lost a significant fraction of their initial mass to tidal stripping. We also find that group members mostly affect disc galaxies indirectly by modifying their common global tidal field. Different initial orbital parameters of group members introduce a significant "scatter" in the evolution of general properties of disc galaxies around a "median" evolution that is similar to when only the effect of the global tidal field is included. Close-encounters introduce a high variability in the properties of disc galaxies, even slowing their evolution in some cases, and could wash out correlations between galaxy properties and the group total mass. The combined effect of the global tidal field and close-encounters appears to be inefficient at forming/enhancing central stellar bulges. This implies that bulges of S0 galaxies should be mostly composed by young stars, which is consistent with recent observations.

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

We present an analysis of the properties of the ICM in an extended set of cosmological hydrodynamical simulations of galaxy clusters and groups performed with the TreePM+SPH GADGET-3 code. Besides a set of non-radiative simulations, we carried out two sets of simulations including radiative cooling, star formation, metal enrichment and feedback from supernovae, one of which also accounts for the effect of feedback from AGN resulting from gas accretion onto super-massive black holes. These simulations are analysed with the aim of studying the relative role played by SN and AGN feedback on the general properties of the diffuse hot baryons in galaxy clusters and groups: scaling relations, temperature, entropy and pressure radial profiles, and ICM chemical enrichment. We find that simulations including AGN feedback produce scaling relations that are in good agreement with X-ray observations at all mass scales. However, our simulations are not able to account for the observed diversity between CC and NCC clusters: unlike for observations, we find that temperature and entropy profiles of relaxed and unrelaxed clusters are quite similar and resemble more the observed behaviour of NCC clusters. As for the pattern of metal enrichment, we find that an enhanced level of iron abundance is produced by AGN feedback with respect to the case of purely SN feedback. As a result, while simulations including AGN produce values of iron abundance in groups in agreement with observations, they over-enrich the ICM in massive clusters. The efficiency of AGN feedback in displacing enriched gas from halos into the inter-galactic medium at high redshift also creates a widespread enrichment in the outskirts of clusters and produces profiles of iron abundance whose slope is in better agreement with observations.

(Abridged) We investigate two modes of coupling the feedback energy from a central AGN to the neighboring gas in galaxy simulations: kinetic - velocity boost, and thermal - heating. We formulate kinetic feedback models for energy-driven wind (EDW) and momentum-driven wind (MDW), using two free parameters: feedback efficiency epsilon_f, and AGN wind velocity v_w. A novel numerical algorithm is implemented in the SPH code GADGET-3, to prevent the expansion of a hole in the gas distribution around the BH. We perform simulations of isolated evolution and merger of disk galaxies, of Milky-Way mass as well as lower and higher masses. We find that in the isolated galaxy BH kinetic feedback generates intermittent bipolar jet-like gas outflows. We infer that current prescriptions for BH subgrid physics in galaxy simulations can grow the BH to observed values even in an isolated disk galaxy. The BH growth is enhanced in a galaxy merger. Comparing the [M_BH - sigma_star] relation obtained in our simulations with observational data, we conclude that it is possible to find parameter sets for a fit in all the models, except for the case with MDW feedback in a galaxy merger, in which the BH is always too massive. The BH thermal feedback implementation of Springel, Di Matteo & Hernquist (2005) within the multiphase star-formation model is found to have negligible impact on gas properties; and the effect claimed in all previous studies is attributed to gas depletion around the BH by the creation of an artificial hole. The BH mass accretion rate in our simulations exhibit heavy fluctuations. The star formation rate is quenched with feedback by removal of gas. The CGM gas at galactocentric distances (20 - 100)/h kpc are found to give the best metallicity observational diagnostic to distinguish between BH models.

We analyze the basic properties of Brightest Cluster Galaxies (BCGs) produced by state of the art cosmological zoom-in hydrodynamical simulations. These simulations have been run with different sub-grid physics included. Here we focus on the results obtained with and without the inclusion of the prescriptions for supermassive black hole (SMBH) growth and of the ensuing Active Galactic Nuclei (AGN) feedback. The latter process goes in the right direction of decreasing significantly the overall formation of stars. However, BCGs end up still containing too much stellar mass, a problem that increases with halo mass, and having an unsatisfactory structure. This is in the sense that their effective radii are too large, and that their density profiles feature a flattening on scales much larger than observed. We also find that our model of thermal AGN feedback has very little effect on the stellar velocity dispersions, which turn out to be very large. Taken together, these problems, which to some extent can be recognized also in other numerical studies typically dealing with smaller halo masses, indicate that on one hand present day sub-resolution models of AGN feedback are not effective enough in diminishing the global formation of stars in the most massive galaxies, but on the other hand they are relatively too effective in their centers. It is likely that a form of feedback generating large scale gas outflows from BCGs precursors, and a more widespread effect over the galaxy volume, can alleviate these difficulties.

We compare three analytical prescriptions for merger times available from the literature to simulations of isolated mergers. We probe three different redshifts, and several halo concentrations, mass ratios, orbital circularities and orbital energies of the satellite. We find that prescriptions available in the literature significantly under-predict long timescales for mergers at high redshift. We argue that these results have not been highlighted previously either because the evolution of halo concentration of satellite galaxies has been neglected (in previous isolated merger simulations), or because long merger times and mergers with high initial orbital circularities are under-represented (for prescriptions based on cosmological simulations). Motivated by the evolution of halo concentration at fixed mass, an explicit dependence on redshift added as t_merger,modified(z) = (1+z)^0.44 t_merger to the prescription based on isolated mergers gives a significant improvement in the predicted merger times up to ~20 t_dyn in the redshift range 0<z<2. When this modified prescription is used to compute galaxy stellar mass functions, we find that it leads up to a 25 per cent increase in the number of low mass galaxies surviving at z=0, and a 10 per cent increase for more massive galaxies. This worsen the known over-prediction in the number of low mass galaxies by hierarchical models of galaxy formation.

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.

We carry out an analysis of a set of cosmological SPH hydrodynamical simulations of galaxy clusters and groups aimed at studying the total baryon budget in clusters, and how this budget is shared between the hot diffuse component and the stellar component. Using the TreePM+SPH GADGET-3 code, we carried out one set of non-radiative simulations, and two sets of simulations including radiative cooling, star formation and feedback from supernovae (SN), one of which also accounting for the effect of feedback from active galactic nuclei (AGN). The analysis is carried out with the twofold aim of studying the implication of stellar and hot gas content on the relative role played by SN and AGN feedback, and to calibrate the cluster baryon fraction and its evolution as a cosmological tool. We find that both radiative simulation sets predict a trend of stellar mass fraction with cluster mass that tends to be weaker than the observed one. However this tension depends on the particular set of observational data considered. Including the effect of AGN feedback alleviates this tension on the stellar mass and predicts values of the hot gas mass fraction and total baryon fraction to be in closer agreement with observational results. We further compute the ratio between the cluster baryon content and the cosmic baryon fraction, Y_b, as a function of cluster-centric radius and redshift. At R_500 we find for massive clusters with M_500>2\times10^{14} h^{-1} M_sun that Y_b is nearly independent of the physical processes included and characterized by a negligible redshift evolution: Y_{b,500}=0.85+/-0.03 with the error accounting for the intrinsic r.m.s. scatter within the set of simulated clusters. At smaller radii, R_2500, the typical value of Y_b slightly decreases, by an amount that depends on the physics included in the simulations, while its scatter increases by about a factor of two.

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

The thick disk rotation--metallicity correlation, \partial V_\phi/\partial[Fe/H] =40\div 50 km s^{-1}dex^{-1} represents an important signature of the formation processes of the galactic disk. We use nondissipative numerical simulations to follow the evolution of a Milky Way (MW)-like disk to verify if secular dynamical processes can account for this correlation in the old thick disk stellar population. We followed the evolution of an ancient disk population represented by 10 million particles whose chemical abundances were assigned by assuming a cosmologically plausible radial metallicity gradient with lower metallicity in the inner regions, as expected for the 10-Gyr-old MW. Essentially, inner disk stars move towards the outer regions and populate layers located at higher |z|. A rotation--metallicity correlation appears, which well resembles the behaviour observed in our Galaxy at a galactocentric distance between 8 kpc and 10 kpc. In particular,we measure a correlation of \partial V_\phi/\partial[Fe/H]\simeq 60 km s^{-1}dex^{-1} for particles at 1.5 kpc < |z| < 2.0 kpc that persists up to 6 Gyr. Our pure N-body models can account for the V_\phi vs. [Fe/H] correlation observed in the thick disk of our Galaxy, suggesting that processes internal to the disk such as heating and radial migration play a role in the formation of this old stellar component. In this scenario, the positive rotation-metallicity correlation of the old thick disk population would represent the relic signature of an ancient "inverse" chemical (radial) gradient in the inner Galaxy, which resulted from accretion of primordial gas.

Forthcoming experiments will enable us to determine tomographic shear spectra at a high precision level. Most predictions about them have until now been biased on algorithms yielding the expected linear and non-linear spectrum of density fluctuations. Even when simulations have been used, so-called Halofit (Smith et al 2003) predictions on fairly large scales have been needed. We wish to go beyond this limitation. We perform N-body and hydrodynamical simulations within a sufficiently large cosmological volume to allow a direct connection between simulations and linear spectra. While covering large length-scales, the simulation resolution is good enough to allow us to explore the high-l harmonics of the cosmic shear (up to l ~ 50000), well into the domain where baryon physics becomes important. We then compare shear spectra in the absence and in presence of various kinds of baryon physics, such as radiative cooling, star formation, and supernova feedback in the form of galactic winds. We distinguish several typical properties of matter fluctuation spectra in the different simulations and test their impact on shear spectra. We compare our outputs with those obtainable using approximate expressions for non--linear spectra, and identify substantial discrepancies even between our results and those of purely N-body results. Our simulations and the treatment of their outputs however enable us, for the first time, to obtain shear results taht are fully independent of any approximate expression, also in the high-l range, where we need to incorporate a non-linear power spectrum of density perturbations, and the effects of baryon physics. This will allow us to fully exploit the cosmological information contained in future high--sensitivity cosmic shear surveys, exploring the physics of cosmic shears via weak lensing measurements.

Forthcoming experiments will enable us to determine high precision tomographic shear spectra. Matter density fluctuation spectra, at various $z$, should then be worked out of them, in order to constrain the model and determine the DE state equation. Available analytical expressions, however, do the opposite, enabling us to derive shear spectra from fluctuation spectra. Here we find the inverse expression, yielding density fluctuation spectra from observational tomographic shear spectra. The procedure involves SVD techniques for matrix inversion. We show in detail how the approach works and provide a few examples.

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 present results from high--resolution cosmological hydrodynamical simulations of a Milky--Way-sized halo, aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter-stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial radius. Gas particles can be accreted through three different channels, depending on the maximum temperature value, $T_{\rm max}$, reached during the particles' past evolution: a cold channel for $T_{\rm max}<2.5 \times 10^5$ K, a hot one for $T>10^6$K, and a warm one for intermediate values of $T_{\rm max}$. We find that the warm channel is at least as important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous findings that the cold mode dominates gas accretion at high redshift. We ascribe this difference to the different supernova feedback scheme implemented in our simulations. While results presented so far in the literature are based on uneffective SN thermal feedback schemes and/or the presence of a kinetic feedback, our simulations include only effective thermal feedback. We argue that observational detections of a warm accretion mode in the high--redshift circum-galactic medium would provide useful constraints on the nature of the feedback that regulates star formation in galaxies.

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 multi-component, collision-less, N-body systems, composed by both dark matter and stars. In our simulations, the evolution of disc galaxies is followed 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. That is, discs on face-on and retrograde orbits are shown to retain longer their disc structures and kinematics, 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.

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]

Magnetic fields in the intra-cluster medium of galaxy clusters have been studied in the past years through different methods. In the next years the up-coming generation of radio telescopes is going to provide new data that have the potential of setting constraints on the properties of magnetic fields in galaxy clusters. Here we present zoomed-in simulations for a set of massive galaxy clusters (M_v > 10^15 M_sun/h). This is an ideal sample to study the evolution of magnetic field during the process of structure formation in detail. Turbulent motions of the gas within the ICM will manifest themselves in a macroscopic magnetic resistivity eta_m, which has to be taken explicitly into account, especially at scales below the resolution limit. We have adapted the MHD GADGET code by Dolag & Stasyszyn (2009) to include the treatment of the magnetic resistivity and for the first time we have included non-ideal MHD equations to better follow the evolution of the magnetic field within galaxy clusters. We investigate which value of the magnetic resistivity eta_m is required to match the magnetic field profile derived from radio observations. We find that a value of eta_m ~ 6 10^27 cm^2/s is necessary to recover the shape of the magnetic field profile inferred from radio observations of the Coma cluster. This value agrees well with the expected level of turbulent motions within the ICM at our resolution limit. The magnetic field profiles can be fitted by a beta-model like profile (Cavaliere & Fusco-Femiano 1976), with small dispersion of the parameters. We find also that that the temperature, density and entropy profiles of the clusters depend on the magnetic resistivity constant, having flatter profiles in the inner regions when the magnetic resistivity increases.

We present results based on an implementation of the Godunov Smoothed Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of the equations of moment and energy conservation, starting from the convolution of these equations with the interpolating kernel. The two most important aspects of the numerical implementation of these equations are (a) the appearance of fluid velocity and pressure obtained from the solution of the Riemann problem between each pair of particles, and (b the absence of an artificial viscosity term. We carry out three different controlled hydrodynamical three-dimensional tests, namely the Sod shock tube, the development of Kelvin-Helmholtz instabilities in a shear flow test, and the "blob" test describing the evolution of a cold cloud moving against a hot wind. The results of our tests confirm and extend in a number of aspects those recently obtained by Cha (2010): (i) GSPH provides a much improved description of contact discontinuities, with respect to SPH, thus avoiding the appearance of spurious pressure forces; (ii) GSPH is able to follow the development of gas-dynamical instabilities, such as the Kevin--Helmholtz and the Rayleigh-Taylor ones; (iii) as a result, GSPH describes the development of curl structures in the shear-flow test and the dissolution of the cold cloud in the "blob" test. We also discuss in detail the effect on the performances of GSPH of changing different aspects of its implementation. The results of our tests demonstrate that GSPH is in fact a highly promising hydrodynamic scheme, also to be coupled to an N-body solver, for astrophysical and cosmological applications. [abridged]

We present a detailed study of scaling relations between total cluster mass and three mass proxies based on X-ray observables: temperature of the intra-cluster medium, gas mass and the product of the two, Y_X. Our analysis is based on two sets of high-resolution hydrodynamical simulations performed with the TreePM-SPH GADGET code. The first set includes about 140 clusters with masses above 5x10^13 M_sun/h (30 having mass above 10^15 M_sun/h), that have been simulated with (i) non-radiative physics and including (ii) cooling, star formation, chemical enrichment and the effect of supernova feedback triggering galactic ejecta. This large statistics is used to quantify the robustness of the scaling relations, to determine their redshift evolution and to calibrate their intrinsic scatter and its distribution. We use a smaller set of clusters including 18 halos with masses above 5x10^13 M_sun/h to test the robustness of mass proxies against changing the physical processes included in simulations (thermal conduction, artificial viscosity, cooling and star formation, galactic winds and AGN feedback). We find the M-Y_X scaling relation to be the least sensitive one to variations of the ICM physics, with its slope and redshift evolution close to the self-similar model predictions. The distribution of the scatter around the best-fitting relations is close to a log-normal one. M_gas has the smallest scatter in mass, with values of sigma_lnM = 0.04-0.06, depending on the physics included in the simulation, and with a mild dependence on redshift. The M-T relation is the one with the largest scatter, with sigma_lnM > 0.1 at z=0, increasing to > 0.15 at z=1. The intrinsic scatter in the M-Y_X relation is slightly larger than in the M-M_gas relation. These results confirm that both Y_X and M_gas mass proxies are well suited for cosmological applications of future large X-ray surveys. [abridged]

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.

When clusters of galaxies are viewed in projection, one cannot avoid picking up foreground/background interlopers (FBIs), that lie within the virial cone (VC), but outside the virial sphere. Structural & kinematic deprojection equations are not known for an expanding Universe, where the Hubble flow (HF) stretches the line-of-sight (LOS) distribution of velocities. We analyze 93 mock relaxed clusters, built from a cosmological simulation. The stacked mock cluster is well fit by an m=5 Einasto DM density profile (but only out to 1.5 virial radii [r_v]), with velocity anisotropy (VA) close to the Mamon-Lokas model with VA radius equal to that of density slope -2. The surface density of FBIs is nearly flat out to r_v, while their LOS velocity distribution shows a dominant gaussian cluster-outskirts component and a flat field component. This distribution of FBIs in projected phase space is nearly universal in mass. A local k=2.7 sigma velocity cut returns the LOS velocity dispersion profile (LOSVDP) expected from the NFW density and VA profiles measured in 3D. The HF causes a shallower outer LOSVDP that cannot be well matched by the Einasto model for any k. After this velocity cut, FBIs still account for 23% of DM particles within the VC (close to the observed fraction of cluster galaxies lying off the Red Sequence). The best-fit projected NFW/Einasto models underestimate the 3D concentration by 6+/-6% (16+/-7%) after (before) the velocity cut, unless a constant background is included in the fit. Assuming the correct mass profile, the VA profile is well recovered from the measured LOSVDP, with a slight bias towards more radial orbits in the outer regions. These small biases are overshadowed by large cluster-cluster variations caused by cosmic variance. An appendix provides an analytical approximation to the surface density, projected mass and tangential shear profiles of the Einasto model.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

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.

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]

Member galaxies within galaxy clusters nowadays can be routinely identified in cosmological, hydrodynamical simulations using methods based on identifying self bound, locally over dense substructures. However, distinguishing the central galaxy from the stellar diffuse component within clusters is notoriously difficult, and in the center it is not even clear if two distinct stellar populations exist. Here, after subtracting all member galaxies, we use the velocity distribution of the remaining stars and detect two dynamically, well-distinct stellar components within simulated galaxy clusters. These differences in the dynamics can be used to apply an un-binding procedure which leads to a spatial separation of the two components into a cD and a diffuse stellar component (DSC). Applying our new algorithm to a cosmological, hydrodynamical simulation we find that -- in line with previous studies -- these two components have clearly distinguished spatial and velocity distributions as well as different star formation histories. We show that the DSC fraction -- which can broadly be associated with the observed intra cluster light -- does not depend on the virial mass of the galaxy cluster and is much more sensitive to the formation history of the cluster. We conclude that the separation of the cD and the DSC in simulations, based on our dynamical criteria, is more physically motivated than current methods which depend on implicit assumptions on a length scale associated with the cD galaxy and therefore represent a step forward in understanding the different stellar components within galaxy clusters. Our results also show the importance of analyzing the dynamics of the DSC to characterize its properties and understand its origin.

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.

We present a study of the effect of AGN feedback on metal enrichment and thermal properties of the intracluster medium (ICM) in hydrodynamical simulations. The cosmological simulations are performed for a set of clusters using a version of the TreePM-SPH Gadget code that follows chemo-dynamical evolution by accounting for metal enrichment by different stellar populations. Besides runs not including any efficient form of energy feedback, we carry out simulations including: (i) kinetic feedback in the form of galactic winds triggered by supernova explosions; (ii) AGN feedback from gas accretion onto super-massive black holes (BHs); (iii) AGN feedback in which a 'radio mode' is included. We find that AGN feedback is able to quench star formation in the brightest cluster galaxies at z<4 and provides correct temperature profiles in the central regions of galaxy groups. However, its effect is not sufficient to create cool cores in massive clusters. AGN feedback creates a widespread enrichment in the outskirts of clusters, thanks to its efficiency in displacing enriched gas from galactic halos to the inter-galactic medium at relatively high redshift. Iron abundance profiles are in better agreement with observations, with a more pristine enrichment of the ICM around and beyond the cluster virial regions. From the pattern of the relative abundances of Silicon and Iron, we conclude that a significant fraction of ICM enrichment in simulations is contributed by a diffuse population of intra-cluster stars. Our simulations also predict that profiles of Z_Si/Z_Fe abundance ratio do not increase at least out to 0.5 R_vir. Our results clearly show that different sources of energy feedback leave distinct imprints in the enrichment pattern of the ICM, that are more evident when looking at cluster external regions.

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

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)

We aim at studying the abundance gradients along the Galactic disk and their dependence upon several parameters: a threshold in the surface gas density regulating star formation, the star formation efficiency, the timescale for the formation of the thin disk and the total surface mass density of the stellar halo. We test a model which considers a cosmological infall law. This law does not predict an inside-out disk formation, but it allows to well fit the properties of the solar vicinity. We study several cases. We find that to reproduce at the same time the abundance, star formation rate and surface gas density gradients along the Galactic disk it is necessary to assume an inside-out formation for the disk. The threshold in the gas density is not necessary and the same effect could be reached by assuming a variable star formation efficiency. A cosmologically derived infall law with an inside-out process for the disk formation and a variable star formation efficiency can indeed well reproduce all the properties of the disk. However, the cosmological model presented here does not have sufficient resolution to capture the requested inside-out formation for the disk.

We present results from two high-resolution hydrodynamical simulations of proto-cluster regions at z~2.1. The simulations have been compared to observational results for the socalled Spiderweb galaxy system, the core of a putative proto-cluster region at z = 2.16, found around a radio galaxy. The simulated regions have been chosen so as to form a poor cluster with M200~10^14 h-1 Msun (C1) and a rich cluster with M200~2x10^15 h-1 Msun (C2) at z = 0. The simulated proto-clusters show evidence of ongoing assembly of a dominating central galaxy. The stellar mass of the brightest cluster galaxy (BCG) of the C2 system is in excess with respect to observational estimates for the Spiderweb galaxy, with a total star formation rate which is also larger than indicated by observations. We find that the projected velocities of galaxies in the C2 cluster are consistent with observations, while those measured for the poorer cluster C1 are too low compared to the observed velocities. We argue that the Spiderweb complex resemble the high-redshift progenitor of a rich galaxy cluster. Our results indicate that the included supernovae feedback is not enough to suppress star formation in these systems, supporting the need of introducing AGN feedback. According to our simulations, a diffuse atmosphere of hot gas in hydrostatic equilibrium should already be present at this redshift, and enriched at a level comparable to that of nearby galaxy clusters. The presence of this gas should be detectable with future deep X-ray observations.

We present the first attempt to analyse the growth of the bar instability in stellar-gaseous disks evolving in a fully consistent cosmological scenario. We explored the role of the cosmology on pure stellar disks with different mass embedded in a cosmological dark matter halo. We deepened such a study by analysing the impact of different gas fractions and of the star formation onset. We found that in all these cases, the stellar bar arising inside the less massive disks, i.e., dark matter (DM)-dominated disks, is still living at redshift zero even if the gas fraction exceeds half of the disk mass. Such a bar is a genuine product of the cosmology. However, in the most massive disks there is a threshold value for their gas percentage and lower limit for the central gas concentration able to destroy the bar when the star formation rate is switched off. On the other hand in the simulations with star formation 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. Even if our results qualitatively agree with the classical ones, i.e. with criteria concerning bar instability derived outside the cosmological framework, the same criteria cannot be validated for the DM-dominated disks.

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.

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.

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 explore the growth and the evolution of the bar instability in stellar-gaseous disks embedded in a suitable dark matter halo evolving in a fully consistent cosmological framework. The aim of this paper is to point out the impact of different gas fractions on the bar formation, inside disks of different disk-to-halo mass ratio, and the role of the cosmological framework. We perform cosmological simulations with the same disk-to-halo mass ratios as in a previous work where the gas was not taken into account. We compare results of the new simulations with the previous ones to investigate the effect of the gas by analysing the morphology of the stellar and gaseous components, the stellar bar strength and the behaviour of its pattern speed. In our cosmological simulations, inside dark-matter dominated disks, a stellar bar, lasting 10 Gyr, is still living at z=0 even if the gaseous fraction exceeds half of the disk mass. However, in the most massive disks we find a threshold value (0.2) of the gas fraction able to destroy the bar. The stellar bar strength is enhanced by the gas and in the more massive disks higher gas fractions increase the bar pattern speed.

We study the origin of the diffuse stellar component (DSC) in 117 galaxy clusters extracted from a cosmological hydrodynamical simulation. We identify all galaxies present in the simulated clusters at 17 output redshifts, and then build the family trees for all the z=0 cluster galaxies. Then for each diffuse star particle identified at z=0, we look for the galaxy to which it once belonged at an earlier redshift, thus linking the presence of the diffuse stellar component to the galaxy formation history. The main results of our analysis are: (i) On average, half of the DSC star particles comes from galaxies associated with the family tree of the most massive galaxy (BCG), one quarter comes from the family trees of other massive galaxies, and the remaining quarter from dissolved galaxies. I.e., the formation of the DSC is parallel to the build-up of the BCG and other massive galaxies. (ii) Most DSC star particles become unbound during mergers in the formation history of the BGCs and of other massive galaxies, independent of cluster mass. Our results suggest that the tidal stripping mechanism is responsible only for a minor fraction of the DSC. (iii) At cluster radii larger than 250 kpc/h, the DSC fraction from the BCG is reduced and the largest contribution comes from the other massive galaxies; in the cluster outskirts, galaxies of all masses contribute to the DSC. (iv) The DSC does not have a preferred redshift of formation: however, most DSC stars are unbound at z<1. (v) The amount of DSC stars at z=0 does not corre late strongly with the global dynamical history of clusters, and increases weakly with cluster mass. [abridged]

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.

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.

We present a study of the galaxy population predicted by hydrodynamical simulations for a set of 19 galaxy clusters based on the GADGET-2 Tree+SPH code. These simulations include gas cooling, star formation, a detailed treatment of stellar evolution and chemical enrichment, as well as SN energy feedback in the form of galactic winds. We compute the spectro-photometric properties of the simulated galaxies. All simulations have been performed for two choices of the stellar initial mass function: a standard Salpeter IMF, and a top-heavier IMF. Several of the observational properties of the galaxy population in nearby clusters are reproduced fairly well by simulations. A Salpeter IMF is successful in accounting for the slope and the normalization of the color-magnitude relation for the bulk of the galaxy population. Simulated clusters have a relation between mass and optical luminosity which generally agrees with observations, both in normalization and slope. We find that galaxies are generally bluer, younger and more star forming in the cluster outskirts, thus reproducing the observational trends. However, simulated clusters have a total number of galaxies which is significantly smaller than the observed one, falling short by about a factor 2-3. Finally, the brightest cluster galaxies are always predicted to be too massive and too blue, when compared to observations, due to gas overcooling in the core cluster regions, even in the presence of a rather efficient SN feedback.

We study the efficiency and reliability of cluster mass estimators that are based on the projected phase-space distribution of galaxies in a cluster region. To this aim, we analyse a data-set of 62 clusters extracted from a concordance LCDM cosmological hydrodynamical simulation. Galaxies (or Dark Matter particles) are first selected in cylinders of given radius (from 0.5 to 1.5 Mpc/h) and ~200 Mpc/h length. Cluster members are then identified by applying a suitable interloper removal algorithm. Two cluster mass estimators are considered: the virial mass estimator (Mvir), and a mass estimator (Msigma) based entirely on the cluster velocity dispersion estimate. Mvir overestimates the true mass by ~10%, and Msigma underestimates the true mass by ~15%, on average, for sample sizes of > 60 cluster members. For smaller sample sizes, the bias of the virial mass estimator substantially increases, while the Msigma estimator becomes essentially unbiased. The dispersion of both mass estimates increases by a factor ~2 as the number of cluster members decreases from ~400 to ~20. The bias in the Mvir estimates is reduced in clusters without significant evidence for subclustering, and when only early-type galaxies are selected. Radially-dependent incompleteness can drastically affect Mvir estimates, but leaves the Msigma estimates almost unaffected. Other observational effects, like centering and velocity errors, and different observational apertures, have little effect on the mass estimates. (Abridged)

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.

We study the properties of the diffuse X-ray background by using the results of a cosmological hydrodynamical simulation of the concordance LambdaCDM model. The simulation follows gravitational and gas dynamics and includes a treatment of physical processes like radiative cooling, star formation and supernova feedback. From the simulation outputs, we produce a set of two-dimensional maps of the intergalactic medium X-ray emission integrated over redshift. We find that the signal in the soft (0.5-2 keV) band is lognormally distributed with a mean intensity of about 4 10^-12 erg s^-1 cm^-2 deg^-2; approximately 40 per cent of the emission originates from warm-hot gas (defined as baryons with 10^5<T<10^7 K), and 90 per cent comes from structures at z<0.9. Since the spectrum is soft, being mostly provided by the intergalactic medium at low temperature, the total mean intensity in the hard (2-10 keV) X-ray band is smaller by a factor of about 4. In order to constrain the physical processes included in our simulation, we compare our results with the observed upper limit (1.2 +/- 0.3) 10^-12 erg s^-1 cm^-2 deg^-2 of the soft X-ray emission due to diffuse gas. To this purpose, we remove the contributions of observable extended objects (groups and clusters of galaxies) from the simulated maps by adopting different detectability criteria which are calibrated on the properties of systems at intermediate redshifts observed by Chandra. We show that the simulated diffuse soft X-ray emission is consistent with the present observed upper limit. However, if future measurements will decrease the level of the unresolved X-ray background by a factor of two, a more efficient feedback mechanism should be required to suppress the soft emission of the gas residing within filaments and group-size haloes.

We discuss an extended set of Tree+SPH simulations of galaxy clusters, with the goal of investigating the interplay between numerical resolution effects and star-formation/feedback processes. The simulated clusters span the mass range (0.1-2.3) 10^{15}Msun/h, with mass resolution varying by several decades. At the highest achieved resolution, we resolve the virial region of a Virgo-like cluster with more than 2 million gas particles and with at least as many dark-matter (DM) particles. Our resolution study confirms that, in the absence of an efficient feedback mechanism, runaway cooling leads to about 35 per cent of baryons in clusters to be locked up in long lived stars at our highest resolution, with no evidence of convergence. However, including feedback causes the fraction of cooled baryons to converge at about 15 per cent already at modest resolution. Feedback also stabilizes other gas-related quantities, such as radial profiles of entropy, gas density and temperature, against variations due to changes in resolution. We also investigate the influence of the gravitational force softening length, and that of numerical heating of the gas induced by two-body encounters between DM and lighter gas particles. We show that simulations where more DM than gas particles are used, show a significantly enhanced efficiency of star formation at z>3. Our results are important for establishing and delineating the regime of numerical reliability of the present generation of hydrodynamical simulations of galaxy clusters.

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.

We study the baryon mass fraction in a set of hydrodynamical simulations of galaxy clusters performed using the Tree+SPH code GADGET-2. We investigate the dependence of the baryon fraction upon the radiative cooling, star formation, feedback through galactic winds, conduction and redshift. Both the cold stellar component and the hot X-ray emitting gas have narrow distributions that, at large cluster-centric distances r>R500, are nearly independent of the physics included in the simulations. Only the non-radiative runs reproduce the gas fraction inferred from observations of the inner regions (r ~ R2500) of massive clusters. When cooling is turned on, the excess star formation is mitigated by the action of galactic winds, but yet not by the amount required by observational data. The baryon fraction within a fixed overdensity increases slightly with redshift, independent of the physical processes involved in the accumulation of baryons in the cluster potential well. In runs with cooling and feedback, the increase in baryons is associated with a larger stellar mass fraction that arises at high redshift as a consequence of more efficient gas cooling. For the same reason, the gas fraction appears less concentrated at higher redshift. We discuss the possible cosmological implications of our results and find that two assumptions generally adopted, (1) mean value of Yb = fb / (Omega_b/Omega_m) not evolving with redshift, and (2) a fixed ratio between f_star and f_gas independent of radius and redshift, might not be valid. In the estimate of the cosmic matter density parameter, this implies some systematic effects of the order of Delta Omega_m/Omega_m < +0.15 for non-radiative runs and Delta Omega_m/Omega_m ~ +0.05 and < -0.05 for radiative simulations.

We perform a detailed analytical and numerical study of the dynamics of perturbations (vortex/aperiodic mode, Rossby and spiral-density waves) in 2D compressible disks with a Keplerian law of rotation. We draw attention to the process of spiral-density wave generation from vortices, discussing, in particular, the initial, most peculiar stages of wave emission. We show that the linear phenomenon of wave generation by vortices in smooth (without inflection points) shear flows found by using the so-called non-modal approach, is directly applicable to the present case. After an analytical non-modal description of the physics and characteristics of the spiral-density wave generation/propagation in the local shearing-sheet model, we follow the process of wave generation by small amplitude coherent circular vortex structures, by direct global numerical simulation, describing the main features of the generated waves.

The fluctuations in the Cosmic Microwave Background (CMB) intensity due to the Sunyaev-Zeldovich (SZ) effect are the sum of a thermal and a kinetic contribution. Separating the two components to measure the peculiar velocity of galaxy clusters requires radio and microwave observations at three or more frequencies, and knowledge of the temperature T_e of the intracluster medium weighted by the electron number density. To quantify the systematics of this procedure, we extract a sample of 117 massive clusters at redshift z=0 from an N-body hydrodynamical simulation, with 2x480^3 particles, of a cosmological volume 192 Mpc/h on a side of a flat Cold Dark Matter model with Omega_0=0.3 and Lambda=0.7. Our simulation includes radiative cooling, star formation and the effect of feedback and galactic winds from supernovae. We find that (1) our simulated clusters reproduce the observed scaling relations between X-ray and SZ properties; (2) bulk flows internal to the intracluster medium affect the velocity estimate by less than 200 km/s in 93 per cent of the cases; (3) using the X-ray emission weighted temperature, as an estimate of T_e, can overestimate the peculiar velocity by 20-50 per cent, if the microwave observations do not spatially resolve the cluster. For spatially resolved clusters, the assumptions on the spatial distribution of the ICM, required to separate the two SZ components, still produce a velocity overestimate of 10-20 per cent, even with an unbiased measure of T_e. Thanks to the large size of our cluster samples, these results set a robust lower limit of 200 km/s to the systematic errors that will affect upcoming measures of cluster peculiar velocities with the SZ effect.

The detection of excess of soft X-ray or Extreme Ultraviolet (EUV) radiation, above the thermal contribution from the hot intracluster medium (ICM), has been a controversial subject ever since the initial discovery of this phenomenon. We use a large--scale hydrodynamical simulation of a concordance $\Lambda$CDM model, to investigate the possible thermal origin for such an excess in a set of 20 simulated clusters having temperatures in the range 1--7 keV. Simulated clusters are analysed by mimicking the observational procedure applied to ROSAT--PSPC data, which for the first time showed evidences for the soft X-ray excess. For cluster--centric distances $0.4< R/R_{\rm vir}< 0.7$ we detect a significant excess in most of the simulated clusters, whose relative amount changes from cluster to cluster and, for the same cluster, by changing the projection direction. In about 30 per cent of the cases, the soft X-ray flux is measured to be at least 50 per cent larger than predicted by the one--temperature plasma model. We find that this excess is generated in most cases within the cluster virialized regions. It is mainly contributed by low--entropy and high--density gas associated with merging sub--halos, rather than to diffuse warm gas. Only in a few cases the excess arises from fore/background groups observed in projection, while no evidence is found for a significant contribution from gas lying within large--scale filaments. We compute the distribution of the relative soft excess, as a function of the cluster--centric distance, and compare it with the observational result by Bonamente et al. (2003) for the Coma cluster. Similar to observations, we find that the relative excess increases with the distance from the cluster center, with no significant excess detected for $R<0.4R_{\rm vir}$. (abridged)

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. In a recent paper Mazzotta et al. show that, if the cluster is thermally complex, T_ew fails at reproducing T_spec, and propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T_sl relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and \beta--model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X-ray temperature function could bias low \sigma_8 by 10-20 per cent. This may alleviate the tension between the value of \sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

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.

(ABRIDGED) We analyze the X-ray properties of a sample of local and high redshift galaxy clusters extracted from a large cosmological hydrodynamical simulation. This simulation has been realized using the Tree+SPH code GADGET-2 for a LambdaCDM model. In our analysis, we consider only objects with T_ew >2 keV and adopt an approach that mimics observations, associating with each measurement an error comparable with recent observations and providing best-fit results via robust techniques. Within the clusters, baryons are distributed among (i) a cold neutral phase, with a relative contribution that increases from less than 1 to 3 per cent at higher redshift, (ii) stars which contribute with about 20 per cent and (iii) the X-ray emitting plasma that contributes by 80 (76) per cent at z=0 (1) to the total baryonic budget. A depletion of the cosmic baryon fraction of ~7 (at z=0) and 5 (at z=1) per cent is measured at the virial radius, R_vir, in good agreement with adiabatic hydrodynamical simulations. We confirm that, also at redshift >0.5, power-law relations hold between gas temperature, T, bolometric luminosity, L, central entropy, S, gas mass, M_gas, and total gravitating mass, M_tot and that these relations are steeper than predicted by simple gravitational collapse. A significant, negative evolution in the L-T and L-M_tot relations and positive evolution in the S-T relation are detected at 0.5 < z < 1 in this set of simulated galaxy clusters. This is partially consistent with recent analyses of the observed properties of z>0.5 X-ray galaxy clusters. By fixing the slope to the values predicted by simple gravitational collapse, we measure at high redshift normalizations lower by 10-40 per cent in the L-T, M_tot-T, M_gas-T, f_gas-T and L-M_tot relations than the observed estimates.

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.

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.

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 show that DM halos, in n--body simulations, have a boundary layer (BL), separating bound from unbound mass. Let T(r) and W(r) be the kinetic and potential energies in shells of halos. We find that, in almost all halos: (i) The virial ratio R=-2T/W has at least one persistent (i.e. resolution independent) minimum r_c, such that, close to it, (ii) the function w=-d log W/d log r has a maximum and (iii) the relation R(r_c)=w(r_c) is almost exactly fulfilled. The radius r_c is the position of BL's, in halos found in simulations of TCDM and LCDM models, run using ART and GADGET codes at various resolutions. We find that 97% of the ~300 halos (per model) with M>4.2 10^14 M_s h^-1 owns a BL. Those with no BL are undergoing major mergings. The mass M_c enclosed in r_c almost coincides with the mass evaluated from velocities, according to virial theorem. Particles at r>r_c are not in virial equilibrium. Using r_c we have a density contrast \Delta_c for each halo. For each mass scale, \Delta_v=178 Omega_m^0.45 is within the range of \Delta_c's, but the spread in \Delta_c is wide and the average \Delta_c is ~25% smaller than the corresponding \Delta_v. The matching of properties derived under the assumption of spherical symmetry is a consequence of violent relaxation destroying features related to previous ellipsoidal non-linear growth. In turn, the spread of final \Delta_c's is an imprint of the different 3-D geometries and of the variable environment during collapses, as suggested by a comparison with Sheth & Tormen analysis.

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.

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.

The cosmological mass function problem is analyzed in full detail in the case of 1D gravity, with analytical, semi-analytical and numerical techniques. The extended Press & Schechter theory is improved by detailing the relation between smoothing radius and mass of the objects. This is done by introducing in the formalism the concept of a growth curve for the objects. The predictions of the extended Press & Schechter theory are compared to large N-body simulations of flat expanding 1D universes with scale-free power spectra of primordial perturbations. The collapsed objects in the simulations are located with a clump-finding algorithm designed to find regions that have undergone orbit crossing or that are in the multi-stream regime (these are different as an effect of the finite size of the multi-stream regions). It is found that the semi-analytical mass function theory, which has no free parameters, is able to recover the properties of collapsed objects both statistically and object by object. In particular, the predictions of regions in orbit crossing are optimized by the use of Gaussian filtering, while the use of sharp k-space filtering apparently allows to reproduce the larger multi-stream regions. The mass function theory does not reproduce well the clumps found with the standard friends-of-friends algorithm; however, the performance of this algorithm has not been thoroughly tested in the 1D cosmology. Our preliminary analyses of the 3D case confirms that the techniques developed in this paper are precious in understanding the cosmological mass function problem in 3D.

We study the cluster mass function in mixed dark matter (MDM) models, using two COBE normalized simulations with Omega_h = 0.26 and n=1.2, and Omega_h = 0.14 and n = 1.05, both with 2 massive nu's (MDM1 and MDM2, respectively). For the sake of comparison, we also simulate a CDM model with spectral index n=0.8 (TCDM), also COBE normalized. We argue that, in our non--hydro simulations, where CDM particles describe both actual CDM and baryons, the galaxy distribution traces CDM particles. Therefore, we use them to define clusters and their velocities to work out cluster masses. As CDM particles are more clustered than HDM and therefore have, in average, greater velocities, this leads to significant differences from PS predictions. Such predictions agree with simulations if both HDM and CDM are used to define clusters. Clusters defined through CDM in MDM models, instead, are less numerous than PS estimates, by a factor ~0.3, at the low mass end; the factor becomes \~0.6-0.8, depending on the mix, on intermediate mass scales (~4-5, h^-1 10^14 Msun) and almost vanishes on the high mass end. Therefore: (i) MDM models expected to overproduce clusters over intermediate scales are viable; (ii) the greater reduction factor at small scales agrees with the observational data dependence on the cluster mass M (which, however, may be partially due to sample incompleteness); (iii) the higher spectral normalization is felt at large scales, where MDM models produce more objects (hence, large clusters) than CDM. MDM1 even exceeds Donahue et al. (1998) findings, while MDM2 is consistent with them. (abridged)

If we treat the galaxies in published redshift catalogues as point sets, we may determine the generalized dimensions of such sets by standard means, outlined here. For galaxy separations up to about 5 Mpc, we find the dimensions of the galaxy set to be about 1.2, with not a strong indication of multifractality. For larger scales, out to about 30 Mpc, there is also good scaling with a dimension of about 1.8. For even larger scales, the data seem too sparse to be conclusive, but we find that the dimension is climbing as the scales increase. We report simulations that suggest a rationalization of such measurements, namely that in the intermediate range the scaling behavior is dominated by flat structures (pancakes) and that the results on the smallest scales are a reflection of the formation of density singularities.

We show that an array of singular density concentrations accounts for the scaling exponents observed in the luminous matter distribution in the Universe for scales below 10Mpc. This model provides a good description of the matter distribution on those smaller scales.

We use the void probability function to compare the redshift--space galaxy distribution in the Perseus--Pisces survey with artificial samples from N-body simulations of standard cold dark matter (CDM) and broken scale invariance (BSI) models. Galaxies are identified as residing in peaks of the evolved density field in such a way as to reproduce both the observed luminosity and two--point correlation functions. Using a similar approach, it was recently shown that the VPF can discriminate between CDM and a cold+hot dark matter (CHDM) model with cold/hot/baryon fractions = 0.6/0.3/0.1. Our main result is that both CDM (as expected from a previous analysis) and BSI fit observational data. The robustness of the result is checked against changing the observer's position in the simulations and the galaxy identification in the evolved density field. Therefore, while the void statistics is sensitive to the passage from CDM to CHDM (different spectrum and different nature of dark matter), it is not to the passage from CDM to BSI (different spectrum but same dark matter). On such a basis, we conjecture that the distribution of voids could be directly sensitive to the nature of dark matter, but scarcely sensitive to the shape of the transfer function.

We use the moments of counts of neighbors as given by the Generalized Correlation Integrals, to study the clustering properties of Dark Matter Halos (DH) in Cold Dark Matter (CDM) and Cold+Hot Dark Matter (CHDM) models. We compare the results with those found in the CfA and SSRS galaxy catalogs. We show that if we apply the analysis in redshift space, both models reproduce equally well the observed clustering of galaxies. Mass segregation is also found in the models: more massive DHs are more clustered compared with less massive ones. In redshift space, this mass segregation is reduced by a factor 2-3 due to the peculiar velocities. Observational catalogs give an indication of luminosity and size segregation, which is consistent with the predictions of the models. Because the mass segregation is smaller in redshift space, it is suggestive that the real luminosity or size segregation of galaxies could be significantly larger than what it is found in redshift catalogs.