We report high quality, H{\alpha} or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs), across the peak of cosmic galaxy evolution (z~0.67-2.45), taken with the ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts with beam-convolved, forward modeling with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from abundance matching scaling relations. Two-thirds or more of the z>1.2 SFGs are baryon dominated within a few Re of typically 5.5 kpc, and have DM fractions less than maximal disks (<fDM (Re)>=0.12). At lower redshift (z<1.2) that fraction is less than one-third. DM fractions correlate inversely with the baryonic angular momentum parameter, baryonic surface density and bulge mass. Inferred low DM fractions cannot apply to the entire disk & halo but more plausibly reflect a flattened, or cored, inner DM density distribution. The typical central 'DM deficit' in these cores relative to NFW distributions is ~30% of the bulge mass. The observations are consistent with rapid radial transport of baryons in the first generation massive gas rich halos forming globally gravitationally unstable disks, and leading to efficient build-up of massive bulges and central black holes. A combination of heating due to dynamical friction and AGN feedback may drive DM out of the initial cusps.

We use deep adaptive optics assisted integral field spectroscopy from SINFONI on the VLT to study the spatially resolved properties of ionized gas outflows driven by active galactic nuclei (AGN) in three galaxies at z~2.2 -- K20-ID5, COS4-11337 and J0901+1814. These systems probe AGN feedback from nuclear to circumgalactic scales, and provide unique insights into the different mechanisms by which AGN-driven outflows interact with their host galaxies. K20-ID5 and COS4-11337 are compact star forming galaxies with powerful $\sim$1500 km s$^{-1}$ AGN-driven outflows that dominate their nuclear H$\alpha$ emission. The outflows do not appear to have any impact on the instantaneous star formation activity of the host galaxies, but they carry a significant amount of kinetic energy which could heat the halo gas and potentially lead to a reduction in the rate of cold gas accretion onto the galaxies. The outflow from COS4-11337 is propagating directly towards its companion galaxy COS4-11363, at a projected separation of 5.4 kpc. COS4-11363 shows signs of shock excitation and recent truncation of star formation activity, which could plausibly have been induced by the outflow from COS4-11337. J0901+1814 is gravitationally lensed, giving us a unique view of a compact (R = 470 $\pm$ 70 pc), relatively low velocity ($\sim$650 km s$^{-1}$) AGN-driven outflow. J0901+1814 has a similar AGN luminosity to COS4-11337, suggesting that the difference in outflow properties is not related to the current AGN luminosity, and may instead reflect a difference in the evolutionary stage of the outflow and/or the coupling efficiency between the AGN ionizing radiation field and the gas in the nuclear regions.

We report on the results of a search for serendipitous sources in CO emission in 110 cubes targeting CO(2-1), CO(3-2), and CO(6-5) at z ~ 1-2 from the second Plateau de Bure High-z Blue-Sequence Survey (PHIBSS2). The PHIBSS2 observations were part of a 4-year legacy program at the IRAM Plateau de Bure Interferometer aimed at studying early galaxy evolution from the perspective of molecular gas reservoirs. We present a catalog of 67 candidate secondary sources from this search, with 45 out of the 110 data cubes showing sources in addition to the primary target that appear to be field detections, unrelated to the central sources. This catalog includes the redshifts, line widths, fluxes, as well as an estimation of their reliability based on their false positive probability. We perform a search in the 3D-HST/CANDELS catalogs for the secondary CO detections and tentatively find that ~64% of these have optical counterparts, which we use to constrain their redshifts. Finally, we use our catalog of candidate CO detections to derive the CO(2-1), CO(3-2), CO(4-3), CO(5-4), and CO(6-5) luminosity functions over a range of redshifts, as well as the molecular gas mass density evolution. Despite the different methodology, these results are in very good agreement with previous observational constraints derived from blind searches in deep fields. They provide an example of the type of "deep field" science that can be carried out with targeted observations.

We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.

We present half-light sizes measured from H${\alpha}$ emission tracing star-formation in 281 star-forming galaxies from the KMOS3D survey at 0.7 < z < 2.7. Sizes are derived by fitting 2D exponential disk models, with bootstrap errors averaging 20%. H${\alpha}$ sizes are a median (mean) of 1.19 (1.26) times larger than those of the stellar continuum, which due to radial dust gradients places an upper limit on the growth in stellar size via star formation, with just 43% intrinsic scatter. At fixed continuum size the H${\alpha}$ size shows no residual trend with stellar mass, star formation rate, redshift or morphology. The only significant residual trend is with the excess obscuration of H${\alpha}$ by dust, at fixed continuum obscuration. The scatter in continuum size at fixed stellar mass is likely driven by the scatter in halo spin parameters. The stability of the ratio of H${\alpha}$ size to continuum size demonstrates a high degree of stability in halo spin and in the transfer of angular momentum to the disk over a wide range of physical conditions and cosmic time. This may require local regulation by feedback processes. The implication of our results, as we demonstrate using a toy model, is that our upper limit on star-formation driven growth is sufficient only to evolve star-forming galaxies approximately along the observed size-mass relation, consistent with the size growth of galaxies at constant cumulative co-moving number density. To explain the observed evolution of the size-mass relation of star-forming disk galaxies other processes, such as the preferential quenching of compact galaxies or galaxy mergers, may be required.

Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $\mu$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular gas in nuclear regions of galaxies. Our study contributes to double the number of OH observations of luminous systems at $z\sim0.2-0.3$, and push the search for molecular outflows based on the OH transition to $z\sim6$. Results. We detect OH high-velocity absorption wings in three of the four ULIRGs. In two cases, IRAS F20036-1547 and IRAS F13352+6402, the blueshifted absorption profiles indicate the presence of powerful and fast molecular gas outflows. Consistent with an inside-out quenching scenario, these outflows are depleting the central reservoir of molecular gas at a similar rate than the intense star formation activity. In the case of the starburst-dominated system IRAS 10091+4704, we detect an inverted P-Cygni profile that is unique among ULIRGs and indicates the presence of a fast ($\sim400$ km s$^{-1}$) inflow of molecular gas at a rate of $\sim100~M_{\odot}~{\rm yr}^{-1}$ towards the central region. Finally, we tentatively detect ($\sim3\sigma$) the OH doublet in absorption in the $z=6.13$ QSO ULAS J131911+095051. The OH feature is blueshifted with a median velocity that suggests the presence of a molecular outflow, although characterized by a modest molecular mass loss rate of $\sim200~M_{\odot}~{\rm yr}^{-1}$. This value is comparable to the small mass outflow rates found in the stacking of the [CII] spectra of other $z\sim6$ QSOs and suggests that ejective feedback in this phase of the evolution of ULAS J131911+095051 has subsided.

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.

We present VLT/SINFONI near-infrared (NIR) integral field spectroscopy of six $z \sim 0.2$ Lyman break galaxy "analogs" (LBAs), from which we detect HI, HeI, and [FeII] recombination lines, and multiple H$_2$ ro-vibrational lines in emission. Pa$\alpha$ kinematics reveal high velocity dispersions and low rotational velocities relative to random motions ($\langle v/\sigma \rangle = 1.2 \pm 0.8$). Matched-aperture comparisons of H$\beta$, H$\alpha$, and Pa$\alpha$ reveal that the nebular color excesses are lower relative to the continuum color excesses than is the case for typical local star-forming systems. We compare observed HeI/HI recombination line ratios to photoionization models to gauge the effective temperatures (T$_{\rm eff}$) of massive ionizing stars, finding the properties of at least one LBA are consistent with extra heating from an active galactic nucleus (AGN) and/or an overabundance of massive stars. We use H$_2$ 1-0 S($\cdot$) ro-vibrational spectra to determine rotational excitation temperature $T_{\rm ex} \sim 2000$ K for warm molecular gas, which we attribute to UV heating in dense photon-dominated regions. Spatially resolved NIR line ratios favor excitation by massive, young stars, rather than supernovae or AGN feedback. Our results suggest that the local analogs of Lyman break galaxies are primarily subject to strong feedback from recent star formation, with evidence for AGN and outflows in some cases.

To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately $\pm$1200 km s$^{-1}$ and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located $\sim$6 kpc north-west from the galaxy moving away at $\sim$900 km s$^{-1}$, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, warm and cold molecular gas phases is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for about $\sim$40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central $\sim$1.5 kpc region of only $\sim3$ Myr, a factor of $\sim2$ shorter than the depletion time by star formation activity.

We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged)

We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.

We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.

We present the completed KMOS$^\mathrm{3D}$ survey $-$ an integral field spectroscopic survey of 739, $\log(M_{\star}/M_{\odot})>9$, galaxies at $0.6<z<2.7$ using the K-band Multi Object Spectrograph (KMOS) at the Very Large Telescope (VLT). KMOS$^\mathrm{3D}$ provides a population-wide census of kinematics, star formation, outflows, and nebular gas conditions both on and off the star-forming galaxy main sequence through the spatially resolved and integrated properties of H$\alpha$, [N II], and [S II] emission lines. We detect H$\alpha$ emission for 91% of galaxies on the main sequence of star-formation and 79% overall. The depth of the survey has allowed us to detect galaxies with star-formation rates below 1 M$_{\odot}$/ yr$^{-1}$, as well as to resolve 81% of detected galaxies with $\geq3$ resolution elements along the kinematic major axis. The detection fraction of H$\alpha$ is a strong function of both color and offset from the main sequence, with the detected and non-detected samples exhibiting different SED shapes. Comparison of H$\alpha$ and UV+IR star formation rates (SFRs) reveal that dust attenuation corrections may be underestimated by 0.5 dex at the highest masses ($\log(M_{\star}/M_{\odot})>10.5$). We confirm our first year results of a high rotation dominated fraction (monotonic velocity gradient and $v_\mathrm{rot}$/$\sigma_0 > \sqrt{3.36}$) of 77% for the full KMOS$^\mathrm{3D}$ H$\alpha$sample. The rotation-dominated fraction is a function of both stellar mass and redshift with the strongest evolution measured over the redshift range of the survey for galaxies with $\log(M_{\star}/M_{\odot})<10.5$. With this paper we include a final data release of all 739 observed objects.

We present 1.4 pc resolution observations of 256 GHz nuclear radio continuum and HCN ($J=3 \to 2$) in the molecular torus of NGC 1068. The integrated radio continuum emission has a flat spectrum consistent with free-free emission and resolves into an X-shaped structure resembling an edge-brightened bicone. HCN is detected in absorption against the continuum, and the absorption spectrum shows a pronounced blue wing that suggests a high-velocity molecular outflow with speeds reaching 450 km/s. Analysis of the off-nucleus emission line kinematics and morphology reveals two nested, rotating disk components. The inner disk, inside $r\sim 1.2$ pc, has kinematics consistent with the nearly edge-on, geometrically thin water megamaser disk in Keplerian rotation around a central mass of $1.66\times 10^7\,\mbox{M}_\odot$. The outer disk, which extends to $\sim 7$~pc radius, counter-rotates relative to the inner disk. The rotation curve of the outer disk is consistent with rotation around the same central mass as the megamaser disk but in the opposite sense. The morphology of the molecular gas is asymmetric around the nuclear continuum source. We speculate that the outer disk formed from more recently introduced molecular gas falling out of the host galaxy or from a captured dwarf satellite galaxy. In NGC 1068, we find direct evidence that the molecular torus consists of counter-rotating and misaligned disks on parsec scales.

We present the $0.6<z<2.6$ evolution of the ionized gas velocity dispersion in 175 star-forming disk galaxies based on data from the full KMOS$^{\rm 3D}$ integral field spectroscopic survey. In a forward-modelling Bayesian framework including instrumental effects and beam-smearing, we fit simultaneously the observed galaxy velocity and velocity dispersion along the kinematic major axis to derive the intrinsic velocity dispersion $\sigma_0$. We find a reduction of the average intrinsic velocity dispersion of disk galaxies as a function of cosmic time, from $\sigma_0\sim45$ km s$^{-1}$ at $z\sim2.3$ to $\sigma_0\sim30$ km s$^{-1}$ at $z\sim0.9$. There is substantial intrinsic scatter ($\sigma_{\sigma_0, {\rm int}}\approx10$ km s$^{-1}$) around the best-fit $\sigma_0-z$-relation beyond what can be accounted for from the typical measurement uncertainties ($\delta\sigma_0\approx12$ km s$^{-1}$), independent of other identifiable galaxy parameters. This potentially suggests a dynamic mechanism such as minor mergers or variation in accretion being responsible for the scatter. Putting our data into the broader literature context, we find that ionized and atomic+molecular velocity dispersions evolve similarly with redshift, with the ionized gas dispersion being $\sim10-15$ km s$^{-1}$ higher on average. We investigate the physical driver of the on average elevated velocity dispersions at higher redshift, and find that our galaxies are at most marginally Toomre-stable, suggesting that their turbulent velocities are powered by gravitational instabilities, while stellar feedback as a driver alone is insufficient. This picture is supported through comparison with a state-of-the-art analytical model of galaxy evolution.

We present ~1" resolution (~2 kpc in the source plane) observations of the CO(1-0), CO(3-2), Halpha, and [N II] lines in the strongly-lensed z=2.26 star-forming galaxy SDSS J0901+1814. We use these observations to constrain the lensing potential of a foreground group of galaxies, and our source-plane reconstructions indicate that SDSS J0901+1814 is a nearly face-on (i~30 degrees) massive disk with r_{1/2}>~4 kpc for its molecular gas. Using our new magnification factors (mu_tot~30), we find that SDSS J0901+1814 has a star formation rate (SFR) of 268^{+63}_{-61} M_sun/yr, M_gas=(1.6^{+0.3}_{-0.2})x10^11x(alpha_CO/4.6) M_sun, and M_star=(9.5^{+3.8}_{-2.8})x10^10 M_sun, which places it on the star-forming galaxy "main sequence." We use our matched high-angular resolution gas and SFR tracers (CO and Halpha, respectively) to perform a spatially resolved (pixel-by-pixel) analysis of SDSS J0901+1814 in terms of the Schmidt-Kennicutt relation. After correcting for the large fraction of obscured star formation (SFR_Halpha/SFR_TIR=0.054^{+0.015}_{-0.014}), we find SDSS J0901+1814 is offset from "normal" star-forming galaxies to higher star formation efficiencies independent of assumptions for the CO-to-H_2 conversion factor. Our mean best-fit index for the Schmidt-Kennicutt relation for SDSS J0901+1814, evaluated with different CO lines and smoothing levels, is n=1.54+/-0.13; however, the index may be affected by gravitational lensing, and we find n=1.24+/-0.02 when analyzing the source-plane reconstructions. While the Schmidt-Kennicutt index largely appears unaffected by which of the two CO transitions we use to trace the molecular gas, the source-plane reconstructions and dynamical modeling suggest that the CO(1-0) emission is more spatially extended than the CO(3-2) emission.

We present the X-ray properties of 108 Dust-Obscured Galaxies (DOGs; F$_{24 \mu m}$/F$_{R} >$ 1000) in the COSMOS field, all of which detected in at least three far-infrared bands with the Herschel Observatory. Out of the entire sample, 22 are individually detected in the hard 2-8 keV X-ray band by the Chandra COSMOS Legacy survey, allowing us to classify them as AGN. Of them, 6 (27%) are Compton Thick AGN candidates with column densities N$_{H}$$>$10$^{24}$ cm$^{-2}$ while 15 are moderately obscured AGNs with 10$^{22}$ $<$ N$_{H}$ $<$ 10$^{24}$ cm$^{-2}$. Additionally, we estimate AGN contributions to the IR luminosity (8-1000$\mu$m rest-frame) greater than 20% for 19 DOGs based on SED decomposition using Spitzer/MIPS 24$\mu$m and the five Herschel bands (100-500 $\mu$m). Only 7 of these are detected in X-rays individually. We performed a X-ray stacking analysis for the 86 undetected DOGs. We find that the AGN fraction in DOGs increases with 24$\mu$m flux and that it is higher than that of the general 24$\mu$m population. However, no significant difference is found when considering only X-ray detections. This strongly motivates the combined use of X-ray and far-IR surveys to successfully probe a wider population of AGNs, particularly for the most obscured ones.

The impact of Active Galactic Nuclei (AGN) on star formation has implications for our understanding of the relationships between supermassive black holes and their galaxies, as well as for the growth of galaxies over the history of the Universe. We report on a high-resolution multi-phase study of the nuclear environment in the nearby Seyfert galaxy NGC 2110 using the Atacama Large Millimeter Array (ALMA), Hubble and Spitzer Space Telescopes, and the Very Large Telescope/SINFONI. We identify a region that is markedly weak in low-excitation CO $2\rightarrow1$ emission from cold molecular gas, but appears to be filled with ionised and warm molecular gas, which indicates that the AGN is directly influencing the properties of the molecular material. Using multiple molecular gas tracers, we demonstrate that, despite the lack of CO line emission, the surface densities and kinematics of molecular gas vary smoothly across the region. Our results demonstrate that the influence of an AGN on star-forming gas can be quite localized. In contrast to widely-held theoretical expectations, we find that molecular gas remains resilient to the glare of energetic AGN feedback.

