The investigation of pulsars between millimetre and optical wavelengths is challenging due to the faintness of the pulsar signals and the relative low sensitivity of the available facilities compared to 100-m class telescopes operating in the centimetre band. The Kinetic Inductance Detector (KID) technology offers large instantaneous bandwidths and a high sensitivity that can help to substantially increase the ability of existing observatories at short wavelengths to detect pulsars and transient emission. To investigate the feasibility of detecting pulsars with KIDs, we observed the anomalous X-ray pulsar XTE J1810-197 with the New IRAM KIDs Array-2 (NIKA2) camera installed at the IRAM 30-m Telescope in Spain. We detected the pulsations from the pulsar with NIKA2 at its two operating frequency bands, 150 and 260 GHz ($\lambda$=2.0 and 1.15 mm, respectively). This is the first time that a pulsar is detected with a receiver based on KID technology in the millimetre band. In addition, this is the first report of short millimetre emission from XTE J1810-197 after its reactivation in December 2018, and it is the first time that the source is detected at 260 GHz, which gives us new insights into the radio emission process of the star.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

NIKA2 is a dual-band millimetric continuum camera of 2900 Kinetic Inductance Detectors (KID), operating at $150$ and $260\,\rm{GHz}$, installed at the IRAM 30-meter telescope. We present the performance assessment of NIKA2 after one year of observation using a dedicated point-source calibration method, referred to as the \emph{baseline} method. Using a large data set acquired between January 2017 and February 2018 that span the whole range of observing elevations and atmospheric conditions encountered at the IRAM 30-m telescope, we test the stability of the performance parameters. We report an instantaneous field of view (FOV) of 6.5' in diameter, filled with an average fraction of $84\%$ and $90\%$ of valid detectors at $150$ and $260\,\rm{GHz}$, respectively. The beam pattern is characterized by a FWHM of $17.6'' \pm 0.1''$ and $11.1''\pm 0.2''$, and a beam efficiency of $77\% \pm 2\%$ and $55\% \pm 3\%$ at $150$ and $260\,\rm{GHz}$, respectively. The rms calibration uncertainties are about $3\%$ at $150\,\rm{GHz}$ and $6\%$ at $260\,\rm{GHz}$. The absolute calibration uncertainties are of $5\%$ and the systematic calibration uncertainties evaluated at the IRAM 30-m reference Winter observing conditions are below $1\%$ in both channels. The noise equivalent flux density (NEFD) at $150$ and $260\,\rm{GHz}$ are of $9 \pm 1\, \rm{mJy}\cdot s^{1/2}$ and $30 \pm 3\, \rm{mJy}\cdot s^{1/2}$. This state-of-the-art performance confers NIKA2 with mapping speeds of $1388 \pm 174$ and $111 \pm 11 \,\rm{arcmin}^2\cdot \rm{mJy}^{-2}\cdot \rm{h}^{-1}$ at $150$ and $260\,\rm{GHz}$. With these unique capabilities of fast dual-band mapping at high (better that 18'') angular resolution, NIKA2 is providing an unprecedented view of the millimetre Universe.

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies.

IRAM operates two observatories - the 30-meter Telescope on Pico Veleta in Spain and NOEMA, an interferometer of ten 15-meter antennas on Plateau de Bure in France. Both observatories allow to observe at millimeter wavelengths. Here, we aim at discussing the complementary between continuum observations with NOEMA and NIKA2 at the 30m and their role at the cutting edge of research in astronomy. In particular, we will review possible synergies of continuum studies from nearby star forming regions to high redshift galaxies at cosmological distances.

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.

The potential similarity of the powering mechanisms of relativistic SNe and GRBs allowed us to make a prediction that relativistic SNe are born in environments similar to those of GRBs, that is, ones which are rich in atomic gas. Here we embark on testing this hypothesis by analysing the properties of the host galaxy NGC 3278 of the relativistic SN 2009bb. This is the first time the atomic gas properties of a relativistic SN host are provided and the first time resolved 21 cm-hydrogen-line (HI) information is provided for a host of an SN of any type in the context of the SN position. We obtained radio observations with ATCA covering the HI line, and optical integral field unit spectroscopy observations with MUSE. The atomic gas distribution of NGC 3278 is not centred on the optical galaxy centre, but instead around a third of atomic gas resides in the region close to the SN position. This galaxy has a few times lower atomic and molecular gas masses than predicted from its SFR. SN 2009bb exploded close to the region with the highest SFR density and the lowest age (~5.5 Myr, corresponding to the initial mass of the progenitor star ~36 Mo). As for GRB hosts, the gas properties of NGC 3278 are consistent with a recent inflow of gas from the intergalactic medium, which explains the concentration of atomic gas close to the SN position and the enhanced SFR. Super-solar metallicity at the position of the SN (unlike for most GRBs) may mean that relativistic explosions signal a recent inflow of gas (and subsequent star formation), and their type (GRBs or SNe) is determined either i) by the metallicity of the inflowing gas, so that metal-poor gas results in a GRB explosion and metal-rich gas results in a relativistic SN explosion without an accompanying GRB, or ii) by the efficiency of gas mixing, or iii) by the type of the galaxy.

Using spectroscopic observations taken for the VIMOS Ultra-Deep Survey (VUDS) we report here on the discovery of PCl J1001+0220, a massive proto-cluster located at $z_{spec}\sim4.57$ in the COSMOS field. The proto-cluster was initially detected as a $\sim12\sigma$ overdensity of typical star-forming galaxies in the blind spectroscopic survey of the early universe ($2<z<6$) performed by VUDS. It was further mapped using a new technique developed that statistically combines spectroscopic and photometric redshifts, the latter derived from a recent compilation of deep multi-band imaging. Through various methods, the descendant halo mass of PCl J1001+0220 is estimated to be $\log(M_{h}/M_{\odot})_{z=0}\sim14.5-15$ with a large amount of mass apparently already in place at $z\sim4.57$. Tentative evidence is found for a fractional excess of older and more massive galaxies within the proto-cluster, an observation which suggests the pervasive early onset of vigorous star formation. No evidence is found for the differences in the star formation rates of member and a matched sample of coeval field galaxies either through rest-frame ultraviolet methods or through stacking extremely deep Very Large Array 3 GHz imaging. Additionally, no evidence for pervasive strong active galactic nuclei (AGN) activity is observed. Analysis of Hubble Space Telescope images provides weak evidence for for an elevated incidence of galaxy-galaxy interaction within the proto-cluster. The spectral properties of the two samples are compared, with a definite suppression of Ly$\alpha$ seen in the average member galaxy relative to the coeval field ($f_{esc,Ly\alpha}=1.8^{+0.3}_{-1.7}$% and $4.0^{+1.0}_{-0.8}$%, respectively). This observation along with other lines of evidence leads us to infer the possible presence of a large, cool diffuse medium within the proto-cluster environment evocative of a nascent intracluster medium.

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 present an analysis of far--infrared (FIR) [CII] and [OI] fine structure line and continuum observations obtained with $Herschel$/PACS, and CO(1-0) observations obtained with the IRAM Plateau de Bure Interferometer, of Lyman Break Analogs (LBAs) at $z\sim 0.2$. The principal aim of this work is to determine the typical ISM properties of $z\sim 1-2$ Main Sequence (MS) galaxies, with stellar masses between $10^{9.5}$ and $10^{11}$ $M_{\odot}$, which are currently not easily detectable in all these lines even with ALMA and NOEMA. We perform PDR modeling and apply different IR diagnostics to derive the main physical parameters of the FIR emitting gas and dust and we compare the derived ISM properties to those of galaxies on and above the MS at different redshifts. We find that the ISM properties of LBAs are quite extreme (low gas temperature, high density and thermal pressure) with respect to those found in local normal spirals and more active local galaxies. LBAs have no [CII] deficit despite having the high specific star formation rates (sSFRs) typical of starbursts. Although LBAs lie above the local MS, we show that their ISM properties are more similar to those of high-redshift MS galaxies than of local galaxies above the main sequence. This data set represents an important reference for planning future ALMA [CII] observations of relatively low-mass MS galaxies at the epoch of the peak of the cosmic star formation.

We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.

(abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

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

We present the serendipitous ALMA detection of a faint submillimeter galaxy (SMG) lensed by a foreground z~1 galaxy. By optimizing the source detection to deblend the system, we accurately build the full spectral energy distribution of the distant galaxy from the I814 band to radio wavelengths. It is extremely red, with a I-K colour larger than 2.5. We estimate a photometric redshift of 3.28 and determine the physical parameters. The distant galaxy turns out to be magnified by the foreground lens by a factor of ~1.5, which implies an intrinsic Ks-band magnitude of ~24.5, a submillimeter flux at 870um of ~2.5 mJy and a SFR of ~150-300Msun/yr, depending on the adopted tracer. These values place our source towards the faint end of the distribution of observed SMGs, and in particular among the still few faint SMGs with a fully characterized spectral energy distribution, which allows us not only to accurately estimate its redshift but also to measure its stellar mass and other physical properties. The galaxy studied in this work is a representative of the population of faint SMGs, of which only few objects are known to date, that are undetected in optical and therefore are not typically accounted for when measuring the cosmic star formation history (SFH). This faint galaxy population thus likely represents an important and missing piece in our understanding of the cosmic SFH. Its observation and characterization is of major importance to achieve a solid picture of galaxy evolution.

We use Herschel 70 to 160um images to study the size of the far-infrared emitting region in about 400 local galaxies and quasar (QSO) hosts. The sample includes normal `main-sequence' star-forming galaxies, as well as infrared luminous galaxies and Palomar-Green QSOs, with different levels and structures of star formation. Assuming Gaussian spatial distribution of the far-infrared (FIR) emission, the excellent stability of the Herschel point spread function (PSF) enables us to measure sizes well below the PSF width, by subtracting widths in quadrature. We derive scalings of FIR size and surface brightness of local galaxies with FIR luminosity, with distance from the star-forming main-sequence, and with FIR color. Luminosities L_FIR ~ 10^11 L_Sun can be reached with a variety of structures spanning 2 dex in size. Ultraluminous L_FIR >~ 10^12 L_Sun galaxies far above the main-sequence inevitably have small Re_70um ~ 0.5 kpc FIR emitting regions with large surface brightness, and can be close to optically thick in the FIR on average over these regions. Compared to these local relations, first ALMA sizes for the dust emission regions in high redshift galaxies, measured at somewhat longer rest wavelengths, suggest larger sizes at the same IR luminosity. We report a remarkably tight relation with 0.15 dex scatter between FIR surface brightness and the ratio of [CII] 158um emission and FIR emission -- the so-called [CII]-deficit is more tightly linked to surface brightness than to FIR luminosity or FIR color. Among 33 z <~ 0.1 PG QSOs with typical L_FIR/L_Bol,AGN ~ 0.1, 19 have a measured 70um half light radius, with median Re_70um = 1.1kpc. This is consistent with the FIR size for galaxies with similar L_FIR but lacking a QSO, in accordance with a scenario where the rest FIR emission of these types of QSOs is, in most cases, due to host star formation.

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

Outflows are invoked in co-evolutionary models to link the growth of SMBH and galaxies through feedback phenomena, and from the analysis of both galaxies and Active Galactic Nuclei (AGN) samples at z$\sim1-3$, it is becoming clear that powerful winds are quite common in AGN hosts. High-resolution and high S/N observations are needed in order to uncover the physical properties of the wind through kinematics analysis. We exploited VIMOS, SINFONI and Subaru/IRCS Adaptive Optics data to study the kinematics properties on the scale the host galaxy of XID5395, a luminous, X-ray obscured Starburst/Quasar merging system at z$\sim1.5$ detected in the XMM-COSMOS field, and associated with an extreme [O II] emitter (EW$\sim200$ \AA). We mapped, for the first time, at high resolution the kinematics of the [O III] and H$\alpha$ line complexes and linked them with the [O II] emission. The high spatial resolution achieved allowed us to resolve all the components of the SB-QSO system. Our analysis with a resolution of few kpc reveals complexities and asymmetries in and around the nucleus of XID5395. The velocity field measured via non parametric analysis reveals different kinematic components, with maximum blueshifted and redshifted velocities up to $\simeq1300$ km s$^{-1}$, not spatially coincident with the nuclear core. These extreme values of the observed velocities and the spatial location can be explained by the presence of fast moving material. We also spectroscopically confirm the presence of a merging system at the same redshift of the AGN host. We propose that EW as large as $>150$ \AA\ in X-ray selected AGN may be an efficient criterion to isolate objects associated to the short, transition phase of "feedback" in the AGN-galaxy co-evolutionary path, which will subsequently evolve in an unobscured QSO, as suggested from the different observational evidences we accumulated for XID5395.

(Abridged) Combining the deepest Herschel extragalactic surveys (PEP, GOODS-H, HerMES), and Monte Carlo mock catalogs, we explore the robustness of dust mass estimates based on modeling of broad band spectral energy distributions (SEDs) with two popular approaches: Draine & Li (2007, DL07) and a modified black body (MBB). As long as the observed SED extends to at least 160-200 micron in the rest frame, M(dust) can be recovered with a >3 sigma significance and without the occurrence of systematics. An average offset of a factor ~1.5 exists between DL07- and MBB-based dust masses, based on consistent dust properties. At the depth of the deepest Herschel surveys (in the GOODS-S field) it is possible to retrieve dust masses with a S/N>=3 for galaxies on the main sequence of star formation (MS) down to M(stars)~1e10 [M(sun)] up to z~1. At higher redshift (z<=2) the same result is achieved only for objects at the tip of the MS or lying above it. Molecular gas masses, obtained converting M(dust) through the metallicity-dependent gas-to-dust ratio delta(GDR), are consistent with those based on the scaling of depletion time, and on CO spectroscopy. Focusing on CO-detected galaxies at z>1, the delta(GDR) dependence on metallicity is consistent with the local relation. We combine far-IR Herschel data and sub-mm ALMA expected fluxes to study the advantages of a full SED coverage.

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.

In order to investigate the FIR properties of radio-active AGN, we have considered three different fields where both radio and FIR observations are the deepest to-date: GOODS-South, GOODS-North and the Lockman Hole. Out of a total of 92 radio-selected AGN, ~64% are found to have a counterpart in Herschel maps. The percentage is maximum in the GOODS-North (72%) and minimum (~50%) in the Lockman Hole, where FIR observations are shallower. Our study shows that in all cases FIR emission is associated to star-forming activity within the host galaxy. Such an activity can even be extremely intense, with star-forming rates as high as ~10^3-10^4 Msun/yr. AGN activity does not inhibit star formation in the host galaxy, just as on-site star-formation does not seem to affect AGN properties, at least those detected at radio wavelengths and for z>~1. Furthermore, physical properties such as the mass and age distributions of the galaxies hosting a radio-active AGN do not seem to be affected by the presence of an ongoing star-forming event. Given the very high rate of FIR detections, we stress that this refers to the majority of the sample: most radio-active AGN are associated with intense episodes of star-formation. However, the two processes proceed independently within the same galaxy, at all redshifts but in the local universe, where powerful enough radio activity reaches the necessary strength to switch off the on-site star formation. Our data also show that for z>~1 the hosts of radio-selected star-forming galaxies and AGN are indistinguishable from each other both in terms of mass and IR luminosity distributions. The two populations only differentiate in the very local universe, whereby the few AGN which are still FIR-active are found in galaxies with much higher masses and luminosities.

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.

Local starbursts have a higher efficiency of converting gas into stars, as compared to typical star-forming galaxies at a given stellar mass, possibly indicative of different modes of star formation. With the peak epoch of galaxy formation occurring at z > 1, it remains to be established whether such an efficient mode of star formation is occurring at high-redshift. To address this issue, we measure the molecular gas content of seven high-redshift (z ~ 1.6) starburst galaxies with the Atacama Large (sub-)Millimeter Array and IRAM/Plateau de Bure Interferometer. Our targets are selected from the sample of Herschel far-infrared detected galaxies having star formation rates (~300-800 Msolar/yr) elevated (>4x) above the star-forming `main sequence', and included in the FMOS-COSMOS near-infrared spectroscopic survey of star-forming galaxies at z ~ 1.6 with Subaru. We detect CO emission in all cases at high levels of significance, indicative of high gas fractions (~30-50%). Even more compelling, we firmly establish with a clean and systematic selection that starbursts, identified as main-sequence outliers, at high redshift generally have a lower ratio of CO to total infrared luminosity as compared to typical main-sequence star-forming galaxies, although with a smaller offset than expected based on past studies of local starbursts. We put forward a hypothesis that there exists a continuous increase in star formation efficiency with elevation from the main sequence with galaxy mergers as a possible physical driver. Along with a heightened star formation efficiency, our high-redshift sample is similar in other respects to local starbursts such as being metal rich and having a higher ionization state of the interstellar medium.

We present new data for four candidate obscured Compton-Thick (CT) quasars at z $\sim$1-2.5 observed with SINFONI VLT spectrograph in AO mode. These sources were selected from a 24$\mu$m Spitzer MIPS survey of the COSMOS field, on the basis of red mid-infrared-to-optical and optical-to-near-infrared colours, with the intention of identifying active galactic nuclei (AGNs) in dust enshrouded environments, where most of the black hole mass is assembled in dust enshrouded environments. Near infrared spectra were analyzed in order to check for emission line features and to search for broad components in the [OIII]-H$\beta$ and H$\alpha$-[NII] regions. X-ray spectral analysis, radio and MIR diagnostics, and SED fitting have also been employed to study the nature of the sources. We successfully identified three objects for which we had only a photometric redshift estimate. Based on their emission line diagnostics and on ancillary multi-wavelength constraints, we find that all four targets harbor obscured AGNs. Broad profiles that could be attributed to the effects of outflows are revealed in only one target, MIRO20581. In particular, we clearly resolved a fast ($\sim$1600 km/s) and extended ($\sim$5 kpc) outflow in the [OIII]5007 emission line. This feature, the commonly used indicator for ionised outflowing gas, was sampled and detected only for this target; hence, we can not exclude the presence of outflows in the other sources. Overall, the constraints we obtain from our targets and from other comparative samples from the literature suggest that these optically faint luminous infrared galaxies, hosting obscured AGNs, may represent a brief evolutionary phase between the post-merger starburst and the unobscured QSO phases.

We present first results from the SXDF-ALMA 1.5 arcmin^2 deep survey at 1.1 mm using Atacama Large Millimeter Array (ALMA). The map reaches a 1sigma depth of 55 uJy/beam and covers 12 Halpha-selected star-forming galaxies at z = 2.19 or z=2.53. We have detected continuum emission from three of our Halpha-selected sample, including one compact star-forming galaxy with high stellar surface density, NB2315-07. They are all red in the rest-frame optical and have stellar masses of log (M*/Msun)>10.9 whereas the other blue, main-sequence galaxies with log(M*/Msun)=10.0-10.8 are exceedingly faint, <290 uJy (2sigma upper limit). We also find the 1.1 mm-brightest galaxy, NB2315-02, to be associated with a compact (R_e=0.7+-0.1 kpc), dusty star-forming component. Given high gas fraction (44^{+20}_{-8}% or 37^{+25}_{-3}%) and high star formation rate surface density (126^{+27}_{-30} Msun yr^{-1}kpc^{-2}), the concentrated starburst can within less than 50^{+12}_{-11} Myr build up a stellar surface density matching that of massive compact galaxies at z~2, provided at least 19+-3% of the total gas is converted into stars in the galaxy centre. On the other hand, NB2315-07, which already has such a high stellar surface density core, shows a gas fraction (23+-8%) and is located in the lower envelope of the star formation main-sequence. This compact less star-forming galaxy is likely to be in an intermediate phase between compact dusty star-forming and quiescent galaxies.

Long gamma-ray bursts (GRBs), among the most energetic events in the Universe, are explosions of massive and short-lived stars, so they pinpoint locations of recent star formation. However, several GRB host galaxies have recently been found to be deficient in molecular gas (H2), believed to be the fuel of star formation. Moreover, optical spectroscopy of GRB afterglows implies that the molecular phase constitutes only a small fraction of the gas along the GRB line-of-sight. Here we report the first ever 21 cm line observations of GRB host galaxies, using the Australia Telescope Compact Array, implying high levels of atomic hydrogen (HI), which suggests that the connection between atomic gas and star formation is stronger than previously thought, with star formation being potentially directly fuelled by atomic gas (or with very efficient HI-to-H2 conversion and rapid exhaustion of molecular gas), as has been theoretically shown to be possible. This can happen in low metallicity gas near the onset of star formation, because cooling of gas (necessary for star formation) is faster than the HI-to-H2 conversion. Indeed, large atomic gas reservoirs, together with low molecular gas masses, stellar and dust masses are consistent with GRB hosts being preferentially galaxies which have very recently started a star formation episode after accreting metal-poor gas from the intergalactic medium. This provides a natural route for forming GRBs in low-metallicity environments. The gas inflow scenario is also consistent with the existence of the companion HI object with no optical counterpart ~19 kpc from the GRB 060505 host, and with the fact that the HI centroids of the GRB 980425 and 060505 hosts do not coincide with optical centres of these galaxies, but are located close to the GRB positions.

We investigate the star formation properties of ~800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~0.1 - 4) and about four orders of magnitude in star formation rate (SFR). It includes both star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs population and the occurrence of AGNs in sources with different level of star formation. Finally we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.

