We present a multi-wavelength study (from X-ray to millimetre) of the nearby low-luminosity active galactic nucleus (LLAGN) NGC 7213. We combine the information from the different bands to characterize the source in terms of contribution from the AGN and the host-galaxy interstellar medium (ISM). This approach allows us to provide a coherent picture of the role of the AGN and its impact, if any, on the star formation and molecular gas properties of the host galaxy. We focused our study on archival ALMA Cycle 1 observations, where the CO(2-1) emission line has been used as a tracer of the molecular gas. Using the 3DBarolo code on ALMA data, we performed the modelling of the molecular gas kinematics traced by the CO(2-1) emission, finding a rotationally dominated pattern. The host-galaxy molecular gas mass was estimated from the integrated CO(2-1) emission line obtained with APEX data, assuming an $\alpha_{CO}$ conversion factor. By using the ALMA data, we would have underestimated the gas masses by a factor $\sim$3, given the filtering out of the large scale emission in interferometric observations. We also performed a complete X-ray spectral analysis on archival observations, revealing a relatively faint and unobscured AGN. The AGN results to be too faint to significantly affect the properties of the host-galaxy, such as star formation activity and molecular gas kinematics and distribution.

We present the first [CII] 158 $\mu$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies.

We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim 4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A} < -20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [C II] detections and 43 are upper limits. With this UV-selected sample, we set stringent lower bounds to the volume density of [C II] emitters at $z\sim4-6$. Our derived LFs are consistent with the $z\sim0$ [C II] LF at $\leq 10^9L_\odot$, but show a hint of an excess of [C II] emitters at the luminous end ($>10^9L_\odot$) compared to $z\sim0$. This excess is consistent with a potential population of [C II]-bright UV-faint sources that are largely excluded from the ALPINE target selection. To further constrain such sources, in a companion analysis, Loiacono et al. (2020) estimate the [C II] LF based on serendipitous [C II] sources in the ALPINE maps. Their resulting single data point is consistent with the excess we observe at high-L. Combined, these studies suggest possible redshift evolution in the [C II] LFs at high luminosity between $z\sim 5$ and $0$. This conclusion is further supported by the inferred [C II] LFs from the far-IR and CO LFs at $z\sim4-6$. Integrating our estimate on the best-fit [C II] LF, we set a constraint on molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$ Mpc$^{-3}$. This is broadly consistent with previous studies, although the uncertainties span more than an order of magnitude. Finally, we find that available model predictions tend to underestimate the number densities of [C II] emitters at $z\sim5$.

We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA Band-7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 16. per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (7 show a counterpart only in the mid-IR and no HST or near-IR identification, while 2 are detected as [CII] emitters at z=5). The 7 HST and near-IR dark galaxies with mid-IR counterpart contribute for about 15 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3.

The ALPINE-ALMA large program targets the [CII] 158 $\mu$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates, we measured the conversion factor from the ALMA 158 $\mu$m rest-frame dust continuum luminosity to the total infrared luminosity (L$_{\rm IR}$) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L$_{\rm IR}$ of 4.4$\times 10^{11}$ L$_\odot$. We also detected 57 additional continuum sources in our ALMA pointings. They are at lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5$\pm$0.2. We measured the 850 $\mu$m number counts between 0.35 and 3.5 mJy, improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40 % of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8$\times$10$^8$ L$_\odot$ and their median full width at half maximum is 252 km/s. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR-L$_{\rm [CII]}$ relations of De Looze et al. (2014) and Lagache et al. (2018).

We present ALMA observations of a merging system at z ~ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII] 158 micron and far-infrared continuum data with multi-wavelength ancillary data we find that the system is composed of two massive (Mstar >~ 10^10 Msun) star-forming galaxies experiencing a major merger (stellar mass ratio r_mass ~ 0.9) at close spatial (~13 kpc; projected) and velocity (delta_v < 300 km/s) separations, and two additional faint narrow [CII]-emitting satellites. The overall system belongs to a larger-scale protocluster environment and is coincident to one of its overdensity peaks. ALMA reveals also the presence of [CII] emission arising from a circumgalactic gas structure, extending out to ~30 kpc. Our morpho-spectral decomposition analysis shows that about 50% of the total flux resides between the individual galaxy components, in a metal-enriched gaseous envelope characterized by a disturbed morphology and complex kinematics. Similarly to observations of shock-excited [CII] emitted from tidal tails in local groups, our results can be interpreted as a possible signature of interstellar gas stripped by strong gravitational interactions, with a possible contribution from material ejected by galactic outflows and emission triggered by star formation in small faint satellites. Our findings suggest that mergers could be an efficient mechanism of gas mixing in the circumgalactic medium around high-z galaxies, and thus play a key role in the galaxy baryon cycle at early epochs.

The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z<2. We discuss three possible effects, namely outflows, halt of gas supplying, and over-efficient star formation, which may jointly contribute to the gas fraction plateau of the latter massive galaxies.

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at redshift ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame $158\,\rm{\mu m}$ obtained with the ALMA large program ALPINE. We study the connection between the UV spectral slope ($\beta$), stellar mass ($M_{\star}$), and infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$). Twenty-three galaxies are individually detected in the continuum at $>3.5\,\sigma$ significance. We perform a stacking analysis using both detections and non-detections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$\beta$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$\beta$ relation as a function of redshift suggests an evolution of dust attenuation properties at $z>4$. Similarly, we find that our galaxies have lower IRX values up to 1 dex on average at fixed mass compared to previously studied IRX-$M_{\star}$ relations at $z\lesssim4$, albeit with significant scatter. This implies a lower obscured fraction of star-formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at $z\gtrsim4$ are characterised by (i) a steeper attenuation curve than at $z\lesssim4$, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies ($\log M_{\star}/M_\odot > 10$) at z~5-6 already exhibit an obscured fraction of star formation of $\sim45\%$, indicating a rapid build-up of dust during the epoch of reionization.

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$\sigma$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($\beta$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($\Delta\beta$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$\lambda$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows.

We report the serendipitous discovery of a bright galaxy (Gal-A) observed as part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). While this galaxy is detected both in line and continuum emission in ALMA Band 7, it is completely dark in UV/optical filters and only presents a marginal detection in the UltraVISTA Ks band. We discuss the nature of the observed ALMA line, i.e. whether the emission comes from [CII] at z~4.6, or from high-J CO transitions at z~2.2. In the first case we find a [CII]-to-FIR luminosity ratio of log(L_[CII]/L_FIR)=-2.5, consistent with the average value for local star-forming galaxies (SFGs); in the second case, instead, the source would lie outside of the empirical relations between L_CO and L_FIR found in the literature. At both redshifts, we derive the star-formation rate (SFR) from the ALMA continuum, and the stellar mass (M*) by using stellar population synthesis models as input for LePHARE spectral energy distribution (SED) fitting. Exploiting our results, we believe that Gal-A is a Main-Sequence (MS), dusty SFG at z=4.6 (i.e. [CII] emitter) with log(SFR/[M/yr])~1.4 and log(M*/M)~9.7. This work underlines the crucial role of the ALPINE survey in making a census of this class of objects, in order to unveil their contribution to the global star-formation rate density (SFRD) of the Universe at the end of the Reionisation epoch.

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

We study the impact of active galactic nuclei (AGN) on the CO Spectral Line Energy Distribution (SLED) of high-$z$ galaxies. In particular, we want to assess if the CO SLED can be used as a probe of AGN activity. To this purpose, we develop a semi-analytical model that takes into account the radiative transfer and the clumpy structure of giant molecular clouds where the CO lines are excited, their distribution in the galaxy disk, and the torus obscuration of the AGN radiation. We study the joint effect on the CO SLED excitation of (i) the X-ray luminosity from the AGN ($L_{X}$), (ii) the size of the molecular disk, (iii) the inclination angle ($\Omega$) of the torus with respect to the molecular disk, and (iv) the GMC mean density. We also discuss the possibility of an enhanced Cosmic Ray Ionization Rate (CRIR). We find that the X-ray Dominated Region (XDR) generated by the AGN in every case enhances the CO SLED for $J>5$, with increasing excitation of high-$J$ CO lines for increasing X-ray luminosity. Because high-$z$ galaxies are compact, the XDR region typically encloses the whole disk, thus its effect can be more important with respect to lower redshift objects. The impact of the XDR can be disentangled from an enhanced CRIR either if $L_X>10^{44} \rm \,erg\, s^{-1}$, or if $\Omega \geq 60^{\circ}$. We finally provide predictions on the CO(7-6)/[CII] and CO(17-16)/[CII] ratios as a function of $L_X$, which can be relevant for ALMA follow up of galaxies and quasars previously detected in [CII].

We derive for the first time the dust mass function (DMF) in a wide redshift range, from z~0.2 up to z~2.5. In order to trace the dust emission, we start from a far-IR (160-um) Herschel selected catalogue in the COSMOS field. We estimate the dust masses by fitting the far-IR data (lam_rest>50um) with a modified black body function and we present a detailed analysis to take into account the incompleteness in dust masses from a far-IR perspective. By parametrizing the observed DMF with a Schechter function in the redshift range 0.1<z<0.25, where we are able to sample faint dust masses, we measure a steep slope (alpha~1.48), as found by the majority of works in the Local Universe. We detect a strong dust mass evolution, with M_d^star at z~2.5 almost one dex larger than in the local Universe, combined with a decrease in their number density. Integrating our DMFs we estimate the dust mass density (DMD), finding a broad peak at z~1, with a decrease by a factor of ~3 towards z~0 and z~2.5. In general, the trend found for the DMD mostly agrees with the derivation of Driver et al. (2018), another DMD determination based also on far-IR detections, and with other measures based on indirect tracers.

We report the detection of [CII]158um emission from a system of three closely-separated sources in the COSMOS field at z~4.56, as part of the ALMA Large Program to INvestigate CII at Early times (ALPINE). The two dominant sources are closely associated, both spatially (1.6"~11kpc) and in velocity (~100km/s), while the third source is slightly more distant (2.8"~18kpc, ~300km/s). The second strongest source features a slight velocity gradient, while no significant velocity gradient is seen in the other two sources. Using the observed [CII] luminosities, we derive a total log(SFR_[CII]/[Msol/year])=2.8+/-0.2, which may be split into contributions of 59%, 31%, and 10% from the central, east, and west sources, respectively. Comparison of these [CII] detections to recent zoom-in cosmological simulations suggests an ongoing major merger. We are thus witnessing a system in a major phase of mass build-up by merging, including an on-going major merger and an upcoming minor merger, which is expected to end up in a single massive galaxy by z~2.5.

The estimate of the number and space density of obscured AGN over cosmic time still represents an open issue. While the obscured AGN population is a key ingredient of the X-ray background synthesis models and is needed to reproduce its shape, a complete census of obscured AGN is still missing. Here we test the selection of obscured sources among the local 12-micron sample of Seyfert galaxies. Our selection is based on a difference up to three orders of magnitude in the ratio between the AGN bolometric luminosity, derived from the spectral energy distribution (SED) decomposition, and the same quantity obtained by the published XMM-Newton 2-10 keV luminosity. The selected sources are UGC05101, NGC1194 and NGC3079 for which the available X-ray wide bandpass, from Chandra and XMM-Newton plus NuSTAR data, extending to energies up to ~30-45 keV, allows us an accurate determination of the column density, and hence of the true intrinsic power. The newly derived NH values clearly indicate heavy obscuration (about 1.2, 2.1 and 2.4 x10^{24} cm-2 for UGC05101, NGC1194 and NGC3079, respectively) and are consistent with the prominent silicate absorption feature observed in the Spitzer-IRS spectra of these sources (at 9.7 micron rest frame). We finally checked that the resulting X-ray luminosities in the 2-10 keV band are in good agreement with those derived from the mid-IR band through empirical L_MIR-L_X relations.

We investigate a sample of 152 dusty sources at 1:5 < z < 2:5 to understand the connection of enhanced Star-Formation-Rate (SFR) and Black-Hole-Accretion-Rate (BHAR). The sources are Herschel-selected, having stellar masses M*>10^10Msun and SFR (100-1000Msun/yr) elevated(>4?) above the star-forming "main sequence", classifying them as Starbursts (SB). Through a multiwavelength fitting approach (including a dusty torus component), we divided the sample into active SBs (dominated by an AGN emission, SBs-AGN, ? 23% of the sample) and purely star-forming SBs (SBs-SFR). We visually inspected their HST/UV-restframe maps: SBs-SFR are generally irregular and composite systems; ? 50% of SBs-AGN are instead dominated by regular compact morphologies. We then found archival ALMA continuum counterparts for 33 galaxies (12 SBs-AGN and 21 SBs-SFR). For these sources we computed dust masses, and, with standard assumptions, we also guessed total molecular gas-masses. SBs turn to be gas rich systems (fgas = Mgas=Mgas/(Mgas+M*) \sim 20%-70%), and the gas fractions of the two SB classes are very similar (fgas = 43 +/-4% and fgas = 42+/- 2%). Our results show that SBs are consistent with a mixture of: 1) highly star-forming merging systems (dominating the SBs-SFR), and 2) primordial galaxies, rapidly growing their M* together with their Black Hole (mainly the more compact SBs-AGN). Anyway, feedback effects have not reduced their fgas yet. Indeed, SBs at z = 2, with relatively low bolometric AGN luminosities in the range 10^44 < Lbol(AGN) < 10^46 erg/s (compared to bright optical and X-ray quasars), are still relatively far from the epoch when the AGN feedback will quench the SFR in the host and will substantially depress the gas fractions.

We aim to study the statistical properties of dusty star-forming galaxies, such as their number counts, luminosity functions (LF) and dust-obscured star-formation rate density (SFRD). We use state-of-the-art de-blended Herschel catalogue in the COSMOS field, generated by combining the Bayesian source extraction tool XID+ and informative prior on the spectral energy distributions, to measure the number counts and LFs at far-infrared (FIR) and sub-millimetre (sub-mm) wavelengths. Thanks to our de-confusion technique and deep multi-wavelength photometry, we are able to achieve more accurate measurements while probing ten times below the confusion limit. Our number counts at 250 microns agree well with previous Herschel studies. However, our counts at 350 and 500 microns are considerably below previous Herschel results. This is due to previous studies suffering from source confusion which worsens towards longer wavelength. Our number counts at 450 and 870 microns show excellent agreement with previous determinations derived from single dish and interferometric observations. Our measurements of the monochromatic LF and the total IR LF agree well with previous results. The increased dynamic range of our measurements allows us to better measure the faint-end slope of the LF and measure the dust-obscured SFRD out to z~6. We find that the fraction of dust obscured star-formation activity is at its highest around z~1 which then decreases towards both low and high redshift. We do not find a shift of balance between z~3 and z~4 in the cosmic star-formation history from being dominated by unobscured star formation at higher redshift to obscured star formation at lower redshift. However, we do find 3<z<4 to be an interesting transition period as the fraction of the total SFRD that is obscured by dust is significantly lower at higher redshifts.

Recent derivations of the galaxy star formation rate density (SFRD) obtained from sub-millimetre (sub-mm) surveys (e.g., SCUBA-2) show a tension with previous works based on Herschel and multi-wavelength data. Some of these works claim that the SFRD derived by pushing the Herschel surveys beyond z~2 are incorrect. However, the current sub-mm surveys obtained from SCUBA-2 data and the methods used to construct the total infrared (IR) luminosity function (LF) and the SFRD could be affected by some limitations. Here we show how these limitations (i.e., selection bias and incompleteness effects) might affect the total IR LF, making the resulting dusty galaxy evolution of difficult interpretation. In particular, we find that the assumed spectral energy distribution (SED) plays a crucial role in the total IR LF derivation, moreover, we confirm that the long-wavelength (e.g., 850-micron) surveys can be incomplete against "warm" SED galaxies, and that the use of a wide spectral coverage of IR wavelengths is crucial to limit the uncertainties and biases.

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

We use two catalogues, a Herschel catalogue selected at 500 mu (HerMES) and an IRAS catalogue selected at 60 mu (RIFSCz), to contrast the sky at these two wavelengths. Both surveys demonstrate the existence of extreme starbursts, with star-formation rates (SFRs) > 5000 Msun/yr. The maximum intrinsic star-formation rate appears to be ~30,000 Msun/yr. The sources with apparent SFR estimates higher than this are in all cases either lensed systems, blazars, or erroneous photometric redshifts. At redshifts of 3 to 5, the time-scale for the Herschel galaxies to make their current mass of stars at their present rate of formation ~ 10^8 yrs, so these galaxies are making a significant fraction of their stars in the current star-formation episode. Using dust mass as a proxy for gas mass, the Herschel galaxies at redshift 3 to 5 have gas masses comparable to their mass in stars. Of the 38 extreme starbursts in our Herschel survey for which we have more complete SED information, over 50% show evidence for QSO-like optical emission, or exhibit AGN dust tori in the mid-infrared SEDs. In all cases however the infrared luminosity is dominated by a starburst component. We derive a mean covering factor for AGN dust as a function of redshift and derive black hole masses and black hole accretion rates. There is a universal ratio of black-hole mass to stellar mass, ~ 10^{-3}, driven by the strong period of star-formation and black-hole growth at z = 1-5.