We present a census of ionized gas outflows in 599 normal galaxies at redshift 0.6<z<2.7, mostly based on integral field spectroscopy of Ha, [NII], and [SII] line emission. The sample fairly homogeneously covers the main sequence of star-forming galaxies with masses 9.0<log(M*/Msun)<11.7, and probes into the regimes of quiescent galaxies and starburst outliers. About 1/3 exhibits the high-velocity component indicative of outflows, roughly equally split into winds driven by star formation (SF) and active galactic nuclei (AGN). The incidence of SF-driven winds correlates mainly with star formation properties. These outflows have typical velocities of ~450 km/s, local electron densities of n_e~380 cm^-3, modest mass loading factors of ~0.1-0.2 at all galaxy masses, and energetics compatible with momentum driving by young stellar populations. The SF-driven winds may escape from log(M*/Msun)<10.3 galaxies but substantial mass, momentum, and energy in hotter and colder outflow phases seem required to account for low galaxy formation efficiencies in the low-mass regime. Faster AGN-driven outflows (~1000-2000 km/s) are commonly detected above log(M*/Msun)~10.7, in up to ~75% of log(M*/Msun)>11.2 galaxies. The incidence, strength, and velocity of AGN-driven winds strongly correlates with stellar mass and central concentration. Their outflowing ionized gas appears denser (n_e~1000 cm^-3), and possibly compressed and shock-excited. These winds have comparable mass loading factors as the SF-driven winds but carry ~10 (~50) times more momentum (energy). The results confirm our previous findings of high duty cycle, energy-driven outflows powered by AGN above the Schechter mass, which may contribute to star formation quenching.

We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star forming galaxies at $z\sim$ 2-2.6, using adaptive optics assisted integral field observations from SINFONI on the VLT. The narrow and broad components of the H$\alpha$ emission are used to simultaneously determine the local star formation rate surface density ($\Sigma_{\rm SFR}$), and the outflow velocity $v_{\rm out}$ and mass outflow rate $\dot{M}_{\rm out}$, respectively. We find clear evidence for faster outflows with larger mass loading factors at higher $\Sigma_{\rm SFR}$. The outflow velocities scale as $v_{\rm out}$ $\propto$ $\Sigma_{\rm SFR}^{0.34 \pm 0.10}$, which suggests that the outflows may be driven by a combination of mechanical energy released by supernova explosions and stellar winds, as well as radiation pressure acting on dust grains. The majority of the outflowing material does not have sufficient velocity to escape from the galaxy halos, but will likely be re-accreted and contribute to the chemical enrichment of the galaxies. In the highest $\Sigma_{\rm SFR}$ regions the outflow component contains an average of $\sim$45% of the H$\alpha$ flux, while in the lower $\Sigma_{\rm SFR}$ regions only $\sim$10% of the H$\alpha$ flux is associated with outflows. The mass loading factor, $\eta$ = $\dot{M}_{\rm out}$/SFR, is positively correlated with $\Sigma_{\rm SFR}$ but is relatively low even at the highest $\Sigma_{\rm SFR}$: $\eta \lesssim$ 0.5 $\times$ (380 cm$^{-3}$/n$_e$). This may be in tension with the $\eta$ $\gtrsim$ 1 required by cosmological simulations, unless a significant fraction of the outflowing mass is in other gas phases and has sufficient velocity to escape the galaxy halos.

Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z=1-1.6) as well as its building-up (z=2-3) and winding-down (z=0.5-0.8) phases. The targets are drawn from the GOODS, COSMOS, and AEGIS deep fields and uniformly sample the MS in the stellar mass (M*) - star formation rate (SFR) plane with log(M*/Msun) = 10-11.8. We describe the survey strategy and sample selection before focusing on the results obtained at z=0.5-0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine their molecular gas masses and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio, gas fraction, and depletion time as well as their dependence with M* and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a wider range of redshifts. The galaxy-averaged Kennicutt-Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS.

The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars. Observations of broad emission lines could potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds. The size of this broad-line region has been estimated by measuring the light travel time delay between the variable nuclear continuum and the emission lines - a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow, whereas in others they are produced by orbiting gas clouds. Directly imaging such regions has not hitherto been possible because of their small angular sizes (< 0.1 milli-arcseconds). Here we report a spatial offset (with a spatial resolution of ten micro-arcseconds or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-{\alpha} line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 300 million solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time.

Nuclear outflows driven by accreting massive black holes are one of the main feedback mechanisms invoked at high-z to reproduce the distinct separation between star-forming, disk galaxies and quiescent spheroidal systems. Yet, our knowledge of feedback at high-z remains limited by the lack of observations of the multiple gas phases in galaxy outflows. In this work we use new deep, high-spatial resolution ALMA CO(3-2) and archival VLT/SINFONI H$\alpha$ observations to study the molecular and ionized components of the AGN-driven outflow in zC400528 ---a massive, main sequence galaxy at z=2.3 in the process of quenching. We detect a powerful molecular outflow that shows a positive velocity gradient and extends for at least ~10 kpc from the nuclear region, about three times the projected size of the ionized wind. The molecular gas in the outflow does not reach velocities high enough to escape the galaxy and is therefore expected to be reaccreted. Keeping in mind the various assumptions involved in the analysis, we find that the mass and energetics of the outflow are dominated by the molecular phase. The AGN-driven outflow in zC400528 is powerful enough to deplete the molecular gas reservoir on a timescale at least twice shorter than that needed to exhaust it by star formation. This suggests that the nuclear outflow is one of the main quenching engines at work in the observed suppression of the central star-formation activity in zC400528.

Gravitationally lensed systems allow a detailed view of galaxies at high redshift. High spatial- and spectral-resolution measurements of arc-like structures can offer unique constraints on the physical and dynamical properties of high-z systems. We present near-infrared spectra centred on the gravitational arcs of six known z ~ 2 lensed star-forming galaxies of stellar masses of 10^9-10^11 Msun and star formation rate (SFR) in the range between 10 and 400 Msun/yr. Ground layer adaptive optics (AO)-assisted observations are obtained at the Large Binocular Telescope (LBT) with the LUCI spectrographs during the commissioning of the ARGOS facility. We used MOS masks with curved slits to follow the extended arched structures and study the diagnostic emission lines. Combining spatially resolved kinematic properties across the arc-like morphologies, emission line diagnostics and archival information, we distinguish between merging and rotationally supported systems, and reveal the possible presence of ejected gas. For galaxies that have evidence for outflows, we derive outflow energetics and mass-loading factors compatible with those observed for stellar winds in local and high-z galaxies. We also use flux ratio diagnostics to derive gas-phase metallicities. The low signal-to-noise ratio in the faint H$\beta$ and nitrogen lines allows us to derive an upper limit of ~ 0.15 dex for the spatial variations in metallicity along the slit for the lensed galaxy J1038. Analysed near-infrared spectra presented here represent the first scientific demonstration of performing AO-assisted multi-object spectroscopy with narrow curved-shape slits. The increased angular and spectral resolution, combined with the binocular operation mode with the 8.4-m-wide eyes of LBT, will allow the characterisation of kinematic and chemical properties of a large sample of galaxies at high-z in the near future.

We present a cross-calibration of CO- and dust-based molecular gas masses at $z \leqslant 0.2$. Our results are based on a survey with the IRAM 30-m telescope collecting CO(1-0) measurements of 78 massive ($\log M_{\star} / M_{\odot} > 10$) galaxies with known gas-phase metallicities, and with IR photometric coverage from WISE (22 $\mu$m ) and Herschel SPIRE (250, 350, 500 $\mu$m). We find a tight relation ($\sim 0.17$ dex scatter) between the gas masses inferred from CO and dust continuum emission, with a minor systematic offset of 0.05 dex. The two methods can be brought into agreement by applying a metallicity-dependent adjustment factor ($\sim 0.13$ dex scatter). We illustrate that the observed offset is consistent with a scenario in which dust traces not only molecular gas, but also part of the ${\rm H \small I}$ reservoir, residing in the ${\rm H_2}$-dominated region of the galaxy. Observations of the CO(2-1) to CO(1-0) line ratio for two thirds of the sample indicate a narrow range in excitation properties, with a median ratio of luminosities $ \left\langle R_{21} \right\rangle \sim 0.64 $. Finally, we find dynamical mass constraints from spectral line profile fitting to agree well with the anticipated mass budget enclosed within an effective radius, once all mass components (stars, gas and dark matter) are accounted for.

Using integral field spectroscopy we investigate the kinematic properties of 35 massive centrally-dense and compact star-forming galaxies (${\log{\overline{M}_*}}=11.1$, $\log{(\Sigma_\mathrm{1kpc})}>9.5$, $\log{(M_\ast/r_e^{1.5})}>10.3$) at $z\sim0.7-3.7$ within the KMOS$^\mathrm{3D}$survey. We spatially resolve 23 compact star-forming galaxies (SFGs) and find that the majority are dominated by rotational motions with velocities ranging from {$95-500$ km s$^{-1}$}. The range of rotation velocities is reflected in a similar range of integrated H$\alpha$ linewidths, $75-400$ km s$^{-1}$, consistent with the kinematic properties of mass-matched extended galaxies from the full KMOS$^\mathrm{3D}$ sample. The fraction of compact SFGs that are classified as `rotation-dominated' or `disk-like' also mirrors the fractions of the full KMOS$^\mathrm{3D}$ sample. We show that integrated line-of-sight gas velocity dispersions from KMOS$^\mathrm{3D}$ are best approximated by a linear combination of their rotation and turbulent velocities with a lesser but still significant contribution from galactic scale winds. The H$\alpha$ exponential disk sizes of compact SFGs are on average $2.5\pm0.2$ kpc, $1-2\times$ the continuum sizes, in agreement with previous work. The compact SFGs have a $1.4\times$ higher AGN incidence than the full KMOS$^\mathrm{3D}$ sample at fixed stellar mass with average AGN fraction of 76\%. Given their high and centrally concentrated stellar masses as well as stellar to dynamical mass ratios close to unity, the compact SFGs are likely to have low molecular gas fractions and to quench on a short time scale unless replenished with inflowing gas. The rotation in these compact systems suggests that their direct descendants are rotating passive galaxies.

We use the Herschel/PACS spectrometer to study the global and spatially resolved far-infrared (FIR) fine-structure line emission in a sample of 52 galaxies that constitute the SHINING survey. These galaxies include star-forming, active-galactic nuclei (AGN), and luminous infrared galaxies (LIRGs). We find an increasing number of galaxies (and kiloparsec size regions within galaxies) with low line-to-FIR continuum ratios as a function of increasing FIR luminosity ($L_{\mathrm{FIR}}$), dust infrared color, $L_{\mathrm{FIR}}$ to molecular gas mass ratio ($L_{\mathrm{FIR}}/M_{\mathrm{mol}}$), and FIR surface brightness ($\Sigma_{\mathrm{FIR}}$). The correlations between the [CII]/FIR or [OI]/FIR ratios with $\Sigma_{\mathrm{FIR}}$ are remarkably tight ($\sim0.3$ dex scatter over almost four orders of magnitude in $\Sigma_{\mathrm{FIR}}$). We observe that galaxies with $L_{\mathrm{FIR}}/M_{\mathrm{mol}} \gtrsim 80\,L_{\odot}\,M_{\odot}^{-1}$ and $\Sigma_{\mathrm{FIR}}\gtrsim10^{11}$ $L_{\odot}$ kpc$^{-2}$ tend to have weak fine-structure line-to-FIR continuum ratios, and that LIRGs with infrared sizes $\gtrsim1$ kpc have line-to-FIR ratios comparable to those observed in typical star-forming galaxies. We analyze the physical mechanisms driving these trends in Paper II (Herrera-Camus et al. 2018). The combined analysis of the [CII], [NII], and [OIII] lines reveals that the fraction of the [CII] line emission that arises from neutral gas increases from 60% to 90% in the most active star-forming regions and that the emission originating in the ionized gas is associated with low-ionization, diffuse gas rather than with dense gas in HII regions. Finally, we report the global and spatially resolved line fluxes of the SHINING galaxies to enable the comparison and planning of future local and high-$z$ studies.

The SHINING survey (Paper I; Herrera-Camus et al. 2018) offers a great opportunity to study the properties of the ionized and neutral media of galaxies from prototypical starbursts and active galactic nuclei (AGN) to heavily obscured objects. Based on Herschel/PACS observations of the main far-infrared (FIR) fine-structure lines, in this paper we analyze the physical mechanisms behind the observed line deficits in galaxies, the apparent offset of luminous infrared galaxies (LIRGs) from the mass-metallicity relation, and the scaling relations between [CII] 158 $\mu$m line emission and star formation rate (SFR). Based on a toy model and the Cloudy code, we conclude that the increase in the ionization parameter with FIR surface brightness can explain the observed decrease in the line-to-FIR continuum ratio of galaxies. In the case of the [CII] line, the increase in the ionization parameter is accompanied by a reduction in the photoelectric heating efficiency and the inability of the line to track the increase in the FUV radiation field as galaxies become more compact and luminous. In the central $\sim$kiloparsec regions of AGN galaxies we observe a significant increase in the [OI] 63 $\mu$m/[CII] line ratio; the AGN impact on the line-to-FIR ratios fades on global scales. Based on extinction-insensitive metallicity measurements of LIRGs we confirm that they lie below the mass-metallicity relation, but the offset is smaller than those reported in studies that use optical-based metal abundances. Finally, we present scaling relations between [CII] emission and SFR in the context of the main-sequence of star-forming galaxies.

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged]

We present deep observations of a $z=1.4$ massive, star-forming galaxy in molecular and ionized gas at comparable spatial resolution (CO 3-2, NOEMA; H$\alpha$, LBT). The kinematic tracers agree well, indicating that both gas phases are subject to the same gravitational potential and physical processes affecting the gas dynamics. We combine the one-dimensional velocity and velocity dispersion profiles in CO and H$\alpha$ to forward-model the galaxy in a Bayesian framework, combining a thick exponential disk, a bulge, and a dark matter halo. We determine the dynamical support due to baryons and dark matter, and find a dark matter fraction within one effective radius of $f_{\rm DM}(\leq$$R_{e})=0.18^{+0.06}_{-0.04}$. Our result strengthens the evidence for strong baryon-dominance on galactic scales of massive $z\sim1-3$ star-forming galaxies recently found based on ionized gas kinematics alone.

We imaged with ALMA and ARGOS/LUCI the molecular gas and the dust and stellar continuum in XID2028, an obscured QSO at z=1.593, where the presence of a massive outflow in the ionized gas component traced by the [O III]5007 emission has been resolved up to 10 kpc. This target represents a unique test case to study QSO 'feedback in action' at the peak epoch of AGN-galaxy coevolution. The QSO has been detected in the CO(5-4) transition and in the 1.3mm continuum, at ~30 and ~20 {\sigma} significance respectively, with both emissions confined in the central (<4 kpc) radius area. Our analysis suggests the presence of a fast rotating molecular disc (v~400 km/s) on very compact scales, and well inside the galaxy extent seen in the rest-frame optical light (~10 kpc, as inferred from the LUCI data). Adding available measurements in additional two CO transitions, CO(2-1) and CO(3-2), we could derive a total gas mass of ~10$^{10}$ M$_\odot$, thanks to a critical assessment of CO excitation and the comparison with Rayleigh-Jeans continuum estimate. This translates into a very low gas fraction (<5%) and depletion time scales of 40-75 Myr, reinforcing the result of atypical gas consumption conditions in XID2028, possibly due to feedback effects on the host galaxy. Finally, we also detect at ~5{\sigma} the presence of high velocity CO gas, which we interpret as a signature of galaxy-scale molecular outflow, spatially coincident with the ionised gas outflow. XID2028 represents therefore a unique case where the measurement of total outflowing mass (~500-800 M$_\odot$/yr) including the molecular and atomic components, in both the ionised and neutral phases, has been attempted for a high-z QSO.

This paper provides an update of our previous scaling relations (Genzel et al.2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M_solar)=9.0-11.8, and star formation rates relative to that on the MS, delta_MS=SFR/SFR(MS), from 10^{-1.3} to 10^{2.2}. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t_depl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z)^{-0.6}x(delta_MS)^{-0.44}, and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass mu_gas depends on (1+z)^{2.5}x (delta_MS)^{0.52}x(M*)^{-0.36}, which tracks the evolution of the specific star formation rate. The redshift dependence of mu_gas requires a curvature term, as may the mass-dependences of t_depl and mu_gas. We find no or only weak correlations of t_depl and mu_gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.

