We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $\mu$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $\mu$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2''$ to $1.5''$ ($1'' \sim6.9$ kpc), our observation clearly separates, and moderately resolves the two lensed images in both continuum and line emission at $\rm S/N>5$ . Our continuum-based lensing model implies an averaged amplification factor of $\sim5.8$ and reveals that the de-lensed continuum image has the S\'ersic index $\simeq 0.95$ and the S\'ersic radius of $\sim0.18'' (\sim 1.24$ kpc). Furthermore, the reconstructed [N II] velocity field in the source plane is dominated by a rotation component with a maximum velocity of $\sim 300$ km/s at large radii, indicating a dark matter halo mass of $\sim 10^{12}M_{\odot}$. This, together with the reconstructed velocity dispersion field being smooth and modest in value ($<100$ km/s) over much of the outer parts of the galaxy, favours the interpretation of ID 141 being a disk galaxy dynamically supported by rotation. The observed [N II]/CO (7-6) and [N II]/[C II] 158 $\mu$m line luminosity ratios, which are consistent with the corresponding line ratio vs. far-infrared color correlation from local luminous infrared galaxies, imply a de-lensed star formation rate of ($1.8\pm 0.6)\times10^3M_\odot$/yr and provide an independent estimate on the size of the star-forming region $0.7^{+0.3}_{-0.3}$ kpc in radius.

We investigate the relation between the detection of the $11.3\,\mu$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,\mu$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,\mu$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN.

Obscured AGN represent a significant fraction of the entire AGN population, especially at high redshift (~70% at z=3--5). They are often characterized by the presence of large gas and dust reservoirs that are thought to sustain and possibly obscure vigorous star formation processes that make these objects shine at far-IR and sub-mm wavelengths. We exploit ALMA Cycle 4 observations of the continuum (~2.1mm) and high-J CO emission of a sample of six X-ray selected SMGs hosting an obscured AGN at z_spec>2.5 in the 7 Ms CDF-S. We measured the masses and sizes of the dust and molecular gas and we derived the gas density and column density on the basis of a uniform sphere geometry. Finally, we compared the measured column densities with those derived from the Chandra X-ray spectra. We detected both the continuum and line emission for three sources for which we measured both the flux density and size. For the undetected sources, we derived an upper limit on the flux density. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range <0.5 - 2.7 x 10^10 M_sun for {\alpha}_CO=0.8 and up to ~2 x 10^11~M_sun for {\alpha}_CO=6.5, and dust mass <0.9 - 4.9 x 10^8 M_sun) and compact (gas major axis 2.1-3.0 kpc, dust major axis 1.4-2.7 kpc). The column densities associated with the ISM are on the order of 10^(23-24) cm-2, which is comparable with those derived from the X-ray spectra. For the detected sources we also derived dynamical masses in the range 0.8 - 3.7 x 10^10 M_sun. We conclude that the ISM of high redshift galaxies can substantially contribute to nuclear obscuration up to the Compton-thick (>10^24 cm-2) regime. In addition, we found that all the detected sources show a velocity gradient reminding one rotating system, even though two of them show peculiar features in their morphology that can be associated with a chaotic, possibly merging, structure.

AGN-driven outflows are believed to play an important role in regulating the growth of galaxies mostly via negative feedback. However, their effects on their hosts are far from clear, especially for low and moderate luminosity Seyferts. To investigate this issue, we have obtained cold molecular gas observations, traced by the CO(2-1) transition, using the NOEMA interferometer of five nearby (distances between 19 and 58 Mpc) Seyfert galaxies. The resolution of approx. 0.3-0.8 arcsec (approx. 30-100 pc) and field of view of NOEMA allowed us to study the CO(2-1) morphology and kinematics in the nuclear regions (approx. 100 pc) and up to radial distances of approx. 900 pc. We have detected CO(2-1) emission in all five galaxies with disky or circumnuclear ring like morphologies. We derived cold molecular gas masses on nuclear (approx. 100 pc) and circumnuclear (approx. 650 pc) scales in the range from $10^6$ to $10^7$M$_{\odot}$ and from $10^7$ to $10^8$ $M_{\odot}$, respectively. In all of our galaxies the bulk of this gas is rotating in the plane of the galaxy. However, non-circular motions are also present. In NGC 4253, NGC 4388 and NGC 7465, we can ascribe the streaming motions to the presence of a large-scale bar. In Mrk 1066 and NGC 4388, the non-circular motions in the nuclear regions are explained as outflowing material due to the interaction of the AGN wind with molecular gas in the galaxy disk. We conclude that for an unambiguous and precise interpretation of the kinematics of the cold molecular gas we need a detailed knowledge of the host galaxy (i.e., presence of bars, interactions, etc) as well as of the ionized gas kinematics and the ionization cone geometry.

In this paper we present our imaging observations on the CO(7-6) line and its underlying continuum emission of the young submillimeter galaxy LESS 073 at redshift 4.755, using the Atacama Large Millimeter/submillimeter Array (ALMA). At the achieved resolution of $\sim$$1^{\prime\prime}.2\times0^{\prime\prime}.9$ ($8\times6$~kpc$^2$), the CO(7-6) emission is largely unresolved (with a deconvolved size of $1^{\prime\prime}.1(\pm0^{\prime\prime}.5) \times 0^{\prime\prime}.9(\pm0^{\prime\prime}.8)$), and the continuum emission is totally unresolved. The CO(7-6) line emission has an integrated flux of $0.86\pm0.08$~Jy km/s, and a line width of $343\pm40$ km/s. The continuum emission has a flux density of 0.51 mJy. By fitting the observed far-infrared (FIR) spectral energy distribution of LESS 073 with a single-temperature modified blackbody function, we obtained a dust temperature $T_{\rm dust}=57.6\pm3.5$ K, 60-to-100 $\mu$m flux density ratio $f_{60}/f_{100}=0.86\pm0.08$, and total infrared luminosity $L_{\rm IR}=(5.8\pm0.9) \times 10^{12}~L_\odot$. The SED-fit-based $f_{60}/f_{100}$ is consistent with those estimated from various line ratios as advocated by our earlier work, indicating that those proposed line-ratio-based method can be used to practically derive $f_{60}/f_{100}$ for high-$z$ sources. The total molecular gas mass of LESS 073 is $(3.3\pm1.7) \times10^{10}~M_\odot$, and the inferred gas depletion time is about 43 Myr.

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.

We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $\alpha_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $\alpha_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $\alpha_{\mathrm{[CI](1-0)}}$ to $\alpha_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $\alpha_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas.

We present the discovery of HLock01-LAB, a luminous and large Lya nebula at z=3.326. Medium-band imaging and long-slit spectroscopic observations with the Gran Telescopio Canarias reveal extended emission in the Lya 1215\AA, CIV1550\AA, and HeII 1640\AA lines over ~100kpc, and a total luminosity L(Lya)=(6.4+/-0.1)x10^44 erg s^-1. HLock01-LAB presents an elongated morphology aligned with two faint radio sources contained within the central ~8kpc of the nebula. The radio structures are consistent to be faint radio jets or lobes of a central galaxy, whose spectrum shows nebular emission characteristic of a type-II active galactic nucleus (AGN). The continuum emission of the AGN at short wavelengths is, however, likely dominated by stellar emission of the host galaxy, for which we derive a stellar mass M* = 2.3x10^11 Msun. The detection of extended emission in CIV and CIII] indicates that the gas within the nebula is not primordial. Feedback may have enriched the halo at at least 50 kpc from the nuclear region. Using rest-frame UV emission-line diagnostics, we find that the gas in the nebula is likely heated by the AGN. Nevertheless, at the center of the nebula we find extreme emission line ratios of Lya/CIV~60 and Lya/HeII~80, one of the highest values measured to date, and well above the standard values of photoionization models (Lya/HeII~30 for case B photoionization). Our data suggest that jet-induced shocks are likely responsible for the increase of the electron temperature and, thus, the observed Lya enhancement in the center of the nebula. This scenario is further supported by the presence of radio structures and perturbed kinematics in this region. The large Lya luminosity in HLock01-LAB is likely due to a combination of AGN photoionization and jet-induced shocks, highlighting the diversity of sources of energy powering Lya nebulae. [abridged]

We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$\mu$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged)

The extreme infrared (IR) luminosity of local luminous and ultra-luminous IR galaxies (U/LIRGs; 11 < log LIR /Lsun < 12 and log LIR /Lsun > 12, respectively) is mainly powered by star-formation processes triggered by mergers or interactions. While U/LIRGs are rare locally, at z > 1, they become more common, they dominate the star-formation rate (SFR) density, and a fraction of them are found to be normal disk galaxies. Therefore, there must be an evolution of the mechanism triggering these intense starbursts with redshift. To investigate this evolution, we present new optical SWIFT integral field spectroscopic H{\alpha}+[NII] observations of a sample of 9 intermediate-z (0.2 < z < 0.4) U/LIRG systems selected from Herschel 250{\mu}m observations. The main results are the following: (a) the ratios between the velocity dispersion and the rotation curve amplitude indicate that 10-25% (1-2 out of 8) might be compatible with being isolated disks while the remaining objects are interacting/merging systems; (b) the ratio between un-obscured and obscured SFR traced by H{\alpha} and LIR, respectively, is similar in both local and these intermediate-z U/LIRGs; and (c) the ratio between 250{\mu}m and the total IR luminosities of these intermediate-z U/LIRGs is higher than that of local U/LIRGs with the same LIR . This indicates a reduced dust temperature in these intermediate-z U/LIRGs. This, together with their already measured enhanced molecular gas content, suggests that the interstellar medium conditions are different in our sample of intermediate-z galaxies when compared to local U/LIRGs.

Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.

[Abridged] We combine new CO(1-0) line observations of 24 intermediate redshift galaxies (0.03 < z < 0.28) along with literature data of galaxies at 0<z<4 to explore scaling relations between the dust and gas content using PAH 6.2 $\mu$m ($L_{6.2}$), CO ($L'_{\rm CO}$), and infrared ($L_{\rm IR}$) luminosities for a wide range of redshifts and physical environments. Our analysis confirms the existence of a universal $L_{6.2}-L_{\rm CO}$ correlation followed by normal star-forming galaxies (SFGs) and starbursts (SBs) at all redshifts. This relation is also followed by local ULIRGs that appear as outliers in the $L_{6.2}-L_{\rm IR}$ and $L_{\rm IR}-L'_{\rm CO}$ relations from the sequence defined by normal SFGs. The emerging tight ($\sigma \approx 0.26$ dex) and linear ($\alpha = 1.03$) relation between $L_{6.2}$ and $L_{\rm CO}$ indicates a $L_{6.2}$ to molecular gas ($M_{\rm H_2}$) conversion factor of $\alpha_{6.2} = M_{\rm H2}/L_{6.2} = (2.7\pm1.3) \times \alpha_{\rm CO}$, where $\alpha_{\rm CO}$ is the $L'_{\rm CO}$ to $M_{\rm H_2}$ conversion factor. We also find that on galaxy integrated scales, PAH emission is better correlated with cold rather than with warm dust emission, suggesting that PAHs are associated with the diffuse cold dust, which is another proxy for $M_{\rm H_2}$. Focusing on normal SFGs among our sample, we employ the dust continuum emission to derive $M_{\rm H_2}$ estimates and find a constant $M_{\rm H_2}/L_{6.2}$ ratio of $\alpha_{6.2} = 12.3 \ M_{\rm H_2}/{\rm L}_{\odot}$ ($\sigma\approx 0.3$ dex). We propose that the presented $L_{6.2}-L'_{\rm CO}$ and $L_{6.2}-M_{\rm H_2}$ relations will serve as useful tools for the determination of the physical properties of high-$z$ SFGs, for which PAH emission will be routinely detected by the James Webb Space Telescope.

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

We present the results from our Atacama Large Millimeter/submillimeter Array (ALMA) imaging observations of the CO (7-6), [CI] 370 um (hereafter [CI]) and [NII] 205 um (hereafter [NI]I) lines and their underlying continuum emission of BRI 1335-0417, an infrared bright quasar at z = 4.407. At the achieved resolutions of 1.1" to 1.2" (or 7.5 to 8.2 kpc), the continuum at 205 and 372 um (rest-frame), the CO (7-6), and the [CI] emissions are at best barely resolved whereas the [NII] emission is well resolved with an ALMA beam de-convolved major axis of 1.3" (+/- 0.3") or 9 (+/-2) kpc. As a warm dense gas tracer, the CO (7-6) emission shows a more compact spatial distribution and a significantly higher peak velocity dispersion than the other two lines that probe lower density gas, a picture favoring a merger-triggered star formation (SF) scenario over an orderly rotating SF disk. The CO (7-6) data also indicate a possible QSO-driven gas outflow that reaches a maximum line-of-sight velocity of 500 to 600 km/s. The far-infrared (FIR) dust temperature (T_dust) of 41.5 K from a gray-body fit to the continuum agrees well with the average T_dust inferred from various line luminosity ratios. The resulting L_CO(7-6)/L_FIR luminosity ratio is consistent with that of local luminous infrared galaxies powered predominantly by SF. The CO(7-6) luminosity-inferred SF rate is 5.1 (+/-1.5) x 10^3 M_solar/yr . The system has an effective star-forming region of 1.7 (+1.7/-0.8) kpc in diameter and a molecular gas reservoir of ~5 x 10^{11} M_solar.

We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies.