Dust-Obscured galaxies (DOGs) are bright 24 um-selected sources with extreme obscuration at optical wavelengths. They are typically characterized by a rising power-law continuum of hot dust (T_D ~ 200-1000K) in the near-IR indicating that their mid-IR luminosity is dominated by an an active galactic nucleus (AGN). DOGs with a fainter 24 um flux display a stellar bump in the near-IR and their mid-IR luminosity appears to be mainly powered by dusty star formation. Alternatively, it may be that the mid-IR emission arising from AGN activity is dominant but the torus is sufficiently opaque to make the near-IR emission from the AGN negligible with respect to the emission from the host component. In an effort to characterize the astrophysical nature of the processes responsible for the IR emission in DOGs, this paper exploits Herschel data (PACS + SPIRE) on a sample of 95 DOGs within the COSMOS field. We derive a wealth of far-IR properties (e.g., total IR luminosities; mid-to-far IR colors; dust temperatures and masses) based on SED fitting. Of particular interest are the 24 um-bright DOGs (F_24um > 1mJy). They present bluer far-IR/mid-IR colors than the rest of the sample, unveiling the potential presence of an AGN. The AGN contribution to the total 8-1000 um flux increases as a function of the rest-frame 8 um-luminosity irrespective of the redshift. This confirms that faint DOGs (L_8 um< 10^12 L_sun) are dominated by star-formation while brighter DOGs show a larger contribution from an AGN.

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

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

We have conducted a detailed object-by-object study of a mass-complete (M*>10^11 M_sun) sample of 56 galaxies at 1.4 < z < 2 in the GOODS-South field, showing that an accurate de-blending in MIPS/24um images is essential to properly assign to each galaxy its own star formation rate (SFR), whereas an automatic procedure often fails. This applies especially to galaxies with SFRs below the Main Sequence (MS) value, which may be in their quenching phase. After that, the sample splits evenly between galaxies forming stars within a factor of 4 of the MS rate (~45%), and sub-MS galaxies with SFRs ~10-1000 times smaller (~55%). We did not find a well defined class of intermediate, transient objects below the MS, suggesting that the conversion of a massive MS galaxy into a quenched remnant may take a relatively short time (<1 Gyr), though a larger sample should be analyzed in the same way to set precise limits on the quenching timescale. X-ray detected AGNs represent a ~30% fraction of the sample, and are found among both star-forming and quenched galaxies. The morphological analysis revealed that ~50% of our massive objects are bulge-dominated, and almost all MS galaxies with a relevant bulge component host an AGN. We also found sub-MS SFRs in many bulge-dominated systems, providing support to the notion that bulge growth, AGN activity and quenching of star formation are closely related to each other.

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

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

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

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

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

[Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.

We study the evolution of the total star formation (SF) activity, total stellar mass and halo occupation distribution in massive halos by using one of the largest X-ray selected sample of galaxy groups with secure spectroscopic identification in the major blank field surveys (ECDFS, CDFN, COSMOS, AEGIS). We provide an accurate measurement of SFR for the bulk of the star-forming galaxies using very deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations. For undetected IR sources, we provide a well-calibrated SFR from SED fitting. We observe a clear evolution in the level of SF activity in galaxy groups. The total SF activity in the high redshift groups (0.5<z<1.1) is higher with respect to the low redshift (0.15<z<0.5) sample at any mass by 0.8+/-0.12 dex. A milder difference (0.35+/-0.1 dex) is observed between the low redshift bin and the groups at z~0. We show that the level of SF activity is declining more rapidly in the more massive halos than in the more common lower mass halos. We do not observe any evolution in the halo occupation distribution and total stellar mass- halo mass relations in groups. The picture emerging from our findings suggests that the galaxy population in the most massive systems is evolving faster than galaxies in lower mass halos, consistently with a "halo downsizing" scenario.

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

We combine multiwavelength data in the AEGIS-XD and C-COSMOS surveys to measure the typical dark matter halo mass of X-ray selected AGN [Lx(2-10keV)>1e42 erg/s] in comparison with far-infrared selected star-forming galaxies detected in the Herschel/PEP survey (PACS Evolutionary Probe; Lir>1e11 solar) and quiescent systems at z~1. We develop a novel method to measure the clustering of extragalactic populations that uses photometric redshift Probability Distribution Functions in addition to any spectroscopy. This is advantageous in that all sources in the sample are used in the clustering analysis, not just the subset with secure spectroscopy. The method works best for large samples. The loss of accuracy because of the lack of spectroscopy is balanced by increasing the number of sources used to measure the clustering. We find that X-ray AGN, far-infrared selected star-forming galaxies and passive systems in the redshift interval 0.6<z<1.4 are found in halos of similar mass, $\log M_{DMH}/(M_{\odot}\,h^{-1})\approx13.0$. We argue that this is because the galaxies in all three samples (AGN, star-forming, passive) have similar stellar mass distributions, approximated by the J-band luminosity. Therefore all galaxies that can potentially host X-ray AGN, because they have stellar masses in the appropriate range, live in dark matter haloes of $\log M_{DMH}/(M_{\odot}\,h^{-1})\approx13.0$ independent of their star-formation rates. This suggests that the stellar mass of X-ray AGN hosts is driving the observed clustering properties of this population. We also speculate that trends between AGN properties (e.g. luminosity, level of obscuration) and large scale environment may be related to differences in the stellar mass of the host galaxies.

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.

Gamma-ray bursts (GRBs) have been proposed as a tool for studying star formation in the Universe, so it is crucial to investigate whether their host galaxies and immediate environments are in any way special compared with other star-forming galaxies. Here we present spatially resolved maps of dust emission of the host galaxy of the closest known GRB 980425 at z=0.0085 using our new high-resolution observations from Herschel, Atacama Pathfinder Experiment (APEX), Atacama Large Millimeter Array (ALMA) and Australia Telescope Compact Array (ATCA). We modelled the spectral energy distributions of the host and of the star-forming region displaying the Wolf-Rayet signatures in the spectrum (WR region), located 800 pc from the GRB position. The host is characterised by low dust content and a high fraction of UV-visible star formation, similar to other dwarf galaxies. These galaxies are abundant in the local universe, so it is not surprising to find a GRB in one of them, assuming the correspondence between the GRB rate and star formation. The WR region contributes substantially to the host emission at the far-infrared, millimetre, and radio wavelengths and we propose that this is a consequence of its high gas density. If dense environments are also found close to the positions of other GRBs, then the ISM density should also be considered, along with metallicity, an important factor influencing whether a given stellar population can produce a GRB.

By exploiting the VLA-COSMOS and the Herschel-PEP surveys, we investigate the Far Infrared (FIR) properties of radio-selected AGN. To this purpose, from VLA-COSMOS we considered the 1537, F[1.4 GHz]>0.06 mJy sources with a reliable redshift estimate, and sub-divided them into star-forming galaxies and AGN solely on the basis of their radio luminosity. The AGN sample is complete with respect to radio selection at all z<~3.5. 832 radio sources have a counterpart in the PEP catalogue. 175 are AGN. Their redshift distribution closely resembles that of the total radio-selected AGN population, and exhibits two marked peaks at z~0.9 and z~2.5. We find that the probability for a radio-selected AGN to be detected at FIR wavelengths is both a function of radio power and redshift, whereby powerful sources are more likely to be FIR emitters at earlier epochs. This is due to two distinct effects: 1) at all radio luminosities, FIR activity monotonically increases with look-back time and 2) radio activity of AGN origin is increasingly less effective at inhibiting FIR emission. Radio-selected AGN with FIR emission are preferentially located in galaxies which are smaller than those hosting FIR-inactive sources. Furthermore, at all z<~2, there seems to be a preferential (stellar) mass scale M ~[10^{10}-10^{11}] Msun which maximizes the chances for FIR emission. We find such FIR (and MIR) emission to be due to processes indistinguishable from those which power star-forming galaxies. It follows that radio emission in at least 35% of the entire AGN population is the sum of two contributions: AGN accretion and star-forming processes within the host galaxy.

We study a hundred of galaxies from the spectroscopic Sloan Digital Sky Survey with individual detections in the Far-Infrared Herschel PACS bands (100 or 160 $\mu$m) and in the GALEX Far-UltraViolet band up to z$\sim$0.4 in the COSMOS and Lockman Hole fields. The galaxies are divided into 4 spectral and 4 morphological types. For the star forming and unclassifiable galaxies we calculate dust extinctions from the UV slope, the H$\alpha$/H$\beta$ ratio and the $L_{\rm IR}/L_{\rm UV}$ ratio. There is a tight correlation between the dust extinction and both $L_{\rm IR}$ and metallicity. We calculate SFR$_{total}$ and compare it with other SFR estimates (H$\alpha$, UV, SDSS) finding a very good agreement between them with smaller dispersions than typical SFR uncertainties. We study the effect of mass and metallicity, finding that it is only significant at high masses for SFR$_{H\alpha}$. For the AGN and composite galaxies we find a tight correlation between SFR and L$_{IR}$ ($\sigma\sim$0.29), while the dispersion in the SFR - L$_{UV}$ relation is larger ($\sigma\sim$0.57). The galaxies follow the prescriptions of the Fundamental Plane in the M-Z-SFR space.