A randomly chosen star in today's Universe is most likely to live in a galaxy with a stellar mass between that of the Milky Way and Andromeda. Yet it remains uncertain how the structural evolution of these bulge-disk systems proceeded. Most of the unobscured star formation we observe building Andromdeda progenitors at 0.7<z<1.5 occurs in disks, but >90% of their star formation is reprocessed by dust and remains unaccounted for. Here we map 500micron dust continuum emission in an Andromeda progenitor at z=1.25 to probe where it is growing through dust-obscured star formation. Combining resolved dust measurements from the NOEMA interferometer with Hubble Space Telescope Halpha maps and multicolor imaging (including new UV data from the HDUV survey), we find a bulge growing by dust-obscured star formation: while the unobscured star formation is centrally suppressed, the dust continuum is centrally concentrated, filling in the ring-like structures evident in the Halpha and UV emission. Reflecting this, the dust emission is more compact than the optical/UV tracers of star formation with r_e(dust)=3.4kpc, r_e(Halpha)/r_e(dust)=1.4, and r_e(UV)/r_e(dust)=1.8. Crucially, however, the bulge and disk of this galaxy are building simultaneously; although the dust emission is more compact than the rest-optical emission (r_e(optical)/r_e(dust)=1.4), it is somewhat less compact than the stellar mass (r_e(M_*)/r_e(dust)=0.9). Taking the 500micron emission as a tracer of star formation, the expected structural evolution of this galaxy can be accounted for by star formation: it will grow in size by Delta(r_e)/Delta(M_*)~0.3 and central surface density by Delta(Sigma_cen)/Delta(M_*)~0.9. Finally, our observations are consistent with a picture in which merging and disk instabilities drive gas to the center of galaxies, boosting global star formation rates above the main sequence and building bulges.

We explore the relationship between X-ray absorption and optical obscuration within the BAT AGN Spectroscopic Survey (BASS) which has been collecting and analyzing the optical and X-ray spectra for 641 hard X-ray selected ($E>14$ keV) active galactic nuclei (AGN). We use the deviation from a linear broad H$\alpha$-to-X-ray relationship as an estimate of the maximum optical obscuration towards the broad line region and compare the $A_{\rm V}$ to the hydrogen column densities ($N_{\rm H}$) found through systematic modeling of their X-ray spectra. We find that the inferred columns implied by $A_{\rm V}$ towards the broad line region (BLR) are often orders of magnitude less than the columns measured towards the X-ray emitting region indicating a small scale origin for the X-ray absorbing gas. After removing 30\% of Sy 1.9s that potentially have been misclassified due to outflows, we find that 86\% (164/190) of the Type 1 population (Sy 1--1.9) are X-ray unabsorbed as expected based on a single obscuring structure. However, 14\% (26/190), of which 70\% (18/26) are classified as Sy 1.9, are X-ray absorbed, suggesting the broad line region itself is providing extra obscuration towards the X-ray corona. The fraction of X-ray absorbed Type 1 AGN remains relatively constant with AGN luminosity and Eddington ratio, indicating a stable broad line region covering fraction.

Using new APEX and JCMT spectroscopy of the CO 2-1 line, we undertake a controlled study of cold molecular gas in moderately luminous Active Galactic Nuclei (AGN) and inactive galaxies from the Luminous Local AGN with Matched Analogs (LLAMA) survey. We use spatially resolved infrared photometry of the LLAMA galaxies from 2MASS, WISE, IRAS & Herschel, corrected for nuclear emission using multi-component spectral energy distribution (SED) fits, to examine the dust-reprocessed star-formation rates (SFRs), molecular gas fractions and star formation efficiencies (SFEs) over their central 1 - 3 kpc. We find that the gas fractions and central SFEs of both active and inactive galaxies are similar when controlling for host stellar mass and morphology (Hubble type). The equivalent central molecular gas depletion times are consistent with the discs of normal spiral galaxies in the local Universe. Despite energetic arguments that the AGN in LLAMA should be capable of disrupting the observable cold molecular gas in their central environments, our results indicate that nuclear radiation only couples weakly with this phase. We find a mild preference for obscured AGN to contain higher amounts of central molecular gas, which suggests a connection between AGN obscuration and the gaseous environment of the nucleus. Systems with depressed SFEs are not found among the LLAMA AGN. We speculate that the processes that sustain the collapse of molecular gas into dense pre-stellar cores may also be a prerequisite for the inflow of material on to AGN accretion disks.

In a complete sample of local 14-195 keV selected AGNs and inactive galaxies, matched by their host galaxy properties, we study the spatially resolved stellar kinematics and luminosity distributions at near-infrared wavelengths on scales of 10-150 pc, using SINFONI on the VLT. In this paper, we present the first half of the sample, which comprises 13 galaxies, 8 AGNs and 5 inactive galaxies. The stellar velocity fields show a disk-like rotating pattern, for which the kinematic position angle is in agreement with the photometric position angle obtained from large scale images. For this set of galaxies, the stellar surface brightness of the inactive galaxy sample is generally comparable to the matched sample of AGN but extends to lower surface brightness. After removal of the bulge contribution, we find a nuclear stellar light excess with an extended nuclear disk structure, and which exhibits a size-luminosity relation. While we expect the excess luminosity to be associated with a dynamically cooler young stellar population, we do not typically see a matching drop in dispersion. This may be because these galaxies have pseudo-bulges in which the intrinsic dispersion increases towards the centre. And although the young stars may have an impact in the observed kinematics, their fraction is too small to dominate over the bulge and compensate the increase in dispersion at small radii, so no dispersion drop is seen. Finally, we find no evidence for a difference in the stellar kinematics and nuclear stellar luminosity excess between these active and inactive galaxies.

We use Herschel data to analyze the size of the far-infrared 70micron emission for z<0.06 local samples of 277 hosts of Swift-BAT selected active galactic nuclei (AGN), and 515 comparison galaxies that are not detected by BAT. For modest far-infrared luminosities 8.5<log(LFIR)<10.5, we find large scatter of half light radii Re70 for both populations, but a typical Re70 <~ 1 kpc for the BAT hosts that is only half that of comparison galaxies of same far-infrared luminosity. The result mostly reflects a more compact distribution of star formation (and hence gas) in the AGN hosts, but compact AGN heated dust may contribute in some extremely AGN-dominated systems. Our findings are in support of an AGN-host coevolution where accretion onto the central black hole and star formation are fed from the same gas reservoir, with more efficient black hole feeding if that reservoir is more concentrated. The significant scatter in the far-infrared sizes emphasizes that we are mostly probing spatial scales much larger than those of actual accretion, and that rapid accretion variations can smear the distinction between the AGN and comparison categories. Large samples are hence needed to detect structural differences that favour feeding of the black hole. No size difference AGN host vs. comparison galaxies is observed at higher far-infrared luminosities log(LFIR)>10.5 (star formation rates >~ 6 Msun/yr), possibly because these are typically reached in more compact regions in the first place.

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

We investigate the stellar mass and baryonic mass Tully-Fisher relations (TFRs) of massive star-forming disk galaxies at redshift z~2.3 and z~0.9 as part of the KMOS^3D integral field spectroscopy survey. Our spatially resolved data allow reliable modelling of individual galaxies, including the effect of pressure support on the inferred gravitational potential. At fixed circular velocity, we find higher baryonic masses and similar stellar masses at z~2.3 as compared to z~0.9. Together with the decreasing gas-to-stellar mass ratios with decreasing redshift, this implies that the contribution of dark matter to the dynamical mass at the galaxy scale increases towards lower redshift. A comparison to local relations reveals a negative evolution of the stellar and baryonic TFR zero-points from z=0 to z~0.9, no evolution of the stellar TFR zero-point from z~0.9 to z~2.3, and a positive evolution of the baryonic TFR zero-point from z~0.9 to z~2.3. We discuss a toy model of disk galaxy evolution to explain the observed, non-monotonic TFR evolution, taking into account the empirically motivated redshift dependencies of galactic gas fractions, and of the relative amount of baryons to dark matter on galaxy and halo scales.

We present a new analysis of the PG quasar sample based on Spitzer and Herschel observations. (I) Assuming PAH-based star formation luminosities (L_SF) similar to Symeonidis et al. (2016, S16), we find mean and median intrinsic AGN spectral energy distributions (SEDs). These, in the FIR, appear hotter and significantly less luminous than the S16 mean intrinsic AGN SED. The differences are mostly due to our normalization of the individual SEDs, that properly accounts for a small number of very FIR-luminous quasars. Our median, PAH-based SED represents ~ 6% increase on the 1-243 micron luminosity of the extended Mor & Netzer (2012, EM12) torus SED, while S16 find a significantly larger difference. It requires large-scale dust with T ~ 20 -- 30 K which, if optically thin and heated by the AGN, would be outside the host galaxy. (II) We also explore the black hole and stellar mass growths, using L_SF estimates from fitting Herschel/PACS observations after subtracting the EM12 torus contribution. We use rough estimates of stellar mass, based on scaling relations, to divide our sample into groups: on, below and above the star formation main sequence (SFMS). Objects on the SFMS show a strong correlation between star formation luminosity and AGN bolometric luminosity, with a logarithmic slope of ~ 0.7. Finally we derive the relative duty cycles of this and another sample of very luminous AGN at z = 2 -- 3.5. Large differences in this quantity indicate different evolutionary pathways for these two populations characterised by significantly different black hole masses.

We explore the H-alpha emission in the massive quiescent galaxies observed by the KMOS-3D survey at 0.7 < z < 2.7. The H-alpha line is robustly detected in 20 out of 120 UVJ-selected quiescent galaxies, and we classify the emission mechanism using the H-alpha line width and the [NII]/H-alpha line ratio. We find that AGN are likely to be responsible for the line emission in more than half of the cases. We also find robust evidence for star formation activity in nine quiescent galaxies, which we explore in detail. The H-alpha kinematics reveal rotating disks in five of the nine galaxies. The dust-corrected H-alpha star formation rates are low (0.2 - 7 Msun/yr), and place these systems significantly below the main sequence. The 24micron-based infrared luminosities, instead, overestimate the star formation rates. These galaxies present a lower gas-phase metallicity compared to star-forming objects with similar stellar mass, and many of them have close companions. We therefore conclude that the low-level star formation activity in these nine quiescent galaxies is likely to be fueled by inflowing gas or minor mergers, and could be a sign of rejuvenation events.

We present spatially resolved ALMA observations of the CO J=3-2 emission line in two massive galaxies at z=2.5 on the star-forming main sequence. Both galaxies have compact dusty star-forming cores with effective radii of Re=1.3 kpc and Re=1.2 kpc in the 870 um continuum emission. The spatial extent of star-forming molecular gas is also compact with Re=1.9 kpc and Re=2.3 kpc, but more extended than the dust emission. Interpreting the observed position-velocity diagrams with dynamical models, we find the starburst cores to be rotation-dominated with the ratio of the maximum rotation velocity to the local velocity dispersion of v/sigma=7.0 (v=386 km/s) and v/sigma_0=4.1 (v=391 km/s). Given that the descendants of these massive galaxies in the local universe are likely ellipticals with v/sigma nearly an order of magnitude lower, the rapidly rotating galaxies would lose significant net angular momentum in the intervening time. The comparisons among dynamical, stellar, gas, and dust mass suggest that the starburst CO-to-H2 conversion factor of alpha_CO=0.8 Msun/(K km/s/pc2) is appropriate in the spatially resolved cores. The dense cores are likely to be formed in extreme environments similar to the central regions of local ultraluminous infrared galaxies. Our work also demonstrates that a combination of medium-resolution CO and high-resolution dust continuum observations is a powerful tool for characterizing the dynamical state of molecular gas in distant galaxies.

We present spectroscopic redshifts of S(870)>2mJy submillimetre galaxies (SMGs) which have been identified from the ALMA follow-up observations of 870um detected sources in the Extended Chandra Deep Field South (the ALMA-LESS survey). We derive spectroscopic redshifts for 52 SMGs, with a median of z=2.4+/-0.1. However, the distribution features a high redshift tail, with ~25% of the SMGs at z>3. Spectral diagnostics suggest that the SMGs are young starbursts, and the velocity offsets between the nebular emission and UV ISM absorption lines suggest that many are driving winds, with velocity offsets up to 2000km/s. Using the spectroscopic redshifts and the extensive UV-to-radio photometry in this field, we produce optimised spectral energy distributions (SEDs) using Magphys, and use the SEDs to infer a median stellar mass of M*=(6+/-1)x10^{10}Msol for our SMGs with spectroscopic redshifts. By combining these stellar masses with the star-formation rates (measured from the far-infrared SEDs), we show that SMGs (on average) lie a factor ~5 above the main-sequence at z~2. We provide this library of 52 template fits with robust and well-sampled SEDs available as a resource for future studies of SMGs, and also release the spectroscopic catalog of ~2000 (mostly infrared-selected) galaxies targeted as part of the spectroscopic campaign.

We exploit the deep resolved Halpha kinematic data from the KMOS^3D and SINS/zC-SINF surveys to examine the largely unexplored outer disk kinematics of star-forming galaxies (SFGs) out to the peak of cosmic star formation. Our sample contains 101 SFGs representative of the more massive (9.3 < log(M*/Msun) < 11.5) main sequence population at 0.6<z<2.6. Through a novel stacking approach we are able to constrain a representative rotation curve extending out to ~4 effective radii. This average rotation curve exhibits a significant drop in rotation velocity beyond the turnover, with a slope of Delta(V)/Delta(R) = $-0.26^{+0.10}_{-0.09}$ in units of normalized coordinates V/V_max and R/R_turn. This result confirms that the fall-off seen previously in some individual galaxies is a common feature of our sample of high-z disks. We show that this outer fall-off strikingly deviates from the flat or mildly rising rotation curves of local spiral galaxies of similar masses. We furthermore compare our data with models including baryons and dark matter demonstrating that the falling stacked rotation curve can be explained by a high mass fraction of baryons relative to the total dark matter halo (m_d>~0.05) in combination with a sizeable level of pressure support in the outer disk. These findings are in agreement with recent studies demonstrating that star-forming disks at high redshift are strongly baryon dominated within the disk scale, and furthermore suggest that pressure gradients caused by large turbulent gas motions are present even in their outer disks. We demonstrate that these results are largely independent of our model assumptions such as the presence of a central stellar bulge, the effect of adiabatic contraction at fixed m_d, and variations in the concentration parameter.

