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.

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.

We report on our efforts to test the Einstein Equivalence Principle by measuring the gravitational redshift with the VLBI spacecraft RadioAstron, in an eccentric orbit around Earth with geocentric distances as small as $\sim$ 7,000 km and up to 350,000 km. The spacecraft and its ground stations are each equipped with stable hydrogen maser frequency standards, and measurements of the redshifted downlink carrier frequencies were obtained at both 8.4 and 15 GHz between 2012 and 2017. Over the course of the $\sim$ 9 d orbit, the gravitational redshift between the spacecraft and the ground stations varies between $6.8 \times 10^{-10}$ and $0.6 \times 10^{-10}$. Since the clock offset between the masers is difficult to estimate independently of the gravitational redshift, only the variation of the gravitational redshift is considered for this analysis. We obtain a preliminary estimate of the fractional deviation of the gravitational redshift from prediction of $\epsilon = -0.016 \pm 0.003_{\rm stat} \pm 0.030_{\rm syst}$ with the systematic uncertainty likely being dominated by unmodelled effects including the error in accounting for the non-relativistic Doppler shift. This result is consistent with zero within the uncertainties. For the first time, the gravitational redshift has been probed over such large distances in the vicinity of Earth. About three orders of magnitude more accurate estimates may be possible with RadioAstron using existing data from dedicated interleaved observations combining uplink and downlink modes of operation.

We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR.

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

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 present an approach to testing the gravitational redshift effect using the RadioAstron satellite. The experiment is based on a modification of the Gravity Probe A scheme of nonrelativistic Doppler compensation and benefits from the highly eccentric orbit and ultra-stable atomic hydrogen maser frequency standard of the RadioAstron satellite. Using the presented techniques we expect to reach an accuracy of the gravitational redshift test of order $10^{-5}$, a magnitude better than that of Gravity Probe A. Data processing is ongoing, our preliminary results agree with the validity of the Einstein Equivalence Principle.

A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths.

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.

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.

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.

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.

By means of the updated PARSEC database of evolutionary tracks of massive stars, we compute the integrated stellar light, the ionizing photon budget and the supernova rates of young simple stellar populations (SSPs), for different metallicities and IMF upper mass limits. Using CLOUDY we compute and include in the SSP spectra the neb- ular emission contribution. We also revisit the thermal and non-thermal radio emission contribution from young stars. Using GRASIL we can thus predict the panchromatic spectrum and the main recombination lines of any type of star-forming galaxy, including the effects of dust absorption and re-emission. We check the new models against the spectral energy distributions (SEDs) of selected well-observed nearby galaxies. From the best-fit models we obtain a consistent set of star formation rate (SFR) calibrations at wavelengths ranging from ultraviolet (UV) to radio. We also provide analytical calibrations that take into account the dependence on metallcity and IMF upper mass limit of the SSPs. We show that the latter limit can be well constrained by combining information from the observed far infrared, 24 {\mu}m, 33 GHz and H{\alpha} luminosities. Another interesting property derived from the fits is that, while in a normal galaxy the attenuation in the lines is significantly higher than that in the nearby continuum, in individual star bursting regions they are similar, supporting the notion that this effect is due to an age selective extinction. Since in these conditions the Balmer decrement method may not be accurate, we provide relations to estimate the attenuation from the observed 24 {\mu}m or 33 GHz fluxes. These relations can be useful for the analysis of young high redshift galaxies.

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

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 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 present a backward approach for the interpretation of the evolution of the near-infrared and the far-infrared luminosity functions across the redshift range 0<z<3. In our method, late-type galaxies are treated by means of a parametric phenomenological method based on PEP/HerMES data up to z~4, whereas spheroids are described by means of a physically motivated backward model. The spectral evolution of spheroids is modelled by means of a single-mass model, associated to a present-day elliptical with K-band luminosity comparable to the break of the local early-type luminosity function. The formation of proto-spheroids is assumed to occurr across the redshift range 1< z < 5. The key parameter is represented by the redshift z_0.5 at which half proto-spheroids are already formed. A statistical study indicates for this parameter values between z_0.5=1.5 and z_0.5=3. We assume as fiducial value z_0.5~2, and show that this assumption allows us to describe accourately the redshift distributions and the source counts. By assuming z_0.5 ~ 2 at the far-IR flux limit of the PEP-COSMOS survey, the PEP-selected sources observed at z>2 can be explained as progenitors of local spheroids caught during their formation. We also test the effects of mass downsizing by dividing the spheroids into three populations of different present-day stellar masses. The results obtained in this case confirm the validity of our approach, i.e. that the bulk of proto-spheroids can be modelled by means of a single model which describes the evolution of galaxies at the break of the present-day early type K-band LF.

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.

Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations match well the evolution of the stellar mass-halo mass relation, have a star forming main sequence that maintains a constant slope out to redshift z $\sim$ 2, and populate projections of the stellar mass - star formation - metallicity plane, similar to observed star forming disc galaxies. We discuss small differences between these projections in observational data and in simulations, and the possible causes for the discrepancies. The light-weighted stellar masses are in good agreement with the simulation values, the differences between the two varying between 0.06 dex and 0.20 dex. We also find a good agreement between the star formation rate tracer and the true (time-averaged) simulation star formation rates. Regardless if we use mass- or light-weighted quantities, our simulations indicate that bursty star formation cycles can account for the scatter in the star forming main sequence.

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.

We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands.

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

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

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 empirically test the relation between the SFR(LIR) derived from the infrared luminosity, LIR, and the SFR(Ha) derived from the Ha emission line luminosity using simple conversion relations. We use a sample of 474 galaxies at z = 0.06 - 0.46 with both Ha detection (from 20k zCOSMOS survey) and new far-IR Herschel data (100 and 160 {\mu}m). We derive SFR(Ha) from the Ha extinction corrected emission line luminosity. We find a very clear trend between E(B - V) and LIR that allows to estimate extinction values for each galaxy even if the Ha emission line measurement is not reliable. We calculate the LIR by integrating from 8 up to 1000 {\mu}m the SED that is best fitting our data. We compare SFR(Ha) with the SFR(LIR). We find a very good agreement between the two SFR estimates, with a slope of m = 1.01 \pm 0.03 in the SFR(LIR) vs SFR(Ha) diagram, a normalization constant of a = -0.08 \pm 0.03 and a dispersion of sigma = 0.28 dex.We study the effect of some intrinsic properties of the galaxies in the SFR(LIR)-SFR(Ha) relation, such as the redshift, the mass, the SSFR or the metallicity. The metallicity is the parameter that affects most the SFR comparison. The mean ratio of the two SFR estimators log[SFR(LIR)/SFR(Ha)] varies by approx. 0.6 dex from metal-poor to metal-rich galaxies (8.1 < log(O/H) + 12 < 9.2). This effect is consistent with the prediction of a theoretical model for the dust evolution in spiral galaxies. Considering different morphological types, we find a very good agreement between the two SFR indicators for the Sa, Sb and Sc morphologically classified galaxies, both in slope and normalization. For the Sd, irregular sample (Sd/Irr), the formal best-fit slope becomes much steeper (m = 1.62 \pm 0.43), but it is still consistent with 1 at the 1.5 sigma level, because of the reduced statistics of this sub-sample.

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.

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.

The cosmic star formation rate (CSFR), is an important clue to investigate the history of the assembly and evolution of galaxies. Here, we develop a method to study the CSFR from a purely theoretical point of view. Starting from detailed models of chemical evolution, we obtain the histories of star formation of galaxies of different morphological types. These histories are then used to determine the luminosity functions of the same galaxies by means of a spectro-photometric code. We obtain the CSFR under different hypothesis. First, we study the hypothesis of a pure luminosity evolution scenario, in which all galaxies are supposed to form at the same redshift and then evolve only in luminosity. Then we consider scenarios in which the number density or the slope of the LFs are assumed to vary with redshift. After comparison with available data we conclude that a pure luminosity evolution does not provide a good fit to the data, especially at very high redshift, although many uncertainties are still present in the data. On the other hand, a variation in the number density of ellipticals and spirals as a function of redshift can provide a better fit to the observed CSFR. We also explore cases of variable slope of the LFs with redshift and variations of number density and slope at the same time. We cannot find any of those cases which can improve the fit to the data respect to the solely number density variation. Finally, we compute the evolution of the average cosmic metallicity in galaxies with redshift.

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.

The spectral energy distribution of galaxies is a complex function of the star formation history and geometrical arrangement of stars and gas in galaxies. The computation of the radiative transfer of stellar radiation through the dust distribution is time-consuming. This aspect becomes unacceptable in particular when dealing with the predictions by semi-analytical galaxy formation models populating cosmological volumes, to be then compared with multi-wavelength surveys. Mainly for this aim, we have implemented an artificial neural network algorithm into the spectro-photometric and radiative transfer code GRASIL in order to compute the spectral energy distribution of galaxies in a short computing time. This allows to avoid the adoption of empirical templates that may have nothing to do with the mock galaxies output by models. The ANN has been implemented to compute the dust emission spectrum (the bottleneck of the computation), and separately for the star-forming molecular clouds and the diffuse dust (due to their different properties and dependencies). We have defined the input neurons effectively determining their emission, which means this implementation has a general applicability and is not linked to a particular galaxy formation model. We have trained the net for the disc and spherical geometries, and tested its performance to reproduce the SED of disc and starburst galaxies, as well as for a semi-analytical model for spheroidal galaxies. We have checked that for this model both the SEDs and the galaxy counts in the Herschel bands obtained with the ANN approximation are almost superimposed to the same quantities obtained with the full GRASIL. We conclude that this method appears robust and advantageous, and will present the application to a more complex SAM in another paper.

The presence of dust strongly affects the way we see galaxies and also the chemical abundances we measure in gas. It is therefore important to study he chemical evolution of galaxies by taking into account dust evolution. We aim at performing a detailed study of abundance ratios of high redshift objects and their dust properties. We focus on Lyman-Break galaxies (LBGs) and Quasar (QSO) hosts and likely progenitors of low- and high-mass present-day elliptical galaxies, respectively. We have adopted a chemical evolution model for elliptical galaxies taking account the dust production from low and intermediate mass stars, supernovae Ia, supernovae II, QSOs and both dust destruction and accretion processes. By means of such a model we have followed the chemical evolution of ellipticals of different baryonic masses. Our model complies with chemical downsizing. We made predictions for the abundance ratios versus metallicity trends for models of differing masses that can be used to constrain the star formation rate, initial mass function and dust mass in observed galaxies. We predict the existence of a high redshift dust mass-stellar mass relationship. We have found a good agreement with the properties of LBGs if we assume that they formed at redshift z=2-4. In particular, a non-negligible amount of dust is needed to explain the observed abundance pattern. We studied the QSO SDSS J114816, one of the most distant QSO ever observed (z=6.4), and we have been able to reproduce the amount of dust measured in this object. The dust is clearly due to the production from supernovae and the most massive AGB stars as well as from the grain growth in the interstellar medium. The QSO dust is likely to dominate only in the very central regions of the galaxies and during the early development of the galactic wind.