We present ALMA observations at 107.291 GHz (band 3) and 214.532 GHz (band 6) of GMASS 0953, a star-forming galaxy at z=2.226 hosting an obscured AGN that has been proposed as a progenitor of compact quiescent galaxies (QG). We measure for the first time the size of the dust and molecular gas emission of GMASS 0953 that we find to be extremely compact ($\sim$1 kpc). This result, coupled with a very high ISM density (n$\sim$10$^{5.5}$ cm$^{-3}$), a low gas mass fraction ($\sim$0.2) and a short gas depletion timescale ($\sim$150 Myr) imply that GMASS 0953 is experiencing an episode of intense star-formation in its central region that will rapidly exhaust its gas reservoirs, likely aided by AGN-induced feedback, confirming its fate as a compact QG. Kinematic analysis of the CO(6-5) line shows evidence of rapidly-rotating gas ($V_{rot}$=320$^{+92}_{-53}$ km s$^{-1}$), as observed also in a handful of similar sources at the same redshift. On-going quenching mechanisms could either destroy the rotation or leave it intact leading the galaxy to evolve into a rotating QG.

The local Seyfert 2 galaxy NGC5135, thanks to its almost face-on appearance, a bulge overdensity of stars, the presence of a large-scale bar, an AGN and a Supernova Remnant, is an excellent target to investigate the dynamics of inflows, outflows, star formation and AGN feedback. Here we present a reconstruction of the gas morphology and kinematics in the inner regions of this galaxy, based on the analysis of Atacama Large Millimeter Array (ALMA) archival data. To our purpose, we combine the available $\sim$100 pc resolution ALMA 1.3 and 0.45 mm observations of dust continuum emission, the spectroscopic maps of two transitions of the CO molecule (tracer of molecular mass in star forming and nuclear regions), and of the CS molecule (tracer of the dense star forming regions) with the outcome of the SED decomposition. By applying the $^{\rm 3D}$BAROLO software (3D-Based Analysis of Rotating Object via Line Observations), we have been able to fit the galaxy rotation curves reconstructing a 3D tilted-ring model of the disk. Most of the observed emitting features are described by our kinematic model. We also attempt an interpretation for the emission in few regions that the axisymmetric model fails to reproduce. The most relevant of these is a region at the northern edge of the inner bar, where multiple velocity components overlap, as a possible consequence of the expansion of a super-bubble.

We present a detailed analysis of the X-ray and molecular gas emission in the nearby galaxy NGC 34, to constrain the properties of molecular gas, and assess whether, and to what extent, the radiation produced by the accretion onto the central black hole affects the CO line emission. We analyse the CO Spectral Line Energy Distribution (SLED) as resulting mainly from Herschel and ALMA data, along with X-ray data from NuSTAR and XMM-Newton. The X-ray data analysis suggests the presence of a heavily obscured AGN with an intrinsic luminosity of L$_{\rm{1-100\,keV}} \simeq 4.0\times10^{42}$ erg s$^{-1}$. ALMA high resolution data ($\theta \simeq 0.2''$) allows us to scan the nuclear region down to a spatial scale of $\approx 100$ pc for the CO(6-5) transition. We model the observed SLED using Photo-Dissociation Region (PDR), X-ray-Dominated Region (XDR), and shock models, finding that a combination of a PDR and an XDR provides the best fit to the observations. The PDR component, characterized by gas density ${\rm log}(n/{\rm cm^{-3}})=2.5$ and temperature $T=30$ K, reproduces the low-J CO line luminosities. The XDR is instead characterised by a denser and warmer gas (${\rm log}(n/{\rm cm^{-3}})=4.5$, $T =65$ K), and is necessary to fit the high-J transitions. The addition of a third component to account for the presence of shocks has been also tested but does not improve the fit of the CO SLED. We conclude that the AGN contribution is significant in heating the molecular gas in NGC 34.

The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic surveys of galaxies that could be performed by SPICA. The surveys are designed to address the question how star formation and black hole accretion activities evolved over cosmic time through spectral diagnostics of the physical conditions of the interstellar/circumnuclear media in galaxies. On the basis of results obtained with Herschel far-IR photometric surveys of distant galaxies and Spitzer and AKARI near- to mid-IR spectroscopic observations of nearby galaxies, we estimate the numbers of the galaxies at redshift z > 0.5, which are expected to be detected in the PAH features or dust continuum by a wide (10 deg^2) or deep (1 deg^2) blind survey, both for a given observation time of 600 hours. As by-products of the wide blind survey, we also expect to detect debris disks, through the mid-IR excess above the photospheric emission of nearby main-sequence stars, and we estimate their number. We demonstrate that the SPICA mid-IR surveys will efficiently provide us with unprecedentedly large spectral samples, which can be studied further in the far-IR with SPICA.

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.

Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

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.