In cold dark matter cosmology, the baryonic components of galaxies are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the total mass of the galaxy and its dark matter halo. In the local Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius. Comparison between the dynamical mass and the sum of stellar and cold gas mass at the peak epoch of galaxy formation, inferred from ancillary data, suggest high baryon factions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (stellar initial mass function, calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of two main factors: first, a large fraction of the massive, high-redshift galaxy population was strongly baryon dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early Universe efficiently condensed at the centres of dark matter halos when gas fractions were high, and dark matter was less concentrated. [Abridged]

We discuss the environment of local hard X-ray selected active galaxies, with reference to two independent group catalogues. We find that the fraction of these AGN in S0 host galaxies decreases strongly as a function of galaxy group size (halo mass) - which contrasts with the increasing fraction of galaxies of S0 type in denser environments. However, there is no evidence for an environmental dependence of AGN in spiral galaxies. Because most AGN are found in spiral galaxies, this dilutes the signature of environmental dependence for the population as a whole. We argue that the differing results for AGN in disk-dominated and bulge-dominated galaxies is related to the source of the gas fuelling the AGN, and so may also impact the luminosity function, duty cycle, and obscuration. We find that there is a significant difference in the luminosity function for AGN in spiral and S0 galaxies, and tentative evidence for some difference in the fraction of obscured AGN.

We present 0".2-resolution Atacama Large Millimeter/submillimeter Array observations at 870 um for 25 Halpha-seleced star-forming galaxies (SFGs) around the main-sequence at z=2.2-2.5. We detect significant 870 um continuum emission in 16 (64%) of these SFGs. The high-resolution maps reveal that the dust emission is mostly radiated from a single region close to the galaxy center. Exploiting the visibility data taken over a wide $uv$ distance range, we measure the half-light radii of the rest-frame far-infrared emission for the best sample of 12 massive galaxies with logM*>11. We find nine galaxies to be associated with extremely compact dust emission with R_{1/2,870um}<1.5 kpc, which is more than a factor of 2 smaller than their rest-optical sizes, R_{1/2,1.6um}=3.2 kpc, and is comparable with optical sizes of massive quiescent galaxies at similar redshifts. As they have an exponential disk with Sersic index of n=1.2 in the rest-optical, they are likely to be in the transition phase from extended disks to compact spheroids. Given their high star formation rate surface densities within the central 1 kpc of Sigma SFR1kpc=40 Msol/yr/kpc^2, the intense circumnuclear starbursts can rapidly build up a central bulge with Sigma M*1kpc>1e10 Msol/kpc^2 in several hundred Myr, i.e. by z~2. Moreover, ionized gas kinematics reveal that they are rotation-supported with an angular momentum as large as that of typical SFGs at z=1-3. Our results suggest bulges are commonly formed in extended rotating disks by internal processes, not involving major mergers.

We exploit deep integral-field spectroscopic observations with KMOS/VLT of 240 star-forming disks at 0.6 < z < 2.6 to dynamically constrain their mass budget. Our sample consists of massive ($\gtrsim 10^{9.8} M_\odot$) galaxies with sizes $R_e \gtrsim 2$ kpc. By contrasting the observed velocity and dispersion profiles to dynamical models, we find that on average the stellar content contributes $32^{+8}_{-7}\%$ of the total dynamical mass, with a significant spread among galaxies (68th percentile range f_star ~ 18 - 62%). Including molecular gas as inferred from CO- and dust-based scaling relations, the estimated baryonic mass adds up to $56^{+17}_{-12}\%$ of total for the typical galaxy in our sample, reaching ~ 90% at z > 2. We conclude that baryons make up most of the mass within the disk regions of high-redshift star-forming disk galaxies, with typical disks at z > 2 being strongly baryon-dominated within $R_e$. Substantial object-to-object variations in both stellar and baryonic mass fractions are observed among the galaxies in our sample, larger than what can be accounted for by the formal uncertainties in their respective measurements. In both cases, the mass fractions correlate most strongly with measures of surface density. High $\Sigma_{star}$ galaxies feature stellar mass fractions closer to unity, and systems with high inferred gas or baryonic surface densities leave less room for additional mass components other than stars and molecular gas. Our findings can be interpreted as more extended disks probing further (and more compact disks probing less far) into the dark matter halos that host them.

We have detected in ALMA observations CO J = 6 - 5 emission from the nucleus of the Seyfert galaxy NGC 1068. The low-velocity (up to +/- 70 km/s relative to systemic) CO emission resolves into a 12x7 pc structure, roughly aligned with the nuclear radio source. Higher-velocity emission (up to +/- 400 km/s) is consistent with a bipolar outflow in a direction nearly perpendicular (roughly 80 degrees) to the nuclear disk. The position-velocity diagram shows that in addition to the outflow, the velocity field may also contain rotation about the disk axis. These observations provide compelling evidence in support of the disk-wind scenario for the AGN obscuring torus.

We use high spectral resolution (R > 8000) data covering 3800-13000\r{A} to study the physical conditions of the broad line region (BLR) of nine nearby Seyfert 1 galaxies. Up to six broad HI lines are present in each spectrum. A comparison - for the first time using simultaneous optical to near-infrared observations - to photoionisation calculations with our devised simple scheme yields the extinction to the BLR at the same time as determining the density and photon flux, and hence distance from the nucleus, of the emitting gas. This points to a typical density for the HI emitting gas of 10$^{11}$cm$^{-3}$ and shows that a significant amount of this gas lies at regions near the dust sublimation radius, consistent with theoretical predictions. We also confirm that in many objects the line ratios are far from case B, the best-fit intrinsic broad-line H$\alpha$/H$\beta$ ratios being in the range 2.5-6.6 as derived with our photoionization modeling scheme. The extinction to the BLR, based on independent estimates from HI and HeII lines, is A$_V$ $\le$ 3 for Seyfert 1-1.5s, while Seyfert 1.8-1.9s have A$_V$ in the range 4-8. A comparison of the extinction towards the BLR and narrow line region (NLR) indicates that the structure obscuring the BLR exists on scales smaller than the NLR. This could be the dusty torus, but dusty nuclear spirals or filaments could also be responsible. The ratios between the X-ray absorbing column N$_H$ and the extinction to the BLR are consistent with the Galactic gas-to-dust ratio if N$_H$ variations are considered.

We analyze the angular momenta of massive star forming galaxies (SFGs) at the peak of the cosmic star formation epoch (z~0.8-2.6). Our sample of ~360 log(M*/Msun) ~ 9.3-11.8 SFGs is mainly based on the KMOS3D and SINS/zC-SINF surveys of H$\alpha$ kinematics, and collectively provides a representative subset of the massive star forming population. The inferred halo scale angular momentum distribution is broadly consistent with that theoretically predicted for their dark matter halos, in terms of mean spin parameter <$\lambda$> ~ 0.037 and its dispersion ($\sigma_{log(\lambda)}$~0.2). Spin parameters correlate with the disk radial scale, and with their stellar surface density, but do not depend significantly on halo mass, stellar mass, or redshift. Our data thus support the long-standing assumption that on average, even at high redshifts, the specific angular momentum of disk galaxies reflects that of their dark matter halos (j_d = j_DM). The lack of correlation between $\lambda$ x (j_d/j_DM) and the nuclear stellar density $\Sigma_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to "compaction" inside massive high-redshift disks. For our sample, the inferred average stellar-to-dark matter mass ratio is ~2%, consistent with abundance matching results. Including the molecular gas, the total baryonic disk-to-dark matter mass ratio is ~5% for halos near $10^{12}$ Msun, which corresponds to 31% of the cosmologically available baryons, implying that high-redshift disks are strongly baryon dominated.

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

We present a complete census of all Herschel-detected sources within the six massive lensing clusters of the HST Frontier Fields (HFF). We provide a robust legacy catalogue of 263 sources with Herschel fluxes, primarily based on imaging from the Herschel Lensing Survey (HLS) and PEP/HerMES Key Programmes. We optimally combine Herschel, Spitzer and WISE infrared (IR) photometry with data from HST, VLA and ground-based observatories, identifying counterparts to gain source redshifts. For each Herschel-detected source we also present magnification factor (mu), intrinsic IR luminosity and characteristic dust temperature, providing a comprehensive view of dust-obscured star formation within the HFF. We demonstrate the utility of our catalogues through an exploratory overview of the magnified population, including more than 20 background sub-LIRGs unreachable by Herschel without the assistance gravitational lensing.

We report a tentative correlation between the outflow characteristics derived from OH absorption at $119\,\mu\text{m}$ and [CII] emission at $158\,\mu\text{m}$ in a sample of 22 local and bright ultraluminous infrared galaxies (ULIRGs). For this sample we investigate whether [CII] broad wings are a good tracer of molecular outflows, and how the two tracers are connected. Fourteen objects in our sample have a broad wing component as traced by [CII], and all of these also show OH119 absorption indicative of an outflow (in 1 case an inflow). The other eight cases, where no broad [CII] component was found, are predominantly objects with no OH outflow or a low-velocity ($\leq 100\,\text{km s}^{-1}$) OH outflow. The full width at half maximum (FWHM) of the broad [CII] component shows a trend with the OH119 blue-shifted velocity, although with significant scatter. Moreover, and despite large uncertainties, the outflow masses derived from OH and broad [CII] show a 1:1 relation. The main conclusion is therefore that broad [CII] wings can be used to trace molecular outflows. This may be particularly relevant at high redshift, where the usual tracers of molecular gas (like low-J CO lines) become hard to observe. Additionally, observations of blue-shifted Na I D $\lambda\lambda 5890,5896$ absorption are available for ten of our sources. Outflow velocities of Na I D show a trend with OH velocity and broad [CII] FWHM. These observations suggest that the atomic and molecular gas phases of the outflow are connected.

We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

We discuss the dense molecular gas in central regions of nearby Seyfert galaxies, and report new arcsec resolution observations of HCN(1-0) and HCO$^+$(1-0) for 3 objects. In NGC 3079 the lines show complex profiles as a result of self-absorption and saturated continuum absorption. H$^{13}$CN reveals the continuum absorption profile, with a peak close to the galaxy's systemic velocity that traces disk rotation, and a second feature with a blue wing extending to $-350$km s$^{-1}$ that most likely traces a nuclear outflow. The morphological and spectral properties of the emission lines allow us to constrain the dense gas dynamics. We combine our kinematic analysis for these 3 objects, as well as another with archival data, with a previous comparable analysis of 4 other objects, to create a sample of 8 Seyferts. In 7 of these, the emission line kinematics imply thick disk structures on radial scales of $\sim$100pc, suggesting such structures are a common occurrence. We find a relation between the circumnuclear LHCN and Mdyn that can be explained by a gas fraction of 10% and a conversion factor {\alpha}HCN $\sim$ 10 between gas mass and HCN luminosity. Finally, adopting a different perspective to probe the physical properties of the gas around AGN, we report on an analysis of molecular line ratios which indicates that the clouds in this region are not self-gravitating.

We study a set of 3319 galaxies in the redshift interval 0.04 < z < 0.15 with far-infrared (FIR) coverage from the Herschel Stripe 82 survey (HerS), and emission-line measurements, redshifts, stellar masses and star-formation rates (SFRs) from the SDSS (DR7) MPA/JHU database. About 40% of the sample are detected in the Herschel/SPIRE 250 micron band. Total infrared (TIR) luminosities derived from HerS and ALLWISE photometry allow us to compare infrared and optical estimates of SFR with unprecedented statistics for diverse classes of galaxies. We find excellent agreement between TIR-derived and emission line-based SFRs for H II galaxies. Other classes, such as active galaxies and evolved galaxies, exhibit systematic discrepancies between optical and TIR SFRs. We demonstrate that these offsets are attributable primarily to survey biases and the large intrinsic uncertainties of the D4000- and colour-based optical calibrations used to estimate the SDSS SFRs of these galaxies. Using a classification scheme which expands upon popular emission-line methods, we demonstrate that emission-line galaxies with uncertain classifications include a population of massive, dusty, metal-rich star-forming systems that are frequently neglected in existing studies. We also study the capabilities of infrared selection of star-forming galaxies. FIR selection reveals a substantial population of galaxies dominated by cold dust which are missed by the long-wavelength WISE bands. Our results demonstrate that Herschel large-area surveys offer the means to construct large, relatively complete samples of local star-forming galaxies with accurate estimates of SFR that can be used to study the interplay between nuclear activity and star-formation.

The optical classification of a Seyfert galaxy and whether it is considered X-ray absorbed are often used interchangeably. But there are many borderline cases and also numerous examples where the optical and X-ray classifications appear to be in conflict. In this article we re-visit the relation between optical obscuration and X-ray absorption in AGNs. We make use of our "dust color" method (Burtscher et al. 2015) to derive the optical obscuration A_V and consistently estimated X-ray absorbing columns using 0.3--150 keV spectral energy distributions. We also take into account the variable nature of the neutral gas column N_H and derive the Seyfert sub-classes of all our objects in a consistent way. We show in a sample of 25 local, hard-X-ray detected Seyfert galaxies (log L_X / (erg/s) ~ 41.5 - 43.5) that there can actually be a good agreement between optical and X-ray classification. If Seyfert types 1.8 and 1.9 are considered unobscured, the threshold between X-ray unabsorbed and absorbed should be chosen at a column N_H = 10^22.3 / cm^2 to be consistent with the optical classification. We find that N_H is related to A_V and that the N_H/A_V ratio is approximately Galactic or higher in all sources, as indicated previously. But in several objects we also see that deviations from the Galactic ratio are only due to a variable X-ray column, showing that (1) deviations from the Galactic N_H/A_V can simply be explained by dust-free neutral gas within the broad line region in some sources, that (2) the dust properties in AGNs can be similar to Galactic dust and that (3) the dust color method is a robust way to estimate the optical extinction towards the sublimation radius in all but the most obscured AGNs.

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

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

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

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

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

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

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by AGN. Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel-PACS observation of an extremely high J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3$\sigma$ upper limit of $2 \times 10^{-17}\,\text{W}\,\text{m}^{-2}$. We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO Spectral Line Energy Distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities $\sim 10^{6}-10^{7}\,\text{cm}^{-3}$, located at $1.6-5\,\text{pc}$ from the AGN, with column densities of at least $10^{25}\,\text{cm}^{-2}$. At such high column densities, however, dust absorbs most of the CO(40-39) line emission at $\lambda = 65.69\, \mu$m. Therefore, even if NGC 1068 has a molecular torus which radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.

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

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

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

We describe a complete volume limited sample of nearby active galaxies selected by their 14-195keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy properties. The active sample appears to have no bias in terms of AGN type, the only difference being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, log L_{14-195keV} [erg/s] = 42.4-43.7, the optically obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.

We combine two approaches to isolate the AGN luminosity at near-infrared wavelengths and relate the near-IR pure AGN luminosity to other tracers of the AGN. Using integral-field spectroscopic data of an archival sample of 51 local AGNs, we estimate the fraction of non-stellar light by comparing the nuclear equivalent width of the stellar 2.3 micron CO absorption feature with the intrinsic value for each galaxy. We compare this fraction to that derived from a spectral decomposition of the integrated light in the central arc second and find them to be consistent with each other. Using our estimates of the near-IR AGN light, we find a strong correlation with presumably isotropic AGN tracers. We show that a significant offset exists between type 1 and type 2 sources in the sense that type 1 sources are 7 (10) times brighter in the near-IR at log L_MIR = 42.5 (log L_X = 42.5). These offsets only becomes clear when treating infrared type 1 sources as type 1 AGNs. All AGNs have very red near-to-mid-IR dust colors. This, as well as the range of observed near-IR temperatures, can be explained with a simple model with only two free parameters: the obscuration to the hot dust and the ratio between the warm and hot dust areas. We find obscurations of A_V (hot) = 5 - 15 mag for infrared type 1 sources and A_V (hot) = 15 - 35 mag for type 2 sources. The ratio of hot dust to warm dust areas of about 1000 is nicely consistent with the ratio of radii of the respective regions as found by infrared interferometry.

Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole growth. Despite their ubiquitous presence, the origin and the acceleration mechanism of such powerful and extended winds is not yet understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess the impact of the outflow on the host galaxy interstellar medium (ISM). We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of CO(3-2) emission in this source allows us to infer the molecular gas content and compare it to the ISM mass derived from the dust emission. We then analyze the results in the context of recent insights on scaling relations, which describe the gas content of the overall population of star-forming galaxies at a similar redshifts. The Star formation efficiency (~100) and gas mass (M_gas=2.1-9.5x10^{10} M_sun) inferred from the CO(3-2) line depend on the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than expected for its observed M$_\star$, sSFR and redshift, based on the most up-to-date calibrations (with gas fraction <20% and depletion time scale <340 Myr). Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback able to remove the gas from the host

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

We report the detection of far-infrared (FIR) CO rotational emission from nearby active galactic nuclei (AGN) and starburst galaxies, as well as several merging systems and Ultra-Luminous Infrared Galaxies (ULIRGs). Using Herschel-PACS, we have detected transitions in the J$_{upp}$ = 14 - 20 range ($\lambda \sim$ 130 - 185 $\mu$m, $\nu \sim$ 1612 - 2300 GHz) with upper limits on (and in two cases, detections of) CO line fluxes up to J$_{upp}$ = 30. The PACS CO data obtained here provide the first well-sampled FIR extragalactic CO SLEDs for this range, and will be an essential reference for future high redshift studies. We find a large range in the overall SLED shape, even amongst galaxies of similar type, demonstrating the uncertainties in relying solely on high-J CO diagnostics to characterize the excitation source of a galaxy. Combining our data with low-J line intensities taken from the literature, we present a CO ratio-ratio diagram and discuss its potential diagnostic value in distinguishing excitation sources and physical properties of the molecular gas. The position of a galaxy on such a diagram is less a signature of its excitation mechanism, than an indicator of the presence (or absence) of warm, dense molecular gas. We then quantitatively analyze the CO emission from a subset of the detected sources with Large Velocity Gradient (LVG) radiative transfer models to fit the CO SLEDs. Using both single-component and two-component LVG models to fit the kinetic temperature, velocity gradient, number density and column density of the gas, we derive the molecular gas mass and the corresponding CO-to-H$_2$ conversion factor, $\alpha_{CO}$, for each respective source. For the ULIRGs we find $\alpha$ values in the canonical range 0.4 - 5 M$_\odot$/(K kms$^{-1}$pc$^2$), while for the other objects, $\alpha$ varies between 0.2 and 14.} Finally, we compare our best-fit LVG model ..

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

Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution.

We present the KMOS^3D survey, a new integral field survey of over 600 galaxies at 0.7<z<2.7 using KMOS at the Very Large Telescope (VLT). The KMOS^3D survey utilizes synergies with multi-wavelength ground and space-based surveys to trace the evolution of spatially-resolved kinematics and star formation from a homogeneous sample over 5 Gyrs of cosmic history. Targets, drawn from a mass-selected parent sample from the 3D-HST survey, cover the star formation-stellar mass ($M_*$) and rest-frame $(U-V)-M_*$ planes uniformly. We describe the selection of targets, the observations, and the data reduction. In the first year of data we detect Halpha emission in 191 $M_*=3\times10^{9}-7\times10^{11}$ Msun galaxies at z=0.7-1.1 and z=1.9-2.7. In the current sample 83% of the resolved galaxies are rotation-dominated, determined from a continuous velocity gradient and $v_{rot}/\sigma>1$, implying that the star-forming 'main sequence' (MS) is primarily composed of rotating galaxies at both redshift regimes. When considering additional stricter criteria, the Halpha kinematic maps indicate at least ~70% of the resolved galaxies are disk-like systems. Our high-quality KMOS data confirm the elevated velocity dispersions reported in previous IFS studies at z>0.7. For rotation-dominated disks, the average intrinsic velocity dispersion decreases by a factor of two from 50 km/s at z~2.3 to 25 km/s at z~0.9 while the rotational velocities at the two redshifts are comparable. Combined with existing results spanning z~0-3, disk velocity dispersions follow an approximate (1+z) evolution that is consistent with the dependence of velocity dispersion on gas fractions predicted by marginally-stable disk theory.

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

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

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