We present the analysis of Spitzer-IRS spectra of four early-type galaxies, NGC 1297, NGC 5044, NGC 6868, and NGC 7079, all classified as LINERs in the optical bands. Their IRS spectra present the full series of H2 rotational emission lines in the range 5--38 microns, atomic lines, and prominent PAH features. We investigate the nature and origin of the PAH emission, characterized by unusually low 6 -- 9/11.3 microns inter-band ratios. After the subtraction of a passive early type galaxy template, we find that the 7 -- 9 microns spectral region requires dust features not normally present in star forming galaxies. Each spectrum is then analyzed with the aim of identifying their components and origin. In contrast to normal star forming galaxies, where cationic PAH emission prevails, our 6--14 microns spectra seem to be dominated by large and neutral PAH emission, responsible for the low 6 -- 9/11.3 microns ratios, plus two broad dust emission features peaking at 8.2 microns and 12 microns. Theses broad components, observed until now mainly in evolved carbon stars and usually attributed to pristine material, contribute approximately 30-50% of the total PAH flux in the 6--14 microns region. We propose that the PAH molecules in our ETGs arise from fresh carbonaceous material which is continuously released by a population of carbon stars, formed in a rejuvenation episode which occurred within the last few Gyr. The analysis of the MIR spectra allows us to infer that, in order to maintain the peculiar size and charge distributions biased to large and neutral PAHs, this material must be shocked, and excited by the weak UV interstellar radiation field of our ETG.

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 use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

We investigate the atomic (HI) and molecular (H_2) Hydrogen content of normal galaxies by combining observational studies linking galaxy stellar and gas budgets to their host dark matter (DM) properties, with a physically grounded galaxy formation model. This enables us to analyse empirical relationships between the virial, stellar, and gaseous masses of galaxies and explore their physical origins. Utilising a semi-analytic model (SAM) to study the evolution of baryonic material within evolving DM halos, we study the effects of baryonic infall and various star formation and feedback mechanisms on the properties of formed galaxies using the most up-to-date physical recipes. We find that in order to significantly improve agreement with observations of low-mass galaxies we must suppress the infall of baryonic material and exploit a two-phase interstellar medium (ISM), where the ratio of HI to H_2 is determined by the galactic disk structure. Modifying the standard Schmidt-Kennicutt star formation law, which acts upon the total cold gas in galaxy discs and includes a critical density threshold, and employing a star formation law which correlates with the H_2 gas mass results in a lower overall star formation rate. This in turn, allows us to simultaneously reproduce stellar, HI and H_2 mass functions of normal galaxies.

We use our model for the formation and evolution of galaxies within a two-phase galaxy formation scenario, showing that the high-redshift domain typically supports the growth of spheroidal systems, whereas at low redshifts the predominant baryonic growth mechanism is quiescent and may therefore support the growth of a disc structure. Under this framework we investigate the evolving galaxy population by comparing key observations at both low and high-redshifts, finding generally good agreement. By analysing the evolutionary properties of this model, we are able to recreate several features of the evolving galaxy population with redshift, naturally reproducing number counts of massive star-forming galaxies at high redshifts, along with the galaxy scaling relations, star formation rate density and evolution of the stellar mass function. Building upon these encouraging agreements, we make model predictions that can be tested by future observations. In particular, we present the expected evolution to z=2 of the super-massive black hole mass function, and we show that the gas fraction in galaxies should decrease with increasing redshift in a mass, with more and more evolution going to higher and higher masses. Also, the characteristic transition mass from disc to bulge dominated system should decrease with increasing redshift.

We study in detail the mid-infrared Spitzer-IRS spectrum of M 87 in the range 5 to 20 micron. Thanks to the high sensitivity of our Spitzer-IRS spectra we can disentangle the stellar and nuclear components of this active galaxy. To this end we have properly subtracted from the M 87 spectrum, the contribution of the underlying stellar continuum, derived from passive Virgo galaxies in our sample. The residual is a clear power-law, without any additional thermal component, with a zero point consistent with that obtained by high spatial resolution, ground based observations. The residual is independent of the adopted passive template. This indicates that the 10 micron silicate emission shown in spectra of M 87 can be entirely accounted for by the underlying old stellar population, leaving little room for a possible torus contribution. The MIR power-law has a slope alpha ~ 0.77-0.82 (S$_\nu\propto\nu^{-\alpha}$), consistent with optically thin synchrotron emission.

We introduce a new technique based on artificial neural networks which allows us to make accurate predictions for the spectral energy distributions (SEDs) of large samples of galaxies, at wavelengths ranging from the far-ultra-violet to the sub-millimetre and radio. The neural net is trained to reproduce the SEDs predicted by a hybrid code comprised of the GALFORM semi-analytical model of galaxy formation, which predicts the full star formation and galaxy merger histories, and the GRASIL spectro-photometric code, which carries out a self-consistent calculation of the SED, including absorption and emission of radiation by dust. Using a small number of galaxy properties predicted by GALFORM, the method reproduces the luminosities of galaxies in the majority of cases to within 10% of those computed directly using GRASIL. The method performs best in the sub-mm and reasonably well in the mid-infrared and the far-ultra-violet. The luminosity error introduced by the method has negligible impact on predicted statistical distributions, such as luminosity functions or colour distributions of galaxies. We use the neural net to predict the overlap between galaxies selected in the rest-frame UV and in the observer-frame sub-mm at z=2. We find that around half of the galaxies with a 850um flux above 5 mJy should have optical magnitudes brighter than R_AB < 25 mag. However, only 1% of the galaxies selected in the rest-frame UV down to R_AB < 25 mag should have 850um fluxes brighter than 5 mJy. Our technique will allow the generation of wide-angle mock catalogues of galaxies selected at rest-frame UV or mid- and far-infrared wavelengths.

We propose and test a scenario for the assembly and evolution of luminous matter in galaxies which substantially differs from that adopted by other semianalytic models. As for the dark matter (DM), we follow the detailed evolution of halos within the canonical LCDM cosmology using standard Montecarlo methods. However, when overlaying prescriptions for baryon evolution, we take into account an effect pointed out in the past few years by a number of studies mostly based on intensive N-body simulations, namely that typical halo growth occurs in two phases: an early, fast collapse phase featuring several major merger events, followed by a late, quiescent accretion onto the halo outskirts. We propose that the two modes of halo growth drive two distinct modes for the evolution of baryonic matter, favoring the development of the spheroidal and disc components of galaxies, respectively. We test this idea using the semianalytic technique. Our galaxy formation model envisages an early coevolution of spheroids and the central supermassive black holes, already tested in our previous works, followed by a relatively quiescent growth of discs around the preformed spheroids. In this exploratory study, we couple our model to the spectrophotometric code GRASIL, and compare our results on several properties of the local galaxy population with observations, Finding an encouraging agreement.

We fit the near-infrared to radio spectral energy distributions of a sample of 30 luminous and ultra-luminous infrared galaxies with models that include both starburst and AGN components. The aim of the work was to determine important physical parameters for this kind of objects such as the optical depth towards the luminosity source, the star formation rate, the star formation efficiency and the AGN fraction. We found that although about half of our sample have best-fit models that include an AGN component, only 30 % have an AGN which accounts for more than 10 % of the infrared luminosity whereas all have an energetically dominant starburst. Our models also determine the mass of dense molecular gas. Assuming that this mass is that traced by the HCN molecule, we reproduce the observed linear relation between HCN luminosity and infrared luminosity found by Gao and Solomon (2004). However, our derived conversion factor between HCN luminosity and the mass of dense molecular gas is a factor of 2 smaller than that assumed by these authors. Finally, we find that the star formation efficiency falls as the starburst ages.

We present photometric evolution models of galaxies, in which, in addition to the stellar component, the effects of an evolving dusty interstellar medium have been included with particular care. Starting from the work of Calura, Pipino & Matteucci (2008), in which chemical evolution models have been used to study the evolution of both the gas and dust components of the interstellar medium in the solar neighbourhood, elliptical and irregular galaxies, it has been possible to combine these models with a spectrophotometric stellar code that includes dust reprocessing (GRASIL) (Silva et al. 1998) to analyse the evolution of the spectral energy distributions (SED) of these galaxies. We test our models against observed SEDs both in the local universe and at high redshift and use them to predict how the percentage of reprocessed starlight evolves for each type of galaxy. The importance of following the dust evolution is investigated by comparing our results with those obtained by adopting simple assumptions to treat this component.

