In the second work of this series, we analyse the connection between the availability of gas and the position of a region with respect to the spatially resolved main sequence (MS) relation. Following the procedure presented in Paper I we obtain 500pc scales estimates of stellar mass and star formation rate surface densities ($\Sigma_{\star}$ and $\Sigma_{\rm{SFR}}$). Our sample consists of five face-on, grand design spiral galaxies located on the MS. Thanks to HI 21cm and $^{12}$CO(2-1) maps, we connect the gas surface densities and gas fractions to the observed star formation properties of each region. We find that the spatially resolved MS ($\sigma=0.23$ dex) is the combination of two relations: the Kennicutt-Schmidt law ($\sigma=0.19$ dex) and the molecular gas MS (MGMS, $\sigma=0.22$ dex); $\Sigma_{\star}$, $\Sigma_{\rm{SFR}}$ and the surface density of the molecular gas, $\Sigma_{\rm{H_2}}$, define a 3D relation as proposed by \citet{2019ApJ...884L..33L}. We find that $\Sigma_{\rm{H_2}}$ steadily increases along the MS relation, varies little towards higher $\Sigma_{\rm{SFR}}$ at fixed stellar surface densities (not enough to sustain the change in SFR), and it is almost constant perpendicular to the relation. The surface density of neutral gas ($\Sigma_{\rm{HI}}$) is constant along the MS, and increases in its upper envelop. $\Sigma_{\rm{SFR}}$ can be expressed as a function of $\Sigma_{\star}$ and $\Sigma_{\rm{HI}}$, following the Equation: $\log\Sigma_{\rm{SFR}}$ = 0.97$\log\Sigma_{\star}$ + 1.99$\log\Sigma_{\rm{HI}}$ - 11.11. Finally, we show that f$_{\rm{gas}}$ increases significantly towards the starburst region in the $\log\Sigma_{\star}$ - $\log\Sigma_{\rm{SFR}}$ plane, accompanied by a slight increase in SFE.

We analyse the spatially resolved relation between stellar mass (M$_{\star}$) and star formation rate (SFR) in disk galaxies (i.e. the Main Sequence, MS). The studied sample includes eight nearby face-on grand-design spirals, e.g. the descendant of high-redshift, rotationally-supported star-forming galaxies. We exploit photometric information over 23 bands, from the UV to the far-IR, from the publicly available DustPedia database to build spatially resolved maps of stellar mass and star formation rates on sub-galactic scales of 0.5-1.5 kpc, by performing a spectral energy distribution fitting procedure that accounts for both the observed and the obscured star formation processes, over a wide range of internal galaxy environments (bulges, spiral arms, outskirts). With more than 30 thousands physical cells, we have derived a definition of the local spatially resolved MS per unit area for disks, $\log(\Sigma_{SFR})$=0.82log$(\Sigma_{*})$-8.69. This is consistent with the bulk of recent results based on optical IFU, using the H$\alpha$ line emission as a SFR tracer. Our work extends the analysis at lower sensitivities in both M$_{\star}$ and SFR surface densities, up to a factor $\sim$ 10. The self consistency of the MS relation over different spatial scales, from sub-galactic to galactic, as well as with a rescaled correlation obtained for high redshift galaxies, clearly proves its universality.