We present a coherent multi-band modelling of the CO Spectral Energy Distribution of the local Seyfert Galaxy NGC7130 to assess the impact of the AGN activity on the molecular gas. We take advantage of all the available data from X-ray to the sub-mm, including ALMA data. The high-resolution (~0.2") ALMA CO(6-5) data constrain the spatial extension of the CO emission down to ~70 pc scale. From the analysis of the archival CHANDRA and NuSTAR data, we infer the presence of a buried, Compton-thick AGN of moderate luminosity, L_2-10keV ~ 1.6x10^{43} ergs-1. We explore photodissociation and X-ray-dominated regions (PDRs and XDRs) models to reproduce the CO emission. We find that PDRs can reproduce the CO lines up to J~6, however, the higher rotational ladder requires the presence of a separate source of excitation. We consider X-ray heating by the AGN as a source of excitation, and find that it can reproduce the observed CO Spectral Energy Distribution. By adopting a composite PDR+XDR model, we derive molecular cloud properties. Our study clearly indicates the capabilities offered by current-generation of instruments to shed light on the properties of nearby galaxies adopting state-of-the art physical modelling.

The coeval AGN and galaxy evolution and the observed local relations between SMBHs and galaxy properties suggest some connection or feedback between SMBH growth and galaxy build-up. We looked for correlations between properties of X-ray detected AGN and their FIR detected host galaxies, to find quantitative evidences for this connection, highly debated in the latest years. We exploit the rich multi-wavelength data set available in the COSMOS field for a large sample (692 sources) of AGN and their hosts, in the redshift range $0.1<z<4$. We use X-ray data to select AGN and determine their properties (intrinsic luminosity and nuclear obscuration), and broad-band SED fitting to derive host galaxy properties (stellar mass $M_*$ and star formation rate SFR). We find that the AGN 2-10 keV luminosity ($L_{\rm X}$) and the host $8-1000~\mu m$ star formation luminosity ($L_{\rm IR}^{\rm SF}$) are significantly correlated. However, the average host $L_{\rm IR}^{\rm SF}$ has a flat distribution in bins of AGN $L_{\rm X}$, while the average AGN $L_{\rm X}$ increases in bins of host $L_{\rm IR}^{\rm SF}$, with logarithmic slope of $\sim0.7$, in the redshifts range $0.4<z<1.2$. We also discuss the comparison between the distribution of these two quantities and the predictions from hydro-dynamical simulations. Finally we find that the average column density ($N_H$) shows a positive correlation with the host $M_*$, at all redshifts, but not with the SFR (or $L_{\rm IR}^{\rm SF}$). This translates into a negative correlation with specific SFR. Our results are in agreement with the idea that BH accretion and SF rates are correlated, but occur with different variability time scales. The presence of a positive correlation between $N_H$ and host $M_*$ suggests that the X-ray $N_H$ is not entirely due to the circum-nuclear obscuring torus, but may also include a contribution from the host galaxy.

The survival of dust grains in galaxies depends on various processes. Dust can be produced in stars, it can grow in the interstellar medium and be destroyed by astration and interstellar shocks. In this paper, we assemble a few data samples of local and distant star-forming galaxies to analyse various dust-related quantities in low and high redshift galaxies, to study how the relations linking the dust mass to the stellar mass and star formation rate evolve with redshift. We interpret the available data by means of chemical evolution models for discs and proto-spheroid (PSPH) starburst galaxies. In particular, we focus on the dust-to-stellar mass (DTS) ratio, as this quantity represents a true measure of how much dust per unit stellar mass survives the various destruction processes in galaxies and is observable. The theoretical models outline the strong dependence of this quantity on the underlying star formation history. Spiral galaxies are characterised by a nearly constant DTS as a function of the stellar mass and cosmic time, whereas PSPHs present an early steep increase of the DTS, which stops at a maximal value and decreases in the latest stages. In their late starburst phase, these models show a decrease of the DTS with their mass, which allows us to explain the observed anti-correlation between the DTS and the stellar mass. The observed redshift evolution of the DTS ratio shows an increase from z~0 to z~1, followed by a roughly constant behaviour at 1<z<2.5. Our models indicate a steep decrease of the global DTS at early times, which implies an expected decrease of the DTS at larger redshift.

We use 3035 Herschel-SPIRE 500$\mu$m sources from 20.3 sq deg of sky in the HerMES Lockman, ES1 and XMM-LSS areas to estimate the star-formation rate density at z = 1-6. 500 mu sources are associated first with 350 and 250 mu sources, and then with Spitzer 24 mu sources from the SWIRE photometric redshift catalogue. The infrared and submillimetre data are fitted with a set of radiative-transfer templates corresponding to cirrus (quiescent) and starburst galaxies. Lensing candidates are removed via a set of colour-colour and colour-redshift constraints. Star-formation rates are found to extend from < 1 to 20,000 Mo/yr. Such high values were also seen in the all-sky IRAS Faint Source Survey. Star-formation rate functions are derived in a series of redshift bins from 0-6, combined with earlier far-infrared estimates, where available, and fitted with a Saunders et al (1990) functional form. The star-formation-rate density as a function of redshift is derived and compared with other estimates. There is reasonable agreement with both infrared and ultraviolet estimates for z < 3, but we find higher star-formation-rate densities than ultraviolet estimates at z = 3-6. Given the considerable uncertainties in the submillimetre estimates, we can not rule out the possibility that the ultraviolet estimates are correct. But the possibility that the ultraviolet estimates have seriously underestimated the contribution of dust-shrouded star-formation can also not be excluded.

We present new estimates of AGN accretion and star-formation luminosity in galaxies obtained for the local 12-$\mu$m sample of Seyfert galaxies (12MGS), by performing a detailed broad-band spectral energy distribution (SED) decomposition including the emission of stars, dust heated by star formation and a possible AGN dusty torus. Thanks to the availability of data from the X-rays to the sub-millimetre, we constrain and test the contribution of the stellar, AGN and star-formation components to the SEDs. The availability of Spitzer-IRS low resolution mid-infrared (mid-IR) spectra is crucial to constrain the dusty torus component at its peak wavelengths. The results of SED-fitting are also tested against the available information in other bands: the reconstructed AGN bolometric luminosity is compared to those derived from X-rays and from the high excitation IR lines tracing AGN activity like [Ne V] and [O IV]. The IR luminosity due to star-formation (SF) and the intrinsic AGN bolometric luminosity are shown to be strongly related to the IR line luminosity. Variations of these relations with different AGN fractions are investigated, showing that the relation dispersions are mainly due to different AGN relative contribution within the galaxy. Extrapolating these local relations between line and SF or AGN luminosities to higher redshifts, by means of recent Herschel galaxy evolution results, we then obtain mid- and far-IR line luminosity functions useful to estimate how many star-forming galaxies and AGN we expect to detect in the different lines at different redshifts and luminosities with future IR facilities (e.g., JWST, SPICA).

We combined the spectroscopic information from the 3D-HST survey with \textit{Herschel} data to characterize the H$\alpha$ dust attenuation properties of a sample of 79 main sequence star-forming galaxies at $z \sim 1$ in the GOODS-S field. The sample was selected in the far-IR, at $\lambda$=100 and/or 160 $\mu$m, and only includes galaxies with a secure H$\alpha$ detection (S/N$>$3). From the low resolution 3D-HST spectra we measured the redshifts and the H$\alpha$ fluxes for the whole sample (a factor of 1/1.2 was applied to the observed fluxes to remove the [NII] contamination). The stellar masses (M$_{\star}$), infrared (L$_{IR}$) and UV luminosities (L$_{UV}$) were derived from the SEDs by fitting multi-band data from GALEX near-UV to SPIRE 500 $\mu$m. We estimated the continuum extinction E$_{star}$(B-V) from both the IRX=L$_{IR}$/L$_{UV}$ ratio and the UV-slope, $\beta$, and found an excellent agreement between the two. The nebular extinction was estimated from comparison of the observed SFR$_{H\alpha}$ and SFR$_{UV}$. We obtained \emph{f}=E$_{star}$(B-V)/E$_{neb}$(B-V)=0.93$\pm$0.06, i.e. higher than the canonical value of \emph{f}=0.44 measured in the local Universe. Our derived dust correction produces good agreement between the H$\alpha$ and IR+UV SFRs for galaxies with SFR$\gtrsim$ 20 M$_{\odot}$/yr and M$_{\star} \gtrsim 5 \times 10^{10}$ M$_{\odot}$, while objects with lower SFR and M$_{\star}$ seem to require a smaller \emph{f}-factor (i.e. higher H$\alpha$ extinction correction). Our results then imply that the nebular extinction for our sample is comparable to that in the optical-UV continuum and suggest that the \emph{f}-factor is a function of both M$_{\star}$ and SFR, in agreement with previous studies.

We derive the CO luminosity function (LF) for different rotational transitions (i.e. (1-0), (3-2), (5-4)) starting from the Herschel LF by Gruppioni et al. and using appropriate $L'_{\rm CO} - L_{\rm IR}$ conversions for different galaxy classes. Our predicted LFs fit the data so far available at $z\approx0$ and $2$. We compare our results with those obtained by semi-analytical models (SAMs): while we find a good agreement over the whole range of luminosities at $z\approx0$, at $z\approx1$ and $z\approx2$ the tension between our LFs and SAMs in the faint and bright ends increases. We finally discuss the contribution of luminous AGN ($L_{X}>10^{44}\,\rm{erg\,s^{-1}}$) to the bright end of the CO LF concluding that they are too rare to reproduce the actual CO luminosity function at $z\approx2$.

We use a spectroscopic sample of 286 star-forming galaxies (SFGs) at 1<z<3 from the GMASS survey to study different star formation rate (SFR) estimators. Infrared (IR) data are used to derive empirical calibrations to correct ultraviolet (UV) and [OII]{\lambda}3727 luminosities for dust extinction and dust-corrected estimates of SFR. In the selection procedure we fully exploit the available spectroscopic information. On the basis of three continuum indices, we are able to identify and exclude from the sample galaxies in which old stellar populations might bring a non-negligible contribution to IR luminosity (LIR) and continuum reddening. Using Spitzer-MIPS and Herschel-PACS data we derive LIR for two-thirds of our sample. The LIR/LUV ratio is used as a probe of effective attenuation (AIRX) to search for correlations with continuum and spectroscopic features. The relation between AIRX and UV continuum slope ({\beta}) was tested for our sample and found to be broadly consistent with the literature results at the same redshift, though with a larger dispersion with respect to UV-selected samples. We find a correlation between the rest-frame equivalent width (EW) of the [OII]{\lambda}3727 line and {\beta}, which is the main result of this work. We therefore propose the [OII]{\lambda}3727 line EW as a dust attenuation probe and calibrate it through AIRX, though the assumption of a reddening curve is still needed to derive the actual attenuation towards the [OII]{\lambda}3727 line. We tested the issue of differential attenuation towards stellar continuum and nebular emission: our results are in line with the traditional prescription of extra attenuation towards nebular lines. A set of relations is provided that allows the recovery of the total unattenuated SFR from UV and [OII]{\lambda}3727 luminosities. (Abridged)

We present predictions for number counts and redshift distributions of galaxies detectable in continuum and in emission lines with the Mid-infrared (MIR) Instrument (SMI) proposed for the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We have considered 24 MIR fine-structure lines, four Polycyclic Aromatic Hydrocarbon (PAH) bands (at 6.2, 7.7, 8.6 and 11.3$\mu$m) and two silicate bands (in emission and in absorption) at 9.7$\mu$m and 18.0$\mu$m. Six of these lines are primarily associated with Active Galactic Nuclei (AGNs), the others with star formation. A survey with the SMI spectrometers of 1 hour integration per field-of-view (FoV) over an area of $1\,\hbox{deg}^2$ will yield $5\,\sigma$ detections of $\simeq 140$ AGN lines and of $\simeq 5.2\times10^{4}$ star-forming galaxies, $\simeq 1.6\times10^{4}$ of which will be detected in at least two lines. The combination of a shallow ($20.0\,\hbox{deg}^{2}$, $1.4\times10^{-1}$ h integration per FoV) and a deep survey ($6.9\times10^{-3}\,\hbox{deg}^{2}$, $635$ h integration time), with the SMI camera, for a total of $\sim$1000 h, will accurately determine the MIR number counts of galaxies and of AGNs over five orders of magnitude in flux density, reaching values more than one order of magnitude fainter than the deepest Spitzer $24\,\mu$m surveys. This will allow us to determine the cosmic star formation rate (SFR) function down to SFRs more than 100 times fainter than reached by the Herschel Observatory.

We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.

Highly obscured active galactic nuclei (AGN) are common in nearby galaxies, but are difficult to observe beyond the local Universe, where they are expected to significantly contribute to the black hole accretion rate density. Furthermore, Compton-thick (CT) absorbers (NH>10^24 cm^-2) suppress even the hard X-ray (2-10 keV) AGN nuclear emission, and therefore the column density distribution above 10^24 cm^-2 is largely unknown. We present the identification and multi-wavelength properties of a heavily obscured (NH>~10^25 cm^-2), intrinsically luminous (L(2-10keV)>10^44 erg s^-1) AGN at z=0.353 in the COSMOS field. Several independent indicators, such as the shape of the X-ray spectrum, the decomposition of the spectral energy distribution and X-ray/[NeV] and X-ray/6{\mu}m luminosity ratios, agree on the fact that the nuclear emission must be suppressed by a 10^25 cm^-2 column density. The host galaxy properties show that this highly obscured AGN is hosted in a massive star-forming galaxy, showing a barred morphology, which is known to correlate with the presence of CT absorbers. Finally, asymmetric and blueshifted components in several optical high-ionization emission lines indicate the presence of a galactic outflow, possibly driven by the intense AGN activity (L(Bol)/L(Edd) = 0.3-0.5). Such highly obscured, highly accreting AGN are intrinsically very rare at low redshift, whereas they are expected to be much more common at the peak of the star formation and BH accretion history, at z~2-3. We demonstrate that a fully multi-wavelength approach can recover a sizable sample of such peculiar sources in large and deep surveys such as COSMOS.

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.

We study the relation of AGN accretion, star formation rate (SFR), and stellar mass (M$_*$) using a sample of $\approx$ 8600 star-forming galaxies up to z=2.5 selected with \textit{Herschel} imaging in the GOODS and COSMOS fields. For each of them we derive SFR and M$_*$, both corrected, when necessary, for emission from an active galactic nucleus (AGN), through the decomposition of their spectral energy distributions (SEDs). About 10 per cent of the sample are detected individually in \textit{Chandra} observations of the fields. For the rest of the sample we stack the X-ray maps to get average X-ray properties. After subtracting the X-ray luminosity expected from star formation and correcting for nuclear obscuration, we derive the average AGN accretion rate for both detected sources and stacks, as a function of M$_{*}$, SFR and redshift. The average accretion rate correlates with SFR and with M$_*$. The dependence on SFR becomes progressively more significant at z$>$0.8. This may suggest that SFR is the original driver of these correlations. We find that average AGN accretion and star formation increase in a similar fashion with offset from the star-forming "main-sequence". Our interpretation is that accretion onto the central black hole and star formation broadly trace each other, irrespective of whether the galaxy is evolving steadily on the main-sequence or bursting.

We investigate the co-evolution of black-hole-accretion-rate (BHAR) and star-formation-rate (SFR) in $1.5<z<2.5$ galaxies displaying a greater diversity of star-forming properties compared to previous studies. We combine X-ray stacking and far-IR photometry of stellar mass-limited samples of normal star-forming, starburst and quiescent/quenched galaxies in the COSMOS field. We corroborate the existence of a strong correlation between BHAR (i.e. the X-ray luminosity, L_X), and stellar mass (M*) for normal star-forming galaxies, although find a steeper relation than previously reported. We find that starbursts show a factor of 3 enhancement in BHAR compared to normal SF galaxies (against a factor of 6 excess in SFR), while quiescents show a deficit of a factor 5.5 at a given mass. One possible interpretation of this is that the starburst phase does not coincide with cosmologically relevant BH growth, or that starburst-inducing mergers are more efficient at boosting SFR than BHAR. Contrary to studies based on smaller samples, we find the BHAR/SFR ratio of main sequence (MS) galaxies is not mass invariant, but scales weakly as M*^(0.43\pm0.09}, implying faster BH growth in more massive galaxies at $z\sim2$. Furthermore, BHAR/SFR during the starburst is a factor of 2 lower than in MS galaxies, at odds with the predictions of hydrodynamical simulations of merger galaxies that foresee a sudden enhancement of L_X/SFR during the merger. Finally, we estimate that the bulk of the accretion density of the Universe at $z\sim2$ is associated with normal star-forming systems, with only 6(+/-1)% and 11(+/-1)% associated with starburst and quiescent galaxies, respectively.

Heavily obscured, Compton Thick (CT, NH>10^24 cm^-2) AGN may represent an important phase in AGN/galaxy co-evolution and are expected to provide a significant contribution to the cosmic X-ray background (CXB). Through direct X-ray spectra analysis, we selected 39 heavily obscured AGN (NH>3x10^23 cm^-2) in the 2 deg^2 XMM-COSMOS survey. After selecting CT AGN based on the fit of a simple absorbed two power law model to the XMM data, the presence of CT AGN was confirmed in 80% of the sources using deeper Chandra data and more complex models. The final sample of CT AGN comprises 10 sources spanning a large range of redshift and luminosity. We collected the multi-wavelength information available for all these sources, in order to study the distribution of SMBH and host properties, such as BH mass (M_BH), Eddington ratio (\lambda_Edd), stellar mass (M*), specific star formation rate (sSFR) in comparison with a sample of unobscured AGN. We find that highly obscured sources tend to have significantly smaller M_BH and higher \lambda_edd with respect to unobscured ones, while a weaker evolution in M* is observed. The sSFR of highly obscured sources is consistent with the one observed in the main sequence of star forming galaxies, at all redshift. We also present optical spectra, spectral energy distribution (SED) and morphology for the sample of 10 CT AGN: all the available optical spectra are dominated by the stellar component of the host galaxy, and a highly obscured torus component is needed in the SED of the CT sources. Exploiting the high resolution Hubble-ACS images available, we conclude that these highly obscured sources have a significantly larger merger fraction with respect to other X-ray selected samples of AGN. Finally we discuss implications in the context of AGN/galaxy co-evolutionary models, and compare our results with the predictions of CXB synthesis models.

The most striking feature of the Cosmic Star Formation History (CSFH) of the Universe is a dramatic drop of the star formation (SF) activity, since z~1. In this work we investigate if the very same process of assembly and growth of structures is one of the major drivers of the observed decline. We study the contribution to the CSFH of galaxies in halos of different masses. This is done by studying the total SFR-halo mass-redshift plane from redshift 0 to redshift z~1.6 in a sample of 57 groups and clusters by using the deepest available mid- and far-infrared surveys conducted with Spitzer MIPS and Herschel PACS and SPIRE. Our results show that low mass groups provide a 60-80% contribution to the CSFH at z~1. Such contribution declines faster than the CSFH in the last 8 billion years to less than 10% at z<0.3, where the overall SF activity is sustained by lower mass halos. More massive systems provide only a marginal contribution (<10%) at any epoch. A simplified abundance matching method shows that the large contribution of low mass groups at z~1 is due to a large fraction (>50%) of very massive, highly star forming Main Sequence galaxies. Below z~1 a quenching process must take place in massive halos to cause the observed faster suppression of their SF activity. Such process must be a slow one though, as most of the models implementing a rapid quenching of the SF activity in accreting satellites significantly underpredicts the observed SF level in massive halos at any redshift. Starvation or the transition from cold to hot accretion would provide a quenching timescale of 1 Gyrs more consistent with the observations. Our results suggest a scenario in which, due to the structure formation process, more and more galaxies experience the group environment and, thus, the associated quenching process. This leads to the progressive suppression of their SF activity shaping the CSFH below z~1.

There is now a large consensus that the current epoch of the Cosmic Star Formation History (CSFH) is dominated by low mass galaxies while the most active phase at 1<z<2 is dominated by more massive galaxies, which undergo a faster evolution. Massive galaxies tend to inhabit very massive halos such as galaxy groups and clusters. We aim to understand whether the observed "galaxy downsizing" could be interpreted as a "halo downsizing", whereas the most massive halos, and their galaxy populations, evolve more rapidly than the halos of lower mass. Thus, we study the contribution to the CSFH of galaxies inhabiting group-sized halos. This is done through the study of the evolution of the Infra-Red (IR) luminosity function of group galaxies from redshift 0 to ~1.6. We use a sample of 39 X-ray selected groups in the Extended Chandra Deep Field South (ECDFS), the Chandra Deep Field North (CDFN), and the COSMOS field, where the deepest available mid- and far-IR surveys have been conducted with Spitzer MIPS and Hersche PACS. Groups at low redshift lack the brightest, rarest, and most star forming IR-emitting galaxies observed in the field. Their IR-emitting galaxies contribute <10% of the comoving volume density of the whole IR galaxy population in the local Universe. At redshift >~1, the most IR-luminous galaxies (LIRGs and ULIRGs) are preferentially located in groups, and this is consistent with a reversal of the star-formation rate vs .density anti-correlation observed in the nearby Universe. At these redshifts, group galaxies contribute 60-80% of the CSFH, i.e. much more than at lower redshifts. Below z~1, the comoving number and SFR densities of IR-emitting galaxies in groups decline significantly faster than those of all IR-emitting galaxies. Our results are consistent with a "halo downsizing" scenario and highlight the significant role of "environment" quenching in shaping the CSFH.

We present new estimates of redshift-dependent luminosity functions of IR lines detectable by SPICA/SAFARI and excited both by star formation and by AGN activity. The new estimates improve over previous work by using updated evolutionary models and dealing in a self consistent way with emission of galaxies as a whole, including both the starburst and the AGN component. New relationships between line and AGN bolometric luminosity have been derived and those between line and IR luminosities of the starburst component have been updated. These ingredients were used to work out predictions for the source counts in 11 mid/far-IR emission lines partially or entirely excited by AGN activity. We find that the statistics of the emission line detection of galaxies as a whole is mainly determined by the star formation rate, because of the rarity of bright AGNs. We also find that the slope of the line integral number counts is flatter than 2 implying that the number of detections at fixed observing time increases more by extending the survey area than by going deeper. We thus propose a wide spectroscopic survey of 1 hour integration per field-of-view over an area of 5 deg$^{2}$ to detect (at 5$\sigma$) $\sim$760 AGNs in [OIV]25.89$\mu$m $-$ the brightest AGN mid-infrared line $-$ out to $z\sim2$. Pointed observations of strongly lensed or hyper-luminous galaxies previously detected by large area surveys such as those by Herschel and by the SPT can provide key information on the galaxy-AGN co-evolution out to higher redshifts.

We compare various star formation rate (SFR) indicators for star-forming galaxies at $1.4<z<2.5$ in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR-stellar mass ($M_*$) relation, compared to that of the dominant main sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR-{$M_*$} relation is $\sim0.8-0.9$, in agreement with that derived from the dust-corrected UV-luminosity. Exploiting deeper 24$\mu$m-Spitzer data we have characterized a sub-sample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the $B$ band. The combination of Herschel with Spitzer data have allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the H$\alpha$ luminosity of a sample of star-forming galaxies at $1.4<z<1.7$. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.

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.

The goal of this work is to compute the number density of far-IR selected galaxies in the comoving frame and along the past lightcone of observationally constrained LTB models and to compare those results with their standard model counterparts. We derived integral number densities and differential number densities using different cosmological distance definitions in the LTB models. Then, we computed selection functions and consistency functions for the luminosity functions in the combined fields of the Herschel/PACS evolutionary probe (PEP) survey in both standard and void cosmologies, from which we derived the observed values of the above-mentioned densities. We used the Kolmogorov-Smirnov statistics to study both the evolution of the consistency functions and its connection to the evolution of the comoving density of sources. Finally, we fitted the power-law behaviour of the densities along the observer's past lightcone. The analysis of the comoving number density shows that the increased flexibility of the LTB models is not enough to fit the observed redshift evolution of the number counts, if it is specialised to a recent best-fit giant void parametrisation. The results for the power-law fits of the densities along the observer's past lightcone show general agreement across both cosmological models studied here around a slope of -2.5 $\pm$ 0.1 for the integral number density on the luminosity-distance volumes. The differential number densities show much bigger slope discrepancies. We conclude that the differential number densities on the observer's past lightcone were still rendered dependent on the cosmological model by the flux limits of the PEP survey. In addition, we show that an intrinsic evolution of the sources must be assumed to fit the comoving number-density redshift evolution in the giant void parametrisation for the LTB models used in this work.

We study a sample of Herschel-PACS selected galaxies within the GOODS-South and the COSMOS fields in the framework of the PACS Evolutionary Probe (PEP) project. Starting from the rich multi-wavelength photometric data-sets available in both fields, we perform a broad-band Spectral Energy Distribution (SED) decomposition to disentangle the possible active galactic nucleus (AGN) contribution from that related to the host galaxy. We find that 37 per cent of the Herschel-selected sample shows signatures of nuclear activity at the 99 per cent confidence level. The probability to reveal AGN activity increases for bright ($L_{\rm 1-1000} > 10^{11} \rm L_{\odot}$) star-forming galaxies at $z>0.3$, becoming about 80 per cent for the brightest ($L_{\rm 1-1000} > 10^{12} \rm L_{\odot}$) infrared (IR) galaxies at $z \geq 1$. Finally, we reconstruct the AGN bolometric luminosity function and the super-massive black hole growth rate across cosmic time up to $z \sim 3$ from a Far-Infrared (FIR) perspective. This work shows general agreement with most of the panchromatic estimates from the literature, with the global black hole growth peaking at $z \sim 2$ and reproducing the observed local black hole mass density with consistent values of the radiative efficiency $\epsilon_{\rm rad}$ ($\sim$0.07).

While continuum imaging data at far-infrared to sub-millimeter wavelengths have provided tight constraints on the population properties of dusty star forming galaxies up to high redshifts, future space missions like the Space Infra-Red Telescope for Cosmology and Astrophysics (SPICA) and ground based facilities like the Cerro Chajnantor Atacama Telescope (CCAT) will allow detailed investigations of their physical properties via their mid-/far-infrared line emission. We present updated predictions for the number counts and the redshift distributions of star forming galaxies spectroscopically detectable by these future missions. These predictions exploit a recent upgrade of evolutionary models, that include the effect of strong gravitational lensing, in the light of the most recent Herschel and South Pole Telescope data. Moreover the relations between line and continuum infrared luminosity are re-assessed, considering also differences among source populations, with the support of extensive simulations that take into account dust obscuration. The derived line luminosity functions are found to be highly sensitive to the spread of the line to continuum luminosity ratios. Estimates of the expected numbers of detections per spectral line by SPICA/SAFARI and by CCAT surveys for different integration times per field of view at fixed total observing time are presented. Comparing with the earlier estimates by Spinoglio et al. (2012) we find, in the case of SPICA/SAFARI, differences within a factor of two in most cases, but occasionally much larger. More substantial differences are found for CCAT.

SPICA is one of the key projects for the future. Not only its instrument suite will open up a discovery window but they will also allow to physically understand some of the phenomena that we still do not understand in the high-redshift universe. Using new homogeneous luminosity functions (LFs) in the Far-Ultraviolet (FUV) from VVDS and in the Far-Infrared (FIR) from Herschel/PEP and Herschel/HerMES, we studied the evolution of the dust attenuation with redshift. With this information, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density (SFRD_TOT). Our main conclusions are that: 1) the dust attenuation A_FUV is found to increase from z = 0 to $z \sim 1.2 and then starts to decrease until our last data point at z = 3.6; 2) the estimated SFRD confirms published results to z \sim 2. At z > 2, we observe either a plateau or a small increase up to z \sim 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT but the origin of this delay is not understood yet, and SPICA instruments will provide clues to move further in the physical understanding of this delay but also on the detection and redshift measurements of galaxies at higher redshifts. This work is further detailed in Burgarella et al. (2013).