Until recently, dust emission has been detected in very few host galaxies of gamma-ray bursts (GRBHs). With Herschel, we have now observed 17 GRBHs up to redshift z~3 and detected seven of them at infrared (IR) wavelengths. This relatively high detection rate (41%) may be due to the composition of our sample which at a median redshift of 1.1 is dominated by the hosts of dark GRBs. Although the numbers are small, statistics suggest that dark GRBs are more likely to be detected in the IR than their optically-bright counterparts. Combining our IR data with optical, near-infrared, and radio data from our own datasets and from the literature, we have constructed spectral energy distributions (SEDs) which span up to 6 orders of magnitude in wavelength. By fitting the SEDs, we have obtained stellar masses, dust masses, star-formation rate (SFR), and extinctions for our sample galaxies. We find that GRBHs are galaxies that tend to have a high specfic SFR (sSFR), and like other star-forming galaxies, their ratios of dust-to-stellar mass are well correlated with sSFR. We incorporate our Herschel sample into a larger compilation of GRBHs, and compare this combined sample to SFR-weighted median stellar masses of the widest, deepest galaxy survey to date. This is done in order to establish whether or not GRBs can be used as an unbiased tracer of cosmic comoving SFR density (SFRD) in the universe. In contrast with previous results, this comparison shows that GRBHs are medium-sized galaxies with relatively high sSFRs; stellar masses and sSFRs of GRBHs as a function of redshift are similar to what is expected for star-forming galaxy populations at similar redshifts. We conclude that there is no strong evidence that GRBs are biased tracers of SFRD; thus they should be able to reliably probe the SFRD to early epochs.

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

Three billion years after the big bang (at redshift z=2), half of the most massive galaxies were already old, quiescent systems with little to no residual star formation and extremely compact with stellar mass densities at least an order of magnitude larger than in low redshift ellipticals, their descendants. Little is known about how they formed, but their evolved, dense stellar populations suggest formation within intense, compact starbursts 1-2 Gyr earlier (at 3<z<6). Simulations show that gas-rich major mergers can give rise to such starbursts which produce dense remnants. Sub-millimeter selected galaxies (SMGs) are prime examples of intense, gas-rich, starbursts. With a new, representative spectroscopic sample of compact quiescent galaxies at z=2 and a statistically well-understood sample of SMGs, we show that z=3-6 SMGs are consistent with being the progenitors of z=2 quiescent galaxies, matching their formation redshifts and their distributions of sizes, stellar masses and internal velocities. Assuming an evolutionary connection, their space densities also match if the mean duty cycle of SMG starbursts is 42 (+40/-29) Myr (consistent with independent estimates), which indicates that the bulk of stars in these massive galaxies were formed in a major, early surge of star-formation. These results suggests a coherent picture of the formation history of the most massive galaxies in the universe, from their initial burst of violent star-formation through their appearance as high stellar-density galaxy cores and to their ultimate fate as giant ellipticals.

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

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

We study the incidence of nuclear obscuration on a complete sample of 1310 AGN selected on the basis of their rest-frame 2-10 keV X-ray flux from the XMM-COSMOS survey, in the redshift range 0.3<z<3.5. We classify the AGN as obscured or un-obscured on the basis of either the optical spectral properties and the overall SED or the shape of the X-ray spectrum. The two classifications agree in about 70% of the objects, and the remaining 30% can be further subdivided into two distinct classes: at low luminosities X-ray un-obscured AGN do not always show signs of broad lines or blue/UV continuum emission in their optical spectra, most likely due to galaxy dilution effects; at high luminosities broad line AGN may have absorbed X-ray spectra, which hints at an increased incidence of small-scale (sub-parsec) dust-free obscuration. We confirm that the fraction of obscured AGN is a decreasing function of the intrinsic X-ray luminosity, while the incidence of absorption shows significant evolution only for the most luminous AGN, which appear to be more commonly obscured at higher redshift. We find no significant difference between the mean stellar masses and star formation rates of obscured and un-obscured AGN hosts. We conclude that the physical state of the medium responsible for obscuration in AGN is complex, and mainly determined by the radiation environment (nuclear luminosity) in a small region enclosed within the gravitational sphere of influence of the central black hole, but is largely insensitive to the wider scale galactic conditions.

We investigate the evolution of the star formation rate (SFR)-density relation in the Extended Chandra Deep Field South (ECDFS) and the Great Observatories Origin Deep Survey (GOODS) fields up to z~1.6. In addition to the "traditional method", in which the environment is defined according to a statistical measurement of the local galaxy density, we use a "dynamical" approach, where galaxies are classified according to three different environment regimes: group, "filament-like", and field. Both methods show no evidence of a SFR-density reversal. Moreover, group galaxies show a mean SFR lower than other environments up to z~1, while at earlier epochs group and field galaxies exhibit consistent levels of star formation (SF) activity. We find that processes related to a massive dark matter halo must be dominant in the suppression of the SF below z~1, with respect to purely density-related processes. We confirm this finding by studying the distribution of galaxies in different environments with respect to the so-called Main Sequence (MS) of star-forming galaxies. Galaxies in both group and "filament-like" environments preferentially lie below the MS up to z~1, with group galaxies exhibiting lower levels of star-forming activity at a given mass. At z>1, the star-forming galaxies in groups reside on the MS. Groups exhibit the highest fraction of quiescent galaxies up to z~1, after which group, "filament-like", and field environments have a similar mix of galaxy types. We conclude that groups are the most efficient locus for star-formation quenching. Thus, a fundamental difference exists between bound and unbound objects, or between dark matter haloes of different masses.

Quasi-stellar objects (QSOs) occur in galaxies in which supermassive black holes (SMBHs) are growing substantially through rapid accretion of gas. Many popular models of the co-evolutionary growth of galaxies and SMBHs predict that QSOs are also sites of substantial recent star formation, mediated by important processes, such as major mergers, which rapidly transform the nature of galaxies. A detailed study of the star-forming properties of QSOs is a critical test of such models. We present a far-infrared Herschel/PACS study of the mean star formation rate (SFR) of a sample of spectroscopically observed QSOs to z~2 from the COSMOS extragalactic survey. This is the largest sample to date of moderately luminous AGNs studied using uniform, deep far-infrared photometry. We study trends of the mean SFR with redshift, black hole mass, nuclear bolometric luminosity and specific accretion rate (Eddington ratio). To minimize systematics, we have undertaken a uniform determination of SMBH properties, as well as an analysis of important selection effects within spectroscopic QSO samples that influence the interpretation of SFR trends. We find that the mean SFRs of these QSOs are consistent with those of normal massive star-forming galaxies with a fixed scaling between SMBH and galaxy mass at all redshifts. No strong enhancement in SFR is found even among the most rapidly accreting systems, at odds with several co-evolutionary models. Finally, we consider the qualitative effects on mean SFR trends from different assumptions about the star-forming properties of QSO hosts and redshift evolution of the SMBH-galaxy relationship. While limited currently by uncertainties, valuable constraints on AGN-galaxy co-evolution can emerge from our approach.

We report the serendipitous detection of a significant overdensity of Herschel-SPIRE 250 micron sources in the vicinity of MRC1138-26. We use an adaptive kernel density estimate to quantify the significance, including a comparison with other fields. The overdensity has a size of ~3.5-4' and stands out at ~5sigma with respect to the background estimate. No features with similar significance were found in four extragalactic control fields: GOODS-North, Lockman, COSMOS and UDS. The chance of having a similar overdensity in a field with the same number but randomly distributed sources is less than 2%. The clump is also visible as a low surface brightness feature in the Planck 857 GHz map. We detect 76 sources at 250 micron (with a signal-to-noise ratio greater than 3), in a region of 4' radius; 43 of those are above a flux density limit of 20 mJy. This is a factor of 3.6 in excess over the average in the four control fields, considering only the sources above 20 mJy. We also find an excess in the number counts of sources with 250 micron flux densities between 30 and 40 mJy, compared to deep extragalactic blank-field number counts. Assuming a fixed dust temperature (30 K) and emissivity (beta=1.5) a crude, blackbody-derived redshift distribution, zBB, of the detected sources is significantly different from the distributions in the control fields and exhibits a significant peak at zBB ~ 1.5, although the actual peak redshift is highly degenerate with the temperature. We tentatively suggest, based on zBB and the similar S250/S350 colours of the sources within the peak, that a significant fraction of the sources in the clump may be at a similar redshift. Since the overdensity lies ~7' south of the z=2.16 Spiderweb protocluster MRC1138-26, an intriguing possibility (that is presently unverifiable given the data in hand) is that it lies within the same large-scale structure.(abridged)

We combine IRAM Plateau de Bure Interferometer and Herschel PACS and SPIRE measurements to study the dust and gas contents of high-redshift star forming galaxies. We present new observations for a sample of 17 lensed galaxies at z=1.4-3.1, which allow us to directly probe the cold ISM of normal star-forming galaxies with stellar masses of ~10^10Msun, a regime otherwise not (yet) accessible by individual detections in Herschel and molecular gas studies. The lensed galaxies are combined with reference samples of sub-millimeter and normal z~1-2 star-forming galaxies with similar far-infrared photometry to study the gas and dust properties of galaxies in the SFR-M*-redshift parameter space. The mean gas depletion timescale of main sequence galaxies at z>2 is measured to be only ~450Myr, a factor of ~1.5 (~5) shorter than at z=1 (z=0), in agreement with a (1+z)^-1 scaling. The mean gas mass fraction at z=2.8 is 40+/-15% (44% after incompleteness correction), suggesting a flattening or even a reversal of the trend of increasing gas fractions with redshift recently observed up to z~2. The depletion timescale and gas fractions of the z>2 normal star-forming galaxies can be explained under the "equilibrium model" for galaxy evolution, in which the gas reservoir of galaxies is the primary driver of the redshift evolution of specific star formation rates. Due to their high star formation efficiencies and low metallicities, the z>2 lensed galaxies have warm dust despite being located on the star formation main sequence. At fixed metallicity, they also have a gas-to-dust ratio 1.7 times larger than observed locally when using the same standard techniques, suggesting that applying the local calibration of the relation between gas-to-dust ratio and metallicity to infer the molecular gas mass of high redshift galaxies may lead to systematic differences with CO-based estimates.

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

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

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$\beta$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