We investigate the nature of star-forming galaxies with reduced specific star formation rate (sSFR) and high stellar masses, those `green valley' objects that seemingly cause a reported bending, or flattening, of the star-forming main sequence. The fact that such objects host large bulges recently led some to suggest that the internal formation of bulges was a late event that induced the sSFRs of massive galaxies to drop in a slow downfall, and thus the main sequence to bend. We have studied in detail a sample of 10 galaxies at $0.45<z<1$ with secure SFR from Herschel, deep Keck optical spectroscopy, and HST imaging from CANDELS allowing us to perform multi-wavelength bulge to disc decomposition, and to derive star formation histories for the separated bulge and disc components. We find that the bulges hosted in these systems below main sequence are virtually all maximally old, with ages approaching the age of the Universe at the time of observation, while discs are young ($\langle$ T$_{50}\rangle \sim 1.5$ Gyr). We conclude that, at least based on our sample, the bending of the main sequence is, for a major part, due to rejuvenation, and we disfavour mechanisms that postulate the internal formation of bulges at late times. The very old stellar ages of our bulges suggest a number density of Early Type Galaxies at $z=1-3$ higher than actually observed. If confirmed, this might represent one of the first direct validations of hierarchical assembly of bulges at high redshifts.

Extreme high-energy peaked BL Lac objects (EHBLs) are an emerging class of blazars with exceptional spectral properties. In blazars, the spectral energy distribution (SED) is dominated by the non-thermal emission of the relativistic jet, and consists of two main broad humps. For the EHBLs, these two components peak in the X-ray and GeV-TeV bands, respectively. Although the number of TeV detected extreme blazars is very limited, recent observations by Imaging Atmospheric Cherenkov Telescopes (IACTs) have revealed that in some of them the energy of the second peak exceeds several TeV (e.g. 1ES 0229+200). Their exceptional hard TeV spectra represent a challenge for the standard leptonic modeling, and a possible hadronic contribution may make these objects high-energy neutrinos producers. Moreover, they are important for the implications on the indirect measurements of the extragalactic background light and of the intergalactic magnetic field. In this contribution, we perform a comparative study of the multi-wavelength spectral energy distributions of a sample of hard X-ray selected EHBL objects. The analysis suggests that the EHBL class is not homogeneous, and a possible sub-classification may be unveiled with TeV gamma-ray observations of the candidates. With the purpose of increasing their number and settle their statistics, we discuss the potential detectability of the currently undetected TeV-emitting EHBLs in our sample by current and next generation of IACTs.

Context. Measurements of the Extragalactic Background Light (EBL) are a fundamental source of information on the collective emission of cosmic sources. Aims. At infrared wavelengths, however, these measurements are precluded by the overwhelming dominance from Interplanetary Dust emission and the Galactic infrared foreground. Only at $\lambda > 300 \ \mu$m, where the foregrounds are minimal, has the Infrared EBL (IR EBL) been inferred from analysis of the COBE maps. The present paper aims to assess the possibility of evaluating the IR EBL from a few $\mu$m up to the peak of the emission at >100 $\mu$m using an indirect method that avoids the foreground problem. Methods. To this purpose we exploit the effect of pair-production from gamma-gamma interaction by considering the highest energy photons emitted by extragalactic sources and their interaction with the IR EBL photons. We simulate observations of a variety of low redshift emitters with the forthcoming Imaging Atmospheric Cherenkov Telescope (IACT) arrays (CTA in particular) and water Cherenkov observatories (LHAASO, HAWC, SWGO) to assess their suitability to constrain the EBL at such long wavelengths. Results. We find that, even under the most extremely favorable conditions of huge emission flares, extremely high-energy emitting blazars are not very useful for our purpose because they are much too distant (>100 Mpc the nearest ones, MKN 501 and MKN 421). Observations of more local Very High Energy (VHE) emitting AGNs, like low-redshift radio galaxies (M87, IC 310, Centaurus A), are better suited and will potentially allow us to constrain the EBL up to $\lambda \simeq 100\ \mu$m.

Extreme high-energy peaked BL Lac objects (EHBLs) are an emerging class of blazars with exceptional spectral properties. The non-thermal emission of the relativistic jet peaks in the spectral energy distribution (SED) plot with the synchrotron emission in X-rays and with the gamma-ray emission in the TeV range or above. These high photon energies may represent a challenge for the standard modeling of these sources. They are important for the implications on the indirect measurements of the extragalactic background light, the intergalactic magnetic field estimate, and the possible origin of extragalactic high-energy neutrinos. In this paper, we perform a comparative study of the multi-wavelength spectra of 32 EHBL objects detected by the Swift-BAT telescope in the hard X-ray band and by the Fermi-LAT telescope in the high-energy gamma-ray band. The source sample presents uniform spectral properties in the broad-band SEDs, except for the TeV gamma-ray band where an interesting bimodality seems to emerge. This suggests that the EHBL class is not homogeneous, and a possible sub-classification of the EHBLs may be unveiled. Furthermore, in order to increase the number of EHBLs and settle their statistics, we discuss the potential detectability of the 14 currently TeV gamma-ray undetected sources in our sample by the Cherenkov telescopes.

Based on HST and MUSE data, we probe the stellar and gas properties (i.e. kinematics, stellar mass, star formation rate) of the radio-loud brightest cluster galaxy (BCG) located at the centre of the X-ray luminous cool core cluster Abell 2667 (z = 0.2343). The bi-dimensional modelling of the BCG surface brightness profile reveals the presence of a complex system of substructures extending all around the galaxy. Clumps of different size and shape plunged into a more diffuse component constitute these substructures, whose intense 'blue' optical colour hints to the presence of a young stellar population. Our results depict the BCG as a massive (M_star ~ 1.38 x 10^11 M_sun) dispersion-supported spheroid (v_star < 150 km/s, sigma_0 ~ 216 km/s) hosting an active supermassive black hole (M_SMBH ~ 3.8 x 10^9 M_sun) whose optical features are typical of low ionisation nuclear emission line regions. Although the velocity pattern of the stars in the BCG is irregular, the stellar kinematics in the regions of the clumps show a positive velocity of ~ 100 km/s, similarly to the gas component. An analysis of the mechanism giving rise to the observed lines in the clumps through empirical diagnostic diagrams points out that the emission is composite, suggesting the contribution from both star formation and AGN. We conclude our analysis describing how scenarios of both chaotic cold accretion and merging with a gas-rich disc galaxy can efficaciously explain the phenomena the BCG is undergoing.

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 quantify the star formation (SF) in the inner cores ($\mathcal{R}$/$R_{200}$$\leq$0.3) of 24 massive galaxy clusters at 0.2$\lesssim$$z$$\lesssim$0.9 observed by the $Herschel$ Lensing Survey and the Cluster Lensing and Supernova survey with $Hubble$. These programmes, covering the rest-frame ultraviolet to far-infrared regimes, allow us to accurately characterize stellar mass-limited ($\mathcal{M}_{*}$$>$$10^{10}$ $M_{\odot}$) samples of star-forming cluster members (not)-detected in the mid- and/or far-infrared. We release the catalogues with the photometry, photometric redshifts, and physical properties of these samples. We also quantify the SF displayed by comparable field samples from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We find that in intermediate-$z$ cluster cores, the SF activity is suppressed with respect the field in terms of both the fraction ($\mathcal{F}$) of star-forming galaxies (SFG) and the rate at which they form stars ($\mathcal{SFR}$ and $s\mathcal{SFR} = \mathcal{SFR}/\mathcal{M}_{*}$). On average, the $\mathcal{F}$ of SFGs is a factor $\sim$$2$ smaller in cluster cores than in the field. Furthermore, SFGs present average $\mathcal{SFR}$ and $s\mathcal{SFR}$ typically $\sim$0.3 dex smaller in the clusters than in the field along the whole redshift range probed. Our results favour long time-scale quenching physical processes as the main driver of SF suppression in the inner cores of clusters since $z$$\sim$0.9, with shorter time-scale processes being very likely responsible for a fraction of the missing SFG population.

We report on optical spectroscopy obtained at the 10.4m Gran Telescopio Canarias of 28 Fermi gamma-ray sources that completes the study of a sample of 60 targets of unidentified objects for which the detection of an X-ray and/or radio source inside the 3FGL error box is available. The observations are aimed to characterize the nature and measure the redshift of these sources. For all optical counterparts, the observations allow us to establish their AGN nature. In particular, we found 24 BL Lac objects, one QSO, one NLSy1 and two objects showing spectral features typical of Seyfert 2 galaxies. For most of them, we determine a spectroscopic redshift, while for five we can set lower limits based on the lack of stellar features from the host galaxy. The global properties of the full sample are briefly discussed.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

For a sample of star forming galaxies in the redshift interval 0.15$<$z$<$0.3, we study how both the relative strength of the AGN infra-red emission, compared to that due to the star formation (SF), and the numerical fraction of AGNs, change as a function of the total stellar mass of the hosting galaxy group (M$^{*}_{\mathrm{group}}$), between $10^{10.25}$ and $10^{11.9}$M$_{\odot}$. Using a multi-component SED fitting analysis, we separate the contribution of stars, AGN torus and star formation to the total emission at different wavelengths. This technique is applied to a new multi-wavelength data-set in the SIMES field (23 not redundant photometric bands), spanning the wavelength range from the UV (GALEX) to the far-IR (Herschel) and including crucial AKARI and WISE mid-IR observations (4.5 \mu m$<\lambda<$24 \mu m), where the BH thermal emission is stronger. This new photometric catalog, that includes our best photo-z estimates, is released through the NASA/IPAC Infrared Science Archive (IRSA). Groups are identified through a friends of friends algorithm ($\sim$62% purity, $\sim$51% completeness). We identified a total of 45 galaxies requiring an AGN emission component, 35 of which in groups and 10 in the field. We find BHAR$\propto ($M$^{*}_{\mathrm{group}})^{1.21\pm0.27}$ and (BHAR/SFR)$\propto ($M$^{*}_{\mathrm{group}})^{1.04\pm0.24}$ while, in the same range of M$^{*}_{\mathrm{group}}$, we do not observe any sensible change in the numerical fraction of AGNs. Our results indicate that the nuclear activity (i.e. the BHAR and the BHAR/SFR ratio) is enhanced when galaxies are located in more massive and richer groups.

We perform lens modelling and source reconstruction of Submillimeter Array (SMA) data for a sample of 12 strongly lensed galaxies selected at 500$\mu$m in the Herschel Astrophysical Terahertz Large Area Survey H-ATLAS. A previous analysis of the same dataset used a single S\`ersic profile to model the light distribution of each background galaxy. Here we model the source brightness distribution with an adaptive pixel scale scheme, extended to work in the Fourier visibility space of interferometry. We also present new SMA observations for seven other candidate lensed galaxies from the H-ATLAS sample. Our derived lens model parameters are in general consistent with previous findings. However, our estimated magnification factors, ranging from 3 to 10, are lower. The discrepancies are observed in particular where the reconstructed source hints at the presence of multiple knots of emission. We define an effective radius of the reconstructed sources based on the area in the source plane where emission is detected above 5$\sigma$. We also fit the reconstructed source surface brightness with an elliptical Gaussian model. We derive a median value $r_{eff}\,\sim 1.77\,$kpc and a median Gaussian full width at half maximum $\sim1.47\,$kpc. After correction for magnification, our sources have intrinsic star formation rates SFR$\,\sim900-3500\,M_{\odot}yr^{-1}$, resulting in a median star formation rate surface density $\Sigma_{SFR}\sim132\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ (or $\sim 218\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ for the Gaussian fit). This is consistent with what observed for other star forming galaxies at similar redshifts, and is significantly below the Eddington limit for a radiation pressure regulated starburst.

We present optical spectroscopy secured at the 10m Gran Telescopio Canarias of the counterparts of 20 extragalactic gamma-ray sources detected by the Fermi satellite. The observations allow us to investigate the nature of these sources and to determine their redshift. We find that all optical counterparts have a spectrum that is consistent with a BL Lac object nature. We are able to determine the redshift for 11 objects and set spectroscopic redshift limits for five targets. Only for four sources the optical spectrum is found featureless. In the latter cases we can set lower limits on the redshift based on the assumption that they are hosted by a typical massive elliptical galaxy whose spectrum is diluted by the non thermal continuum. The observations allow us to unveil the nature of these gamma-ray sources and provide a sanity check of a tool to discover the counterparts of gamma-ray emitters/blazars based on their multiwavelength emission.

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.

In addition to its relevant astrophysical and cosmological significance, the Extragalactic Background Light (EBL) is a fundamental source of opacity for cosmic high energy photons, as well as a limitation for the propagation of high-energy particles in the Universe. We review our previously published determinations of the EBL photon density in the Universe and its evolution with cosmic time, in the light of recent surveys of IR sources at long wavelengths. We exploit deep survey observations by the Herschel Space Observatory and the Spitzer telescope, matched to optical and near-IR photometric and spectroscopic data, to re-estimate number counts and luminosity functions longwards of a few microns, and the contribution of resolved sources to the EBL. These new data indicate slightly lower photon densities in the mid- and far-infrared and sub-millimeter compared to previous determinations. This implies slightly lower cosmic opacity for photon-photon interactions. The new data do not modify previously published EBL modeling in the UV-optical and near-IR up to several microns, while reducing the photon density at longer wavelengths. This improved model of the EBL alleviates some tension that had emerged in the interpretation of the highest-energy TeV observations of local \textit{blazar}s, reducing the case for new physics beyond the standard model (like violations of the Lorenz Invariance, LIV, at the highest particle energies), or for exotic astrophysics, that had sometimes been called for to explain it. Applications of this improved EBL model on current data are considered, as well as perspectives for future instrumentation, the Cherenkov Telescope Array (CTA) in particular.

We discuss a new method for unveiling the possible blazar AGN nature among the numerous population of Unassociated Gamma-ray sources (UGS) in the Fermi catalogues. Our tool relies on positional correspondence of the Fermi object with X-ray sources (mostly from Swift-XRT), correlated with other radio, IR and optical data in the field. We built a set of Spectral Energy Distributions (SED) templates representative of the various blazar classes, and we quantitatively compared them to the observed multi-wavelength flux density data for all Swift-XRT sources found within the Fermi error-box, by taking advantage of some well-recognised regularities in the broad-band spectral properties of the objects. We tested the procedure by comparison with a few well-known blazars, and tested the chance for false positive recognition of UGS sources against known pulsars and other Galactic and extragalactic sources. Based on our spectral recognition tool, we find the blazar candidate counterparts for 14 2FGL UGSs among 183 selected at high galactic latitudes. Further our tool also allows us rough estimates of the redshift for the candidate blazar. In a few cases in which this has been possible (i.e. when the counterpart was a SDSS object), we verified that our estimate is consistent with the measured redshift. The estimated redshifts of the proposed UGS counterparts are larger, on average, than those of known Fermi blazars, a fact that might explain the lack of previous association or identification in published catalogues.

We present rest-frame optical spectra from the FMOS-COSMOS survey of twelve $z \sim 1.6$ \textit{Herschel} starburst galaxies, with Star Formation Rate (SFR) elevated by $\times$8, on average, above the star-forming Main Sequence (MS). Comparing the H$\alpha$ to IR luminosity ratio and the Balmer Decrement we find that the optically-thin regions of the sources contain on average only $\sim 10$ percent of the total SFR whereas $\sim90$ percent comes from an extremely obscured component which is revealed only by far-IR observations and is optically-thick even in H$\alpha$. We measure the [NII]$_{6583}$/H$\alpha$ ratio, suggesting that the less obscured regions have a metal content similar to that of the MS population at the same stellar masses and redshifts. However, our objects appear to be metal-rich outliers from the metallicity-SFR anticorrelation observed at fixed stellar mass for the MS population. The [SII]$_{6732}$/[SII]$_{6717}$ ratio from the average spectrum indicates an electron density $n_{\rm e} \sim 1,100\ \mathrm{cm}^{-3}$, larger than what estimated for MS galaxies but only at the 1.5$\sigma$ level. Our results provide supporting evidence that high-$z$ MS outliers are the analogous of local ULIRGs, and are consistent with a major merger origin for the starburst event.

We present the Spitzer-IRAC/MIPS Extragalactic survey (SIMES) in the South Ecliptic Pole (SEP) field. The large area covered (7.7 deg$^2$), together with one of the lowest Galactic cirrus emissions in the entire sky and a very extensive coverage by Spitzer, Herschel, Akari, and GALEX, make the SIMES field ideal for extragalactic studies. The elongated geometry of the SIMES area ($\approx$4:1), allowing for a significant cosmic variance reduction, further improves the quality of statistical studies in this field. Here we present the reduction and photometric measurements of the Spitzer/IRAC data. The survey reaches a depth of 1.93 and 1.75 $\mu$Jy (1$\sigma$) at 3.6 and 4.5 $\mu$m, respectively. We discuss the multiwavelength IRAC--based catalog, completed with optical, mid-- and far--IR observations. We detect 341,000 sources with F$_{3.6\mu m} \geq 3\sigma$. Of these, 10% have an associated 24 $\mu$m counterpart, while 2.7% have an associated SPIRE source. We release the catalog through the NASA/IPAC Infrared Science Archive (IRSA). Two scientific applications of these IRAC data are presented in this paper: first we compute integral number counts at 3.6 $\mu$m. Second, we use the [3.6]--[4.5] color index to identify galaxy clusters at z$>$1.3. We select 27 clusters in the full area, a result consistent with previous studies at similar depth.

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 used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 {\mu}m at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.

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.

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.

We complement our previous analysis of a sample of z~1-2 luminous and ultra-luminous infrared galaxies ((U)LIRGs), by adding deep VLA radio observations at 1.4 GHz to a large data-set from the far-UV to the sub-mm, including Spitzer and Herschel data. Given the relatively small number of (U)LIRGs in our sample with high S/N radio data, and to extend our study to a different family of galaxies, we also include 6 well sampled near IR-selected BzK galaxies at z~1.5. From our analysis based on the radiative transfer spectral synthesis code GRASIL, we find that, while the IR luminosity may be a biased tracer of the star formation rate (SFR) depending on the age of stars dominating the dust heating, the inclusion of the radio flux offers significantly tighter constraints on SFR. Our predicted SFRs are in good agreement with the estimates based on rest-frame radio luminosity and the Bell (2003) calibration. The extensive spectro-photometric coverage of our sample allows us to set important constraints on the SF history of individual objects. For essentially all galaxies we find evidence for a rather continuous SFR and a peak epoch of SF preceding that of the observation by a few Gyrs. This seems to correspond to a formation redshift of z~5-6. We finally show that our physical analysis may affect the interpretation of the SFR-M* diagram, by possibly shifting, with respect to previous works, the position of the most dust obscured objects to higher M* and lower SFRs.

We present the most complete study to date of the X-ray emission from star-formation in high redshift (median z=0.7; z<1.5), IR-luminous (L_IR=10^10-10^13 L_sun) galaxies detected by Herschel's PACS and SPIRE instruments. For our purpose we take advantage of the deepest X-ray data to date, the Chandra deep fields (North and South). Sources which host AGN are removed from our analysis by means of multiple AGN indicators. We find an AGN fraction of 18+/-2 per cent amongst our sample and note that AGN entirely dominate at values of log[L_X/L_IR]>-3 in both hard and soft X-ray bands. From the sources which are star-formation dominated, only a small fraction are individually X-ray detected and for the bulk of the sample we calculate average X-ray luminosities through stacking. We find an average soft X-ray to infrared ratio of log[L_SX/L_IR]=-4.3 and an average hard X-ray to infrared ratio of log[L_HX/L_IR]=-3.8. We report that the X-ray/IR correlation is approximately linear through the entire range of L_IR and z probed and, although broadly consistent with the local (z<0.1) one, it does display some discrepancies. We suggest that these discrepancies are unlikely to be physical, i.e. due to an intrinsic change in the X-ray properties of star-forming galaxies with cosmic time, as there is no significant evidence for evolution of the L_X/L_IR ratio with redshift. Instead they are possibly due to selection effects and remaining AGN contamination. We also examine whether dust obscuration in the galaxy plays a role in attenuating X-rays from star-formation, by investigating changes in the L_X/L_IR ratio as a function of the average dust temperature. We conclude that X-rays do not suffer any measurable attenuation in the host galaxy.

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.

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

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

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.

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.

In this work we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in these same galaxies. To do this we apply a multi-component, multi-band spectral synthesis technique to a sample of 250 micron selected galaxies of the Herschel Multi-tiered Extragalactic Survey, with IRS spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fit hot dust model (torus). We find a correlation between the obscured star formation rate (SFR) derived from the IR luminosity of the starburst component, SFR_IR and SFR_PAH, derived from the luminosity of the PAH features, L_PAH, with SFR_FIR taking higher values than SFR_PAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of L_PAH to that of the starburst component, L_PAH/L_SB, is almost constant for AGN-dominated objects but decreases with increasing L_SB for starburst-dominated objects. SFR_FIR increases with the accretion luminosity, L_acc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot and cold dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

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.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

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

Dust emission at sub-millimetre wavelengths allows us to trace the early phases of star formation in the Universe. In order to understand the physical processes involved in this mode of star formation, it is essential to gain knowledge about the dark matter structures - most importantly their masses - that sub-millimetre galaxies live in. Here we use the magnification effect of gravitational lensing to determine the average mass and dust content of sub-millimetre galaxies with 250mu flux densities of S_250>15mJy selected using data from the Herschel Multi-tiered Extragalactic Survey. The positions of hundreds of sub-millimetre foreground lenses are cross-correlated with the positions of background Lyman-break galaxies at z~3-5 selected using optical data from the Canada-France Hawaii Telescope Legacy Survey. We detect a cross-correlation signal at the 7-sigma level over a sky area of one square degree, with ~80% of this signal being due to magnification, whereas the remaining ~20% comes from dust extinction. Adopting some simple assumptions for the dark matter and dust profiles and the redshift distribution enables us to estimate the average mass of the halos hosting the sub-millimetre galaxies to be log(M_200/M_sun)=13.17+0.05-0.08(stat.) and their average dust mass fraction (at radii of >10kpc) to be M_dust/M_200~6x10^-5. This supports the picture that sub-millimetre galaxies are dusty, forming stars at a high rate, reside in massive group-sized halos, and are a crucial phase in the assembly and evolution of structure in the Universe.

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $\mu$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $\Delta$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $\Delta$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

The Planck Early Release Compact Source Catalog (ERCSC) has offered the first opportunity to accurately determine the luminosity function of dusty galaxies in the very local Universe (i.e. distances <~ 100 Mpc), at several (sub-)millimetre wavelengths, using blindly selected samples of low redshift sources, unaffected by cosmological evolution. This project, however, requires careful consideration of a variety of issues including the choice of the appropriate flux density measurement, the separation of dusty galaxies from radio sources and from Galactic sources, the correction for the CO emission, the effect of density inhomogeneities, and more. We present estimates of the local luminosity functions at 857 GHz (350 microns), 545 GHz (550 microns) and 353 GHz (850 microns) extending across the characteristic luminosity L_star, and a preliminary estimate over a limited luminosity range at 217 GHz (1382 microns). At 850 microns and for luminosities L >~ L_star our results agree with previous estimates, derived from the SCUBA Local Universe Galaxy Survey (SLUGS), but are higher than the latter at L <~ L_star. We also find good agreement with estimates at 350 and 500 microns based on preliminary Herschel survey data.

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

The infrared (IR) emission of M_* galaxies (10^{10.4} < M_{star} < 10^{11.0} M_\sun) in galaxy pairs, derived using data obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement induced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies (SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation enhancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z=0 and z=1, this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by a factor of ~10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic medium IGM in the dark matter halos (DMH) hosting these pairs.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

The populations of both quiescent and actively star-forming galaxies at 1<z<2 are still under-represented in our spectroscopic census of galaxies throughout the history of the Universe. In the light of galaxy formation models, however, the evolution of galaxies at these redshifts is of pivotal importance and merits further investigation. We therefore designed a spectroscopic observing campaign of a sample of both massive, quiescent and star-forming galaxies at z>1.4, called Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS). To determine redshifts and physical properties, such as metallicity, dust content, dynamical masses, and star formation history, we performed ultra-deep spectroscopy with the red-sensitive optical spectrograph FORS2 at the VLT. Our sample consists of objects, within the CDFS/GOODS area, detected at 4.5 micron, to be sensitive to stellar mass rather than star formation intensity. The spectroscopic targets were selected with a photometric redshift constraint (z>1.4) and magnitude constraints (B(AB)<26, I(AB)<26.5), which should ensure that these are faint, distant, and fairly massive galaxies. We present the sample selection, survey design, observations, data reduction, and spectroscopic redshifts. Up to 30 hours of spectroscopy of 174 spectroscopic targets and 70 additional objects enabled us to determine 210 redshifts, of which 145 are at z>1.4. From the redshifts and photometry, we deduce that the BzK selection criteria are efficient (82%) and suffer low contamination (11%). Several papers based on the GMASS survey show its value for studies of galaxy formation and evolution. We publicly release the redshifts and reduced spectra. In combination with existing and on-going additional observations in CDFS/GOODS, this data set provides a legacy for future studies of distant galaxies.

The geometry of the dust distribution within the inner regions of Active Galactic Nuclei (AGN) is still a debated issue and relates directly with the AGN unified scheme. Traditionally, models discussed in the literature assume one of two distinct dust distributions in what is believed to be a toroidal region around the Supermassive Black Holes: a continuous distribution, customarily referred to as smooth, and a concentration of dust in clumps or clouds, referred to as clumpy. In this paper we perform a thorough comparison between two of the most popular models in the literature, namely the smooth models by Fritz. et al. 2006 and the clumpy models by Nenkova et al. 2008a, in their common parameters space. Particular attention is paid to the silicate features at ~9.7 and ~18 micron, the width of the infrared bump, the near-infrared index and the luminosity at 12.3 micron, all previously reported as possible diagnostic tools to distinguish between the two dust distributions. We find that, due to the different dust chemical compositions used in the two models, the behaviour of the silicate features at 9.7 and 18 micron is quite distinct between the two models. The width of the infrared bump and the peak of the infrared emission can take comparable values, their distributions do, however, vary. The near-infrared index is also quite different, due partly to the primary sources adopted by the two models. Models with matched parameters do not produce similar SEDs and virtually no random parameter combinations can result in seemingly identical SEDs.

We discuss the clustering properties of galaxies with signs of ongoing star formation detected by the Spitzer Space Telescope at 24mum band in the SWIRE Lockman Hole field. The sample of mid-IR-selected galaxies includes ~20,000 objects detected above a flux threshold of S24mum=310muJy. We adopt optical/near-IR color selection criteria to split the sample into the lower-redshift and higher-redshift galaxy populations. We measure the angular correlation function on scales of theta=0.01-3.5 deg, from which, using the Limber inversion along with the redshift distribution established for similarly selected source populations in the GOODS fields (Rodighiero et al. 2010), we obtain comoving correlation lengths of r0=4.98+-0.28 h^-1 Mpc and r0 =8.04+-0.69 h^-1 Mpc for the low-z (<z>=0.7) and high-z (<z>=1.7) subsamples, respectively. Comparing these measurements with the correlation functions of dark matter halos identified in the Bolshoi cosmological simulation (Klypin et al. 2011}, we find that the high-redshift objects reside in progressively more massive halos reaching Mtot>3e12 h^-1 Msun, compared to Mtot>7e11 h^-1 Msun for the low-redshift population. Approximate estimates of the IR luminosities based on the catalogs of 24mum sources in the GOODS fields show that our high-z subsample represents a population of "distant ULIRGs" with LIR>10^12Lsun, while the low-z subsample mainly consists of "LIRGs", LIR~10^11Lsun. The comparison of number density of the 24mum selected galaxies and of dark matter halos with derived minimum mass Mtot shows that only 20% of such halos may host star-forming galaxies.

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $\mu$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Bo\"{o}tes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $\mu$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $\mu$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We use rest-frame UV spectroscopy to investigate the properties related to large-scale gas outflows, and to the dust extinction and star-formation rates of a sample of z ~ 2 star-forming galaxies from the Galaxy Mass Assembly ultradeep Spectroscopic Survey (GMASS). Dust extinction is estimated from the rest-frame UV continuum slope and used to obtain dust-corrected star-formation rates for the galaxies of the sample. For the entire sample, a mean value of the continuum slope <\beta> = -1.11 \pm 0.44 (r.m.s.) was derived, while the average SFR was found to be <SFR> = 52 \pm 48 M_sun/yr (r.m.s.). A positive correlation between SFR and stellar mass was observed, in agreement with other works, the logarithmic slope of the relation being 1.10 \pm 0.10. Low-ionization absorption lines, associated with the interstellar medium, were found to be blueshifted, with respect to the rest frame of the system, which indicates that there is outflowing gas with typical velocities of the order of ~ 100 km/s. Finally, investigating correlations between galaxy UV spectral characteristics and galaxy general properties, we report a possible correlation between the equivalent width of the interstellar absorption lines and SFR, stellar mass, and colour excess similar to that seen to hold at different redshifts.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

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.

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

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.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

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.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

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.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

[Abridged]We exploit a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent 24, 15, 12 and 8micron monochromatic and bolometric IR luminosity functions (LF) over the 0<z<2.5 redshift range. Our analysis is based on the combination of data from deep Spitzer surveys in the VVDS-SWIRE and GOODS areas. To our limiting flux of S(24)=400microJy our derived sample in VVDS-SWIRE includes 1494 sources, and 666 and 904 sources brighter than S(24)=80microJy are catalogued in GOODS-S and GOODS-N, respectively, for a total area of ~0.9 square degs. We obtain reliable optical identifications and redshifts, providing us a rich and robust dataset for our luminosity function determination. Based on the multi-wavelength information available, we constrain the LFs at 8, 12, 15 and 24micron. We also extrapolate total IR luminosities from our best-fit to the observed SEDs of each source, and use this to derive the bolometric LF and comoving volume emissivity up to z~2.5. In the 0<z<1 interval, the bolometric IR luminosity density evolves as (1+z)^3.8+/-0.4. Although more uncertain at higher-z, our results show a flattening of the IR luminosity density at z>1. The mean redshift of the peak in the source number density shifts with luminosity: the brighest IR galaxies appear to be forming stars earlier in cosmic time (z>1.5), while the less luminous ones keep doing it at more recent epochs (z~1 for L(IR)<10^11L_sun). Our results suggest a rapid increase of the galaxy IR comoving volume emissivity back to z~1 and a constant average emissivity at z>1. We also seem to find a difference in the evolution rate of the source number densities as a function of luminosity, a downsizing evolutionary pattern similar to that reported from other samples of cosmic sources.

[abridged] At the end of the Spitzer cryogenic mission and the onset of the Herschel era, we review our current knowledge on galaxy evolution at IR wavelengths. We also develop new tools for the analysis of background fluctuations to constrain source counts in regimes of high confusion. We analyse a large variety of new data on galaxy evolution and high-z source populations from Spitzer surveys, as well as complementary data from sub-mm (BLAST) and millimetric ground-based observations. These data confirm earlier indications about a very rapid increase of galaxy volume emissivity up to z~1. This is the fastest evolution rate observed for galaxies at any wavelengths. The observed Spitzer counts at 24 micron require a combination of fast evolution for the dominant population and a bumpy spectrum with substantial PAH emission at z~1 to 2. Confusion-limited number counts at longer wavelengths confirm these results. All the present data require that the fast observed evolution from z=0 to 1 flattens around redshift 1 and then keeps approximately constant up to z~2.5 at least. The present-day IR data provide clear evidence for the existence of a distinct population of very luminous galaxies becoming dominant at z > 1. Their cosmological evolution, peaking around z~2, shows a faster decay with cosmic time than lower luminosity systems, whose maximal activity is set around z~1, then supporting an earlier phase of formation for the most luminous and massive galaxies. From a comparison of our results on the comoving IR emissivity with recent estimates of the redshift-dependent stellar mass functions of galaxies, we infer that a large fraction (>=50%) of the IR activity at z > 1 should be due to obscured AGN accretion and that concomitant SF in high-z luminous sources should follow a top-heavy IMF.

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.

We present a detailed analysis of the far-IR properties of the bright, lensed, z = 2.3, SMG, SMM J2135-0102, using new observations with Herschel, SCUBA-2 and the VLA. These data allow us to constrain the galaxy's SED and show that it has an intrinsic rest-frame 8-1000um luminosity, L(bol), of (2.3 +/- 0.2) x 10^12 L(sun) and a likely SFR of ~400 M(sun)/yr. The galaxy sits on the far-IR/radio correlation for far-IR-selected galaxies. At ~>70um, the SED can be described adequately by dust components with T(d) ~ 30 and 60K. Using SPIRE's Fourier Transform Spectrometer we report a detection of the [CII] 158um cooling line. If the [CII], CO and far-IR continuum arise in photo-dissociation regions, we derive a characteristic gas density, n ~ 10^3 cm^-3, and a far-UV radiation field, G_0, 10^3x stronger than the Milky Way. L([CII])/L(bol) is significantly higher than in local ULIRGs but similar to the values found in local star-forming galaxies and starburst nuclei. This is consistent with SMM J2135-0102 being powered by starburst clumps distributed across ~2 kpc, evidence that SMGs are not simply scaled-up ULIRGs. Our results show that SPIRE's FTS has the ability to measure the redshifts of distant, obscured galaxies via the blind detection of atomic cooling lines, but it will not be competitive with ground-based CO-line searches. It will, however, allow detailed study of the integrated properties of high-redshift galaxies, as well as the chemistry of their ISM, once more suitably bright candidates have been found.

We present a new semi-analytical model of galaxy formation, GECO (Galaxy Evolution COde), aimed at a better understanding of when and how the two processes of star formation and galaxy assembly have taken place. Our model is structured into a Monte Carlo algorithm based on the Extended Press-Schechter theory, for the representation of the merging hierarchy of dark matter halos, and a set of analytic algorithms for the treatment of the baryonic physics, including classical recipes for the gas cooling, the star formation time-scales, galaxy mergers and SN feedback. Together with the galaxies, the parallel growth of BHs is followed in time and their feedback on the hosting galaxies is modelled. We set the model free parameters by matching with data on local stellar mass functions and the BH-bulge relation at z=0. Based on such local boundary conditions, we investigate how data on the high-redshift universe constrain our understanding of the physical processes driving the evolution, focusing in particular on the assembly of stellar mass and on the star formation history. Since both processes are currently strongly constrained by cosmological near- and far-IR surveys, the basic physics of the Lambda CDM hierarchical clustering concept of galaxy formation can be effectively tested by us by comparison with the most reliable set of observables. Our investigation shows that when the time-scales of the stellar formation and mass assembly are studied as a function of dark matter halo mass and the single galaxy stellar mass, the 'downsizing' fashion of star formation appears to be a natural outcome of the model, reproduced even in the absence of the AGN feedback. On the contrary, the stellar mass assembly history turns out to follow a more standard hierarchical pattern progressive in cosmic time, with the more massive systems assembled at late times mainly through dissipationless mergers.

The paper has been deeply reviewed and compeltely re-written.

The SAGACE experiment consists of a mm/sub-mm telescope with a 3-m diameter primary mirror, coupled to a cryogenic multi-beam differential spectrometer. SAGACE explores the sky in the 100-760 GHz frequency range, using four diffraction-limited bolometer arrays. The instrument is designed to perform spectroscopic surveys of the Sunyaev-Zeldovich effects of thousands of galaxy clusters, of the spectral energy distribution of active galactic nuclei, and of the [CII] line of a thousand galaxies in the redshift desert. In 2008 a full phase-A study for a national small mission was completed and delivered to the Italian Space Agency (ASI). We have shown that taking advantage of the differential operation of the Fourier Transform Spectrometer, this ambitious instrument can operate from a Molniya orbit, and can be built and operated within the tight budget of a small mission.

We report a moderate-depth (70 ksec), contiguous 0.7 sq.deg, Chandra survey, in the Lockman Hole Field of the Spitzer/SWIRE Legacy Survey coincident with a completed, ultra-deep VLA survey with deep optical and near-infrared imaging in-hand. The primary motivation is to distinguish starburst galaxies and AGN, including the significant, highly obscured (log N_H >23) subset. Chandra has detected 775 X-ray sources to a limiting broad band (0.3-8 keV) flux ~4E-16 erg/cm^2/s. We present the X-ray catalog, fluxes, hardness ratios and multi-wavelength fluxes. The log N vs. log S agrees with those of previous surveys covering similar flux ranges. The Chandra and Spitzer flux limits are well matched: 771 (99%) of the X-ray sources have IR or optical counterparts, and 333 have MIPS 24 um detections. There are 4 optical-only X-ray sources and 4 with no visible optical/IR counterpart. The very deep (~2.7 uJy rms) VLA data yields 251 (> 4 sigma) radio counterparts, 44% of the X-ray sources in the field. We confirm that the tendency for lower X-ray flux sources to be harder is primarily due to absorption. As expected, there is no correlation between observed IR and X-ray flux. Optically bright, Type 1 and red AGN lie in distinct regions of the IR vs X-ray flux plots, demonstrating the wide range of SEDs in this sample and providing the potential for classification/source selection. The classical definition of radio-loudness, R_L, is not effective at selecting strong radio sources for red and relatively optically faint AGN. Using the 24 um to radio flux ratio (q_24) results in the usual ~10% radio-loud AGN (13 of 147).

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.

We measured stellar velocity dispersions sigma and dynamical masses of 9 massive (M~10^11 Msun) early-type galaxies (ETG) from the GMASS sample at redshift 1.4<z<2.0. The sigma are based on individual spectra for two galaxies at z~1.4 and on a stacked spectrum for 7 galaxies with 1.6<z<2.0, with 202-h of exposure at the ESO Very Large Telescope. We constructed detailed axisymmetric dynamical models for the objects, based on the Jeans equations, taking the observed surface brightness (from deep HST/ACS observations), PSF and slit effects into account. Our dynamical masses M_Jeans agree within ~30% with virial estimates M_vir=5*Re*sigma^2/G, although the latter tend to be smaller. This suggests that sizes are not underestimated by more than a similar fraction. Our M_Jeans also agrees within a factor <2 with the M_pop previously derived using stellar population models and 11 bands photometry. This confirms that the galaxies are intrinsically massive. The inferred mass-to-light ratios M/L_U in the very age-sensitive rest frame U-band are consistent with passive evolution in the past ~1 Gyr (formation redshift z_f~3). A 'bottom-light' stellar Initial Mass Function (IMF) appears to be required to ensure close agreement between M_Jeans and M_pop at z~2, as it does at z~0. The GMASS ETGs are on average more dense than their local counterpart. However a few percent of local ETGs of similar dynamical masses also have comparable sigma and mass surface density Sigma_50 inside Re.

The BL Lac object 1ES 0229+200 (z=0.14) has been detected by HESS during observations taking place in 2005-2006. The TeV spectrum, when corrected for the absorption of gamma-ray photons through the interaction with the extragalactic background light, is extremely hard, even if the most conservative level for the background is considered. The case of 1ES 0229+200 is very similar to that of 1ES 1101-232, for which a possible explanation, in the framework of the standard one-zone synchrotron-self Compton model, is that the high-energy emission is synchrotron-self Compton radiation of electrons distributed as a power law with a large value of the minimum energy. In this scenario the hard TeV spectrum is accompanied by a very hard synchrotron continuum below the soft X-ray band. We will show that recent Swift observations of 1ES 0229+200 in the critical UV-X-ray band strongly support this model, showing the presence of the expected spectral break and hard continuum between the UV and the X-ray bands.

We present a study of a galaxy overdensity at z=1.6, Cl 0332-2742, uncovered in the Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS). We conservatively estimate the overdensity in redshift space for the spike, containing 42 galaxies at z=1.6, to be 8.3+/-1.5 and find a velocity dispersion of 500 km/s. A map of the surface density of galaxies at z=1.6 in the GMASS field shows that its structure is irregular. The differences in the physical properties of Cl 0332-2742 member and field galaxies agree with the latest hierarchical galaxy formation models: for overdensity members, the star formation rate (SFR), and specific SFR, is approximately 50% lower than for the field galaxies; overdensity galaxies are twice the age, on average, of field galaxies; and there is a higher proportion of both massive (M>10^{10.7} M_sol), and early-type galaxies, inside Cl 0332-2742 than in the field. Among the 42 members, 7 have spectra consistent with being passively evolving, massive galaxies. In a z-J colour-magnitude diagram, the photometric data of these early-type galaxies are in close agreement with a theoretical red sequence of a galaxy cluster at redshift z=1.6, which formed most of its stars in a short burst of star formation at z~3. We conclude that the redshift spike at z=1.6 in the GMASS field represents a sheet-like structure in the cosmic web and the area with the highest surface density within this structure, containing already seven passively evolving galaxies, will evolve into a cluster of galaxies at a later time.

We present submillimetre observations of the z=3.653 quasar SDSS160705+533558 together with data in the optical and infrared. The object is unusually bright in the far-IR and submm with an IR luminosity of ~10^14 L_sun. We ascribe this luminosity to a combination of AGN and starburst emission, with the starburst forming stars at a rate of a few thousand solar masses per year. Submillimetre Array (SMA) imaging observations with a resolution ~1" show that the submm (850 micron) emission is extended on scales of 10--35kpc and is offset from the optical position by ~10 kpc. This morphology is dissimilar to that found in submm galaxies, which are generally un- or marginally resolved on arcsecond scales, or submm-luminous AGN where the AGN lies at the peak of the submm or molecular emission. The simplest explanation is that the object is in the early stages of a merger between a gas rich galaxy, which hosts the starburst, and a gas-poor AGN-host galaxy, which is responsible for the quasar emission. It is also possible that jet induced star formation might contribute to the unusual morphology.

The properties of dust attenuation at rest-frame UV wavelengths are inferred from very high-quality FORS2 spectra of 78 galaxies from the GMASS survey at 1<z<2.5. These objects complement a previously investigated sample of 108 UV-luminous galaxies at similar redshifts, selected from the FDF spectroscopic survey, the K20 survey, and the GDDS. The shape of the UV extinction curve is constrained by a parametric description of the rest-frame UV continuum. The UV bump is further characterised by fitting Lorentzian-like profiles. Spectra exhibit a significant 2175A feature in at least 30% of the cases. If attenuation is dominated by dust ejected from the galaxy main body via superwinds, UV extinction curves in-between those of the SMC and LMC characterise the sample galaxies. The fraction of galaxies with extinction curves differing from the SMC one increases, if more dust resides in the galactic plane or dust attenuation depends on stellar age. On average, the width of the manifested UV bumps is about 60% of the values typical of the LMC and Milky Way. This suggests the presence of dust similar to that found in the LMC2 supershell close to 30Dor. The presence of the carriers of the UV bump at 1<z<2.5 argues for outflows from AGB stars being copious then. Consistent with their higher SFRs, the GMASS galaxies with a manifested UV bump are more luminous at rest-frame 8mum, where the emission is dominated by PAHs (also products of AGB stars). In addition, they exhibit stronger UV absorption features, mostly of interstellar origin, which indicates overall more evolved stellar populations. We conclude that diversification of the small-size dust component has already started in the most evolved star-forming systems at 1<z<2.5.

We describe the construction and the properties of the SWIRE-SDSS database, a preliminary derivation of the Far-Infrared Local Luminosity Functions at 24/70/160 micron based on such a database and ways in which VO tools will allow to refine and extend such work.

We report on--off pointed MAMBO observations at 1.2 mm of 61 Spitzer-selected star-forming galaxies from the SWIRE survey. The sources are selected on the basis of bright 24um fluxes (f_24um>0.4mJy) and of stellar dominated near-infrared spectral energy distributions in order to favor z~2 starburst galaxies. The average 1.2mm flux for the whole sample is 1.5+/-0.2 mJy. Our analysis focuses on 29 sources in the Lockman Hole field where the average 1.2mm flux (1.9+/-0.3 mJy) is higher than in other fields (1.1+/-0.2 mJy). The analysis of the sources multi-wavelength spectral energy distributions indicates that they are starburst galaxies with far-infrared luminosities ~10^12-10^13.3 Lsun, and stellar masses of ~0.2-6 x10^11 M_sun. Compared to sub-millimeter selected galaxies (SMGs), the SWIRE-MAMBO sources are among those with the largest 24um/millimeter flux ratios. The origin of such large ratios is investigated by comparing the average mid-infrared spectra and the stacked far-infrared spectral energy distributions of the SWIRE-MAMBO sources and of SMGs. The mid-infrared spectra exhibit strong PAH features, and a warm dust continuum. The warm dust continuum contributes to ~34% of the mid-infrared emission, and is likely associated with an AGN component. This constribution is consistent with what is found in SMGs. The large 24um/1.2mm flux ratios are thus not due to AGN emission, but rather to enhanced PAH emission compared to SMGs. The analysis of the stacked far-infrared fluxes yields warmer dust temperatures than typically observed in SMGs. Our selection favors warm ultra-luminous infrared sources at high-z, a class of objects that is rarely found in SMG samples. Our sample is the largest Spitzer-selected sample detected at millimeter wavelengths currently available.

The science case for FIR/Submm surveys of the extragalactic sky to be carried from Dome C are reviewed. The main questions concerning the formation and evolution of galaxies making up the CIRB are outlined and opportunities to exploit Dome C unique observing conditions through single-dish observations are discussed.

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.

(abridged) The ESIS survey is the optical follow up of the SWIRE/Spitzer in the ELAIS-S1 region of the sky. In the era of observational cosmology, the main efforts are focused on the study of galaxy evolution and its environmental dependence. Wide area, multiwavelength, extragalactic surveys are needed in order to probe sufficiently large volumes, minimize cosmic variance and find significant numbers of rare objects. We present VIMOS I and z band imaging belonging to the ESIS survey. A total of ~4 deg2 were targeted in I and ~1 deg2 in z. More than 300000 galaxies have been detected in the I band and ~50000 in the z band. Object coordinates are defined within an uncertainty of ~0.2 arcsec r.m.s., with respect to GSC 2.2. We reach a 90% average completeness at 23.1 and 22.5 mag (Vega) in the I and z bands, respectively. On the basis of IRAC colors, we identified galaxies having the 1.6 um stellar peak shifted to z=1-3. The new I, z band data provide reliable constraints to avoid low-redshift interlopers and reinforce this selection. Roughly 1000 galaxies between z=2-3 were identified over the ESIS ~4 deg2, at the SWIRE 5.8 um depth (25.8 uJy at 3sigma). These are the best galaxy candidates to dominate the massive tail (M>1e11 Msun) of the z>2 mass function.

The background radiations in the optical and the infrared constitute a relevant cause of energy loss in the propagation of high energy particles through space. In particular, TeV observations with Cherenkov telescopes of extragalactic sources are influenced by the opacity effects due to the interaction of the very high-energy source photons with the background light. With the aim of assessing with the best possible detail these opacity terms, we have modelled the extragalactic optical and IR backgrounds using available information on cosmic sources in the universe from far-UV to sub-mm wavelengths over a wide range of cosmic epochs. We have exploited the relevant cosmological survey data - including number counts, redshift distributions, luminosity functions - from ground-based observatories in the optical, near-IR, and sub-mm, as well as multi-wavelength information coming from space telescopes, HST, ISO and Spitzer. Additional constraints have been used from direct measurements or upper limits on the extragalactic backgrounds by dedicated missions (COBE). All data were fitted and interpolated with a multi-wavelength backward evolutionary model, allowing us to estimate the background photon density and its redshift evolution. From the redshift-dependent background spectrum, the photon-photon opacities for sources of high-energy emission at any redshifts were then computed. The same results can also be used to compute the optical depths for any kind of processes in the intergalactic space involving interactions with background photons (like scattering of cosmic-ray particles). We have applied our photon-photon opacity estimates to the analysis of spectral data at TeV energies on a few BLAZARs of particular interest. [abridged]

We derive the properties of dusty tori in Active Galactic Nuclei (AGN) from the comparison of observed Spectral Energy Distributions (SEDs) of SDSS quasars and a precomputed grid of torus models. The observed SEDs comprise SDSS photometry, 2MASS J, H, and K data, whenever available and mid-Infrared (MIR) data from the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. The adopted model is that of Fritz et al., 2006. The fit is performed by standard chi^2 minimisation, the model however can be multi-component comprising a stellar and a starburst components, whenever necessary. Models with low equatorial optical depth, tau_9.7, were allowed as well as ``traditional'' models with tau_9.7 > 1.0, corresponding to A_V > 22 and the results were compared. Fits using high optical depth tori models only produced dust more compactly distributed than in the configuration where all tau_9.7 models were permitted. Tori with decreasing dust density with the distance from the centre were favoured while there was no clear preference for models with or without angular variation of the dust density. The computed outer radii of the tori are of some tens of parsecs large but can reach, in a few cases, a few hundreds of parsecs. The mass of dust, M_Dust, and infrared luminosity, L_IR, integrated in the wavelength range between 1 and 1000 micron, do not show significant variations with redshift, once the observational biases are taken into account. Objects with 70 micron detections, representing 25% of the sample, are studied separately and the starburst contribution (whenever present) to the IR luminosity can reach, in the most extreme but very few cases, 80%.

The Galaxy Mass Assembly ultra-deep Spectroscopic Survey samples a part of the CDFS to unprecedented depth. The resulting distribution of 150 z>1.4 redshifts reveals a significant peak at z=1.6, part of a larger overdensity found at this redshift. The 42 spectroscopic members of this structure, called Cl 0332-2742, form an overdensity in redshift of a factor 11+/-3 and have a velocity dispersion of 450 km/s. We derive a total mass for Cl 0332-2742 of ~7x10^14 Msol. The colours of its early-type galaxies are consistent with a theoretical red sequence of galaxies with stars formed at z=3.0. In addition, there are more massive, passive and older, but less star forming galaxies in CL 0332-2742 than in the field. We conclude that this structure is a cluster under assembly at z=1.6.

The aim of this work is to study the evolution of the rest frame color distribution of galaxies with the redshift, in particular in the critical interval 1.4<z<3. We combine ultradeep spectroscopy from the GMASS project (Galaxy Mass Assembly ultradeep Spectroscopic Survey) with GOODS multi-band photometry (from optical to mid-infrared) to study a sample of 1021 galaxies up to m(4.5$\mu$m)=23. We find that the distribution of galaxies in the (U-B) color vs stellar mass plane is bimodal up to at least redshift z=2. We define a mass complete sample of galaxies residing on the red-sequence, selecting objects with log(M/M_{\odot})>10.1, and we study their morphological and spectro-photometric properties. We show that the contribution to this sample of early-type galaxies, defined as galaxies with a spheroidal morphology and no star formation, decreases from 60-70% at z<0.5 down to ~50% at redshift z=2. At z>2 we still find red galaxies in the mass complete sample, even if the bimodality is not seen any more. About 25% of these red galaxies at z>2 are passively evolving, with the bulk of their stars formed at redshift z>`3.