We use 16 micron, Spitzer-IRS, blue peakup photometry of 50 early-type galaxies in the Coma cluster to define the mid-infrared colour-magnitude relation. We compare with recent simple stellar population models that include the mid-infrared emission from the extended, dusty envelopes of evolved stars. The Ks-[16] colour in these models is very sensitive to the relative population of dusty Asymptotic Giant Branch (AGB) stars. We find that the passively evolving early-type galaxies define a sequence of approximately constant age (~10 Gyr) with varying metallicity. Several galaxies that lie on the optical/near-infrared colour-magnitude relation do not lie on the mid-infrared relation. This illustrates the sensitivity of the Ks-[16] colour to age. The fact that a colour-magnitude relation is seen in the mid-infrared underlines the extremely passive nature of the majority (68%) of early-type galaxies in the Coma cluster. The corollary of this is that 32% of the early-type galaxies in our sample are not `passive', insofar as they are either significantly younger than 10 Gyr or they have had some rejuvenation episode within the last few Gyr.

We present initial results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z>3 by means of deep near-IR spectroscopy. Gas metallicities are measured, for an initial sample of nine star forming galaxies at z~3.5, by means of optical nebular lines redshifted into the near-IR. Stellar masses are accurately determined by using Spitzer-IRAC data, which sample the rest-frame near-IR stellar light in these distant galaxies. When compared with previous surveys, the mass-metallicity relation inferred at z~3.5 shows an evolution much stronger than observed at lower redshifts. The evolution is prominent even in massive galaxies, indicating that z~3 is an epoch of major action in terms of star formation and metal enrichment also for massive systems. There are also indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.5 is difficult to reconcile with the predictions of some hierarchical evolutionary models. Such discrepancies suggest that at z>3 galaxies are assembled mostly with relatively un-evolved sub-units, i.e. small galaxies with low star formation efficiency. The bulk of the star formation and metallicity evolution probably occurs once small galaxies are already assembled into bigger systems.

In a previous work, three bright MIR/radio sources were discovered in the nuclear region of NGC1365. We here confirm that these sources are young and massive ``embedded'' clusters, and derive their parameters, such as extinction, age and mass. Using ISAAC and VISIR at the VLT we obtained NIR and MIR maps and LR spectra. The dataset is first interpreted by comparing the observations with images and spectra of the close-by young cluster R136 in the LMC and then by using model predictions for both the nebular emission lines and the spectral energy distribution of the sources. We produce maps of the region containing the three sources in the R, J, Ks, L' bands and at 12.8micro. We also provide spectra in K, L and N. The spectral energy distribution of the three sources rises with wavelength. Emission lines from ionised hydrogen and molecular hydrogen are detected, as well as PAH emission. Conspicuous [NeII]12.8 line is also present, while neither the [ArIII] 8.9 nor the [SIV] 10.4 lines are detected. This provides a stringent constraint on the age of the sources: they are relatively evolved young clusters (6-8 Myr). Owing to their ionising photon rates and ages, they must be extremely massive clusters (around 10^7 solar masses). Their MIR spectral energy distribution suggests the presence of two components: (1) an optically thin component, with a continuum comparable to that of R136, and (2) an optically thick component which might be related to subsequent or on-going episodes of star formation. These sources are good candidates for evolving according to a bi-modal hydrodynamical regime, in which matter is trapped at the centre of a compact and massive cluster.

Motivated by the recent observational and theoretical evidence that long Gamma-Ray Bursts (GRBs) are likely associated with low metallicity, rapidly rotating massive stars, we examine the cosmological star formation rate (SFR) below a critical metallicity Z_crit Z_sun/10 - Z_sun/5, to estimate the event rate of high-redshift long GRB progenitors. To this purpose, we exploit a galaxy formation scenario already successfully tested on a wealth of observational data on (proto)spheroids, Lyman break galaxies, Lyman alpha emitters, submm galaxies, quasars, and local early-type galaxies. We find that the predicted rate of long GRBs amounts to about 300 events/yr/sr, of which about 30 per cent occur at z>~6. Correspondingly, the GRB number counts well agree with the bright SWIFT data, without the need for an intrinsic luminosity evolution. Moreover, the above framework enables us to predict properties of the GRB host galaxies. Most GRBs are associated with low mass galaxy halos M_H<~10^11 M_sun, and effectively trace the formation of small galaxies in such halos. The hosts are young, with age smaller than 5*10^7 yr, gas rich, but poorly extincted (A_V<~0.1) because of their chemical immaturity; this also implies high specific SFR and quite extreme alpha-enhancement. Only the minority of hosts residing in large halos with M_H>~10^12 M_sun have larger extinction (A_V~0.7-1), SFRs exceeding 100 M_sun/yr and can be detected at submm wavelengths. Most of the hosts have UV magnitudes in the range -20 <~M_1350<~ -16, and Lyman alpha luminosity in the range 2*10^40 <~L_Lya<~2*10^42 erg/s. GRB hosts are thus tracing the faint end of the luminosity function of Lyman break galaxies and Lyman alpha emitters.

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.

We present predictions for the evolution of the galaxy luminosity function, number counts and redshift distributions in the IR based on the Lambda-CDM cosmological model. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF seems to be necessary to explain the number counts of faint sub-mm galaxies. We compare the model with observational data from the SPITZER Space Telescope, with the model parameters fixed at values chosen before SPITZER data became available. We find that the model matches the observed evolution in the IR remarkably well over the whole range of wavelengths probed by SPITZER. In particular, the SPITZER data show that there is strong evolution in the mid-IR galaxy luminosity function over the redshift range z ~ 0-2, and this is reproduced by our model without requiring any adjustment of parameters. On the other hand, a model with a normal IMF in starbursts predicts far too little evolution in the mid-IR luminosity function, and is therefore excluded.

We present preliminary results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z~3 by means of deep near-IR spectroscopy. Gas metallicities and stellar masses are measured for an initial sample of nine star forming galaxies at z~3.3. When compared with previous surveys, the mass-metallicity relation inferred at z~3.3 shows an evolution significantly stronger than observed at lower redshifts. There are also some indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered as the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.3 is difficult to reconcile with the predictions of some hierarchical evolutionary models. We shortly discuss the possible implications of such discrepancies.

We fit the near-infrared to radio spectral energy distributions of 30 luminous and ultra-luminous infrared galaxies with pure starburst models or models that include both starburst and AGN components to determine important physical parameters for this population of objects. In particular we constrain the optical depth towards the luminosity source, the star formation rate, the star formation efficiency and the AGN fraction. We find that although about half of our sample have best-fit models that include an AGN component, only 30% have an AGN which accounts for more than 10% of the infrared luminosity, whereas all have an energetically dominant starburst. Our derived AGN fractions are generally in good agreement other measurements based in the mid-infrared line ratios measured by Spitzer IRS, but lower than those derived from PAH equivalent widths or the mid-infrared spectral slope. Our models determine the mass of dense molecular gas via the extinction required to reproduce the SED. Assuming that this mass is that traced by HCN, we reproduce the observed linear relation between HCN and infrared luminosities found by Gao & Solomon. We also find that the star formation efficiency, defined as the current star formation rate per unit of dense molecular gas mass, is enhanced in the ULIRGs phase. If the evolution of ULIRGs includes a phase in which an AGN contributes an important fraction to the infrared luminosity, this phase should last an order of magnitude less time than the starburst phase. Because the mass of dense molecular gas which we derive is consistent with observations of the HCN molecule,it should be possible to estimate the mass of dense, star-forming molecular gas in such objects when molecular line data are not available.

As a constraint for new starburst/AGN models of IRAS bright galaxies we determine the radio spectra of 31 luminous and ultraluminous IRAS galaxies (LIRGs/ULIRGs). We construct the radio spectra using both new and archival data. From our sample of radio spectra we find that very few have a straight power-law slope. Although some sources show a flattening of the radio spectral slope at high frequencies the average spectrum shows a steepening of the radio spectrum from 1.4 to 22.5 GHz. This is unexpected because in sources with high rates of star formation we expect flat spectrum, free-free emission to make a significant contribution to the radio flux at higher radio frequencies. Despite this trend the radio spectral indices between 8.4 and 22.5 GHz are flatter for sources with higher values of the FIR-radio flux density ratio q, when this is calculated at 8.4 GHz. Therefore, sources that are deficient in radio emission relative to FIR emission (presumably younger sources) have a larger thermal component to their radio emission. However, we find no correlation between the radio spectral index between 1.4 and 4.8 GHz and q at 8.4 GHz. Because the low frequency spectral index is affected by free-free absorption, and this is a function of source size for a given mass of ionized gas, this is evidence that the ionized gas in ULIRGs shows a range of densities. The youngest LIRGs and ULIRGs are characterized by a larger contribution to their high-frequency radio spectra from free-free emission. However, the youngest sources are not those that have the greatest free-free absorption at low radio frequencies. The sources in which the effects of free-free absorption are strongest are instead the most compact sources. Although these have the warmest FIR colours, they are not necessarily the youngest sources.

We have studied the submillimetre (submm) properties of the following classes of near-infrared (NIR)-selected massive galaxies at high redshifts: BzK-selected star-forming galaxies (BzKs); distant red galaxies (DRGs); and extremely red objects (EROs). We used the SCUBA HAlf Degree Extragalactic Survey (SHADES), the largest uniform submm survey to date. Partial overlap of SIRIUS/NIR images and SHADES in SXDF has allowed us to identify 4 submm-bright NIR-selected galaxies, which are detected in the mid-infrared, 24 micron, and the radio, 1.4 GHz. We find that all of our submm-bright NIR-selected galaxies satisfy the BzK selection criteria, except for one galaxy whose B-z and z-K colours are however close to the BzK colour boundary. Two of the submm-bright NIR-selected galaxies satisfy all of the selection criteria we considered, i.e. they belong to the BzK-DRG-ERO overlapping population, or `extremely red' BzKs. Although these extremely red BzKs are rare (0.25 arcmin^{-2}), up to 20 % of this population could be submm galaxies. This fraction is significantly higher than that found for other galaxy populations studied here. Via a stacking analysis, we have detected the 850 micron flux of submm-faint BzKs and EROs in our SCUBA maps. While the contribution of z~2 BzKs to the submm background is about 10--15 % and similar to that from EROs typically at z~1, BzKs have a higher fraction (~30 %) of submm flux in resolved sources compared with EROs and submm sources as a whole. From the SED fitting analysis for both submm-bright and submm-faint BzKs, we found no clear signature that submm-bright BzKs are experiencing a specifically luminous evolutionary phase, compared with submm-faint BzKs. An alternative explanation might be that submm-bright BzKs are more massive than submm-faint ones.

We explore the onset of star formation in the early Universe, exploiting the observations of high-redshift Lyman-break galaxies (LBGs) and Lyman alpha emitters (LAEs), in the framework of the galaxy formation scenario elaborated by Granato et al. (2004) already successfully tested against the wealth of data on later evolutionary stages. Complementing the model with a simple, physically plausible, recipe for the evolution of dust attenuation in metal poor galaxies we reproduce the luminosity functions (LFs) of LBGs and of LAEs at different redshifts. This recipe yields a much faster increase with galactic age of attenuation in more massive galaxies, endowed with higher star formation rates. These objects have therefore shorter lifetimes in the LAE and LBG phases, and are more easily detected in the dusty submillimeter bright phase. The short UV bright lifetimes of massive objects strongly mitigate the effect of the fast increase of the massive halo density with decreasing redshift, thus accounting for the weaker evolution of the LBG LF, compared to that of the halo mass function, and the even weaker evolution between z~6 and z~3 of the LAE LF. LAEs are on the average expected to be younger, with lower stellar masses, and associated to less massive halos than LBGs. Finally, we show that the intergalactic medium can be completely reionized at redshift z~6-7 by massive stars shining in protogalactic spheroids with halo masses from a few 10^10 to a few 10^11 M_sun, showing up as faint LBGs with magnitude in the range -17<M_1350<-20, without resorting to any special stellar initial mass function.

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.