In the local Universe, galaxy properties show a strong dependence on environment. In cluster cores, early type galaxies dominate, whereas star-forming galaxies are more and more common in the outskirts. At higher redshifts and in somewhat less dense environments (e.g. galaxy groups), the situation is less clear. One open issue is that of whether and how the star formation rate (SFR) of galaxies in groups depends on the distance from the centre of mass. To shed light on this topic, we have built a sample of X-ray selected galaxy groups at 0<z<1.6 in various blank fields (ECDFS, COSMOS, GOODS). We use a sample of spectroscopically confirmed group members with stellar mass M >10^10.3 M_sun in order to have a high spectroscopic completeness. As we use only spectroscopic redshifts, our results are not affected by uncertainties due to projection effects. We use several SFR indicators to link the star formation (SF) activity to the galaxy environment. Taking advantage of the extremely deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations, we have an accurate, broad-band measure of the SFR for the bulk of the star-forming galaxies. We use multi-wavelength SED fitting techniques to estimate the stellar masses of all objects and the SFR of the MIPS and PACS undetected galaxies. We analyse the dependence of the SF activity, stellar mass and specific SFR on the group-centric distance, up to z~1.6, for the first time. We do not find any correlation between the mean SFR and group-centric distance at any redshift. We do not observe any strong mass segregation either, in agreement with predictions from simulations. Our results suggest that either groups have a much smaller spread in accretion times with respect to the clusters and that the relaxation time is longer than the group crossing time.