We combine ultradeep optical spectroscopy from the GMASS project ("Galaxy Mass Assembly ultradeep Spectroscopic Survey") with GOODS multi-band photometry and HST imaging to study a sample of passive galaxiesat 1.39<z<1.99 selected at 4.5 microns. A stacked spectrum with an equivalent integration time of ~500 hours was obtained is publicly released. The spectral and photometric SED properties indicate very weak or absent star formation, moderately old stellar ages of ~1 Gyr (for solar metallicity) and stellar masses in the range of 10^{10-11} solar masses, thus implying that the major star formation and assembly processes for these galaxies occurred at z>2. These galaxies have morphologies that are predominantly compact and spheroidal.However, their sizes (R_e <~ 1 kpc) are much smaller than those of spheroids in the present--day Universe. Their stellar mass surface densities are consequently higher by ~1 dex if compared to spheroids at z~0 with the same mass. Their rest-frame B-band surface brightness scales with the effective radius, but the offset with respect to the surface brightness of the local Kormendy relation is too large to be explained by simple passive evolution. At z~1, a larger fraction of passive galaxies follows the z~0 size -- mass relation. Superdense relics with R_e~1 kpc are extremely rare at z~0 with respect to z>1, and absent if R_e<1 kpc. Because of the similar sizes and mass densities, we suggest that the superdense passive galaxies at 1<z<2 are the remnants of the powerful starbursts occurring in submillimeter--selected galaxies at z>2. The results are compared with theoretical models and the main implications discussed in the framework of massive galaxy formation and evolution.