We present multiline CO observations of the complex submillimeter galaxy SMM J00266+1708. Using the Zpectrometer on the Green Bank Telescope, we provide the first precise spectroscopic measurement of its redshift (z=2.742). Based on followup CO(1-0), CO(3-2), and CO(5-4) mapping, SMM J00266+1708 appears to have two distinct components separated by ~500 km/s that are nearly coincident along our line of sight. The two components show hints of different kinematics, with the blue-shifted component dispersion-dominated and the red-shifted component showing a clear velocity gradient. CO line ratios differ slightly between the two components, indicating that the physical conditions in their molecular gas may not be alike. We tentatively infer that SMM J00266+1708 is an ongoing merger with a mass ratio of (7.8+/-4.0)/sin^2(i), with its overall size and surface brightness closely resembling that of other merging systems. We perform large velocity gradient modeling of the CO emission from both components and find that each component's properties are consistent with a single phase of molecular gas (i.e., a single temperatures and density); additional multi-phase modeling of the red-shifted component, although motivated by a CO(1-0) size larger than the CO(3-2) size, is inconclusive. SMM J00266+1708 provides evidence of early stage mergers within the submillimeter galaxy population. Continuum observations of J00266 at the ~1" resolution of our observations could not have distinguished between the two components due to their separation (0.73" +/- 0.06"), illustrating that the additional velocity information provided by spectral line studies is important for addressing the prevalence of unresolved galaxy pairs in low-resolution submillimeter surveys.

We present X-shooter at the Vewry Large Telescope observations of a sample of 10 luminous, X-ray obscured QSOs at z$\sim1.5$ from the XMM-COSMOS survey, expected to be caught in the transitioning phase from starburst to active galactic nucleus (AGN)-dominated systems. The main selection criterion is X-ray detection at bright fluxes (L$>=10^{44}$ erg s$^{-1}$) coupled to red optical-to-NIR-to-MIR colors. Thanks to its large wavelength coverage, X-shooter allowed us to determine accurate redshifts from the presence of multiple emission lines for five out of six targets for which we had only a photometric redshift estimate, with an 80 percent success rate, significantly larger than what is observed in similar programs of spectroscopic follow-up of red QSOs. We report the detection of broad and shifted components in the [OIII]$\lambda\lambda$5007,4959 complexes for six out of eight sources with these lines observable in regions free from strong atmospheric absorptions. The full width half-maximum (FWHM) associated with the broad components are in the range FWHM$\sim900-1600$ km s$^{-1}$, larger than the average value observed in SDSS Type 2 AGN samples at similar observed [OIII] luminosity, but comparable to those observed for QSO/ultraluminous infrared galaxies systems for which the presence of kpc scale outflows have been revealed through integral field unit spectroscopy. Although the total outflow energetics (inferred under reasonable assumptions) may be consistent with winds accelerated by stellar processes, we favour an AGN origin for the outflows given the high outflow velocities oberved (v$>1000$ km s$^{-1}$) and the presence of strong winds also in objects undetected in the far infrared.

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

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

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

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

We present near-infrared integral field spectroscopy of ten Herschel selected galaxies at z~1.5. From detailed mapping of the H$\alpha$ and [NII] emission lines we trace the dynamics, star formation rates, metallicities and also investigate gas fractions for these galaxies. For a few galaxies the distribution of star formation as traced by H$\alpha$ only traces a small fraction of the stellar disc, which could be tracing recent minor merging events. The rest-frame ultraviolet (UV) continuum emission often has a distribution completely different from H$\alpha$, which warns about the use of UV-SED based star formation tracers in these systems. Our analysis of galaxy dynamics shows that minor dynamical disruptions (e.g. minor merging) are generally not enough to cause a deviation from the established 'Main Sequence' relation. Most galaxies are found to follow the fundamental metallicity relation (FMR), although with large scatter. One galaxy, (a small satellite galaxy of a massive companion) is found to deviate strongly from the FMR. This deviation is in nice agreement with the correlation recently discovered in local galaxies between gas metallicity and environment, which has been ascribed to enriched inter-galactic medium (IGM) in dense environments, and therefore suggests that here the IGM was already significantly enriched by z~1.5.

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

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

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

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

We present the activities carried out to calibrate and characterise the performance of the elements of attitude control and measurement on board the Herschel spacecraft. The main calibration parameters and the evolution of the indicators of the pointing performance are described, from the initial values derived from the observations carried out in the performance verification phase to those attained in the last year and half of mission, an absolute pointing error around or even below 1 arcsec, a spatial relative pointing error of some 1 arcsec and a pointing stability below 0.2 arsec. The actions carried out at the ground segment to improve the spacecraft pointing measurements are outlined. On-going and future developments towards a final refinement of the Herschel astrometry are also summarised. A brief description of the different components of the attitude control and measurement system (both in the space and in the ground segments) is also given for reference. We stress the importance of the cooperation between the different actors (scientists, flight dynamics and systems engineers, attitude control and measurement hardware designers, star-tracker manufacturers, etc.) to attain the final level of performance.

Using a sample of 25683 star-forming and 2821 passive galaxies at $z\sim2$, selected in the COSMOS field following the BzK color criterion, we study the hosting halo mass and environment of galaxies as a function of their physical properties. Spitzer and Herschel provide accurate SFR estimates for starburst galaxies. We measure the auto- and cross-correlation functions of various galaxy sub-samples and infer the properties of their hosting halos using both an HOD model and the linear bias at large scale. We find that passive and star-forming galaxies obey a similarly rising relation between the halo and stellar mass. The mean host halo mass of star forming galaxies increases with the star formation rate between 30 and 200 M$_\odot$.yr$^{-1}$, but flattens for higher values, except if we select only main-sequence galaxies. This reflects the expected transition from a regime of secular co-evolution of the halos and the galaxies to a regime of episodic starburst. We find similar large scale biases for main-sequence, passive, and starburst galaxies at equal stellar mass, suggesting that these populations live in halos of the same mass. We detect an excess of clustering on small scales for passive galaxies and showed, by measuring the large-scale bias of close pairs, that this excess is caused by a small fraction ($\sim16%$) of passive galaxies being hosted by massive halos ($\sim 3 \times 10^{13}$ M$_\odot$) as satellites. Finally, extrapolating the growth of halos hosting the z$\sim$2 population, we show that M$_\star \sim 10^{10}$ M$_\odot$ galaxies at z$\sim$2 will evolve, on average, into massive (M$_\star \sim 10^{11}$ M$_\odot$), field galaxies in the local Universe and M$_\star \sim 10^{11}$ M$_\odot$ galaxies at z=2 into local, massive, group galaxies. The most massive main-sequence galaxies and close pairs of massive, passive galaxies end up in today's clusters.

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

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth.

We have investigated the gas content of a sample of several hundred AGN host galaxies at z$<$1 and compared it with a sample of inactive galaxies, matched in bins of stellar mass and redshift. Gas masses have been inferred from the dust masses, obtained by stacked Herschel far-IR and sub-mm data in the GOODS and COSMOS fields, under reasonable assumptions and metallicity scaling relations for the dust-to-gas ratio. We find that AGNs are on average hosted in galaxies much more gas rich than inactive galaxies. In the vast majority of stellar mass bins, the average gas content of AGN hosts is higher than in inactive galaxies. The difference is up to a factor of ten higher in low stellar mass galaxies, with a significance of 6.5$\sigma$. In almost half of the AGN sample the gas content is three times higher than in the control sample of inactive galaxies. Our result strongly suggests that the probability of having an AGN activated is simply driven by the amount of gas in the host galaxy; this can be explained in simple terms of statistical probability of having a gas cloud falling into the gravitational potential of the black hole. The increased probability of an AGN being hosted by a star-forming galaxy, identified by previous works, may be a consequence of the relationship between gas content and AGN activity, found in this paper, combined with the Schmidt-Kennicutt law for star formation.

We report the detection of ubiquitous powerful nuclear outflows in massive (> 10^11 Msun) z~2 star-forming galaxies (SFGs), which are plausibly driven by an Active Galactic Nucleus (AGN). The sample consists of the eight most massive SFGs from our SINS/zC-SINF survey of galaxy kinematics with the imaging spectrometer SINFONI, six of which have sensitive high-resolution adaptive optics (AO) assisted observations. All of the objects are disks hosting a significant stellar bulge. The spectra in their central regions exhibit a broad component in Halpha and forbidden [NII] and [SII] line emission, with typical velocity FWHM ~ 1500 km/s, [NII]/Halpha ratio ~ 0.6, and intrinsic extent of 2 - 3 kpc. These properties are consistent with warm ionized gas outflows associated with Type 2 AGN, the presence of which is confirmed via independent diagnostics in half the galaxies. The data imply a median ionized gas mass outflow rate of ~ 60 Msun/yr and mass loading of ~ 3. At larger radii, a weaker broad component is detected but with lower FWHM ~ 485 km/s and [NII]/Halpha ~ 0.35, characteristic for star formation-driven outflows as found in the lower-mass SINS/zC-SINF galaxies. The high inferred mass outflow rates and frequent occurrence suggest the nuclear outflows efficiently expel gas out of the centers of the galaxies with high duty cycles, and may thus contribute to the process of star formation quenching in massive galaxies. Larger samples at high masses will be crucial to confirm the importance and energetics of the nuclear outflow phenomenon, and its connection to AGN activity and bulge growth.

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

Roughly half of the radiation from evolving galaxies in the early universe reaches us in the far-infrared and submillimeter wavelength range. Recent major advances in observing capabilities, in particular the launch of the Herschel Space Observatory in 2009, have dramatically enhanced our ability to use this information in the context of multiwavelength studies of galaxy evolution. Near its peak, three quarters of the cosmic infrared background is now resolved into individually detected sources. The use of far-infrared diagnostics of dust-obscured star formation and of interstellar medium conditions has expanded from rare extreme high-redshift galaxies to more typical main sequence galaxies and hosts of active galactic nuclei, out to z>~2. These studies shed light on the evolving role of steady equilibrium processes and of brief starbursts, at and since the peak of cosmic star formation and black hole accretion. This review presents a selection of recent far-infrared studies of galaxy evolution, with an emphasis on Herschel results

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

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

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

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

We present a study of the morphology and intensity of star formation in the host galaxies of eight Palomar-Green quasars using observations with the Hubble Space Telescope. Our observations are motivated by recent evidence for a close relationship between black hole growth and the stellar mass evolution in its host galaxy. We use narrow-band [O II] $\lambda$3727, H$\beta$, [O III] $\lambda$5007 and Pa$\alpha$ images, taken with the WFPC2 and NICMOS instruments, to map the morphology of line-emitting regions, and, after extinction corrections, diagnose the excitation mechanism and infer star-formation rates. Significant challenges in this type of work are the separation of the quasar light from the stellar continuum and the quasar-excited gas from the star-forming regions. To this end, we present a novel technique for image decomposition and subtraction of quasar light. Our primary result is the detection of extended line-emitting regions with sizes ranging from 0.5 to 5 kpc and distributed symmetrically around the nucleus, powered primarily by star formation. We determine star-formation rates of order a few tens of M$_\odot$/yr. The host galaxies of our target quasars have stellar masses of order $10^{11}$ M$_\odot$ and specific star formation rates on a par with those of M82 and luminous infrared galaxies. As such they fall at the upper envelope or just above the star-formation mass sequence in the specific star formation vs stellar mass diagram. We see a clear trend of increasing star formation rate with quasar luminosity, reinforcing the link between the growth of the stellar mass of the host and the black hole mass found by other authors.

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

We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excitation, presumably spatially extended. Particularly prominent, the blue wing of the absorption detected in the in-ladder 2Pi_{3/2} J=9/2-7/2 OH doublet at 65 um, with E_lower=290 K, indicates that the excited outflowing gas is generated in a compact and warm (circum)nuclear region. Because the excited, outflowing OH gas in Mrk 231 is associated with the warm, far-IR continuum source, it is likely more compact (diameter of 200-300 pc) than that probed by CO and HCN. Nevertheless, its mass-outflow rate per unit of solid angle as inferred from OH is similar to that previously derived from CO, >~70x(2.5x10^{-6}/X_{OH}) Msun yr^{-1} sr^{-1}, where X_{OH} is the OH abundance relative to H nuclei. In spherical symmetry, this would correspond to >~850x(2.5x10^{-6}/X_{OH}) Msun yr^{-1}, though significant collimation is inferred from the line profiles. The momentum flux of the excited component attains ~15 L_{AGN}/c, with an OH column density of (1.5-3)x10^{17} cm^-2 and a mechanical luminosity of ~10^{11} Lsun. The detection of very excited OH peaking at central velocities indicates the presence of a nuclear reservoir of gas rich in OH, plausibly the 130-pc scale circumnuclear torus previously detected in OH megamaser emission, that may be feeding the outflow. An exceptional ^{18}OH enhancement, with OH/^{18}OH<~30 at both central and blueshifted velocities, is likely the result of interstellar-medium processing by recent starburst/SNe activity.

We analyze the resolved stellar populations of 473 massive star-forming galaxies at 0.7 < z < 1.5, with multi-wavelength broad-band imaging from CANDELS and Halpha surface brightness profiles at the same kiloparsec resolution from 3D-HST. Together, this unique data set sheds light on how the assembled stellar mass is distributed within galaxies, and where new stars are being formed. We find the Halpha morphologies to resemble more closely those observed in the ACS I band than in the WFC3 H band, especially for the larger systems. We next derive a novel prescription for Halpha dust corrections, which accounts for extra extinction towards HII regions. The prescription leads to consistent SFR estimates and reproduces the observed relation between the Halpha/UV luminosity ratio and visual extinction, both on a pixel-by-pixel and on a galaxy-integrated level. We find the surface density of star formation to correlate with the surface density of assembled stellar mass for spatially resolved regions within galaxies, akin to the so-called 'main sequence of star formation' established on a galaxy-integrated level. Deviations from this relation towards lower equivalent widths are found in the inner regions of galaxies. Clumps and spiral features, on the other hand, are associated with enhanced Halpha equivalent widths, bluer colors, and higher specific star formation rates compared to the underlying disk. Their Halpha/UV luminosity ratio is lower than that of the underlying disk, suggesting the ACS clump selection preferentially picks up those regions of elevated star formation activity that are the least obscured by dust. Our analysis emphasizes that monochromatic studies of galaxy structure can be severely limited by mass-to-light ratio variations due to dust and spatially inhomogeneous star formation histories.

As part of the SINS/zC-SINF surveys of high-z galaxy kinematics, we derive the radial distributions of H-alpha surface brightness, stellar mass surface density, and dynamical mass at ~2 kpc resolution in 19 z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. From these data we infer the radial distribution of the Toomre Q-parameter for these main-sequence star forming galaxies (SFGs), covering almost two decades of stellar mass (10^9.6 to 10^11.5 solar masses). In more than half of our SFGs, the H-alpha distributions cannot be fit by a centrally peaked distribution, such as an exponential, but are better described by a ring, or the combination of a ring and an exponential. At the same time the kinematic data indicate the presence of a mass distribution more centrally concentrated than a single exponential distribution for 5 of the 19 galaxies. The resulting Q-distributions are centrally peaked for all, and significantly exceed unity there for three quarters of the SFGs. The occurrence of H-alpha rings and of large nuclear Q-values is strongly correlated, and is more common for the more massive SFGs. While our sample is small and there remain substantial uncertainties and caveats, our observations are consistent with a scenario in which cloud fragmentation and global star formation are secularly suppressed in gas rich high-z disks from the inside out, as the central stellar mass density of the disks grows.

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

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

In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.

This paper provides an overview of the PACS photometer flux calibration concept, in particular for the principal observation mode, the scan map. The absolute flux calibration is tied to the photospheric models of five fiducial stellar standards (alpha Boo, alpha Cet, alpha Tau, beta And, gamma Dra). The data processing steps to arrive at a consistent and homogeneous calibration are outlined. In the current state the relative photometric accuracy is around 2% in all bands. Starting from the present calibration status, the characterization and correction for instrumental effects affecting the relative calibration accuracy is described and an outlook for the final achievable calibration numbers is given. After including all the correction for the instrumental effects, the relative photometric calibration accuracy (repeatability) will be as good as 0.5% in the blue and green band and 2% in the red band. This excellent calibration starts to reveal possible inconsistencies between the models of the K-type and the M-type stellar calibrators. The absolute calibration accuracy is therefore mainly limited by the 5% uncertainty of the celestial standard models in all three bands. The PACS bolometer response was extremely stable over the entire Herschel mission and a single, time-independent response calibration file is sufficient for the processing and calibration of the science observations. The dedicated measurements of the internal calibration sources were needed only to characterize secondary effects. No aging effects of the bolometer or the filters have been found. Also, we found no signs of filter leaks. The PACS photometric system is very well characterized with a constant energy spectrum nu*Fnu = lambda*Flambda = const as a reference. Colour corrections for a wide range of sources SEDs are determined and tabulated.

In the era of deep, large-area far-infrared (FIR) surveys from the Herschel Space Telescope, the bulk of the star-formation in distant galaxies, once hidden by dust, is now being revealed. The FIR provides probably the cleanest view of SF in the host galaxies of Active Galactic Nuclei (AGNs) over cosmic time. We report results from studies of the relationships between SF, AGN activity and AGN obscuration out to z=2.5, which employ some of the deepest Herschel and X-ray datasets currently available, while spanning orders of magnitude in the dynamic range of AGN properties. We highlight the role of gaseous supply in modulating both SF and AGN activity without necessarily implying a direct causal connection between these phenomenon. The role of starburst- or major merger-fueled AGN activity at low and high redshifts is discussed in the context of our results.