We present ALMA Band 6 $^{12}$CO(2--1) line and rest-frame 232GHz continuum observations of the nearby Compton-thick Seyfert galaxy NGC5643 with angular resolutions 0.11-0.26arcsec (9-21pc). The CO(2--1) integrated line map reveals emission from the nuclear and circumnuclear region with a two-arm nuclear spiral extending 10arcsec on each side. The circumnuclear CO(2--1) kinematics can be fitted with a rotating disk, although there are regions with large residual velocities and/or velocity dispersions. The CO(2--1) line profiles of these regions show two different velocity components. One is ascribed to the circular component and the other to the interaction of the AGN outflow, as traced by the [O III]5007AA emission, with molecular gas in the disk a few hundred parsecs from the AGN. On nuclear scales, we detected an inclined CO(2--1) disk (diameter 26 pc, FWHM) oriented almost in a north-south direction. The CO(2--1) nuclear kinematics can be fitted with a rotating disk which appears to be tilted with respect to the large scale disk. There are strong non-circular motions in the central 0.2-0.3 arcsec with velocities of up to 110km/s. In the absence of a nuclear bar, these motions could be explained as radial outflows in the nuclear disk. We estimate a total molecular gas mass for the nuclear disk of $M({\rm H}_2)=1.1\times 10^7\,M_\odot$ and an H$_2$ column density toward the location of the AGN of $N({\rm H}_2)\sim 5 \times 10^{23}\,{\rm cm}^{-2}$, for a standard CO-to-H$_2$ conversion factor. We interpret this nuclear molecular gas disk as the obscuring torus of NGC5643 as well as the collimating structure of the ionization cone.

We present new CO(2-1) observations of 3 low-z (~350 Mpc) ULIRG systems (6 nuclei) observed with ALMA at high-spatial resolution (~500 pc). We detect massive cold molecular gas outflows in 5 out of 6 nuclei (0.3-5)x10^8 Msun. These outflows are spatially resolved with deprojected radii of 0.25-1 kpc although high-velocity molecular gas is detected up to ~0.5-1.8 kpc (1-6 kpc deprojected). The mass outflow rates are 12-400 Msun/yr and the inclination corrected average velocity of the outflowing gas 350-550 km/s (v_max = 500-900 km/s). The origin of these outflows can be explained by the nuclear starbursts although the contribution of an obscured AGN can not be completely ruled out. The position angle (PA) of the outflowing gas along the kinematic minor axis of the nuclear molecular disk suggests that the outflow axis is perpendicular to the disk for three of these outflows. Only in one case, the outflow PA is clearly not along the kinematic minor axis. The outflow depletion times are 15-80 Myr which are slightly shorter than the star-formation (SF) depletion times (30-80 Myr). However, we estimate that only 15-30% of the outflowing gas will escape the gravitational potential of the nucleus. The majority of the outflowing gas will return to the disk after 5-10 Myr and become available to form new stars. Therefore, these outflows will not likely quench the nuclear starbursts. These outflows would be consistent with being driven by radiation pressure (momentum-driven) only if the coupling between radiation and dust increases with increasing SF rates. This can be achieved if the dust optical depth is higher in objects with higher SF. The relatively small sizes (<1 kpc) and dynamical times (<3 Myr) of the cold molecular outflows suggests that molecular gas cannot survive longer in the outflow environment or that it cannot form efficiently beyond these distances or times. (Abridged)

The submillimeter galaxy (SMG) HERMES J105751.1+573027 (hereafter HLock01) at z = 2.9574 +/- 0.0001 is one of the brightest gravitationally lensed sources discovered in the Herschel Multi-tiered Extragalactic Survey. Apart from the high flux densities in the far-infrared, it is also extremely bright in the rest-frame ultraviolet (UV), with a total apparent magnitude m_UV = 19.7 mag. We report here deep spectroscopic observations with the Gran Telescopio Canarias of the optically bright lensed images of HLock01. Our results suggest that HLock01 is a merger system composed of the Herschel-selected SMG and an optically bright Lyman break-like galaxy (LBG), separated by only 3.3 kpc in projection. While the SMG appears very massive (M* = 5x10^11 Msun), with a highly extinguished stellar component (A_V = 4.3), the LBG is a young, lower-mass (M* = 1x10^10 Msun), but still luminous (10xL*_UV) satellite galaxy. Detailed analysis of the high signal-to-noise (S/N) rest-frame UV spectrum of the LBG shows complex kinematics of the gas, exhibiting both blueshifted and redshifted absorption components. While the blueshifted component is associated with strong galactic outflows from the massive stars in the LBG, as is common in most star-forming galaxies, the redshifted component may be associated with gas inflow seen along a favorable sightline to the LBG. We also find evidence of an extended gas reservoir around HLock01 at an impact parameter of 110 kpc, through the detection of C II 1334$\AA$ absorption in the red wing of a bright Ly-alpha emitter at z = 3.327. The data presented here highlight the power of gravitational lensing in high S/N studies to probe deeply into the physics of high-z star forming galaxies.

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60${\mu}$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500${\mu}$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $\beta$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $\beta$ distribution (median $\beta$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850${\mu}$m. These latter sources generally have cooler dust and/or redder 100-to-250 ${\mu}$m colours than our 60${\mu}$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population.

In an earlier paper we modeled the far-infrared emission from a star-forming galaxy using the photoionisation code CLOUDY and presented metallicity sensitive diagnostics based on far-infrared fine structure line ratios. Here, we focus on the applicability of the [OIII]88/[NII]122 microns line ratio as a gas phase metallicity indicator in high redshift submillimetre luminous galaxies. The [OIII]88/[NII]122 microns ratio is strongly dependent on the ionization parameter (which is related to the total number of ionizing photons) as well as the gas electron density. We demonstrate how the ratio of 88/$122 continuum flux measurements can provide a reasonable estimate of the ionization parameter while the availability of the [NII]205 microns line can constrain the electron density. Using the [OIII]88/[NII]122 microns line ratios from a sample of nearby normal and star-forming galaxies we measure their gas phase metallicities and find that their mass metallicity relation is consistent with the one derived using optical emission lines. Using new, previously unpublished, Herschel spectroscopic observations of key far-infrared fine structure lines of the z~3 galaxy HLSW-01 and additional published measurements of far-infrared fine structure lines of high-z submillimetre luminous galaxies we derive gas phase metallicities using their [OIII]88/[NII]122 microns line ratio. We find that the metallicities of these z~3 submm luminous galaxies are consistent with solar metallicities and that they appear to follow the mass-metallicity relation expected for z~3 systems.

We present stacked average far-infrared spectra of a sample of 197 dusty, star-forming galaxies (DSFGs) at $0.005 < z < 4$ using close to 90% of the SPIRE Fourier Transform Spectrometer (FTS) extragalactic data archive from the Herschel Space Observatory based on 3.5 years of science operations. These spectra explore an observed-frame $\rm 447\,GHz-1568\,GHz$ ($\rm 191\,\mu m-671\,\mu m$) frequency (wavelength) range allowing us to observe the main atomic and molecular lines emitted by gas in the interstellar medium. The sample is sub-divided into five redshift bins at $0.005 < z < 0.05$, $0.05 < z < 0.2$, $0.2 < z < 0.5$, $0.8 < z <2$, and $2 < z < 4$. To study the dependence of observed spectral lines on total infrared luminosity, the sources in a subset of the redshift bins are stacked in luminosity bins. These stacked spectra are used to determine the average properties of the interstellar medium and dense molecular gas properties of DSFGs, in particular, the fine-structure line ([CII] 158 $\mu$m and [OI] 63 $\mu$m) luminosity ratios, and the line to far-IR luminosity ratios are used to model the gas density and radiation field strength in the photodissociation regions (PDRs). For the low-redshift sample, we additionally present the average spectral line energy distributions (SLED) of CO and $\rm{H_2O}$ rotational transitions and also consider PDR conditions based on a combination of [CI] 370 $\mu$m and 609 $\mu$m and $\rm CO (7-6)$ lines. For the high-z ($0.8 < z < 4$) sample PDR models suggest a molecular gas distribution in the presence of a radiation field that is at least a factor of 10$^3$ larger than the Milky-Way and with a neutral gas density of roughly 10$^3$ to 10$^5$ cm$^{-3}$. The corresponding PDR models for the low-z sample suggest a UV radiation field and gas density comparable to those at high-z.

We report the detection of ADFS-27, a dusty, starbursting major merger at a redshift of z=5.655, using the Atacama Large Millimeter/submillimeter Array (ALMA). ADFS-27 was selected from Herschel/SPIRE and APEX/LABOCA data as an extremely red "870 micron riser" (i.e., S_250<S_350<S_500<S_870), demonstrating the utility of this technique to identify some of the highest-redshift dusty galaxies. A scan of the 3mm atmospheric window with ALMA yields detections of CO(5-4) and CO(6-5) emission, and a tentative detection of H2O(211-202) emission, which provides an unambiguous redshift measurement. The strength of the CO lines implies a large molecular gas reservoir with a mass of M_gas=2.5x10^11(alpha_CO/0.8)(0.39/r_51) Msun, sufficient to maintain its ~2400 Msun/yr starburst for at least ~100 Myr. The 870 micron dust continuum emission is resolved into two components, 1.8 and 2.1 kpc in diameter, separated by 9.0 kpc, with comparable dust luminosities, suggesting an ongoing major merger. The infrared luminosity of L_IR~=2.4x10^13Lsun implies that this system represents a binary hyper-luminous infrared galaxy, the most distant of its kind presently known. This also implies star formation rate surface densities of Sigma_SFR=730 and 750Msun/yr/kpc2, consistent with a binary "maximum starburst". The discovery of this rare system is consistent with a significantly higher space density than previously thought for the most luminous dusty starbursts within the first billion years of cosmic time, easing tensions regarding the space densities of z~6 quasars and massive quiescent galaxies at z>~3.

The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.

We present the results from Atacama Large Millimeter/submillimeter Array (ALMA) imaging in the [NII] 205 micron fine-structure line (hereafter [NII]) and the underlying continuum of BRI 1202-0725, an interacting galaxy system at $z =$ 4.7, consisting of an optical QSO, a sub-millimeter galaxy (SMG) and two Lyman-$\alpha$ emitters (LAEs), all within $\sim$25 kpc of the QSO. We detect the QSO and SMG in both [NII] and continuum. At the $\sim$$1"$ (or 6.6 kpc) resolution, both QSO and SMG are resolved in [NII], with the de-convolved major axes of $\sim$9 and $\sim$14 kpc, respectively. In contrast, their continuum emissions are much more compact and unresolved even at an enhanced resolution of $\sim$$0.7"$. The ratio of the [NII] flux to the existing CO (7$-$6) flux is used to constrain the dust temperature ($T_{\rm dust}$) for a more accurate determination of the FIR luminosity $L_{\rm FIR}$. Our best estimated $T_{\rm dust}$ equals $43 (\pm 2)$ K for both galaxies (assuming an emissivity index $\beta = 1.8$). The resulting $L_{\rm CO(7-6)}/L_{\rm FIR}$ ratios are statistically consistent with that of local luminous infrared galaxies, confirming that $L_{\rm CO(7-6)}$ traces the star formation (SF) rate (SFR) in these galaxies. We estimate that the on-going SF of the QSO (SMG) has a SFR of 5.1 $(6.9) \times 10^3 M_{\odot}$ yr$^{-1}$ ($\pm$ 30%) assuming Chabrier initial mass function, takes place within a diameter (at half maximum) of 1.3 (1.5) kpc, and shall consume the existing 5 $(5) \times 10^{11} M_{\odot}$ of molecular gas in 10 $(7) \times 10^7$ years.

The abundance of metals in galaxies is a key parameter which permits to distinguish between different galaxy formation and evolution models. Most of the metallicity determinations are based on optical line ratios. However, the optical spectral range is subject to dust extinction and, for high-z objects (z > 3), some of the lines used in optical metallicity diagnostics are shifted to wavelengths not accessible to ground based observatories. For this reason, we explore metallicity diagnostics using far-infrared (IR) line ratios which can provide a suitable alternative in such situations. To investigate these far-IR line ratios, we modeled the emission of a starburst with the photoionization code CLOUDY. The most sensitive far-IR ratios to measure metallicities are the [OIII]52$\mu$m and 88$\mu$m to [NIII]57$\mu$m ratios. We show that this ratio produces robust metallicities in the presence of an AGN and is insensitive to changes in the age of the ionizing stellar. Another metallicity sensitive ratio is the [OIII]88$\mu$m/[NII]122$\mu$m ratio, although it depends on the ionization parameter. We propose various mid- and far-IR line ratios to break this dependency. Finally, we apply these far-IR diagnostics to a sample of 19 local ultraluminous IR galaxies (ULIRGs) observed with Herschel and Spitzer. We find that the gas-phase metallicity in these local ULIRGs is in the range 0.7 < Z_gas/Z_sun < 1.5, which corresponds to 8.5 < 12 + log (O/H) < 8.9. The inferred metallicities agree well with previous estimates for local ULIRGs and this confirms that they lie below the local mass-metallicity relation.

We present millimetre dust emission measurements of two Lyman Break Galaxies at z~3 and construct for the first time fully sampled infrared spectral energy distributions (SEDs), from mid-IR to the Rayleigh-Jeans tail, of individually detected, unlensed, UV-selected, main sequence (MS) galaxies at $z=3$. The SED modelling of the two sources confirms previous findings, based on stacked ensembles, of an increasing mean radiation field <U> with redshift, consistent with a rapidly decreasing gas metallicity in z > 2 galaxies. Complementing our study with CO[3-2] emission line observations, we measure the molecular gas mass (M_H2) reservoir of the systems using three independent approaches: 1) CO line observations, 2) the dust to gas mass ratio vs metallicity relation and 3) a single band, dust emission flux on the Rayleigh-Jeans side of the SED. All techniques return consistent M_H2 estimates within a factor of ~2 or less, yielding gas depletion time-scales (tau_dep ~ 0.35 Gyrs) and gas-to-stellar mass ratios (M_H2/M* ~ 0.5-1) for our z~3 massive MS galaxies. The overall properties of our galaxies are consistent with trends and relations established at lower redshifts, extending the apparent uniformity of star-forming galaxies over the last 11.5 billion years.