Context: Galaxy metallicities have been measured to redshift z~2 by gas-phase oxygen abundances of the interstellar medium using the R23 and N2 methods. Galaxy stellar metallicities provide crucial data for chemical evolution models but have not been assessed reliably much outside the local Universe. Aims: We determine the iron-abundance, stellar metallicity of star-forming galaxies (SFGs) at redshift z~2, observed as part of the Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS). Methods: We compute the equivalent width of a rest-frame mid-ultraviolet, photospheric absorption-line index, the 1978 index found to vary monotonically with stellar metallicity by Rix and collaborators. We normalise and combine 75 SFG spectra from the GMASS survey to produce a spectrum corresponding to a total integration time 1652.5 hours (and a signal-to-noise ratio ~100 for our 1.5 angstrom binning) of FORS2 spectroscopic observations at the Very Large Telescope. Results: We measure an iron-abundance, stellar metallicity of log (Z/Zsolar) = -0.574+/-0.159 for our spectrum representative of a galaxy of stellar mass 9.4 x 10^9 Msolar assuming a Chabrier IMF. We find that the R04 model SFG spectrum for log (Z/Zsolar) = -0.699 solar metallicity provides the best description of our GMASS coadded spectrum. For similar galaxy stellar mass, our stellar metallicity is ~0.25 dex lower than the oxygen-abundance, gas-phase metallicity quantified by Erb and collaborators for UV-selected star-forming galaxies at z=2. Conclusions: We conclude that we are witnessing the establishment of a light-element overabundance in galaxies as they are being formed at redshift z~2. Our measurements are reminiscent of the alpha-element enhancement seen in low-redshift, galactic bulges and early-type galaxies. (Abridged)

This paper presents galaxy source counts at 24 microns in the six Spitzer Wide-field InfraRed Extragalactic (SWIRE) fields. The source counts are compared to counts in other fields, and to model predictions that have been updated since the launch of Spitzer. This analysis confirms a very steep rise in the Euclidean-normalized differential number counts between 2 mJy and 0.3 mJy. Variations in the counts between fields show the effects of sample variance in the flux range 0.5-10 mJy, up to 100% larger than Poisson errors. Nonetheless, a "shoulder" in the normalized counts persists at around 3 mJy. The peak of the normalized counts at 0.3 mJy is higher and narrower than most models predict. In the ELAIS N1 field, the 24 micron data are combined with Spitzer-IRAC data and five-band optical imaging, and these bandmerged data are fit with photometric redshift templates. Above 1 mJy the counts are dominated by galaxies at z less than 0.3. By 300 microJy, about 25% are between z ~ 0.3-0.8, and a significant fraction are at z ~ 1.3-2. At low redshifts the counts are dominated by spirals, and starbursts rise in number density to outnumber the spirals' contribution to the counts below 1 mJy.

We investigate the clustering of galaxies selected in the 3.6 micron band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for eleven samples with flux limits of S_3.6 > 4-400 mujy, over an 8 square degree field. The angular clustering strength is measured at >5-sigma significance at all flux limits, with amplitudes of A=(0.49-29)\times10^{-3} at one degree, for a power-law model, A\theta^{-0.8}. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the GalICS (Galaxies In Cosmological Simulations) models of galaxy evolution and with parameterized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z<0.5) the comoving correlation length is approximately constant, r_0=6.1\pm0.5h^{-1} Mpc, and then decreases with increasing redshift to a value of 2.9\pm0.3h^{-1} Mpc for the faintest sample, for which the median redshift is z=1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over 0.5<z<1.

We use a simple optical/infrared (IR) photometric selection of high-redshift QSOs that identifies a Lyman Break in the optical photometry and requires a red IR color to distinguish QSOs from common interlopers. The search yields 100 z~3 (U-dropout) QSO candidates with 19<r'<22 over 11.7 deg^2 in the ELAIS-N1 (EN1) and ELAIS-N2 (EN2) fields of the Spitzer Wide-area Infrared Extragalactic (SWIRE) Legacy Survey. The z~3 selection is reliable, with spectroscopic follow-up of 10 candidates confirming they are all QSOs at 2.83<z<3.44. We find that our z~4$ (g'-dropout) sample suffers from both unreliability and incompleteness but present 7 previously unidentified QSOs at 3.50<z<3.89. Detailed simulations show our z~3 completeness to be ~80-90% from 3.0<z<3.5, significantly better than the ~30-80% completeness of the SDSS at these redshifts. The resulting luminosity function extends two magnitudes fainter than SDSS and has a faint end slope of beta=-1.42 +- 0.15, consistent with values measured at lower redshift. Therefore, we see no evidence for evolution of the faint end slope of the QSO luminosity function. Including the SDSS QSO sample, we have now directly measured the space density of QSOs responsible for ~70% of the QSO UV luminosity density at z~3. We derive a maximum rate of HI photoionization from QSOs at z~3.2, Gamma = 4.8x10^-13 s^-1, about half of the total rate inferred through studies of the Ly-alpha forest. Therefore, star-forming galaxies and QSOs must contribute comparably to the photoionization of HI in the intergalactic medium at z~3.

(Abridged) We selected high-z massive galaxies at 5.8 microns, in the SWIRE ELAIS-S1 field (1 sq. deg.). Galaxies with the 1.6 microns stellar peak redshifted into the IRAC bands (z~1-3, called ``IR-peakers'') were identified. Stellar masses were derived by means of spectro-photometric fitting and used to compute the stellar mass function (MF) at z=1-2 and 2-3. A parametric fit to the MF was performed, based on a Bayesian formalism, and the stellar mass density of massive galaxies above z=2 determined. We present the first systematic study of the very-massive tail of the galaxy stellar mass function at high redshift. A total of 326 sources were selected. The majority of these galaxies have stellar masses in excess of 1e11 Msun and lie at z>1.5. The availability of mid-IR data turned out to be a valuable tool to constrain the contribution of young stars to galaxy SEDs, and thus their M(stars)/L ratio. The influence of near-IR data and of the chosen stellar library on the SED fitting are also discussed. A significant evolution is found not only for galaxies with M~1e11 Msun, but also in the highest mass bins considered. The comoving number density of these galaxies was lower by more than a factor of 10 at z=2-3, with respect to the local estimate. SWIRE 5.8 micron peakers more massive than 1.6x1e11 Msun provide 30-50% of the total stellar mass density in galaxies at z=2-3.

We present the results of infrared L-band (3-4 micron) and M-band (4-5 micron) VLT-ISAAC spectroscopy of five bright Ultraluminous InfraRed Galaxies (ULIRGs) hosting an AGN. From our analysis we distinguish two types of sources: ULIRGs where the AGN is unobscured (with a flat continuum and no absorption features at 3.4 micron and 4.6 micron), and those with highly obscured AGNs (with a steep, reddened continuum and absorption features due to hydrocarbons and CO). Starburst activity is also present in all of the sources as inferred from the 3.3 micron PAH emission line. A strong correlation is found between continuum slope and CO optical depth, which suggests that deep carbon monoxide absorption is a common feature of highly obscured ULIRG AGN. Finally we show that the AGN dominates the 3-4 micron emission, even if its contribution to the bolometric luminosity is small.

(Abridged) This work explores the existence of high redshift massive galaxies unveiled with Spitzer+IRAC, but missed by conventional selection techniques based on optical and near-infrared observations. To this end, we use the multi-wavelength imaging data available for the GOODS-South field, and select a flux-limited sample from the IRAC 3.6um image to m(AB)<23.26. We confine our study to the galaxies undetected by the optical HST+ACS imaging and close to the detection limit of the K-band image (K>23.5 AB). Our selection unveiled 20 galaxies on which we performed a detailed photometric analysis. For each galaxy, we built an SED based on optical-to-8um photometry to estimate the photo-z and to derive the main galaxies physical properties. The majority of the sample sources show degenerate/bimodal solutions for the photometric redshifts (Abridged). These can either be heavily dust-enshrouded (Av~2-4) starbursts at 2<z<3 with bolometric luminosities L(IR)>10^12 Lsun, or massive post-starburst galaxies in the redshift interval 4<z<9 with stellar masses of 10^11 Msun. One galaxy, the only source in our sample with both an X-ray and a 24um detection, might be an extremely massive object at z~8 detected during a post-starburst phase with concomitant QSO activity (although a lower-z solution is not excluded). Our investigation of Spitzer-selected galaxies with very red SEDs and completely undetected in the optical reveals a potential population of massive galaxies at z>4 which appear to include significant AGN emissions. These sources may be the oldest stellar systems at z~4. These, previously unrecognized, optically obscured objects might provide an important contribution to the massive-end (M>10^11 sun) of the high-z stellar mass function and they would almost double it (Abridged).