We present a quantitative analysis of the astrophysical and cosmological information that can be extracted from the many important wide-area, shallow surveys that will be carried out in the next few years. Our calculations combine the predictions of the physical model by Granato et al. (2004) for the formation and evolution of spheroidal galaxies with up-to-date phenomenological models for the evolution of starburst and normal late-type galaxies and of radio sources. We compute the expected number counts and the redshift distributions of these source populations separately and then focus on proto-spheroidal galaxies. For the latter objects we predict the counts and redshift distributions of strongly lensed sources at 250, 350, 500, and 850 micron, the angular correlation function of sources detected in the surveys considered, the angular power spectra due to clustering of sources below the detection limit in Herschel and Planck surveys. An optimal survey for selecting strongly lensed proto-spheroidal galaxies is described, and it is shown how they can be easily distinguished from the other source populations. We also discuss the detectability of the imprints of the 1-halo and 2-halo regimes on angular correlation functions and clustering power spectra, as well as the constraints on cosmological parameters that can be obtained from the determinations of these quantities. The novel data relevant to derive the first sub-millimeter estimates of the local luminosity functions of starburst and late-type galaxies, and the constraints on the properties of rare source populations, such as blazars, are also briefly described.

The mid infrared emission of early type galaxies traces the presence of intermediate age stellar populations as well as even tiny amounts of ongoing star formation. Here we discuss high S/N Spitzer IRS spectra of a sample of Virgo early type galaxies, with particular reference to NGC 4435. We show that, by combining mid infrared spectroscopic observations with existing broad band fluxes, it is possible to obtain a very clean picture of the nuclear activity in this galaxy.

We are performing a systematic study of the properties of early-type galaxies in the mid infrared spectral region with the Spitzer Space Telescope. We present here high S/N Spitzer IRS spectra of 17 Virgo early-type galaxies. Thirteen objects of the sample (76%) show a pronounced broad feature (above 10 microns) which is spatially extended and likely of stellar origin. We argue that this feature is (mostly) due to silicate emission from circumstellar envelopes of asymptotic giant branch (AGB) stars. The remaining 4 objects, namely NGC 4486, NGC 4636, NGC 4550 and NGC 4435, are characterized by various levels and type of activity.

We analyse the attenuation properties of a sample of UV selected galaxies, with the use of the spectrophotometric model Grasil. In particular, we focus on the relation between dust attenuation and the reddening in the UV spectral region. We show that a realistic modelling of geometrical distribution of dust and of the different population of stars can explain the UV reddening of normal spiral galaxies also with a standard Milky Way dust. Our results clearly underline that it is fundamental to take into account that younger stars suffer a higher attenuation than older stars (the age-dependent extinction) because stars are born in more-than-average dusty environments. In this work we also find that the concentration of young stars on the galactic plane of spirals has a relevant impact on the expected UV colours, impact that has not been explored before this paper. Finally, we discuss the role of IMF in shaping the relation between UV reddening and dust attenuation, and we show that a Kroupa IMF is more consistent with observed data than the classical Salpeter IMF.

The ~800 optically unseen (R>25.5) 24mum-selected sources in the complete Spitzer First Look Survey sample (Fadda et al. 2006) with F[24mum]>0.35 mJy are found to be very strongly clustered. If, as indicated by several lines of circumstantial evidence, they are ultraluminous far-IR galaxies at z ~ [1.6-2.7], the amplitude of their spatial correlation function is very high. The associated comoving clustering length is estimated to be r_0=14.0_{-2.4}^{+2.1} Mpc, value which puts these sources amongst the most strongly clustered populations of our known universe. Their 8mum-24mum colours suggest that the AGN contribution dominates above F[24mum] ~ 0.8 mJy, consistent with earlier analyses. The properties of these objects (number counts, redshift distribution, clustering amplitude) are fully consistent with those of proto-spheroidal galaxies in the process of forming most of their stars and of growing their active nucleus, as described by the Granato et al. (2004) model. In particular, the inferred space density of such galaxies at z ~ 2 is much higher than what expected from most semi-analytic models. Matches of the observed projected correlation function w(\theta) with models derived within the so-called Halo Occupation Scenario show that these sources have to be hosted by haloes more massive than ~10^{13.4} M_\odot. This value is significantly higher than that for the typical galactic haloes hosting massive elliptical galaxies, suggesting a duration of the starburst phase of massive high-redshift dusty galaxies of T_B ~ 0.5 Gyr.

We present new HI observations of the nearby dwarf galaxy NGC 3741. This galaxy has an extremely extended HI disk, which allows us to trace the rotation curve out to unprecedented distances in terms of the optical disk: we reach 42 B-band exponential scale lengths or about 7 kpc. The HI disk is strongly warped, but the warp is very symmetric. The distribution and kinematics are accurately derived by building model data cubes, which closely reproduce the observations. In order to account for the observed features in the data cube, radial motions of the order of 5-13 km/s are needed. They are consistent with an inner bar of several hundreds of pc and accretion of material in the outer regions. The observed rotation curve was decomposed into its stellar, gaseous and dark components. The Burkert dark halo (with a central constant density core) provides very good fits. The dark halo density distribution predicted by the LambdaCDM theory fails to fit the data, unless NGC 3741 is a 2.5-sigma exception to the predicted relation between concentration parameter and virial mass and has at the same time a high value of the virial mass (though poorly constrained) of 10$^{11}$ solar masses. Noticeably, MOND seems to be consistent with the observed rotation curve. Scaling up the contribution of the gaseous disk also gives a good fit.

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.

We present a population synthesis study of NGC4435, an early-type Virgo galaxy interacting with NGC4438. We combine new spectroscopic observations obtained with the Spitzer Space Telescope IRS instrument with IRAC archival data and broad band data from the literature. The IRS spectrum shows prominent PAH features, low ionization emission lines and H_2 rotational lines arising from the dusty circumnuclear disk characterizing this galaxy. The central SED, from X-ray to radio, is well fitted by a model of an exponential burst superimposed on an old simple stellar population. From the lack of high excitation nebular lines, the [NeIII]15.5/[NeII]12.8 ratio, the temperature of molecular hydrogen, and the fit to the full X-ray to radio SED we argue that the present activity of the galaxy is driven by star formation alone. The AGN contribution to the ionizing flux is constrained to be less than 2%. The age of the burst is found to be around 190 Myr and it is fully consistent with the notion that the star formation process has been triggered by the interaction with NGC4438. The mass involved in the rejuvenation episode turns out to be less than 1.5% of the stellar galaxy mass sampled in a 5" central aperture. This is enough to render NGC4435 closely similar to a typical interacting early-type galaxy with inverted CaII[H+K] lines that will later turn into a typical cluster E+A galaxy and enforces the notion that these objects are the result of a recent rejuvenation episode rather than a genuine delayed formation.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

We show that our previously proposed anti-hierarchical baryon collapse scenario for the joint evolution of black holes and host galaxies predicts quasar luminosity functions at redshifts 1.5<z<6 and local properties in nice agreement with observations. In our model the quasar activity marks and originates the transition between an earlier phase of violent and heavily dust-enshrouded starburst activity promoting rapid black hole growth, and a later phase of almost passive evolution; the former is traced by the submillimeter-selected sources, while the latter accounts for the high number density of massive galaxies at substantial redshifts z>1.5, the population of Extremely Red Objects, and the properties of local ellipticals.

We have obtained high S/N Spitzer IRS spectra of 17 Virgo early-type galaxies that lie on the colour-magnitude relation of passively evolving galaxies in the cluster. To flux calibrate these extended sources we have devised a new procedure that allows us to obtain the intrinsic spectral energy distribution and to disentangle resolved and unresolved emission within the same object. Thirteen objects of the sample (76%) show a pronounced broad silicate feature (above 10micron) which is spatially extended and likely of stellar origin, in agreement with model predictions. The other 4 objects (24%) are characterized by different levels of activity. In NGC 4486 (M87) the line emission and the broad silicate emission are evidently unresolved and, given also the typical shape of the continuum, they likely originate in the nuclear torus. NGC 4636 show emission lines superimposed to extended silicate emission (i.e. likely of stellar origin, pushing the percentage of galaxies with silicate emission to 82%). Finally NGC 4550 and NGC 4435 are characterized by PAH and line emission, arising from a central unresolved region.

We present a hierarchical galaxy formation model which can account for the number counts of sources detected through their emission at sub-millimetre wavelengths. The first stage in our approach is an ab initio calculation of the star formation histories for a representative sample of galaxies, which is carried out using the semi-analytical galaxy formation model GALFORM. These star formation histories are then input into the spectro-photometric code GRASIL, to produce a spectral energy distribution for each galaxy. Dust extinction and emission are treated self consistently in our model, without having to resort to ad-hoc assumptions about the amount of attenuation by dust or the temperature at which the dust radiates. We argue that it is necessary to modify the form of the stellar initial mass function in starbursts in order to match the observed number of sub-mm sources, if we are to retain the previous good matches enjoyed between observations and model predictions in the local universe. We also list some other observational tests that have been passed by our model.

We show that the ABC scenario we proposed for the co-evolution of spheroids and QSOs predicts accretion rates and masses of supermassive black holes in sub-mm galaxies in keeping with recent X-ray determinations. These masses are well below the local values, and those predicted by alternative models. The observed column densities may be mostly due to ISM in the galaxy. The contribution of the associated nuclear activity to the X-ray background is likely negligible, while they may contribute a sizeable fraction $\sim 10 %$ to hard-X cumulative counts at the faintest observed fluxes.

At variance with most semi-analytic models, in the Anti-hierarchical Baryon Collapse scenario (Granato et al. 2001, 2004) the main driver of the galaxy formation and evolution is not the merging sequence but are baryon processes. This approach emphasizes, still in the framework of the hierarchical clustering paradigm for dark matter halos, feedback processes from supernova explosions and from active nuclei, that tie together star formation in spheroidal galaxies and the growth of black holes at their centers. We review some recent results showing the remarkably successful predictive power of this scenario, which allows us to account for the evolution with cosmic time of a broad variety of properties of galaxies and active nuclei, which proved to be very challenging for competing models.

We present high signal to noise ratio Spitzer Infrared Spectrograph observations of 17 Virgo early-type galaxies. The galaxies were selected from those that define the colour-magnitude relation of the cluster, with the aim of detecting the silicate emission of their dusty, mass-losing evolved stars. To flux calibrate these extended sources we have devised a new procedure that allows us to obtain the intrinsic spectral energy distribution and to disentangle resolved and unresolved emission within the same object. We have found that thirteen objects of the sample (76%) are passively evolving galaxies with a pronounced broad silicate feature which is spatially extended and likely of stellar origin, in agreement with model predictions. The other 4 objects (24%) are characterized by different levels of activity. In NGC 4486 (M 87) the line emission and the broad silicate emission are evidently unresolved and, given also the typical shape of the continuum, they likely originate in the nuclear torus. NGC 4636 shows emission lines superimposed on extended (i.e. stellar) silicate emission, thus pushing the percentage of galaxies with silicate emission to 82%. Finally, NGC 4550 and NGC 4435 are characterized by polycyclic aromatic hydrocarbon (PAH) and line emission, arising from a central unresolved region. A more detailed analysis of our sample, with updated models, will be presented in a forthcoming paper.

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.