We explore the nature of the long-wavelength mid-infrared (MIR) emission of a sample of 13000 local Type II (narrow-line) Active Galactic Nuclei (AGNs) from the Sloan Digital Sky Survey (SDSS) using 12 and 22 micron photometry from the WISE all-sky survey. In combination with FIRST 1.4 GHz measurements, we show that AGNs divide into two relatively distinct populations or "branches" in the plane of MIR and radio luminosity. Seyfert galaxies lie almost exclusively on a MIR-bright branch (Branch A), while low-ionization nuclear emission line galaxies (LINERs) are split evenly into Branch A and the MIR-faint Branch B. We devise various tests to constrain the processes that define the branches, including a comparison to the properties of pure star-forming (SF) inactive galaxies on the MIR-Radio plane. We demonstrate that the total MIR emission of objects on Branch A, including most Seyfert galaxies, is governed primarily by host star-formation, with about 15% of the 22 micron luminosity coming from AGN-heated dust. This implies that on-going dusty star-formation is a general property of Seyfert host galaxies. We show that the 12 micron broad-band luminosity of AGNs on Branch A is suppressed with respect to star-forming galaxies, possibly due to the destruction of PAHs or deeper 10 microns Si absorption in AGNs. We uncover a correlation between the MIR luminosity and [O III] luminosity in AGNs. This suggests a relationship between the SFR and nuclear luminosity in the AGN population, but we caution on the importance of selection effects inherent to such AGN-dominated emission-line galaxies in driving such a correlation. We highlight the MIR-radio plane as a useful tool in comparative studies of SF and nuclear activity in AGN.

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

We report the results from a systematic search for molecular (OH-119 um) outflows with Herschel-PACS in a sample of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. We find that the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7-um silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the detection of OH absorption profiles with median velocities more blueshifted than -50 km/sec, is seen in 26 (70%) of the 37 OH-detected targets, suggesting a wide-angle (~145 degrees) outflow geometry. Conversely, unambiguous evidence for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted than +50 km/sec, is seen in only 4 objects, suggesting a planar or filamentary geometry for the inflowing gas. Terminal outflow velocities of ~-1000 km/sec are measured in several objects, but median outflow velocities are typically ~-200 km s^{-1}. While the outflow velocities show no statistically significant dependence on the star formation rate, they are distinctly more blueshifted among systems with large AGN fractions and luminosities [log (L_AGN / L_sun) > 11.8 +/- 0.3]. The quasars in these systems play a dominant role in driving the molecular outflows. In contrast, the most AGN dominated systems, where OH is seen purely in emission, show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular outflows subside once the quasar has cleared a path through the obscuring material.

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

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

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

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

Based on high-resolution, spatially resolved data of 10 z~2 star-forming galaxies from the SINS/zC-SINF survey and LUCI data for 12 additional galaxies, we probe the excitation properties of high-z galaxies and the impact of active galactic nuclei (AGN), shocks and photoionization. We explore how these spatially-resolved line ratios can inform our interpretation of integrated emission line ratios obtained at high redshift. Many of our galaxies fall in the `composite' region of the z~0 [NII]/Ha versus [OIII]/Hb diagnostic (BPT) diagram, between star-forming galaxies and those with AGN. Based on our resolved measurements, we find that some of these galaxies likely host an AGN, while others appear to be affected by the presence of shocks possibly caused by an outflow or from enhanced ionization parameter as compared with HII regions in normal local star-forming galaxies. We find that the Mass-Excitation (MEx) diagnostic, which separates purely star-forming and AGN hosting local galaxies in the [OIII]/Hb versus stellar mass plane, does not properly separate z~2 galaxies classified according to the BPT diagram. However, if we shift the galaxies based on the offset between the local and z~2 mass-metallicity relation (i.e. to the mass they would have at z~0 with the same metallicity), we find better agreement between the MEx and BPT diagnostics. Finally, we find that metallicity calibrations based on [NII]/Ha are more biased by shocks and AGN at high-z than the [OIII]/Hb/[NII]/Ha calibration.

We present F110W (~J) and F160W (~H) observations of ten submillimeter galaxies (SMGs) obtained with the Hubble Space Telescope's (HST's) NICMOS camera. Our targets have optical redshifts in the range 2.20<z<2.81 confirmed by millimeter CO or mid-IR spectroscopy, guaranteeing that the two bands sample the rest-frame optical with the Balmer break falling between them. Eight of ten are detected in both bands, while two are detected in F160W only. We study their F160W morphologies, applying a maximum-deblending detection algorithm to distinguish multiple- from single-component configurations, leading to reassessments for several objects. Based on our NICMOS imaging and/or previous dynamical evidence we identify five SMGs as multiple sources, which we interpret as merging systems. Additionally, we calculate morphological parameters asymmetry (A) and Gini coefficient (G); thanks to our sample's limited redshift range we recover the trend that multiple-component, merger-like morphologies are reflected in higher asymmetries. We analyze the stellar populations of nine objects with F110W/F160W photometry, using archival HST optical data when available. For multiple systems, we are able to model the individual components that build up an SMG. With the available data we cannot discriminate among star formation histories, but we constrain stellar masses and mass ratios for merger-like SMG systems, obtaining a mean log(M_*/M_sun)=10.9+/-0.2 for our full sample, with individual values log(M_*/M_sun)~9.6-11.8. The morphologies and mass ratios of the least and most massive systems match the predictions of the major-merger and cold accretion SMG formation scenarios, respectively, suggesting that both channels may have a role in the population's origin.

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

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

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

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

Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us to study the star formation efficiency on subgalactic scales around redshift z = 1.2. We present a method for obtaining the gas and star formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding scales are not resolved. This method is based on identifying these structures in position-velocity diagrams corresponding to slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star forming distant galaxies - EGS13003805, EGS13004291, EGS12007881, and EGS13019128 in the AEGIS terminology - to determine the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [OII] line. We identify the ensembles of clumps with the angular resolution available in both CO and [OII] spectroscopy; i.e., 1-1.5". SFR and gas surface densities are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation on a scale of ~8 kpc. The data generally indicates an average depletion time of 1.9 Gyr, but with significant variations from point to point within the galaxies.

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

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

We report matched resolution, imaging spectroscopy of the CO J=3-2 line (with the IRAM Plateau de Bure millimeter interferometer) and of the H-alpha line (with LUCI at the Large Binocular Telescope)in the massive z=1.53 main-sequence galaxy EGS 13011166, as part of the "Plateau de Bure high-z, blue sequence survey (PHIBSS). We combine these data with HST V-J-J-H-band maps to derive spatially resolved distributions of stellar surface density, star formation rate, molecular gas surface density, optical extinction and gas kinematics. The spatial distribution and kinematics of the ionized and molecular gas are remarkably similar and are well modeled by a turbulent, globally Toomre unstable rotating disk. The stellar surface density distribution is smoother than the clumpy rest-frame UV/optical light distribution, and peaks in an obscured, star forming massive bulge near the dynamical center. The molecular gas surface density and the effective optical screen extinction track each other and are well modeled by a 'mixed' extinction model. The inferred slope of the spatially resolved molecular gas to star formation rate relation depends strongly on the adopted extinction model and can vary from 0.8 to 1.7. For the preferred mixed dust-gas model we find a near linear slope.

We analyze the spectra, spatial distributions and kinematics of Ha, [NII] and [SII] emission in a sample of 42, z~2.2 UV/optically selected star forming galaxies (SFGs) from the SINS & zC-SINF surveys, 35 of which were observed in the adaptive optics mode of SINFONI. This is supplemented by kinematic data from 48 z~1-2.5 galaxies from the literature. We find that the kinematic classification of the high-z SFGs as `dispersion dominated' or `rotation dominated' correlates most strongly with their intrinsic sizes. Smaller galaxies are more likely `dispersion-dominated' for two main reasons: 1) The rotation velocity scales linearly with galaxy size but intrinsic velocity dispersion does not depend on size, and as such, their ratio is systematically lower for smaller galaxies, and 2) Beam smearing strongly decreases large-scale velocity gradients and increases observed dispersion much more for galaxies with sizes at or below the resolution. Dispersion dominated SFGs may thus have intrinsic properties similar to `rotation dominated' SFGs, but are primarily more compact, lower mass, less metal enriched and may have higher gas fractions, plausibly because they represent an earlier evolutionary state.

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

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

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

Active Galactic Nuclei (AGN) play a decisive role in galaxy evolution, particularly so when operating in a radiatively inefficient mode, where they launch powerful jets that reshape their surroundings. However, identifying them is difficult, since radio observations commonly have resolutions of between 1 arcsec and 10 arcsec, which is equally sensitive to radio emission from star-forming activity and from AGN. Very Long Baseline Interferometry (VLBI) observations allow one to filter out all but the most compact non-thermal emission from radio survey data. The observational and computational demands to do this in large surveys have been, until recently, too high to make such undertakings feasible. Only the recent advent of wide-field observing techniques have facilitated such observations, and we here present the results from a survey of 217 radio sources in the Lockman Hole/XMM field. We describe in detail some new aspects of the calibration, including primary beam correction, multi-source self-calibration, and mosaicing. As a result, we detected 65 out of the 217 radio sources and were able to construct, for the first time, the source counts of VLBI-detected AGN. They indicate that at least 15%-25% of the sub-mJy radio sources are AGN-driven, consistent with recent findings using other AGN selection techniques. We have used ancillary data to investigate the AGN hosts. We find that among the sources nearby enough to be resolved in the optical images, 88% (23/26) could be classified as early-type or bulge-dominated galaxies. While 50% of these sources are correctly represented by the SED of an early-type galaxy, for the rest the best fit was obtained with a heavily extinct starburst template, an effect we ascribe to a degeneracy in the fit. Overall, the typical hosts of VLBI-detected sources are in good agreement with being early-type or bulge-dominated galaxies.

Mrk 231, the nearest (z = 0.0422) quasar, hosts both a galactic-scale wind and a nuclear-scale iron low-ionization broad absorption line (FeLoBAL) outflow. We recently obtained a far-ultraviolet (FUV) spectrum of this object covering ~1150 - 1470 A with the Cosmic Origins Spectrograph on board the Hubble Space Telescope. This spectrum is highly peculiar, highlighted by the presence of faint (~< 2% of predictions based on H-alpha), broad (>~ 10,000 km/s at the base), and highly blueshifted (centroid at ~ -3500 km/s) Ly-alpha emission. The FUV continuum emission is slightly declining at shorter wavelengths (consistent with F_lambda ~ lambda^1.7) and does not show the presence of any obvious photospheric or wind stellar features. Surprisingly, the FUV spectrum also does not show any unambiguous broad absorption features. It thus appears to be dominated by the AGN, rather than hot stars, and virtually unfiltered by the dusty FeLoBAL screen. The observed Ly-alpha emission is best explained if it is produced in the outflowing BAL cloud system, while the Balmer lines arise primarily from the standard broad emission line region seen through the dusty (A_V ~ 7 mag.) broad absorption line region. Two possible geometric models are discussed in the context of these new results.

Based on SINFONI Ha, [NII] and [SII] AO data of 30 z \sim 2 star-forming galaxies (SFGs) from the SINS and zcSINF surveys, we find a strong correlation of the Ha broad flux fraction with the star formation surface density of the galaxy, with an apparent threshold for strong outflows occurring at 1 Msun yr^-1 kpc^-2. Above this threshold, we find that SFGs with logm_\ast>10 have similar or perhaps greater wind mass loading factors (eta = Mdotout/SFR) and faster outflow velocities than lower mass SFGs. This trend suggests that the majority of outflowing gas at z \sim 2 may derive from high-mass SFGs, and that the z \sim 2 mass-metallicity relation is driven more by dilution of enriched gas in the galaxy gas reservoir than by the efficiency of outflows. The mass loading factor is also correlated with the SFR and inclination, such that more star-forming and face-on galaxies launch more powerful outflows. For galaxies that have evidence for strong outflows, we find that the broad emission is spatially extended to at least the half-light radius (\sim a few kpc). We propose that the observed threshold for strong outflows and the observed mass loading of these winds can be explained by a simple model wherein break-out of winds is governed by pressure balance in the disk. Using the ratio of the [SII] doublet in a broad and narrow component, we find that outflowing gas has a density of \sim10-100 cm^-3, significantly less than that of the star forming gas (600 cm^-3).

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 PHIBSS, the IRAM Plateau de Bure high-z blue sequence CO 3-2 survey of the molecular gas properties in normal star forming galaxies (SFGs) near the cosmic star formation peak. PHIBSS provides 52 CO detections in two redshift slices at z~1.2 and 2.2, with log(M*(M_solar))>10.4 and log(SFR(M_solar/yr))>1.5. Including a correction for the incomplete coverage of the M*-SFR plane, we infer average gas fractions of ~0.33 at z~1.2 and ~0.47 at z~2.2. Gas fractions drop with stellar mass, in agreement with cosmological simulations including strong star formation feedback. Most of the z~1-3 SFGs are rotationally supported turbulent disks. The sizes of CO and UV/optical emission are comparable. The molecular gas - star formation relation for the z=1-3 SFGs is near-linear, with a ~0.7 Gyrs gas depletion timescale; changes in depletion time are only a secondary effect. Since this timescale is much less than the Hubble time in all SFGs between z~0 and 2, fresh gas must be supplied with a fairly high duty cycle over several billion years. At given z and M*, gas fractions correlate strongly with the specific star formation rate. The variation of specific star formation rate between z~0 and 3 is mainly controlled by the fraction of baryonic mass that resides in cold gas.

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

We report on Herschel/PACS observations of absorption lines of OH+, H2O+ and H3O+ in NGC 4418 and Arp 220. Excited lines of OH+ and H2O+ with E_lower of at least 285 and \sim200 K, respectively, are detected in both sources, indicating radiative pumping and location in the high radiation density environment of the nuclear regions. Abundance ratios OH+/H2O+ of 1-2.5 are estimated in the nuclei of both sources. The inferred OH+ column and abundance relative to H nuclei are (0.5-1)x10^{16} cm-2 and \sim2x10^{-8}, respectively. Additionally, in Arp 220, an extended low excitation component around the nuclear region is found to have OH+/H2O+\sim5-10. H3O+ is detected in both sources with N(H3O+)\sim(0.5-2)x10^{16} cm-2, and in Arp 220 the pure inversion, metastable lines indicate a high rotational temperature of ~500 K, indicative of formation pumping and/or hot gas. Simple chemical models favor an ionization sequence dominated by H+ - O+ - OH+ - H2O+ - H3O+, and we also argue that the H+ production is most likely dominated by X-ray/cosmic ray ionization. The full set of observations and models leads us to propose that the molecular ions arise in a relatively low density (\gtrsim10^4 cm-3) interclump medium, in which case the ionization rate per H nucleus (including secondary ionizations) is zeta>10^{-13} s-1, a lower limit that is severalx10^2 times the highest rate estimates for Galactic regions. In Arp 220, our lower limit for zeta is compatible with estimates for the cosmic ray energy density inferred previously from the supernova rate and synchrotron radio emission, and also with the expected ionization rate produced by X-rays. In NGC 4418, we argue that X-ray ionization due to an AGN is responsible for the molecular ion production.

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

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

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]

[Abridged] We present maps of the main cooling lines of the neutral atomic gas ([OI] at 63 and 145 micron and [CII] at 158 micron) and in the [OIII] 88 micron line of the starburst galaxy M82, carried out with the PACS spectrometer on board the Herschel satellite. By applying PDR modeling we derive maps of the main ISM physical parameters, including the [CII] optical depth, at unprecedented spatial resolution (~300 pc). We can clearly kinematically separate the disk from the outflow in all lines. The [CII] and [OI] distributions are consistent with PDR emission both in the disk and in the outflow. Surprisingly, in the outflow, the atomic and the ionized gas traced by the [OIII] line both have a deprojected velocity of ~75 km/s, very similar to the average velocity of the outflowing cold molecular gas (~ 100 km/s) and several times smaller than the outflowing material detected in Halpha (~ 600 km/s). This suggests that the cold molecular and neutral atomic gas and the ionized gas traced by the [OIII] 88 micron line are dynamically coupled to each other but decoupled from the Halpha emitting gas. We propose a scenario where cold clouds from the disk are entrained into the outflow by the winds where they likely evaporate, surviving as small, fairly dense cloudlets (n_H\sim 500-1000 cm^-3, G_0\sim 500- 1000, T_gas\sim300 K). We show that the UV photons provided by the starburst are sufficient to excite the PDR shells around the molecular cores. The mass of the neutral atomic gas in the outflow is \gtrsim 5-12x 10^7 M_sun to be compared with that of the molecular gas (3.3 x 10^8 M_sun) and of the Halpha emitting gas (5.8 x 10^6 M_sun). The mass loading factor, (dM/dt)/SFR, of the molecular plus neutral atomic gas in the outflow is ~ 2. Energy and momentum driven outflow models can explain the data equally well, if all the outflowing gas components are taken into account.

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

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

