Search

Context. Green valley (GV) galaxies are objects defined on a color-magnitude, or color-mass, diagram which are associated with a transition from a star-forming to a quiescent state (quenching), or vice versa (rejuvenation). Aims. We study the sub-mm emission of galaxies in the GV and link it with their physical evolutionary properties. Methods. We exploit a semi-analytic model (SAM) for galaxy evolution which includes a detailed treatment of dust production and evolution in galactic contexts. We model the observational properties of simulated galaxies by post-processing the SAM catalogs with the spectral synthesis and radiative transfer code GRASIL. Results. Our model produces a clear bimodality (and thus a GV) in the color-mass diagram, although some tensions arise when compared to observations. After having introduced a new criterion for identifying the GV in any dataset, we find that most blue and GV galaxies in our model are emitters at $250\,{\mu \rm m}$ (with fluxes $\geq 25\,{\rm mJy}$), in keeping with GAMA results. In the color-$S_{250\, \mu \rm m}$ diagram, GV galaxies exhibit intermediate colors yet maintain significant sub-mm fluxes. This occurs because, while specific star formation rates drop sharply during quenching, the dust content remains relatively high during the GV transition, powering sub-mm emission. Rejuvenating galaxies in the GV, which were previously red, have experienced a star formation burst that shifts their color to green, but their $S_{250\, \mu \rm m}$ fluxes remain low due to still low dust mass. Conclusions. Our galaxy evolution model highlights the delay between star formation and dust evolution, showing that sub-mm emission is not always a safe indicator of star formation activity, with quenching (rejuvenating) GV galaxies featuring relatively large (low) sub-mm emission., Comment: 12 pages, 13 figures. Submitted to A&A

Context: Molecular hydrogen ($\rm{H_2}$) is crucial in galaxy formation and evolution, serving as the main fuel for star formation (SF). In metal-enriched environments, $\rm{H_2}$ primarily forms on interstellar dust grain surfaces. However, due to the complexities of modelling this process, SF in cosmological simulations often relies on empirical or theoretical frameworks validated only in the Local Universe to estimate the abundance of $\rm{H_2}$. Aims: This study aims to model the connection between star, dust, and $\rm{H_2}$ formation processes in cosmological simulations. Methods: We include $\rm{H_2}$ formation on dust grain surfaces and account for molecule destruction and radiation shielding into the SF and feedback model MUPPI. Results: The model reproduces key properties of observed galaxies for stellar, dust, and $\rm{H_2}$ components. The cosmic density of $\rm{H_2}$ ($\rho_{\rm{H2}}$) peaks around $z=1.5$, then decreases by half towards $z=0$, showing milder evolution than observed. The $\rm{H_2}$ mass function since $z=2$ also shows gentler evolution. Our model successfully recovers the integrated molecular Kennicutt-Schmidt (mKS) law between surface star formation rate ($\Sigma_{\rm SFR}$) and surface $\rm{H_2}$ density ($\Sigma_{\rm H2}$) at $z=0$, already evident at $z=2$ with a higher normalization. We find hints of a broken power law with a steeper slope at higher $\Sigma_{\rm H2}$, aligning with some observational findings. Additionally, the $\rm{H_2}$-to-dust mass ratio in galaxies shows a decreasing trend with gas metallicity and stellar mass. The $\rm{H_2}$-to-dust mass fraction for the global galaxy population is higher at higher redshifts. The analysis of the atomic-to-molecular transition on a particle-by-particle basis suggests that gas metallicity cannot reliably substitute the dust-to-gas ratio in models simulating dust-promoted $\rm{H_2}$., Comment: Accepted for publication in A&A

The large-scale environment of the cosmic web is believed to impact galaxy evolution, but there is still no consensus regarding the mechanisms. We use a semi-analytic model (SAM) galaxy catalog to study the star formation and dust content of local galaxies in different cosmic environments of the cosmic web, namely voids, filaments, walls, and nodes. We find a strong impact of the environment only for galaxies with $M_{\rm stars}\lesssim10^{10.8}\, M_\odot$: the less dense the environment, the larger the star formation rate and dust content at fixed stellar mass. This is attributed to the fact that galaxies in less dense environments typically feature younger stellar populations, a slower evolution of their stellar mass and a delayed star formation compared to galaxies in denser environments. As for galaxies with $M_{\rm stars}\gtrsim 10^{10.8}\, M_\odot$ differences among environments are milder due to the disc instability (DI) driven supermassive black hole (SMBH) growth implemented in the SAM, which makes SMBH growth, and thus galaxy quenching, environment insensitive. We qualitatively test our predictions against observations by identifying environments in the SDSS-DR16 using dust masses derived from the GAMA survey. The agreement is encouraging, particularly at ${\rm log} \, M_{\rm stars}/M_\odot\gtrsim 10.5-11$, where sSFRs and dust masses appear quite environment-insensitive. This result confirms the importance of in situ growth channels of SMBHs., Comment: 21 pages, 12 figures, 3 tables. Accepted for publication in ApJ