The main energy-generating mechanisms in galaxies are black hole (BH) accretion and star formation (SF) and the interplay of these processes is driving the evolution of galaxies. MIR/FIR spectroscopy are able to distinguish between BH accretion and SF, as it was shown in the past by infrared spectroscopy from the space by the Infrared Space Observatory and Spitzer. Spitzer and Herschel spectroscopy together can trace the AGN and the SF components in galaxies, with extinction free lines, almost only in the local Universe, except for a few distant objects. One of the major goals of the study of galaxy evolution is to understand the history of the luminosity source of galaxies along cosmic time. This goal can be achieved with far-IR spectroscopic cosmological surveys. SPICA in combination with ground based large single dish submillimeter telescopes, such as CCAT, will offer a unique opportunity to do this. We use galaxy evolution models linked to the observed MIR-FIR counts (including Herschel) to predict the number of sources and their IR lines fluxes, as derived from observations of local galaxies. A shallow survey in an area of 0.5 square degrees, with a typical integration time of 1 hour per pointing, will be able to detect thousands of galaxies in at least three emission lines, using SAFARI, the far-IR spectrometer onboard of SPICA.

This is the first paper of a series aiming at investigating galaxy formation and evolution in the giant-void class of the Lemaitre-Tolman-Bondi (LTB) models that best fits current cosmological observations. Here we investigate the Luminosity Function (LF) methodology, and how its estimates would be affected by a change on the cosmological model assumed in its computation. Are the current observational constraints on the allowed Cosmology enough to yield robust LF results? We use the far-infrared source catalogues built on the observations performed with the Herschel/PACS instrument, and selected as part of the PACS evolutionary probe (PEP) survey. Schechter profiles are obtained in redshift bins up to z approximately 4, assuming comoving volumes in both the standard model, that is, Friedmann-Lemaitre-Robertson-Walker metric with a perfect fluid energy-momentum tensor, and non-homogeneous LTB dust models, parametrized to fit the current combination of results stemming from the observations of supernovae Ia, the cosmic microwave background, and baryonic acoustic oscillations. We find that the luminosity functions computed assuming both the standard model and LTB void models show in general good agreement. However, the faint-end slope in the void models shows a significant departure from the standard model up to redshift 0.4. We demonstrate that this result is not artificially caused by the used LF estimator which turns out to be robust under the differences in matter-energy density profiles of the models. The differences found in the LF slopes at the faint end are due to variation in the luminosities of the sources, which depend on the geometrical part of the model. It follows that either the standard model is over-estimating the number density of faint sources or the void models are under-estimating it.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

Using new homogeneous LFs in the FUV and in the FIR Herschel/PEP and Herschel/HerMES, we study the evolution of the dust attenuation with redshift. With this information in hand, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density SFRD_TOT. By integrating SFRD_TOT, we follow the mass building and analyze the redshift evolution of the stellar mass density (SMD). This letter aims at providing a complete view of star formation from the local universe to z = 4 and, using assumptions on earlier star formation history, compares this evolution to what was known before in an attempt to draw a homogeneous picture of the global evolution of star formation in galaxies. The main conclusions of this letter are: 1) the dust attenuation A_FUV is found to increase from z = 0 to z \sim 1.2 and then starts to decrease up to our last data point at z = 3.6; 2) the estimated SFRD confirms published results up to z = 2. At z > 2, we observe either a plateau or a small increase up to z = 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT and a likely origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) using assumptions (namely exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to z_form = 10, we integrate SFRD_TOT and find a good agreement with the published SMDs.

By making use of Herschel-PEP observations of the COSMOS and Extended Groth Strip fields, we have estimated the dependence of the clustering properties of FIR-selected sources on their 100um fluxes. Our analysis shows a tendency for the clustering strength to decrease with limiting fluxes: r0(S100um >8 mJy)~4.3 Mpc and r0(S100um >5 mJy)~5.8 Mpc. These values convert into minimum halo masses Mmin~10^{11.6} Msun for sources brighter than 8 mJy and Mmin~10^{12.4} Msun for S100um > 5 mJy galaxies. We show such an increase of the clustering strength to be due to an intervening population of z~2 sources, which are very strongly clustered and whose relative contribution, equal to about 10% of the total counts at S100um > 2 mJy, rapidly decreases for brighter flux cuts. By removing such a contribution, we find that z <~ 1 FIR galaxies have approximately the same clustering properties, irrespective of their flux level. The above results were then used to investigate the intrinsic dependence on cosmic epoch of the clustering strength of dusty star-forming galaxies between z~0 and z~2.5. In order to remove any bias in the selection process, the adopted sample only includes galaxies observed at the same rest-frame wavelength, lambda~60 um, which have comparable luminosities and therefore star-formation rates (SFR>~100 Msun/yr). Our analysis shows that the same amount of (intense) star forming activity takes place in extremely different environments at the different cosmological epochs. For z<~1 the hosts of such star forming systems are small, Mmin~10^{11} Msun, isolated galaxies. High (z~2) redshift star formation instead seems to uniquely take place in extremely massive/cluster-like halos, Mmin~10^{13.5} Msun, which are associated with the highest peaks of the density fluctuation field at those epochs. (abridged)

We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths and to determine the infrared luminosity function of galaxies down to LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun at z~2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3-sigma) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

[Abridged] We present a comprehensive investigation of the cosmological evolution of the luminosity function (LF) of galaxies and active galactic nuclei (AGN) in the infrared (IR). Based on the observed dichotomy in the ages of stellar populations of early-type galaxies on one side and late-type galaxies on the other, the model interprets the epoch-dependent LFs at z \geq 1.5 using a physical model for the evolution of proto-spheroidal galaxies and of the associated AGNs, while IR galaxies at z<1.5 are interpreted as being mostly late-type 'cold' (normal) and 'warm' (starburst) galaxies. As for proto-spheroids, in addition to the epoch-dependent LFs of stellar and AGN components separately, we have worked out the evolving LFs of these objects as a whole (stellar plus AGN component). The model provides a physical explanation for the observed positive evolution of both galaxies and AGNs up to z \simeq 2.5 and for the negative evolution at higher redshifts, for the sharp transition from Euclidean to extremely steep counts at (sub-)mm wavelengths, as well as the (sub-)mm counts of strongly lensed galaxies, that are hard to account for by alternative, physical or phenomenological, approaches. The evolution of late-type galaxies and of z<1.5 AGNs is described using a parametric phenomenological approach. The modeled AGN contributions to the counts and to the cosmic infrared background (CIB) are always subdominant with maximal at mid-IR wavelengths. The model provides a good fit to the multi-wavelength (from the mid-IR to millimeter waves) data on LFs at different redshifts and on number counts (both global and per redshift slices). A prediction of the present model is a systematic variation with wavelength of the populations dominating the counts and the contributions to the CIB intensity. The implied specific trend for cross-wavelength CIB power spectra is found to be in good agreement with the data.

(abridged) Far-infrared Herschel photometry from the PEP and HerMES programs is combined with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields. Based on this rich dataset, we reproduce the restframe UV to FIR ten-colors distribution of galaxies using a superposition of multi-variate Gaussian modes. The median SED of each mode is then fitted with a modified version of the MAGPHYS code that combines stellar light, emission from dust heated by stars and a possible warm dust contribution heated by an AGN. The defined Gaussian grouping is also used to identify rare sources. The zoology of outliers includes Herschel-detected ellipticals, very blue z~1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eight other popular methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent with each other within a median absolute deviation of 10-20%, the scatter being dominated more by fine tuning of the codes, rather than by the choice of SED templates. Finally, the library is used to classify 24 micron detected sources in PEP GOODS fields. AGN appear to be distributed in the stellar mass (M*) vs. star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with the tendency to lie on the high SFR side of the "main sequence". The incidence of warmer star-forming sources grows for objects with higher specific star formation rates (sSFR), and they tend to populate the "off-sequence" region of the M*-SFR-z space.

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 use Herschel-PACS far-infrared data, combined with previous multi-band information and mid-IR spectra, to properly account for the presence of an active nucleus and constrain its energetic contribution in luminous infrared (IR) sources at z\sim2. The sample is composed of 24 sources in the GOODS-South field, with typical IR luminosity of 10^{12} Lo. Data from the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN emission. We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy and an AGN component. A smooth-torus model for circum-nuclear dust is used to account for the direct and re-processed contribution from the AGN. We confirm that galaxies with typical L_{8-1000um}\sim10^{12}Lo at z\sim2 are powered predominantly by star-formation. An AGN component is present in nine objects (\sim35% of the sample) at the 3sigma confidence level, but its contribution to the 8-1000 um emission accounts for only \sim5% of the energy budget. The AGN contribution rises to \sim23% over the 5-30 um range (in agreement with Spitzer IRS results) and to \sim60% over the narrow 2-6 um range. The presence of an AGN is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating the presence of significant absorption, i.e. NH\sim10^{23}-10^{24} cm^{-2}. An additional source shows indications of obscured AGN emission from X-ray data. The comparison between the mid-IR--derived X-ray luminosities and those obtained from X-ray data suggests that obscuration is likely present also in the remaining six sources that harbour an AGN according to the SED-fitting analysis.

Star formation and accretion onto supermassive black holes in the nuclei of galaxies are the two most energetic processes in the Universe, producing the bulk of the observed emission throughout its history. We simulated the luminosity functions of star-forming and active galaxies for spectral lines that are thought to be good spectroscopic tracers of either phenomenon, as a function of redshift. We focused on the infrared (IR) and sub-millimeter domains, where the effects of dust obscuration are minimal. Using three different and independent theoretical models for galaxy formation and evolution, constrained by multi-wavelength luminosity functions, we computed the number of star-forming and active galaxies per IR luminosity and redshift bin. We converted the continuum luminosity counts into spectral line counts using relationships that we calibrated on mid- and far-IR spectroscopic surveys of galaxies in the local universe. Our results demonstrate that future facilities optimized for survey-mode observations, i.e., the Space Infrared Telescope for Cosmology and Astrophysics (SPICA) and the Cerro Chajnantor Atacama Telescope (CCAT), will be able to observe thousands of z>1 galaxies in key fine-structure lines, e.g., [SiII], [OI], [OIII], [CII], in a half-square-degree survey, with one hour integration time per field of view. Fainter lines such as [OIV], [NeV] and H_2 (0-0)S1 will be observed in several tens of bright galaxies at 1<z<2, while diagnostic diagrams of active-nucleus vs star-formation activity will be feasible even for normal z~1 galaxies. We discuss the new parameter space that these future telescopes will cover and that strongly motivate their construction.

Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (\Sigma(SFR/M)) in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift ~1.6. We use deep far-infrared PACS data at 100 and 160 um to accurately estimate the total star formation rate of the Luminous Infrared Galaxy population of 9 clusters with mass ~10^{15} M_{\odot}, and 9 groups/poor clusters with mass ~ 5 x 10^{13} M_{\odot}. Estimates of the field \Sigma(SFR/M) are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the universe. The field \Sigma(SFR/M) increases with redshift up to z~1 and it is constant thereafter. The evolution of the \Sigma(SFR/M)-z relation in galaxy systems is much faster than in the field. Up to redshift z~0.2, the field has a higher \Sigma(SFR/M) than galaxy groups and galaxy clusters. At higher redshifts, galaxy groups and the field have similar \Sigma(SFR/M), while massive clusters have significantly lower \Sigma(SFR/M) than both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z~1.6, where the group \Sigma(SFR/M) lies above the field \Sigma(SFR/M)-z relation. We discuss possible interpretations of our results in terms of the processes of downsizing, and star-formation quenching.

Two main modes of star formation are know to control the growth of galaxies: a relatively steady one in disk-like galaxies, defining a tight star formation rate (SFR)-stellar mass sequence, and a starburst mode in outliers to such a sequence which is generally interpreted as driven by merging. Such starburst galaxies are rare but have much higher SFRs, and it is of interest to establish the relative importance of these two modes. PACS/Herschel observations over the whole COSMOS and GOODS-South fields, in conjunction with previous optical/near-IR data, have allowed us to accurately quantify for the first time the relative contribution of the two modes to the global SFR density in the redshift interval 1.5<z<2.5, i.e., at the cosmic peak of the star formation activity. The logarithmic distributions of galaxy SFRs at fixed stellar mass are well described by Gaussians, with starburst galaxies representing only a relatively minor deviation that becomes apparent for SFRs more than 4 times higher than on the main sequence. Such starburst galaxies represent only 2% of mass-selected star forming galaxies and account for only 10% of the cosmic SFR density at z~2. Only when limited to SFR>1000M(sun)/yr, off-sequence sources significantly contribute to the SFR density (46+/-20%). We conclude that merger-driven starbursts play a relatively minor role for the formation of stars in galaxies, whereas they may represent a critical phase towards the quenching of star formation and morphological transformation in galaxies.

We have studied the evolution of high redshift quiescent galaxies over an effective area of ~1.7 deg^2 in the COSMOS field. Galaxies have been divided according to their star-formation activity and the evolution of the different populations has been investigated in detail. We have studied an IRAC (mag_3.6 < 22.0) selected sample of ~18000 galaxies at z > 1.4 with multi-wavelength coverage. We have derived accurate photometric redshifts (sigma=0.06) and other important physical parameters through a SED-fitting procedure. We have divided our sample into actively star-forming, intermediate and quiescent galaxies depending on their specific star formation rate. We have computed the galaxy stellar mass function of the total sample and the different populations at z=1.4-3.0. We have studied the properties of high redshift quiescent galaxies finding that they are old (1-4 Gyr), massive (log(M/M_sun)~10.65), weakly star forming stellar populations with low dust extinction (E(B-V) < 0.15) and small e-folding time scales (tau ~ 0.1-0.3 Gyr). We observe a significant evolution of the quiescent stellar mass function from 2.5 < z < 3.0 to 1.4 < z < 1.6, increasing by ~ 1 dex in this redshift interval. We find that z ~ 1.5 is an epoch of transition of the GSMF. The fraction of star-forming galaxies decreases from 60% to 20% from z ~ 2.5-3.0 to z ~ 1.4-1.6 for log(M/M_sun) > 11, while the quiescent population increases from 10% to 50% at the same redshift and mass intervals. We compare the fraction of quiescent galaxies derived with that predicted by theoretical models and find that the Kitzbichler & White (2007) model is the one that better reproduces the data. Finally, we calculate the stellar mass density of the star-forming and quiescent populations finding that there is already a significant number of quiescent galaxies at z > 2.5 (rho~6.0 MsunMpc^-3).

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high redshift galaxy populations, thus testing and superseeding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

ABRIDGED-This paper presents the first direct estimate of the 3D clustering properties of far-infrared sources up to z~3. This has been possible thanks to the Pacs Evolutionary Probe (PEP) survey of the GOODS South field performed with the PACS instrument onboard the Herschel Satellite. An analysis of the two-point correlation function over the whole redshift range spanned by the data reports for the correlation length, r_0~6.3 Mpc and r_0~6.7 Mpc, respectively at 100um and 160um, corresponding to dark matter halo masses M>~10^{12.4} M_sun. Objects at z~2 instead seem to be more strongly clustered, with r_0~19 Mpc and r_0~17 Mpc in the two considered PACS channels. This dramatic increase of the correlation length between z~1 and z~2 is connected with the presence of a wide, M>~10^{14} M_sun, filamentary structure which includes more than 50% of the sources detected at z~2. An investigation of the properties of such sources indicates the possibility for boosted star-forming activity in those which reside within the overdense environment with respect of more isolated galaxies found in the same redshift range. Lastly, we also present our results on the evolution of the relationship between luminous and dark matter in star-forming galaxies between z~1 and z~2. We find that the increase of (average) stellar mass in galaxies <M*> between z~1 and z~2 is about a factor 10 lower than that of the dark matter haloes hosting such objects (<M*>[z~1]/<M*>[z~2] ~ 0.4 vs M_{halo}[z~1]/M_{halo}[z~2] ~ 0.04). Our findings agree with the evolutionary picture of downsizing whereby massive galaxies at z~2 were more actively forming stars than their z~1 counterparts, while at the same time contained a lower fraction of their mass in the form of luminous matter.