We study the near-infrared properties of the super star cluster NGC1750-1 in order to constrain its spatial extent, its stellar population and its age. We use adaptive optics assisted integral field spectroscopy with SINFONI on the VLT. We estimate the spatial extent of the cluster and extract its K-band spectrum from which we constrain the age of the dominant stellar population. Our observations have an angular resolution of about 0.11", providing an upper limit on the cluster radius of 2.85+/-0.50 pc depending on the assumed distance. The K-band spectrum is dominated by strong CO absorption bandheads typical of red supergiants. Its spectral type is equivalent to a K4-5I star. Using evolutionary tracks from the Geneva and Utrecht groups, we determine an age of 12+/-6 Myr. The large uncertainty is rooted in the large difference between the Geneva and Utrecht tracks in the red supergiants regime. The absence of ionized gas lines in the K-band spectrum is consistent with the absence of O and/or Wolf-Rayet stars in the cluster, as expected for the estimated age.

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

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

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

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

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

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

Most theoretical models invoke quasar driven outflows to quench star formation in massive galaxies, this feedback mechanism is required to account for the population of old and passive galaxies observed in the local universe. The discovery of massive, old and passive galaxies at z=2, implies that such quasar feedback onto the host galaxy must have been at work very early on, close to the reionization epoch. We have observed the [CII]158um transition in SDSSJ114816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most distant massive outflow ever detected and is likely tracing the long sought quasar feedback, already at work in the early Universe. The outflow is marginally resolved on scales of about 16 kpc, implying that the outflow can really affect the whole galaxy, as required by quasar feedback models. The inferred outflow rate, dM/dt > 3500 Msun/yr, is the highest ever found. At this rate the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few million years.

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

We report the detection of far-IR CO rotational emission from the prototypical Seyfert 2 galaxy NGC 1068. Using Herschel-PACS, we have detected 11 transitions in the J_upper=14-30 (E_upper/k_B = 580-2565 K) range, all of which are consistent with arising from within the central 10" (700 pc). The detected transitions are modeled as arising from 2 different components: a moderate excitation (ME) component close to the galaxy systemic velocity, and a high excitation (HE) component that is blueshifted by ~80 km s^{-1}. We employ a large velocity gradient (LVG) model and derive n_H2~10^{5.6} cm^{-3}, T_kin~170 K, and M_H2~10^{6.7} M_sun for the ME component, and n_H2~10^{6.4} cm^{-3}, T_kin~570 K, and M_H2~10^{5.6} M_sun for the HE component, although for both components the uncertainties in the density and mass are plus/minus (0.6-0.9) dex. We compare the CO line profiles with those of other molecular tracers observed at higher spatial and spectral resolution, and find that the ME transitions are consistent with these lines arising in the ~200 pc diameter ring of material traced by H_2 1-0 S(1) observations. The blueshift of the HE lines may also be consistent with the bluest regions of this H_2 ring, but a better kinematic match is found with a clump of infalling gas ~40 pc north of the AGN. We discuss the prospects of placing the HE component near the AGN, and conclude that while the moderate thermal pressure precludes an association with the ~1 pc radius H_2O maser disk, the HE component could potentially be located only a few parsecs more distant from the AGN, and might then provide the N_H~10^{25} cm^{-2} column obscuring the nuclear hard X-rays. Finally, we also report sensitive upper limits extending up to J_upper=50, which place constraints on a previous model prediction for the CO emission from the X-ray obscuring torus. [Abridged]

We perform a detailed analysis of the resolved colors and stellar populations of a complete sample of 323 star-forming galaxies at 0.5 < z < 1.5, and 326 star-forming galaxies at 1.5 < z < 2.5 in the ERS and CANDELS-Deep region of GOODS-South. Galaxies were selected to be more massive than 10^10 Msun and have specific star formation rates above 1/t_H. We model the 7-band optical ACS + near-IR WFC3 spectral energy distributions of individual bins of pixels, accounting simultaneously for the galaxy-integrated photometric constraints available over a longer wavelength range. We analyze variations in rest-frame color, stellar surface mass density, age, and extinction as a function of galactocentric radius and local surface brightness/density, and measure structural parameters on luminosity and stellar mass maps. We find evidence for redder colors, older stellar ages, and increased dust extinction in the nuclei of galaxies. Big star-forming clumps seen in star formation tracers are less prominent or even invisible on the inferred stellar mass distributions. Off-center clumps contribute up to ~20% to the integrated SFR, but only 7% or less to the integrated mass of all massive star-forming galaxies at z ~ 1 and z ~ 2, with the fractional contributions being a decreasing function of wavelength used to select the clumps. The stellar mass profiles tend to have smaller sizes and M20 coefficients, and higher concentration and Gini coefficients than the light distribution. Our results are consistent with an inside-out disk growth scenario with brief (100 - 200 Myr) episodic local enhancements in star formation superposed on the underlying disk. Alternatively, the young ages of off-center clumps may signal inward clump migration, provided this happens efficiently on the order of an orbital timescale.

We present the results from a survey for 12CO emission in 40 luminous sub-millimetre galaxies (SMGs), with 850um fluxes of S850 = 4 - 20 mJy, conducted with the Plateau de Bure Interferometer. We detect 12CO emission in 32 SMGs at z~1.2 - 4.1, including 16 SMGs not previously published. Using multiple 12CO line (J_up =2 - 7) observations, we derive a median spectral line energy distribution for luminous SMGs and use this to estimate a mean gas mass of (5.3 +/- 1.0) \times 10^10 Msun. We report the discovery of a fundamental relationship between 12CO FWHM and 12CO line luminosity in high-redshift starbursts, which we interpret as a natural consequence of the baryon-dominated dynamics within the regions probed by our observations. We use far-infrared luminosities to assess the star-formation efficiency in our SMGs, finding a steepening of the L'CO-LFIR relation as a function of increasing 12CO J_up transition. We derive dynamical masses and molecular gas masses, and use these to determine the redshift evolution of the gas content of SMGs, finding that they do not appear to be significantly more gas rich than less vigorously star-forming galaxies at high redshifts. Finally, we collate X-ray observations, and study the interdependence of gas and dynamical properties of SMGs with their AGN activity and supermassive black hole masses (MBH), finding that SMGs lie significantly below the local M_BH-sigma relation. We conclude that SMGs represent a class of massive, gas-rich ultraluminous galaxies with somewhat heterogeneous properties, ranging from starbursting disc-like systems with L~10^12 L_sun, to the most highly star-forming mergers in the Universe.

We have obtained high-resolution data of the z 2 ring-like, clumpy star-forming galaxy (SFG) ZC406690 using the VLT/SINFONI with AO (in K-band) and in seeing-limited mode (in H- and J-band). Our data includes all of the main strong optical emission lines: [OII], [OIII], Ha, Hb, [NII] and [SII]. We find broad, blueshifted Ha and [OIII] emission line wings in the spectra of the galaxy's massive, star-forming clumps (sigma \sim 85 km s^-1) and even broader wings (up to 70% of the total Ha flux, with sigma \sim 290 km s^-1) in regions spatially offset from the clumps by \sim 2 kpc. The broad emission likely originates from large-scale outflows with mass outflow rates from individual clumps that are 1-8x the SFR of the clumps. Based on emission line ratio diagnostics ([NII]/Ha and [SII]/Ha) and photoionization and shock models, we find that the emission from the clumps is due to a combination of photoionization from the star-forming regions and shocks generated in the outflowing component, with 5-30% of the emission deriving from shocks. In terms of the ionization parameter (6x10^7-10^8 cm/s, based on both the SFR and the O32 ratio), density (local electron densities of 300-1800 cm^-3 in and around the clumps, and ionized gas column densities of 1200-8000 Msol/pc^2), and SFR (10-40 Msol/yr), these clumps more closely resemble nuclear starburst regions of local ULIRGs and dwarf irregulars than HII regions in local galaxies. However, the star-forming clumps are not located in the nucleus as in local starburst galaxies but instead are situated in a ring several kpc from the center of their high-redshift host galaxy, and have an overall disk-like morphology. The two brightest clumps are quite different in terms of their internal properties, energetics and relative ages, and thus we are given a glimpse at two different stages in the formation and evolution of rapidly star-forming giant clumps at high-z.

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

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

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

Herschel/PACS spectroscopy of the luminous infrared galaxies NGC4418 and Arp220 reveals high excitation in H2O, OH, HCN, and NH3. In NGC4418, absorption lines were detected with E_low>800 K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), while in Arp220 the excitation is somewhat lower. While outflow signatures in moderate excitation lines are seen in Arp220 as reported in previous studies, in NGC4418 the lines tracing its outer regions are redshifted relative to the nucleus, suggesting an inflow with Mdot<~12 Msun yr^{-1}. Both galaxies have warm (Tdust>~100 K) nuclear continuum components, together with a more extended component that is much more prominent and massive in Arp220. A chemical dichotomy is found in both sources: on the one hand, the nuclear regions have high H2O abundances, ~10^{-5}, and high HCN/H2O and HCN/NH3 column density ratios of 0.1-0.4 and 2-5, respectively, indicating a chemistry typical of evolved hot cores where grain mantle evaporation has occurred. On the other hand, the high OH abundance, with OH/H2O ratios of ~0.5, indicates the effects of X-rays and/or cosmic rays. The nuclear media have surface brightnesses >~10^{13} Lsun/kpc^2 and are estimated to be thick (N_H>~10^{25} cm^{-2}). While NGC4418 shows weak absorption in H2^{18}O and ^{18}OH, with a ^{16}O-to-^{18}O ratio of >~250-500, the strong absorption of the rare isotopologues in Arp220 indicates ^{16}O-to-^{18}O of 70-130. Further away from the nuclear regions, the H2O abundance decreases to <~10^{-7} and the OH/H2O ratio is 2.5-10. Despite the different scales of NGC4418, Arp220, and Mrk231, preliminary evidence is found for an evolutionary sequence from infall, hot-core like chemistry, and solar oxygen isotope ratio to high velocity outflow, disruption of the hot core chemistry and cumulative high mass stellar processing of 18O.

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

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