In this paper we analyze the evolution of actively star forming galaxies in the mid-infrared (MIR). This spectral region, characterized by continuum emission by hot dust and by the presence of strong emission features generally ascribed to polycyclic aromatic hydrocarbon (PAH) molecules, is the most strongly affected by the heating processes associated with star formation and/or active galactic nuclei (AGN). Following the detailed observational characterization of galaxies in the MIR by ISO, we have updated the modelling of this spectral region in our spectro-photometric model GRASIL (Silva et al. 1998). In the diffuse component we have updated the treatment of PAHs according to the model by Li & Draine (2001). As for the dense phase of the ISM associated with the star forming regions, the molecular clouds, we strongly decrease the abundance of PAHs as compared to that in the cirrus, basing on the observational evidences of the lack or weakness of PAH bands close to the newly formed stars, possibly due to the destruction of the molecules in strong UV fields. The robustness of the model is checked by fitting near infrared to radio broad band spectra and the corresponding detailed MIR spectra of a large sample of galaxies (Lu et al. 2003), at once. With this model, we have analyzed the larger sample of actively star forming galaxies by Dale et al. (2000). We show that the observed trends of galaxies in the ISO-IRAS-Radio color-color plots can be interpreted in terms of different evolutionary phases of star formation activity, and the consequent different dominance in the spectral energy distribution (SED) of the diffuse or dense phase of the ISM.

We present new results from a multi-wavelength model of galaxy formation, which combines a semi-analytical treatment of the formation of galaxies within the CDM framework with a sophisticated treatment of absorption and emission of radiation by dust. We find that the model, which incorporates a top-heavy IMF in bursts, agrees well with the evolution of the rest-frame far-UV luminosity function over the range z=0-6, with the IR number counts in all bands measured by SPITZER, and with the observed evolution of the mid-IR luminosity function for z=0-2.

We show that the observed Velocity Dispersion Function of E/S0 galaxies matches strikingly well the distribution function of virial velocities of massive halos virializing at z > 1.5, as predicted by the standard hierarchical clustering scenario in a \LambdaCDM cosmology, for a constant ratio sigma/V_vir = 0.55 \pm 0.05, close to the value expected at virialization if it typically occurred at z > 3. This strongly suggests that dissipative processes and later merging events had little impact on the matter density profile. Adopting the above sigma/V_vir ratio, the observed relationships between photometric and dynamical properties which define the fundamental plane of elliptical galaxies, such as the luminosity-sigma (Faber-Jackson) and the luminosity-effective radius relations, as well as the M_BH-sigma relation, are nicely reproduced. Their shapes turn out to be determined by the mutual feedback of star-formation (and supernova explosions)and nuclear activity, along the lines discussed by Granato et al. (2004). To our knowledge, this is the first semi-analytic model for which simultaneous fits of the fundamental plane relations and of the epoch-dependent luminosity function of spheroidal galaxies have been presented.

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.

The possible contributions of the various classes of extragalactic sources (including, in addition to the canonical radio sources, GHz Peaked Spectrum sources, advection-dominated sources, starburst galaxies, high-redshift proto-spheroidal galaxies) to the arcminute scale fluctuations measured by the CBI, BIMA, and ACBAR experiments are discussed. At 30 GHz, fluctuations due to radio sources undetected by ancillary low-frequency surveys may be higher than estimated by the CBI and BIMA groups. High-redshift dusty galaxies, whose fluctuations may be strongly enhanced by the effect of clustering, could contribute to the BIMA excess signal, and dominate at 150 GHz (the ACBAR frequency). Moreover, in the present data situation, the dust emission of these high-redshift sources set an unavoidable limit to the detection of primordial CMB anisotropies at high multipoles, even at frequencies as low as $\simeq 30$ GHz. It is concluded that the possibility that the excess power at high multipoles is dominated by unsubtracted extragalactic sources cannot be ruled out. On the other hand, there is room for a contribution from the Sunyaev-Zeldovich effect within clusters of galaxies, with a density fluctuation amplitude parameter $\sigma_8$ consistent with the values preferred by current data.

We present here a model for simulating the panchromatic spectral energy distribution of galaxies, which aims to be a complete tool to study the complex multi-wavelength picture of the universe. The model take into account all important components that concur to the SED of galaxies at wavelengths from X-rays to the radio. We review the modeling of each component and provide several applications, interpreting observations of galaxy of different types at all the wavelengths.

We present predictions for the abundance of sub-mm galaxies (SMGs) and Lyman-break galaxies (LBGs) in the $\Lambda$CDM cosmology. A key feature of our model is the self-consistent calculation of the absorption and emission of radiation by dust. The new model successfully matches the LBG luminosity function, as well reproducing the properties of the local galaxy population in the optical and IR. The model can also explain the observed galaxy number counts at $850\mum$, but only if we assume a top-heavy IMF for the stars formed in bursts. The predicted redshift distribution of SMGs depends relatively little on their flux over the range 1-$10\mjy$, with a median value of $z\approx 2.0$ at a flux of $5\mjy$, in very good agreement with the recent measurement by Chapman et al The counts of SMGs are predicted to be dominated by ongoing starbursts. However, in the model these bursts are responsible for making only a few per cent of the stellar mass locked up in massive ellipticals at the present day.

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]

We present a physical model for the coevolution of massive spheroidal galaxies and active nuclei at their centers. Supernova heating is increasingly effective in slowing down the star formation and in driving gas outflows in smaller and smaller dark matter halos. Thus the more massive protogalaxies virializing at early times are the sites of faster star formation. The correspondingly higher radiation drag causes a faster angular momentum loss by the gas and induces a larger accretion rate onto the central black hole. In turn, the kinetic energy of the outflows powered by the active nuclei can unbind the residual gas in a time shorter for larger halos. The model accounts for a broad variety of dynamical, photometric and metallicity properties of early-type galaxies, for the M_BH -- \sigma relation and for the local supermassive black-hole mass function.

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.

We have performed a detailed analysis of the local super-massive black-hole (SMBH) mass function based on both kinematic and photometric data and derived an accurate analytical fit in the range 10^6 <= (M_BH/M_sun) <= 5*10^9. We find a total SMBH mass density of (4.2+/-1.1)*10^5 M_sun/Mpc^3, about 25% of which is contributed by SMBHs residing in bulges of late type galaxies. Exploiting up-to-date luminosity functions of hard X-ray and optically selected AGNs, we have studied the accretion history of the SMBH population. If most of the accretion happens at constant \dot{M_BH}/M_BH the local SMBH mass function is fully accounted for by mass accreted by X-ray selected AGNs, with bolometric corrections indicated by current observations and a standard mass-to-light conversion efficiency \epsilon ~10%. The analysis of the accretion history highlights that the most massive BHs accreted their mass faster and at higher redshifts (z>1.5), while the lower mass BHs responsible for most of the hard X-ray background have mostly grown at z<1.5. The accreted mass function matches the local SMBH mass function if \epsilon ~0.09(+0.04,-0.03) and the Eddington ratio \lambda=L/L_Edd \~0.3(+0.3,-0.1) (68% confidence errors). The visibility time, during which AGNs are luminous enough to be detected by the currently available X-ray surveys, ranges from ~0.1 Gyr for present day BH masses M_BH(z=0) ~10^6 M_sun to ~0.3 Gyr for M_BH(z=0) >= 10^9 M_sun. The mass accreted during luminous phases is >= 25-30% even if we assume extreme values of \epsilon (\epsilon \~0.3-0.4). An unlikely fine tuning of the parameters would be required to account for the local SMBH mass function accomodating a dominant contribution from 'dark' BH growth (due, e.g., to BH coalescence).

We present adaptive optics assisted spectroscopy of three quasars obtained with NACO at VLT. The high angular resolution achieved with the adaptive optics (~0.08"), joined to the diagnostic power of near-IR spectroscopy, allow us to investigate the properties of the innermost 100 pc of these quasars. In the quasar with the best adaptive optics correction, PG1126-041, we spatially resolve the Pa-alpha emission within the nuclear 100 pc. The comparison with higher excitation lines suggests that the narrow Pa-alpha emission is due to nuclear star formation. The inferred intensity of the nuclear star formation (13 M(sun)/yr) may account for most of the far-IR luminosity observed in this quasar.

We present detailed predictions for the confusion noise due to extragalactic sources in the far-IR/(sub)-millimeter channels of ESA/ISO, NASA/Spitzer, ESA/Herschel and ESA/Planck satellites, including the contribution from clustering of unresolved SCUBA galaxies. Clustering is found to increase the confusion noise, compared to the case of purely Poisson fluctuations, by 10-15% for the lowest frequency (i.e. lowest angular resolution) Spitzer and Herschel channels, by 25-35% for the 175 micron ISOPHOT channel, and to dominate in the case of Planck/HFI channels at nu>143GHz. Although our calculations make use of a specific evolutionary model (Granato et al. 2004), the results are strongly constrained by the observed counts and by data on the redshift distribution of SCUBA sources, and therefore are not expected to be heavily model dependent. The main uncertainty arises from the poor observational definition of the source clustering properties. Two models have been used for the latter: a power-law with constant slope and a redshift-independent comoving correlation length,r_0, and the standard theoretical model for clustering evolution in a LambdaCDM universe, with a redshift-dependent bias factor. In both cases, the clustering amplitude has been normalized to yield a unit angular correlation function at theta_0=1-2 arcsec for 850 micron sources fainter than 2 mJy, consistent with the results by Peacock et al. (2000). This normalization yields, for the first model, r_0=8.3$ Mpc/h, and, for the second model, an effective mass of dark matter haloes in which these sources reside of M_halo=1.8*10^{13} M_sun/h. These results are consistent with independent estimates for SCUBA galaxies and for other, likely related, sources.

We have performed 23 GHz VLA observations of 7 compact, luminous infrared galaxies, selected to have evidence of starburst activity. New and published multi-frequency data are combined to obtain the spectral energy distributions of all 7 galaxies from the near-infrared to the radio (at 1.4 GHz). These SEDs are compared with new models, for dust enshrouded galaxies, which account for both starburst and AGN components. In all 7 galaxies the starburst provides the dominant contribution to the infrared luminosity; in 4 sources no contribution from an AGN is required. Although AGN may contribute up to 50 percent of the total far--infrared emission, the starbursts always dominate in the radio. The SEDs of most of our sources are best fit with a very high optical depth of (>=50) at 1 micron. The scatter in the far-infrared/radio correlation, found among luminous IRAS sources, is due mainly to the different evolutionary status of their starburst components. The short time-scale of the star formation process amplifies the delay between the far-infrared and radio emission. This becomes more evident at low radio frequencies (below about 1 GHz) where synchrotron radiation is the dominant process. In the far-infrared (at wavelengths shorter than 100 micron) an additional source of scatter is provided by AGN, where present. AGN may be detected in the near-infrared by the absence of the knee, typical of stellar photospheres. However, near-infrared data alone cannot constrain the level at which AGN contribute because the interpretation of their observed properties, in this wave-band, depends strongly on model parameters.

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.

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]

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.