We present a new backward evolution model for galaxies and AGNs in the infrared (IR). What is new in this model is the separate study of the evolutionary properties of the different IR populations (i.e. spiral galaxies, starburst galaxies, low-luminosity AGNs, "unobscured" type 1 AGNs and "obscured" type 2 AGNs) defined through a detailed analysis of the spectral energy distributions (SEDs) of large samples of IR selected sources. The evolutionary parameters have been constrained by means of all the available observables from surveys in the mid- and far-IR (source counts, redshift and luminosity distributions, luminosity functions). By decomposing the SEDs representative of the three AGN classes into three distinct components (a stellar component emitting most of its power in the optical/near-IR, an AGN component due to hot dust heated by the central black hole peaking in the mid-IR, and a starburst component dominating the far-IR spectrum) we have disentangled the AGN contribution to the monochromatic and total IR luminosity emitted by the different populations considered in our model from that due to star-formation activity. We have then obtained an estimate of the total IR luminosity density (and star-formation density - SFD - produced by IR galaxies) and the first ever estimate of the black hole mass accretion density (BHAR) from the IR. The derived evolution of the BHAR is in agreement with estimates from X-rays, though the BHAR values we derive from IR are slightly higher than the X-ray ones. Finally, we have simulated source counts, redshift distributions and SFD and BHAR that we expect to obtain with the future cosmological Surveys in the mid-/far-IR that will be performed with JWST-MIRI and SPICA-SAFARI.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

We present a search for Herschel-PACS counterparts of dust-obscured, high-redshift objects previously selected at submillimeter and millimeter wavelengths in the Great Observatories Origins Deep Survey North field. We detect 22 of 56 submillimeter galaxies (39%) with a SNR of >=3 at 100 micron down to 3.0 mJy, and/or at 160 micron down to 5.7 mJy. The fraction of SMGs seen at 160 micron is higher than that at 100 micron. About 50% of radio-identified SMGs are associated with PACS sources. We find a trend between the SCUBA/PACS flux ratio and redshift, suggesting that these flux ratios could be used as a coarse redshift indicator. PACS undetected submm/mm selected sources tend to lie at higher redshifts than the PACS detected ones. A total of 12 sources (21% of our SMG sample) remain unidentified and the fact that they are blank fields at Herschel-PACS and VLA 20 cm wavelength may imply higher redshifts for them than for the average SMG population (e.g., z>3-4). The Herschel-PACS imaging of these dust-obscured starbursts at high-redshifts suggests that their far-infrared spectral energy distributions have significantly different shapes than template libraries of local infrared galaxies.

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.

Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 microns observations toward the cluster lens Abell 2218, to penetrate the Herschel confusion limit. We derive source counts down to a flux density of 1 mJy at 100 microns and 2 mJy at 160 microns, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 beams/source in the two bands, approaching source density confusion at 160 microns. The slope of the counts below the turnover of the Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94-35 mJy at 100 microns and 1.47-35 mJy at 160 microns, we retrieve a cosmic infrared background (CIB) surface brightness of ~8.0 and ~9.9 nW m^-2 sr^-1, in the respective bands. These values correspond to 55% (+/- 24%) and 77% (+/- 31%) of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower fields, these figures increase to 62% (+/- 25%) and 88% (+/- 32%) CIB total fractions, resolved at 100 and 160 microns, disregarding the high uncertainties of DIRBE absolute values.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Sensitive Herschel far-infrared observations can break degeneracies that were inherent to previous studies of star formation in high-z AGN hosts. Combining PACS 100 and 160um observations of the GOODS-N field with 2Msec Chandra data, we detect ~20% of X-ray AGN individually at >3sig. The host far-infrared luminosity of AGN with L2-10~10^43erg/s increases with redshift by an order of magnitude from z=0 to z~1. In contrast, there is little dependence of far-infrared luminosity on AGN luminosity, for L2-10<~10^44erg/s AGN at z>~1. We do not find a dependence of far-infrared luminosity on X-ray obscuring column, for our sample which is dominated by L2-10<10^44erg/s AGN. In conjunction with properties of local and luminous high-z AGN, we interpret these results as reflecting the interplay between two paths of AGN/host coevolution. A correlation of AGN luminosity and host star formation is traced locally over a wide range of luminosities and also extends to luminous high z AGN. This correlation reflects an evolutionary connection, likely via merging. For lower AGN luminosities, star formation is similar to that in non-active massive galaxies and shows little dependence on AGN luminosity. The level of this secular, non-merger driven star formation increasingly dominates over the correlation at increasing redshift.

We exploit the deepest existing far-infrared (FIR) data obtained so far by Herschel at 100 and 160 um in the GOODS-N, as part of the PACS Evolutionary Probe (PEP) survey, to derive for the first time the evolution of the rest-frame 60-um, 90-um, and total IR luminosity functions (LFs) of galaxies and AGNs from z=0 to unprecedented high redshifts (z~2-3). The PEP LFs were computed using the 1/Vmax method. The FIR sources were classified by means of a detailed broad- band SED-fitting analysis and spectral characterisation. Based on the best-fit model results, k-correction and total IR (8-1000 um) luminosity were obtained for each source. LFs (monochromatic and total) were then derived for various IR populations separately in different redshift bins and compared to backward evolution model predictions. We detect strong evolution in the LF to at least z~2. Objects with SEDs similar to local spiral galaxies are the major contributors to the star formation density (SFD) at z< 0.3, then, as redshift increases, moderate SF galaxies - most likely containing a low-luminosity AGN - start dominating up to z ~= 1.5. At >1.5 the SFD is dominated by the contributions of starburst galaxies. In agreement with previous findings, the comoving IR LD derived from our data evolves approximately as (1 + z)^(3.8+/-0.3) up to z~1, there being some evidence of flattening up to z~2.

We use deep observations obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel space observatory to study the far-infrared (FIR) properties of submm and optically faint radio galaxies (SMGs and OFRGs). From literature we compiled a sample of 35 securely identified SMGs and nine OFRGs located in the GOODS-N and the A2218 fields. This sample is cross-matched with our PACS 100 um and 160 um multi-wavelength catalogs. About half of the galaxies in our sample are detected with PACS. The dust temperatures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 um (only the upper limits for the OFRGs) flux densities with a single modified (beta=1.5) black body function. The median dust temperature of our SMG sample is T=36+/-8K while for our OFRG sample it is T=47+/-3K. For both samples, median dust temperatures derived from Herschel data agree well with previous estimates. In particular, Chapman et al. (2005) found a dust temperature of T=36+/-7K for a large sample of SMGs assuming the validity of the FIR/radio correlation. The agreement between our studies confirms that the local FIR/radio correlation effectively holds at high redshift even though we find <q>=2.17+/-0.19, a slightly lower value than that observed in local systems. The median IR luminosities of SMGs and OFRGs are 4.6*10^12 Lsun and 2.6*10^12 Lsun, respectively. We note that for both samples the IR luminosity estimates from the radio part of the spectral energy distribution are accurate, while estimates from the mid-IR are considerably (x3) more uncertain. Our observations confirm the remarkably high luminosities of SMGs and thus imply median star-formation rates of 960Msun yr^-1 for SMGs with S(850um)>5mJy and 460Msun yr^-1 for SMGs with S(850um)>2mJy, assuming a Chabrier IMF and no dominant AGN contribution to the far-infrared luminosity.

We exploit deep observations of the GOODS-N field taken with PACS, on board of Herschel, as part of the PEP guaranteed time, to study the link between star formation and stellar mass in galaxies to z~2. Starting from a stellar mass-selected sample of ~4500 galaxies with mag[4.5mu]<23 (AB), we identify ~350 objects with a PACS detection at 100 or 160mu and ~1500 with only Spitzer 24 mu counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the SEDs. Consistently with other Herschel results, we find that LIR based only on 24 mu data is overestimated by a median factor ~1.8 at z~2, whereas it is underestimated (with our approach) up to a factor ~1.6 at 0.5<z<1.0. We then exploit this calibration to correct LIR based on the MIPS fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the SFR, of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of the relations between the star formation rate (SFR), the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, being almost flat at z<1.0 and reaching a slope of alpha=-0.50^(+0.13)_(-0.16) at z~2, at odds with recent works based on radio-stacking analysis at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ~15 for massive M>10^(11) (M_sun) galaxies from z=0 to z=2, and seems to flatten at z>1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, implying that they have formed their stars earlier and more rapidly than their low mass counterparts.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 um) by exploiting Herschel/PACS observations of the GOODS-N, Lockman Hole, and COSMOS fields in the PACS Evolutionary Probe (PEP) guaranteed-time survey. The GOODS-N data reach 3 sigma depths of ~3.0 mJy at 100 um and ~5.7 mJy at 160 um. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 um. Differential number counts extend from a few mJy up to 100-200 mJy, and are approximated as a double power law, with the break lying between 5 and 10 mJy. The available ancillary information allows us to split number counts into redshift bins. At z<=0.5 we isolate a class of luminous sources (L(IR)~1e11 Lsun), whose SEDs resemble late-spiral galaxies, peaking at ~130 um restframe and significantly colder than what is expected on the basis of pre-Herschel models. By integrating number counts over the whole covered flux range, we obtain a surface brightness of 6.36 +/- 1.67 and 6.58 +/-1.62 [nW m^-2 sr^-1] at 100 and 160 um, resolving ~45% and ~52% of the CIB, respectively. When stacking 24 um sources, the inferred CIB lies within 1.1 sigma and 0.5 sigma from direct measurements in the two bands, and fractions increase to 50% and 75%.Most of this resolved CIB fraction was radiated at z<=1.0, with 160 um sources found at higher redshift than 100 um ones.

We combine near-to-mid-IR Spitzer data with shorter wavelength observations (optical to X-rays) to get insights on the properties of a sample of luminous, obscured Active Galactic Nuclei (AGN). We aim at modeling their broad-band Spectral Energy Distributions (SEDs) in order to estimate the main parameters related to the dusty torus. The sample comprises 16 obscured high-redshift (0.9<z<2.1) xray luminous quasars (L_2-10 ~ 10^44 erg s-1) selected from the HELLAS2XMM survey. The SEDs are described by a multi-component model including a stellar component, an AGN component and a starburst. The majority (~80%) of the sources show moderate optical depth (tau_9.7um<3) and the derived column densities N_H are consistent with the xray inferred values (10^22 <N_H< 3x10^23 cm-2) for most of the objects, confirming that the sources are moderately obscured Compton-thin AGN. Accretion luminosities in the range 5x10^44 < Lbol < 4x10^46 erg s-1 are inferred. We compare model luminosities with those obtained by integrating the observed SED, finding that the latter are lower by a factor of ~2 in the median. The discrepancy can be as high as an order of magnitude for models with high optical depth (tau_9.7um=10). The ratio between the luminosities obtained by the fitting procedure and from the observed SED suggest that, at least for Type~2 AGN, observed bolometric luminosities are likely to underestimate intrinsic ones and the effect is more severe for highly obscured sources. Bolometric corrections from the hard X-ray band are computed and have a median value of k_2-10kev ~ 20. The obscured AGN in our sample are characterized by relatively low Eddington ratios (median lambda_Edd~0.08). On average, they are consistent with the Eddington ratio increasing at increasing bolometric correction (e.g. Vasudevan & Fabiam 2009).

We show the feasibility of spectroscopic cosmological surveys with the SAFARI instrument onboard of SPICA. The work is done through simulations that make use of both empirical methods, i.e. the use of observed luminosity functions and theoretical models for galaxy formation and evolution. The relations assumed between the line emission to trace AGN and star formation activity have been derived from the observations of local samples of galaxies. The results converge to indicate the use of blind spectroscopy with the SAFARI FTS at various resolutions to study galaxy evolution from the local to the distant (z~3) Universe. Specifically, two different and independent galaxy evolution models predict about 7-10 sources to be spectroscopically detected in more than one line in a 2'x 2'SAFARI field of view, down to the expected flux limits of SAFARI, with about 20% of sources to be detected at z>2. SPICA-SAFARI will be therefore excellent at detecting high-z sources and at assessing in a direct way their nature (e.g whether mainly AGN or Star Formation powered) thanks to blind spectroscopy.

We present a catalog of optical spectroscopic identifications of sources detected by Spitzer at 3.6 or 24 micron down to 10 and 280 microJy, respectively, in the SWIRE/XMM-Newton/ELAIS-S1 field and classified via line width analysis and diagnostic diagrams. A total of 1376 sources down to R~24.2 mag have been identified (1362 detected at 3.6 micron, 419 at 24 micron, and 405 at both) by low-resolution optical spectroscopy carried out with FORS2, VIMOS, and EFOSC2 at the Very Large Telescope and 3.6m ESO telescopes. The spectroscopic campaigns have been carried out over the central 0.6 square degrees area of ELAIS-S1 which, in particular, has also been observed by XMM-Newton and Chandra. We find the first direct optical spectroscopic evidence that the fraction of active galactic nuclei (AGN; mostly AGN2) increases with increasing F(24 micron)/F(R) ratio, reaching values of 70(+/-20)% in the range 316<F(24 micron)F(R)<1000. We present an IRAC-MIPS color-color diagram able to separate AGN1 from obscured AGN2 candidates. After having corrected for the spectroscopic incompleteness of our sample, it results that the AGN fraction at F(24 micron)=0.8 mJy is ~22(+/-7)% and decreases slowly to ~19(+/-5)% down to F(24 micron)=0.3 mJy.

Recent works have suggested that selection criteria based on MIR colors can be used to reveal a population of dust-enshrouded, extremely luminous quasars at z>1. However the X-ray spectral properties of these intriguing sources still remain largely unexplored. We report on an X-ray spectroscopic study of a sample of 44 very bright mid-IR galaxies with extreme mid-IR to optical flux ratios (MIR/O>2000). The X-ray coverage of the sample is highly inhomogeneous (from snap-shot 5 ks Chandra observations to medium-deep XMM exposures of 70 ks) and, consequently, a sizable fraction of them (~43%) remains undetected in the 0.5-10 keV band. The vast majority (95%) of the detected sources (23) show an absorption column density NH>10e22 cm-2 and, remarkably, we also find that 50% of them can be classified as Type 2 quasars on the basis of their absorption properties and X-ray luminosity. Moreover, most of the X-ray undetected sources show extreme mid-IR colors, consistent with being luminous AGN-powered objects, suggesting they might host heavily obscured (possibly Compton-thick) quasars in X-rays. This demonstrates that our selection criteria applied to a wide area survey is very efficient in finding a large number of Type 2 quasars at z > 1. The existence of this class of very powerful, obscured quasars at high z could have important implications in the context of the formation and cosmological evolution of accreting supermassive black holes and their host galaxies.

We present multi-wavelength observations (from optical to sub-millimeter, including Spitzer and SCUBA) of H2XMMJ 003357.2-120038 (also GD158_19), an X-ray selected, luminous narrow-line (Type 2) quasar at z=1.957 selected from the HELLAS2XMM survey. Its broad-band properties can be reasonably well modeled assuming three components: a stellar component to account for the optical and near-IR emission, an AGN component (i.e., dust heated by an accreting active nucleus), dominant in the mid-IR, with an optical depth at 9.7 micron along the line of sight (close to the equatorial plane of the obscuring matter) of tau(9.7)=1 and a full covering angle of the reprocessing matter (torus) of 140 degrees, and a far-IR starburst component (i.e., dust heated by star formation) to reproduce the wide bump observed longward of 70 micron. The derived star-formation rate is about 1500 solar masses per year. The overall modeling indicates that GD158_19 is a high-redshift X-ray luminous, obscured quasar with coeval powerful AGN activity and intense star formation. It is probably caught before the process of expelling the obscuring gas has started, thus quenching the star formation.

We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey (ELAIS) with S15 > 0.8 mJy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. We classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of source: those dominated by an unobscured AGN (QSOs), obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of the 7.7-micron PAH band into the selection 15 micron band, while AGN spread over the 0.5 < z < 3.1 range. Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Mo yr^-1, but the implied FIR luminosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of the PAH carriers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, but could instead indicate a connection between obscuration and star formation. However, the observed correlation between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of the AGN emission in these sources being produced in a dust torus.

We present optical identifications and a multi-band catalogue of a sample of 478 X-ray sources in the XMM and Chandra surveys of the central 0.6 deg^2 of the ELAIS-S1 field. The optical/infrared counterpart of each X-ray source was identified using R and IRAC 3.6 um bands. This method was complemented by the precise positions obtained through Chandra observations. Approximately 94% of the counterparts are detected in the R band, while the remaining are blank fields in the optical down to R~24.5, but have a near-infrared counterpart detected by IRAC within 6 arcsec from the XMM centroid. The multi-band catalogue contains photometry in ten photometric bands (B to the MIPS 24 um). We determined redshift and classification for 237 sources (~50% of the sample) brighter than R=24. We classified 47% of the sources with spectroscopic redshift as broad-line active galactic nuclei (BL AGNs) with z=0.1-3.5, while sources without broad-lines are about 46% of the spectroscopic sample and are found up to z=2.6. We identified 11 type 2 QSOs among the sources with X/O>8, with z=0.9-2.6, high 2-10 keV luminosity (log(L2-10keV)>=43.8 erg/s) and hard X-ray colors suggesting large absorbing columns at the rest frame (logN_H up to 23.6 cm^-2). BL AGNs show on average blue optical-to-near-infrared colors, softer X-ray colors and X-ray-to-optical colors typical of optically selected AGNs. Conversely, narrow-line sources show redder optical colors, harder X-ray flux ratio and span a wider range of X-ray-to-optical colors. On average the SEDs of high-luminosity BL AGNs resemble the power-law typical of unobscured AGNs. The SEDs of NOT BL AGNs are dominated by the galaxy emission in the optical/near-infrared, and show a rise in the mid-infrared which suggests the presence of an obscured active nucleus.