Our goal is to probe the populations of obscured and unobscured AGN investigating their optical-IR and X-ray properties as a function of X-ray flux, luminosity and redshift within a hard X-ray selected sample of 136 X-ray sources in the XMM Medium Deep Survey (XMDS) with wide multiwavelength coverage. The XMDS area is covered with optical photometry from the VVDS and CFHTLS surveys and infrared Spitzer data. Based on the X-ray luminosity and X-ray to optical ratio, 132 sources are likely AGN, of which 122 have unambiguous optical - IR identification. The observed optical and IR spectral energy distributions of sources are fitted with AGN/galaxy templates in order to classify them and compute photometric redshifts. 70% of the AGN are fitted by a type 2 AGN or a star forming galaxy template and are grouped together in a single class of ``optically obscured'' AGN. They have ``red'' optical colors and generally show significant X-ray absorption from X-ray spectra or hardness ratios (N$_H > 10^{22}$ cm$^{-2}$). Sources with SEDs typical of type 1 AGN have ``blue'' optical colors and exhibit X-ray absorption in about 30% of cases. We performed a stacking analysis for obscured and type 1 AGN. The stacked X-ray spectrum of obscured AGN is flatter than that of type 1 AGN and has an average spectral slope of Gamma = 1.6. The subsample of objects fitted by a galaxy template has an even harder stacked spectrum, with Gamma = 1.2 - 1.3. The obscured fraction is larger at lower fluxes, lower redshifts and lower luminosities. X-ray absorption is less common than ``optical'' obscuration and its incidence is nearly constant with redshift and luminosity. This implies that X-ray absorption is not necessarily related to optical obscuration.

We present the spectral energy distributions (SEDs) of a hard X-ray selected sample. The sample contains 136 sources with F(2-10 keV)>10^-14 erg/cm^2/s and 132 are AGNs. The sources are detected in a 1 square degree area of the XMM-Newton-Medium Deep Survey where optical data from the VVDS, CFHTLS surveys, and infrared data from the SWIRE survey are available. Based on a SED fitting technique we derive photometric redshifts with sigma(1+z)=0.11 and 6% of outliers and identify AGN signatures in 83% of the objects. This fraction is higher than derived when a spectroscopic classification is available. The remaining 17+9-6% of AGNs shows star-forming galaxy SEDs (SF class). The sources with AGN signatures are divided in two classes, AGN1 (33+6-1%) and AGN2 (50+6-11). The AGN1 and AGN2 classes include sources whose SEDs are fitted by type 1 and type 2 AGN templates, respectively. On average, AGN1s show soft X-ray spectra, consistent with being unabsorbed, while AGN2s and SFs show hard X-ray spectra, consistent with being absorbed. The analysis of the average SEDs as a function of X-ray luminosity shows a reddening of the IR SEDs, consistent with a decreasing contribution from the host galaxy at higher luminosities. The AGNs in the SF classes are likely obscured in the mid-infrared, as suggested by their low L(3-20micron)/Lcorr(0.5-10 keV) ratios. We confirm the previously found correlation for AGNs between the radio luminosity and the X-ray and the mid-infrared luminosities. The X-ray-radio correlation can be used to identify heavily absorbed AGNs. However, the estimated radio fluxes for the missing AGN population responsible for the bulk of the background at E>10 keV are too faint to be detected even in the deepest current radio surveys.

(Abridged) High-redshift ultra luminous infrared galaxies contribute the bulk of the cosmic IR background and are the best candidates for very massive galaxies in formation at z>1.5. We present Keck/LRIS optical spectroscopy of 35 z>1.4 luminous IR galaxies in the Spitzer Wide-area Infra-Red Extragalactic survey (SWIRE) northern fields (Lockman Hole, ELAIS-N1, ELAIS-N2). The primary targets belong to the ``IR-peak'' class of galaxies, having the 1.6 micron (restframe) stellar feature detected in the IRAC Spitzer channels.The spectral energy distributions of the main targets are thoroughly analyzed, by means of spectro-photometric synthesis and multi-component fits (stars + starburst dust + AGN torus). The IR-peak selection technique is confirmed to successfully select objects above z=1.4, though some of the observed sources lie at lower redshift than expected. Among the 16 galaxies with spectroscopic redshift, 62% host an AGN component, two thirds being type-1 and one third type-2 objects. The selection, limited to r'<24.5, is likely biased to optically-bright AGNs. The SEDs of non-AGN IR-peakers resemble those of starbursts (SFR=20-500 Msun/yr) hosted in massive (M>1e11 Msun) galaxies. The presence of an AGN component provides a plausible explanation for the spectroscopic/photometric redshift discrepancies, as the torus produces an apparent shift of the peak to longer wavelengths. These sources are analyzed in IRAC and optical-IR color spaces. In addition to the IR-peak galaxies, we present redshifts and spectral properties for 150 objects, out of a total of 301 sources on slits.

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

We study the enviromental dependence and the morphological composition of the galaxy color-magnitude diagram at z~0.7, using a pilot sub-sample of ~2000 galaxies from the COSMOS surve, with I_AB<24, photometric redshift within 0.61<z<0.85 and galaxy morphologies based on the HST-ACS data. The C-M diagram shows a clear red-sequence dominated by early-type galaxies and also a remarkably well-defined "blue sequence'" described by late-type objects. While the percentage of objects populating the two sequences is a function of environment, with a clear morphology/color-density relation, their normalization and slope are independent of local density. We identify and study a number of objects with "anomalous'' colors, given their morphology. Red late-type galaxies are found to be mostly highly-inclined spiral galaxies, with colors dominated by internal reddening (33% contamination with respect to truly passive spheroidals). Conversely, the population of blue early-type galaxies is composed by objects of moderate luminosity and mass, concurring to only ~5% of the mass in spheroidal galaxies. The majority of them (~70%) occupy a position in the surface-brightness/effective-radius plane not consistent with them being precursors of current epoch elliptical galaxies. In a color-mass diagram, color sequences are even better defined, with red galaxies covering in general a wider range of masses at nearly constant color, and blue galaxies showing a more pronounced dependence of color on mass. While the red sequence is adequately reproduced by models of passive evolution, the blue sequence is better interpreted as a specific star-formation sequence. Its substantial invariance with respect to local density suggests that the overall, "secular'' star formation is driven more by galaxy mass than by environment. [...]

We have identified a large-scale structure at z~0.73 in the COSMOS field, coherently described by the distribution of galaxy photometric redshifts, an ACS weak-lensing convergence map and the distribution of extended X-ray sources in a mosaic of XMM observations. The main peak seen in these maps corresponds to a rich cluster with Tx= 3.51+0.60/-0.46 keV and Lx=(1.56+/-0.04) x 10^{44} erg/s ([0.1-2.4] keV band). We estimate an X-ray mass within $r500$ corresponding to M500~1.6 x 10^{14} Msun and a total lensing mass (extrapolated by fitting a NFW profile) M(NFW)=(6+/-3) x 10^15 Msun. We use an automated morphological classification of all galaxies brighter than I_AB=24 over the structure area to measure the fraction of early-type objects as a function of local projected density Sigma_10, based on photometric redshifts derived from ground-based deep multi-band photometry. We recover a robust morphology-density relation at this redshift, indicating, for comparable local densities, a smaller fraction of early-type galaxies than today. Interestingly, this difference is less strong at the highest densities and becomes more severe in intermediate environments. We also find, however, local "inversions'' of the observed global relation, possibly driven by the large-scale environment. In particular, we find direct correspondence of a large concentration of disk galaxies to (the colder side of) a possible shock region detected in the X-ray temperature map and surface brightness distribution of the dominant cluster. We interpret this as potential evidence of shock-induced star formation in existing galaxy disks, during the ongoing merger between two sub-clusters.

The Spitzer-COSMOS survey (S-COSMOS) is a Legacy program (Cycles 2+3) designed to carry out a uniform deep survey of the full 2 sq deg COSMOS field in all seven Spitzer bands (3.6, 4.5, 5.6, 8.0, 24.0, 70.0, 160.0 u). This paper describes the survey parameters, mapping strategy, data reduction procedures, achieved sensitivities to date, and the complete data set for future reference. We show that the observed infrared backgrounds in the S-COSMOS field are within 10% of the predicted background levels. The fluctuations in the background at 24u have been measured and do not show any significant contribution from cirrus, as expected. In addition, we report on the number of asteroid detections in the low galactic latitude COSMOS field. We use the Cycle 2 S-COSMOS data to determine preliminary number counts, and compare our results with those from previous Spitzer Legacy surveys (e.g. SWIRE, GOODS). The results from this "first analysis" confirm that the S-COSMOS survey will have sufficient sensitivity with IRAC to detect ~ L* disks and spheroids out to z ~ 3, and with MIPS to detect ultraluminous starbursts and AGN out to z ~3 at 24u and out to z ~1.5-2 at 70u and 160u.

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 present the first set of XMM-Newton EPIC observations in the 2 square degree COSMOS field. The strength of the COSMOS project is the unprecedented combination of a large solid angle and sensitivity over the whole multiwavelength spectrum. The XMM-Newton observations are very efficient in localizing and identifying active galactic nuclei (AGN) and clusters as well as groups of galaxies. One of the primary goals of the XMM-Newton Cosmos survey is to study the co-evolution of active galactic nuclei as a function of their environment in the Cosmic web. Here we present the log of observations, images and a summary of first research highlights for the first pass of 25 XMM-Newton pointings across the field. In the existing dataset we have detected 1416 new X-ray sources in the 0.5-2, 2-4.5 and 4.5-10 keV bands to an equivalent 0.5-2 keV flux limit of 7x10-16 erg cm-2 s-1. The number of sources is expected to grow to almost 2000 in the final coverage of the survey. From an X-ray color color analysis we identify a population of heavily obscured, partially leaky or reflecting absorbers, most of which are likely to be nearby, Compton-thick AGN.

The Spitzer Space Telescope is devoting a significant fraction of the observing time to multi-wavelength cosmological surveys of different depths in various low-background sky regions. Several tens of thousand mid-IR galaxies have been detected over a wide interval of redshifts. A progressively clearer picture of galaxy evolution is emerging, which emphasizes populations of luminous galaxies at z>1,likely corresponding to the main phases of stellar formation and galaxy assembly. These results are entirely consistent with previous outcomes from ISO, SCUBA and COBE observations, and provide valuable constraints of high statistical and photometric precision. We briefly report here on our attempt to extract from statistical data some general properties of galaxy evolution and describe evidence that a population of very luminous objects at z>1.5 share different properties from those of starbursts at lower redshifts, indicating some seemingly anti-hierarchical behavior of galaxy evolution in the IR. We warn, however, that these results are based on large, still uncertain, extrapolations of the observed mid-IR to bolometric fluxes, for measuring which the forthcoming far-IR and submillimetre Herschel Space Observatory will be needed. We finally comment, based on our present understanding, about Herschel capabilities for investigating the early phases of galaxy evolution.

SPIRE, the Spectral and Photometric Imaging Receiver, is one of three instruments to fly on ESA's Herschel Space Observatory. It contains a three-band imaging photometer operating at 250, 360 and 520 microns, and an imaging Fourier Transform Spectrometer (FTS) covering 200-670 microns. It will be used for many extragalactic science programmes, a number of which will be implemented as Herschel Key Projects. The SPIRE consortium's Guaranteed Time (GT) programme will devote more than 1000 hours to Key Projects covering the high-z universe and local galaxies. It is also expected that substantial amounts of Herschel Open Time will be used for further extragalactic investigations. The high-z part of the SPIRE GT programme will focus on blank-field surveys with a range of depths and areas optimised to sample the luminosity-redshift plane and characterise the bolometric luminosity density of the universe at high z. Fields will be selected that are well covered at other wavelengths to facilitate source identifications and enable detailed studies of the redshifts, spectral energy distributions, and infrared properties of detected galaxies. The local galaxies programme will include a detailed spectral and photometric study of a sample of well resolved nearby galaxies, a survey of more than 300 local galaxies designed to provide a statistical survey of dust in the nearby universe, and a study of the ISM in low-metallicity environments, bridging the gap between the local universe and primordial galaxies.

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.

We construct rest-frame luminosity functions at 3.6, 4.5, 5.8, 8 and 24 microns over the redshift range 0<z<2 for galaxies and 0<z<4 for optical QSOs, using optical and infrared data from the Spitzer Wide-area InfraRed Extragalactic survey. The 3.6 and 4.5 micron galaxy LFs show evidence for moderate positive luminosity evolution up to z~1.5, consistent with the passive ageing of evolved stellar populations. Their comoving luminosity density was found to evolve passively, gradually increasing out to z~0.5-1 but flattening, or even declining, at higher redshift. Conversely, the 24 micron galaxy LF, which is more sensitive to obscured star formation and/or AGN activity, undergoes strong positive evolution, with the derived IR energy density and SFR density proportional to (1+z)^gamma with gamma=4.5^{+0.7}_{-0.6} and the majority of this evolution occurring since z~1. Optical QSOs, however, show positive luminosity evolution in all bands, out to the highest redshifts (3<z<4). Modelling as L* proportional to (1+z)^gamma gave gamma=1.3^{+0.1}_{-0.1} at 3.6 micron, gamma=1.0^{+0.1}_{-0.1} at 4.5 micron and stronger evolution at the longer wavelengths (5.8, 8 and 24 micron), of gamma~3. Comparison of the galaxy LFs to predictions from a semi-analytic model based on CDM indicate that an IMF skewed towards higher mass star formation in bursts compared to locally is preferred. As a result the currently inferred massive star formation rates in distant sub-mm sources may require substantial downwards revision.

One of the main questions of galaxy formation is how and when galaxies assembled their mass. Observations suggest that most mass assembly and star formation occured at 1 < z < 3. However, in the same redshift range, fully assembled old and massive galaxies co-exist with star-bursting systems still in their assembling stage. In order to explore galaxy evolution in this critical redshift range, we started a project called GMASS (Galaxy Mass Assembly ultra-deep Spectroscopic Survey) based on an 145 hours ESO Large Program and aimed at doing ultradeep spectroscopy with the VLT of a sample of high redshift galaxies slected with IRAC/Spitzer m(4.5mu) < 23, z(phot) > 1.4) from the GOODS-South/HUDF region. The selection at 4.5 micron is crucial to select high redshift galaxies based on the stellar mass or objects affected by strong dust extinction. The spectroscopic deep integration times (15 - 50 hours) are needed to identify unambiguously the redshift and nature of these galaxies, which redshifts are mostly in the "desert" where spectroscopy is traditionally difficult due to the absence of strong emission lines at optical wavelengths. The survey is almost complete and some preliminary results are presented here.

We discuss spectral energy distributions, photometric redshifts, redshift distributions, luminosity functions, source-counts and the far infrared to optical luminosity ratio for sources in the SWIRE Legacy Survey. The spectral energy distributions of selected SWIRE sources are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 5%. While further spectroscopic data are needed to confirm this result, the prospect of determining good photometric redshifts for the 2 million extragalactic objects in SWIRE is excellent. The distribution of the different infrared sed types in the L{ir}/L{opt} versus L{ir} plane, where L{ir} and L{opt} are the infrared and optical bolometric luminosities, is discussed. Source-counts at 24, 70 and 160 mu are discussed, and luminosity functions at 3.6 and 24 mu are presented.

This work focuses on the properties of dusty tori in active galactic nuclei (AGN) derived from the comparison of SDSS type 1 quasars with mid-Infrared (MIR) counterparts and a new, detailed torus model. The infrared data were taken by the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. Basic model parameters are constraint, such as the density law of the graphite and silicate grains, the torus size and its opening angle. A whole variety of optical depths is supported. The favoured models are those with decreasing density with distance from the centre, while there is no clear tendency as to the covering factor, ie small, medium and large covering factors are almost equally distributed. Based on the models that better describe the observed SEDs, properties such as the accretion luminosity, the mass of dust, the inner to outer radius ratio and the hydrogen column density are computed.