The Sunyaev-Zeldovich (SZ) effect has proven to be an extremely powerful tool to study the physical and evolutionary properties of rich clusters of galaxies. Upcoming SZ experiments, with their much improved sensitivity and angular resolution, will provide unique information also on phases of galaxy evolution characterized by the presence of large amounts of hot proto-galactic gas. We present a preliminary analysis of the SZ signals that can be expected at the collapse of the proto-galaxy when, according to the standard scenario, the gas is heated at its virial temperature, and during episodes of strong energy injections from the active nucleus. The contributions of such signals to excess power on arc-minute scales recently found by CBI and BIMA experiments are briefly discussed.

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.

We summarize our detailed, physically grounded, model for the early co-evolution of spheroidal galaxies and of active nuclei at their centers (astro-ph/0307202). Our predictions are 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. Also, the black hole mass function and the relationship between the black hole mass and the velocity dispersion in the galaxy are nicely reproduced.

We present a model for nebular emission in star forming galaxies, which takes into account the effects of dust reprocessing. The nebular emissions have been computed with CLOUDY and then included into GRASIL, our spectrophotometric code specifically developed for dusty galaxies. The interface between nebular emission and population synthesis is based on a set of pre-computed HII region emission models covering a wide range of physical quantities. Concerning the extinction properties of normal star forming galaxies, we are able to interpret the observed lack of correlation between the attenuation measured at Halpha and in the UV band as a consequence of age selective extinction. We also find that, for these galaxies with modest SFR, the ratio FIR/UV provides the best constraints on the UV attenuation. Our model also allows to deal with different SFR estimators in a consistent way, from the UV to radio wavelengths, and to discuss the uncertainties arising from the different physical conditions encountered in star forming galaxies. We provide our best estimates of SFR/luminosity calibrations, together with their expected range of variation. It results that SFR derived through Halpha, even when corrected for extinction using the Balmer decrement, is affected by important uncertainties due to age selective extinction. Another remarkable result is that SFR from UV luminosity corrected by means of the ratio FIR/UV has a small uncertainty. Finally, our model provides a calibration of SFR from radio luminosity; we are also able to reproduce the observed FIR/radio ratio.

We revisit in this paper the location of the various components observed in the AGN of NGC1068. Discrepancies between previously published studies are explained, and a new measurement for the absolute location of the K-band emission peak is provided. It is found to be consistent with the position of the central engine as derived by Gallimore (1997), Capetti (1997) and Kishimoto (1999). A series of map overlays is then presented and discussed. Model predictions of dusty tori show that the nuclear unresolved NIR-MIR emission is compatible with a broad range of models: the nuclear SED alone does not strongly constrain the torus geometry, while placing reasonable constraints on its size and thickness. The extended MIR emission observed within the ionizing cone is shown to be well explained by the presence of optically thick dust clouds exposed to the central engine radiation and having a small covering factor. Conversely, a distribution of diffuse dust particles within the ionizing cone is discarded. A simple model for the H2 and CO emission observed perpendicularly to the axis of the ionizing cone is proposed. We show that a slight tilt between the molecular disc and the Compton thick central absorber naturally reproduces the observed distribution of H2 of CO emission.

In view of the extensive evidence of tight inter-relationships between spheroidal galaxies (and galactic bulges) with massive black holes hosted at their centers, a consistent model must deal jointly with the evolution of the two components. We describe one such model, which successfully accounts for the local luminosity function of spheroidal galaxies, for their photometric and chemical properties, for 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, for 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.

Our multi-wavelength model GRASIL for the SED of galaxies is described, in particular the recent extension to the radio and X-ray range. With our model we can study different aspects of galaxy evolution by exploiting all available spectral observations, where different emission components dominate.

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.

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.

We present results from a deep mid-IR survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15um with the CAM instrument on board ISO. We found reliable optical/near-IR associations for 32 of the 35 sources detected in this field by Oliver et al. (2002, Paper I): eight of them were identified as stars, one is definitely an AGN, a second seems likely to be an AGN, too, while the remaining 22 appear to be normal spiral or starburst galaxies. Using model spectral energy distributions (SEDs) of similar galaxies, we compare methods for estimating the star formation rates (SFRs) in these objects, finding that an estimator based on integrated (3-1000um) IR luminosity reproduces the model SFRs best. Applying this estimator to model fits to the SEDs of our 22 spiral and starburst galaxies, we find that they are forming stars at rates of ~1-100 M_sol/yr, with a median value of ~40M_sol/yr, assuming an Einstein - de Sitter universe with a Hubble constant of 50 km/s/Mpc, and star formation taking place according to a Salpeter (1955) IMF across the mass range 0.1-100M_sol. We split the redshift range 0.0<z<0.6 into two equal-volume bins to compute raw estimates of the star formation rate density contributed by these sources, assuming the same cosmology and IMF as above and computing errors based on estimated uncertainties in the SFRs of individual galaxies. We compare these results with other estimates of the SFR density made with the same assumptions, showing them to be consistent with the results of Flores et al. (1999) from their ISO survey of the CFRS 1415+52 field. However, the relatively small volume of our survey means that our SFR density estimates suffer from a large sampling variance, implying that our results, by themselves, do not place tight constraints on the global mean SFR density.

We present results from a deep mid-infrared survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15 micron with the CAM instrument on board the Infrared Space Observatory (ISO). The final map in each band was constructed by the coaddition of four independent rasters, registered using bright sources securely detected in all rasters, with the absolute astrometry being defined by a radio source detected at both 7 and 15 micron. We sought detections of bright sources in a circular region of radius 2.5 arcmin at the centre of each map, in a manner that simulations indicated would produce highly reliable and complete source catalogues using simple selection criteria. Merging source lists in the two bands yielded a catalogue of 35 distinct sources, which we calibrated photometrically using photospheric models of late-type stars detected in our data. We present extragalactic source count results in both bands, and discuss the constraints they impose on models of galaxy evolution models, given the volume of space sampled by this galaxy population.

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.

In view of the extensive evidence of tight inter-relationships between spheroidal galaxies (and galactic bulges) with massive black holes hosted at their centers, a consistent model must deal jointly with the evolution of the two components. We describe one such model, which successfully accounts for the local luminosity function of spheroidal galaxies, for their photometric and chemical properties, for 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, for 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.

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.

We show that the predictions of the recent physical model by Granato et al. (2001) are fully consistent with the available estimates of the clustering properties of galaxies at (sub)-mm wavelengths. The model is then exploited to compute the expected power spectrum of fluctuations due to clustering at the frequencies of the High Frequency Instrument (HFI) of ESA's {\sc Planck} satellite. The clustering signal is found to be detectable in regions of low interstellar dust emission. A further distinctive prediction of the adopted model is a remarkably high fraction of gravitationally lensed sources at bright mm/sub-mm fluxes. In fact, according to the adopted model, only thanks to gravitational lensing, a significant number of high-$z$ forming spheroidal galaxies will be detectable by large area, shallow surveys at mm/sub-mm wavelengths, such as those carried out by P{\sc lanck}/HFI. Allowing for other source populations, we find that the fraction of gravitationally lensed sources at bright mm/sub-mm fluxes is expected to be up to $\simeq 40%$.

The clustering properties of SCUBA-selected galaxies are investigated within the framework of a unifying scheme relating the formation of QSOs and spheroids. The theoretical angular correlation function is derived for different bias functions, corresponding to different values of the ratio $M_{halo}/M_{sph}$ between the mass of the dark halos hosting such galaxies and the mass in stars produced at the end of the major star-formation burst. SCUBA sources are predicted to be strongly clustered, with a clustering strength increasing with mass. Comparisons with the best available measurements show better fits for $M_{halo}/M_{sph}\simeq 100$. The model can also account for the clustering of Lyman-break galaxies, seen as the optical counterpart of low- to intermediate-mass primeval spheroidal galaxies. Best agreement is once again obtained for high values of the $M_{halo}/M_{sph}$ ratio. We also discuss implications for small scale fluctuations observed at different wavelengths by forthcoming experiments such as the Planck mission aimed at mapping the Cosmic Microwave Background (CMB).

The intensity of the Cosmic Far-IR Background and the strong evolution of galaxies in the far-IR to mm wavelength range demonstrate that the bulk of starlight emitted during the early phases of galaxy evolution was reprocessed by dust. Therefore, the optical view of the galaxy formation process is highly incomplete and biased, and must be complemented with far-IR/sub-mm observations. We review the impact of sub-mm surveys on our understanding of the evolutionary history of spheroidal galaxies. A recent model, bringing into play also the inter-relationships between formation and early evolution of spheroidal galaxies and quasars, is described and some implications are outlined.

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.

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.

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.

We examine the tight correlation observed between the far infrared (FIR) and radio emission in normal and star burst galaxies. We show that significant deviations from the average relation are to be expected in young star burst or post star burst galaxies, due to the different fading times of FIR and Radio power.

This paper investigates the clustering properties of SCUBA-selected galaxies in the framework of a unifying scheme relating the formation of QSOs and spheroids (Granato et al. 2000). The theoretical angular correlation function is derived for different bias functions, corresponding to different values of the ratio $M_{\rm halo}/M_{\rm sph}$ between the mass of the dark halo and the final mass in stars. SCUBA sources are predicted to be strongly clustered, with a clustering strength increasing with mass. We show that the model accounts for the clustering of Lyman-break Galaxies, seen as the optical counterpart of low- to intermediate-mass primeval spheroidal galaxies and is also consistent with the observed angular correlation function of Extremely Red Objects. Best agreement is obtained for $M_{\rm halo}/M_{\rm sph}=100$. We also consider the implications for small scale fluctuations observed at sub-mm wavelengths by current or forthcoming experiments aimed at mapping the Cosmic Microwave Background (CMB). The predicted amplitude of the clustering signal in the 350 GHz channel of the Planck mission strongly depends on the halo-to-bulge mass ratio and may be of comparable amplitude to primary CMB anisotropies for multipole numbers $l\simgt 50$.