We present a measurement of the spatial clustering of submillimetre galaxies (SMGs) at z = 1-3. Using data from the 870 micron LESS survey, we employ a novel technique to measure the cross-correlation between SMGs and galaxies, accounting for the full probability distributions for photometric redshifts of the galaxies. From the observed projected two-point cross-correlation function we derive the linear bias and characteristic dark matter (DM) halo masses for the SMGs. We detect clustering in the cross-correlation between SMGs and galaxies at the > 4 sigma level. For the SMG autocorrelation we obtain r_0 = 7.7 (+1.8,-2.3) h^-1 Mpc, and derive a corresponding DM halo mass of log(M_halo [h^-1 M_sun]) = 12.8 (+0.3,-0.5). Based on the evolution of DM haloes derived from simulations, we show that that the z = 0 descendants of SMGs are typically massive (~2-3 L*) elliptical galaxies residing in moderate- to high-mass groups (log(M_halo [h^-1 M_sun]) = 13.3 (+0.3,-0.5). From the observed clustering we estimate an SMG lifetime of ~100 Myr, consistent with lifetimes derived from gas consumption times and star-formation timescales, although with considerable uncertainties. The clustering of SMGs at z ~ 2 is consistent with measurements for optically-selected quasi-stellar objects (QSOs), supporting evolutionary scenarios linking starbursts and QSOs. Given that SMGs reside in haloes of characteristic mass ~ 6 x 10^12 h^-1 M_sun, we demonstrate that the redshift distribution of SMGs can be described remarkably well by the combination of two effects: the cosmological growth of structure and the evolution of the molecular gas fraction in galaxies. We conclude that the powerful starbursts in SMGs likely represent a short-lived but universal phase in massive galaxy evolution, associated with the transition between cold gas-rich, star-forming galaxies and passively evolving systems. [Abridged]

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

We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor {\alpha}CO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log({\alpha}CO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.

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

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

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.

The zCOSMOS SINFONI project is aimed at studying the physical and kinematical properties of a sample of massive z~1.4-2.5 star-forming galaxies, through SINFONI near-IR integral field spectroscopy (IFS), combined with the multi-wavelength information from the zCOSMOS (COSMOS) survey. The project is based on 1 hour of natural-seeing observations per target, and Adaptive Optics (AO) follow-up for a major part of the sample, which includes 30 galaxies selected from the zCOSMOS/VIMOS spectroscopic survey. This first paper presents the sample selection, and the global physical characterization of the target galaxies from multicolor photometry, i.e., star formation rate (SFR), stellar mass, age, etc. The Halpha integrated properties such as, flux, velocity dispersion, and size, are derived from the natural-seeing observations, while the follow up AO observations will be presented in the next paper of this series. Our sample appears to be well representative of star-forming galaxies at z~2, covering a wide range in mass and SFR. The Halpha integrated properties of the 25 Halpha detected galaxies are similar to those of other IFS samples at the same redshifts. Good agreement is found among the SFRs derived from Halpha luminosity and other diagnostic methods, provided the extinction affecting the Halpha luminosity is about twice that affecting the continuum. A preliminary kinematic analysis, based on the maximum observed velocity difference across the source, and on the integrated velocity dispersion, indicates that the sample splits nearly 50-50 into rotation-dominated and velocity dispersion-dominated galaxies, in good agreement with previous surveys.

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

We present a study of the infrared properties of X-ray selected, moderate luminosity (Lx=10^{42}-10^{44}ergs/s) active galactic nuclei (AGNs) up to z~3, to explore the links between star formation in galaxies and accretion onto their central black holes. We use 100um and 160um fluxes from GOODS-Herschel -the deepest survey yet undertaken by the Herschel telescope- and show that in >94 per cent of cases these fluxes are dominated by the host. We find no evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs at any redshift, suggesting that star-formation is decoupled from nuclear (AGN) activity. The star formation rates of AGN hosts increase strongly with redshift; by a factor of 43 from z<0.1 to z=2-3 for AGNs with the same X-ray luminosities. This increase is consistent with the factor of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (main-sequence) galaxies. Indeed, the average SSFRs of AGN hosts are only marginally (20 per cent) lower than those of main-sequence galaxies, with this small deficit being due to a fraction of AGNs residing in quiescent (low-SSFR) galaxies. We estimate 79+/-10 per cent of moderate AGNs are hosted in main-sequence galaxies, 15+/-7 per cent in quiescent galaxies and <10 per cent in strongly starbursting galaxies. The fractions of all main sequence galaxies at z<2 experiencing a period of moderate nuclear activity is strongly dependent on galaxy stellar mass (Mstars); rising from a few per cent at Mstars~10^{10}Msun to >20 per cent at Mstars>10^{11}Msun. Our results indicate that it is galaxy stellar mass that is most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than violent mergers, suggesting that disk instabilities could be an important AGN feeding mechanism.

[Abridged] In order to understand which process (e.g. galactic winds, cold accretion) is responsible for the cool (T~10^4 K) halo gas around galaxies, we embarked on a program to study the star-formation properties of galaxies selected by their MgII absorption signature in quasar spectra. Specifically, we searched for the H-alpha line emission from galaxies near very strong z=2 MgII absorbers (with rest-frame equivalent width EW>2 \AA) because these could be the sign-posts of outflows or inflows. Surprisingly, we detect H-alpha from only 4 hosts out of 20 sight-lines (and 2 out of the 19 HI-selected sight-lines), despite reaching a star-formation rate (SFR) sensitivity limit of 2.9 M/yr (5-sigma) for a Chabrier initial mass function. This low success rate is in contrast with our z=1 survey where we detected 66%\ (14/21) of the MgII hosts. Taking into account the difference in sensitivity between the two surveys, we should have been able to detect >11.4 of the 20 z=2 hosts whereas we found only 4 galaxies. Interestingly, all the z=2 detected hosts have observed SFR greater than 9 M/yr, well above our sensitivity limit, while at z=1 they all have SFR less than 9 M/yr, an evolution that is in good agreement with the evolution of the SFR main sequence. Moreover, we show that the z=2 undetected hosts are not hidden under the quasar continuum after stacking our data and that they also cannot be outside our surveyed area. Hence, strong MgII absorbers could trace star-formation driven winds in low-mass halos (Mhalo < 10^{10.6} Msun). Alternatively, our results imply that z=2 galaxies traced by strong MgII absorbers do not form stars at a rate expected (3--10 M/yr) for their (halo or stellar) masses, supporting the existence of a transition in accretion efficiency at Mhalo ~ 10^{11} Msun. This scenario can explain both the detections and the non-detections.

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

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

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high redshift galaxies in the optical/near infrared (NIR) domain. This can lead to the contamination of high redshift galaxy searches based on broad band optical/NIR photometry by lower redshift dusty galaxies as both kind of galaxies cannot be distinguished. The contamination rate could be as high as 50%. {This work shows how the far infrared (FIR) domain can help to recognize likely low-z interlopers in an optical/NIR search for high-z galaxies.} We analyse the FIR SEDs of two galaxies proposed as very high redshift ($z>7$) dropout candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160\,$\mu$m, with SPIRE/Herschel at 250, 350 and 500\,$\mu$m and with LABOCA/APEX at 870\,$\mu$m. We find that redshifts $>7$ would imply extreme FIR SEDs (with dust temperatures $>100$\,K and FIR luminosities $>10^{13}$\,$L_{\odot}$). At z$\sim$2, instead, the SEDs of both sources would be compatible with that of typical ULIRGs/SMGs. Considering all the data available for these sources from visible to FIR we re-estimate the redshifts and we find $z\sim$1.6--2.5. Due to the strong spectral breaks observed in these galaxies, standard templates from the literature fail to reproduce the visible-near IR part of the SEDs even when additional extinction is included. These sources resemble strongly dust obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at $z=4$. Galaxies with similar SEDs could contaminate other high redshift surveys.

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.

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

The kinematics of distant galaxies, from z=0.1 to z>2, play a key role in our understanding of galaxy evolution from early times to the present. One of the important parameters is the intrinsic, or local, velocity dispersion of a galaxy, which allows one to quantify the degree of non-circular motions such as pressure support. However, this is difficult to measure because the observed dispersion includes the effects of (often severe) beam smearing on the velocity gradient. Here we investigate four methods of measuring the dispersion that have been used in the literature, to assess their effectiveness at recovering the intrinsic dispersion. We discuss the biasses inherent in each method, and apply them to model disk galaxies in order to determine which methods yield meaningful quantities, and under what conditions. All the mean weighted dispersion estimators are affected by (residual) beam smearing. In contrast, the dispersion recovered by fitting a spatially and spectrally convolved disk model to the data is unbiassed by the beam smearing it is trying to compensate. Because of this, and because the bias it does exhibit depends only on the signal-to-noise, it can be considered reliable. However, at very low signal-to-noise, all methods should be used with caution.

We study the properties of luminous stellar clumps identified in deep, high resolution HST/NIC2 F160W imaging at 1.6um of six z~2 star-forming galaxies with existing near-IR integral field spectroscopy from SINFONI at the VLT. Individual clumps contribute ~0.5%-15% of the galaxy-integrated rest-frame ~5000A emission, with median of about 2%; the total contribution of clump light ranges from 10%-25%. The median intrinsic clump size and stellar mass are ~1kpc and log(Mstar[Msun])~9, in the ranges for clumps identified in rest-UV or line emission in other studies. The clump sizes and masses in the subset of disks are broadly consistent with expectations for clump formation via gravitational instabilities in gas-rich, turbulent disks given the host galaxies' global properties. By combining the NIC2 data with ACS/F814W imaging available for one source, and AO-assisted SINFONI Halpha data for another, we infer modest color, M/L, and stellar age variations within each galaxy. In these two objects, sets of clumps identified at different wavelengths do not fully overlap; NIC2-identified clumps tend to be redder/older than ACS- or Halpha-identified clumps without rest-frame optical counterparts. There is evidence for a systematic trend of older ages at smaller galactocentric radii among the clumps, consistent with scenarios where inward migration of clumps transports material towards the central regions. From constraints on a bulge-like component at radii <1-3kpc, none of the five disks in our sample appears to contain a compact massive stellar core, and we do not discern a trend of bulge stellar mass fraction with stellar age of the galaxy. Further observations are necessary to probe the build-up of stellar bulges and the role of clumps in this process.

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

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

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

We present deep 1.2 mm continuum mapping of a 566 arcmin^2 area within the Lockman Hole North field, previously a target of the Spitzer Wide-area Infrared Extragalactic (SWIRE) survey and extremely deep 20 cm mapping with the Very Large Array, which we have obtained using the Max-Planck millimeter bolometer (MAMBO) array on the IRAM 30 m telescope. After filtering, our full map has an RMS sensitivity ranging from 0.45 to 1.5 mJy/beam, with an average of 0.75 mJy/beam. Using the pixel flux distribution in a map made from our best data, we determine the shape, normalization, and approximate flux density cutoff for 1.2 mm number counts well below our nominal sensitivity and confusion limits. After validating our full dataset through comparison with this map, we successfully detect 41 1.2 mm sources with S/N > 4.0 and S(1.2 mm)\simeq 2-5 mJy. We use the most significant of these detections to directly determine the integral number counts down to 1.8 mJy, which are consistent with the results of the pixel flux distribution analysis. 93% of our 41 individual detections have 20 cm counterparts, 49% have Spitzer/MIPS 24 micron counterparts, and one may have a significant Chandra X-ray counterpart. We resolve \simeq 3% of the cosmic infrared background (CIB) at 1.2 mm into significant detections, and directly estimate a 0.05 mJy faint-end cutoff for the counts that is consistent with the full intensity of the 1.2 mm CIB. The median redshift of our 17 detections with spectroscopic or robust photometric redshifts is z(median)=2.3, and rises to z(median)=2.9 when we include redshifts estimated from the radio/far-infrared spectral index. By using a nearest neighbor and angular correlation function analysis, we find evidence that our S/N>4.0 detections are clustered at the 95% confidence level.

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

We derive the extinction curve towards the Galactic Center from 1 to 19 micron. We use hydrogen emission lines of the minispiral observed by ISO-SWS and SINFONI. The extinction free flux reference is the 2 cm continuum emission observed by the VLA. Towards the inner 14" * 20" we find an extinction of A(2.166 micron)=2.62 +/- 0.11, with a power-law slope of alpha=-2.11 +/- 0.06 shortward of 2.8 micron, consistent with the average near infrared slope from the recent literature. At longer wavelengths, however, we find that the extinction is grayer than shortward of 2.8 micron. We find it is not possible to fit the observed extinction curve with a dust model consisting of pure carbonaceous and silicate grains only, and the addition of composite particles, including ices, is needed to explain the observations. Combining a distance dependent extinction with our distance independent extinction we derive the distance to the GC to be R_0=7.94 +/- 0.65 kpc. Towards Sgr A* (r<0.5") we obtain A_H=4.21 +/- 0.10, A_Ks=2.42 +/- 0.10 and A_L'=1.09 +/- 0.13.

Mass outflows driven by stars and active galactic nuclei are a key element in many current models of galaxy evolution. They may produce the observed black hole-galaxy mass relation and regulate and quench both star formation in the host galaxy and black hole accretion. However, observational evidence of such feedback processes through outflows of the bulk of the star forming molecular gas is still scarce. Here we report the detection of massive molecular outflows, traced by the hydroxyl molecule (OH), in far-infrared spectra of ULIRGs obtained with Herschel-PACS as part of the SHINING key project. In some of these objects the (terminal) outflow velocities exceed 1000 km/s, and their outflow rates (up to ~1200 M_sol/yr) are several times larger than their star formation rates. We compare the outflow signatures in different types of ULIRGs and in starburst galaxies to address the issue of the energy source (AGN or starburst) of these outflows. We report preliminary evidence that ULIRGs with a higher AGN luminosity (and higher AGN contribution to L_IR) have higher terminal velocities and shorter gas depletion time scales. The outflows in the observed ULIRGs are able to expel the cold gas reservoirs from the centres of these objects within ~1E6-1E8 years.

We use deep 70, 100 and 160 um observations taken with PACS, the Photodetector Array Camera and Spectrometer on board of Herschel, as part of the PACS Evolutionary Probe (PEP) guaranteed time, to study the relation between star formation rate and environment at redshift ~ 1 in the GOODS-S and GOODS-N fields. We use the SDSS spectroscopic catalog to build the local analog and study the evolution of the star formation activity dependence on the environment. At z ~ 1 we observe a reversal of the relation between star formation rate and local density, confirming the results based on Spitzer 24 um data. However, due to the high accuracy provided by PACS in measuring the star formation rate also for AGN hosts, we identify in this class of objects the cause for the reversal of the density-SFR relation. Indeed, AGN hosts favor high stellar masses, dense regions and high star formation rates. Without the AGN contribution the relation flattens consistently with respect to the local analog in the same range of star formation rates. As in the local universe, the specific star formation rate anti-correlates with the density. This is due to mass segregation both at high and low redshift. The contribution of AGN hosts does not affect this anti-correlation, since AGN hosts exhibit the same specific star formation rate as star forming galaxies at the same mass. The same global trends and AGN contribution is observed once the relations are studied per morphological type. We study the specific star formation rate vs stellar mass relation in three density regimes. Our data provides an indication that at M/M_{\odot} > 10^{11} the mean specific star formation rate tends to be higher at higher density, while the opposite trend is observed in the local SDSS star forming sample.

We have studied the properties of giant star forming clumps in five z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. The clumps reside in disk regions where the Toomre Q-parameter is below unity, consistent with their being bound and having formed from gravitational instability. Broad H{\alpha}/[NII] line wings demonstrate that the clumps are launching sites of powerful outflows. The inferred outflow rates are comparable to or exceed the star formation rates, in one case by a factor of eight. Typical clumps may lose a fraction of their original gas by feedback in a few hundred million years, allowing them to migrate into the center. The most active clumps may lose much of their mass and disrupt in the disk. The clumps leave a modest imprint on the gas kinematics. Velocity gradients across the clumps are 10-40 km/s/kpc, similar to the galactic rotation gradients. Given beam smearing and clump sizes, these gradients may be consistent with significant rotational support in typical clumps. Extreme clumps may not be rotationally supported; either they are not virialized, or they are predominantly pressure supported. The velocity dispersion is spatially rather constant and increases only weakly with star formation surface density. The large velocity dispersions may be driven by the release of gravitational energy, either at the outer disk/accreting streams interface, and/or by the clump migration within the disk. Spatial variations in the inferred gas phase oxygen abundance are broadly consistent with inside-out growing disks, and/or with inward migration of the clumps.

We report initial results from the far-infrared fine structure line observations of a sample of 44 local starbursts, Seyfert galaxies and infrared luminous galaxies obtained with the PACS spectrometer on board Herschel. We show that the ratio between the far-infrared luminosity and the molecular gas mass, Lfir/MH2, is a much better proxy for the relative brightness of the far-infrared lines than Lfir alone. Galaxies with high Lfir/MH2 ratios tend to have weaker fine structure lines relative to their far-infrared continuum than galaxies with Lfir/MH2 < 80 Lsun/Msun. A deficit of the [C II] 158 micron line relative to Lfir was previously found with the ISO satellite, but now we show for the first time that this is a general aspect of all far-infrared fine structure lines, regardless of their origin in the ionized or neutral phase of the interstellar medium. The Lfir/MH2 value where these line deficits start to manifest is similar to the limit that separates between the two modes of star formation recently found in galaxies on the basis of studies of their gas-star formation relations. Our finding that the properties of the interstellar medium are also significantly different in these regimes provides independent support for the different star forming relations in normal disk galaxies and major merger systems. We use the spectral synthesis code Cloudy to model the emission of the lines. The expected increase of the ionization parameter with Lfir/MH2 can simultaneously explain the line deficits in the [C II], [N II] and [O I] lines.

The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fraction. The galaxies segregate into dispersion-dominated systems with 1<vmax/sigma<3, maximum rotational velocities vmax<200 km/s and disk half-light radii rd=1-3 kpc and rotation-dominated systems with vmax>200 km/s, vmax/sigma>3 and rd=4-8 kpc. For the dispersion-dominated sample, radial pressure gradients partly compensate the gravitational force, reducing the rotational velocities. Including this pressure effect in the MMW model, dispersion-dominated galaxies can be fitted well with spin parameters lf lambda_d=0.03-0.05 for high disk mass fractions of m_d=0.2 and with lambda_d=0.01-0.03 for m_d=0.05. These values are in good agreement with cosmological expectations. For the rotation-dominated sample however pressure effects are small and better agreement with theoretically expected disk spin parameters can only be achieved if the dark halo mass contribution in the visible disk regime (2-3*rd) is smaller than predicted by the MMW model. We argue that these galaxies can still be embedded in standard cold dark matter halos if the halos did not contract adiabatically in response to disk formation. It is shown that the observed high turbulent gas motions of the galaxies are consistent with a Toomre instability parameter Q=1 which is equal to the critical value, expected for gravitational disk instability to be the major driver of turbulence. The dominant energy source of turbulence is then the potential energy of the gas in the disk.

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 low-resolution, rest-frame ~ 5 - 12 micron Spitzer/IRS spectra of two lensed z ~ 2 UV-bright star-forming galaxies, SDSS J120602.09+514229.5 and SDSS J090122.37+181432.3. Using the magnification boost from lensing, we are able to study the physical properties of these objects in greater detail than is possible for unlensed systems. In both targets, we detect strong PAH emission at 6.2, 7.7, and 11.3 microns, indicating the presence of vigorous star formation. For J1206, we find a steeply rising continuum and significant [S IV] emission, suggesting that a moderately hard radiation field is powering continuum emission from small dust grains. The strength of the [S IV] emission also implies a sub-solar metallicity of ~ 0.5 Z_{Sun}, confirming published rest-frame optical measurements. In J0901, the PAH lines have large rest-frame equivalent widths (> 1 micron) and the continuum rises slowly with wavelength, suggesting that any AGN contribution to L_{IR} is insignificant, in contrast to the implications of optical emission-line diagnostics. Using [O III] line flux as a proxy for AGN strength, we estimate that the AGN in J0901 provides only a small fraction of its mid-infrared continuum flux. By combining the detection of [Ar II] with an upper limit on [Ar III] emission, we infer a metallicity of > 1.3 Z_{Sun}. This work highlights the importance of combining rest-frame optical and mid-IR spectroscopy in order to understand the detailed properties of star-forming galaxies at high redshift.

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

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

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

We present ultra-deep mid-IR spectra of 48 infrared-luminous galaxies in the GOODS-South field obtained with the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. These galaxies are selected among faint infrared sources (0.14 - 0.5 mJy at 24 um) in two redshift bins (0.76-1.05 and 1.75-2.4) to sample the major contributors to the cosmic infrared background at the most active epochs. We estimate redshifts for 92% of the sample using PAH and Si absorption features. Only few of these galaxies (5% at z~1 and 12% at z~2) have their total infrared luminosity dominated by emission from AGN. The averaged mid-IR spectra of the z~1 LIRGs and of the z~2 ULIRGs are very similar to the averaged spectrum of local starbursts and HII-like ULIRGs, respectively. We find that 6.2um PAH equivalent widths reach a plateau of ~1 um for L(24 mu) < 1E11 L(sun). At higher luminosities, EW (6.2 mu) anti-correlates with L(24 um). Intriguingly, high-z ULIRGs and SMG lie above the local EW (6.2 um) - L(24 um) relationship suggesting that, at a given luminosity, high-z ULIRGs have AGN contributions to their dust emission lower than those of local counterparts. A quantitative analysis of their morphology shows that most of the luminous IR galaxies have morphologies similar to those of IR-quiet galaxies at the same redshift. All z~2 ULIRGs of our sample are IR-excess BzK galaxies and most of them have L(FIR)/L(1600A) ratios higher than those of starburst galaxies at a given UV slope. The ``IR excess'' (Daddi et al. 2007) is mostly due to strong 7.7 um PAH emission and under-estimation of UV dust extinction. On the basis of the AGN-powered L (6 um) continuum measured directly from the mid-IR spectra, we estimate an average intrinsic X-ray AGN luminosity of L(2-10 keV) = (0.1 +/- 0.6) 1E43 erg/s, a value substantially lower than the prediction by Daddi et al. (2007).

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

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

The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210\mu\ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85\mu\ m or 85-125\mu\m and 125-210\mu\ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions.

We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

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

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

We present Herschel-PACS observations of rest-frame mid-infrared and far-infrared spectral line emissions from two lensed, ultra-luminous infrared galaxies at high redshift: MIPS J142824.0+352619 (MIPS J1428), a starburst-dominated system at z = 1.3, and IRAS F10214+4724 (F10214), a source at z = 2.3 hosting both star-formation and a luminous AGN. We have detected [OI]63 micron and [OIII]52 micron in MIPS J1428, and tentatively [OIII]52 micron in F10214. Together with the recent ZEUS-CSO [CII]158 micron detection in MIPS J1428 we can for the first time combine [OI], [CII] and far-IR (FIR) continuum measurements for photo-dissociation (PDR) modeling of an ultra-luminous (L_IR > 10^12 L_sun) star forming galaxy at the peak epoch of cosmic star formation. We find that MIPS J1428, contrary to average local ULIRGs, does not show a deficit in [OI] relative to FIR. The combination of far-UV flux G_0 and gas density n (derived from the PDR models), as well as the star formation efficiency (derived from CO and FIR) is similar to normal or starburst galaxies, despite the high infrared luminosity of this system. In contrast, F10214 has stringent upper limits on [OIV] and [SIII], and an [OIII]/FIR ratio at least an order of magnitude lower than local starbursts or AGN, similar to local ULIRGs.

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

In this first paper on the results of our Herschel PACS survey of local Ultraluminous Infrared Galaxies (ULIRGs), as part of our SHINING survey of local galaxies, we present far-infrared spectroscopy of Mrk 231, the most luminous of the local ULIRGs, and a type 1 broad absorption line AGN. For the first time in a ULIRG, all observed far-infrared fine-structure lines in the PACS range were detected and all were found to be deficient relative to the far infrared luminosity by 1 - 2 orders of magnitude compared with lower luminosity galaxies. The deficits are similar to those for the mid-infrared lines, with the most deficient lines showing high ionization potentials. Aged starbursts may account for part of the deficits, but partial covering of the highest excitation AGN powered regions may explain the remaining line deficits. A massive molecular outflow, discovered in OH and 18OH, showing outflow velocities out to at least 1400 km/sec, is a unique signature of the clearing out of the molecular disk that formed by dissipative collapse during the merger. The outflow is characterized by extremely high ratios of 18O / 16O suggestive of interstellar medium processing by advanced starbursts.

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

We use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.

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

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