Using the large multi-wavelength data set in the chandra/SWIRE Survey (0.6 square degrees in the Lockman Hole), we show evidence for the existence of highly obscured (Compton-thick) AGN, estimate a lower limit to their surface density and characterize their multi-wavelength properties. Two independent selection methods based on the X-ray and infrared spectral properties are presented. The two selected samples contain 1) 5 X-ray sources with hard X-ray spectra and column densities > 10^24 cm-2, and 2) 120 infrared sources with red and AGN-dominated infrared spectral energy distributions (SEDs). We estimate a surface density of at least 25 Compton-thick AGN per square degree detected in the infrared in the chandra/SWIRE field of which ~40% show distinct AGN signatures in their optical/near-infrared SEDs, the remainings being dominated by the host-galaxy emission. Only ~33% of all Compton-thick AGN are detected in the X-rays at our depth (F(0.3-8 keV)>10^-15 erg/cm2/s. We report the discovery of two sources in our sample of Compton-thick AGN, SWIRE_J104409.95+585224.8 (z=2.54) and SWIRE_J104406.30+583954.1 (z=2.43), which are the most luminous Compton-thick AGN at high-z currently known. The properties of these two sources are discussed in detail with an analysis of their spectra, SEDs, luminosities and black-hole masses.

We compare the history of the galaxy mass build-up, as inferred from near-IR observations, and the Star Formation Rate of massive stars in the comoving volume traced by deep extensive far-IR surveys, both possible now with the Spitzer Space Telescope. These two independent and complementary approaches to the history of galaxy formation consistently indicate that a wide interval of cosmic epochs between z~0.7 to z~2 brackets the main evolutionary phases. The rate of the integrated galaxy mass growth indicated by the IR-based comoving SFR appears consistent with the observed decrease of the stellar mass densities with redshift. There are also indications that the evolution with z of the total population depends on galaxy mass, being stronger for moderate-mass, but almost absent up to z=1.4 for high-mass galaxies, thus confirming previous evidence for a "downsizing" effect in galaxy formation. The most massive galaxies appear already mostly in place by z~1. Although a precise matching of this galaxy build-up with the growth of nuclear super-massive black-holes is not possible with the present data (due to difficulties for an accurate census of the obscured AGN phenomenon), some preliminary indications reveal a similar mass/luminosity dependence for AGN evolution as for the hosting galaxies.

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 constrain the evolution of the galaxy mass and luminosity functions from the analysis of (public) multi-wavelength data in the Chandra Deep Field South (CDFS) area, obtained from the GOODS and other projects, and including very deep high-resolution imaging by HST/ACS. Our reference catalogue of faint high-redshift galaxies, which we have thoroughly tested for completeness and reliability, comes from a deep (S(3.6micron)>1 microJy) image by IRAC on the Spitzer Observatory. These imaging data in the field are complemented with extensive optical spectroscopy by the ESO VLT/FORS2 and VIMOS spectrographs, while deep K-band VLT/ISAAC imaging is also used to derive further complementary statistical constraints and to assist the source identification and SED analysis. We have selected a highly reliable IRAC 3.6micron sub-sample of 1478 galaxies with S(3.6)>10microJy, 47% of which have spectroscopic redshift, while for the remaining objects both COMBO-17 and Hyperz are used to estimate the photometric redshift. This very extensive dataset is exploited to assess evolutionary effects in the galaxy luminosity and stellar mass functions, while luminosity/density evolution is further constrained with the number counts and redshift distributions. The deep ACS imaging allows us to differentiate these evolutionary paths by morphological type, which our simulations show to be reliable at least up to z=1.5 for the two main early- (E/S0) and late-type (Sp/Irr) classes. These data, as well as our direct estimate of the stellar mass function above M=10^(10)M_sun for the spheroidal subclass, consistently evidence a progressive dearth of such objects to occur starting at z=0.7, paralleled by an increase in luminosity. (abridged)

We present 3-5 micron spectroscopy of the interacting system NGC 6240, showing the presence of two active galactic nuclei. The brightest (southern) nucleus shows up with a starburst-like emission, with a prominent 3.3 micron emission feature. However, the presence of an AGN is revealed by the detection of a broad Br alpha emission line, with a width of ~1,800 km/s. The spectrum of the faintest (northern) nucleus shows typical AGN features, such as a steep continuum and broad absorption features in the M-band. We discuss the physical properties of the dusty absorbers/emitters, and show that in both nuclei the AGN is dominant in the 3-5 micron band, but its contribution to the total luminosity is small (a few percent of the starburst emission).

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 describe improved modelling of the emission by dust in a toroidal--like structure heated by a central illuminating source within Active Galactic Nuclei (AGN). We chose a simple but realistic torus geometry, a flared disc, and a dust grain distribution function including a full range of grain sizes. The optical depth within the torus is computed in detail taking into account the different sublimation temperatures of the silicate and graphite grains, which solves previously reported inconsistencies in the silicate emission feature in type-1 AGN. We exploit this model to study the spectral energy distributions (SEDs) of 58 extragalactic (both type-1 and type-2) sources using archival optical and infrared (IR) data. We find that both AGN and starburst contributions are often required to reproduce the observed SEDs, although in a few cases they are very well fitted by a pure AGN component. The AGN contribution to the far-IR luminosity is found to be higher in type-1 sources, with all the type-2 requiring a substantial contribution from a circum-nuclear starburst. Our results appear in agreement with the AGN Unified Scheme, since the distributions of key parameters of the torus models turn out to be compatible for type-1 and type-2 AGN. Further support to the unification concept comes from comparison with medium-resolution IR spectra of type-1 AGN by the Spitzer observatory, showing evidence for a moderate silicate emission around 10 \mums which our code reproduces. From our analysis we infer accretion flows in the inner nucleus of local AGN characterized by high equatorial optical depths ($A_V\simeq 100$), moderate sizes ($R_{max}<100 pc$) and very high covering factors ($f\simeq 80$ per cent) on average.

We present ground-based optical and Spitzer infrared imaging observations of the interacting galaxy UGC 10214, the "Tadpole Galaxy" (z = 0.0310), focusing on the star formation activity in the nuclear, disk, spiral arms and tidal tail regions. The major findings of this study are that the Tadpole is actively forming stars in the main disk outside of the nucleus and in the tidal plume, with an estimated mean star formation rate of ~2 to 4 M_sun/yr. The most prominent sites of mid-infrared emission define a "ring" morphology that, combined with the overall morphology of the system, suggest the interaction may belong to the rare class of off-center collisional ring systems that form both shock-induced rings of star formation and tidal plumes. The nuclear emission is solely powered by older stars, with little evidence for ongoing star formation at the center of the Tadpole. Extra-nuclear star formation accounts for >50% of the total star formation in the disk and spiral arms, featuring infrared-bright 'hot spots' that exhibit strong PAH emission, whose band strength is comparable to that of late-type star-forming disk galaxies. The tidal tail, which extends 2 arcmin (~75 kpc) into the intergalactic medium, is populated by super massive star clusters likely triggered by the galaxy-galaxy interaction that has distorted UGC 10214 into its current "tadpole" shape.

We present the results of L-band spectroscopical observations of local bright ULIRGs, performed with ISAAC at the VLT. We discuss the L-band spectral features of seven single sources, and the statistical properties of a complete sample of 15 sources obtained combining our observations with other published 3-4 micron spectra. Our main results are: 1. When a spectral indicator suggesting the presence of an AGN (low equivalent width of the 3.3 micron emission line, steep lambda-f_lambda spectrum, presence of an absorption feature at 3.4 micron) is found, the AGN is always confirmed by independent analysis at other wavelengths. Conversely, in all known AGNs at least one of the above indicators is present. 2. Two new diagnostic diagrams are proposed combining the above indicators, in which starbursts and AGNs are clearly and completely separated. 3. The above diagnostic techniques are possible with spectra of relatively low quality, which can be obtained for several tens of ULIRGs with currently available telescopes. This makes L-band spectroscopy the current best tool to disentangle AGNs and starbursts contributions in ULIRGs. 4. The L-band properties of ULIRGs are heterogeneous. However, we show that all the spectral differences among ULIRGs can be reproduced starting from pure intrinsic AGN and SB spectra and varying two parameters: the amount of dust extinction of the AGN component, and the relative AGN/SB contribution to the bolometric luminosity. 5. Using the above decomposition model, we show that AGNs in ULIRGs have a low dust-to-gas ratio and a dust extinction curve different from Galactic. 6. The estimate of the contribution of AGNs in a complete sample show that AGNs are hosted by ~2/3 of ULIRGs, but their energetic contribution is relevant only in ~20% of the sample.

We present results of a morphological analysis of a small subset of the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) galaxy population. The analysis is based on public ACS data taken inside the SWIRE N1 field, which are the deepest optical high-resolution imaging available within the SWIRE fields as of today. Our reference sample includes 156 galaxies detected by both ACS and SWIRE. Among the various galaxy morphologies, we disentangle two main classes, spheroids (or bulge-dominated galaxies) and disk-dominated ones, for which we compute the number counts as a function of flux. We then limit our sample to objects with IRAC fluxes brighter than 10 microJy, estimated ~90% completeness limit of the SWIRE catalogues, and compare the observed counts to model predictions. We find that the observed counts of the spheroidal population agree with the expectations of a hierarchical model while a monolithic scenario predicts steeper counts. Both scenaria, however, under-predict the number of late-type galaxies. These observations show that the large majority (close to 80 per cent) of the 3.6 and 4.5 micron galaxy population, even at these moderately faint fluxes, is dominated by spiral and irregular galaxies or mergers.

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.

We present low resolution multi-object spectroscopy of an I-band magnitude limited (I_{AB} ~ 23--23.5) sample of galaxies located in an area centered on the Hubble Deep Field-South (HDFS). The observations were obtained using the Focal Reducer low dispersion Spectrograph (FORS) on the ESO Very Large Telescope. Thirty-two primary spectroscopic targets in the HST-WFPC2 HDFS were supplemented with galaxies detected in the Infrared Space Observatory's survey of the HDFS and the ESO Imaging Deep Survey to comprise a sample of 100 galaxies for spectroscopic observations. Based on detections of several emission lines, such as [OII]3727, H_beta and [OIII]5007, or other spectroscopic features, we have measured accurate redshifts for 50 objects in the central HDFS and flanking fields. The redshift range of the current sample of galaxies is 0.6--1.2, with a median redshift of 1.13 (at I ~ 23.5 not corrected for completeness). The sample is dominated by starburst galaxies with only a small fraction of ellipticals (~10%). For the emission line objects, the extinction corrected [OII]3727 line strengths yield estimates of star formation rates in the range 0.5--30 M_solar/yr. We have used the present data to derive the [OII]3727 luminosity function up to redshift of 1.2. When combined with [OII]3727 luminosity densities for the local and high redshift Universe, our results confirm the steep rise in the star formation rate (SFR) to z ~ 1.3.

We report on near-infrared medium-resolution spectroscopy of a sample of luminous and ultra luminous infrared galaxies (LIRGs-ULIRGs), carried out with SOFI at the ESO 3.5m New Technology Telescope. Because of wavelength dependence of the attenuation, the detection of the Pa_alfa or Br_gamma line in the Ks band should provide relevant constraints on SFR and the contribution of an AGN. We find, however, that the intensities of the Pa_alfa and Br_gamma lines, even corrected for slit losses, are on average only 10% and 40%, respectively, of that expected from a normal starburst of similar bolometric luminosity. The corresponding star formation rates, after correcting for the attenuation derived from the NIR-optical emission line ratios, are 14% and 60% of that expected if the far infrared luminosity were entirely powered by the starburst. This confirms the existence of a recombination photon deficit, particularly in the case of the Pa_alfa line, already found in the Br_gamma line in other infrared galaxies of similar luminosity. In discussing the possible causes of the discrepancy, we find unlikely that it is due to the presence of an AGN, though two objects show evidence of broadening of the Pa_alfa line and of the presence of coronal line emission. In fact, from our own observations and data collected from the literature we argue that the studied galaxies appear to be predominantly powered by a nuclear starburst. Two scenarios compatible with the present data are that either there exists a highly attenuated nuclear star forming region, and/or that a significant fraction of the ionizing photons are absorbed by dust within the HII regions. We suggest that observations in the Br_alpha spectral region could constitute a powerful tool to disentangle these two possibilities.

We exploit deep combined observations with Spitzer and Chandra of the SWIRE survey in the ELAIS-N1 region, to investigate the nature of the faint X-ray and IR sources in common, to identify AGN/starburst diagnostics, and to study the sources of the X-ray and IR cosmic backgrounds. In the 17'x17' area of the Chandra ACIS-I image there are 3400 SWIRE near-IR sources with 4 sigma detections in at least 2 IRAC bands and 988 sources detected at 24micron with MIPS brighter than 0.1 mJy. Of these, 102 IRAC and 59 MIPS sources have Chandra counterparts, out of a total of 122 X-ray sources present in the area with S(0.5-8 kev)>10^(-15) erg/cm^2/s. We have constructed SEDs for each source using data from the 4 IRAC wavebands, Chandra fluxes, and optical follow-up data in the wavebands U, g', r', i', Z, and H. We fit a number of spectral templates to the SEDs at optical and infrared wavelengths to determine photometric redshifts and spectral categories, and also make use of diagnostics based on the X-ray luminosities, hardness ratios, X-ray to infrared spectral slopes and optical morphologies. Although we have spectroscopic redshifts for only a minority of the Chandra sources, the available SEDs constrain the redshifts for most of the sample sources, which turn out to be typically at 0.5<z<2. We find that 39% of the Chandra sources are dominated by type-1 AGN emission, 23% display optical/IR spectra typical of type-2 AGNs, while the remaining 38% fraction show starburst-like or even normal galaxy spectra. Since we prove that all these galaxies are dominated by AGN emission in X-rays this brings the fraction of type-1 AGNs to be 80% of the type-2: even assuming that all the Chandra sources undetected by Spitzer are type-2 AGNs, the type-1 fraction would exceed 1/3 of the total population (abridged).

We discuss optical associations, spectral energy distributions and photometric redshifts for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The band-merged IRAC (3.6, 4.5, 5.8 and 8.0 mu) and MIPS (24, 70, 160 mu) data have been associated with optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, and with our own optical Ugri data in Lockman-VF. The spectral energy distributions of selected ELAIS sources in N1 detected by SWIRE, most with spectroscopic redshifts, are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 7%. The photometric redshifts are used to derive the 3.6 and 24 mu redshift distribution and to compare this with the predictions of models. For those sources with a clear mid infrared excess, relative to the galaxy starlight model used for the optical and near infrared, the mid and far infrared data are modelled in terms of the same dust emission templates. The proportions found of each template type are: cirrus 31%, M82 29%, Arp 220 10%, AGN dust tori 29%. The distribution of the different infrared sed types in the L_{ir}/L_{opt} versus L_{ir} plane, where L_{ir} and L_{opt} are the infrared and optical bolometric luminosities, is discussed.

We present a detailed analysis of the evolution of the rest-frame B-band morphology of K-selected galaxies with 0 < z < 2.5. This work is based on the K20 spectroscopic sample (Ks<20) located within the Chandra Deep Field South area, coupled with the public deep GOODS HST+ACSmulti-band optical imaging available in that field. Thanks to the spectroscopic completeness of this catalog reaching 94%, we can compare the morphological and spectroscopic properties of galaxies with unprecedented detail. Our morphological analysis includes visual inspection and automatic procedures using both parametric (e.g. the Sersic indices treated by the GALFIT and GASPHOT packages) and non-parametric (the Concentration, Asymmetry and clumpineSs, CAS) methods. Our analysis classifies the K20 galaxies as: 60/300 (20%, class 1) normal ellipticals/S0; 14/300 (4%, class 2) perturbed or peculiar ellipticals; 80/300 (27%, class 3) normal spirals; 48/300 (16%, class 4) perturbed or actively star-forming spirals; 98/300 (33%, class 5) irregulars. The evolution of the merging fraction is constrained up to z~2, by carefully accounting effects of morphological k-correction: both asymmetry criterion and pair statistic show an increasing merging fraction as a function of the redshift. We finally analyse the redshift-dependence of the effective radii for early- and late-type galaxies and find some mild evidence for a decrease with z of the early-type galaxy sizes, while disks and irregulars remain constant. Altogether, this analysis of the K20 sample indicates the large predominance of spirals and irregulars at 0.5<z<1.5 in K-band selected samples at even moderate depths.