We describe the far-infrared (FIR; rest-frame 8--1000\mu m) properties of a sample of 443 H\alpha-selected star-forming galaxies in the COSMOS and UDS fields detected by the HiZELS imaging survey. Sources are identified using narrow-band filters in combination with broad-band photometry to uniformly select H\alpha\ (and [OII] if available) emitters in a narrow redshift slice at z = 1.47+/-0.02. We use a stacking approach in Spitzer, Herschel (from PEP and HerMES surveys) and AzTEC images to describe their typical FIR properties. We find that HiZELS galaxies with observed H\alpha\ luminosities of ~ 10^{8.1-9.1} Lo have bolometric FIR luminosities of typical LIRGs, L_FIR ~ 10^{11.48+/-0.05} Lo. Combining the H\alpha\ and FIR luminosities, we derive median SFR = 32+/-5 Mo/yr and H\alpha\ extinctions of A(H\alpha) = 1.0+/-0.2 mag. Perhaps surprisingly, little difference is seen in typical HiZELS extinction levels compared to local star-forming galaxies. We confirm previous empirical stellar mass (M*) to A(H\alpha) relations and the little or no evolution up to z = 1.47. For HiZELS galaxies, we provide an empirical parametrisation of the SFR as a function of (u-z)_rest colours and 3.6\mu m photometry. We find that the observed H\alpha\ luminosity is a dominant SFR tracer when (u-z)_rest ~< 0.9 mag or when 3.6\mu m photometry > 22 mag (Vega) or when M* < 10^9.7 Mo. We do not find any correlation between the [OII]/H\alpha\ and FIR luminosity, suggesting that this emission line ratio does not trace the extinction of the most obscured star-forming regions. The luminosity-limited HiZELS sample tends to lie above of the so-called `main sequence' for star-forming galaxies, especially at low M*. This work suggests that obscured star formation is linked to the assembly of M*, with deeper potential wells in massive galaxies providing dense, heavily obscured environments in which stars can form rapidly.

We use deep far-infrared data from the PEP/GOODS-Herschel surveys and rest frame ultraviolet photometry to study the evolution of the molecular gas mass function of normal star forming galaxies. Computing the molecular gas mass, M(mol), by scaling star formation rates (SFR) through depletion timescales, or combining IR luminosity and obscuration properties as in Nordon et al., we obtain M(mol) for roughly 700, z=0.2-3.0 galaxies near the star forming "main sequence". The number density of galaxies follows a Schechter function of M(mol). The characteristic mass M* is found to strongly evolve up to z~1, and then to flatten at earlier epochs, resembling the infrared luminosity evolution of similar objects. At z~1, our result is supported by an estimate based on the stellar mass function of star forming galaxies and gas fraction scalings from the PHIBSS survey. We compare our measurements to results from current models, finding better agreement with those that are treating star formation laws directly rather than in post-processing. Integrating the mass function, we study the evolution of the M(mol) density and its density parameter Omega(mol).

We explore the question of whether low and moderate luminosity Active Galactic Nuclei (AGNs) are preferentially found in galaxies that are undergoing a transition from active star formation to quiescence. This notion has been suggested by studies of the UV-to-optical colors of AGN hosts, which find them to be common among galaxies in the so-called "Green Valley", a region of galaxy color space believed to be composed mostly of galaxies undergoing star-formation quenching. Combining the deepest current X-ray and Herschel. PACS far-infrared (FIR) observations of the two Chandra Deep Fields (CDFs) with redshifts, stellar masses and rest-frame photometry derived from the extensive and uniform multi-wavelength data in these fields, we compare the rest-frame U-V color distributions and SFR distributions of AGNs and carefully constructed samples of inactive control galaxies. The UV-to-optical colors of AGNs are consistent with equally massive inactive galaxies at redshifts out to z~2, but we show that such colors are poor tracers of star formation. While the FIR distributions of both star-forming AGNs and star-forming inactive galaxies are statistically similar, we show that AGNs are preferentially found in star-forming host galaxies, or, in other words, AGNs are less likely to be found in weakly star-forming or quenched galaxies. We postulate that, among X-ray selected AGNs of low and moderate accretion luminosities, the supply of cold gas primarily determines the accretion rate distribution of the nuclear black holes.

Star-formation in the galaxy populations of local massive clusters is reduced with respect to field galaxies, and tends to be suppressed in the core region. Indications of a reversal of the star-formation--density relation have been observed in a few z >1.4 clusters. Using deep imaging from 100-500um from PACS and SPIRE onboard Herschel, we investigate the infrared properties of spectroscopic and photo-z cluster members, and of Halpha emitters in XMMU J2235.3-2557, one of the most massive, distant, X-ray selected clusters known. Our analysis is based mostly on fitting of the galaxies spectral energy distribution in the rest-frame 8-1000um. We measure total IR luminosity, deriving star formation rates (SFRs) ranging from 89-463 Msun/yr for 13 galaxies individually detected by Herschel, all located beyond the core region (r >250 kpc). We perform a stacking analysis of nine star-forming members not detected by PACS, yielding a detection with SFR=48 Msun/yr. Using a color criterion based on a star-forming galaxy SED at the cluster redshift we select 41 PACS sources as candidate star-forming cluster members. We characterize a population of highly obscured SF galaxies in the outskirts of XMMU J2235.3-2557. We do not find evidence for a reversal of the SF-density relation in this massive, distant cluster.

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

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

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

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.

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) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

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

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

We use the infrared excess (IRX) FIR/UV luminosity ratio to study the relation between the effective UV attenuation (A_IRX) and the UV spectral slope (beta) in a sample of 450 1<z<2.5 galaxies. The FIR data is from very deep Herschel observations in the GOODS fields that allow us to detect galaxies with SFRs typical of galaxies with log(M)>9.3. Thus, we are able to study galaxies on and even below the main SFR-stellar mass relation (main sequence). We find that main sequence galaxies form a tight sequence in the IRX--beta plane, which has a flatter slope than commonly used relations. This slope favors a SMC-like UV extinction curve, though the interpretation is model dependent. The scatter in the IRX-beta plane, correlates with the position of the galaxies in the SFR-M plane. Using a smaller sample of galaxies with CO gas masses, we study the relation between the UV attenuation and the molecular gas content. We find a very tight relation between the scatter in the IRX-beta plane and the specific attenuation (S_A), a quantity that represents the attenuation contributed by the molecular gas mass per young star. S_A is sensitive to both the geometrical arrangement of stars and dust, and to the compactness of the star forming regions. We use this empirical relation to derive a method for estimating molecular gas masses using only widely available integrated rest-frame UV and FIR photometry. The method produces gas masses with an accuracy between 0.12-0.16 dex in samples of normal galaxies between z~0 and z~1.5. Major mergers and sub-millimeter galaxies follow a different S_A relation.

We investigate the effect of the "high-pass filter" data reduction technique on the Herschel PACS PSF and noise of the PACS maps at the 70, 100 and 160 um bands and in medium and fast scan speeds. This branch of the PACS Photometer pipeline is the most used for cosmological observations and for point-source observations.The calibration of the flux loss due to the median removal applied by the PACS pipeline (high-pass filter) is done via dedicated simulations obtained by "polluting" real PACS timelines with fake sources at different flux levels. The effect of the data reduction parameter settings on the final map noise is done by using selected observations of blank fields with high data redundancy. We show that the running median removal can cause significant flux losses at any flux level. We analyse the advantages and disadvantages of several masking strategies and suggest that a mask based on putting circular patches on prior positions is the best solution to reduce the amount of flux loss. We provide a calibration of the point-source flux loss for several masking strategies in a large range of data reduction parameters, and as a function of the source flux. We also show that, for stacking analysis, the impact of the high-pass filtering effect is to reduce significantly the clustering effect. The analysis of the global noise and noise components of the PACS maps shows that the dominant parameter in determining the final noise is the high-pass filter width. We also provide simple fitting functions to build the error map from the coverage map and to estimate the cross-correlation correction factor in a representative portion of the data reduction parameter space.

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

Deep Herschel imaging and 12CO(2-1) line luminosities from the IRAM PdBI are combined for a sample of 17 galaxies at z>1 from the GOODS-N field. The sample includes galaxies both on and above the main sequence (MS) traced by star-forming galaxies in the SFR-M* plane. The far-infrared data are used to derive dust masses, Mdust. Combined with an empirical prescription for the dependence of the gas-to-dust ratio on metallicity (GDR), the CO luminosities and Mdust values are used to derive for each galaxy the CO-to-H2 conversion factor, alpha_co. Like in the local Universe, the value of alpha_co is a factor of ~5 smaller in starbursts compared to normal star-forming galaxies (SFGs). We also uncover a relation between alpha_co and dust temperature (Tdust; alpha_co decreasing with increasing Tdust) as obtained from modified blackbody fits to the far-infrared data. While the absolute normalization of the alpha_co(Tdust) relation is uncertain, the global trend is robust against possible systematic biases in the determination of Mdust, GDR or metallicity. Although we cannot formally distinguish between a step and a smooth evolution of alpha_co with the dust temperature, we can conclude that in galaxies of near-solar metallicity, a critical value of Tdust=30K can be used to determine whether the appropriate alpha_co is closer to the starburst value (1.0 Msun(K kms pc^2)^-1, if Tdust>30K) or closer to the Galactic value (4.35 Msun (K kms pc^2)^-1, if Tdust<30K). This indicator has the great advantage of being less subjective than visual morphological classifications of mergers/SFGs, which can be difficult at high z because of the clumpy nature of SFGs. In the absence of far-infrared data, the offset of a galaxy from the main sequence (i.e., log[SSFR(galaxy)/SSFR_MS(M*,z)]) can be used to identify galaxies requiring the use of an alpha_co conversion factor lower than the Galactic value.

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.

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.

Using the wide multi-band photometry available in the COSMOS field we explore the host galaxy properties of a large sample of Active Galactic Nuclei (AGN) obtained by combining X-ray and optical spectroscopic selections. Based on a careful study of their Spectral Energy Distribution (SED), which has been parametrized using a 2-component (AGN+galaxy) model fit, we derived dust-corrected rest-frame magnitudes, colors, stellar masses and star formation rates (SFRs). We find that AGN hosts span a large range of stellar masses and SFRs. No color-bimodality is seen at any redshift in the AGN hosts, which are found to be mainly massive, red galaxies. Once accounting for the color-mass degeneracy in well defined mass-matched samples, we find a residual marginal enhancement of AGN incidence in redder galaxies with lower specific star formation rates, and we argue that this result might emerge because of our ability to properly account for AGN light contamination and dust extinction. Interestingly, we find that the probability for a galaxy to host a black hole growing at any given "specific accretion rate" (i.e. the ratio of X-ray luminosity to the host stellar mass) is almost independent of the host galaxy mass, while it decreases as a power-law with Lx/M. By analyzing the normalization of such probability distribution, we show how the incidence of AGN increases with redshift as rapidly as (1+z)^4, in close resemblance with the overall evolution of the specific star formation rate of the entire galaxy population. Although AGN activity and star formation in galaxies do appear to have a common triggering mechanism, at least in a statistical sense, within the COSMOS sample we do not find strong evidence of any 'smoking gun' signaling powerful AGN influence on the star-forming properties of their hosts galaxies.

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

We report results from a deep Jansky Very Large Array (JVLA) search for CO 1-0 line emission from galaxies in a candidate galaxy cluster at z~1.55 in the COSMOS field. We target 4 galaxies with optical spectroscopic redshifts in the range z=1.47-1.59. Two of these 4 galaxies, ID51613 and ID51813, are nominally detected in CO line emission at the 3-4 sigma level. We find CO luminosities of 2.4x10^10 K km/s pc^2 and 1.3x10^10 K km/s pc^2, respectively. Taking advantage from the clustering and 2-GHz bandwidth of the JVLA, we perform a search for emission lines in the proximity of optical sources within the field of view of our observations. We limit our search to galaxies with K<23.5 (AB) and z_phot=1.2-1.8. We find 2 bright optical galaxies to be associated with significant emission line peaks (>4 sigma) in the data cube, which we identify with the CO line emission. To test the reliability of the line peaks found, we performed a parallel search for line peaks using a Bayesian inference method. Monte Carlo simulations show that such associations are statistically significant, with probabilities of chance association of 3.5% and 10.7% for ID 51207 and ID 51380, respectively. Modeling of their optical/IR SEDs indicates that the CO detected galaxies and candidates have stellar masses and SFRs in the range (0.3-1.1)x10^11 M_sun and 60-160 M_sun/yr, with SFEs comparable to that found in other star-forming galaxies at similar redshifts. By comparing the space density of CO emitters derived from our observations with the space density derived from previous CO detections at z~1.5, and with semi-analytic predictions for the CO luminosity function, we suggest that the latter tend to underestimate the number of CO galaxies detected at high-redshift. Finally, we argue about the benefits of future blind CO searches in clustered fields with upcoming submm/radio facilities.

We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 galaxies with log(M_*/M_sun) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deepest H-band Hubble Space Telescope survey of its type prior to the installation of WFC3. We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Survey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median log(M_*/M_sun) = 10.8 and z = 2.0, associated with 250 um sources detected with SNR > 3. We derive mean total IR luminosity log L_IR (L_sun) = 12.36 +/- 0.05 and corresponding star formation rate SFR_(IR+UV) = (280 +/- 40) M_sun/yr for these objects, and find them to have mean dust temperature T_dust ~ 35 K. We find that the SFR derived from the far-IR photometry combined with UV-based estimates of unobscured SFR for these galaxies is on average more than a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, although we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis to measure the mean SFR in bins of stellar mass. We obtain detections at the 2-4 sigma level at SPIRE wavelengths for samples with log(M_*/M_sun) > 10. In contrast to the Herschel detected GNS galaxies, we find that estimates of SFR_(IR+UV) for the stacked samples are comparable to those derived from extinction corrected UV emission, although the uncertainties are large. We find evidence for an increasing fraction of dust obscured star formation with stellar mass, finding SFR_IR/SFR_UV \propto M_*^{0.7 +/- 0.2}, which is likely a consequence of the mass--metallicity relation.

We present the results of a Herschel-PACS study of a sample of 97 LINERs at redshift z\sim 0.3 selected from the zCOSMOS survey. Of these sources, 34 are detected in a least one PACS band, enabling reliable estimates of the far-infrared L(FIR) luminosities, and a comparison to the FIR luminosities of local LINERs. Many of our PACS-detected LINERs are also UV sources detected by Galex. Assuming that the FIR is produced in young dusty star-forming regions, the typical star-formation rates (SFRs) for the host galaxies in our sample is \sim 10 M_Sun yr-1, 1-2 orders of magnitude larger than in many local LINERs. Given stellar masses inferred from optical/NIR photometry of the (unobscured) evolved stellar populations, we find that the entire sample lies close to the star-forming "main sequence" for galaxies at redshift 0.3. For young star-forming regions, the H\alpha- and UV-based estimates of the SFRs are much smaller than the FIR-based estimates, by factors \sim 30, even assuming that all of the H\alpha\ emission is produced by O-star ionization rather than by the AGNs. These discrepancies may be due to large (and uncertain) extinctions towards the young stellar systems. Alternatively, the H\alpha\ and UV emissions could be tracing residual star-formation in an older less obscured population with decaying star formation. We also compare L(SF) and L(AGN) in local LINERs and in our sample and comment on the problematic use of several line diagnostic diagrams in cases similar to the sample under study.

We use the photometric redshift method of Chakrabarti & McKee (2008) to infer photometric redshifts of submillimeter galaxies with far-IR (FIR) $\it{Herschel}$ data obtained as part of the PACS Evolutionary Probe (PEP) program. For the sample with spectroscopic redshifts, we demonstrate the validity of this method over a large range of redshifts ($ 4 \ga z \ga 0.3$) and luminosities, finding an average accuracy in $(1+z_{\rm phot})/(1+z_{\rm spec})$ of 10%. Thus, this method is more accurate than other FIR photometric redshift methods. This method is different from typical FIR photometric methods in deriving redshifts from the light-to-gas mass ($L/M$) ratio of infrared-bright galaxies inferred from the FIR spectral energy distribution (SED), rather than dust temperatures. Once the redshift is derived, we can determine physical properties of infrared bright galaxies, including the temperature variation within the dust envelope, luminosity, mass, and surface density. We use data from the GOODS-S field to calculate the star formation rate density (SFRD) of sub-mm bright sources detected by AzTEC and PACS. The AzTEC-PACS sources, which have a threshold $850 \micron$ flux $\ga 5 \rm mJy$, contribute 15% of the SFRD from all ULIRGs ($L_{\rm IR} \ga 10^{12} L_{\odot}$), and 3% of the total SFRD at $z \sim 2$.

Bolometric luminosities and Eddington ratios of both X-ray selected broad-line (Type-1) and narrow-line (Type-2) AGN from the XMM-Newton survey in the COSMOS field are presented. The sample is composed by 929 AGN (382 Type-1 AGN and 547 Type-2 AGN) and it covers a wide range of redshifts, X-ray luminosities and absorbing column densities. About 65% of the sources are spectroscopically identified as either Type-1 or Type-2 AGN (83% and 52% respectively), while accurate photometric redshifts are available for the rest of the sample. The study of such a large sample of X-ray selected AGN with a high quality multi-wavelength coverage from the far-infrared (now with the inclusion of Herschel data at 100 micron and 160 micron) to the optical-UV allows us to obtain accurate estimates of bolometric luminosities, bolometric corrections and Eddington ratios. The kbol-Lbol relations derived in the present work are calibrated for the first time against a sizable AGN sample, and rely on observed redshifts, X-ray luminosities and column density distributions. We find that kbol is significantly lower at high Lbol with respect to previous estimates by Marconi et al. (2004) and Hopkins et al. (2007). Black hole masses and Eddington ratios are available for 170 Type-1 AGN, while black hole masses for Type-2 AGN are computed for 481 objects using the black hole mass-stellar mass relation and the morphological information. We confirm a trend between kbol and lambda_Edd, with lower hard X-ray bolometric corrections at lower Eddington ratios for both Type-1 and Type-2 AGN. We find that, on average, Eddington ratio increases with redshift for all Types of AGN at any given Mbh, while no clear evolution with redshift is seen at any given Lbol.

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.

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

The analysis of the physical properties of low-redshift Ly$\alpha$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$\alpha$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$\alpha$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$\alpha$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.

We compare the average star formation (SF) activity in X-ray selected AGN hosts with mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are broadly consistent with the locus of their `main sequence'. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.