We present the broad-band Spectral Energy Distributions (SEDs) of the largest available highly (72%) complete spectroscopic sample of mid-infrared (MIR) selected galaxies and AGN at intermediate redshift. The sample contains 203 extragalactic sources from the 15-micron survey in the ELAIS-SWIRE field S1, all with measured spectroscopic redshift. Most of these sources have full multi-wavelength coverage from the far-UV to the far-infrared and lie in the redshift range 0.1<z<1.3. Due to its size, this sample allows us for the first time to characterise the spectral properties of the sources responsible for the strong evolution observed in the MIR. Based on SED-fitting technique we have classified the MIR sources, identifying AGN signatures in about 50% of them. This fraction is significantly higher than that derived from optical spectroscopy (~29%) and is due in particular to the identification of AGN activity in objects spectroscopically classified as galaxies. It is likely that in most of our objects, the AGN is either obscured or of low-luminosity, and thus it does not dominate the energetic output at any wavelength, except in the MIR, showing up just in the range where the host galaxy SED has a minimum. The fraction of AGN strongly depends on the flux density, with that derived through the SED-fitting being about 20% at S(15)~0.5-1 mJy and gradually increasing up to 100% at S(15)>10 mJy, while that obtained from optical spectroscopy never being >30%, even at the higher flux densities. The results of this work will be very useful for updating all the models aimed at interpreting the deep infrared survey data and, in particular, for constraining the nature and the role of dust-obscured systems in the intermediate/high-redshift Universe.

A large population of heavily obscured, Compton Thick AGNs is predicted by models of galaxy formation, models of Cosmic X-ray Background and by the ``relic'' super-massive black-hole mass function measured from local bulges. However, so far only a handful of Compton thick AGNs have been possibly detected using even the deepest Chandra and XMM surveys. Compton-thick AGNs can be recovered thanks to the reprocessing of the AGN UV emission in the infrared by selecting sources with AGN luminosity's in the mid-infrared and faint near-infrared and optical emission. To this purpose, we make use of deep HST, VLT, Spitzer and Chandra data on the Chandra Deep Field South to constrain the number of Compton thick AGN in this field. We show that sources with high 24$\mu$m to optical flux ratios and red colors form a distinct source population, and that their infrared luminosity is dominated by AGN emission. Analysis of the X-ray properties of these extreme sources shows that most of them (80$\pm15%$) are indeed likely to be highly obscured, Compton thick AGNs. The number of infrared selected, Compton thick AGNs with 5.8$\mu$m luminosity higher than $10^{44.2}$ erg s$^{-1}$ turns out to be similar to that of X-ray selected, unobscured and moderately obscured AGNs with 2-10 keV luminosity higher than $10^{43}$ erg s$^{-1}$ in the redshift bin 1.2-2.6. This ``factor of 2'' source population is exactly what it is needed to solve the discrepancies between model predictions and X-ray AGN selection.

Context: The counts of galaxies and AGN in the mid infra-red (MIR) bands are important instruments for studying their cosmological evolution. However, the classic spectral line ratios techniques can become misleading when trying to properly separate AGN from starbursts or even from apparently normal galaxies. Aims: We use X-ray band observations to discriminate AGN activity in previously classified MIR-selected starburst galaxies and to derive updated AGN1 and (Compton thin) AGN2 counts at 15 um. Methods: XMM observations of the ELAIS-S1 15um sample down to flux limits ~2x10^-15 erg cm^-2 s^-1 (2-10 keV band) were used. We classified as AGN all those MIR sources with a unabsorbed 2-10 keV X-ray luminosity higher that ~10^42 erg/s. Results: We find that at least about 13(+/-6) per cent of the previously classified starburst galaxies harbor an AGN. According to these figures, we provide an updated estimate of the counts of AGN1 and (Compton thin) AGN2 at 15 um. It turns out that at least 24% of the extragalactic sources brighter than 0.6 my at 15 um are AGN (~13% contribution to the extragalactic background produced at fluxes brighter than 0.6 mJy).

Aims: We aim at estimating the spectral energy distributions (SEDs) and the physical parameters related to the black holes harbored in eight high X-ray-to-optical (F_X/F_R>10) obscured quasars at z>0.9 selected in the 2--10 keV band from the HELLAS2XMM survey. Methods: We use IRAC and MIPS 24 micron observations, along with optical and Ks-band photometry, to obtain the SEDs of the sources. The observed SEDs are modeled using a combination of an elliptical template and torus emission (using the phenomenological templates of Silva et al. 2004) for six sources associated with passive galaxies; for two point-like sources, the empirical SEDs of red quasars are adopted. The bolometric luminosities and the M_BH-L_K relation are used to provide an estimate of the masses and Eddington ratios of the black holes residing in these AGN. Results: All of our sources are detected in the IRAC and MIPS (at 24 micron) bands. The SED modeling described above is in good agreement with the observed near- and mid-infrared data. The derived bolometric luminosities are in the range ~10^45-10^47 erg s^-1, and the median 2--10 keV bolometric correction is ~25, consistent with the widely adopted value derived by Elvis et al. (1994). For the objects with elliptical-like profiles in the K_s band, we derive high stellar masses (0.8-6.2)X10^11 Mo, black hole masses in the range (0.2-2.5)X10^9 Mo, and Eddington ratios L/L_Edd<0.1, suggesting a low-accretion phase.

Recent X-ray surveys have provided a large number of high-luminosity, obscured Active Galactic Nuclei (AGN), the so-called Type 2 quasars. Despite the large amount of multi-wavelength supporting data, the main parameters related to the black holes harbored in such AGN are still poorly known. Here we present the results obtained for a sample of eight Type 2 quasars in the redshift range 0.9-2.1 selected from the HELLAS2XMM survey, for which we used Ks-band, Spitzer IRAC and MIPS data at 24 micron to estimate bolometric corrections, black hole masses, and Eddington ratios.

We exploit very deep mid-IR (MIR) and X-ray observations by Spitzer and Chandra in the GOODS North to identify signs of hidden (either starburst or AGN) activity in spheroidal galaxies between z~0.3 and 1. Our reference is a complete sample of 168 morphologically classified spheroidal galaxies with z[AB]<22.5 selected from GOODS ACS imaging. Nineteen of these have 24 micron detections in the GOODS catalogue, half of which have an X-ray counterpart in the 2 Ms Chandra catalogue, while about 25% have 1.4 GHz fluxes larger than 40 microJy. Traces of hidden activity in the spheroidal population are also searched for in the deep X-ray images and 14 additional galaxies are detected in X-rays only. The nature of the observed MIR emissions is investigated by modelling their SEDs based on the available multi-wavelength photometry, including X-ray, UV, optical, near-IR, MIR and radio fluxes, and optical spectroscopy. The amount of dust derived from the IR emission observed by Spitzer appears in excess of that expected by mass loss from evolved stars. Our analysis of the X-ray and MIR properties leads us to conclude that at least 8 of the 19 24-um bright sources should hide an obscured AGN, while the X-ray undetected sources are more likely dominated by star formation. [abridged] We conclude that ~30 objects (~20%) of the original flux-limited sample of 168 spheroidal galaxies in the GOODS-North are detected during phases of prominent activity, of both stellar and quasar origin. Due to the short expected lifetimes of the IR and X-ray emissions, this fraction might imply a significant level of activity in this class of galaxies during the relatively recent cosmic epochs, z~0.3 to z~1, under investigation.

We report results from a Spitzer GO-1 program of IRS spectroscopy of a large sample of Luminous Infrared Galaxies and quasars selected from the European Large Area ISO Survey (ELAIS) which have a wide multiwavelength coverage, including ISOCAM, ISOPHOT, IRAC and MIPS (from SWIRE), and optical photometry. We present the sample selection and results from the IRS spectroscopy.

A proper analysis of the evolution of sources emitting in the Mid-Infrared is strongly dependent on their broad-band spectral properties (SEDs) at different redshifts and luminosities and on a reliable classification allowing to disentangle AGN from star-formation activity. The diagnostic diagrams based on the optical line ratios are often ambiguous and/or misleading not allowing a proper separation of the galaxy/AGN populations. Thanks to the combination of deep Spitzer and X-rays data a much better census of the hidden AGN activity and dust-obscured star-forming galaxies can be obtained, constraining galaxy and AGN evolutionary models.

We have surveyed with XMM-Newton the central ~0.6 deg2 region of the ELAIS-S1 field down to flux limits of ~5.5X10-16 cgs (0.5-2 keV, S band), ~2X10-15 cgs (2-10 keV, H band), and ~4X10-15 cgs (5-10 keV, HH band). We detect a total of 478 sources, 395 and 205 of which detected in the S and H bands respectively. We identified 7 clearly extended sources and estimated their redshift through X-ray spectral fits with thermal models. In four cases the redshift is consistent with z=0.4. We have computed the angular correlation function of the sources in the S and H bands, finding best fit correlation angles theta_0=5.2+/-3.8 arcsec and theta_0=12.8+/-7.8 arcsec respectively. A rough estimate of the present-day correlation length r_0 can be obtained inverting the Limber equation and assuming an appropriate redshift distribution dN/dz. The results range between 12.8 and 9.8 h-1 Mpc in the S band and between 17.9 and 13.4 h-1 Mpc in the H band, with 30-40% statistical errors, assuming either smooth redshift distributions or redshift distributions with spikes accounting for the presence of a structure at z=0.4. The relative density of the S band sources is higher near the clusters and groups at z~0.4 and extends toward East and toward South/West. This suggests that the structure is complex, with a size comparable to the full XMM-Newton field. Conversely, the highest relative source densities of the H band sources are located in the central-west region of the field.

The ESO-Spitzer extragalactic Imaging Survey (ESIS) is the optical follow up of the Spitzer Wide-Area InfraRed Extragalactic (SWIRE) survey in the ELAIS-S1 area. This paper presents B, V, R Wide Field Imager observations of the first 1.5 square degree of the ESIS survey. Data reduction is described including astrometric calibration, illumination and color corrections, completeness and photometric accuracy estimates. Number counts and color distributions are compared to literature observational and theoretical data, including non-evolutionary, PLE, evolutionary and semi-analytic Lambda-CDM galaxy models, as well as Milky Way stellar predictions. ESIS data are in good agreement with previous works and are best reproduced by evolutionary and hierarchical Lambda-CDM scenarios. The ELAIS-S1 area benefits from extensive follow-up from X-ray to radio frequencies: some potential uses of the multi-wavelength observations are illustrated. Optical-Spitzer color-color plots promise to be very powerful tools to disentangle different classes of sources (e.g. AGNs, starbursts, quiescent galaxies). Ultraviolet GALEX data are matched to optical and Spitzer samples, leading to a discussion of galaxy properties in the UV-to-24 microns color space. The spectral energy distribution of a few objects, from the X-rays to the far-IR are presented as examples of the multi-wavelength study of galaxy emission components in different spectral domains.

We study the evolution of the luminosity function (LF) of type-1 and type-2 AGN in the mid-infrared, and derive their contribution to the Cosmic InfraRed Background (CIRB) and the expected deep source counts to be observed by Spitzer at 24 micron. The sample of type-1 and type-2 AGN was selected at 15 micron (ISO) and 12 micron (IRAS), and classified on the basis of their optical spectra. Local templates of type-1 and type-2 AGN have been used to derive the intrinsic 15 micron luminosities. We adopted an evolving smooth two-power law shape of the LF, whose parameters have been derived using an un-binned maximum likelihood method. We find that the LF of type-1 AGN is compatible with a pure luminosity evolution (L(z)=L(0)(1+z)^k_L) model where k_L~2.9. A small flattening of the faint slope of the LF with increasing redshift is favoured by the data. A similar evolutionary scenario is found for the type-2 population with a rate k_L ranging from ~1.8 to 2.6, depending significantly on the adopted mid-infrared spectral energy distribution. Also for type-2 AGN a flattening of the LF with increasing redshift is suggested by the data, possibly caused by the loss of a fraction of type-2 AGN hidden within the optically classified starburst and normal galaxies. The type-1 AGN contribution to the CIRB at 15 micron is (4.2-12.1) x 10e-11 W m^-2 sr^-1, while the type-2 AGN contribution is (5.5-11.0) x 10e-11 W m^-2 sr^-1. We expect that Spitzer will observe, down to a flux limit of S_24 = 0.01 mJy, a density of ~1200 deg^-2 type-1 and ~1000 deg^-2 type-2 optically classified AGN. The derived total contribution of the AGN galaxies to the CIRB (4-10%) and Spitzer counts should be considered as lower limits, because of a possible loss of type-2 sources caused by the optical classification.

We address the question of how to deal with confusion limited surveys in the mid-infrared domain by using informations from higher frequency observations over the same sky regions. Such informations, once applied to apparently extended mid-infrared sources, which are indeed ``blends'' of two or more different sources, allow us to disentangle the single counterparts and to split the measured flux density into different components. We present the application of this method to the 24 micron Spitzer archival data in the GOODS-EN1 test field, where apparently extended, ``blended'' sources constitute about 20% of a reliable sample of 983 sources detected above the 5-sigma threshold down to 23 microJy. As higher frequency data-set we have considered the public IRAC images and catalogues on the same field. We show that the 24 micron sample is almost unbiased down to ~40 microJy and the careful application of the deblending procedure does not require any statistical completeness correction (at least at the flux level considered). This is probed by direct comparison of our results with those of Chary et al. (2004), who analysed the same data-set through extensive Monte Carlo simulations. The deblending procedure reduces of about 30% the confusion limit of the MIPS 24 micron survey, allowing one to obtain reliable source counts down to ~40 microJy. The extrapolation of the source counts down to fainter fluxes suggests that our 24 micron sample is able to resolve ~62% of the cosmic background down to a flux level of 38 microJy.

We investigate the clustering properties of galaxies in the recently completed ELAIS-S1 redshift survey through their spatial two point autocorrelation function. We used a sub-sample of the ELAIS-S1 catalog covering approximately 4 deg^2 and consisting of 148 objects selected at 15 micron with a flux >0.5 mJy and redshift z<0.5. We detected a positive signal in the correlation function that, in the range of separations 1-10 h mpc is well approximated by a power law with a slope gamma = 1.4 and a correlation length s =5.4 h mpc. This result is in good agreement with the redshift-space correlation function measured in more local samples of mid infrared selected galaxies like the IRAS PSC$z$ redshift survey. This suggests a lack of significant clustering evolution of infrared selected objects out to $z=0.5$ that is further confirmed by the consistency found between the correlation functions measured in a local (z<0.2) and a distant (0.2<z<0.5) subsample of ELAIS-S1 galaxies. We also confirm that optically selected galaxies in the local redshift surveys, especially those of the SDSS sample, are significantly more clustered than infrared objects.

We present the Final Analysis of the European Large Area ISO Survey (ELAIS) 15 $\mu$m observations, carried out with the ISOCAM instrument on board the Infrared Space Observatory (ISO). The data reduction method, known as LARI method, is based on a mathematical model of the detector's behaviour and was specifically designed for the detection of faint sources in ISO-CAM/PHOT data. The method is fully interactive and leads to very reliable and complete source lists. The resulting catalogue includes 1923 sources detected with $S/N > 5$ in the {0.5 -- 100 mJy} flux range and over an area of 10.85 {deg$^2$} split into four fields, making it the largest non-serendipitous extragalactic source catalogue obtained to date from ISO data. This paper presents the concepts underlying the data reduction method together with its latest enhancements. The data reduction process, the production and basic properties of the resulting catalogue are discussed. The catalogue quality is assessed by means of detailed simulations, optical identifications and comparison with previous analyses.

The comparison between the new Spitzer data at 24 micron and the previous ISOCAM data at 15 micron is a key tool to understand galaxy properties and evolution in the infrared and to interpret the observed number counts, since the combination of Spitzer with the ISO cosmological surveys provides for the first time the direct view of the Universe in the Infrared up to z~2. We present the prediction in the Spitzer 24-micron band of a phenomenological model for galaxy evolution derived from the 15-micron data. Without any ``a posteriori'' update, the model predictions seem to agree well with the recently published 24-micron extragalactic source counts, suggesting that the peak in the 24-micron counts is dominated by ``starburst'' galaxies like those detected by ISOCAM at 15 micron, but at higher redshifts (1 < z < 2 instead of 0.5 < z < 1.5).