Observations suggest that the amount of galactic dust in the Universe decreased by a factor $\sim 2-3$ during the last $\sim 8$ Gyr. However, cosmological models of galaxy evolution usually struggle to explain this decrease. Here we use the semi-analytic model (SAM) L-Galaxies2020 to show that this drop may be reproduced assuming standard prescriptions for dust production and evolution. We extend the SAM with i) a state-of-the-art dust model which adopts the two-size approximation and ii) a new disc instability criterion which triggers bulge and central black hole growth. The model reproduces some fundamental properties of the local galaxy population, such as the fraction of spheroid-dominated galaxies and some scaling relations involving dust. Moreover, the model predicts a galactic dust drop from $z \sim 1 \rightarrow 0$, which becomes closer to the observed one when adopting the new treatment of disc instabilities. This result is related to the newly implemented super-massive black hole growth during disc instabilities, which enhances the quenching of massive galaxies. Consequently, these objects feature a lower gas and dust content. We provide a census of the contribution of all the processes affecting the galactic dust content. Accretion is the dominant dust mass growth process. Destruction by supernovae, astration and ejection by winds have all a non-negligible role in decreasing the overall dust content in galaxies below $z \sim 1$. We also discuss predictions concerning extra-galactic dust, confirming that a sputtering efficiency lower than the canonical one is required to match the few available observations., Comment: 20 pages, 20 figures, 1 table. Accepted for publication in MNRAS

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

We study the evolution of dust in a cosmological volume using a hydrodynamical simulation in which the dust production is coupled with the MUPPI (MUlti Phase Particle Integrator) sub-resolution model of star formation and feedback. As for the latter, we keep as reference the model setup calibrated previously to match the general properties of Milky Way like galaxies in zoom-in simulations. However, we suggest that an increase of the star formation efficiency with the local dust to gas ratio would better reproduce the observed evolution of the cosmic star formation density. Moreover, the paucity of quenched galaxies at low redshift demands a stronger role of AGN feedback. We tune the parameters ruling direct dust production from evolved stars and accretion in the inter stellar medium to get scaling relations involving dust, stellar mass and metallicity in good agreement with observations. In low mass galaxies the accretion process is inefficient. As a consequence, they remain poorer in silicate and small grains than higher mass ones. We reproduce reasonably well the few available data on the radial distribution of dust outside the galactic region, supporting the assumption that the dust and gas dynamics are well coupled at galactic scales., Comment: MNRAS accepted

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

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

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

We present a detailed study of the shapes and alignments of different galaxy cluster components using hydrodynamical simulations. We compute shape parameters from the Dark Matter (DM) distribution, the galaxy members and the intra-cluster light (ICL). We assess how well the DM cluster shape can be constrained by means of the identified galaxy member positions and the ICL. Further, we address the dilution factor introduced when estimating the cluster elongation using weak-lensing stacking techniques, which arises due to the misalignment between the total surface mass distribution and the distribution of luminous tracers. The dilution is computed considering the alignment between the DM and the Brightest Cluster Galaxy, the galaxy members and the ICL. Our study shows that distributions of galaxy members and ICL are less spherical than the DM component, although both are well aligned with the semi-major axis of the later. We find that the distribution of galaxy members hosted in more concentrated subhalos is more elongated than the distribution of the DM. Moreover, these galaxies are better aligned with the dark matter component compared to the distribution of galaxies hosted in less concentrated subhalos. We conclude that the positions of galaxy members can be used as suitable tracers to estimate the cluster surface density orientation, even when a low number of members is considered. Our results provide useful information for interpreting the constraints on the shapes of galaxy clusters in observational studies., Comment: Submitted MNRAS. 12 pages, 7 figures (+ appendix)

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

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., Comment: 10 pages, 8 figures. MNRAS accepted. Comments welcome

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

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

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

The spectral energy distributions (SEDs) of low-mass low-metallicity (dwarf) galaxies are a challenging piece of the puzzle of galaxy formation in the near Universe. These SEDs show some particular features in the submillimeter to far-infrared wavelength range compared to normal, larger, and metal-richer galaxies that cannot be explained by the current models. These can be summarized as: a broadening of the IR peak, which implies a warmer dust component; an excess of emission in the submm ($\sim$500 $\mu$m), that causes a flattening of the submm/FIR slope; and a very low intensity of PAH emission features. With the aim of explaining these emission patterns, the SEDs of a sample of 27 simulated dwarf galaxies were calculated using the GRASIL-3D radiation transfer code. This code has the particularity that it separately treats the radiative transfer through dust grains within molecular clouds and within the cirrus, the dense and diffuse components of the gas phase, respectively. The simulated galaxies have stellar masses ranging from 10$^6$-10$^9$ M$_\odot$, and were obtained from a single simulation run within a Local Group environment with initial conditions from the CLUES project. We report a study of the IRAS, Spitzer and Herschel bands luminosities, and of the SFRs, dust, and gas (HI and H$_2$) mass contents. We find a satisfactory agreement with observational data, with GRASIL-3D naturally reproducing the specific spectral features mentioned above. We conclude that the GRASIL-3D two-component dust model gives a physical interpretation of the emission of dwarf galaxies: molecular clouds and cirrus represent the warm and cold dust components, respectively, needed to reproduce observational data., Comment: 22 pages, 18 figures. Accepted for publication in Astronomy & Astrophysics. Final version for journal