We present the first detection of the ortho-H2O 4_23-3_30 transition at 448 GHz in space. We observed this transition in the local (z = 0.010) luminous infrared (IR) galaxy ESO 320-G030 (IRAS F11506-3851) using the Atacama Large Millimeter/submillimeter Array (ALMA). The water 4_23-3_30 emission, which originates in the highly obscured nucleus of this galaxy, is spatially resolved over a region of ~65 pc in diameter and shows a regular rotation pattern compatible with the global molecular and ionized gas kinematics. The line profile is symmetric and well fitted by a Gaussian with an integrated flux of 37.0 +- 0.7 Jy km s-1 . Models predict this water transition as a potential collisionally excited maser transition. On the contrary, in this galaxy, we find that the 4_23-3_30 emission is primarily excited by the intense far-IR radiation field present in its nucleus. According to our modeling, this transition is a probe of deeply buried galaxy nuclei thanks to the high dust optical depths (tau_100{\mu}m > 1, N_H > 1e24 cm-2) required to efficiently excite it.

The European Far-Infrared (FIR) Space Roadmap focuses on fundamental, yet still unresolved, astrophysical questions that can only be answered through a far-infrared space mission and gives an overview of the technology required to answer them. The document discusses topics ranging from Solar System and Planet Formation, Our Galaxy and nearby Galaxies and Distant Galaxies and Galaxy Evolution. The FIR Roadmap was open to comments from the wider astronomical community following a presentation during EWASS 2016.

Submillimeter galaxies (SMGs) at $z\gtrsim1$ are luminous in the far-infrared and have star-formation rates, SFR, of hundreds to thousands of solar masses per year. However, it is unclear whether they are true analogs of local ULIRGs or whether the mode of their star formation is more similar to that in local disk galaxies. We target these questions by using Herschel-PACS to examine the conditions in the interstellar medium (ISM) in far-infrared luminous SMGs at z~1-4. We present 70-160 micron photometry and spectroscopy of the [OIV]26 micron, [FeII]26 micron, [SIII]33 micron, [SiII]34 micron, [OIII]52 micron, [NIII]57 micron, and [OI]63 micron fine-structure lines and the S(0) and S(1) hydrogen rotational lines in 13 lensed SMGs identified by their brightness in early Herschel data. Most of the 13 targets are not individually spectroscopically detected and we instead focus on stacking these spectra with observations of an additional 32 SMGs from the \herschel\ archive -- representing a complete compilation of PACS spectroscopy of SMGs. We detect [OI]63 micron, [SiII]34 micron, and [NIII]57 micron at >3sigma in the stacked spectra, determining that the average strengths of these lines relative to the far-IR continuum are $(0.36\pm0.12)\times10^{-3}$, $(0.84\pm0.17)\times10^{-3}$, and $(0.27\pm0.10)\times10^{-3}$, respectively. Using the [OIII]52/[NIII]57 emission line ratio we show that SMGs have average gas-phase metallicities $\gtrsim Z_{\rm sun}$. By using PDR modelling and combining the new spectral measurements with integrated far-infrared fluxes and existing [CII]158 micron data we show that SMGs have average gas densities, n, of $\sim10^{1-3}{\rm cm^{-3}}$ and FUV field strengths, $G_0\sim10^{2.2-4.5}$ (in Habing units: $1.6\times10^{-3}{\rm erg~cm^{-2}~s^{-1}}$), consistent with both local ULIRGs and lower luminosity star-forming galaxies.

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671\mu$m with photometry in the SPIRE bands at 250$\mu$m, 350$\mu$m and 500$\mu$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $\mu$m, [CI] 609 $\mu$m, and [NII] 205 $\mu$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

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

The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.

We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.

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

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst.

The Far-InfraRed Spectroscopic Explorer (FIRSPEX) is a candidate mission in response to a bi-lateral Small-mission call issued by the European Space Agency (ESA) and the Chinese Academy of Sciences (CAS). FIRSPEX is a small satellite (~1m telescope) operating from Low Earth Orbit (LEO). It consists of a number of heterodyne detection bands targeting key molecular and atomic transitions in the terahertz (THz) and Supra-Terahertz (>1 THz) frequency range. The FIRSPEX bands are: [CII] 158 microns (1.9 THz), [NII] 205 microns (1.46 THz), [CI] 370 microns (0.89 THz), CO(6-5) 433 microns (0.69 THz). The primary goal of FIRSPEX is to perform an unbiased all sky spectroscopic survey in four far-infrared lines delivering the first 3D-maps (high spectral resolution) of the Galaxy. The spectroscopic surveys will build on the heritage of Herschel and complement the broad-band all-sky surveys carried out by the IRAS and AKARI observatories. In addition FIRSPEX will enable targeted observations of nearby and distant galaxies allowing for an in-depth study of the ISM components.

The Herschel Multi-tiered Extragalactic Survey (HerMES) has identified large numbers of dusty star-forming galaxies (DSFGs) over a wide range in redshift. A detailed understanding of these DSFGs is hampered by the limited spatial resolution of Herschel. We present 870um 0.45" resolution imaging from the Atacama Large Millimeter/submillimeter Array (ALMA) of 29 HerMES DSFGs with far-infrared (FIR) flux densities in between the brightest of sources found by Herschel and fainter DSFGs found in ground-based sub-millimeter (sub-mm) surveys. We identify 62 sources down to the 5-sigma point-source sensitivity limit in our ALMA sample (sigma~0.2mJy), of which 6 are strongly lensed (showing multiple images) and 36 experience significant amplification (mu>1.1). To characterize the properties of the ALMA sources, we introduce and make use of uvmcmcfit, a publicly available Markov chain Monte Carlo analysis tool for interferometric observations of lensed galaxies. Our lens models tentatively favor intrinsic number counts for DSFGs with a steep fall off above 8mJy at 880um. Nearly 70% of the Herschel sources comprise multiple ALMA counterparts, consistent with previous research indicating that the multiplicity rate is high in bright sub-mm sources. Our ALMA sources are located significantly closer to each other than expected based on results from theoretical models as well as fainter DSFGs identified in the LABOCA ECDFS Submillimeter Survey. The high multiplicity rate and low projected separations argue in favor of interactions and mergers driving the prodigious emission from the brightest DSFGs as well as the sharp downturn above S_880=8mJy.

[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 present Herschel far-IR photometry and spectroscopy as well as ground based CO observations of an intermediate redshift (0.21 < z < 0.88) sample of Herschel-selected (ultra)-luminous infrared galaxies (L_IR > 10^11.5L_sun). With these measurements we trace the dust continuum, far-IR atomic line emission, in particular [CII]\,157.7microns, as well as the molecular gas of z~0.3 (U)LIRGs and perform a detailed investigation of the interstellar medium of the population. We find that the majority of Herschel-selected intermediate redshift (U)LIRGs have L_CII/L_FIR ratios that are a factor of about 10 higher than that of local ULIRGs and comparable to that of local normal and high-$z$ star forming galaxies. Using our sample to bridge local and high-z [CII] observations, we find that the majority of galaxies at all redshifts and all luminosities follow a L_CII-L_FIR relation with a slope of unity, from which local ULIRGs and high-z AGN dominated sources are clear outliers. We also confirm that the strong anti-correlation between the L_CII/L_FIR ratio and the far-IR color L_60/L_100 observed in the local Universe holds over a broad range of redshifts and luminosities, in the sense that warmer sources exhibit lower L_CII/L_FIR at any epoch. Intermediate redshift ULIRGs are also characterised by large molecular gas reservoirs and by lower star formation efficiencies compared to that of local ULIRGs. The high L_CII/L_FIR ratios, the moderate star formation efficiencies (L_LIR/L_CO or L_IR/M_gas) and the relatively low dust temperatures of our sample (which are also common characteristics of high-z star forming galaxies with ULIRG-like luminosities) indicate that the evolution of the physical properties of (U)LIRGs between the present day and z > 1 is already significant by z ~ 0.3.

We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3.

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies at $z>$4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their Spectral Energy Distributions (SED) and a composite model to fit the SEDs. The model includes a stellar emission component at $\lambda_{\rm rest} <$ 3.5$ \mu$m; a hot dust component peaking at $\lambda_{rest} \sim$ 5$\,\mu$m; and cold dust component which becomes significant for $\lambda_{\rm rest} >$ 50$\,\mu$m. Six objects in the sample are detected at 250 and 350$ \mu$m. The dust temperatures for the sources in this sample are in the range of 40$-$80 K, and their $L_{\rm FIR}$ $\sim$ 10$^{13}$ L$_{\odot}$ qualifies them as Hyper$-$Luminous Infrared Galaxies (HyperLIRGs). The mean FIR-radio index for this sample is around $< q > = 2.2$ indicating no radio excess in their radio emission. Most sources in the sample have 24$ \mu$m detections corresponding to a rest-frame 4.5$ \mu$m luminosity of Log$_{10}$(L$_{4.5}$ / L$_{\odot}$) = 11 $\sim$ 11.5. Their L$_{\rm 4.5}$/$L_{\rm FIR}$ ratios are very similar to those of starburst dominated submillimeter galaxies at $z \sim$ 2. The $L_{\rm CO}-L_{\rm FIR}$ relation for this sample is consistent with that determined for local ULIRGs and SMGs at $z \sim$ 2. We conclude that submillimeter galaxies at $z >$ 4 are hotter and more luminous in the FIR, but otherwise very similar to those at $z \sim$ 2. None of these sources show any sign of the strong QSO phase being triggered.

We report the first results from a spectroscopic survey of the [CII] 158um line from a sample of intermediate redshift (0.2<z<0.8) (ultra)-luminous infrared galaxies, (U)LIRGs (LIR>10^11.5 Lsun), using the SPIRE-Fourier Transform Spectrometer (FTS) on board the Herschel Space Observatory. This is the first survey of [CII] emission, an important tracer of star-formation, at a redshift range where the star-formation rate density of the Universe increases rapidly. We detect strong [CII] 158um line emission from over 80% of the sample. We find that the [CII] line is luminous, in the range (0.8-4)x10^(-3) of the far-infrared continuum luminosity of our sources, and appears to arise from photodissociation regions on the surface of molecular clouds. The L[CII]/LIR ratio in our intermediate redshift (U)LIRGs is on average ~10 times larger than that of local ULIRGs. Furthermore, we find that the L[CII]/LIR and L[CII]/LCO(1-0) ratios in our sample are similar to those of local normal galaxies and high-z star-forming galaxies. ULIRGs at z~0.5 show many similarities to the properties of local normal and high-z star forming galaxies. Our findings strongly suggest that rapid evolution in the properties of the star forming regions of luminous infrared galaxies is likely to have occurred in the last 5 billion years.

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

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

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

We study the link between observed ultraviolet luminosity, stellar mass, and dust attenuation within rest-frame UV-selected samples at z~ 4, 3, and 1.5. We measure by stacking at 250, 350, and 500 um in the Herschel/SPIRE images from the HerMES program the average infrared luminosity as a function of stellar mass and UV luminosity. We find that dust attenuation is mostly correlated with stellar mass. There is also a secondary dependence with UV luminosity: at a given UV luminosity, dust attenuation increases with stellar mass, while at a given stellar mass it decreases with UV luminosity. We provide new empirical recipes to correct for dust attenuation given the observed UV luminosity and the stellar mass. Our results also enable us to put new constraints on the average relation between star formation rate and stellar mass at z~ 4, 3, and 1.5. The star formation rate-stellar mass relations are well described by power laws (SFR~ M^0.7), with the amplitudes being similar at z~4 and z~3, and decreasing by a factor of 4 at z~1.5 at a given stellar mass. We further investigate the evolution with redshift of the specific star formation rate. Our results are in the upper range of previous measurements, in particular at z~3, and are consistent with a plateau at 3<z<4. Current model predictions (either analytic, semi-analytic or hydrodynamic) are inconsistent with these values, as they yield lower predictions than the observations in the redshift range we explore. We use these results to discuss the star formation histories of galaxies in the framework of the Main Sequence of star-forming galaxies. Our results suggest that galaxies at high redshift (2.5<z<4) stay around 1 Gyr on the Main Sequence. With decreasing redshift, this time increases such that z=1 Main Sequence galaxies with 10^8<M_*/Msun<10^10 stay on the Main Sequence until z=0.

We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter.

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

We present $^{12}$CO \textit{J=3$\rightarrow$2} maps of NGC 2976 and NGC 3351 obtained with the James Clerk Maxwell Telescope (JCMT), both early targets of the JCMT Nearby Galaxy Legacy Survey (NGLS). We combine the present observations with $^{12}$CO \textit{J=1$\rightarrow$0} data and find that the computed $^{12}$CO \textit{J=3$\rightarrow$2} to $^{12}$CO \textit{J=1$\rightarrow$0} line ratio ($R_{31}$) agrees with values measured in other NGLS field galaxies. We compute the M$_{\rm{H_2}}$ value and find that it is robust against the value of $R_{31}$ used. Using \textsc{Hi} data from the The \textsc{Hi} Nearby Galaxy Survey (THINGS) survey, we find a tight correlation between surface density of H$_2$ and star formation rate density for NGC 3351 when $^{12}$CO \textit{J=3$\rightarrow$2} data are used. Finally, we compare the $^{12}$CO \textit{J=3$\rightarrow$2} intensity with the PAH 8 $\mu$m surface brightness and find a good correlation in the high surface brightness regions. We extend this study to include all 25 \textit{Spitzer} Infrared Nearby Galaxies Survey (SINGS) galaxies within the NGLS sample and find a tight correlation at large spatial scales. We suggest that both PAH 8 $\mu$m and $^{12}$CO \textit{J=3$\rightarrow$2} are likely to originate in regions of active star formation.

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

Dust-obscured galaxies (DOGs) are a UV-faint, IR-bright galaxy population that reside at z~2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far-IR observations of a complete sample of DOGs in the 2 deg2 of COSMOS. The 3077 DOGs have <z>=1.9+/-0.3 and are selected from 24um and r+ observations using a color cut of r+ - [24]>=7.5 (AB mag) and S24>=100uJy. Based on the near-IR SEDs, 47% are star-formation dominated and 10% are AGN-dominated. We use SPIRE far-IR photometry from HerMES to calculate the IR luminosity and characteristic dust temperature for the 1572 (51%) DOGs that are detected at 250um (>=3{\sigma}). For the remaining 1505 (49%) that are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and undetected DOGs have average IR luminosities of (2.8+/-0.4) x 1012 LSun and (0.77+-0.08)x10^12LSun, and dust temperatures of 34+/-7 K and 37+/-3 K, respectively. The IR luminosity function of DOGs with S24>=100uJy is calculated, using far-IR observations and stacking. DOGs contribute 10-30% to the total star formation rate density of the Universe at z=1.5-2.5, dominated by 250um detected and bump DOGs. For comparison, DOGs contribute 30% to the star-formation rate density of the Universe for all galaxies at z=1.5-2.5 with S24>=100uJy. DOGs have a large scatter about the star-formation main sequence and their specific star-formation rates show that the observed phase of star-formation could be responsible for their observed stellar mass at z~2.

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 on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km/s, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km/s (1000 km/s). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km/s are thought to require an active galactic nucleus (AGN) to drive them, about 2/3 of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources which are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less obscured AGN. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles implying neutral gas outflow masses of at least 2-4.5 x 10^8 Msun.

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

We present observations of the nearby spiral galaxy IC342 with the Herschel Spectral and Photometric Imaging Receiver (SPIRE) Fourier Transform Spectrometer. The spectral range afforded by SPIRE, 196-671 microns, allows us to access a number of 12CO lines from J=4--3 to J=13--12 with the highest J transitions observed for the first time. In addition we present measurements of 13CO, [CI] and [NII]. We use a radiative transfer code coupled with Bayesian likelihood analysis to model and constrain the temperature, density and column density of the gas. We find two 12CO components, one at 35 K and one at 400 K with CO column densities of 6.3x10^{17} cm^{-2} and 0.4x10^{17} cm^{-2} and CO gas masses of 1.26x10^{7} Msolar and 0.15x10^{7} Msolar, for the cold and warm components, respectively. The inclusion of the high-J 12CO line observations, indicate the existence of a much warmer gas component (~400 K) confirming earlier findings from H_{2} rotational line analysis from ISO and Spitzer. The mass of the warm gas is 10% of the cold gas, but it likely dominates the CO luminosity. In addition, we detect strong emission from [NII] 205microns and the {3}P_{1}->{3}P_{0} and {3}P_{2} ->{3}P_{1} [CI] lines at 370 and 608 microns, respectively. The measured 12CO line ratios can be explained by Photon-dominated region (PDR) models although additional heating by e.g. cosmic rays cannot be excluded. The measured [CI] line ratio together with the derived [C] column density of 2.1x10^{17} cm^{-2} and the fact that [CI] is weaker than CO emission in IC342 suggests that [CI] likely arises in a thin layer on the outside of the CO emitting molecular clouds consistent with PDRs playing an important role.

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

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.

We perform a stacking analysis of H-ATLAS data in order to obtain isothermal dust temperatures and rest-frame luminosities at 250um (L250), for 1599 radio sources over the H-ATLAS P1 GAMA area. The radio sample is generated using a combination of NVSS data and K-band UKIDSS-LAS data, over 0.01<z<0.8. The FIR properties of the sample are investigated as a function of L_1.4 GHz, z, projected radio-source size and radio spectral index. In order to search for stellar-mass dependent relations, we split the parent sample into those sources which are below and above 1.5 L_{K}^{*}. Correcting for stellar mass and redshift, we find no relation between the L250 and L_1.4GHz of radio AGN. This implies that a galaxy's nominal radio luminosity has little or no bearing on the SFR and/or dust mass content of the host system. The L250 of both the radio detected and non radio-detected galaxies (defined as those sources not detected at 1.4GHz but detected in the SDSS with r'<22) rises with increasing z. For mass and colour matched samples, sub-1.5 L_{K}^{*} and super-1.5 L_{K}^{*} radio-detected galaxies have 0.89\pm0.18 and 0.49\pm0.12 times the L250 of their non radio-detected counterparts. We explain these results in terms of the hotter, denser and richer halo environments massive radio-galaxies maintain and are embedded in. These results indicate that all massive radio galaxies (>1.5 L_{K}^{*}) may have systematically lower FIR luminosities (~25%) than their colour-matched non radio-detected counterparts. Compact radio sources (<30kpc) are associated with higher L250s and dust temperatures than their more extended (>30kpc) counterparts. The higher dust temperature suggests that this may be attributed to enhanced SFRs, but whether this is directly or indirectly due to radio activity (e.g. jet induced or merger-driven SF) is as yet unknown.

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.