We make predictions for the cosmological surveys to be conducted by the future Herschel mission operating in the far-infrared. The far-infrared bands match the peak of the CIRB, the brightest background of astrophysical origin. Therefore, surveys in these bands will provide essential information on the evolutionary properties of Luminous and Ultra-Luminous Infrared Galaxies (LIGs and ULIGs), starburst and normal galaxies. Our predictions are based on a new phenomenological model obtained from the 15-micron luminosity function of galaxies and AGN, fitting all the ISOCAM observables (source counts and redshift distributions) and also the recently published Spitzer source counts in the 24-micron band. We discuss the confusion noise due to extragalactic sources, depending strongly on the shape of the source counts and on the telescope parameters. We derive the fraction of the CIRB expected to be resolved by Herschel in the different wavebands and we discuss extragalactic surveys that could be carried on by Herschel for different scientific puropouses (i.e. ultra-deep, deep and shallow).

We present the 90 micron luminosity function of the Final Analysis of the European Large Area ISO Survey (ELAIS), extending the sample size of our previous analysis (paper IV) by about a factor of 4. Our sample extends to z=1.1, around 50 times the comoving volume of paper IV, and 10^{7.7} < h^{-2}L/Lsun < 10^{12.5}. From our optical spectroscopy campaigns of the northern ELAIS 90 mircon survey (7.4 deg^2 in total, to S(90um)>70mJy), we obtained redshifts for 61% of the sample (151 redshifts) to B<21 identified at 7 microns, 15 microns, 20cm or with bright (B<18.5) optical identifications. The selection function is well-defined, permitting the construction of the 90 micron luminosity function of the Final Analysis catalogue in the ELAIS northern fields, which is in excellent agreement with our Preliminary Analysis luminosity function in the ELAIS S1 field from paper IV. The luminosity function is also in good agreement with the IRAS-based prediction of Serjeant & Harrison (2004), which if correct requires luminosity evolution of (1+z)^{3.4 +/- 1.0} for consistency with the source counts. This implies an evolution in comoving volume averaged star formation rate at z<~1 consistent with that derived from rest-frame optical and ultraviolet surveys.

We present the final band-merged ELAIS catalogue at 6.7, 15, 90, and 175 $\mu$m, and the associated data at u,g,r,i,z,J,H,K, and 20cm. The origin of the survey, infrared and radio observations, data-reduction and optical identifications are briefly reviewed, and a summary of the area covered, and completeness limit for each infrared band is given. A detailed discussion of the band-merging and optical association strategy is given. The total catalogues consists of 2860 sources. For extragalactic sources observed in 3 or more infrared bands, colour-colour diagrams are presented and discussed in terms of the contributing infrared populations. Spectral energy distributions are shown for selected sources and compared with cirrus, M82 and Arp220 starburst, and AGN dust torus models. Spectroscopic redshifts are tabulated, where available. For the N1 and N2 areas, the INT ugriz Wide Field Survey, permits photometric redshifts to be estimated for galaxies and quasars. These agree well with the spectroscopic redshifts, within the uncertainty of the photometric method. The redshift distribution is given for selected ELAIS bands and colour-redshift diagrams are discussed. There is a high proportion of ultraluminous infrared galaxies in the ELAIS Catalogue ($> 10 %$ of 15 $\mu$m sources), many with Arp220-like colours. 10 hyperluminous infrared galaxies and 10 EROs are found in the survey. The large numbers of ultraluminous galaxies imply very strong evolution in the star-formation rate between z = 0 and 1.

We present a multi-wavelength catalog (15 um, R, K-band, 1.4 GHz flux) plus spectroscopic identifications for 406 15 um sources detected in the ELAIS region S1, over the flux density range 0.5<S<150 mJy. 332 (~82%) sources are optically identified down to R~23.0. Spectra or bona fide stellar identifications are obtained for 290 objects. The areal coverage, MIR and optical completeness of the sample are discussed in order to allow statistical and evolutionary analyses. Two main spectroscopic classes have been found to dominate the MIR extragalactic population: z<0.5 star-forming galaxies (from absorbed to extreme starbursts: nuL_nu~10^8-10^11 L_odot), which account for \~75% of the sources, and AGN (both type 1 and 2), which account for ~25% of the sources. About 20% of the extragalactic sources are dust-enshrouded starburst galaxies [e(a) spectra], and all the starburst galaxies appear more dust extincted in the optical than nearby normal galaxies. We also identified 91 stellar objects (~22% of the MIR sources). The counts for starburst galaxies and AGN down to 0.6 mJy have been derived. A general trend is found in the optical-MIR SED of the galaxies, where the MIR-luminous objects have larger MIR to optical luminosity ratios.

We present the first determination of the 15micron luminosity functio galaxies from the European Large Area ISO survey (ELAIS) southern fields. We have adopted a new criterion to separate the quiescent, non-evolving and the starburst, evolving populations based on the ratio of mid-infrared to optical luminosities. Strong evolution is suggested by our data for the starburst galaxy population, while normal spiral galaxies are consistent with no evolution. The starburst population must evolve both in luminosity and in density with rates of the order L(z) (1+z)^(3.5) and rho(z) (1+z)^3.8 up to z~1. The evolutionary parameters of our model have been tested by comparing the model predictions with other observables, like source counts at all flux density levels (from 0.1 to 300 mJy) and redshift distributions and luminosity functions at high-z (0.7 < z < 1.0 from HDF-N data). The agreement between our model predictions and the observed data is remarkably good. We use our data to estimate the star-formation density of the Universe up to z=0.4 and we use the luminosity function model to predict the trend of the star-formation history up to z=1.

We present an analysis of the two-point angular correlation function of the ELAIS S1 survey. The survey covers 4 deg$^2$ and contains 462 sources detected at 15$\mu$m to a 5$\sigma$ flux limit of 0.45 mJy. Using the 329 extragalactic sources not repeated in different observations, we detect a significant clustering signal; the resulting angular correlation function can be fitted by a exponential law $w(\theta) = A \theta^{1-\gamma}$ with $A = 0.014\pm0.005$ and $\gamma = 2.04\pm0.18$. Assuming a redshift distribution of the objects, we invert Limber's equation and deduce a spatial correlation length $r_0=4.3^{+0.4}_{-0.7} h^{-1}$ Mpc. This is smaller than that obtained from optical surveys but it is in agreement with results from IRAS. This extends to higher redshift the observational evidence that infrared selected surveys show smaller correlation lengths (i.e. reduced clustering amplitudes) than optical surveys.

We have carried out observations at 18.5 GHz with the Australia Telescope Compact Array (ATCA) of 250 out of the 258 Southern extragalactic sources in the complete 5 GHz 1 Jy sample by K\"uhr et al. (1981). In this paper we focus on the polarization properties of this sample, while other properties will be addressed in a future paper. In our analysis we subdivide the sample into flat and steep spectrum sources following Stickel et al. (1994) classification, where spectral indices were measured between 2.7 and 5 GHz. The polarized flux has been measured with a S/N > 5 for 170 sources (114 flat-spectrum and 56 steep-spectrum) and upper limits have been set for additional 27 sources (12 flat-spectrum and 15 steep-spectrum). The median polarization degree at 18.5 GHz for the flat-spectrum sub-sample is \Pi_{18.5}\simeq 2.7%, about a factor of 2 higher than at 1.4 GHz (\Pi_{1.4}\simeq 1.4%, based on NVSS data). For flat-spectrum sources we find a weak correlation between \Pi_{18.5} and the high frequency (5--18.5 GHz) spectral index. No evidences of significant correlations of the polarization degree with other source properties are found. The median value of \Pi_{18.5} for the steep spectrum sources is \simeq 4.8%, but our sample might be biased against extended sources. We find indications of a correlation between \Pi_{18.5} and both the low frequency (1.4--5 GHz) and the high frequency (5--18.5 GHz) spectral indices. An important application of this work is the possibility to estimate the contamination of CMB polarization maps by extragalactic radio sources. Our results indicate that such contamination is within the range of estimates given by Mesa et al. (2002).

We have studied the optical, near-IR and radio properties of a complete sample of 43 sources detected at 15-micron in one of the deeper ELAIS repeatedly observed region. The extragalactic objects in this sample have 15-micron flux densities in the range 0.4-10 mJy, where the source counts start diverging from no evolution models. About 90% of the sources (39 out of 43) have optical counterparts brighter than I=21 mag. Eight of these 39 sources have been identified with stars on the basis of imaging data, while for another 22 sources we have obtained optical spectroscopy, reaching a high identification percentage (30/43, ~70%). All but one of the 28 sources with flux density > 0.7 mJy are identified. Most of the extragalactic objects are normal spiral or starburst galaxies at moderate redshift (z_med~0.2); four objects are Active Galactic Nuclei. We have used the 15-micron, H_alpha and 1.4-GHz luminosities as indicators of star-formation rate and we have compared the results obtained in these three bands. While 1.4-GHz and 15-micron estimates are in good agreement, showing that our galaxies are forming stars at a median rate of ~40 Mo/yr, the raw H_alpha-based estimates are a factor ~5-10 lower and need a mean correction of ~2 mag to be brought on the same scale as the other two indicators. A correction of ~2 mag is consistent with what suggested by the Balmer decrements H_alpha/H_beta and by the optical colours. Moreover, it is intermediate between the correction found locally for normal spirals and the correction needed for high-luminosity 15-micron objects, suggesting that the average extinction suffered by galaxies increases with infrared luminosity.

The radio counterparts to the 15-micron sources in the European Large Area ISO Survey southern fields are identified in 1.4-GHz maps down to about 80 microJy. The radio-MIR correlation is investigated and derived for the first time at these flux densities for a sample of this size. Our results show that radio and MIR luminosities correlate almost as well as radio and FIR, at least up to z ~ 0.6. Using the derived relation and its spread together with the observed 15-micron counts, we have estimated the expected contribution of the 15-micron extragalactic populations to the radio source counts and the role of MIR starburst galaxies in the well known 1.4-GHz source excess observed at sub-mJy levels. Our analysis demonstrates that IR emitting starburst galaxies do not contribute significantly to the 1.4-GHz counts for strong sources, but start to become a significant fraction of the radio source population at flux densities < 0.5 - 0.8 mJy. They are expected to be responsible for more than 60% of the observed radio counts at < 0.05 mJy. These results are in agreement with the existing results on optical identifications of faint radio sources.

The techniques and software tools developed for the reduction and analysis of ISO-CAM/PHOT data with the LARI method are presented. The method, designed for the detection of faint sources in ISO raster observations, is based on the assumption of the existence of two different time scales in the detectors' transient behaviour, accounting either for fast or slow detectors' response. The specifically developed IDL software includes: a reduction pipeline performing basic operations such as deglitching and background determination; the fitting procedures proper, modelling the time history of individual pixels and detecting any flux excess with respect to the local background ascribable to potential sources; mapping, source extraction and flux estimation procedures; simulation procedures allowing one to estimate the errors arising from different instrumental and reduction effects. Moreover, an easy-to-use graphical user interface allows one to quickly browse the data and carry out the substantial amount of interactive analysis required when the automatic fit fails and to check the reliability of detected sources. This method provides source lists of great reliability and completeness and an outstanding photometric accuracy, particularly at low redundancy levels, where the reliability of ISO-CAM/PHOT source lists, even at moderately bright flux levels, has been a long-standing issue. In this work a description of the techniques and software tools that have been developed is given, alongside with some highlights from the results obtained thanks to their application to different fields.

We present J and K-band near-infrared photometry of a sample of mid-infrared sources detected by the Infrared Space Observatory (ISO) as part of the European Large Area ISO-Survey (ELAIS) and study their classification and star-forming properties. We have used the Preliminary ELAIS Catalogue for the 6.7 micron (LW2) and 15 micron (LW3) fluxes. All of the high-reliability LW2 sources and 80 per cent of the LW3 sources are identified in the near-IR survey reaching K = 17.5 mag. The near- to mid-IR flux ratios can effectively be used to separate stars from galaxies in mid-IR surveys. The stars detected in our survey region are used to derive a new accurate calibration for the ELAIS ISOCAM data in both the LW2 and LW3 filters. We show that near to mid-IR colour-colour diagrams can be used to further classify galaxies, as well as study star-formation. The ISOCAM ELAIS survey is found to mostly detect strongly star-forming late-type galaxies, possibly starburst powered galaxies, and it also picks out obscured AGN. The ELAIS galaxies yield an average mid-IR flux ratio LW2/LW3 = 0.67 +/- 0.27. We discuss this [6.7/15] ratio as a star formation tracer using ISO and IRAS data of a local comparison sample. We find that the [K/15] ratio is also agood indicator of activity level in galaxies and conclude that the drop in the [6.7/15] ratio seen in strongly star-forming galaxies is a result of both an increase of 15 mic emission and an apparent depletion of 6.7 mic emission. Near-IR data together with the mid-IR give the possibility to estimate the relative amount of interstellar matter in the galaxies.

We present the 15 micron extragalactic source counts from the Final Analysis Catalogue of the European Large Area ISO Survey southern hemisphere field S1, extracted using the Lari method. The large number of extragalactic sources (about 350) detected over this area between 0.5 and 100 mJy guarantee a high statistical significance of the source counts in the previously poorly covered flux density range between IRAS and the Deep ISOCAM Surveys. The bright counts in S1 (> 2mJy) are consistent with a flat, Euclidean slope, suggesting the dominance of a non-evolving population. In contrast, at fainter fluxes (< 2 mJy) our counts show a strong departure from no evolution models, with a very steep super-Euclidean slope down to our flux limit. Strong luminosity and density evolution is needed at least for the population of star-forming galaxies, together with a luminosity break around 10^10.8 Lsun in their local luminosity function in order to reproduce our sharp increase of the counts below 2 mJy and the redshift distributions observed for 15 micron sources at different flux levels. The contribution of the strongly evolving starburst population (down to 50 microJy) to the 15 micron cosmic background is estimated to be 2.2 nW m^-2 sr^-1 (67% of the total mid-infrared background estimate).

We present the 15 micron luminosity function of type 1 AGN (QSO + Seyfert 1). Our sample of 21 high-redshift sources is selected from the Preliminary Analysis catalogue in the S1 field of the European Large Area ISO Survey (ELAIS). To study the cosmological evolution of the AGN1 luminosity function, our sample has been combined with a local sample of 41 sources observed by IRAS. We find that the luminosity function of AGN1 at 15 micron is fairly well represented by a double-power-law-function. There is evidence for significant cosmological evolution consistent with a pure luminosity evolution (PLE) model L(z) (1+z)^k, with k=3.0-3.3. The value of k depends on the existence or not of an evolution cut-off at redshift ~2, and on the adopted cosmology. From the luminosity function and its evolution we estimate a contribution of AGN1 to the Cosmic Infrared Background (CIRB) of nuI_nu ~ 6 x 10^{-11}W m^{-2} sr^{-1} at 15 micron. This corresponds to ~2-3% of the total observed CIRB at this wavelength. Under the usual assumptions of unified models for AGN, the expected contribution of the whole AGN population to the CIRB at 15 micron is 10-15%.

We present a new determination of the local volume-averaged star formation rate from the 1.4 GHz luminosity function of star forming galaxies. Our sample, taken from the B<=12 Revised Shapley-Ames catalogue (231 normal spiral galaxies over effective area 7.1 sr) has ~100% complete radio detections and is insensitive to dust obscuration and cirrus contamination. After removal of known active galaxies, the best-fit Schechter function has a faint-end slope of -1.27+/-0.07 in agreement with the local Halpha luminosity function, characteristic luminosity L*=(2.6+/-0.7)*10^{22} W/Hz and density phi* = (4.8 +/-1.1)*10^{-4} / Mpc^3. The inferred local radio luminosity density of (1.73+/-0.37+/-0.03)*10^{19} W/Hz/Mpc^3 (Poisson noise, large scale structure fluctuations) implies a volume averaged star formation rate ~2 x larger than the Gallego et al. Halpha estimate, i.e. rho(1.4 GHz} = (2.10+/-0.45+/-0.04) *10^{-2}$ Msun/yr/Mpc^3 for a Salpeter initial mass function from 0.1-125 Msun and Hubble constant of 50 km/s/Mpc. We demonstrate that the Balmer decrement is a highly unreliable extinction estimator, and argue that optical-UV SFRs are easily underestimated, particularly at high redshift.