In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$\alpha$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$\mu$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$\alpha$ emitting galaxies whose rest-frame Ly$\alpha$ equivalent widths (Ly$\alpha$ EW$_{rest-frame}$) are above 20\AA, the typical threshold in narrow-band searches. Ly$\alpha$ emitting galaxies with Ly$\alpha$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$\alpha$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$\alpha$ emission.

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep IRS mid-infrared spectra, to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7<z<2.5 normal star forming galaxies around the main sequence (the redshift-dependent relation of star formation rate and stellar mass). Our deep data confirm from individual far-infrared detections that z~2 star formation rates are overestimated if based on 24 um fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame nuLnu(8) to 8-1000 um infrared luminosity (LIR) tend to lie along lines of constant offset from the main sequence. We explore the relation between SED shape and offset in specific star formation rate (SSFR) from the redshift-dependent main sequence. Main sequence galaxies tend to have a similar nuLnu(8)/LIR regardless of LIR and redshift, up to z~2.5, and nuLnu(8)/LIR decreases with increasing offset above the main sequence in a consistent way at the studied redshifts. We provide a redshift-independent calibration of SED templates in the range of 8--60 um as a function of log(SSFR) offset from the main sequence. Redshift dependency enters only through the evolution of the main sequence with time. Ultra deep IRS spectra match these SED trends well and verify that they are mostly due to a change in ratio of PAH to LIR rather than continua of hidden AGN. Alternatively, we discuss the dependence of nuLnu(8)/LIR on LIR. Same nuLnu(8)/LIR is reached at increasingly higher LIR at higher redshift, with shifts relative to local by 0.5 and 0.8 dex in log(LIR) at redshifts z~1 and z~2. Corresponding SED template calibrations are provided for use if no stellar masses are in hand. For most of those z~2 star forming galaxies that also host an AGN, the mid-infrared is dominated by the star forming component.

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.

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

We analyze the dependence of galaxy structure (size and Sersic index) and mode of star formation (\Sigma_SFR and SFR_IR/SFR_UV) on the position of galaxies in the SFR versus Mass diagram. Our sample comprises roughly 640000 galaxies at z~0.1, 130000 galaxies at z~1, and 36000 galaxies at z~2. Structural measurements for all but the z~0.1 galaxies were based on HST imaging, and SFRs are derived using a Herschel-calibrated ladder of SFR indicators. We find that a correlation between the structure and stellar population of galaxies (i.e., a 'Hubble sequence') is already in place since at least z~2.5. At all epochs, typical star-forming galaxies on the main sequence are well approximated by exponential disks, while the profiles of quiescent galaxies are better described by de Vaucouleurs profiles. In the upper envelope of the main sequence, the relation between the SFR and Sersic index reverses, suggesting a rapid build-up of the central mass concentration in these starbursting outliers. We observe quiescent, moderately and highly star-forming systems to co-exist over an order of magnitude or more in stellar mass. At each mass and redshift, galaxies on the main sequence have the largest size. The rate of size growth correlates with specific SFR, and so does \Sigma_SFR at each redshift. A simple model using an empirically determined SF law and metallicity scaling, in combination with an assumed geometry for dust and stars is able to relate the observed \Sigma_SFR and SFR_IR/SFR_UV, provided a more patchy dust geometry is assumed for high-redshift galaxies.

Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.

One remaining open question regarding the physical properties of Ly$\alpha$ emitters (LAEs) is their dust content and its evolution with redshift. The variety of results is large and with those reported by now is difficult to establish clear relations between dust, other fundamental parameters of galaxies (star-formation rate, metallicity or age) and redshift. In this Letter, we report \emph{Herschel} PACS-100$\mu$m, PACS-160$\mu$m and \emph{Spitzer} MIPS-24$\mu$m detections of a sample of spectroscopically GALEX-selected LAEs at z$\sim$0.3 and $\sim$1.0. Five out of ten and one out of two LAEs are detected in, at least, one PACS band at z$\sim$0.3 and $\sim$1.0, respectively. These measurements have a great importance given that they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. MIPS-24$\mu$m detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain that mid-IR/FIR detected star-forming (SF) LAEs at z$\sim$0.3 have dust content within 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.0, with a median value of A$_{1200\textrm{\AA}}$$\sim$1.1. This range broadens out to 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.5 when considering those LAEs at z$\sim$1.0. Only one SF LAE is undetected both in MIPS-24$\mu$m and PACS, with $A_{1200\AA}$ $\lesssim$0.75. These results seem to be larger than those reported for high-redshift LAEs and, therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic offset in the results.

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

We report on the detection of bright CO(4-3) line emission in two powerful, obscured quasars discovered in the SWIRE survey, SW022513 and SW022550 at z>3.4. We analyze the line strength and profile to determine the gas mass, dynamical mass and the gas dynamics for both galaxies. In SW022513 we may have found the first evidence for a molecular, AGN-driven wind in the early Universe. The line profile in SW022513 is broad (FWHM = 1000 km/s) and blueshifted by -200 km/s relative to systemic (where the systemic velocity is estimated from the narrow components of ionized gas lines, as is commonly done for AGN at low and high redshifts). SW022550 has a more regular, double-peaked profile, which is marginally spatially resolved in our data, consistent with either a merger or an extended disk. The molecular gas masses, 4x10^10 Msun, are large and account for <30% of the stellar mass, making these obscured QSOs as gas rich as other powerful CO emitting galaxies at high redshift, i.e., submillimeter galaxies. Our sources exhibit relatively lower star-formation efficiencies compared to other dusty, powerful starburst galaxies at high redshift. We speculate that this could be a consequence of the AGN perturbing the molecular gas.

We compare multi-wavelength SFR indicators out to z~3 in GOODS-South. Our analysis uniquely combines U-to-8um photometry from FIREWORKS, MIPS 24um and PACS 70, 100, and 160um photometry from the PEP survey, and Ha spectroscopy from the SINS survey. We describe a set of conversions that lead to a continuity across SFR indicators. A luminosity-independent conversion from 24um to total infrared luminosity yields estimates of LIR that are in the median consistent with the LIR derived from PACS photometry, albeit with significant scatter. Dust correction methods perform well at low to intermediate levels of star formation. They fail to recover the total amount of star formation in systems with large SFR_IR/SFR_UV ratios, typically occuring at the highest SFRs (SFR_UV+IR \gtrsim 100 Msun/yr) and redshifts (z \gtrsim 2.5) probed. Finally, we confirm that Ha-based SFRs at 1.5<z<2.6 are consistent with SFR_SED and SFR_UV+IR provided extra attenuation towards HII regions is taken into account (Av,neb = Av,continuum / 0.44). With the cross-calibrated SFR indicators in hand, we perform a consistency check on the star formation histories inferred from SED modeling. We compare the observed SFR-M relations and mass functions at a range of redshifts to equivalents that are computed by evolving lower redshift galaxies backwards in time. We find evidence for underestimated stellar ages when no stringent constraints on formation epoch are applied. We demonstrate how resolved SED modeling, or alternatively deep UV data, may help to overcome this bias. The age bias is most severe for galaxies with young stellar populations, and reduces towards older systems. Finally, our analysis suggests that SFHs typically vary on timescales that are long (at least several 100 Myr) compared to the galaxies' dynamical time.

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

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z<=1. The 50%-light redshifts lie at z=0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).

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

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

We want to study the mid-infrared properties and the starburst and AGN contributions, of 24um sources at z~2, through analysis of mid-infrared spectra combined with millimeter, radio, and infrared photometry. Mid-infrared spectroscopy allows us to recover accurate redshifts. A complete sample of 16 Spitzer-selected sources (ULIRGs) believed to be starbursts at z~2 ("5.8um-peakers") was selected in the (0.5 sq.deg.) J1064+56 SWIRE Lockman Hole field. These sources have S(24um)>0.5mJy, a stellar emission peak redshifted to 5.8um, and r'(Vega)>23. The entire sample was observed with the low resolution units of the Spitzer/IRS infrared spectrograph. These sources have 1.2mm observations with IRAM 30m/MAMBO and very deep 20cm observations from the VLA. Nine of our sources also benefit from 350um observation and detection from CSO/SHARC-II. The entire sample shows good quality IRS spectra dominated by strong PAH features. The main PAH features at 6.2, 7.7, 8.6, and 11.3um have high S/N average luminosities of 2.90, 10.38, 3.62, and 2.29x10^{10}Lsun, respectively. We derived accurate redshifts spanning from 1.75 to 2.28. The average of these redshifts is 2.017. This result confirms that the selection criteria of "5.8um-peakers" associated with a strong detection at 24um are reliable to select sources at z~2. We have analyzed the different correlations between PAH emission and infrared, millimeter, and radio emission. Practically all our sources are strongly dominated by starburst emission. We have also defined two subsamples based on the equivalent width at 7.7um to investigate AGN contributions. Our sample contains strong starbursts and represents a particularly 24um-bright class of SMGs. The very good correlation between PAH and far-IR luminosities is now confirmed in high-z starburst ULIRGs. These sources show a small AGN contribution to the mid-IR, around ~20% in most cases.

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

We present first insights into the far-IR properties for a sample of IRAC and MIPS-24um detected Lyman Break Galaxies (LBGs) at z ~ 3, as derived from observations in the northern field of the Great Observatories Origins Survey (GOODS-N) carried out with the PACS instrument on board the Herschel Space Observatory. Although none of our galaxies are detected by Herschel, we employ a stacking technique to construct, for the first time, the average spectral energy distribution of infrared luminous LBGs from UV to radio wavelengths. We derive a median IR luminosity of L_{IR} = 1.6 x 10^12 Lo, placing the population in the class of ultra luminous infrared galaxies (ULIRGs). Complementing our study with existing multi-wavelength data, we put constraints on the dust temperature of the population and find that for their L_{IR}, MIPS-LBGs are warmer than submm-luminous galaxies while they fall in the locus of the L_{IR}-T_{d} relation of the local ULIRGs. This, along with estimates based on the average SED, explains the marginal detection of LBGs in current sub-mm surveys and suggests that these latter studies introduce a bias towards the detection of colder ULIRGs in the high-z universe, while missing high-z ULIRGS with warmer dust.

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

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

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

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

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z~2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 micron and UV SFRs. We find that SFRs from 24 micron alone are higher by a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24um)>12.2 L_sun. The SFGs and AGNs tend to exhibit the same 24 micron excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions.

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

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

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

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

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