Many massive star forming disc galaxies in the redshift range 3 to 0.5 are observed to have a clumpy morphology showing giant clumps of size $\sim$1 kpc and masses of about $10^7M_{\odot}$ to $10^{10} M_{\odot}$. The nature and fate of these giant clumps is still under debate. In this work we use 19 high-resolution simulations of disc galaxies from the NIHAO sample to study the formation and the evolution of clumps in the discs of high redshift galaxies. We use mock HST - CANDELS observations created with the radiative transfer code GRASIL-3D to carry out, for the first time, a quantitative comparison of the observed fraction of clumpy galaxies and its evolution with redshift with simulations. We find a good agreement between the observed clumpy fraction and the one of the NIHAO galaxies. We find that dust attenuation can suppress intrinsically bright clumps and enhance less luminous ones. In our galaxy sample we only find clumps in light (u-band) from young stars but not in stellar mass surface density maps. This means that the NIHAO sample does not show clumpy stellar discs but rather a clumpy light distribution originating from clumpy star formation events. The clumps found in the NIHAO sample match observed age/color gradients as a function of distance from the galaxy center but they show no sign of inward migration. Clumps in our simulations disperse on timescales of a about a hundred Myr and their contribution to bulge growth is negligible., Comment: 17 pages, 18 figures, accepted by MNRAS

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

We 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., Comment: 22 pages, 11 figures, version accepted for publication in MNRAS, references and minor revisions added

Comparisons between observed and predicted strong lensing properties of galaxy clusters have been routinely used to claim either tension or consistency with $\Lambda$CDM cosmology. However, standard approaches to such cosmological tests are unable to quantify the preference for one cosmology over another. We advocate approximating the relevant Bayes factor using a marginal likelihood that is based on the following summary statistic: the posterior probability distribution function for the parameters of the scaling relation between Einstein radii and cluster mass, $\alpha$ and $\beta$. We demonstrate, for the first time, a method of estimating the marginal likelihood using the X-ray selected $z>0.5$ MACS clusters as a case in point and employing both N-body and hydrodynamic simulations of clusters. We investigate the uncertainty in this estimate and consequential ability to compare competing cosmologies, that arises from incomplete descriptions of baryonic processes, discrepancies in cluster selection criteria, redshift distribution, and dynamical state. The relation between triaxial cluster masses at various overdensities provide a promising alternative to the strong lensing test., Comment: 15 pages, 6 figures, 1 table, accepted to MNRAS

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., Comment: 14 pages, 11 figures, MNRAS in press. Minor editorial improvements

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., Comment: 13 pages, 12 figures, 4 tables

We have investigated the gas content of a sample of several hundred AGN host galaxies at z$<$1 and compared it with a sample of inactive galaxies, matched in bins of stellar mass and redshift. Gas masses have been inferred from the dust masses, obtained by stacked Herschel far-IR and sub-mm data in the GOODS and COSMOS fields, under reasonable assumptions and metallicity scaling relations for the dust-to-gas ratio. We find that AGNs are on average hosted in galaxies much more gas rich than inactive galaxies. In the vast majority of stellar mass bins, the average gas content of AGN hosts is higher than in inactive galaxies. The difference is up to a factor of ten higher in low stellar mass galaxies, with a significance of 6.5$\sigma$. In almost half of the AGN sample the gas content is three times higher than in the control sample of inactive galaxies. Our result strongly suggests that the probability of having an AGN activated is simply driven by the amount of gas in the host galaxy; this can be explained in simple terms of statistical probability of having a gas cloud falling into the gravitational potential of the black hole. The increased probability of an AGN being hosted by a star-forming galaxy, identified by previous works, may be a consequence of the relationship between gas content and AGN activity, found in this paper, combined with the Schmidt-Kennicutt law for star formation., Comment: Accepted for pubblication in MNRAS. 7 pages, 4 figures

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., Comment: 17 pages, 14 figures, accepted for publication on MNRAS, comments welcome