We present the optical identifications of a 95 micron ISOPHOT sample in the Lockman Hole over an area of about half square degree. The catalogue (Rodighiero et al. 2003) includes 36 sources, making up a complete flux-limited sample for fluxes > 100 mJy. Reliable sources were detected, with decreasing but well-controlled completeness, down to 20 mJy. We have combined mid-IR and radio catalogues in this area to identify the potential optical counterparts of the far-IR sources. We found 14 radio and 13 15 micron associations, 10 of which have both associations. For the 11 sources with spectroscopic redshift, we have performed a spectrophotometric analysis of the observed Spectral Energy Distributions. Four of these 95 micron sources have been classified as faint IR galaxies (L_FIR<1.e11 solar luminosity), six as LIRGs and only one ULIRG. We have discussed the redshift distribution of these objects, comparing our results with evolutionary model predictions 95 and 175 micron. Given their moderate distances (the bulk of the closest spectroscopically identified objects lying at z<0.2), their luminosities and star formation rates (median value 10 solar masses/yr), the sources unveiled by ISOPHOT at 95 micron seem to correspond to the low redshift (z<0.3) FIRBACK 175 micron population, composed of dusty, star-forming galaxies with moderate star formation rates. We computed and compared different SFR estimators, and found that the SF derived from the bolometric IR luminosity is well correlated with that computed from the radio and mid-IR fluxes.

We present a mid-infrared analysis of 35 quasars with spectroscopic redshifts selected from the Spitzer Wide-area InfraRed Extragalactic Survey (SWIRE). We discuss their optical and mid-infrared (MIR) colors, and show that these quasars occupy well defined regions in MIR color-color space. We examine the issue of type-I AGN candidate selection in detail and propose new selection methods based on mid-IR colors. The available multi-band data allows us to construct two new, well-sampled quasar templates, covering wavelengths from the ultraviolet to the MIR.

We present a re--analysis of the European Large Area ISO Survey (ELAIS) 90mum observations carried out with ISOPHOT, an instrument on board the ESA's Infrared Space Observatory (ISO). With more than 12 sq. deg., the ELAIS survey is the largest area covered by ISO in a single program and is about one order of magnitude deeper than the IRAS 100mum survey. The data analysis is presented and was mainly performed with the Phot Interactive Analysis software but using the pairwise method of Stickel et al. (2003) for signal processing from ERD (Edited Raw Data) to SCP (Signal per Chopper Plateau). The ELAIS 90mum catalogue contains 229 reliable sources with fluxes larger than 70 mJy and is available at www.blackwell-synergy.com. Number counts are presented and show an excess above the no-evolution model prediction. This confirms the strong evolution detected at shorter(15mum) and longer (170mum) wavelengths in other ISO surveys. The ELAIS counts are in agreement with previous works at 90mum and in particular with the deeper counts extracted from the Lockman hole observations. Comparison with recent evolutionary models show that the models of Franceschini et al. and Guiderdoni et al. which includes a heavily-extinguished population of galaxies give the best fit to the data. Deeper observations are nevertheless required to better discriminate between the model predictions in the far-infrared and are scheduled with the Spitzer Space Telescope which already started operating and will also be performed by ASTRO-F.

We present a cross-correlation of the European Large Area ISO survey (ELAIS) with the ELAIS Deep X-ray Survey of the N1 and N2 fields. There are 7 Chandra point sources with matches in the ELAIS Final Analysis 15um catalogue, out of a total of 28 extragalactic ISO sources present in the Chandra fields. Five of these are consistent with AGN giving an AGN fraction of ~19 per cent in the 15um flux range 0.8-6 mJy. We have co-added the hard X-ray fluxes of the individually-undetected ISO sources and find a low significance detection consistent with star formation in the remaining population. We combine our point source cross-correlation fraction with the XMM-Newton observations of the Lockman Hole and Chandra observations of the Hubble Deep Field North to constrain source count models of the mid-infrared galaxy population. The low dust-enshrouded AGN fraction in ELAIS implied by the number of cross-identifications between the ELAIS mid-infrared sample and the Chandra point sources is encouraging for the use of mid-infrared surveys to constrain the cosmic star formation history, provided there are not further large undetected populations of Compton-thick AGN.

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 a new analysis of the ISOCAM 14.3 micron deep survey in a 20x20 square arcmins area in the Lockman Hole. This survey is intermediate between the ultra-deep surveys and the shallow surveys in the ELAIS fields. The data have been analyzed with the method presented by Lari et al. (2001). We have produced a catalogue of 283 sources detected above the 5-sigma threshold, with fluxes in the interval 0.1-8 mJy. The catalogue is 90% complete at 1 mJy. The positional accuracy, estimated from the cross-correlation of infrared and optical sources, is around 1.5 arcsec. The search for the optical counterparts of the sources in the survey is performed on a medium-deep r' band optical image (5-sigma depth of r'=25), making use of the radio detections when available. The photometry has been checked through simulations and by comparing the data with those presented in a shallower and more extended ISOCAM survey in the Lockman Hole, that we have presented in a companion paper. Only 15% of the 14.3 micron sources do not have an optical counterpart down to r'=25 mag. We use the 6.7/14.3 micron colour as a star/galaxy separator, together with a visual inspection of the optical image and an analysis of the observed Spectral Energy Distribution of the ISOCAM sources. The stars in the sample turn out to be only 6% of the sample. We discuss the 14.3 micron counts of extragalactic sources, combining our catalogue with that obtained from the shallower ISOCAM survey. The data in the two surveys are consistent, and our results fully support the claims in previous works for the existence of an evolving population of infrared galaxies, confirming the evident departure from non-evolutionary model predictions.

We present the image and catalogue of the 14.3 micron shallow survey of 0.55 square degrees in the region of the Lockman Hole (10h52m03s +57d21m46s, J2000) with the Infrared Space Observatory (ISO). The data have been analyzed with the recent algorithm by Lari et al. (2001) conceived to exploit ISO data in an optimal way, especially in the case of shallow surveys with low redundancy. Photometry has been accurately evaluated through extensive simulations and also the absolute calibration has been checked using a set of 21 stars detected at 14.3 micron, optical, and near-IR bands. On the basis of simulations, we evaluate that the survey is 80%, 50%, and 20% complete at 0.8, 0.6, and 0.45 mJy, respectively. Below the 20% completeness limit, fluxes are generally overestimated since the sources are preferentially detected if their positions correspond to positive oscillations of the noise. Moreover, from a comparison with the deep survey, we estimate that only sources brighter than 0.45 mJy are highly reliable. Only 5% of these sources do not have optical counterparts down to r'=25. Since none of the Spitzer imaging bands cover the 14.3 micron wavelength range, this data set will remain unique until the advent of the James Webb Space Telescope.

We report the discovery of a Compton-thick AGN and of intense star-formation activity in the nucleus and disk, respectively, of the nearly edge-on superwind galaxy NGC 4666. Spatially unresolved emission is detected by BeppoSAX only at energies <10 keV, whereas spatially resolved emission from the whole disk is detected by XMM-Newton. A prominent (EW ~ 1-2 keV) emission line at ~6.4 keV is detected by both instruments. From the XMM-Newton data alone the line is spectrally localized at E ~ 6.42 +/- 0.03 keV, and seems to be spatially concentrated in the nuclear region of NGC 4666. This, together with the presence of a flat (Gamma ~ 1.3) continuum in the nuclear region, suggests the existence of a strongly absorbed (i.e., Compton-thick) AGN, whose intrinsic 2-10 keV luminosity is estimated to be L_{2-10} > 2 x 10^{41} erg/s. At energies <1 keV the integrated (BeppoSAX) spectrum is dominated by a ~0.25 keV thermal gas component distributed throughout the disk (resolved by XMM-Newton). At energies ~2-10 keV, the integrated spectrum is dominated by a steep (G > 2) power-law (PL) component. The latter emission is likely due to unresolved sources with luminosity L ~ 10^{38} - 10^{39} erg/s that are most likely accreting binaries (with BH masses <8 M_sun). Such binaries, which are known to dominate the X-ray point-source luminosity in nearby star-forming galaxies, have Gamma ~ 2 PL spectra in the relevant energy range. A Gamma ~ 1.8 PL contribution from Compton scattering of (the radio-emitting) relativistic electrons by the ambient FIR photons may add a truly diffuse component to the 2-10 keV emission.

We characterize the SWIRE galaxy populations in the SWIRE validation field within the Lockman Hole, based on the 3.6-24$\mu$ Spitzer data and deep U,g',r',r' optical imaging within an area ~1/3 sq. deg for ~16,000 Spitzer-SWIRE sources. The entire SWIRE survey will discover over 2.3 million galaxies at 3.6$\mu$m and almost 350,000 at 24$\mu$m; ~70,000 of these will be 5-band 3.6-24$\mu$ detections. The colors cover a broad range, generally well represented by redshifted spectral energy distributions of known galaxy populations, however significant samples of unusually blue objects in the [3.6-4.5]$\mu$m color are found, as well as many objects very red in the 3.6-24$\mu$m mid-IR. Nine of these are investigated and are interpreted as star-forming systems, starbursts and AGN from z=0.37 to 2.8, with luminosities from L$_{IR}$=10$^{10.3}$ to 10$^{13.7}$ L$_{\odot}$

We report results from a new analysis of a deep 95 micron imaging survey with ISOPHOT on board the Infrared Space Observatory, over a ~1 square degree area within the Lockman Hole, which extends the statistics of our previous study (Rodighiero et al. 2003). Within the survey area we detect sixty-four sources with S/N>3 (roughly corresponding to a flux limit of 16 mJy). Extensive simulations indicate that the sample is almost complete at fluxes > 100 mJy, while the incompleteness can be quantified down to ~30 mJy. The 95 micron galaxy counts reveal a steep slope below 100 mJy (alpha~1.6), in excess of that expected for a non-evolving source population. In agreement with counts data from ISO at 15 and 175 micron, this favours a model where the IR populations evolve both in number and luminosity densities. We finally comment on some differences found with other ISO results in this area.

Analysing mid--infrared ISOCAM images of the cluster of galaxies J1888.16CL, we identified among its members several particularly active galaxies with total infrared luminosities well above 10^{11} Lsun. If powered by dust enshrouded starbursts, as suggested by their optical spectra, these Luminous Infrared Galaxies would exhibit star formation rates surprisingly high in the cluster environment. The triggering mechanism is unclear but could be tidal collisions within sub-structures or infalling groups.

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

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

Based on recent work on spectral decomposition of the emission of star-forming galaxies, we assess whether the integrated 2-10 keV emission from high-mass X-ray binaries (HMXBs), L_{2-10}^{HMXB}, can be used as a reliable estimator of ongoing star formation rate (SFR). Using a sample of 46 local (z < 0.1) star forming galaxies, and spectral modeling of ASCA, BeppoSAX, and XMM-Newton data, we demonstrate the existence of a linear SFR-L_{2-10}^{HMXB} relation which holds over ~5 decades in X-ray luminosity and SFR. The total 2-10 keV luminosity is not a precise SFR indicator because at low SFR (i.e., in normal and moderately-starbursting galaxies) it is substantially affected by the emission of low-mass X-ray binaries, which do not trace the current SFR due to their long evolution lifetimes, while at very high SFR (i.e., for very luminous FIR-selected galaxies) it is frequently affected by the presence of strongly obscured AGNs. The availability of purely SB-powered galaxies - whose 2-10 keV emission is mainly due to HMXBs - allows us to properly calibrate the SFR-L_{2-10}^{HMXB} relation. The SFR-L_{2-10}^{HMXB} relation holds also for distant (z ~ 1) galaxies in the Hubble Deep Field North sample, for which we lack spectral information, but whose SFR can be estimated from deep radio data. If confirmed by more detailed observations, it may be possible to use the deduced relation to identify distant galaxies that are X-ray overluminous for their (independently estimated) SFR, and are therefore likely to hide strongly absorbed AGNs.

We discuss XMM-Newton and BeppoSAX observations of MKN 231, the lowest-redshift Broad Absorption Line (BAL) QSO known so far and one of the best-studied Ultra Luminous Infrared Galaxies. By combining the XMM-Newton spectral resolution and the high-energy bandpass of BeppoSAX we have been able to study in more detail than previously possible its 0.2--50 keV spectral properties. The BeppoSAX PDS data unveiled, for the first time, a highly absorbed ($N_H\sim 2\times 10^{24}$cm$^{-2}$) power-law component. We find that: a) below 10 keV we are seeing only reprocessed radiation through reflection and/or scattering; b) the intrinsic 2-10 keV luminosity of MKN 231 is $1^{+1.0}_{-0.5} \times 10^{44}$ ergs s$^{-1}$, i.e. more than an order of magnitude greater than previous measurements; c) the starburst activity significantly contributes to the soft ($E<2$ keV) X-ray emission.

We present a new tool for the photometric estimate of stellar masses in distant galaxies. The observed SEDs are fitted by combining single stellar populations, with different SFRs and amounts of dust extinction. This approach gives us the best flexibility when dealing with the widest variety of physical situation for the target galaxies. In particular we tested the code on three classes of sources: dusty ISO-selected starbursts, K-band selected ellipticals/S0s and z=2-3 Lyman-break galaxies. We pay particular attention in evaluating the uncertainties in the stellar mass estimate, due to degeneracies in the physical parameters, different SFHs or metallicities. Based on optical-NIR photometric data, the stellar masses are found to have overall uncertainties of a factor of ~2 for E/S0s, ~2-5 for the starbursts population, and up to 10 for Ly-break galaxies. In any case the latter appear to correspond to a galaxy population significantly less massive than those observed at lower redshifts, possibly indicating substantial stellar build-up at z~1-2 in the field galaxy population. Using simulated deep SIRTF/IRAC observations of starbursts and Lyman-break galaxies, we investigate how an extension of the wavelength dynamic range will decrease the uncertainties in the stellar mass estimate, and find that they will reduce for both classes to factors of 2-3, good enough for statistically reliable determinations of the galaxy evolutionary mass functions.

Deep VLT optical spectroscopy, HST+ACS (GOODS) imaging and VLA observations are used to unveil the nature of a complete sample of 47 EROs with R-Ks>5 and Ks<20. The spectroscopic redshift completeness is 62%. Morphological classification was derived for each ERO through visual inspection and surface brightness profile fitting. Three main ERO morphological types are found: E/S0 galaxies (30-37%), spiral-like (24-46%) and irregular systems (17-39%). The only ERO detected in the radio is likely to host an obscured AGN. The average radio luminosity of the star-forming EROs undetected in the radio implies star formation rates of the order of 33 Msun yr^-1. The colors, redshifts and masses of the E/S0 galaxy subsample imply a minimum formation redshift z_f \sim 2. With this z_f there is enough time to have old and massive stellar spheroids already assembled at z \sim 1. We verify that the R-Ks vs. J-Ks color diagram is efficient in segregating old and dusty-star-forming EROs.

We present an L-band spectrum of the Ultraluminous Infrared Galaxy IRAS 19254-7245 (the Superantennae), obtained with VLT-ISAAC. The high signal to noise ratio allows a study of the main spectral features with unprecedented detail for an extragalactic source. We argue that the main energy source in the IR is an obscured AGN. This is indicated by the low equivalent width of the 3.3 micron PAH feature, the broad absorption feature at 3.4 um, and the steep continuum at lambda>3.7 um (f_lambda ~lambda^(2.7)). The substructure of the 3.4 um absorption feature indicates that the absorption is due to hydrocarbon chains of 6-7 carbon atoms.

We present BeppoSAX and XMM-Newton observations of MKN 231. These observations and in particular the BeppoSAX PDS data allowed us to unveil, for the first time, the highly absorbed (N$_H\sim 2\times10^{24}$ cm$^{-2}$) AGN component. We find that: a) the AGN powering MKN231 has an intrinsic 2-10 keV luminosity of $1^{+1}_{-0.5}\times 10^{44}$ erg/s; b) the strong starburst activity contributes significantly in the 0.1-10 keV energy range. We propose that the starburst activity strongly contributes to the far infrared luminosity of MKN 231; this is also suggested by the multiwavelength properties of MKN 231.