Local luminous infrared (IR) galaxies (LIRGs) have both high star formation rates (SFR) and a high AGN (Seyfert and AGN/starburst composite) incidence. Therefore, they are ideal candidates to explore the co-evolution of black hole (BH) growth and star formation (SF) activity, not necessarily associated with major mergers. Here, we use Spitzer/IRS spectroscopy of a complete volume-limited sample of local LIRGs (distances of <78Mpc). We estimate typical BH masses of 3x10^7 M_sun using [NeIII]15.56micron and optical [OIII]5007A gas velocity dispersions and literature stellar velocity dispersions. We find that in a large fraction of local LIRGs the current SFR is taking place not only in the inner nuclear ~1.5kpc region, as estimated from the nuclear 11.3micron PAH luminosities, but also in the host galaxy. We next use the ratios between the SFRs and BH accretion rates (BHAR) to study whether the SF activity and BH growth are contemporaneous in local LIRGs. On average, local LIRGs have SFR to BHAR ratios higher than those of optically selected Seyferts of similar AGN luminosities. However, the majority of the IR-bright galaxies in the RSA Seyfert sample behave like local LIRGs. Moreover, the AGN incidence tends to be higher in local LIRGs with the lowest SFRs. All this suggests that in local LIRGs there is a distinct IR-bright star forming phase taking place prior to the bulk of the current BH growth (i.e., AGN phase). The latter is reflected first as a composite and then as a Seyfert, and later as a non-LIRG optically identified Seyfert nucleus with moderate SF in its host galaxy.

We present a multi-wavelength study of a 3.6 $\mu$m-selected galaxy sample in the Extended Groth strip. The sample is complete for galaxies with stellar mass $>10^{9.5}$ \Msun and redshift $0.4<z<1.2$. In this redshift range, the IRAC 3.6 $\mu$m band measures the rest-frame near-infrared band, permitting nearly unbiased selection with respect to both quiescent and star-forming galaxies. The numerous spectroscopic redshifts available in the EGS are used to train an Artificial Neural Network to estimate photometric redshifts. The distribution of photometric redshift errors is Gaussian with standard deviation ${\sim}0.025(1+z)$, and the fraction of redshift failures (${>}3\sigma$ errors) is about 3.5%. A new method of validation based on pair statistics confirms the estimate of standard deviation even for galaxies lacking spectroscopic redshifts. Basic galaxy properties measured include rest-frame $U-B$ colors, $B$- and $K$-band absolute magnitudes, and stellar masses. We divide the sample into quiescent and star-forming galaxies according to their rest-frame $U-B$ colors and 24 to 3.6 \micron\ flux density ratios and derive rest $K$-band luminosity functions and stellar mass functions for quiescent, star forming, and all galaxies. The results show that massive, quiescent galaxies were in place by $z\approx1$, but lower mass galaxies generally ceased their star formation at later epochs.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

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

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.

Using data from the mid-infrared to millimeter wavelengths for individual galaxies and for stacked ensembles at 0.5<z<2, we derive robust estimates of dust masses (Mdust) for main sequence (MS) galaxies, which obey a tight correlation between star formation rate (SFR) and stellar mass (M*), and for star-bursting galaxies that fall outside that relation. Exploiting the correlation of gas to dust mass with metallicity (Mgas/Mdust -Z), we use our measurements to constrain the gas content, CO-to-H2 conversion factors (a_co) and star formation efficiencies (SFE) of these distant galaxies. Using large statistical samples, we confirm that a_co and SFE are an order of magnitude higher and lower, respectively, in MS galaxies at high redshifts compared to the values of local galaxies with equivalently high infrared luminosities. For galaxies within the MS, we show that the variations of specific star formation rates (sSFR=SFR/M*) are driven by varying gas fractions. For relatively massive galaxies like those in our samples, we show that the hardness of the radiation field, <U>, which is proportional to the dust mass weighted luminosity (LIR/Mdust), and the primary parameter defining the shape of the SED, is equivalent to SFE/Z. For MS galaxies we measure this quantity, <U>, showing that it does not depend significantly on either the stellar mass or the sSFR. This is explained as a simple consequence of the existing correlations between SFR-M*, M*-Z and Mgas-SFR. Instead, we show that <U> (or LIR/Mdust) does evolve, with MS galaxies having harder radiation fields and thus warmer temperatures as redshift increases from z=0 to 2, a trend which can also be understood based on the redshift evolution of the M*-Z and SFR-M* relations. These results motivate the construction of a universal set of SED templates for MS galaxies which vary as a function of redshift with only one parameter, <U>.

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

We present observations of the CO[3-2] emission towards two massive and infrared luminous Lyman Break Galaxies at z = 3.21 and z = 2.92, using the IRAM Plateau de Bure Interferometer, placing first constraints on the molecular gas masses (Mgas) of non-lensed LBGs. Their overall properties are consistent with those of typical (Main-Sequence) galaxies at their redshifts, with specific star formation rates ~1.6 and ~2.2 Gyr^(-1), despite their large infrared luminosities L_IR ~2-3 x 10^12 Lsun derived from Herschel. With one plausible CO detection (spurious detection probability of 10^(-3)) and one upper limit, we investigate the evolution of the molecular gas-to-stellar mass ratio (Mgas/M*) with redshift. Our data suggest that the steep evolution of Mgas/M* of normal galaxies up to z~2 is followed by a flattening at higher redshifts, providing supporting evidence for the existence of a plateau in the evolution of the specific star formation rate at z > 2.5.

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

We report Herschel SPIRE (250, 350, and 500 micron) detections of 32 quasars with redshifts 0.5 < z < 3.6 from the Herschel Multi-tiered Extragalactic Survey. These sources are from a MIPS 24 micron flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame Spectral Energy Distributions (SEDs)from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25K to 60K, with a mean of 34K. The FIR luminosities range from 10^{11.3} to 10^{13.5} Lsun, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3--20 micron, rest-frame), and the bolometric luminosities derived using the 5100 A index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities.

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

Mrk938 is a luminous infrared galaxy in the local Universe believed to be the remnant of a galaxy merger. It shows a Seyfert 2 nucleus and intense star formation according to optical spectroscopic observations. We have studied this galaxy using new Herschel far-IR imaging data in addition to archival X-ray, UV, optical, near-IR and mid-IR data. Mid- and far-IR data are crucial to characterise the starburst contribution, allowing us to shed new light on its nature and to study the coexistence of AGN and starburst activity in the local Universe. The decomposition of the mid-IR Spitzer spectrum shows that the AGN bolometric contribution to the mid-IR and total infrared luminosity is small (Lbol(AGN)/LIR~0.02), which agrees with previous estimations. We have characterised the physical nature of its strong infrared emission and constrained it to a relatively compact emitting region of <2kpc. It is in this obscured region where most of the current star formation activity is taking place as expected for LIRGs. We have used Herschel imaging data for the first time to constrain the cold dust emission with unprecedented accuracy. We have fitted the integrated far-IR spectral energy distribution and derived the properties of the dust, obtaining a dust mass of 3x10^7Msun. The far-IR is dominated by emission at 35K, consistent with dust heated by the on-going star formation activity.

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