A serious concern for semi-analytical galaxy formation models, aiming to simulate multi-wavelength surveys and to thoroughly explore the model parameter space, is the extremely time consuming numerical solution of the radiative transfer of stellar radiation through dusty media. To overcome this problem, we have implemented an artificial neural network algorithm in the radiative transfer code GRASIL, in order to significantly speed up the computation of the infrared SED. The ANN we have implemented is of general use, in that its input neurons are defined as those quantities effectively determining the shape of the IR SED. Therefore, the training of the ANN can be performed with any model and then applied to other models. We made a blind test to check the algorithm, by applying a net trained with a standard chemical evolution model (i.e. CHE_EVO) to a mock catalogue extracted from the SAM MORGANA, and compared galaxy counts and evolution of the luminosity functions in several near-IR to sub-mm bands, and also the spectral differences for a large subset of randomly extracted models. The ANN is able to excellently approximate the full computation, but with a gain in CPU time by $\sim 2$ orders of magnitude. It is only advisable that the training covers reasonably well the range of values of the input neurons in the application. Indeed in the sub-mm at high redshift, a tiny fraction of models with some sensible input neurons out of the range of the trained net cause wrong answer by the ANN. These are extreme starbursting models with high optical depths, favorably selected by sub-mm observations, and difficult to predict a priori., Comment: 11 pages, 7 figures, MNRAS, accepted

We perform controlled numerical experiments to asses the effect of baryon mass loss on the inner structure of large galactic dark matter haloes. This mass expulsion is intended to mimic both the supernovae and AGN feedbacks, as well as the evolution of stellar populations. This study is meant in particular for precursors of massive Early Type Galaxies, wherein strong AGN feedback (often dubbed "QSO mode" in galaxy formation models) has been proposed to remove on a short timescale, of the order of a few dynamical times, a substantial fraction of their baryons. In a previous paper we evaluated the observational consequences (size increase) of this process on the galactic structure (Ragone-Figueroa & Granato 2011). Here we focus on the distribution of dark matter in the galactic region. It is shown that the inner region of the DM halo expands and its density profile flattens by a sizeable amount, with little dependence on the expulsion timescale. We also evaluate the effect of the commonly made approximation of treating the baryonic component as a potential that changes in intensity without any variation in shape. This approximation leads to some underestimates of the halo expansion and its slope flattening. We conclude that cuspy density profiles in ETGs could be difficult to reconcile with an effective AGN (or stellar) feedback during the evolution of these systems., Comment: 9 Pages, 5 figures, MNRAS accepted. Final version with minor revision

Observations performed in the last few years indicate that most massive early-type galaxies (ETGs) observed at redshift z>1 exhibit sizes smaller by a factor of a few than local ETGs of analogous stellar mass. We present numerical simulations of the effect of baryonic mass loss on the structure of a spheroidal stellar system, embedded in a dark matter halo. This process, invoked as a possible explanation of the observed size increase of ETGs since z \sim 2, could be caused either by QSO/starburst driven galactic winds, promptly ejecting from Early Type Galaxies (ETGs) the residual gas and halting star formation (galactic winds), or by stellar mass returned to the ISM in the final stages of stellar evolution. Indeed, we find that a conceivable loss of about 50% of the baryonic mass can produce a significant size increase. However, the puffing up due to galactic winds occurs when the stellar populations are much younger than the estimated ages >0.5 Gyr of compact high-z ETGs. Therefore, while it may have had a role in deciding the final structure of ETGs, it cannot explain the evolution observed so far of their size-mass relation; its signature should be searched for in much younger systems. Conversely, the mass loss due to stellar evolution could cause a relatively modest expansion of passively evolving stellar systems later on, contributing to, without dominating, the observed evolution of their mass-size relationship., Comment: MNRAS, accepted 2011 March 07. Received 2011 February 25; in original form 2011 January 20