A complete sample of 33 sources believed to be starbursts ("5.8micron-peakers") was selected in the (0.5 sq. deg.) J1046+56 field with selection criteria F_(24micron)>400muJy, the presence of a redshifted stellar emission peak at 5.8um, and r'(Vega)>23. The field, part of the SWIRE Lockman Hole field, benefits from very deep VLA/GMRT 20cm, 50cm, and 90cm radio data (all 33 sources are detected at 50cm), and deep 160micron and 70micron Spitzer data. The 33 sources, with photometric redshifts ~1.5-2.5, were observed at 1.2mm with IRAM-30m/MAMBO to an rms ~0.7-0.8mJy in most cases. Their millimeter, radio, 7-band Spitzer, and near-IR properties were jointly analyzed. The entire sample of 33 sources has an average 1.2mm flux density of 1.56+/-0.22mJy and a median of 1.61mJy, so the majority of the sources can be considered SMGs. Four sources have confirmed 4sigma detections, and nine were tentatively detected at the 3sigma level. Because of its 24micron selection, our sample shows systematically lower F_(1.2mm)/F_(24micron) flux ratios than classical SMGs, probably because of enhanced PAH emission. A median FIR SED was built by stacking images at the positions of 21 sources in the region of deepest Spitzer coverage. Its parameters are T_(dust)=37+/-8K, L_(FIR)=2.5x10^{12}Lo, and SFR=450Mo/yr. The FIR-radio correlation provides another estimate of L_(FIR) for each source, with an average value of 4.1x10^{12}Lo; however, this value may be overestimated because of some AGN contribution. Most of our targets are also luminous star-forming BzK galaxies which constitute a significant fraction of weak SMGs at 1.7<z<2.3.

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 used the large binocular camera (LBC) mounted on the large binocular telescope (LBT) to observe the Lockman Hole in the U, B, and V bands. Our observations cover an area of 925 sq.arcmin. We reached depths of 26.7, 26.3, and 26.3 mag(AB) in the three bands, respectively, in terms of 50% source detection efficiency, making this survey the deepest U-band survey and one of the deepest B and V band surveys with respect to its covered area. We extracted a large number of sources (~89000), detected in all three bands and examined their surface density, comparing it with models of galaxy evolution. We find good agreement with previous claims of a steep faint-end slope of the luminosity functions, caused by late-type and irregular galaxies at z>1.5. A population of dwarf star-forming galaxies at 1.5<z<2.5 is needed to explain the U-band number counts. We also find evidence of strong supernova feedback at high redshift. This survey is complementary to the r, i, and z Lockman Hole survey conducted with the Subaru telescope and provides the essential wavelength coverage to derive photometric redshifts and select different types of sources from the Lockman Hole for further study.

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.

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 report the discovery of two sources at z=3.867 and z=3.427 that exhibit powerful starburst and AGN activities. They benefit from data from radio to X rays from the CFHTLS-D1/SWIRE/XMDS surveys. Follow-up optical and near-infrared spectroscopy, and millimeter IRAM/MAMBO observations are also available. We performed an analysis of their spectral energy distributions to understand the origin of their emission and constrain their luminosities. A comparison with other composite systems at similar redshifts from the literature is also presented. The AGN and starburst bolometric luminosities are ~10^13 Lsun. The AGN emission dominates at X ray, optical, mid-infrared wavelengths, and probably in the radio. The starburst emission dominates in the far-infrared. The estimated star formation rates range from 500 to 3000Msun/yr. The AGN near-infrared and X ray emissions are heavily obscured in both sources with an estimated dust extinction Av>4, and Compton-thick gas column densities. The two sources are the most obscured and most luminous AGNs detected at millimeter wavelengths currently known. The sources presented in this work are heavily obscured QSOs, but their properties are not fully explained by the standard AGN unification model. In one source, the ultraviolet and optical spectra suggest the presence of outflowing gas and shocks, and both sources show emission from hot dust, most likely in the vicinity of the nucleus. Evidence of moderate AGN-driven radio activity is found in both sources. The two sources lie on the local M_BH-M_bulge relation. To remain on this relation, their star formation rate has to decrease. Our results support evolutionary models that invoke radio feedback as star formation quenching mechanism, and suggest that such a mechanism might play a major role also in powerful AGNs.

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.

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.

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

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 present the SWIRE Photometric Redshift Catalogue, 1025119 redshifts of unprecedented reliability and accuracy. Our method is based on fixed galaxy and QSO templates applied to data at 0.36-4.5 mu, and on a set of 4 infrared emission templates fitted to infrared excess data at 3.6-170 mu. The code involves two passes through the data, to try to optimize recognition of AGN dust tori. A few carefully justified priors are used and are the key to supression of outliers. Extinction, A_V, is allowed as a free parameter. We use a set of 5982 spectroscopic redshifts, taken from the literature and from our own spectroscopic surveys, to analyze the performance of our method as a function of the number of photometric bands used in the solution and the reduced chi^2. For 7 photometric bands the rms value of (z_{phot}-z_{spec})/(1+z_{spec}) is 3.5%, and the percentage of catastrophic outliers is ~1%. We discuss the redshift distributions at 3.6 and 24 mu. In individual fields, structure in the redshift distribution corresponds to clusters which can be seen in the spectroscopic redshift distribution. 10% of sources in the SWIRE photometric redshift catalogue have z >2, and 4% have z>3, so this catalogue is a huge resource for high redshift galaxies. A key parameter for understanding the evolutionary status of infrared galaxies is L_{ir}/L_{opt}, which can be interpreted as the specific star-formation rate for starbursts. For dust tori around Type 1 AGN, L_{tor}/L_{opt} is a measure of the torus covering factor and we deduce a mean covering factor of 40%.

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)

We have conducted sensitive (1 sigma<30 uJy) 1.4 GHz radio observations with the Australia Telescope Compact Array of a field largely coincident with infrared observations of the Spitzer Wide-Area Extragalactic Survey. The field is centred on the European Large Area ISO Survey S1 region and has a total area of 3.9 deg. We describe the observations and calibration, source extraction, and cross-matching to infrared sources. Two catalogues are presented; one of the radio components found in the image and one of radio sources with counterparts in the infrared and extracted from the literature. 1366 radio components were grouped into 1276 sources, 1183 of which were matched to infrared sources. We discover 31 radio sources with no infrared counterpart at all, adding to the class of Infrared-Faint Radio Sources.

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.

(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) Among the various issues that remain open in the field of accretion onto black hole X-ray binaries (BHBs) is the way the gas accretes at very low Eddington ratios, in the so-called quiescent regime. While there is general agreement that the X-rays are produced by a population of high-energy electrons near to the BH, the controversy comes about in modeling the contribution from inflowing vs. outflowing particles, and their relative energy budget. Recent Spitzer observations of three quiescent BHBs have shown evidence for excess emission with respect to the tail of the companion star between 8-24 micron. We suggest that synchrotron emission from a partially self-absorbed outflow might be responsible for the observed mid-IR excess, in place of, or in addition to, thermal emission from circumbinary material. If so, then the jet synchrotron luminosity exceeds the measured 2-10 keV luminosity by a factor of a few in these systems. In turn, the mechanical power stored in the jet exceeds the bolometric X-ray luminosity at least by 4 orders of magnitude. We then compile the broadband spectral energy distribution (SED) of A0620-00, the lowest Eddington-ratio stellar mass BH with a known radio counterpart, by means of simultaneous radio, optical and X-ray observations, and the archival Spitzer data. We are able to fit the SED of A0620-00 with a `maximally jet-dominated' model in which the radio through the soft X-rays are dominated by synchrotron emission, while the hard X-rays are dominated by inverse Compton at the jet base. The fitted parameters land in a range of values that is reminiscent of the Galactic Center super-massive BH Sgr A*. Most notably, the inferred ratio of the jet acceleration rate to local cooling rates is two orders of magnitude weaker with respect to higher luminosity, hard state sources.

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

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 present the rest-frame optical and infrared colours of a complete sample of 1114 z<0.3 galaxies from the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) and the Sloan Digital Sky Survey (SDSS). We discuss the optical and infrared colours of our sample and analyse in detail the contribution of dusty star-forming galaxies and AGN to optically selected red sequence galaxies. We propose that the optical (g-r) colour and infrared log(L_{24}/L_{3.6}) colour of galaxies in our sample are determined primarily by a bulge-to-disk ratio. The (g-r) colour is found to be sensitive to the bulge-to-disk ratio for disk-dominated galaxies, whereas the log(L_{24}/L_{3.6}) colour is more sensitive for bulge-dominated systems. We identify ~18% (195 sources) of our sample as having red optical colours and infrared excess. Typically, the infrared luminosities of these galaxies are found to be at the high end of star-forming galaxies with blue optical colours. Using emission line diagnostic diagrams, 78 are found to have an AGN contribution, and 117 are identified as star-forming systems. The red (g-r) colour of the star-forming galaxies could be explained by extinction. However, their high optical luminosities cannot. We conclude that they have a significant bulge component. The number densities of optically red star-forming galaxies are found to correspond to ~13% of the total number density of our sample. In addition, these systems contribute ~13% of the total optical luminosity density, and 28% of the total infrared luminosity density of our SWIRE/SDSS sample. These objects may reduce the need for "dry-mergers".

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.

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 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 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 present a combined galaxy-quasar approach to template-fitting photometric redshift techniques and show the method to be a powerful one. The code (ImpZ) is presented, developed and applied to two spectroscopic redshift catalogues, namely the Isaac Newton Telescope Wide Angle Survey ELAIS N1 and N2 fields and the Chandra Deep Field North. In particular, optical size information is used to improve the redshift determination. The success of the code is shown to be very good with Delta z/(1+z) constrained to within 0.1 for 92 per cent of the galaxies in our sample. The extension of template-fitting to quasars is found to be reasonable with Delta z/(1+z) constrained to within 0.25 for 68 per cent of the quasars in our sample. Various template extensions into the far-UV are also tested.

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

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