We present the estimate of the evolution of type 1 AGNs in the hard (2-10 keV) X-rays drawn from the HELLAS survey, and in the IR (15um) obtained from the ELAIS survey. We find that the local luminosity function (LF) of AGN1 in the 2-10 keV band is fairly well represented by a double-power-law-function. There is evidence for significant cosmological evolution according to a pure luminosity evolution model Lx(z)~(1+z)^k, with k=2.12 and k=2.19 in a (Omega_m,Omega_lambda)=(1.0,0.0) and in a (Omega_m,Omega_lambda)=(0.3,0.7) cosmology respectively. In a (Omega_m,Omega_lambda)=(1.0,0.0) Universe the data show an excess of faint high redshift type 1 AGN which is well modeled by a luminosity dependent density evolution, similarly to what observed in the soft X-rays. In the IR band, with a (Omega_m,Omega_lambda)=(1.0,0.0) cosmology, the evolution found is similar to what observed in other wavebands, the LF is a double-power-law-function with a bright slope 2.9 and a faint slope 1.1, following a luminosity evolution model L(z)~(1+z)^3.0.

We have developed a new data reduction technique for ISOCAM LW data and have applied it to the European Large Area ISO Survey (ELAIS) LW3 (15 micron) observations in the southern hemisphere (S1). This method, known as LARI technique and based on the assumption of the existence of two different time scales in ISOCAM transients (accounting either for fast or slow detector response), was particularly designed for the detection of faint sources. In the ELAIS S1 field we obtained a catalogue of 462 15 micron sources with signal-to-noise ratio >= 5 and flux densities in the range 0.45 - 150 mJy (filling the whole flux range between the Deep ISOCAM Surveys and the IRAS Faint Source Survey). The completeness at different flux levels and the photometric accuracy of this catalogue have been tested with simulations. Here we present a detailed description of the method and discuss the results obtained by its application to the S1 LW3 data.

We present BeppoSAX observations of the Southern S1 region in the European Large Area ISO Survey (ELAIS). These observations cover an area of ~1.7 sq. deg. and reach an on-axis (~0.7 sq. deg) 2-10 keV (HX) sensitivity of ~10E-13 cgs. This is the first HX analysis of an ISOCAM survey. We detect 9 sources with a signal to noise ratio SNR_HX>3, 4 additional sources with a 1.3 to 10 keV (T) SNR_T>3 and 2 additional sources which appear to be associated with QSOs with SNR_T>2.9. The number densities of the SNR_HX>3 sources are consistent with the ASCA and BeppoSAX logN-logS functions. Six BeppoSAX sources have reliable ISOCAM 15 micron counterparts within ~60 arcsec. All these ISOCAM sources have optical counterparts of R<20 mags. Five of these sources have been previously optically classified giving 4 QSOs and 1 BALQSO at z=2.2. The remaining unclassified source has X-ray and photometric properties consistent with that of a nearby Seyfert galaxy. One further HX source has a 2.6sigma ISOCAM counterpart associated with a galaxy at z=0.325. If this ISOCAM source is real, the HX/MIR properties suggest either an unusual QSO or a cD cluster galaxy. We have constructed MIR and HX spectral energy distributions to compute the expected HX/MIR ratios for these classes of objects up to z=3.2 and assess the HX/MIR survey depth. The BALQSO has an observed X-ray softness ratio and HX/MIR flux ratio similar to QSOs but different to those found for low redshift BALQSOs. This difference can be explained in terms of absorption, and suggests that high redshift BALQSOs should be comparatively easy to detect in the HX band, allowing their true fraction in the high redshift QSO population to be determined (ABRIDGED).

We present the luminosity function of 90um selected galaxies from the European Large Area ISO Survey (ELAIS), extending to z=0.3. Their luminosities are in the range 10^9 < h_65^-2 L/Lsun < 10^12, i.e. non-ultraluminous. From our sample of 37 reliably detected galaxies in the ELAIS S1 region from the Efstathiou et al. (2000) S_90 >= 100mJy database, we found optical, 15um or 1.4GHz identifications for 24 (65%). We have obtained 2dF and UK Schmidt FLAIR spectroscopy of 89% of IDs to rigid multivariate flux limits. We construct a luminosity function assuming (a) our spectroscopic subset is an unbiased sparse sample, and (b) there are no galaxies which would not be represented in our spectroscopic sample at {\it any} redshift. We argue that we can be confident of both assumptions. We find the luminosity function is well-described by the local 100um luminosity function of Rowan-Robinson, Helou & Walker (1987). {\it Assuming} this local normalisation, we derive luminosity evolution of (1+z)^{2.45\pm0.85} (95% confidence). We argue that star formation dominates the bolometric luminosities of these galaxies and we derive comoving star formation rates in broad agreement with the Flores et al. (1999) and Rowan-Robinson et al. (1997) mid-IR-based estimates.

We present the first results of our HST WFPC2 F814W snapshot imaging survey, targeting virtually all sub-mJy decimetric radio-selected star-forming galaxies. The radio selection at ~1 GHz is free from extinction effects and the radio luminosities are largely unaffected by AGN contamination, making these galaxies ideal tracers of the cosmic star formation history. A sub-sample of 4 targets is presented here, selected at 1.4 GHz from the spectroscopically homogenous and complete samples of Benn et al. (1993) and Hopkins et al. (1999). The redshifts are confined to a narrow range around z~0.2, to avoid differential evolution, with a radio luminosity close to L* where the galaxies dominate the comoving volume-averaged star formation rate. We find clearly disturbed morphologies resembling those of ultraluminous infrared galaxies, indicating that galaxy interactions may be the dominant mechanism for triggering star formation at these epochs. The morphologies are also clearly different from coeval quasars and radiogalaxies, as found in star-forming galaxies selected at other wavelengths. This may prove challenging for models which propose direct causal links between AGN evolution and the cosmic star formation history at these epochs. The asymmetries are typically much larger than seen in the CFRS at similar redshifts, optical luminosities and Halpha-derived star formation rates, indicating the possible existence of an obscuration-related morphological bias in such samples.

We present the first preliminary estimate of the evolution of type 1 AGNs at 15 um. A new sample of sources selected from the ELAIS-S1 survey has been used together with the RMS sample. The ELAIS-S1 survey covers an effective area of 2.2 deg^2 (total: 4 deg^2) down to 1.2 mJy at 15 um. Spectroscopic identifications have been carried out in the range 17.2<R<20.0 at the 3.6m/ESO, NTT/ESO and 2dF/AAT. The luminosity function is steeper at bright luminosities in comparison to the previous estimates of Rush, Malkan & Spinoglio (1993). The evolution is compatible with a Pure Luminosity Evolution model with L(z)=L(0)(1+z)^3.4.

We present some recent results from the European Large Area ISO Survey (ELAIS). This survey was the largest non-serendipitous ISO field survey. A preliminary reduction has recently been completed and catalogues of sources released to the community. Early results show strongly evolving source counts. A comprehensive identification programme is underway and a number of extremely luminous objects have already been discovered. This survey provides an exciting legacy from the ISO mission and (amongst many goals) will allow us to provide important constraints on the obscured star-formation history of the Universe.

We present preliminary source counts at 6.7um and 15um from the Preliminary Analysis of the European Large Area ISO survey, with limiting flux densities of \~2mJy at 15um & ~1mJy at 6.7um. We separate the stellar contribution from the extragalactic using identifications with APM sources made with the likelihood ratio technique. We quantify the completeness & reliability of our source extraction using (a) repeated observations over small areas, (b) cross-IDs with stars of known spectral type, (c) detections of the PSF wings around bright sources, (d) comparison with independent algorithms. Flux calibration at 15um was performed using stellar IDs; the calibration does not agree with the pre-flight estimates, probably due to effects of detector hysteresis and photometric aperture correction. The 6.7um extragalactic counts are broadly reproduced in the Pearson & Rowan-Robinson model, but the Franceschini et al. (1997) model underpredicts the observed source density by ~0.5-1 dex, though the photometry at 6.7um is still preliminary. At 15um the extragalactic counts are in excellent agreement with the predictions of the Pearson & Rowan-Robinson (1996), Franceschini et al. (1994), Guiderdoni et al. (1997) and the evolving models of Xu et al. (1998), over 7 orders of magnitude in 15um flux density. The counts agree with other estimates from the ISOCAM instrument at overlapping flux densities (Elbaz et al. 1999), provided a consistent flux calibration is used. Luminosity evolution at a rate of (1+z)^3, incorporating mid-IR spectral features, provides a better fit to the 15um differential counts than (1+z)^4 density evolution. No-evolution models are excluded, and implying that below around 10mJy at 15um the source counts become dominated by an evolving cosmological population of dust-shrouded starbursts and/or active galaxies.

We have derived the relationship between polarization and intensity fluctuations due to point sources. In the case of a Poisson distribution of a population with uniform evolution properties and constant polarization degree, polarization fluctuations are simply equal to intensity fluctuations times the average polarization degree. Conservative estimates of the polarization degree of the classes of extragalactic sources contributing to fluctuations in the frequency ranges covered by the forthcoming space missions MAP and Planck Surveyor indicate that extragalactic sources will not be a strong limiting factor to measurements of the polarization of the Cosmic Microwave Background.

The European Large Area ISO Survey (ELAIS) has surveyed 12 square degrees of the sky at 15 and 90 microns, and subsets of this area at 6.75 and 175 microns, using the Infrared Space Observatory (ISO). This project was the largest single open time programme executed by ISO, taking 375 hours of data. A preliminary catalogue of more than 1000 galaxies has been produced. In this talk we describe the goals of the project, describe the follow-up programmes that are in progress, and present some first scientific results including a provisional number count analysis at 15 and 90 microns.

We present the X-ray data and the optical identifications for a deep ROSAT PSPC observation in the "Marano field". In the inner region of the ROSAT field (15' radius) we detected 50 X-ray sources with Sx >= 3.7x10^(-15) erg/cm^2/s. When corrected for the different sensitivity over the field, the estimated observed surface density at Sx >= 4x10^(-15) erg/cm^2/s is 272+/-40 sources/sq.deg. Four X-ray sources, corresponding to 8% of the total sample, have been detected in radio images with a flux limit of about 0.2 mJy. Careful statistical analysis of multicolour CCD data in the error boxes of the 50 X-ray sources has led to the identification of 42 sources, corresponding to 84% of the X-ray sample. These 42 reliable identifications are 33 AGNs (including two radio galaxies and one BL Lac candidate; 79% of the identified sources), 2 galaxies, 3 groups or clusters of galaxies and 4 stars. We also show that it is likely that a few of the 8 unidentified sources are such because the derived X-ray positions may be offset with respect to the real ones due to confusion effects. The unidentified sources have a large ratio of X-ray to optical fluxes and most of them have harder than average X-ray spectra. Since most of the identified objects with these characteristics in our field and in the Lockman field are AGNs, we conclude that also most of these sources are likely to be AGNs. Finally, comparing the optically and X-ray selected samples of AGNs in this field, we estimate the ``efficiency'' of AGN selection with deep X-ray exposures and with purely optical criteria.

A deep survey of the European Large Area ISO Survey (ELAIS) field in the southern celestial hemisphere (hereinafter S1) has been carried out with the Australia Telescope Compact Array (ATCA) at 1.4 GHz. The S1 region, covering about 4 square degrees, has been also surveyed in the mid- and far-infrared (5-200 micron) with the Infrared Space Observatory (ISO). The radio survey provides uniform coverage of the entire S1 region, with a sensitivity (5 sigma) of 0.4 mJy over the whole area and 0.2 mJy in the center. To this sensitivity, virtually all the radio counterparts of the far-infrared extragalactic ISO sources should be detected. This paper presents a radio sample - complete at the 5 sigma level - consisting of 581 sources detected at 1.4 GHz. Of these radio sources, 349 have peak flux density in the range 0.2-1 mJy, forming a new homogeneous sample of sub-mJy radio sources. Due to its size, depth and multi-waveband coverage, the sample will allow us to study in greater detail the sub-mJy radio source population. The full catalogue will be available from http://athena.ph.ac.uk/

Deep imaging and spectroscopy have been carried out for optical counterparts of a sample of 68 faint radio sources (S > 0.2 mJy) in the ``Marano Field''. About 60% of the sources have been optically identified on deep CCD exposures (limit R ~ 24.0) or ESO 3.6-m plates (limit bJ ~ 22.5). Thirty-four spectra (50% of the total radio sample) were obtained with the ESO 3.6-m telescope and 30 redshifts were determined. In addition to a few broad line active galactic nuclei, three main spectroscopic classes have been found to dominate the faint radio galaxy population: (1) Early-type galaxies with 0.1 < z < 0.8. (2) Late-type galaxies at moderate redshift (z < 0.4), with relatively bright magnitudes (B < 22.5) and sub-mJy radio fluxes. (3) A group of bright high-redshift (z > 0.8) radio galaxies with moderate-to-strong [OII] emission. They have spectra, colours and absolute magnitudes similar to those of the classical bright elliptical radio galaxies found in surveys carried out at higher radio fluxes. Star-forming galaxies do not constitute the main population of our radio sources identified with galaxies. In fact, even at sub-mJy level the majority of our radio sources are identified with early-type galaxies. This apparent discrepancy with previous results is due to the fainter magnitude limit reached in our spectroscopic identifications. Moreover, using mainly the large radio-to-optical ratio and the information from the available limits on the optical magnitudes of the unidentified radio sources, we conclude that the great majority of them are likely to be early-type galaxies, at z > 1. If correct, it would suggest that the evolution of the radio luminosity function of spiral galaxies, including starbursts, might not be as strong as suggested in previous evolutionary models.

We present a new determination of the co-moving star formation density at redshifts z~<0.35 from the 1.4 GHz luminosity function of sub-mJy star-burst galaxies. Our sample, taken from Benn et al. (1993), is insensitive to dust obscuration. The shape of the Luminosity function of this sample is indistinguishable from a number of reasonable a prior models of the luminosity function. Using these shapes we calculate the modest corrections (typically ~< 20 per cent to the observed 1.4 GHz luminosity density. We find that the cosmic variance in our estimate of this luminosity density is large. We find a luminosity density in broad agreement with that from the RSA sample by Condon et al. (1987). We infer a co-moving star formation rate surprisingly similar to coeval estimates from the Canada-France Redshift Survey, in both ultraviolet and Halpha, although the later may also be affected by Cosmic variance. We conclude that the intermediate 0.05<z<0.3 star-formation rate is not yet well determined partially due to uncertain extinction corrections, but also partially due to cosmic variance. We suggest that deep moderate area radio surveys will improve this situation considerably.

We have used the Very Large Array(VLA) in C configuration to carry out a sensitive 20cm radio survey of regions of sky that have been surveyed in the Far Infra-Red over the wavelength range 5-200 microns with ISO as part of the European Large Area ISO Survey(ELAIS). As usual in surveys based on a relatively small number of overlapping VLA pointings the flux limit varies over the area surveyed. The survey has a flux limit that varies from a 5$\sigma$ limit of 0.135mJy over an area of 0.12deg$^2$ to a 5$\sigma$ limit of 1.15mJy or better over the whole region covered of 4.22 deg$^2$. In this paper we present the radio catalogue of 867 sources. These regions of sky have previously been surveyed to shallow flux limits at 20cm with the VLA as part of the VLA D configuration NVSS(FWHM=45 arcsec) and VLA B configuration FIRST(FWHM=5 arcsec) surveys. We have carried out a a detailed comparison of the reliability of our own survey and these two independent surveys in order to assess the reliability and completeness of each survey.

Radio surveys with the Australia Telescope Compact Array have been carried out at 1.4 and 2.4 GHz with a limiting flux of 0.2 mJy at each frequency in the Marano Field, in which deep optical and X-ray (ROSAT) data are also available. In this paper we present the two radio samples, complete at the 5 sigma local level, consisting of 63 and 48 sources respectively at 1.4 and 2.4 GHz. The 1.4 GHz normalized differential source counts show a flattening below about one mJy, in agreement with the results from previous surveys. The 2.4 GHz counts, which are the deepest at this or similar frequencies (e.g. 2.7 GHz), agree well with the 2.7 GHz counts at higher fluxes and with the extrapolations down to 2 mJy based on fluctuation analyses. The 2.4-1.4 GHz spectral index distributions are presented for both the complete samples in two flux density intervals. The median spectral index for the 1.4 GHz sample remains constant at about 0.8, down to the lowest fluxes (S(1.4)=0.2 mJy), while for the higher frequency sample the spectral index distribution flattens in the lower flux density interval (S(2.4) < 0.8 mJy). A significant number of sources with inverted spectrum (< 0) does appear in both samples, at low flux level (< 2 mJy). These sources, which are about 25% of the complete sample at 2.4 GHz, are probably the ``bright'' counterpart of the inverted spectrum sources which appear to be almost 50% of the sources in the even deeper radio surveys (at about 20-40 microJy).