The SIRTF Wide-area InfraRed Extragalactic survey (SWIRE), the largest SIRTF Legacy program, is a wide-area, imaging survey to trace the evolution of dusty, star-forming galaxies, evolved stellar populations, and AGN as a function of environment, from redshifts z~3 to the current epoch. SWIRE will survey 7 high-latitude fields, totaling 60 - 65 sq. deg. in all 7 SIRTF bands: IRAC 3.6, 4.5, 5.6, 8 microns and MIPS 24, 70, 160 microns. The Legacy Extragalactic Catalog may contain in excess of 2 million IR-selected galaxies, dominated by (1) ~150,000 luminous infrared galaxies (LIRGs: L{FIR}>10^11 L_sun), ~7000 of these with z>2; (2) 1 million early-type galaxies, ~10,000 with z>2; and (3) \~20,000 classical AGN, plus significantly more dust-obscured QSO/AGN among the LIRGs. SWIRE will provide an unprecedented view of the evolution of galaxies, structure, and AGN. The key scientific goals of SWIRE are: (1) to determine the evolution of actively star-forming and passively evolving galaxies in order to understand the history of galaxy formation in the context of cosmic structure formation; (2) to determine the evolution of the spatial distribution and clustering of evolved galaxies, starbursts and AGN in the key redshift range, 0.5<z<3, over which much of cosmic evolution has occurred; (3) to determine the evolutionary relationship between ``normal galaxies'' and AGN, and the contribution of AGN accretion energy vs stellar nucleosynthesis to the cosmic backgrounds. SWIRE's large area is important to create statistically significant population samples over enough volume cells that we can resolve the star formation history as a function of epoch and environment The large volume is also optimised for finding rare objects.

XMM-Newton observations of 10 ULIRGs are reported. The aim is to investigate in hard X-rays a complete ULIRG sample selected from the bright IRAS 60$\mu$m catalogue. All sources are detected in X-rays, 5 of which for the first time. These observations confirm that ULIRGs are intrinsically faint X-rays sources, their observed X-ray luminosities being typically L(2-10 keV)<1E42-43 erg/s, whereas their bolometric luminosities are L>1E45 erg/s. In all sources we find evidence for thermal emission from hot plasma with kT~0.7keV, dominating the X-ray spectra below 1keV, and likely associated with a nuclear or circumnuclear starburst. This thermal emission appears uncorrelated with the FIR luminosity, suggesting that,in addition to the ongoing rate of star formation, other parameters may also affect it. The soft X-ray emission appears to be extended on a scale of ~30kpc for Mkn231 and IRAS19254-7245, possible evidence of galactic superwinds. In these 2 sources, in IRAS20551-4250 and IRAS23128-5919 we find evidence for the presence of hidden AGNs, while a minor AGN contribution may be suspected also in IRAS20100-4156. In particular, we have detected a strong Fe line at 6.4keV in the spectrum of IRAS19254-7245 and a weaker one in Mkn231, suggestive of deeply buried AGNs. For the other sources, the X-ray luminosities and spectral shapes are consistent with hot thermal plasma and X-ray binary emissions of mainly starburst origin. We find that the 2-10keV luminosities in these sources, most likely due to high-mass X-ray binaries, are correlated with L_FIR: both luminosities are good indicators of the current global SFR in the galaxy. The composite nature of ULIRGs is then confirmed, with hints for a predominance of the starburst over the AGN phenomenon in these objects.

We report results of a new analysis of a deep 95 mu imaging survey with the photo-polarimeter ISOPHOT on board the Infrared Space Observatory, over a half square degree area within the Lockman Hole. To this end we exploit a newly developed parametric algorithm able to identify and clean spurious signals induced by cosmic-rays impacts and by transient effects and non-linearities in the detectors. These results provide us with the currently deepest -- to our knowledge -- far-IR image of the extragalactic sky. Within the survey area we detect thirty-six sources with S/N>3 (corresponding to a flux of 16 mJy), making up a complete flux-limited sample around 100$ mJy. Reliable sources are detected, with decreasing but well-controlled completeness, down to 20 mJy. The source extraction process and the completeness, the photometric and astrometric accuracies of this catalogue have been tested by us with extensive simulations accounting for all the details of the procedure. We estimate source counts down to a flux of ~30 mJy, at which limit we evaluate that from 10% to 20% of the cosmic IR background has been resolved into sources. The 95mu galaxy counts reveal a steep slope at fluxes lower than 100 mJy (alpha~1.6), in excess of that expected for a non-evolving source population. The shape of these counts agrees with those determined by ISO at 15 and 175 mu, and starts setting strong constraints on the evolution models for the far-IR galaxy populations.

We have performed spatially-resolved spectroscopy of the double-nucleated Ultra-Luminous Infrared Galaxy IRAS 19254-7245, ``the Superantennae'', along the line connecting the two nuclei. These data are analysed with a spectral synthesis code, to derive the star formation and extinction properties of the galaxy. The star formation history (SFH) of the two nuclei is similarly characterized by two different main episodes: a recent burst, responsible of the observed emission lines, and an older one, occurred roughly 1 Gyr ago. We tentatively associate this bimodal SFH with a double encounter in the dynamical history of the merger. We have complemented our study with a detailed analysis of the broad band spectral energy distribution of the Superantennae, built from published photometry, providing the separate optical-to-mm SEDs of the two nuclei. Our analysis shows that: a) the southern nucleus is responsible for about 80% of the total infrared luminosity of the system, b) the L-band luminosity in the southern nucleus is dominated by the emission from an obscured AGN, providing about 40 to 50% of the bolometric flux between 8 and 1000 microns; c) the northern nucleus does not show evidence for AGN emission and appears to be in a post-starburst phase. As for the relative strengths of the AGN and starburst components, we find that, while they are comparable at FIR and sub-mm wavelengths, in the radio the Sy2 emission dominates by an order of magnitude the starburst.

We present a spectroscopic analysis of 21 galaxies with z=0.2-1.5 drawn from a 25 sq. arcmin ultra-deep 15 micron ISOCAM survey centered in the HDFS. VLT/ISAAC spectra are reported for 18 sources, aimed at detecting the redshifted Ha. Optical data come from VLT/FORS2 and NTT/EMMI. Ha line emission is detected in all the observed objects. Our analysis of emission lines, morphology, and SEDs shows evidence for AGN in only two of these sources. The Ha luminosities indicate star formation rates (SFR) between 0.5 and 20 Msun, without extinction corrections. We find good correlations between the mid-IR, the radio and Ha luminosities, confirming the mid-IR light as a good tracer of star formation (while SFR based on Ha flux show large scatter and offset, still to be understood). We have estimated the baryonic masses in stars by fitting the optical-IR SED, and found that the host galaxies of ISO sources are massive members of groups with SFR=10 to 300 Msun/yr. We have finally compared this ongoing SF activity with the already formed stellar masses to estimate the timescales t(SF) for the stellar build-up, which turn-out to be widely spread between 0.1 Gyrs to more than 10 Gyr. The faint ISOCAM galaxies appear to form a composite population, including moderately active but very massive spiral-like galaxies, and very luminous ongoing starbursts, in a continuous sequence. From the observed t(SF), assuming typical starburst timescales, we infer that, with few exceptions, only a fraction of the galactic stars can be formed in any single starburst, while several of such episodes during a protracted SF history are required for the whole galactic build-up.

We describe the first results of a programme to obtain rotation curves of z~1 disc galaxies in the near-infrared using the Ha emission line in order to study the Tully-Fisher relation. To put any observed evolution into perspective and to investigate any possible selection biases, we constructed a control sample of low redshift galaxies that had rotation velocities and images available for measuring their dynamical, photometric, and morphological properties. Compared to local objects with isophotal sizes similar to the high redshift targets, we find that our sample of galaxies with spatially resolved rotation curves, the most distant sample so far (<z>~0.9), clearly reveals a brightening of ~1.1 mag in the rest-frame B-band. The observed offset can be explained by a combination of increasing surface brightness, decreasing rotation speeds, and slightly smaller disc scale lengths of the high redshift galaxies.

Phenomenological models for evolution of dusty galaxies and E/S0 galaxies, respectively, are developed to address two major questions concenring galaxy populations in deep infrared (IR) surveys: (1) Do normal late-type galaxies or starburst galaxies (including galaxies with obscured AGNs) dominate among sources in deep IR surveys? (2) How much do E/S0 galaxies contribute to the counts in deep mid-infrared (MIR: 3 -- 20$\mu m$) surveys? Among three new models for evolution of dusty galaxies, it is assumed in Model S1 that starburst galaxies are the dominant population, and in Model S2 that normal galaxies dominate. Model S3 is an intermediate model. Comparing the model predictions with a wide range of observational data collected from the literature, we find that none of these models can be ruled out, given the uncertainties of the data. We show that the most direct method to distinguish these models is to compare the predicted color distributions of IR galaxies with observations, which will soon be available from the SWIRE survey. The models for E/S0 galaxies follow a simple passive evolution approach. Among the three E/S0 models (E1, E2 and E3) investigated in this paper, Model E2 which is specified by a peak formation redshift $z_{peak}=2$, and an e-folding formation time scale $\omega=2$ Gyr, fits the data best. This suggests a synchronization between the evolution of E/S0 galaxies and of starburst galaxies, in the sense that the peak of the formation function of E/S0s ($z_{peak}=2$) is close to the peak of the evolution functions of starburst galaxies ($z_{peak} = 1.4$). New predictions for counts and confusion limits in the SIRTF bands are presented.

We present the XMM-Newton observation of the merging system IRAS 19254-7245, also known as The Superantennae, whose southern nucleus is classified as a Seyfert 2 galaxy. The XMM-Newton data have allowed us to perform a detailed X-ray imaging and spectral analysis of this system. We clearly detect, for the first time in this system, a strong EW ~ 1.4 keV Fe emission line at 6.49+/-0.1 keV (rest-frame). The X-ray spectrum requires a soft thermal component (kT~0.9 keV; L(0.5-2) ~ 4E41 cgs), likely associated with the starburst, and a hard power-law continuum above 2 keV (observed L(2-10) ~ 4E42 cgs). We confirm the flatness of this latter component, already noted in previous ASCA data. This flatness, together with the detection of the strong Fe-Kalpha line and other broad band indicators, suggest the presence of a Compton-thick AGN with intrinsic luminosity > 1E44 cgs. We show that a Compton-thick model can perfectly reproduce the X-ray spectral properties of this object.

Using a long (=150 ksec), broad-band (0.1--40 keV) BeppoSAX observation of the merging starburst system Arp 299 (=IC 694 + NGC 3690) we found the first unambiguous evidence of the presence of a deeply buried (N_H = 2.5E24 cm-2) AGN having an intrinsic luminosity of L(0.5-100 keV) = 1.9E43 cgs. The X-ray spectral properties of this AGN are discussed in detail as well as the thermal component detected at soft X-ray energies which, most likely, is associated with the starburst.

We report on a spectro-polarimetric analysis of four southern Ultra-Luminous Infrared Galaxies (ULIRGs), aimed at constraining the presence of hidden broad AGN lines. For IRAS 19254--7245 ({\sl The Superantennae}) we find evidence for a significant level of polarized light in the H$\alpha$ line with FWHM$>$2300 km/s. Some degree of polarization is also detected in IRAS 20551-4250, though with lower significance. In the two other sources (IRAS 20100--4156 and IRAS 22491--1808) no polarized signals are detected. Although it is unclear from the present data if the origin of polarization is due to reflected light from an AGN or more simply to dichroic transmission by a dust slab, we find interesting correlation between the presence of polarized components in the optical spectra and independent evidence for AGN emissions in hard X-rays and the far-IR.

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 report the first results from a study of the internal kinematics, based on spatially resolved H_alpha velocity profiles, of three galaxies at redshift z~0.6 and one at redshift z~0.8, detected by ISOCAM in the Hubble Deep Field South. The kinematics are derived from high resolution near-infrared VLT spectroscopy. One of the galaxies is a massive spiral which possesses a very large rotational velocity of 460 km/s and contains a mass of 10^12 M_solar (within 20 kpc), significantly higher than the dynamical masses measured in most other local and high redshift spirals. Two of the galaxies comprise a counter-rotating interacting system, while the fourth is also a large spiral. The observed galaxies are representative examples of the morphologies encountered among ISOCAM galaxies. The mass-to-light (M /L_bol) ratios of ISOCAM galaxies lie between those of local luminous IR galaxies and massive spirals. We measure an offset of 1.6+/-0.3 mag in the rest frame B-band and of 0.7+/-0.3 mag in the rest frame I-band when we compare the four ISOCAM galaxies to the local Tully-Fisher B and I-band relations. We conclude that the large IR luminosity of the ISOCAM population results from a combination of large mass and efficient triggering of star formation. Since ISOCAM galaxies contribute significantly to the Cosmic Infrared Background our results imply that a relatively small number of very massive and IR luminous objects contribute significantly to the IR background and star formation activity near z~0.7.

We critically review the basic assumptions of the standard model for the synthesis of the XRB in the light of new data from ultradeep surveys by Chandra and XMM, resolving major parts of it. Important constraints come in particular from the observed z-distributions of faint hard X-ray sources -- showing large excesses at redshifts (z~0.8) much lower than expected by the synthesis models -- and from their X/optical/IR SEDs combined with the IR counts of type-II AGNs. Our analysis, although supporting the general scheme which interprets the XRB as due to absorbed AGNs with broad N_H distributions, requires major revision of the other postulate of the XRB synthesis models: the AGN unification. We argue that the unification scheme based on a simple orientation effect fails at high redshifts, where galaxy activity is induced by strong interactions and mergers among gas-rich systems. This helps explaining the observational evidence that type-I and II AGNs follow different evolutionary patterns, with type-I quasars providing a very biased trace of this activity. Combined deep X-ray and IR surveys consistently find that the universe has experienced a violent phase of galaxy activity around z~1, probably related with the assembly of massive galaxies, involving both SF and obscured AGN fueling. Roughly 10 to 20% of this activity has involved substantial AGN emission, this fraction likely reflecting the AGN/starburst duty cycle during the activation phase.

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 report the discovery of an unusual low--ionization broad absorption line quasar at z=1.776 which exhibits absorption lines from many excited states of Feii. This member of the rare class of 'FeLoBAL' quasars was serendipitously found in a mid-infrared (MIR) survey of distant clusters carried out with the ISOCAM camera. ISO J005645.1-273816 has a high MIR to UV luminosity ratio, suggesting a strong dust obscuration plus emission from very hot dust. This characteristic makes MIR surveys particularly efficient at detecting LoBAL quasars.

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

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

The SIRTF Wide-area Infrared Extragalactic (SWIRE) survey is a "Legacy Program" using 851 hours of SIRTF observing time to conduct a set of large-area (67 sq. deg. split into 7 fields) high Galactic latitude imaging surveys, achieving 5-sigma sensitivities of 0.45/2.75/17.5 mJy at 24/70/160 micron with MIPS and of 7.3/9.7/27.5/32.5 microJy at 3.6/4.5/5.8/8.0 micron with IRAC. These data will yield highly uniform source catalogs and high-resolution calibrated images, providing an unprecedented view of the universe on co-moving scales up to several hundreds Mpc and to substantial cosmological depths (z\simeq 2.5 for luminous sources). SWIRE will, for the first time, study evolved stellar systems (from IRAC data) versus active star-forming systems and AGNs (from MIPS data) in the same volume, generating catalogues with of order of 2 million infrared-selected galaxies. These fields will have extensive data at other wavebands, particularly in the optical, near-IR and X-rays. SWIRE will provide a complement to smaller and deeper observations in the SIRTF Guaranteed Time and the Legacy Program GOODS, by allowing the investigation of the effect of environment on galaxy evolution. We expand here on capabilities of SWIRE to study with IRAC the evolution of the bright end of the galaxy mass function as a function of cosmic time.

We present preliminary results of XMM-Newton observations of 5 Ultra-luminous Infrared Galaxies (ULIRGs), part of a mini-survey program dedicated to 10 ULIRGs selected from the bright IRAS sample. For 3 of them (IRAS 20551-4250, IRAS 19254-7245 and Mkn 231) we find strong evidence for the presence of a hidden AGN, while for two others (IRAS 20110-4156 IRAS 22491-1808) the S/N ratio of the data does not allow us to be conclusive. In particular, we have detected a strong Fe-K line in the X-ray spectra of IRAS19254-7245, with an equivalent width (~2 keV) suggestive that most of the energy source in this object is due to a deeply buried AGN.