We present a large set of theoretical isochrones, whose distinctive features mostly reside on the greatly improved treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase. Essentially, we have coupled the TP-AGB tracks described in Paper I, at their stages of pre-flash quiescent H-shell burning, with the evolutionary tracks for the previous evolutionary phases from Girardi et al. (2000). Theoretical isochrones for any intermediate value of age and metallicity are then derived by interpolation in the grids. We take care that the isochrones keep, to a good level of detail, the several peculiarities present in these TP-AGB tracks. Theoretical isochrones are then converted to about 20 different photometric systems -- including traditional ground-based systems, and those of recent major wide-field surveys such as SDSS, OGLE, DENIS, 2MASS, UKIDSS, etc., -- by means of synthetic photometry applied to an updated library of stellar spectra, suitably extended to include C-type stars. Finally, we correct the predicted photometry by the effect of circumstellar dust during the mass-losing stages of the AGB evolution, which allows us to improve the results for the optical-to-infrared systems, and to simulate mid- and far-IR systems such as those of Spitzer and AKARI. Access to the data is provided both via a web repository of static tables (http://stev.oapd.inaf.it/dustyAGB07 and CDS), and via an interactive web interface (http://stev.oapd.inaf.it/cmd) that provides tables for any intermediate value of age and metallicity, for several photometric systems, and for different choices of dust properties., Comment: 25 pages, accepted for publication in A&A, revised according to the latest referee's indications, isochrones are available at http://stev.oapd.inaf.it/cmd

Massive tori with $\approx 10^{8-9}M_{\odot}$ are predicted to extend on $\sim $100 pc scale around the centre of elliptical galaxy progenitors by a model of a supermassive black hole (SMBH) growth coeval to the spheroidal population of the host galaxy. Direct detection of such massive tori would cast light on a key physical condition that allows the rapid growth of SMBHs and the appearance of QSOs at high redshift. For this reason, we examine the detectability of such structures at substantial redshift with the Atacama Large Millimeter Array (ALMA). We propose that submillimeter galaxies (SMGs) are the best targets to test our predictions. In order to assess the observational feasibility, we estimate the expected number counts of SMGs with massive tori and check the detectability with the ALMA instrument, the unique facility which can resolve the central region of high redshift objects. Our work shows that ALMA will be able to resolve and detect high-$J$ ($J >$ 4) CO emissions from $\sim$100 pc scale extended massive tori up to $z\approx2$. Observations of lensed SMGs will yield excellent spatial resolution, allowing even to resolve their massive tori at higher redshift. We discuss further the detectability of the HCN molecule, as a better tracer of the high density gas expected in such tori. The final goal of these kind of observations is to pinpoint possible physical mechanisms that storage in the very central galactic regions very large amount gas on timescale of several 10$^{8}$ yr., Comment: 22 pages, 5 figures, accepted for publication in ApJ

In this invited contribution I review the justifications for the attempts, currently very popular, to include in semi-analytic models of galaxy formation prescriptions to describe the mutual link between the star formation and nuclear activity in galaxies, which has been for surprisingly long time neglected., Comment: Review in the Proceedings of the IAU Symposium 235, 2006, Eds F. Combes & J. Palous. 4 pages

We point out that in Granato & Danese 1994 and Granato et al. 1997 we predicted maximum observable sizes for the putative torus in NGC1068 of 10-20 pc, not "hundreds of parsecs" as stated by M. Elitzur in astro-ph/0512025., Comment: WEB note, 2 pages, no figures

We summarize our physical model for the early co-evolution of spheroidal galaxies and of active nuclei at their centers. Our predictions are in excellent agreement with a number of observables which proved to be extremely challenging for all the current semi-analytic models, including the sub-mm counts and the corresponding redshift distributions, and the epoch-dependent K-band luminosity function of spheroidal galaxies. Also, the black hole mass function and the relationship between the black hole mass and the velocity dispersion in the galaxy are nicely reproduced. The mild AGN activity revealed by X-ray observations of SCUBA sources is in keeping with our scenario, and testify the build up of SMBH triggered by intense star formation., Comment: 6 Pages, 5 figures, Invited talk at "Multiwavelength AGN surveys", Cozumel, Mexico, December 2003, eds. R. Mujica & R. Maiolino, Singapore: World Scientific, 2004

Granato et al (2004) have elaborated a physically grounded model exploiting the mutual feedback between star-forming spheroidal galaxies and the active nuclei growing in their cores to overcome, in the framework of the hierarchical clustering scenario for galaxy formation, one of the main challenges facing such scenario, the fact that massive spheroidal galaxies appear to have formed much earlier and faster than predicted by previous semi-analytical models. After having assessed the values of the two parameters that control the effect of the complex radiative transfer processes on the time-dependent SEDs we have compared the model predictions with a variety of IR to mm data. Our results support a rather strict continuity between objects where stars formed and evolved massive early-type galaxies, indicating that large spheroidal galaxies formed most of their stars when they were already assembled as single objects. The model is successful in reproducing the observed z distribution of K\le20 galaxies at z>1, in contrast with both the classical monolithic and the semi-analytic models, the ratio of star-forming to passively evolving spheroids and the counts and z distributions of EROs. The model also favourably compares with the ISOCAM 6.7 mu counts, with the corresponding z distribution, and with IRAC counts, which probe primarily the passive evolution phase, and with the submm SCUBA and MAMBO data, probing the active star-formation phase. The observed fraction of 24mu selected sources with no detectable emission in either the 8mu or R band nicely corresponds to the predicted surface density of star-forming spheroids with 8mu fluxes below the detection limit. Predictions for the z distributions of 24mu sources detected by MIPS surveys are pointed out. [Abridged], Comment: 12 pages, MNRAS, accepted

We estimate the contribution of AGNs and of their host galaxies to the infrared background. We use the luminosity function and evolution of AGNs recently determined by the hard X-ray surveys, and new Spectral Energy Distributions connecting the X-ray and the infrared emission, divided in intervals of absorption. These two ingredients allow us to determine the contribution of AGNs to the infrared background by using mostly observed quantities, with only minor assumptions. We obtain that AGN emission contributes little to the infrared background ($<$5% over most of the infrared bands), implying that the latter is dominated by star formation. However, AGN host galaxies may contribute significantly to the infrared background, and more specifically 10--20% in the 1--20$\mu$m range and $\sim$5% at $\lambda<60\mu m$. We also give the contribution of AGNs and of their host galaxies to the source number counts in various infrared bands, focusing on those which will be observed with Spitzer. We also report a significant discrepancy between the expected contribution of AGN hosts to the submm background and bright submm number counts with the observational constraints. We discuss the causes and implications of this discrepancy and the possible effects on the Spitzer far-IR bands., Comment: to appear in MNRAS, replaced with accepted version, paper shortened, results unchanged

Granato et al(2004) have elaborated a physically grounded model exploiting the mutual feedback between star-forming spheroidal galaxies and the active nuclei growing in their cores to overcome, within the hierarchical clustering scenario for galaxy formation, one of the main challenges facing such scenario, the fact that massive spheroidal galaxies appear to have formed earlier and faster than predicted by previous models. Adopting the choice by Granato et al for the parameters governing the history of the SF,of chemical abundances and of the gas and dust content of galaxies, we are left with only two parameters affecting the time and mass dependent SED of spheroidal galaxies. After complementing the model with a simple description of evolutionary properties of starburst, normal late-type galaxies and AGNs we have successfully compared the model with a broad variety of observational data, deep K-band, ISO, IRAS, SCUBA, radio counts, the corresponding redshift distributions, the IR background spectrum, and also with data for EROs. We also present detailed predictions for the GOODS and SWIRE surveys with the Spitzer Space Telescope. We find that the GOODS deep survey at 24$\mu$m and the SWIRE surveys at 70 and 160$\mu$m are likely to be severely confusion limited. The GOODS surveys in the IRAC channels are expected to resolve most of the background, to explore the full passive evolution phase of spheroidal galaxies and most of their active star-forming phase, detecting galaxies up to z\simeq 4. A substantial number of high z star-forming spheroidal galaxies should also be detected by the 24\mum SWIRE and GOODS surveys, while the 70 and 160\mum will be particularly useful to study the evolution of such galaxies in the range 1 \lsim z \lsim 2.[abridged], Comment: submitted to A&A

We determined C, N and alpha-element relative abundances in the gas surrounding six QSOs at an average redshift of < z > ~ 2.4, by studying six narrow associated absorption systems in UVES high-resolution spectra. We found five systems with a metallicity (measured by C/H) consistent or above the solar value. The ionization structure observed in the associated systems is clearly different from that of the intervening ones, indicating that the associated systems are influenced by the strong UV flux from the QSO. There is a possible correlation (anticorrelation) between [N/C] ([Si/C]) and [C/H] of the studied associated systems, and [N/C] >= 0 when [C/H] >= 0. We have compared these observational results with the predictions of a model simulating the joint evolution of QSOs and their spheroidal hosts. The agreement turns out to be very good, in particular, the case envisaging massive haloes and high star-formation rates recovers both the correlation between [N/C] and [C/H] and the anticorrelation for [Si/C] vs. [C/H]. Narrow associated absorption systems prove to be powerful tracers of the chemical abundances in gas belonging to high redshift spheroidal galaxies. The outflow of this same gas, triggered by the QSO feedback, is probably going to contribute to the early enrichment of the surrounding intergalactic medium. A larger statistics, possibly increasing the number of ionisation stages, chemical elements and the redshift range, would allow us to put firm constraints on detailed chemical evolution models of galaxies at high redshifts., Comment: 13 pages, 10 figures, submitted to MNRAS. Revised version after referee comments. Some changes in the paper structure, critical revision of errors, conclusions unchanged

The mutual feedback between star formation and nuclear activity in large spheroidal galaxies may be a key ingredient to overcome several difficulties plaguing current semi-analytic models for galaxy formation. We discuss some observational implications of the model by Granato et al. (2003) for the co-evolution of galaxies and active nuclei at their centers and stress the potential of the forthcoming surveys of the Sunyaev-Zeldovich effect on arcminute scales, down to $\mu$K levels, to investigate the early galaxy formation phases, difficult to access by other means., Comment: 6, pages, 1 figure, to appear in proc. of the meeting "Baryons on Cosmic Structures", Roma, October 20-21, 2003

We present a physically motivated model for the early co-evolution of massive spheroidal galaxies and active nuclei at their centers. Within dark matter halos, forming at the rate predicted by the canonical hierarchical clustering scenario, the gas evolution is controlled by gravity, radiative cooling, and heating by feedback from supernovae and from the growing active nucleus. Supernova heating is increasingly effective with decreasing binding energy in slowing down the star formation and in driving gas outflows. The more massive proto-galaxies virializing at earlier times are thus the sites of the faster star-formation. The correspondingly higher radiation drag fastens the angular momentum loss by the gas, resulting in a larger accretion rate onto the central black-hole. In turn, the kinetic energy carried by outflows driven by active nuclei can unbind the residual gas, thus halting both the star formation and the black-hole growth, in a time again shorter for larger halos. For the most massive galaxies the gas unbinding time is short enough for the bulk of the star-formation to be completed before type Ia supernovae can substantially increase the $Fe$ abundance of the interstellar medium, thus accounting for the $\alpha$-enhancement seen in the largest galaxies. The feedback from supernovae and from the active nucleus also determines the relationship between the black-hole mass and the mass, or the velocity dispersion, of the host galaxy, as well as the black-hole mass function. These and other model predictions are in excellent agreement with observations for a number of observables which proved to be extremely challenging for all the current semi-analytic models, including the sub-mm counts and the corresponding redshift distributions, and the epoch-dependent K-band luminosity function of spheroidal galaxies., Comment: 23 pages, 15 figures, revised after ApJ referee report

We revisit the nature of the FIR/Radio correlation by means of the most recent models for star forming galaxies. We model the IR emission with our population synthesis code, GRASIL (Silva et al. 1998). As for the radio emission, we revisit the simple model of Condon & Yin (1990). We find that a tightFIR/Radio correlation is natural when the synchrotron mechanism dominates over the inverse Compton, and the electrons cooling time is shorter than the fading time of the supernova rate. Observations indicate that both these conditions are met in star forming galaxies. However since the radio non thermal emission is delayed, deviations are expected both in the early phases of a starburst, when the radio thermal component dominates, and in the post-starburst phase, when the bulk of the NT component originates from less massive stars. This delay allows the analysis of obscured starbursts with a time resolution of a few tens of Myrs, unreachable with other star formation indicators. We suggest to complement the analysis of the deviations from the FIR/Radio correlation with the radio slope to obtain characteristic parameters of the burst. The analysis of a sample of compact ULIRGs shows that they are intense but transient starbursts, to which one should not apply usual SF indicators devised for constant SF rates. We also discuss the possibility of using the q- radio slope diagram to asses the presence of obscured AGN. A firm prediction of the models is an apparent radio excess during the post-starburst phase, which seems to be typical of a class of star forming galaxies in rich cluster cores. We discuss how deviations from the correlation, due to the evolutionary status of the starburst, affect the technique of photometric redshift determination widely used for high-z sources., Comment: accepted by A&A, 16 pages

Surveys at mm wavelengths emphasize rare classes of extragalactic radio sources characterized by spectra keeping flat or inverted up to high frequencies (such as blazars, GPS sources, advection dominated sources, or even gamma-ray burst afterglows), which are are difficult to single out at lower frequencies where the counts are dominated by far more numerous populations. Such surveys may lead to the discovery of populations with extreme properties, such as sources of the general GPS or blazar type but peaking at mm wavelengths, whose existence is hinted by currently available data. At relatively lower flux levels mm counts are dominated by high-redshift dusty galaxies, such as those detected by SCUBA and MAMBO surveys. These data have profound implications for our understanding of the formation and early evolution of galaxies. High sensitivity surveys at cm/mm wavelengths may also detect galaxy-scale Sunyaev-Zeldovich effects due to gas heated by central active nuclei., Comment: 3 figures, Proceedings conference "SRT: the impact of large antennas on Radioastronomy and Space Science", Cagliari 7-10 Novermber 2001

In view of the extensive evidence of a tight inter-relationship between spheroidal galaxies (and galactic bulges) and massive black holes hosted at their centers, a consistent model must deal jointly with the evolution of the two components. We describe one viable model, which successfully accounts for the local luminosity function of spheroidal galaxies, their photometric and chemical properties, deep galaxy counts in different wavebands, including those in the (sub)-mm region which proved to be critical for current semi-analytic models stemming from the standard hierarchical clustering picture, clustering properties of SCUBA galaxies, of EROs, and of LBGs, as well as for the local mass function of massive black holes and for quasar evolution. Predictions that can be tested by surveys carried out by SIRTF are presented., Comment: 6 pages, 2 fig. To appear in proceedings of "The Mass of Galaxies at Low and High Redshift", held Oct 2001 in Venice, eds. R. Bender and A. Renzini (ESO Astrophysics Symposia, Springer-Verlag)

We analyse very deep X-ray and mid-IR surveys in common areas of the Lockman Hole and the HDF North to study the sources of the X-ray background (XRB) and to test the standard obscured accretion paradigm. We detect with ISO a rich population of X-ray luminous sources with red optical colours, including a fraction identified with Extremely Red Objects (R-K > 5) and galaxies with SEDs typical of normal massive ellipticals or spirals at z ~ 1. The high 0.5-10 keV X-ray luminosities of these objects (1E43-1E45 erg/s) indicate that the ultimate energy source is gravitational accretion, while the X-ray to IR flux ratios and the X-ray spectral hardness show evidence of photoelectric absorption at low X-ray energies. An important hint on the physics comes from the mid-IR data at 6.7 and 15 um, well reproduced by model spectra of completely obscured quasars under standard assumptions and l.o.s. optical depths tau ~ 30-40. Other predictions of the standard XRB picture, like the distributions of intrinsic bolometric luminosities and the relative fractions of type-I and -II objects (1:3), are also consistent with our results. Obscured gravitational accretion is then confirmed as being responsible for the bulk of the X-ray background, since we detect in the IR the down-graded energy photoelectrically absorbed in X-rays: 63% of the faint 5-10 keV XMM sources are detected in the mid-IR by Fadda et al. (2001). However, although as much as 90% of the X-ray energy production could be converted to IR photons, no more than 20% of the Cosmic IR Background can be attributed to X-ray loud AGNs., Comment: 7 pages, 5 postscript figures, ApJ submitted

We present our model to interpret the SED of galaxies. The model for the UV to sub-mm SED is already well established (Silva et al 1998). We remind here its main features and show some applications. Recently we have extended the model to the radio range (Bressan et al 2001), and we have started to include the X-ray emission from the stellar component., Comment: 4 pages, to be published in "The link between stars and cosmology", 26-30 March, 2001, Puerto Vallarta, Mexico, by Kluwer, eds. M. Chavez, A. Bressan, A. Buzzoni, and D. Mayya

Star forming galaxies exhibit a variety of physical conditions, from quiescent normal spirals to the most powerful dusty starbursts. In order to study these complex systems, we need a suitable tool to analyze the information coming from observations at all wavelengths. We present a new spectro-photometric model which considers in a consistent way starlight as reprocessed by gas and dust. We discuss preliminary results to interpret some observed properties of VLIRGs., Comment: 8 pages, to be published in "The link between stars and cosmology", 26-30 March, 2001, Puerto Vallarta, Mexico, by Kluwer, eds. M. Chavez, A. Bressan, A. Buzzoni, and D. Mayya

We summarize recent advances in the determination of the cosmic history of star formation and other properties of high-z galaxies, and the relevance of this information in our understanding of the formation of structures. We emphasize the importance of dust reprocessing in the high--z universe, as demonstrated in particular by IR and sub-mm data. This demand a panchromatic approach to observations and suitable modelling tools. We spend also some words on expectations from future instruments., Comment: 8 pages, to be published in "The link between stars and cosmology", 26-30 March, 2001, Puerto Vallarta, Mexico, by Kluwer, eds. M. Chavez, A. Bressan, A. Buzzoni, and D. Mayya

We present detailed calculations of the magnification distribution, including both weak and strong lensing, using very recent solutions of the Dyer-Roeder (1973) equation for light propagation in a inhomogeneous universe with a cosmological constant and up-to-date models for the evolving cosmological distribution of dark matter halos. We address in particular the effect of a non-zero cosmological constant, of different density profiles of lenses, and of finite sizes of lensed sources. We show that, if dark matter fluctuations are normalized to the local cluster abundance, in the presence of a cosmological constant the optical depth for lensing {\em decreases} compared to the case of an Einstein-de Sitter universe, because halos in the relevant mass range are less abundant over a large redshift interval. We also discuss the differences in the magnification probability distributions produced by Navarro, Frenk & White (NFW) and by Singular Isothermal Sphere (SIS) density profiles of lenses. We find that the NFW lens is more efficient for moderate magnifications ($2\lsim A\lsim 4$), and less efficient for larger magnifications. Moreover, we discuss quantitatively the maximum magnification, $A_{\rm max}$, that can be achieved in the case of extended sources (galaxies) with realistic luminosity profiles, taking into account the possible ellipticity of the lens potential. Finally, we apply our results to a class of sources following the luminosity evolution typical for a unified scheme of QSO formation. We find that the lensed source counts at $850 \mu$m can be larger than the unlensed ones by several orders of magnitude at flux densities $\gta100 $mJy., Comment: 13 pages including 14 figures, MNRAS submitted

I review recent advances in the determination of the cosmic history of star formation, and its relevance in our understanding of the formation of structures. I emphasize the importance of dust reprocessing in the high--z universe, as demonstrated in particular by IR and sub-mm data. This demand a panchromatic approach to observations and suitable modelling tools. I summarize the basic requirements for these models and to what extent they are satisfied by models published so far., Comment: six pages, no figures, invited review for Stars, Gas, and Dust in Galaxies (eds. D. Alloin, K. Olsen, and G. Galaz)

We present models of photometric evolution of galaxies in which the effects of a dusty interstellar medium have been included with particular care. A chemical evolution code follows the star formation rate, the gas fraction and the metallicity, basic ingredients for the stellar population synthesis. The latter is performed with a grid of integrated spectra of simple stellar populations (SSP) of different ages and metallicities, in which the effects of dusty envelopes around asymptotic giant branch (AGB) stars are included. The residual fraction of gas in the galaxy is divided into two phases: the star forming molecular clouds and the cirrus. The relative amount is a model parameter. The molecular gas is sub--divided into clouds of given mass and radius: it is supposed that each SSP is born within the cloud and progressively escapes it. The emitted spectrum of the star forming molecular clouds is computed with a radiative transfer code. The cirrus emission is derived by describing the galaxy as an axially symmetric system, in which the local dust emissivity is consistently calculated as a function of the local field intensity due to the stellar component. Effects of very small grains, subject to temperature fluctuations, as well as PAHs are included. The model is compared and calibrated with available data of normal and starburst galaxies in the local universe, in particular new broad--band and spectroscopic ISO observations. It will be a powerful tool to investigate the star formation, the IMF, supernovae rate in nearby starbursts and normal galaxies, as well as to predict the evolution of luminosity functions of different types of galaxies at wavelengths covering four decades., Comment: 25 pages, 13 figures, accepted for publication on ApJ main journal