We present results about the distribution and kinematics of the molecular environment of the AGN in NGC1068, over a 1.5'' x 3.5'' region around the central engine in NGC1068, derived from H2 line emission detected with ISAAC at VLT/ANTU on ESO/Paranal. The H2 emitting molecular gas is found to be distributed along the East-West direction and with two main peak emission (knots) located at a distance of about 70 pc from the central engine. The eastern H2 knot is more intense than the western one. The line profiles mapped across the entire 1.5'' x 3.5'' region, at a spatial resolution of 0.3'' x 0.45'', appear to be quite complex with either a blue or red wing. At first order, we find a velocity difference of 140 km/s between the two knots; if interpreted as quasi-keplerian velocity, this implies a central enclosed mass of 10^8 solar masses.

Direct and indirect observational evidence leads to the conclusion that high redshift QSOs did shine in the core of early type proto-galaxies during their main episode of star formation. Exploting this fact, we derive the rate of formation of this kind of stellar systems at high redshift by using the QSO Luminosity Function. The elemental proportions in elliptical galaxies, the descendents of the QSO hosts, suggest that the star formation was more rapid in more massive objects. We show that this is expected to occur in Dark Matter haloes, when the processes of cooling and heating is considered. This is also confirmed by comparing the observed sub-mm counts to those derived by coupling the formation rate and the star formation rate of the spheroidal galaxies with a detailed model for their SED evolution. In this scenario SCUBA galaxies and Lyman Break Galaxies are early type proto-galaxies forming the bulk of their stars before the onset of QSO activity.

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.

We analyze a morphologically-selected complete sample of 52 late-type (spiral and irregular) galaxies in the Hubble Deep Field North with total K-magnitudes brighter than K=20.47. This sample exploits in particular the ultimate imaging quality achieved by HST in this field, allowing us to clearly disentangle the galaxy morphologies, based on accurate profiles of the surface brightness distributions. Our purpose was to investigate systematic differences between the two classes, as for colours, redshift distributions and ages of the dominant stellar populations. This sample appears to miss significantly galaxies above z=1.4 (in a similar way as an early-type galaxy sample previously studied by us), a fact which may be explained as a global decline of the underlying mass function for galaxies at these high redshifts. Differences between early and late-types are apparent -particularly in the colour distributions and the evolutionary star-formation (SF) rates per unit volume-, although the complication in spectro-photometric modelling introduced by dust-extinction in the gas-rich systems prevents us to reach conclusive results on the single sources. However, we find that an integrated quantity like the comoving star-formation rate density (SFR(t)) as a function of redshift is much less affected by these uncertainties: by combining this with the previously studied early-type galaxy sample, we find a shallower dependence of SFR(t) on z between z=0.2 and z=1.5 than found by Lilly et al. (1995).

Recently (Granato, Lacey, Silva et al. 2000, astro-ph/0001308) we have combined our spectrophotometric galaxy evolution code which includes dust reprocessing (GRASIL, Silva et al. 1998) with semi-analytical galaxy formation models (GALFORM, Cole et al. 1999). One of the most characteristic features of the former is that the dust is divided in two main phases: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. As a consequence, stellar populations of different ages have different geometrical relationships with the two phases, which is essential in understanding several observed properties of galaxies, in particular those undergoing major episodes of star formation at any redshift. Indeed, our merged GRASIL+GALFORM model reproduces fairly well the SEDs of normal spirals and starbursts from the far-UV to the sub-mm and their internal extinction properties. In particular in the model the observed starburst attenuation law (Calzetti 1999) is accounted for as an effect of geometry of stars and dust, and has nothing to do with the optical properties of dust grains.

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.

We present a new evolutionary model for the far-UV to sub-mm properties of the galaxy population. This combines a semi-analytic galaxy formation model based on hierarchical clustering (GALFORM) with a spectro-photometric code which includes dust reprocessing (GRASIL). The former provides the star formation and metal enrichment histories, together with the gas mass and various geometrical parameters, for a representative sample of galaxies formed in different density environments. These quantities allow us to model the SEDs of galaxies, taking into account stellar emission and also dust extinction and re-emission. Two phases are considered for the dust: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. The model includes both galaxies forming stars quiescently in disks, and starbursts triggered by galaxy mergers. We test our models against the observed spectro-photometric properties of galaxies in the local Universe. The models reproduce fairly well the SEDs of normal spirals and starbursts, and their internal extinction properties. The starbursts follow the observed relationship between the FIR to UV luminosity ratio and the slope of the UV continuum. They also reproduce the observed starburst attenuation law (Calzetti et al 99). This result is remarkable, because we use a dust mixture which reproduces the Milky Way extinction law. It suggests that the observed attenuation law is related to the geometry of the stars and dust. We compute galaxy luminosity functions over our wide range of wavelengths, which turn out to be in good agreement with observational data. The UV continuum turns out to be a poor star formation indicator for our models, whilst the infrared luminosity is much more reliable.

We have explored, in the general framework of the ``young source'' scenario, evolutionary models for GHz Peaked Spectrum (GPS) galaxies and quasars which reproduce the observed counts, redshift and peak frequency distributions of currently available samples. Substantially different cosmological evolution properties are found for the two populations: the quasar luminosity function must evolve strongly up to $z\sim 1$, while the data on galaxies may be consistent with no evolution. The models show that GPS sources (mostly quasars) may comprise quite a significant fraction of bright ($S> 1 $Jy) radio sources at $\nu \geq 30 $GHz if the intrinsic distribution of peak frequencies extends up to $\sim 1000 $GHz. In any case, however, their fraction decreases rapidly with decreasing flux and their contribution to small scale fluctuations in the frequency range covered by the forthcoming space missions MAP and Planck Surveyor is expected to be minor.

We studied with BeppoSAX the infrared luminous galaxy IRAS 09104+4109 over a very wide X-ray band from 0.1 to 80 keV. Our observations indicate the dominance of a thermal component at energies below 8 keV, which we attribute to the free-free emission from the intracluster (IC) plasma surrounding the source. Above 10 keV we find evidence for the existence of flux in excess with respect to the free-free IC plasma emission. This, together with the marginal detection of a neutral iron line at $\sim 6.4$ keV, gives a strong indication for the presence of an AGN deeply buried within the source. This component is best modelled by a strongly absorbed (N_H > 5 10^{24} cm^{-2}) power-law plus unabsorbed reflection spectrum (0.15<R<0.3), which is also responsible for the cold iron line. The unabsorbed broad-band (2-100 keV) X-ray emission of this AGN is 2.5 10$^{46} erg/s, well within the range of quasar luminosities. Our results indicate that IRAS 09104+4109 is indeed the prototype of a rare class of sources, the luminous type-2 QSOs. The association of this source with a huge cooling-flow of $\sim 1000 M_\odot$ in the cluster, as indicated by the X-ray data, might suggest that such condition of extremely fast mass accumulation could favour the survival of a thick obscuring envelope, which would otherwise be quickly destroyed by the very luminous central source.

We present the basic features and preliminary results of the interface between our spectro-photometric model GRASIL (that calculates galactic SED from the UV to the sub-mm with a detailed computation of dust extinction and thermal reemission) with the semi-analytical galaxy formation model GALFORM (that computes galaxy formation and evolution in the hierarchical scenario, providing the star formation history as an input to our model). With these two models we are able to synthetize simulated samples of a few thousands galaxies suited for statistical studies of galaxy properties to investigate on galaxy formation and evolution. We find good agreement with the available data of SED and luminosity functions.

We summarize our modelling of galaxy photometric evolution (GRASIL code). By including the effects of dust grains and PAHs molecules in a two phases clumpy medium, where clumps are associated to star forming regions, we reproduce the observed UV to radio SEDs of galaxies with star formation rates from zero to several hundreds $M_\odot/yr$. GRASIL is a powerful tool to investigate the star formation, the initial mass function and the 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 six decades. It may be interfaced with any device providing the star formation and metallicity histories of a galaxy. As an application, we have investigated the properties of early-type galaxies in the HDF, tracking the contribution of this population to the cosmic star-formation history, which has a broad peak between z=1.5 and 4. To explain the absence of objects at z $\gsim$ 1.3, we suggest a sequence of dust-enshrouded merging-driven starbursts in the first few Gyrs of galaxies lifetime. We are at present working on a complementary sample of late type objects selected in a similar way.

We test emission models of circum-nuclear dust torii around quasars, at low and high redshifts, by using a large collection of photometric data for an unbiased sample of 120 optically-selected objects with millimetric and sub-millimetric fluxes, including new unpublished data. Under the assumption that the dust is heated by a point-like source with a power-law primary spectrum, as defined by the observed optical-UV continuum, we infer the basic model parameters, such as dust masses, temperature distributions and torus sizes, by numerically solving the radiative transfer equation in the dust distribution. In addition to the substantiated statistics, an essential improvement over previous analyses comes from the use of optical-UV data to constrain the primary illuminating continuum, which is needed to estimate dust temperatures and sizes. Dependences of the best-fit parameters on luminosity and redshift are studied and the contribution of dust in the host galaxy to the observed fluxes is briefly mentioned. This analysis constrains the properties of the enriched interstellar medium in the galaxies hosting the quasars. The dust abundance does not display appreciable trends as a function of redshift, from z ~ 1 to almost 5, and shows that dust and metals are at least as, and often more, abundant at these early epochs than they are in local galactic counterparts. This evolutionary pattern is remarkably at variance with respect to what is expected for disk galaxies, like the Milky Way, slowly building metals during the whole Hubble time. It rather points in favour of a much more active phase of star-formation at early epochs, probably related to the formation of massive spheroidal galaxies.

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.

We have investigated the properties of a complete K-band selected sample of 35 elliptical and S0 galaxies brighter than K=20.15 in the Hubble Deep Field, as representative of the field galaxy population. This sample has been derived from deep K-band image by the KPNO-IRIM camera, by applying a rigorous morphological classification scheme based on quantitative analyses of the surface brightness profiles. The broad-band spectra of the sample galaxies allow us to date their dominant stellar populations. The majority of bright early-types in this field are found at redshifts z<1.3 to display colors indicative of a fairly wide range of ages (typically 1.5 to 3 Gyrs). We find that the major episodes of star-formation building up typical M_star galaxies have taken place during a wide redshift interval 1<z<4 for q_0=0.5, which becomes 1<z<3 for q_0=0.15. Our estimated galactic masses, for a Salpeter IMF, are found in the range from a few billion to a few 10^{11} solar masses already at z=1. So the bright end of the E/S0 population is mostly in place by that cosmic epoch, with space densities, masses and luminosities consistent with those of the local field E/S0 population. What distinguishes the present sample is a remarkable absence of objects at z>1.3, which should be detectable during the luminous star-formation phase expected to happen at these redshifts. Obvious solutions are a) that the merging events imply perturbed morphologies which prevent selecting them by our morphological classification filter, or b) that a dust-polluted ISM obscures the (either continuous or episodic) events of star-formation. We conclude that the likely solution is a combination thereof, i.e. a set of dust-enshrouded merging-driven starbursts occurring during the first few Gyrs of the galaxy's lifetime.