We report detection of four IRAC sources in the GOODS-South field with an extremely red color of H$-$[3.6]$>$4.5. The four sources are not detected in the deep HST WFC3 H-band image with H$_{limit}$=28.3 mag. We find that only 3 types of SED templates can produce such a red H$-$[3.6] color: a very dusty SED with the Calzetti extinction of A$_V$=16 mag at z=0.8; a very dusty SED with the SMC extinction of A$_V$=8 mag at z=2.0$\sim$2.2; and an 1Gyr SSP with A$_V \sim$0.8 at z=5.7. We argue that these sources are unlikely dusty galaxies at z$\leq$2.2 based on absent strong MIPS 24$\mu$m emission. The old stellar population model at z$>$4.5 remains a possible solution for the 4 sources. At z$>$4.5, these sources have stellar masses of Log(M$_*$/M$_{\odot}$)=10.6$\sim$11.2. One source, ERS-1, is also a type-II X-ray QSO with L$_{2-8keV}$=1.6$\times 10^{44}$ erg s$^{-1}$. One of the four sources is an X-ray QSO and another one is a HyperLIRG, suggesting a galaxy-merging scenario for the formation of these massive galaxies at high redshifts.

We present broad band photometry and photometric redshifts for 187611 sources located in ~0.5deg^2 in the Lockman Hole area. The catalog includes 389 X-ray detected sources identified with the very deep XMM-Newton observations available for an area of 0.2 deg^2. The source detection was performed on the Rc, z' and B band images and the available photometry is spanning from the far ultraviolet to the mid infrared, reaching in the best case scenario 21 bands. Astrometry corrections and photometric cross-calibrations over the entire dataset allowed the computation of accurate photometric redshifts. Special treatment is undertaken for the X-ray sources, the majority of which is active galactic nuclei. Comparing the photometric redshifts to the available spectroscopic redshifts we achieve for normal galaxies an accuracy of \sigma_{\Delta z/(1+z)}=0.036, with 12.7% outliers, while for the X-ray detected sources the accuracy is \sigma_{\Delta z/(1+z)}=0.069, with 18.3% outliers, where the outliers are defined as sources with |z_{phot}-z_{spec}|>0.15 (1+z_{spec})}. These results are a significant improvement over the previously available photometric redshifts for normal galaxies in the Lockman Hole, while it is the first time that photometric redshifts are computed and made public for AGN for this field.

We quantify the active galactic nucleus (AGN) contribution to the mid-infrared (mid-IR) and the total infrared (IR, 8-1000micron) emission in a complete volume-limited sample of 53 local luminous infrared galaxies (LIRGs). We decompose the Spitzer Infrared Spectrograph (IRS) low-resolution 5-38micron spectra of the LIRGs into AGN and starburst components using clumpy torus models and star-forming galaxy templates, respectively. We find that 50% (25/50) of local LIRGs have an AGN component detected with this method. There is good agreement between these AGN detections through mid-IR spectral decomposition and other AGN indicators, such as the optical spectral class, mid-IR spectral features and X-ray properties. Taking all the AGN indicators together, the AGN detection rate in the individual nuclei of LIRGs is ~62%. The derived AGN bolometric luminosities are in the range L_bol(AGN)=0.4 -50x10^{43} erg/s. The AGN bolometric contribution to the IR luminosities of the galaxies is generally small, with 70% of LIRGs having L_bol(AGN)/L_IR<0.05. Only ~8% of local LIRGs have a significant AGN bolometric contribution L_bol(AGN)/L_IR > 0.25. From the comparison of our results with literature results of ultraluminous infrared galaxies, we confirm that in the local universe the AGN bolometric contribution to the IR luminosity increases with the IR luminosity of the galaxy/system. If we add up the AGN bolometric luminosities we find that AGNs only account for 5%^{+8%}_{-3%} of the total IR luminosity produced by local LIRGs (with and without AGN detections). This proves that the bulk of the IR luminosity of local LIRGs is due to star formation activity. Taking the newly determined IR luminosity density of LIRGs in the local universe, we then estimate an AGN IR luminosity density of Omega_IR(AGN) = 3x10^5 L_sun Mpc^{-3}$ in LIRGs.

We explore the gas-to-dust mass ratio (G/D) and the CO luminosity-to-Mgas conversion factor (a_co) of two well studied galaxies in the GOODS-N field, that are expected to have different star forming modes, the starburst GN20 at z=4.05 and the normal star-forming galaxy BzK-21000 at z=1.52. Detailed sampling is available for their Rayleigh-Jeans emission via ground based mm interferometry (1.1-6.6mm) along with Herschel, PACS and SPIRE data that probe the peak of their infrared emission. Using the physically motivated Draine & Li (2007) models, as well as a modified black body function, we measure the dust mass (Md) of the sources and find 2.0^{+0.7}_{-0.6} x 10^{9} Msun for GN20 and 8.6^{+0.6}_{-0.9} x 10^{8} Msun for BzK-21000. The addition of mm data reduces the uncertainties of the derived Md by a factor of ~2, allowing the use of the local G\D vs metallicity relation to place constraints on the a_co values of the two sources. For GN20 we derive a conversion factor of a_co < 1.0 Msun pc^{-2}(K km s^{-1})^{-1}, consistent with that of local ULIRGs, while for BzK-21000 we find a considerably higher value, a_co ~4.0 Msun pc^{-2}(K km s^{-1})^{-1}, in agreement with an independent kinematic derivation reported previously. The implied star formation efficiency is ~25 Lsun/Msun for BzK-21000, a factor of ~5-10 lower than that of GN20. The findings for these two sources support the existence of different disk-like and starburst star-formation modes in distant galaxies, although a larger sample is required to draw statistically robust results

Ultra Steep Spectrum (USS) radio sources have been successfully used to select powerful radio sources at high redshifts (z>~2). Typically restricted to large-sky surveys and relatively bright radio flux densities, it has gradually become possible to extend the USS search to sub-mJy levels, thanks to the recent appearance of sensitive low-frequency radio facilities. Here a first detailed analysis of the nature of the faintest USS sources is presented. By using Giant Metrewave Radio Telescope and Very Large Array radio observations of the Lockman Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz integrated fluxes above 100 microJy, is assembled. Deep infrared data at 3.6 and 4.5 micron from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used to reliably identify counterparts for 48 (83%) of these sources, showing an average total magnitude of [3.6](AB)=19.8 mag. Spectroscopic redshifts for 14 USS sources, together with photometric redshift estimates, improved by the use of the deep SERVS data, for a further 19 objects, show redshifts ranging from z=0.1 to z=2.8, peaking at z~0.6 and tailing off at high redshifts. The remaining 25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, with 10 sources undetected at 3.6 and 4.5 micron (typically [3.6]>22-23 mag, from local measurements), which suggests the likely existence of higher redshifts among the sub-mJy USS population. The comparison with the Square Kilometre Array Design Studies Simulated Skies models indicate that Fanaroff-Riley type I radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of the faintest USS population, and raises the possibility that the high efficiency of the USS technique for the selection of high redshift sources remains even at the sub-mJy level.

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

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

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

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

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

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

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

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

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

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.

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

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

SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 {\mu}m, and an imaging Fourier transform spectrometer (FTS) covering 194-671 {\mu}m (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The spectrometer wavelength accuracy is determined to be better than 1/10th of the line FWHM. The astrometric accuracy in SPIRE maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models of astronomical standards.

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

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

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

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

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

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

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

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

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

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

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

We present a multi-wavelength, UV-to-radio analysis for a sample of massive (M$_{\ast}$ $\sim$ 10$^{10}$ M$_\odot$) IRAC- and MIPS 24$\mu$m-detected Lyman Break Galaxies (LBGs) with spectroscopic redshifts z$\sim$3 in the GOODS-North field (L$_{\rm UV}$$>1.8\times$L$^{\ast}_{z=3}$). For LBGs without individual 24$\mu$m detections, we employ stacking techniques at 24$\mu$m, 1.1mm and 1.4GHz, to construct the average UV-to-radio spectral energy distribution and find it to be consistent with that of a Luminous Infrared Galaxy (LIRG) with L$\rm_{IR}$=4.5$^{+1.1}_{-2.3}$$\times 10^{11}$ L$_{\odot}$ and a specific star formation rate (SSFR) of 4.3 Gyr$^{-1}$ that corresponds to a mass doubling time $\sim$230 Myrs. On the other hand, when considering the 24$\mu$m-detected LBGs we find among them galaxies with L$\rm_{IR}> 10^{12}$ L$_{\odot}$, indicating that the space density of $z\sim$3 UV-selected Ultra-luminous Infrared Galaxies (ULIRGs) is $\sim$(1.5$\pm$0.5)$\times 10^{-5}$ Mpc$^{-3}$. We compare measurements of star formation rates (SFRs) from data at different wavelengths and find that there is tight correlation (Kendall's $\tau >$ 99.7%) and excellent agreement between the values derived from dust-corrected UV, mid-IR, mm and radio data for the whole range of L$\rm_{IR}$ up to L$\rm_{IR}$ $\sim$ 10$^{13}$ L$_{\odot}$. This range is greater than that for which the correlation is known to hold at z$\sim$2, possibly due to the lack of significant contribution from PAHs to the 24$\mu$m flux at $z\sim$3. The fact that this agreement is observed for galaxies with L$\rm_{IR}$ $>$ 10$^{12}$ L$_{\odot}$ suggests that star-formation in UV-selected ULIRGs, as well as the bulk of star-formation activity at this redshift, is not embedded in optically thick regions as seen in local ULIRGs and submillimeter-selected galaxies at $z=2$.

We present results of a large survey of the mid--IR properties of 248 Lyman Break Galaxies with confirmed spectroscopic redshift using deep Spitzer/IRAC observations in six cosmological fields. We model the Spectral Energy Distributions (SEDs) employing a revised version of the Bruzual and Charlot synthesis population code that incorporates a new treatment of the TP--AGB phase (CB07). Our primary aim is to investigate the impact of the AGB phase in the stellar masses of the LBGs, and compare our new results with previous stellar mass estimates. Based on the new CB07 code we find that the stellar masses of LBGs are smaller on average by a factor of ~1.4 compared to previous estimates. LBGs with 8um and/or 24um detections show higher masses (M~10^11 Mo) than LBGs faint in the IRAC bands (M~10^9 Mo). The ages of these massive LBGs are considerably higher than the rest of the population, indicating that they have been star-forming for at least ~1 Gyr. We also show how the addition of the IRAC bands, improves the accuracy of the estimated stellar masses and reduced the scatter on the derived M/L ratios. In particular, we present a tight correlation between the 8um IRAC band (rest-frame K for galaxies at z~3) and the stellar mass. We calculate the number density of massive (M > 10^11 Mo) LBGs and find it to be $\Phi$= (1.12 $\pm$ 0.4) x 10^(-5) Mpc^(-3), ~1.5 times lower than that found by previous studies. Finally, based on UV-corrected SFRs we investigate the SFR-stellar mass correlation at z~3, find it similar to the one observed at other redshifts and show that our data place the peak of the evolution of the specific star formation rate at z~3.

We present new mid-infrared and submillimetre observations for a sample of eight high redshift type-2 QSOs located in the Chandra Deep Field South. The sources are X-ray absorbed with luminosities in excess of 10^44 erg/s. Two of the targets have robust detections, S/N > 4, while a further three targets are marginally detected with S/N > =2.5. All sources are detected in multiple mid-infrared bands with the Spitzer Space Telescope. The multiwavelength spectral energy distributions (SEDs) of the type-2 QSOs are compared to those of two local ultraluminous galaxies (Arp220 and IR22491) in order to assess contributions from a star-forming component in various parts of the SED. We suggest that their submillimetre emission is possibly due to a starburst while a large fraction of the mid-infrared energy is likely to originate in the obscured central quasar. Using the mid-infrared and submm observations we derive infrared luminosities which are found to be in excess of L>10^12Lsun. The submillimetre (850micron) to X-ray (2 keV) spectral indices (alpha_SX) span a wide range. About half of the type-2 QSOs have values typical for a Compton-thick AGN with only 1 per cent of the nuclear emission seen through scattering and, the remaining with values typical of submm-bright galaxies. Combining the available observational evidence we outline a possible scenario for the early stages of evolution of these sources.

We analyze a sample of galaxies chosen to have F24>0.5 mJy and satisfy a certain IRAC color criterion. IRS spectra yield redshifts, spectral types, and PAH luminosities, to which we add multi-wavelength broadband photometry. Stellar population modeling and IRS spectra together demonstrate that the double criteria used to select this sample have efficiently isolated massive star-forming galaxies at z~1.9. This is the first starburst-dominated ULIRG sample at high redshift with total infrared luminosity measured directly from FIR and millimeter photometry, and as such gives us the first accurate view of broadband SEDs for starburst galaxies at extremely high luminosity and at all wavelengths. Similar broadband data are assembled for three other galaxy samples -- local starburst galaxies, local AGN/ULIRGS, and a second 24mu-luminous z~1.9 sample dominated by AGN. L(PAH)/L(IR) for the new z~1.9 starburst sample is the highest ever seen, some three times higher than in local starbursts, whereas in AGNs this ratio is depressed below the starburst trend, often severely. Several pieces of evidence imply that AGNs exist in this starburst dominated sample, except two of which even host very strong AGN, while they still have very strong PAH emission. The ACS images show most objects have very extended morphologies in the rest-frame UV band, thus extended distribution of PAH molecules. Such an extended distribution prevents further destruction PAH molecules by central AGNs. We conclude that objects in this sample are ULIRGs powered mainly by starburst; and the total infrared luminosity density contributed by this type of objects is 0.9-2.6x 10^7 Lsol/Mpc^3.

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

We present millimetre observations of a sample of 12 high redshift ultraluminous infrared galaxies (ULIRGs) in the Extended Growth Strip (EGS). These objects were initially selected on the basis of their observed mid--IR colours (0.0 < [3.6]-[4.5] < 0.4 and -0.7 < [3.6]-[8.0] < 0.5) to lie at high redshift 1.5 < z < 3, and subsequent 20-38 micron mid-IR spectroscopy confirms that they lie in a narrow redshift window centered on z=2. We detect 9/12 of the objects in our sample at high significance (>3 sigma) with a mean 1200\micron flux of <F_1200> = 1.6+/-0.1 mJy. Our millimetre photometry, combined with existing far-IR photometry from the Far-IR Deep Extragalactic Legacy (FIDEL) Survey and accurate spectroscopic redshifts, places constraints both sides of the thermal dust peak. This allows us to estimate the dust properties, including the far--IR luminosity, dust temperature, and dust mass. We find that our sample is similar to other high-z and intermediate-z ULIRGs, and local systems, but has a different dust selection function than submillimeter-selected galaxies. Finally, we use existing 20cm radio continuum imaging to test the far-IR/radio correlation at high redshift. We find that our sample is consistent with the local relation, implying little evolution. Furthermore, this suggests that our sample selection method is efficient at identifying ultraluminous, starburst--dominated systems within a very narrow redshift range centered at z~2.

We present photometric analysis of deep mid-infrared observations obtained by Spitzer/IRAC covering the fields Q1422+2309, Q2233+1341, DSF2237a,b, HDFN, SSA22a,b and B20902+34, giving the number counts and the depths for each field. In a sample of 751 LBGs lying in those fields, 443, 448, 137 and 152 are identified at 3.6microns, 4.5microns, 5.8microns, 8.0microns IRAC bands respectively, expanding their spectral energy distribution to rest-near-infrared and revealing that LBGs display a variety of colours. Their rest-near-infrared properties are rather inhomogeneous, ranging from those that are bright in IRAC bands and exhibit [R]-[3.6] > 1.5 colours to those that are faint or not detected at all in IRAC bands with [R]-[3.6] < 1.5 colours and these two groups of LBGs are investigated. We compare the mid-IR colours of the LBGs with the colours of star-forming galaxies and we find that LBGs have colours consistent with star-foming galaxies at z~3. The properties of the LBGs detected in the 8microns IRAC band (rest frame K-band) are examined separately, showing that they exhibit redder [R]-[3.6] colours than the rest of the population and that IRAC 8microns band can be used as a diagnostic tool, to separate AGN dominated objects from normal star-forming galaxies at z~3

We examine the infrared properties of 43 high redshift (0.1<z<1.2), infrared-luminous galaxies in the Extended Groth Strip (EGS), selected by a deep 70 micron survey with the Multiband Imaging Photometer on Spitzer (MIPS). In addition and with reference to starburst-type Spectral Energy Distributions (SEDs), we derive a set of equations for estimating the total infrared luminosity (L_IR) in the range 8-1000 microns using photometry from at least one MIPS band. 42 out of 43 of our sources' optical/infrared SEDs (lambda_observed < 160 microns) are starburst-type, with only one object displaying a prominent power-law near-infrared continuum. For a quantitative analysis, models of radiation transfer in dusty media are fit onto the infrared photometry, revealing that the majority of galaxies are represented by high extinction, A_v>35 and for a large fraction (~50 per cent) the SED turns over into the Rayleigh-Jeans regime at wavelengths longward of 90 microns. For comparison, we also fit semi-empirical templates based on local galaxy data, however, these underestimate the far-infrared SED shape by a factor of at least 2 and in extreme cases up to 10 for the majority (~70 per cent) of the sources. Further investigation of SED characteristics reveals that the mid-infrared (70/24 microns) continuum slope is decoupled from various galaxy properties such as the total infrared luminosity and far-infrared peak, quantified by the L_160/L_70 ratio. In view of these results, we propose that these high-redshift galaxies have different properties to their local counterparts, in the sense that large amounts of dust cause heavy obscuration and are responsible for an additional cold emissive component, appearing as a far-infrared excess in their SEDs.

We study the stellar and star formation (SF) properties of the host galaxies of 58 X-ray selected AGN in the GOODS portion of the Chandra Deep Field South (CDF-S) region at z~0.5-1.4. The AGN are selected such that their rest-frame UV to near-IR SEDs are dominated by stellar emission, i.e., they show a prominent 1.6micron bump, thus minimizing the AGN emission 'contamination'. This AGN population comprises approximately 50% of the X-ray selected AGN at these redshifts. Using models of stellar and dust emission we model their SEDs to derive stellar masses (M*) and total (UV+IR) star formation rates (SFR). AGN reside in the most massive galaxies at the redshifts probed here. Their characteristic stellar masses (M*~7.8x10^10 Msun and M*~1.2x10^11 Msun at median z of 0.67 and 1.07, respectively) appear to be representative of the X-ray selected AGN population at these redshifts, and are intermediate between those of local type 2 AGN and high redshift (z~2) AGN. The inferred black hole masses (MBH~ 2x10^8 Msun) of typical AGN are similar to those of optically identified quasars at similar redshifts. Since the AGN in our sample are much less luminous (L(2-10keV)<10^44 erg/s) than quasars, typical AGN have low Eddington ratios (eta~0.01-0.001). This suggests that, at least at intermediate z, the cosmic AGN 'downsizing' is due to both a decrease in the characteristic stellar mass of typical host galaxies, and less efficient accretion. Finally there is no strong evidence in AGN host galaxies for either highly suppressed SF (expected if AGN played a role in quenching SF) or elevated SF when compared to mass selected (i.e., IRAC-selected) galaxies of similar stellar masses and redshifts. This may be explained by the fact that galaxies with M*~5x 10^10 - 5x10^11 Msun are still being assembled at the redshifts probed here.

Using a combination of deep MID-IR observations obtained by IRAC, MIPS and IRS on board Spitzer we investigate the MID-IR properties of Lyman Break Galaxies (LBGs) at z~3, establish a better understanding of their nature and attempt a complete characterisation of the population. With deep mid-infrared and optical observations of ~1000 LBGs covered by IRAC/MIPS and from the ground respectively, we extend the spectral energy distributions (SEDs) of the LBGs to mid-infrared. Spitzer data reveal for the first time that the mid-infrared properties of the population are inhomogeneous ranging from those with marginal IRAC detections to those with bright rest-frame near-infrared colors and those detected at 24mu MIPS band revealing the newly discovered population of the Infrared Luminous Lyman Break Galaxies (ILLBGs). To investigate this diversity, we examine the photometric properties of the population and we use stellar population synthesis models to probe the stellar content of these galaxies. We find that a fraction of LBGs have very red colors and large estimated stellar masses M* > 5x10^10Mo. We discuss the link between these LBGs and submm-luminous galaxies and we report the detection of rest frame 6.2 and 7.7mu emission features arising from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of an infrared-luminous Lyman break galaxy at z=3.01.

A SST survey in the NOAO Deep-Wide Field in Bo\"otes provides a complete, 8-micron-selected sample of galaxies to a limiting (Vega) magnitude of 13.5. In the 6.88 deg$^2$ field sampled, 79% of the 4867 galaxies have spectroscopic redshifts, allowing an accurate determination of the local (z<0.3) galaxy luminosity function. Stellar and dust emission can be separated on the basis of observed galaxy colors. Dust emission (mostly PAH) accounts for 80% of the 8 micron luminosity, stellar photospheres account for 19%, and AGN emission accounts for roughly 1 %. A sub-sample of the 8 micron-selected galaxies have blue, early-type colors, but even most of these have significant PAH emission. The luminosity functions for the total 8 micron luminosity and for the dust emission alone are both well fit by Schechter functions. For the 8 micron luminosity function, the characteristic luminosity is \nu L_{\nu}^*(8.0 \micron) = 1.8 \times 10^{10}$ \Lsun while for the dust emission alone it is 1.6 x 10^{10}$ \Lsun\null. The average 8 \micron luminosity density at z<0.3 is 3.1 x 10^7 \Lsun Mpc^{-3}, and the average luminosity density from dust alone is 2.5 x 10^7 \Lsun Mpc^{-3}. This luminos ity arises predominantly from galaxies with 8 \micron luminosities ($\nu L_{\nu}$) between $2\times 10^9$ and $2 x 10^{10}$ \Lsun, i.e., normal galaxies, not LIRGs or ULIRGs.

We present 0.5 -160 micron Spectral Energy Distributions (SEDs) of galaxies, detected at 70microns with the Multiband Imaging Photometer for Spitzer (MIPS), using broadband imaging data from Spitzer and ground-based telescopes. Spectroscopic redshifts, in the range 0.2<z<1.5, have been measured as part of the Deep Extragalactic Evolutionary Probe2 (DEEP2) project. Based on the SEDs we explore the nature and physical properties of the sources. Using the optical spectra we derive Hbeta and [OII]-based Star Formation Rates (SFR) which are 10-100 times lower than SFR estimates based on IR and radio. The median offset in SFR between optical and IR is reduced by a factor of ~3 when we apply a typical extinction corrections. We investigate mid-to-far infrared correlations for low redshift (>0.5) and high redshift (0.5<z<1.2) bins. Using this unique ``far-infrared'' selected sample we derive an empirical mid to far-infrared relationship that can be used to estimate the infrared energy budget of galaxies in the high-redshift universe. Our sample can be used as a template to translate far-infrared luminosities into bolometric luminosities for high redshift objects.

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

We report the detection of rest--frame 6.2 and 7.7 \micron emission features arising from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of an infrared-luminous Lyman break galaxy at z=3.01. This is currently the highest redshift galaxy where these PAH emission features have been detected. The total infrared luminosity inferred from the MIPS 24 \micron and radio flux density is 2$\times10^{13}$ L$_{\odot}$, which qualifies this object as a so--called hyperluminous infrared galaxy (HyLIRG). However, unlike local HyLIRGs which are generally associated with QSO/AGNs and have weak or absent PAH emission features, this HyLIRG has very strong 6.2 and 7.7 \micron PAH emission. We argue that intense star formation dominates the infrared emission of this source, although we cannot rule out the presence of a deeply obscured AGN. This LBG appears to be a distorted system in the HST ACS F606W and F814W images, possibly indicating that a significant merger or interaction is driving the large IR luminosity.

We study 87 Extremely Red Objects (EROs), selected both to have color redder than R-[3.6] = 4.0, and to have confirmed spectroscopic redshifts. Together, these two constraints result in this sample populating a fairly narrow redshift range at 0.76 < z < 1.42. The key new ingredient included here is deep Spitzer Space Telescope InfraRed Array Camera (IRAC) data. Based on [3.6]-[8.0] color, we demonstrate that it is possible to classify EROs into early-type, dusty starburst, or power-law (AGN) types. We present ultraviolet to mid-infrared spectral energy distributions (SEDs) and Advanced Camera for Surveys (ACS) images, both of which support our simple IRAC color classification.

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

We describe the spectral energy distributions (SEDs) of Lyman Break Galaxies (LBGs) at z~3 using deep mid-infrared and optical observations of the Extended Groth Strip, obtained with IRAC and MIPS on board Spitzer and from the ground, respectively. We focus on LBGs with detections at all four IRAC bands, in particular the 26 galaxies with IRAC 8 micron band (rest--frame K-band) detections. We use stellar population synthesis models and probe the stellar content of these galaxies. Based on best--fit continuous star-formation models we derive estimates of the stellar mass for these LBGs. As in previous studies, we find that a fraction of LBGs have very red colors and large estimated stellar masses (M_stellar > 5 x 10 ^{10} M_solar): the present Spitzer data allow us, for the first time, to study these massive LBGs in detail. We discuss the link between these LBGs and submm-luminous galaxies.We find that the number density of these massive LBGs at high redshift is higher than predicted by current semi-analytic models of galaxy evolution.

We identify 17 possible 8.0 micron-selected counterparts to the submillimeter galaxies in the CUDSS 14-hour field, derived from deep imaging carried out with the IRAC and MIPS instruments aboard the Spitzer Space Telescope. Ten of the 17 counterparts are not the same as those previously identified at shorter wavelengths. We argue that 8.0 micron selection offers a better means for identifying counterparts to submillimeter galaxies than near-infrared or optical selection. Based on the panchromatic SEDs, most counterparts appear to be powered by ongoing star formation. Power-law fits to the SEDs suggest that five objects in the 8.0 micron-selected counterpart sample harbor dominant AGNs; a sixth object is identified as a possible AGN. The 3.6 to 8.0 micron colors of the infrared-selected counterparts are significantly redder than the general IRAC galaxy population in the CUDSS 14-hour field.

We investigate the nature of a sample of 92 Spitzer/MIPS 24 micron selected galaxies in the CDFS, showing power law-like emission in the Spitzer/IRAC 3.6-8 micron bands. The main goal is to determine whether the galaxies not detected in X-rays (47% of the sample) are part of the hypothetical population of obscured AGN not detected even in deep X-ray surveys. The majority of the IR power-law galaxies are ULIRGs at z>1, and those with LIRG-like IR luminosities are usually detected in X-rays. The optical to IR spectral energy distributions (SEDs) of the X-ray detected galaxies are almost equally divided between a BLAGN SED class (similar to an optically selected QSO) and a NLAGN SED (similar to the BLAGN SED but with an obscured UV/optical continuum). A small fraction of SEDs resemble warm ULIRG galaxies (e.g., Mrk231). Most galaxies not detected in X-rays have SEDs in the NLAGN+ULIRG class as they tend to be optically fainter, and possibly more obscured. Moreover, the IR power-law galaxies have SEDs significantly different from those of high-z (z_sp>1) IR (24 micron) selected and optically bright (VVDS I_AB<=24) star-forming galaxies whose SEDs show a very prominent stellar bump at 1.6 micron. The galaxies detected in X-rays have 2-8 keV rest-frame luminosities typical of AGN. The galaxies not detected in X-rays have global X-ray to mid-IR SED properties that make them good candidates to contain IR bright X-ray absorbed AGN. If all these sources are actually obscured AGN, we would observe a ratio of obscured to unobscured 24 micron detected AGN of 2:1, whereas models predict a ratio of up to 3:1. Additional studies using Spitzer to detect X-ray-quiet AGN are likely to find more such obscured sources.

We present low resolution near-infrared spectroscopy of an unbiased sample of 24 ultraluminous infrared galaxies (ULIRGs), selected from samples previously observed spectroscopically in the mid-infrared with the Infrared Space Observatory (ISO). Qualitatively, the near-infrared spectra resemble those of starbursts. Only in one ULIRG, IRAS 04114-5117E, do we find spectroscopic evidence for AGN activity. The spectroscopic classification in the near-infrared is in very good agreement with the mid-infrared one. For a subset of our sample for which extinction corrections can be derived from Pa-alpha and Br-gamma, we find rather high Pa-alpha luminosities, in accordance with the powering source of these galaxies being star formation.[Fe] emission is strong in ULIRGs and may be linked to starburst and superwind activity. Additionally, our sample includes two unusual objects. The first, IRAS F00183-7111, exhibits extreme [Fe] emission and the second, IRAS F23578-5307, is according to our knowledge one of the most luminous infrared galaxies in H2 rotation-vibration emission.

A deep mid- and far-infrared survey in the Extended Groth Strip (EGS) area gives 3.6 to 8\micron flux densities or upper limits for 253 Lyman Break Galaxies (LBGs). The LBGs are a diverse population but with properties correlated with luminosity. The LBGs show a factor of 30 range in indicated stellar mass and a factor of 10 in apparent dust content relative to stellar mass. About 5% of LBGs are luminous at all wavelengths with powerful emission at rest 6\micron. In the rest 0.9 to 2\micron spectral range these galaxies have stellar spectral slopes with no sign of an AGN power law component, suggesting that their emission is mainly powered by intensive star formation. Galaxies in this luminous population share the infrared properties of cold SCUBA sources: both are massive and dusty starburst galaxies at $2<z<3$; their stellar mass is larger than $10^{11} M_{\odot}$. We suggest that these galaxies are the progenitors of present-day giant elliptical galaxies, with a substantial fraction of their stars already formed at $z \approx 3$.

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 report on the first SCUBA detection of a Type 2 QSO at z=3.660 in the Chandra Deep Field South. This source is X-ray absorbed, shows only narrow emission lines in the optical spectrum and is detected in the sub-mm: it is the ideal candidate in an evolution scheme for AGN (e.g. Fabian (1999); Page et al. (2004)) of an early phase corresponding to the main growth of the host galaxy and formation of the central black hole. The overall photometry (from the radio to the X-ray energy band) of this source is well reproduced by the spectral energy distribution (SED) of NGC 6240, while it is incompatible with the spectrum of a Type 1 QSO (3C273) or a starburst galaxy (Arp 220). Its sub-mm (850 \mu m) to X-ray (2 keV) spectral slope (alpha_SX) is close to the predicted value for a Compton-thick AGN in which only 1% of the nuclear emission emerges through scattering. Using the observed flux at 850 \mu m we have derived a SFR=550--680 M$_odot/yr and an estimate of the dust mass, M$_dust=4.2 10^8 M_odot

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

We present Spitzer IRAC (3.6, 4.5, 5.8, 8.0 um) and MIPS (24 um) observations of the SCUBA submillimeter sources and uJy VLA radio sources in a 5'x5' area in the Lockman Hole East region. Out of the ~40 SCUBA/VLA sources in the field, Spitzer counterparts were detected for nearly all except for the few low-weight SCUBA detections. We show that the majority (80-90%) of the detected sources are cold (i.e., starburst-like) infrared-luminous galaxies (L_{IR} > 10^{11} Lsun) at redshift 0.5<z<3.5, whose star-formation rate density (SFRD) is comparable to that of the optically-selected star-forming galaxies.

We investigate the luminous X-ray sources in the Lockman Hole (LH) and the Extended Groth Strip (EGS) detected at 24microns using MIPS and also with IRAC on board Spitzer. We assemble optical/infrared spectral energy distributions (SEDs) for 45 X-ray/24micron sources in the EGS and LH. Only about 1/4 of the hard X-ray/24micron sources show pure type 1 AGN SEDs. More than half of the X-ray/24micron sources have stellar-emission-dominated or obscured SEDs, similar to those of local type 2 AGN and spiral/starburst galaxies. One-third of the sources detected in hard X-rays do not have a 24micron counterpart. Two such sources in the LH have SEDs resembling those of S0/elliptical galaxies. The broad variety of SEDs in the optical-to-Spitzer bands of X-ray selected AGN means that AGN selected according to the behavior in the optical/infrared will have to be supplemented by other kinds of data (e.g., X-ray) to produce unbiased samples of AGN.

We present Spitzer observations in five wavebands between 3.6 and 24um of an unbiased sample of 9 luminous, dusty galaxies selected at 1200um by the MAMBO camera on the IRAM 30-m telescope, a population akin to the well-known submm or `SCUBA' galaxies (hereafter SMGs). Owing to the coarse resolution of submm/mm instrumentation, SMGs have traditionally been difficult to identify at other wavelengths. We compare our multi-wavelength catalogs to show that the overlap between 24 and 1200um must be close to complete at these flux levels. We find that all (4/4) of the most secure >=4sigma SMGs have robust >=4sigma counterparts at 1.4GHz, while the fraction drops to 7/9 using all >=3sigma SMGs. We show that combining mid-IR and marginal (>=3sigma) radio detections provides plausible identifications in the remaining cases, enabling us to identify the complete sample. Accretion onto an obscured central engine is betrayed by the shape of the mid-IR continuum emission for several sources, confirming Spitzer's potential to weed out active galaxies. We demonstrate the power of a S(24um)/S(8um) vs S(8um)/S(4.5um) color-color plot as a diagnostic for this purpose. However, we conclude that the majority (~75%) of SMGs have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. Sensitive 24-um observations are clearly a useful route to identify and characterize reliable counterparts to high-redshift far-IR-bright galaxies, complementing what is possible via deep radio imaging.

We present preliminary results on 24micron detections of luminous infrared galaxies at z>1 with the Multiband Imaging Photometer for Spitzer (MIPS). Observations were performed in the Lockman Hole and the Extended Groth Strip (EGS), and were supplemented by data obtained with the Infrared Array Camera (IRAC) between 3 and 9microns. The positional accuracy of ~2arcsec for most MIPS/IRAC detections provides unambiguous identifications of their optical counterparts. Using spectroscopic redshifts from the Deep Extragalactic Evolutionary Probe survey, we identify 24micron sources at z>1 in the EGS, while the combination of the MIPS/IRAC observations with $BVRIJHK$ ancillary data in the Lockman Hole also shows very clear cases of galaxies with photometric redshifts at 1<z<2.5. The observed 24micron fluxes indicate infrared luminosities greater than 10^11 L_sol, while the data at shorter wavelengths reveal rather red and probably massive (M>=M*) galaxy counterparts. It is the first time that this population of luminous objects is detected up to z~2.5 in the infrared. Our work demonstrates the ability of the MIPS instrument to probe the dusty Universe at very high redshift, and illustrates how the forthcoming Spitzer deep surveys will offer a unique opportunity to illuminate a dark side of cosmic history not explored by previous infrared experiments.

IRAC imaging of a 4'7x4'7 area in the Lockman Hole detected over 400 galaxies in the IRAC 3.6 micron and 4.5 micron bands, 120 in the 5.8 micron, and 80 in the 8 micron bandin 30 minutes of observing time. Color-color diagrams suggest that about half of these galaxies are at redshifts 0.6<z<1.3 with about a quarter at higher redshifts (z>1.3). We also detect IRAC counterparts for 6 of the 7 SCUBA sources and all 9 XMM sources in this area. The detection of the counterparts of the SCUBA sources and galaxies at z>1.3 demonstrates the ability of IRAC to probe the universe at very high redshifts.

We present sub-millimeter statistical detections of galaxies discovered in the 5'x5' Spitzer Early Release Observations (to 4-15 microJy 5 sigma at 3.6-8 microns, 170 microJy at 24 microns) through a stacking analysis of our reanalysed SCUBA 8mJy survey maps, and a Spitzer identification of a new sub-millimeter point source in the 8mJy survey region. For sources detected at 5.8 or 8 microns (154 and 111 sources respectively), we detect positive skews in the sub-millimeter flux distributions at 99.2-99.8% confidence using Kolmogorov-Smirnov tests, at both 850 microns and 450 microns. We also marginally detect the Spitzer 24 micron galaxies at 850 microns at 97% confidence, and place limits on the mean sub-millimeter fluxes of the 3.6 and 4.5 micron sources. Integrating the sub-millimeter fluxes of the Spitzer populations, we find the 5.8 micron galaxies contribute 0.12 +/- 0.05 nW/m^2/sr to the 850 micron background, and 2.4 +/- 0.7 nW/m^2/sr to the 450 micron background; similar contributions are made by the 8 micron-selected sample. We infer that the populations dominating the 5.8 and 8 micron extragalactic background light also contribute around a quarter of the 850 micron background and the majority of the 450 micron background.

We investigate Extremely Red Objects (EROs) using near- and mid-infrared observations in five passbands (3.6 to 24 micron) obtained from the Spitzer Space Telescope, and deep ground-based R and K imaging. The great sensitivity of the IRAC camera allows us to detect 64 EROs in only 12 minutes of IRAC exposure time, by means of an R-[3.6] color cut (analogous to the traditional red R-K cut). A pure infrared K-[3.6] red cut detects a somewhat different population and may be more effective at selecting z > 1.3 EROs. We find 17% of all galaxies detected by IRAC at 3.6 or 4.5 micron to be EROs. These percentages rise to about 40% at 5.8 micron, and about 60% at 8.0 micron. We utilize the spectral bump at 1.6 micron to divide the EROs into broad redshift slices using only near-infrared colors (2.2/3.6/4.5 micron). We conclude that two-thirds of all EROs lie at redshift z > 1.3. Detections at 24 micron imply that at least 11% of 0.6 < z < 1.3 EROs and at least 22% of z > 1.3 EROs are dusty star-forming galaxies.

We present ISO PHOT-S spectra of four galaxies known or suspected to host a central AGN, selected from the initial Iras/Rosat sample of Boller et al. (1992). Two of them had no obvious Seyfert features in their previous optical spectra: IRAS 14201+2956, and IRAS 21582+1018. The latter was bright enough to also allow SWS observations around selected neon lines, to establish its excitation. While both PHOT-S spectra are characteristic of starburst-dominated galaxies, the neon line ratios in IRAS 21582+1018 indicate the presence of a hard excitation source. New, high-resolution, optical spectra show only a weak, broad component around Halpha, classifying now these two objects as Sey 1.9 galaxies. The two other galaxies observed are the NLS1 galaxies Mrk 359 and Mrk 1388. Their ISO spectra however do not reveal the typical, strong PAH features found in the starburst galaxies and are more like those of standard Seyferts. These results show therefore that, although IR observations were expected to be able to always reveal the presence of an active nucleus by piercing through the central obscuration, the result may be ambiguous: the broad band IR energy distribution can still be dominated by starburts located in a circumnuclear region, and the AGN appear only in specific observations (high-excitation lines in the IR, or optical spectra with better quality than classification spectra). The obscuration needs however to be patchy rather than complete, to explain the detection of the high-excitation lines or broad Balmer wings. Only high-energy observations can then establish the strength of the central AGN and the amount of extinction with certainty.

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 preliminary results of our H-band survey of the Chandra Deep Field South (CDFS). The observations were made using SofI on the NTT, and cover 0.027 square degrees with a 50% completeness limit of H = 20.5, and 0.17 square degrees with a 50% completeness limit of H =19.8. We used SExtractor to extract sources from our fields. In total we have detected 4819 objects. Star-galaxy separation was performed using the SExtractor parameter ``stellarity index''. All objects with an index of 0.5 or lower were classified as galaxies. According to this criterion, 80 % of our detections are galaxies. We then compare our results with previous observations of the CDFS. Our astrometric solutions are in good agreement with the Las Campanas Infrared Survey (LCIRS), the COMBO-17 and the ESO-EIS surveys. The photometry of our catalog compares satisfactorily with the results of the LCIRS, as well as with the GOODS data. Galaxy number counts are presented and compared with the LCIRS results. The present data are intended to complement the recent and future multi-wavelenghth observations of the CDFS and will be used, in conjuction with additional multiband photometry, to find counterparts of the upcoming mid-infrared surveys with SIRTF.

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

We present new near-infrared Keck and VLT spectroscopic data on the stellar dynamics in late stage, ultra-luminous infrared galaxy (ULIRG) mergers . We now have information on the structural and kinematic properties of 18 ULIRGs, 8 of which contain QSO-like active galactic nuclei. The host properties (velocity dispersion, effective radius, effective surface brightness, M_K) of AGN-dominated and star formation dominated ULIRGs are similar. ULIRGs fall remarkably close to the fundamental plane of early type galaxies. They populate a wide range of the plane, are on average similar to L*-rotating ellipticals, but are well offset from giant ellipticals and optically/UV bright, low-z QSOs/radio galaxies. ULIRGs and local QSOs/radio galaxies are very similar in their distributions of bolometric and extinction corrected near-IR luminosities, but ULIRGs have smaller effective radii and velocity dispersions than the local QSO/radio galaxy population. Hence, their host masses and inferred black hole masses are correspondingly smaller. The latter are more akin to those of local Seyfert galaxies. ULIRGs thus resemble local QSOs in their near-IR and bolometric luminosities because they are (much more) efficiently forming stars and/or feeding their black holes, and not because they have QSO-like, very massive black holes. We conclude that ULIRGs as a class cannot evolve into optically bright QSOs. They will more likely become quiescent, moderate mass field ellipticals or, when active, might resemble the X-ray bright, early type galaxies that have recently been found by the Chandra Observatory.

We present results from a survey of molecular hydrogen emission from a sample of Starburst and Seyfert galaxies carried out with the Infrared Space Observatory (ISO). Pure rotational H_2 emission has been detected in a number of extragalactic objects and a variety of environments. A number of transitions from S(7) to S(0) are detected in both Starbursts and Seyferts. Using excitation diagrams we derive temperatures and masses of the ``warm'' molecular hydrogen. We find that the temperature of the ``warm'' gas is similar in Starbursts and Seyferts (those Seyferts for which we have firm detections of the S(0) line) with a value of around T~150 K. This ``warm'' gas accounts for as much as 10% of the total galactic mass (as probed by CO molecular observations) in Starbursts. The fraction of ``warm'' gas is overall higher in Seyferts, ranging between 2--35%. We then investigate the origin of the warm H_2 emission. Comparison with published theoretical models and Galactic templates implies that although emission from photodissociation regions (PDR) alone could explain the emission from Starbursts and Seyferts, most likely a combination of PDR, shock emission and gas heated by X-rays (mostly for the Seyferts) is responsible for H_2 excitation in extragalactic environments. Finally, we find that although PAH and H_2 line emission correlate well in Starbursts and the large scale emission in AGN, H_2 emission is much stronger compared to PAH emission in cases where a ``pure'' AGN dominates the energy output.

With the advent of efficient near-IR spectrometers on 10m-class telescopes, exploiting the new generation of low readout noise, large format detectors, OH avoidance and sub-arcsecond seeing, 1-2.4 micron spectroscopy is becoming a key means of obtaining detailed galaxy dynamics and for studies of high-z galaxies. In the following we present the results of three recent IR spectroscopy studies on the dynamics of ULIRG mergers, super star clusters in the Antennae, and on the properties of the rotation curves of z~1 disk galaxies, carried out with ISAAC on the VLT and NIRSPEC on the Keck.

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

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

The IRAM Plateau de Bure mm interferometer and deep K-band imaging have been used to identify the brightest submm source detected in the Lockman field of the UK 8mJy SCUBA survey. The near infrared counterpart is an extended (20-30kpc), clumpy, and extremely red object. The spectral energy distribution suggests it to be a dusty star forming object at a redshift of about 3 (2-4). Its star formation rate and near-infrared properties are consistent with Lockman850.1 being a massive elliptical in formation.

We report high quality near-infrared spectroscopy of 12 ultra-luminous infrared galaxy mergers (ULIRGs). Our new VLT and Keck data provide ~0.5" resolution, stellar and gas kinematics of these galaxies most of which are compact systems in the last merger stages. We confirm that ULIRG mergers are 'ellipticals-in-formation'. Random motions dominate their stellar dynamics, but significant rotation is common. Gas and stellar dynamics are decoupled in most systems. ULIRGs fall on or near the fundamental plane of hot stellar systems, and especially on its less evolution sensitive, r(eff)-sigma projection. The ULIRG velocity dispersion distribution, their location in the fundamental plane and their distribution of v(rot)*sin(i)/sigma closely resemble those of intermediate mass (~L*), elliptical galaxies with moderate rotation. As a group ULIRGs do not resemble giant ellipticals with large cores and little rotation. Our results are in good agreement with other recent studies indicating that disky ellipticals with compact cores or cusps can form through dissipative mergers of gas rich, disk galaxies while giant ellipticals with large cores have a different formation history.

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

We present low resolution mid infrared spectra of 16 ultraluminous infrared galaxies (ULIRGs) obtained with the CVF spectroscopy mode of ISOCAM on board the Infrared Space Observatory ISO. Our sample completes previous ISO spectroscopy of ultra- and hyperluminous infrared galaxies towards higher luminosities. The combined samples cover an infrared luminosity range of \~10^{12 -13.1} Lo. For about half of the high luminosity ULIRGs studied here, strong aromatic emission bands suggest starburst dominance. Other spectra are dominated by a strong AGN-related continuum with weak superposed emission features of uncertain nature. An improved method to quantitatively characterize the relative contribution of star formation and AGN activity to the mid-infrared emission of ULIRGs is presented. As dominant source of the bolometric luminosity, starbursts prevail at the lower end and AGNs at the higher end of this range. The transition between mostly starburst and mostly AGN powered occurs at ~10^{12.4}$ to $10^{12.5} Lo, and individual luminous starbursts are found up to ~10^{12.65} Lo.

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

The prominent interacting galaxy pair NGC 4038/9 contains many active star-forming regions and is continuously forming new star clusters. We present a self-consistent n-body model for this system which includes an SPH gas component. The model qualitatively explains the apparent concentration of gas in the so-called overlap region between the two nuclei as a bridge of gas connecting the two galaxies. Projected on the sky, the bridge appears as a dense spot of gas. We discuss some implications for the evolution of Ultra-luminous infrared galaxies.

We present results and source counts at 90micron extracted from the Preliminary Analysis of the European Large Area ISO Survey (ELAIS). The survey covered about 11.6 square degrees of the sky in four main areas and was carried out with the PHOT instrument onboard the Infrared Space Observatory (ISO). The survey is at least an order of magnitude deeper than the IRAS 100micron survey and is expected to provide constraints on the formation and evolution of galaxies. The majority of the detected sources are associated with galaxies on optical images. In some cases the optical associations are interacting pairs or small groups of galaxies suggesting the sample may include a significant fraction of luminous infrared galaxies. The source counts extracted from a reliable subset of the detected sources are in agreement with strongly evolving models of the starburst galaxy population.

We report the results of near-infrared VLT-ISAAC spectroscopy of a sample of 12 galaxies at z = 0.4-1.4, drawn from the ISOCAM survey of the Hubble Deep Field South. We find that the rest frame R-band spectra of the ISOCAM galaxies resemble those of powerful dust-enshrouded starbursts. Halpha emission is detected in 11 out of 12 objects down to a flux limit of 7x10^(-17) erg/cm^2/s, corresponding to a luminosity limit of 10^41 erg/s at z = 0.6, (for an Ho = 50 and Omega = 0.3 cosmology). From the Halpha luminosities in these galaxies we derive estimates of the star formation rate in the range 2--50 Mo/yr for stellar masses 1--100 Mo. The raw Halpha-based star formation rates are an order of magnitude or more lower than SFR(FIR) estimates based on ISOCAM LW3 fluxes. If the Halpha emission is corrected for extinction the median offset is reduced to a factor of 3. The sample galaxies are part of a new population of optically faint but infrared--luminous active starburst galaxies, which are characterized by an extremely high rate of evolution with redshift up to z~1.5 and expected to contribute significantly to the cosmic far-IR extragalactic background.

A new population of faint galaxies characterized by an extremely high rate of evolution with redshift up to z~1.5 has recently been discovered by ISO. These sources are likely to contribute significantly to the cosmic far-IR extragalactic background. We have carried out near-infrared VLT-ISAAC spectroscopy of a sample of ISOCAM galaxies from the Hubble Deep Field South. The rest-frame R-band spectral properties of the ISO population resembles that of powerful dust-enshrouded active starburst galaxies.

We present a low resolution mid-infrared spectroscopic survey of an unbiased sample of 62 Ultraluminous Infrared Galaxies (ULIRGs) using ISOPHOT-S on board ISO. For comparison we also present ISOPHOT-S spectra for 23 AGN and 15 starburst and normal galaxies. The line-to-continuum ratio of the 7.7 $\mu$m PAH emission feature is used as a discriminator between starburst and AGN activity in ULIRGs. We find that the majority of ULIRGs are predominantly powered by starbursts. The ratio L$_{PAH}/$L$_{IR}$ for starburst dominated ULIRGs is similar to the one found for template starbursts.The shapes of the PAH features are sometimes unusual.Extinction has a noticeable effect on the PAH spectra of ULIRG starbursts. We have obtained high resolution near-infrared imaging for the Southern ISOPHOT-S ULIRGs to investigate their evolution stage. The majority (68%) of the ULIRGs imaged are double systems and all show distorted morphologies. Using separations for our work and from the literature, we have examined whether ULIRGs that are advanced mergers are more AGN-like. We have found no such evidence contrary to what is postulated by the classical evolutionary scenario. No correlation is found between the stage of merger in ULIRGs and their infrared luminosity.We also find that the total mass of interstellar gas, as estimated from the CO (1-0) luminosity, does not decrease with decreasing merger separation. When both an AGN and a starburst occur concurrently in ULIRGs, we find that the starburst dominates the luminosity output. We propose that the available gas reservoir and the individual structure of the interacting galaxies plays a major role in the evolution of the system.

The IRAM interferometer has been used to detect the submm source HDF 850.1 found by Hughes et al. (1998) in the Hubble Deep Field. The flux density measured at 1.3mm is 2.2 mJy, in agreement with the flux density measured at the JCMT. The flux densities and upper limits measured at 3.4, 2.8, 1.3, 0.85, and 0.45 mm show that the emission is from dust. We suggest that the 1.3mm dust source is associated with the optical arc-like feature, 3-593.0, that has a photometric redshift of about 1.7. If HDF 850.1 is at this redshift and unlensed, its spectral energy distribution, combined with that of 3-593.0, matches closely that of the ultraluminous galaxy VII Zw 31. Another possibility is that the dust source may be gravitationally lensed by the elliptical galaxy 3-586.0 at a redshift of 1.

We report the first results of a low resolution mid-infrared spectroscopic survey of an unbiased, far-infrared selected sample of 60 ultraluminous infrared galaxies, using ISOPHOT-S on board ISO. We use the ratio of the 7.7um `PAH' emission feature to the local continuum as a discriminator between starburst and AGN activity. About 80% of all the ULIRGs are found to be predominantly powered by star formation but the fraction of AGN powered objects increases with luminosity. Observed ratios of the PAH features in ULIRGs differ slightly from those in lower luminosity starbursts. This can be plausibly explained by the higher extinction and/or different physical conditions in the interstellar medium of ULIRGs. The PAH feature-to-continuum ratio is anticorrelated with the ratio of feature-free 5.9um continuum to the IRAS 60um continuum, confirming suggestions that strong mid-IR continuum is a prime AGN signature. The location of starburst-dominated ULIRGs in such a diagram is consistent with previous ISO-SWS spectroscopy which implies significant extinction even in the mid-infrared. We have searched for indications that ULIRGs which are advanced mergers might be more AGN-like, as postulated by the classical evolutionary scenario. No such trend has been found amongst those objects for which near infrared images are available to assess their likely merger status.

We have obtained ISOPHOT-S low resolution mid-infrared spectra of a sample of 60 Ultraluminous Infrared Galaxies (ULIRGs). We use the strength of the `PAH' mid-infrared features as a discriminator between starburst and AGN activity, and to probe for evolutionary effects. We focus on the fact that observed ratios of PAH features in ULIRGs differ slightly from those in lower luminosity starbursts. We suggest that such PAH ratio changes relate to the conditions in the interstellar medium in these galaxies, and in particular to extinction.

We present an ISO SWS and ISOPHOT-S, mid-infrared spectroscopic survey of 15 ultra-luminous IRAS galaxies. We combine the survey results with a detailed case study, based on near-IR and mm imaging spectroscopy, of one of the sample galaxies (UGC 5101). We compare the near- and mid-IR characteristics of these ultra-luminous galaxies to ISO and literature data of thirty starburst and active galactic nuclei (AGN), template galaxies. We find that 1) 70-80% of the ultra-luminous IRAS galaxies in our sample are predominantly powered by recently formed massive stars. 20-30% are powered by a central AGN. These conclusions are based on a new infrared 'diagnostic diagram' involving the ratio of high to low excitation mid-IR emission lines on the one hand, and on the strength of the 7.7um PAH feature on the other hand. 2) at least half of the sources probably have simultaneously an active nucleus and starburst activity in a 1-2 kpc diameter circum-nuclear disk/ring. 3) the mid-infrared emitting regions are highly obscured. After correction for these extinctions, we estimate that the star forming regions in ULIRGs have ages between 10^7 and 10^8 years, similar to but somewhat larger than those found in lower luminosity starburst galaxies. 4) in the sample we have studied there is no obvious trend for the AGN component to dominate in the most compact, and thus most advanced mergers. Instead, at any given time during the merger evolution, the time dependent compression of the circum-nuclear interstellar gas, the accretion rate onto the central black hole and the associated radiation efficiency may determine whether star formation or AGN activity dominates the luminosity of the system.