We review estimates of small scale fluctuations due to extragalactic point sources in the Planck Surveyor frequency bands. While our undestanding of the spectral and evolutionary properties of these sources is far from complete, conservative estimates allow us to confidently conclude that, in the frequency range 100--200 GHz, their contaminating effect is well below the expected anisotropy level of the cosmic microwave background (CMB), down to angular scales of at least $\simeq 10'$. Hence, an accurate subtraction of foreground fluctuations is not critical for the determination of the CMB power spectrum up to multipoles $\ell \simeq 1000$. In any case, Planck's wide frequency coverage will allow to carefully control foreground contributions.On the other hand, the all sky surveys at 9 frequencies, spanning the range 30--900 GHz, will be unique in providing complete samples comprising from several hundreds to many thousands of extragalactic sources, selected in an essentially unexplored frequency region. New classes of sources may be revealed in these data. The familiar ``flat''-spectrum radio sources should show spectral features carrying essential information on their physical properties. Crucial information will be provided to understand the nature of radio sources with strongly inverted spectra. Scenarios for the cosmological evolution of galaxies will be extensively tested.

For an extended sample of Seyfert galaxies we compile from the literature the infrared fluxes in the four IRAS bands, the ground-based small-beam (~ 5-10") fluxes in the standard Q, N, M, L (or L') bands, and the nuclear (non-stellar) estimated fluxes in the JHK bands. We estimate nuclear fluxes in the L band by applying a correction for stellar light. From the statistical study of the infrared colors and luminosities, we derive the typical SEDs of Seyfert 1 and 2 nuclei and the typical differences in luminosities between the two types of objects in the mid- and near-infrared spectral ranges. The observational data appear to severely challenge many models of dusty tori, which hardly account for the shapes of the SEDs and the degree of anisotropy observed in Seyfert galaxies. In particular, at variance with some earlier claims, very thick and compact tori are basically inconsistent with these observational constraints. The most successful models, though having problems in accounting for several details, can fit the major infrared observational data of both Seyfert 1 and Seyfert 2 nuclei with tori which extend up to several hundreds pc and have fairly low optical thickness.

We have investigated the spectro-photometric properties of the Asymptotic Giant Branch (AGB) stars and their contribution to the integrated infrared emission in simple stellar populations (SSP). Adopting analytical relations describing the evolution of these stars in the HR diagram and empirical relations for the mass-loss rate and the wind terminal velocity, we were able to model the effects of the dusty envelope around these stars, with a minimal number of parameters. We computed isochrones at different age and initial metal content. We compare our models with existing infrared colors of M giants and Mira stars and with IRAS PSC data. Contrary to previous models, in the new isochrones the mass-loss rate, which establishes the duration of the AGB phase, also determines the spectral properties of the stars. The contribution of these stars to the integrated light of the population is thus obtained in a consistent way. We find that the emission in the mid infrared is about one order of magnitude larger when dust is taken into account in an intermediate age population, irrespective of the particular mixture adopted. The dependence of the integrated colors on the metallicity and age is discussed, with particular emphasis on the problem of age-metallicity degeneracy. We show that, contrary to the case of optical or near infrared colors, the adoption of a suitable pass-band in the mid infrared allows a fair separation of the two effects. We suggest intermediate redshift elliptical galaxies as possible targets of this method of solving the age-metallicity dilemma. The new SSP models constitute a first step in a more extended study aimed at modelling the spectral properties of the galaxies from the ultraviolet to the far infrared.

Our present understanding of the extragalactic source counts and background radiation at infrared and sub-mm wavelengths is reviewed. Available count data are used to constrain evolutionary models of galaxies and Active Nuclei. The CIRB, on the other hand, provides crucial information on the integrated past IR emissivity, including sources so faint to be never accessible. The near-IR and sub-mm are suited to search for the CIRB. These cosmological windows are ideal to detect redshifted photons from the two most prominent broad emission features in galaxy spectra: the stellar photospheric 1 micron peak and the one at 100 micron due to dust re-radiation. The recently claimed detection of an isotropic diffuse component at 100-200 microns would support the evidence for strong cosmic evolution from $z=0$ to $z\sim 2$ of faint gas-rich late-type galaxies, as inferred from direct long-$\lambda$ counts. Furthermore, an equally intense CIRB flux measured in the wavelength range 200 to 500 $\mu m$ may be in support of models envisaging a dust-enshrouded phase during formation of spheroidal galaxies (e.g. Franceschini et al., 1994): the background spectrum at such long wavelengths would indicate an high redshift ($z>2-3$) for this active phase of star formation, implying a relevant constraint on structure formation scenarios. We argue that, given the difficulties for sky surveys in this spectral domain, only sporadic tests of these ideas will be achieved in the next several years. Only a mission like FIRST, in combination with large mm arrays and optical telescopes, will allow to thoroughly investigate, via photometric and spectroscopic surveys, these IR-bright early phases in the evolution of galaxies and AGNs.

We present the first results of a detailed modeling of chemical and photometric evolution of galaxies including the effects of a dusty interstellar medium. A chemical evolution code follows the SF rate, the gas fraction and the metallicity, basic ingredients for the stellar population synthesis. The latter is performed with a grid of integrated spectra for SSP of different ages and metallicities, in which the effects of dusty envelopes around AGB stellar are included. The residual fraction of gas in the galaxy is divided into two phases: the star forming molecular clouds and the diffuse medium. The molecular gas is sub-divided into clouds of given mass and radius: it is supposed that any SSP is borne within the cloud and progressively escapes it. The emitted spectrum of the star forming molecular clouds is computed with a radiative transfer code. The diffuse dust emission (cirrus) 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 PAH are included. The model is compared and calibrated with available data of starburst galaxies in the local universe, in particular new broad-band and spectroscopic ISO observations. It will be a fundamental tool to extract information on global evolution of the stellar component and of ISM of galaxies from observations covering four decades in $\lambda$. The combination of observations from HST, Keck, ISO and ground-based optical, IR and sub-mm telescopes already provided several objects at relevant $z$ with spectral information on this large $\lambda$ range. Their number will increase after the completion of ISO surveys, and will burst when SIRFT, FIRST, Planck

We have investigated the properties of a bright K-band selected sample of early-type galaxies in the Hubble Deep Field, as representative of the field galaxy population. This dataset is unique as for the morphological information on faint high-z sources, and for complete photometric and spectroscopic coverage. The majority of bright early-type galaxies in this field are found at redshifts $z \leq 1.3$ to share common properties with those of high-z cluster samples, as for the age and mass of the dominant stellar population -- which are found to be as old as 3-5 Gyr and as massive as $10^{11} M_\odot$ already at $z\simeq 1$. There is no evidence of a steepening of the mass function from $z$=0 to $z$=1, as inferred by some authors from analyses of optically-selected samples and favoured by hierarchical clustering models forming most of the E/S0s at $z<1$. What distinguishes this sample is a remarkable absence of objects at $z>1.3$, which would be expected as clearly detectable above the flux limits, given the aged properties of the lower redshift counterparts. So, something hide them at high redshifts. Merging could be an explanation, but it would require that already formed and old stellar systems would assemble on short timescales of 1 Gyr or less. Only a negligible fraction of the galactic mass in young stars has to be added, a rather contrived situation from the dynamical point of view. An alternative interpretation could be that a dust-polluted ISM obscures the first 3 to 4 Gyr of the evolution of these massive systems, after which a galactic wind and gas consumption makes them transparent. The presence of dust would have relevant (and testable) implications for the evaluation of the global energetics from galaxy formation and for the visibility of the early phases by current and future infrared

We use photometric data, from the optical to the mm for a large sample of optically selected radio-quiet quasars, at low and high redshifts, to test emission models from circum-nuclear dusty torii around them. Model parameters, such as dust mass, temperature distributions and torus sizes are inferred, under the assumption that the dust is heated only by the central nuclear source. Dependences of best-fit parameters on luminosity and redshift are studied and the contribution of dust in the host galaxy to the observed fluxes is briefly mentioned. Only tentative conclusions can be drawn from this yet large dataset and detailed modelling, because crucial spectral information is missing at wavelengths where the QSO dust emission is more prominent, i.e. at $\lambda \simeq 50$ to 600 $\mu m$. Poor information in this wavelength range affects our understanding not only of the high-redshift objects, but also of the low-redshift sample. We then expect fundamental contributions by the FIRST mission in this field, which is obviously related to the more general one of galaxy formation.

Two families of models of dusty tori in AGNs (moderately thick and extended versus very thick and compact) are tested against available observations. The confrontation suggests that the former class better explains the IR broad-band spectra of both broad and narrow line AGNs, the anisotropy of the emission deduced by comparing IR properties of Seyfert 1 and 2 nuclei, the results of IR spectroscopy and those of high spatial resolution observations. There is however clear evidence for a broad distribution of optical depths. We also examine the relationship between IR and X-ray emission. The data support a view in which the matter responsible for the X-ray absorption is mostly dust free, lying inside the dust sublimation radius. The consequences of these results for the hard X-ray background as well as IR counts and background are discussed.

We investigate models for the power supply and broad-band spectral energy distribution (SED) of hyperluminous IR galaxies, recently discovered at high redshifts, in terms of the emission from an active nucleus embedded in a torus-like dusty structure. We find consistent solutions in terms of a simple torus model extended several hundreds of parsecs, with $A_V$ in the equatorial plane of a few hundreds and a typical covering factor of over $50\%$. Objects as different as the prototype high-z galaxy \fd, the z=0.93 IR object \fq, \fz\ found in a high-z cooling flow, and the optically selected BAL "Cloverleaf" quasar, are all fitted by the same solution, for decreasing values of the polar angle to the line-of-sight and proper scaling of the luminosities. We suggest that such luminous high-z IR objects are heavily buried quasars surrounded by large amounts of dust with high covering factors and large optical depths. Comparison with UVX QSOs suggests that they are observed during a transient phase. Forthcoming observations in the far-IR will soon allow probing this phase and its relationship with the -- possibly concomitant -- formation of the nuclear black hole and of the host galaxy.