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.

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

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.

Over the past decade increasingly robust estimates of the dense molecular gas content in galaxy populations between redshift 0 and the peak of cosmic galaxy/star formation from redshift 1-3 have become available. This rapid progress has been possible due to the advent of powerful ground-based, and space telescopes for combined study of several millimeter to far-IR, line or continuum tracers of the molecular gas and dust components. The main conclusions of this review are: 1. Star forming galaxies contained much more molecular gas at earlier cosmic epochs than at the present time. 2. The galaxy integrated depletion time scale for converting the gas into stars depends primarily on z or Hubble time, and at a given z, on the vertical location of a galaxy along the star-formation rate versus stellar mass "main-sequence" (MS) correlation. 3. Global rates of galaxy gas accretion primarily control the evolution of the cold molecular gas content and star formation rates of the dominant MS galaxy population, which in turn vary with the cosmological expansion. A second key driver may be global disk fragmentation in high-z, gas rich galaxies, which ties local free-fall time scales to galactic orbital times, and leads to rapid radial matter transport and bulge growth. Third, the low star formation efficiency inside molecular clouds is plausibly set by super-sonic streaming motions, and internal turbulence, which in turn may be driven by conversion of gravitational energy at high-z, and/or by local feedback from massive stars at low-z. 4. A simple 'gas regulator' model is remarkably successful in predicting the combined evolution of molecular gas fractions, star formation rates, galactic winds, and gas phase metallicities.

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.

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare.

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.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

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 study the r31=L'CO(3-2)/L'CO(1-0) luminosity line ratio in a sample of nearby (z < 0.05) galaxies: 25 star-forming galaxies (SFGs) from the xCOLD GASS survey, 36 hard X-ray selected AGN host galaxies from BASS and 37 infrared luminous galaxies from SLUGS. We find a trend for r31 to increase with star-formation efficiency (SFE). We model r31 using the UCL-PDR code and find that the gas density is the main parameter responsible for variation of r31, while the interstellar radiation field and cosmic ray ionization rate play only a minor role. We interpret these results to indicate a relation between SFE and gas density. We do not find a difference in the r31 value of SFGs and AGN host galaxies, when the galaxies are matched in SSFR (<r31>= 0.52 +/- 0.04 for SFGs and <r31> = 0.53 +/- 0.06 for AGN hosts). According to the results of UCL-PDR models, the X-rays can contribute to the enhancement of the CO line ratio, but only for strong X-ray fluxes and for high gas density (nH > 10$^4$ cm-3). We find a mild tightening of the Kennicutt-Schmidt relation when we use the molecular gas mass surface density traced by CO(3-2) (Pearson correlation coefficient R=0.83), instead of the molecular gas mass surface density traced by CO(1-0) (R=0.78), but the increase in correlation is not statistically significant (p-value=0.06). This suggests that the CO(3-2) line can be reliably used to study the relation between SFR and molecular gas for normal SFGs at high redshift, and to compare it with studies of low-redshift galaxies, as is common practice.

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 investigate the fueling and the feedback of nuclear activity in the Seyfert 2 galaxy NGC1068, by studying the distribution and kinematics of molecular gas in the torus and its connections. We use ALMA to image the emission of a set of molecular gas tracers in the circumnuclear disk (CND) and the torus of the galaxy using the CO(2-1), CO(3-2) and HCO+(4-3) lines with spatial resolutions ~0.03"-0.09"(2-6pc). ALMA resolves the CND as an asymmetric ringed disk of D~400pc-size and mass of ~1.4x10^8 Msun. The inner edge of the ring is associated with edge-brightened arcs of NIR polarized emission identified with the working surface of the AGN ionized wind. ALMA proves the existence of a molecular torus of M_torus ~ 3x10^5Msun, which extends over a large range of spatial scales D=10-30pc around the central engine. The new observations evidence the density radial stratification of the torus: the HCO+(4-3) torus, with a full size D=11pc, is a factor of 2-3 smaller than its CO(2-1) and CO(3-2) counterparts, which have full-sizes D=26pc and D=28pc, respectively. The torus is connected to the CND through a network of gas streamers. The kinematics of molecular gas show strong departures from circular motions in the torus, the gas streamers, and the CND. These velocity distortions are interconnected and are part of a 3D outflow that reflects the effects of AGN feedback on the kinematics of molecular gas across a wide range of spatial scales. We conclude that a wide-angle AGN wind launched from the accretion disk is impacting a sizeable fraction of the gas inside the torus (~0.4-0.6 x M_torus). However, a large gas reservoir (~1.2-1.8 x 10^5Msun) close to the equatorial plane of the torus remains unaffected by the AGN wind and can continue fueling the AGN for ~1-4Myr. AGN fueling seems nevertheless thwarted on intermediate scales (15pc < r < 50pc).

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

The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (`feedback') are the main uncertainties in modern cosmological simulations of galaxy formation and evolution. These processes determine the properties of galaxies, but are poorly understood on the $\lesssim$100 pc scale of individual giant molecular clouds (GMCs) resolved in modern galaxy formation simulations. The key question is why the timescale for depleting molecular gas through star formation in galaxies ($t_{\rm dep}\approx 2$ Gyr) exceeds the dynamical timescale of GMCs by two orders of magnitude. Either most of a GMC's mass is converted into stars over many dynamical times, or only a small fraction turns into stars before the GMC is dispersed on a dynamical timescale. Here we report our observation that molecular gas and star formation are spatially de-correlated on GMC scales in the nearby flocculent spiral galaxy NGC300, contrary to their tight correlation on galactic scales. We demonstrate that this de-correlation implies rapid evolutionary cycling between GMCs, star formation, and feedback. We apply a novel statistical method to quantify the evolutionary timeline and find that star formation is regulated by efficient stellar feedback, driving GMC dispersal on short timescales (~1.5 Myr) due to radiation and stellar winds, prior to supernova explosions. This feedback limits GMC lifetimes to about one dynamical timescale (~10 Myr), with integrated star formation efficiencies of only 2-3%. Our findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven lifecycles, that vary with the galactic environment and collectively define how galaxies form stars. Systematic applications of this multi-scale analysis to large galaxy samples will provide key input for a predictive, bottom-up theory of galaxy formation and evolution.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

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.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

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.

We derive two-dimensional dust attenuation maps at $\sim1~\mathrm{kpc}$ resolution from the UV continuum for ten galaxies on the $z\sim2$ Star-Forming Main Sequence (SFMS). Comparison with IR data shows that 9 out of 10 galaxies do not require further obscuration in addition to the UV-based correction, though our sample does not include the most heavily obscured, massive galaxies. The individual rest-frame $V$-band dust attenuation (A$_{\rm V}$) radial profiles scatter around an average profile that gently decreases from $\sim1.8$ mag in the center down to $\sim0.6$ mag at $\sim3-4$ half-mass radii. We use these maps to correct UV- and H$\alpha$-based star-formation rates (SFRs), which agree with each other. At masses $<10^{11}~M_{\rm sun}$, the dust-corrected specific SFR (sSFR) profiles are on average radially constant at a mass-doubling timescale of $\sim300~\mathrm{Myr}$, pointing at a synchronous growth of bulge and disk components. At masses $>10^{11}~M_{\rm sun}$, the sSFR profiles are typically centrally-suppressed by a factor of $\sim10$ relative to the galaxy outskirts. With total central obscuration disfavored, this indicates that at least a fraction of massive $z\sim2$ SFMS galaxies have started their inside-out star-formation quenching that will move them to the quenched sequence. In combination with other observations, galaxies above and below the ridge of the SFMS relation have respectively centrally-enhanced and centrally-suppressed sSFRs relative to their outskirts, supporting a picture where bulges are built due to gas `compaction' that leads to a high central SFR as galaxies move towards the upper envelope of SFMS.

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.

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 introduce xCOLD GASS, a legacy survey providing a census of molecular gas in the local Universe. Building upon the original COLD GASS survey, we present here the full sample of 532 galaxies with CO(1-0) measurements from the IRAM-30m telescope. The sample is mass-selected in the redshift interval $0.01<z<0.05$ from SDSS, and therefore representative of the local galaxy population with M$_{\ast}>10^9$M$_{\odot}$. The CO(1-0) flux measurements are complemented by observations of the CO(2-1) line with both the IRAM-30m and APEX telescopes, HI observations from Arecibo, and photometry from SDSS, WISE and GALEX. Combining the IRAM and APEX data, we find that the CO(2-1) to CO(1-0) luminosity ratio for integrated measurements is $r_{21}=0.79\pm0.03$, with no systematic variations across the sample. The CO(1-0) luminosity function is constructed and best fit with a Schechter function with parameters {$L_{\mathrm{CO}}^* = (7.77\pm2.11) \times 10^9\,\mathrm{K\,km\,s^{-1}\, pc^{2}}$, $\phi^{*} = (9.84\pm5.41) \times 10^{-4} \, \mathrm{Mpc^{-3}}$ and $\alpha = -1.19\pm0.05$}. With the sample now complete down to stellar masses of $10^9$M$_{\odot}$, we are able to extend our study of gas scaling relations and confirm that both molecular gas fraction and depletion timescale vary with specific star formation rate (or offset from the star-formation main sequence) much more strongly than they depend on stellar mass. Comparing the xCOLD GASS results with outputs from hydrodynamic and semi-analytic models, we highlight the constraining power of cold gas scaling relations on models of galaxy formation.

We study the feedback of star formation and nuclear activity on the chemistry of molecular gas in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing if the abundances of key molecular species like ethynyl (C2H), a classical tracer of PDR, change in the different environments of the disk of the galaxy. We have used ALMA to map the emission of the hyperfine multiplet of C2H(N=1-0) and its underlying continuum emission in the central r~35"(2.5kpc)-region of the disk of NGC1068 with a spatial resolution 1.0"x0.7"(50-70pc). We have developed a set of time-dependent chemical models to determine the origin of the C2H gas. A sizeable fraction of the total C2H line emission is detected from the r~1.3kpc starburst (SB) ring. However, the brightest C2H emission originates from a r~200pc off-centered circumnuclear disk (CND), where evidence of a molecular outflow has been previously found in other molecular tracers imaged by ALMA. We also detect significant emission that connects the CND with the outer disk. We derived the fractional abundances of C2H (X(C2H)) assuming LTE conditions. Our estimates range from X(C2H)~a few 10^-8 in the SB ring up to X(C2H)~ a few 10^-7 in the outflow region. PDR models that incorporate gas-grain chemistry are able to account for X(C2H) in the SB ring for moderately dense (n(H2)>10^4 cm^-3) and moderately UV-irradiated gas (UV-field<10xDraine field) in a steady-state regime. However, the high fractional abundances estimated for C2H in the outflow region can only be reached at very early times (T< 10^2-10^3 yr) in models of UV/X-ray irradiated dense gas (n(H2)>10^4-10^5) cm^-3). We interpret that the transient conditions required to fit the high values of X(C2H) in the outflow are likely due to UV/X-ray irradiated non-dissociative shocks associated with the highly turbulent interface between the outflow and the molecular gas in NGC1068.

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 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 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 present Herschel PACS observations of the [CII] 158 micron emission line in a sample of 24 intermediate mass (9<logM$_\ast$/M$_\odot$<10) and low metallicity (0.4< Z/Z$_\odot$<1.0) galaxies from the xCOLD GASS survey. Combining them with IRAM CO(1-0) measurements, we establish scaling relations between integrated and molecular region [CII]/CO(1-0) luminosity ratios as a function of integrated galaxy properties. A Bayesian analysis reveals that only two parameters, metallicity and offset from the star formation main sequence, $\Delta$MS, are needed to quantify variations in the luminosity ratio; metallicity describes the total dust content available to shield CO from UV radiation, while $\Delta$MS describes the strength of this radiation field. We connect the [CII]/CO luminosity ratio to the CO-to-H$_2$ conversion factor and find a multivariate conversion function $\alpha_{CO}$, which can be used up to z~2.5. This function depends primarily on metallicity, with a second order dependence on $\Delta$MS. We apply this to the full xCOLD GASS and PHIBSS1 surveys and investigate molecular gas scaling relations. We find a flattening of the relation between gas mass fraction and stellar mass at logM$_\ast$/M$_\odot$<10. While the molecular gas depletion time varies with sSFR, it is mostly independent of mass, indicating that the low L$_{CO}$/SFR ratios long observed in low mass galaxies are entirely due to photodissociation of CO, and not to an enhanced star formation efficiency.

We present 90 mas (37 pc) resolution ALMA imaging of Arp 220 in the CO (1-0) line and continuum at $\lambda = 2.6$ mm. The internal gas distribution and kinematics of both galactic nuclei are well-resolved for the first time. In the West nucleus, the major gas and dust emission extends out to 0.2\arcsec radius (74 pc); the central resolution element shows a strong peak in the dust emission but a factor 3 dip in the CO line emission. In this nucleus, the dust is apparently optically thick ($\tau_{\rm 2.6mm} \sim1$) at $\lambda = 2.6$ mm with a dust brightness temperature $\sim147$ K. The column of ISM at this nucleus is $\rm N_{H2} \geq 2\times10^{26}$ cm$^{-2}$, corresponding to $\sim$900 gr cm$^{-2}$. The East nucleus is more elongated with radial extent 0.3\arcsec or $\sim111$ pc. The derived kinematics of the nuclear disks provide a good fit to the line profiles, yielding the emissivity distributions, the rotation curves and velocity dispersions. In the West nucleus, there is evidence of a central Keplerian component requiring a central mass of $8\times10^8$ \msun. The intrinsic widths of the emission lines are $\Delta \rm v (FWHM)$ = 250 (West) and 120 (East) \kms. Given the very short dissipation timescales for turbulence ($\lesssim10^5$ yrs), we suggest that the line widths may be due to semi-coherent motions within the nuclear disks. The symmetry of the nuclear disk structures is impressive -- implying the merger timescale is significantly longer than the rotation period of the disks.

We present the ALMA survey of CO(2-1) emission from the 1/5 solar metallicity, Local Group dwarf galaxy NGC 6822. We achieve high (0.9 arcsec ~ 2 pc) spatial resolution while covering large area: four 250 pc x 250 pc regions that encompass ~2/3 of NGC 6822's star formation. In these regions, we resolve ~150 compact CO clumps that have small radii (~2-3 pc), narrow line width (~1 km/s), and low filling factor across the galaxy. This is consistent with other recent studies of low metallicity galaxies, but here shown with a 15 times larger sample. At parsec scales, CO emission correlates with 8 micron emission better than with 24 micron emission and anti-correlates with Halpha, so that PAH emission may be an effective tracer of molecular gas at low metallicity. The properties of the CO clumps resemble those of similar-size structures in Galactic clouds except of slightly lower surface brightness and CO-to-H2 ratio ~1-2 times the Galactic value. The clumps exist inside larger atomic-molecular complexes with masses typical for giant molecular cloud. Using dust to trace H2 for the entire complex, we find CO-to-H2 to be ~20-25 times the Galactic value, but with strong dependence on spatial scale and variations between complexes that may track their evolutionary state. The H2-to-HI ratio is low globally and only mildly above unity within the complexes. The SFR-to-H2 ratio is ~3-5 times higher in the complexes than in massive disk galaxies, but after accounting for the bias from targeting star-forming regions, we conclude that the global molecular gas depletion time may be as long as in massive disk galaxies.

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 an analysis of the relationship between the CO-H$_{2}$ conversion factor ($\alpha_{\rm CO}$) and total mass surface density ($\Sigma_{\rm tot}$) in star-forming galaxies at $z < 1.5$. Our sample, which is drawn from the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) and the CO Legacy Database for GASS (COLD GASS), includes 'normal,' massive star-forming galaxies that dominate the evolution of the cosmic star formation rate (SFR) at this epoch and probe the $\Sigma_{\rm tot}$ regime where the strongest variation in $\alpha_{\rm CO}$ is observed. We constrain $\alpha_{\rm CO}$ via existing CO observations, measurements of the star formation rate, and an assumed molecular gas depletion time ($t_{\rm dep}$=$M_{\rm gas}$/SFR) --- the latter two of which establish the total molecular gas mass independent of the observed CO luminosity. For a broad range of adopted depletion times, we find that $\alpha_{\rm CO}$ is independent of total mass surface density, with little deviation from the canonical Milky Way value. This runs contrary to a scenario in which $\alpha_{\rm CO}$ decreases as surface density increases within the extended clouds of molecular gas that potentially fuel clumps of star formation in $z\sim1$ galaxies, similar to those observed in local ULIRGs. Instead, our results suggest that molecular gas, both at $z\sim0$ and $z\sim1$, is primarily in the form of self-gravitating molecular clouds. While CO observations suggest a factor of $\sim3$ reduction in the average molecular gas depletion time between $z \sim 0$ and $z\sim1$, we find that, for typical galaxies, the structure of molecular gas and the process of star formation at $z \sim 1$ is otherwise remarkably similar to that observed in local star-forming systems.

We make publicly available a catalog of calibrated environmental measures for galaxies in the five 3D-HST/CANDELS deep fields. Leveraging the spectroscopic and grism redshifts from the 3D-HST survey, multi wavelength photometry from CANDELS, and wider field public data for edge corrections, we derive densities in fixed apertures to characterize the environment of galaxies brighter than $JH_{140} < 24$ mag in the redshift range $0.5<z<3.0$. By linking observed galaxies to a mock sample, selected to reproduce the 3D-HST sample selection and redshift accuracy, each 3D-HST galaxy is assigned a probability density function of the host halo mass, and a probability that is a central or a satellite galaxy. The same procedure is applied to a $z=0$ sample selected from SDSS. We compute the fraction of passive central and satellite galaxies as a function of stellar and halo mass, and redshift, and then derive the fraction of galaxies that were quenched by environment specific processes. Using the mock sample, we estimate that the timescale for satellite quenching is $t_{\rm quench} \sim 2-5$ Gyr; longer at lower stellar mass or lower redshift, but remarkably independent of halo mass. This indicates that, in the range of environments commonly found within the 3D-HST sample, satellites are quenched by exhaustion of their gas reservoir in absence of cosmological accretion. We find that the quenching times can be separated into a delay phase during which satellite galaxies behave similarly to centrals at fixed stellar mass, and a phase where the star formation rate drops rapidly ($\sim 0.4-0.6$ Gyr), as shown previously at $z=0$. We conclude that this scenario requires satellite galaxies to retain a large reservoir of multi-phase gas upon accretion, even at high redshift, and that this gas sustains star formation for the long quenching times observed.

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 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 use spectra from the ALFALFA, GASS and COLD GASS surveys to quantify variations in the mean atomic and molecular gas mass fractions throughout the SFR-M* plane and along the main sequence (MS) of star-forming galaxies. Although galaxies well below the MS tend to be undetected in the Arecibo and IRAM observations, reliable mean atomic and molecular gas fractions can be obtained through a spectral stacking technique. We find that the position of galaxies in the SFR-M* plane can be explained mostly by their global cold gas reservoirs as observed in the HI line, with in addition systematic variations in the molecular-to-atomic ratio and star formation efficiency. When looking at galaxies within +/-0.4 dex of the MS, we find that as stellar mass increases, both atomic and molecular gas mass fractions decrease, stellar bulges become more prominent, and the mean stellar ages increase. Both star formation efficiency and molecular-to-atomic ratios vary little for massive main sequence galaxies, indicating that the flattening of the MS is due to the global decrease of the cold gas reservoirs of galaxies rather than to bottlenecks in the process of converting cold atomic gas to stars.

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 have used the Atacama Large Millimeter Array (ALMA) to map the emission of the CO(6-5) molecular line and the 432 {\mu}m continuum emission from the 300 pc-sized circumnuclear disk (CND) of the nearby Seyfert 2 galaxy NGC 1068 with a spatial resolution of ~4 pc. These observations spatially resolve the CND and, for the first time, image the dust emission, the molecular gas distribution, and the kinematics from a 7-10 pc-diameter disk that represents the submillimeter counterpart of the putative torus of NGC 1068. We fitted the nuclear spectral energy distribution of the torus using ALMA and near and mid-infrared (NIR/MIR) data with CLUMPY models. The mass and radius of the best-fit solution for the torus are both consistent with the values derived from the ALMA data alone: Mgas_torus=(1+-0.3)x10^5 Msun and Rtorus=3.5+-0.5 pc. The dynamics of the molecular gas in the torus show non-circular motions and enhanced turbulence superposed on the rotating pattern of the disk. The kinematic major axis of the CO torus is tilted relative to its morphological major axis. By contrast with the nearly edge-on orientation of the H2O megamaser disk, we have found evidence suggesting that the molecular torus is less inclined (i=34deg-66deg) at larger radii. The lopsided morphology and complex kinematics of the torus could be the signature of the Papaloizou-Pringle instability, long predicted to likely drive the dynamical evolution of active galactic nuclei (AGN) tori.

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.

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 present the first observations of H$^{13}$CN$(1-0)$, H$^{13}$CO$^+(1-0)$ and SiO$(2-1)$ in NGC\,6240, obtained with the IRAM PdBI. Combining a Markov Chain Monte Carlo (MCMC) code with Large Velocity Gradient (LVG) modelling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass$-$luminosity conversion factors. We find $\sim10^{10}M_\odot$ of dense molecular gas in cold, dense clouds ($T_{\rm k}\sim10$\,K, $n_{{\rm H}_2}\sim10^6$\,cm$^{-3}$) with a volume filling factor $<0.002$, embedded in a shock heated molecular medium ($T_{\rm k}\sim2000$\,K, $n_{{\rm H}_2}\sim10^{3.6}$\,cm$^{-3}$), both surrounded by an extended diffuse phase ($T_{\rm k}\sim200$\,K, $n_{{\rm H}_2}\sim10^{2.5}$\,cm$^{-3}$). We derive a global $\alpha_{\rm CO}=1.5^{7.1}_{1.1}$ with gas masses $\log_{10}\left(M / [M_\odot]\right)=10.1_{10.0}^{10.8}$, dominated by the dense gas. We also find $\alpha_{\rm HCN} = 32^{89}_{13}$, which traces the cold, dense gas. The [$^{12}$C]/[$^{13}$C] ratio is only slightly elevated ($98^{230}_{65}$), contrary to the very high [CO]/[$^{13}$CO] ratio (300-500) reported in the literature. However, we find very high [HCN]/[H$^{13}$CN] and [HCO$^+$]/[H$^{13}$CO$^+$] abundance ratios $(300^{500}_{200})$ which we attribute to isotope fractionation in the cold, dense clouds.

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 Karl G. Jansky Very Large Array observations of the CO J=1-0 transition in a sample of four $z\sim2$ main sequence galaxies. These galaxies are in the blue sequence of star-forming galaxies at their redshift, and are part of the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) which imaged them in CO J=3-2. Two galaxies are imaged here at high signal-to-noise, allowing determinations of their disk sizes, line profiles, molecular surface densities, and excitation. Using these and published measurements, we show that the CO and optical disks have similar sizes in main-sequence galaxies, and in the galaxy where we can compare CO J=1-0 and J=3-2 sizes we find these are also very similar. Assuming a Galactic CO-to-H$_2$ conversion, we measure surface densities of $\Sigma_{mol}\sim1200$ M$_\odot$pc$^{-2}$ in projection and estimate $\Sigma_{mol}\sim500-900$ M$_\odot$pc$^{-2}$ deprojected. Finally, our data yields velocity-integrated Rayleigh-Jeans brightness temperature line ratios $r_{31}$ that are approximately unity. In addition to the similar disk sizes, the very similar line profiles in J=1-0 and J=3-2 indicate that both transitions sample the same kinematics, implying that their emission is coextensive. We conclude that in these two main sequence galaxies there is no evidence for significant excitation gradients or a large molecular reservoir that is diffuse or cold and not involved in active star-formation. We suggest that $r_{31}$ in very actively star-forming galaxies is likely an indicator of how well mixed the star formation activity and the molecular reservoir are.

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.

Galaxy evolution scenarios predict that the feedback of star formation and nuclear activity (AGN) can drive the transformation of gas-rich spiral mergers into ULIRGs, and, eventually, lead to the build-up of QSO/elliptical hosts. We study the role that star formation and AGN feedback have in launching and maintaining the molecular outflows in two starburst-dominated advanced mergers, NGC1614 and IRAS17208-0014, by analyzing the distribution and kinematics of their molecular gas reservoirs. We have used the PdBI array to image with high spatial resolution (0.5"-1.2") the CO(1-0) and CO(2-1) line emissions in NGC1614 and IRAS17208-0014, respectively. The velocity fields of the gas are analyzed and modeled to find the evidence of molecular outflows in these sources and characterize the mass, momentum and energy of these components. While most (>95%) of the CO emission stems from spatially-resolved (~2-3kpc-diameter) rotating disks, we also detect in both mergers the emission from high-velocity line wings that extend up to +-500-700km/s, well beyond the estimated virial range associated with rotation and turbulence. The kinematic major axis of the line wing emission is tilted by ~90deg in NGC1614 and by ~180deg in IRAS17208-0014 relative to their respective rotating disk major axes. These results can be explained by the existence of non-coplanar molecular outflows in both systems. In stark contrast with NGC1614, where star formation alone can drive its molecular outflow, the mass, energy and momentum budget requirements of the molecular outflow in IRAS17208-0014 can be best accounted for by the existence of a so far undetected (hidden) AGN of L_AGN~7x10^11 L_sun. The geometry of the molecular outflow in IRAS17208-0014 suggests that the outflow is launched by a non-coplanar disk that may be associated with a buried AGN in the western nucleus.

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.

Most present-day galaxies with stellar masses $\geq10^{11}$ solar masses show no ongoing star formation and are dense spheroids. Ten billion years ago, similarly massive galaxies were typically forming stars at rates of hundreds solar masses per year. It is debated how star formation ceased, on which timescales, and how this "quenching" relates to the emergence of dense spheroids. We measured stellar mass and star-formation rate surface density distributions in star-forming galaxies at redshift 2.2 with $\sim1$ kiloparsec resolution. We find that, in the most massive galaxies, star formation is quenched from the inside out, on timescales less than 1 billion years in the inner regions, up to a few billion years in the outer disks. These galaxies sustain high star-formation activity at large radii, while hosting fully grown and already quenched bulges in their cores.

We present the analysis of HST $J$- and $H$-band imaging for 29 galaxies on the star-forming main sequence at $z\sim2$, which have Adaptive Optics VLT SINFONI integral field spectroscopy from our SINS/zC-SINF program. The SINFONI H$\alpha$ data resolve the on-going star-formation and the ionized gas kinematics on scales of $1-2$ kpc; the near-IR images trace the galaxies' rest-frame optical morphologies and distributions of stellar mass in old stellar populations at a similar resolution. The global light profiles of most galaxies show disk-like properties well described by a single S\'ersic profile with $n\sim1$, with only $\sim15%$ requiring a high $n>3$ S\'ersic index, all more massive than $10^{10}M_\odot$. In bulge+disk fits, about $40%$ of galaxies have a measurable bulge component in the light profiles, with $\sim15%$ showing a substantial bulge-to-total ratio $B/T\ge0.3$. This is a lower limit to the frequency of $z\sim2$ massive galaxies with a developed bulge component in stellar mass because it could be hidden by dust and/or outshined by a thick actively star-forming disk component. The galaxies' rest-optical half-light radii range between $1-7$ kpc, with a median of 2.1 kpc, and lie slightly above the size-mass relation at these epochs reported in the literature. This is attributed to differences in sample selection and definitions of size and/or mass measurements. The $(u-g)_{rest}$ color gradient and scatter within individual $z\sim2$ massive galaxies with $\ge10^{11}M_\odot$ are as high as in $z=0$ low-mass, late-type galaxies, and are consistent with the high star-formation rates of massive $z\sim2$ galaxies being sustained at large galactocentric distances.

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.

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 report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

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.

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 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 a detailed analysis of ALMA Bands 7 and 9 data of CO, HCO+, HCN and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~ 200 pc circumnuclear disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and whether the different line intensity ratios we find within the CND as well as between the CND and the SB ring are due to excitation effects (gas density and temperature differences) or to a different chemistry. We estimate the column densities of each species in Local Thermodynamic Equilibrium (LTE). We then compute large one-dimensional non-LTE radiative transfer grids (using RADEX) by using first only the CO transitions, and then all the available molecules in order to constrain the densities, temperatures and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN and CO/CS column density ratios. RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. The LTE, RADEX and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND.

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 investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.

Most massive galaxies are thought to have formed their dense stellar cores at early cosmic epochs. However, cores in their formation phase have not yet been observed. Previous studies have found galaxies with high gas velocity dispersions or small apparent sizes but so far no objects have been identified with both the stellar structure and the gas dynamics of a forming core. Here we present a candidate core in formation 11 billion years ago, at z=2.3. GOODS-N-774 has a stellar mass of 1.0x10^11 Msun, a half-light radius of 1.0 kpc, and a star formation rate of 90[+45-20]Msun/yr. The star forming gas has a velocity dispersion 317+-30 km/s, amongst the highest ever measured. It is similar to the stellar velocity dispersions of the putative descendants of GOODS-N-774, compact quiescent galaxies at z~2 and giant elliptical galaxies in the nearby Universe. Galaxies such as GOODS-N-774 appear to be rare; however, from the star formation rate and size of the galaxy we infer that many star forming cores may be heavily obscured, and could be missed in optical and near-infrared surveys.

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 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 report ALMA observations of CO(3-2) emission in the Seyfert 2 double-barred galaxy NGC1433, at the unprecedented spatial resolution of 0.5"=24 pc. Our aim is to probe AGN feeding and feedback phenomena through the morphology and dynamics of the gas inside the central kpc. The CO map, which covers the whole nuclear region (nuclear bar and ring), reveals a nuclear gaseous spiral structure, inside the nuclear ring encircling the nuclear stellar bar. This gaseous spiral is well correlated with the dusty spiral seen in Hubble Space Telescope images. The nuclear spiral winds up in a pseudo-ring at 200 pc radius, which might correspond to the inner ILR. Continuum emission is detected at 0.87 mm only at the very centre, and its origin is more likely thermal dust emission than non-thermal emission from the AGN. It might correspond to the molecular torus expected to exist in this Seyfert 2 galaxy. The HCN(4-3) and HCO+(4-3) lines were observed simultaneously, but only upper limits are derived, with a ratio to the CO(3-2) line lower than 1/60 at 3sigma, indicating a relatively low abundance of very dense gas. The kinematics of the gas over the nuclear disk reveal rather regular rotation only slightly perturbed by streaming motions due to the spiral; the primary and secondary bars are too closely aligned with the galaxy major or minor axis to leave a signature in the projected velocities. Near the nucleus, there is an intense high-velocity CO emission feature redshifted to 200 km/s (if located in the plane), with a blue-shifted counterpart, at 2" (100 pc) from the centre. While the CO spectra are quite narrow in the centre, this wide component is interpreted as an outflow involving a molecular mass of 3.6 10^6 Mo and a flow rate 7 Mo/yr. The flow could be in part driven by the central star formation, but is mainly boosted by the AGN through its wind or radio jets.

We report ALMA observations of CO(3-2) emission in the Seyfert 1 galaxy NGC 1566, at a spatial resolution of 25 pc. Our aim is to investigate the morphology and dynamics of the gas inside the central kpc, and to probe nuclear fueling and feedback phenomena. NGC 1566 has a nuclear bar of 1.7 kpc radius and a conspicuous grand design spiral starting from this radius. The ALMA field of view, of diameter 0.9 kpc, lies well inside the nuclear bar and reveals a molecular trailing spiral structure from 50 to 300~pc in size, which is contributing to fuel the nucleus, according to its negative gravity torques. The spiral starts with a large pitch angle from the center and then winds up in a pseudo-ring at the inner Lindblad resonance (ILR) of the nuclear bar. This is the first time that a trailing spiral structure is clearly seen driving the gas inwards inside the ILR ring of the nuclear bar. This phenomenon shows that the massive central black hole has a significant dynamical influence on the gas, triggering its fueling. The gaseous spiral is well correlated with the dusty spiral seen through extinction in HST images, and also with a spiral feature emitting 0.87mm continuum. This continuum emission must come essentially from cold dust heated by the interstellar radiation field. The HCN(4-3) and HCO+(4-3) lines were simultaneously mapped and detected in the nuclear spiral. The HCO+(4-3) line is 3 times stronger than the HCN(4-3), as expected when star formation excitation dominates over active galactic nucleus (AGN) heating. The CO(3-2)/HCO+(4-3) integrated intensity ratio is \sim 100. The molecular gas is in remarkably regular rotation, with only slight non-circular motions at the periphery of the nuclear spiral arms. These perturbations are quite small, and no outflow nor AGN feedback is detected.

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

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

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.

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.

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

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.

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

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

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

Using observations from the GASS and COLD GASS surveys and complementary data from SDSS and GALEX, we investigate the nature of variations in gas depletion time observed across the local massive galaxy population. The large and unbiased COLD GASS sample allows us to assess the relative importance of galaxy interactions, bar instabilities, morphologies and the presence of AGN in regulating star formation efficiency. Both the H2 mass fraction and depletion time vary as a function of the distance of a galaxy from the main sequence in the SFR-M* plane. The longest gas depletion times are found in below-main sequence bulge-dominated galaxies that are either gas-poor, or else on average less efficient than disk-dominated galaxy at converting into stars any cold gas they may have. We find no link between AGN and these long depletion times. The galaxies undergoing mergers or showing signs of morphological disruptions have the shortest molecular gas depletion times, while those hosting strong stellar bars have only marginally higher global star formation efficiencies as compared to matched control samples. Our interpretation is that depletion time variations are caused by changes in the ratio between the gas mass traced by the CO(1-0) observations, and the gas mass in high density star-forming cores, with interactions, mergers and bar instabilities able to locally increase pressure and raise the ratio of efficiently star-forming gas to CO-detected gas. Building a sample representative of the local massive galaxy population, we derive a global Kennicutt-Schmidt relation of slope 1.18+/-0.24, and observe structure within the scatter around this relation, with galaxies having low (high) stellar mass surface densities lying systematically above (below) the mean relation, suggesting that gas surface density is not the only parameter driving the global star formation ability of a galaxy.

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.

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 present new ~1" resolution data of the dense molecular gas in the central 50-100 pc of four nearby Seyfert galaxies. PdBI observations of HCN and, in 2 of the 4 sources, simultaneously HCO+ allow us to carefully constrain the dynamical state of the dense gas surrounding the AGN. Analysis of the kinematics shows large line widths of 100-200 km/s FWHM that can only partially arise from beam smearing of the velocity gradient. The observed morphological and kinematic parameters (dimensions, major axis position angle, red and blue channel separation, and integrated line width) are well reproduced by a thick disk, where the emitting dense gas has a large intrinsic dispersion (20-40 km/s), implying that it exists at significant scale heights (25-30% of the disk radius). To put the observed kinematics in the context of the starburst and AGN evolution, we estimate the Toomre Q parameter. We find this is always greater than the critical value, i.e. Q is above the limit such that the gas is stable against rapid star formation. This is supported by the lack of direct evidence, in these 4 Seyfert galaxies, for on-going star formation close around the AGN. Instead, any current star formation tends to be located in a circumnuclear ring. We conclude that the physical conditions are indeed not suited to star formation within the central ~100 pc.

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 analyze the role of bars in the build-up of central mass concentrations in massive, disk galaxies. Our parent sample consists of 3757 face-on disk galaxies with redshifts between 0.01 and 0.05, selected from the seventh Data Release of the Sloan Digital Sky Survey. 1555 galaxies with bars are identified using position angle and ellipticity profiles of the $i$-band light. We compare the ratio of the specific star formation rate measured in the 1-3 kpc central region of the galaxy to that measured for the whole galaxy. Galaxies with strong bars have centrally enhanced star formation; the degree of enhancement depends primarily on the ellipticity of the bar, and not on the size of the bar or on the mass or structure of the host galaxy. The fraction of galaxies with strong bars is highest at stellar masses greater than $3 \times 10^{10} M_{\odot}$, stellar surface densities less than $3 \times 10^8 M_{\odot}$ and concentration indices less than 2.5. In this region of parameter space, galaxies with strong bars either have enhanced central star formation rates, or star formation that is {\em suppressed} compared to the mean. This suggests that bars may play a role in the eventual quenching of star formation in galaxies. Only 50% of galaxies with strongly concentrated star formation have strong bars, indicating that other processes such as galaxy interactions also induce central star-bursts. We also find that the ratio of the size of the bar to that of the disk depends mainly on the colour of the galaxy, suggesting that the growth and destruction of bars are regulated by gas accretion, as suggested by simulations.

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.

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.

We compare the semi-analytic models of galaxy formation of Fu et al. (2010), which track the evolution of the radial profiles of atomic and molecular gas in galaxies, with gas fraction scaling relations derived from the COLD GASS survey (Saintonge et al 2011). The models provide a good description of how condensed baryons in galaxies with gas are partitioned into stars, atomic and molecular gas as a function of galaxy stellar mass and surface density. The models do not reproduce the tight observed relation between stellar surface density and bulge-to-disk ratio for this population. We then turn to an analysis of the"quenched" population of galaxies without detectable cold gas. The current implementation of radio-mode feedback in the models disagrees strongly with the data. In the models, gas cooling shuts down in nearly all galaxies in dark matter halos above a mass of 10**12 M_sun. As a result, stellar mass is the observable that best predicts whether a galaxy has little or no neutral gas. In contrast, our data show that quenching is largely independent of stellar mass. Instead, there are clear thresholds in bulge-to-disk ratio and in stellar surface density that demarcate the location of quenched galaxies. We propose that processes associated with bulge formation play a key role in depleting the neutral gas in galaxies and that further gas accretion is suppressed following the formation of the bulge, even in dark matter halos of low mass.

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 have obtained long-slit spectra of 174 star-forming galaxies with stellar masses greater than 10^10 M_\odot from the GALEX Arecibo SDSS (GASS) survey. These galaxies have both HI and H_2 mass measurements. The average metallicity profile is strikingly flat out to R_90, the radius enclosing 90% of the r-band light. Metallicity profiles which decline steadily with radius are found primarily for galaxies in our sample with low stellar mass (Log(M_*)<10.2), concentration, and/or mean stellar mass density. Beyond ~R_90, however, around 10 percent of the galaxies in our sample exhibit a sharp downturn in metallicity. Remarkably, we find that the magnitude of the outer metallicity drop is well correlated with the total HI content of the galaxy (measured as f_HI=M_HI/M_*). We examine the radial profiles of stellar population ages and star formation rate densities, and conclude that the galaxies with largest outer metallicity drops are actively growing their stellar disks, with mass doubling times across the whole disk only one third as long as a typical GASS galaxy. We also describe a correlation between local stellar mass density and metallicity, which is valid across all galaxies in our sample. We argue that much of the recent stellar mass growth at the edges of these galaxies can be linked to the accretion or radial transport of relatively pristine gas from beyond the galaxies' stellar disks.

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.

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.

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

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

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

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.

As part of an extensive study of the physical properties of active galactic nuclei (AGN) we report high spatial resolution near-IR integral-field spectroscopy of the narrow-line region (NLR) and coronal-line region (CLR) of seven Seyfert galaxies. These measurements elucidate for the first time the two-dimensional spatial distribution and kinematics of the recombination line Br{\gamma} and high-ionization lines [Sivi], [Alix] and [Caviii] on scales <300 pc from the AGN. The observations reveal kinematic signatures of rotation and outflow in the NLR and CLR. The spatially resolved kinematics can be modeled as a combination of an outflow bicone and a rotating disk coincident with the molecular gas. High-excitation emission is seen in both components, suggesting it is leaking out of a clumpy torus. While NGC 1068 (Seyfert 2) is viewed nearly edge-on, intermediate-type Seyferts are viewed at intermediate angles, consistent with unified schemes. A correlation between the outflow velocity and the molecular gas mass in r<30 pc indicates that the accumulation of gas around the AGN increases the collimation and velocity of the outflow. The outflow rate is 2-3 orders of magnitude greater than the accretion rate, implying that the outflow is mass-loaded by the surrounding interstellar medium (ISM). In half of the observed AGN the kinetic power of the outflow is of the order of the power required by two-stage feedback models to be thermally coupled to the ISM and match the M-{\sigma}* relation. In these objects the radio jet is clearly interacting with the ISM, indicative of a link between jet power and outflow power.

Narrow-Line Seyfert 1 (NLS1) galaxies show extreme properties with respect to the other Seyfert galaxies. Indeed, they are thought to be accreting at Eddington rates and to possess low mass black holes. Therefore, they may represent a key class of objects for understanding the co-evolution of black holes and their host galaxies. We propose that NLS1s represent a class of AGN in which the black hole growth is, and has always been, dominated by secular evolution. Firstly, by looking at the NLS1 host galaxy properties in the literature, we show that the evolution of NLS1s is presently driven by secular processes, much more so than for Broad-Line Seyfert 1s (BLS1s). Secondly, we study the bulges of NLS1 and BLS1 galaxies. Our results demonstrate that NLS1 host bulges are pseudo-bulges and are statistically different from BLS1 bulges. This difference points to the particular importance of secular processes in the past evolution of their hosts. We build on this result to understand the implications on their evolution and the duration of their duty cycle. We show that NLS1s are not necessarily in a special phase of black hole growth and that several Gyr are required for their black hole masses to become similar to BLS1s. Finally, in the light of our results, we discuss the location of NLS1 galaxies on the MBH-\sigma plane and speculate about the connection between the NLS1 galaxy properties and their black hole spin.

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

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.

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 study the relation between molecular gas and star formation in a volume-limited sample of 222 galaxies from the COLD GASS survey, with measurements of the CO(1-0) line from the IRAM 30m telescope. The galaxies are at redshifts 0.025<z<0.05 and have stellar masses in the range 10.0<log(M*/Msun)<11.5. The IRAM measurements are complemented by deep Arecibo HI observations and homogeneous SDSS and GALEX photometry. A reference sample that includes both UV and far-IR data is used to calibrate our estimates of star formation rates from the seven optical/UV bands. The mean molecular gas depletion timescale, tdep(H2), for all the galaxies in our sample is 1 Gyr, however tdep(H2) increases by a factor of 6 from a value of ~0.5 Gyr for galaxies with stellar masses of 10^10 Msun to ~3 Gyr for galaxies with masses of a few times 10^11 Msun. In contrast, the atomic gas depletion timescale remains contant at a value of around 3 Gyr. This implies that in high mass galaxies, molecular and atomic gas depletion timescales are comparable, but in low mass galaxies, molecular gas is being consumed much more quickly than atomic gas. The strongest dependences of tdep(H2) are on the stellar mass of the galaxy (parameterized as log tdep(H2)= (0.36+/-0.07)(log M* - 10.70)+(9.03+/-0.99)), and on the specific star formation rate. A single tdep(H2) versus sSFR relation is able to fit both "normal" star-forming galaxies in our COLD GASS sample, as well as more extreme starburst galaxies (LIRGs and ULIRGs), which have tdep(H2) < 10^8 yr. Normal galaxies at z=1-2 are displaced with respect to the local galaxy population in the tdep(H2) versus sSFR plane and have molecular gas depletion times that are a factor of 3-5 times longer at a given value of sSFR due to their significantly larger gas fractions.

We are conducting COLD GASS, a legacy survey for molecular gas in nearby galaxies. Using the IRAM 30m telescope, we measure the CO(1-0) line in a sample of ~350 nearby (D=100-200 Mpc), massive galaxies (log(M*/Msun)>10.0). The sample is selected purely according to stellar mass, and therefore provides an unbiased view of molecular gas in these systems. By combining the IRAM data with SDSS photometry and spectroscopy, GALEX imaging and high-quality Arecibo HI data, we investigate the partition of condensed baryons between stars, atomic gas and molecular gas in 0.1-10L* galaxies. In this paper, we present CO luminosities and molecular hydrogen masses for the first 222 galaxies. The overall CO detection rate is 54%, but our survey also uncovers the existence of sharp thresholds in galaxy structural parameters such as stellar mass surface density and concentration index, below which all galaxies have a measurable cold gas component but above which the detection rate of the CO line drops suddenly. The mean molecular gas fraction MH2/M* of the CO detections is 0.066+/-0.039, and this fraction does not depend on stellar mass, but is a strong function of NUV-r colour. Through stacking, we set a firm upper limit of MH2/M*=0.0016+/-0.0005 for red galaxies with NUV-r>5.0. The average molecular-to-atomic hydrogen ratio in present-day galaxies is 0.3, with significant scatter from one galaxy to the next. The existence of strong detection thresholds in both the HI and CO lines suggests that "quenching" processes have occurred in these systems. Intriguingly, atomic gas strongly dominates in the minority of galaxies with significant cold gas that lie above these thresholds. This suggests that some re-accretion of gas may still be possible following the quenching event.

We analyse new spatially resolved integral field spectroscopic H- and K-band data at a resolution of 0.3arcsec (100pc) and re-analyse interferometric CO(2-1) line observations of the prototypical merging system Arp220. We find that the majority of the K-band luminosity is due to a 10Myr old starburst, with a significant contribution from an underlying ca 1Gyr old stellar population, and a small contribution from stars less than 8Myr old. The Calzetti (2000) reddening law provides the best fit to photometric datapoints spanning 0.45um -- 2.12um. Furthermore, estimates of the bolometric luminosity from IRAS fluxes in conjunction with our stellar population analysis indicate that we observe less than 10% of the emitted K-band light. The stellar and CO(2-1) kinematic centre of the western nucleus coincides with the compact hot dust emission, indicating that the latter marks the centre of the gravitational potential. In the eastern nucleus, the CO(2-1) data are well matched by a model in which the gas orbits around the peak of the dust emission. This, and the similarity of the K-band tracer kinematics, shows that despite the irregular morphology, the eastern nucleus is also a kinematically coherent structure. Comparison of the extinction map and EW_CO and EW_Br-gamma maps indicates that the lower half of the eastern nucleus is significantly more extincted than the upper half, suggesting that the lower half is buried in the larger-scale gas disk.

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 analyse subarcsecond resolution interferometric CO line data for twelve sub-millimetre-luminous (S850um > 5mJy) galaxies with redshifts between 1 and 3, presenting new data for four of them. Morphologically and kinematically most of the twelve systems appear to be major mergers. Five of them are well-resolved binary systems, and seven are compact or poorly resolved. Of the four binary systems for which mass measurements for both separate components can be made, all have mass ratios of 1:3 or closer. Furthermore, comparison of the ratio of compact to binary systems with that observed in local ULIRGs indicates that at least a significant fraction of the compact SMGs must also be late-stage mergers. In addition, the dynamical and gas masses we derive are most consistent with the lower end of the range of stellar masses published for these systems, favouring cosmological models in which SMGs result from mergers. These results all point to the same conclusion, that likely most of the bright SMGs with L_IR > 5x10e12L_sun are major mergers.

We present spatially resolved integral field spectroscopic K-band data at a resolution of 0.13" (60pc) and interferometric CO(2-1) line observations of the prototypical merging system NGC6240. Despite the clear rotational signature, the stellar kinematics in the two nuclei are dominated by dispersion. We use Jeans modelling to derive the masses and the mass-to-light ratios of the nuclei. Combining the luminosities with the spatially resolved Br-gamma equivalent width shows that only 1/3 of the K-band continuum from the nuclei is associated with the most recent star forming episode; and that less than 30% of the system's bolometric luminosity and only 9% of its stellar mass is due to this starburst. The star formation properties, calculated from typical merger star formation histories, demonstrate the impact of different assumptions about the star formation history. The properties of the nuclei, and the existence of a prominent old stellar population, indicate that the nuclei are remnants of the progenitor galaxies' bulges.

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

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

We report new observations of the galaxy UGC8802 obtained through GASS, the GALEX Arecibo SDSS Survey, which show this galaxy to be in a remarkable evolutionary state. UGC8802 (GASS35981) is a disk galaxy with stellar mass M*=2x10^10 Msolar which appears to contain an additional 2.1x10^10 Msolar of HI gas. New millimeter observations with the IRAM 30m telescope indicate a molecular gas mass only a tenth this large. Using deep long-slit spectroscopy, we examine the spatially resolved star formation rate and metallicity profiles of GASS35981 for clues to its history. We find that the star formation surface density in this galaxy is low (Sigma_SFR=0.003 Msolar/yr/kpc^2) and that the star formation is spread remarkably evenly across the galaxy. The low molecular gas masses measured in our three IRAM pointings are largely consistent with the total star formation measured within the same apertures. Our MMT long-slit spectrum reveals a sharp drop in metallicity in the outer disk of GASS35981. The ratio of current star formation rate to existing stellar mass surface density in the outer disk is extremely high, implying that all the stars must have formed within the past ~1Gyr. At current star formation rates, however, GASS35981 will not consume its HI reservoir for another 5-7 Gyr. Despite its exceptionally large gas fraction for a galaxy this massive, GASS35981 has a regular rotation curve and exhibits no sign of a recent interaction or merger. We speculate that GASS35981 may have acquired its gas directly from the inter-galactic medium, and that it and other similar galaxies identified in the GASS survey may provide rare local glimpses of gas accretion processes that were more common during the prime epoch of disk galaxy formation at z~1.

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

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

We use measurements of the HI content, stellar mass and star formation rates in ~190 massive galaxies with stellar masses greater than 10^10 Msun, obtained from the Galex Arecibo SDSS Survey (GASS) described in Paper I (Catinella et al. 2010) to explore the global scaling relations associated with the bin-averaged ratio of the star formation rate over the HI mass, which we call the HI-based star formation efficiency (SFE). Unlike the mean specific star formation rate, which decreases with stellar mass and stellar mass surface density, the star formation efficiency remains relatively constant across the sample with a value close to SFE = 10^-9.5 yr^-1 (or an equivalent gas consumption timescale of ~3 Gyr). Specifically, we find little variation in SFE with stellar mass, stellar mass surface density, NUV-r color and concentration. We interpret these results as an indication that external processes or feedback mechanisms that control the gas supply are important for regulating star formation in massive galaxies. An investigation into the detailed distribution of SFEs reveals that approximately 5% of the sample shows high efficiencies with SFE > 10^-9 yr^-1, and we suggest that this is very likely due to a deficiency of cold gas rather than an excess star formation rate. Conversely, we also find a similar fraction of galaxies that appear to be gas-rich for their given specific star-formation rate, although these galaxies show both a higher than average gas fraction and lower than average specific star formation rate. Both of these populations are plausible candidates for "transition" galaxies, showing potential for a change (either decrease or increase) in their specific star formation rate in the near future. We also find that 36+/-5% of the total HI mass density and 47+/-5% of the total SFR density is found in galaxies with stellar mass greater than 10^10 Msun. [abridged]

Using the cosmological baryonic accretion rate and normal star formation efficiencies, we present a very simple model for star-forming galaxies (SFGs) that accounts for the mass and redshift dependencies of the SFR-Mass and Tully-Fisher relations from z=2 to the present. The time evolution follows from the fact that each modelled galaxy approaches a steady state where the SFR follows the (net) cold gas accretion rate. The key feature of the model is a halo mass floor M_{min}~10^{11} below which accretion is quenched in order to simultaneously account for the observed slopes of the SFR-Mass and Tully-Fischer relations. The same successes cannot be achieved via a star-formation threshold (or delay) nor by varying the SF efficiency or the feedback efficiency. Combined with the mass ceiling for cold accretion due to virial shock heating, the mass floor M_{min} explains galaxy "downsizing", where more massive galaxies formed earlier and over a shorter period of time. It turns out that the model also accounts for the observed galactic baryon and gas fractions as a function of mass and time, and the cosmic SFR density from z~6 to z=0, which are all resulting from the mass floor M_{min}. The model helps to understand that it is the cosmological decline of accretion rate that drives the decrease of cosmic SFR density between z~2 and z=0 and the rise of the cosmic SFR density allows us to put a constraint on our main parameter M_{min}~10^{11} solar masses. Among the physical mechanisms that could be responsible for the mass floor, we view that photo-ionization feedback (from first in-situ hot stars) lowering the cooling efficiency is likely to play a large role.

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.

We describe several projects addressing the growth of galaxies and massive black holes, for which adaptive optics is mandatory to reach high spatial resolution but is also a challenge due to the lack of guide stars and long integrations. In each case kinematics of the stars and gas, derived from integral field spectroscopy, plays a key role. We explain why deconvolution is not an option, and that instead the PSF is used to convolve a physical model to the required resolution. We discuss the level of detail with which the PSF needs to be known, and the ways available to derive it. We explain how signal-to-noise can limit the resolution achievable and show there are many science cases that require high, but not necessarily diffraction limited, resolution. Finally, we consider what requirements astrometry and photometry place on adaptive optics performance and design.

NGC 2992 is an intermediate Seyfert 1 galaxy showing outflows on kilo parsec scales which might be due either to AGN or starburst activity. We therefore aim at investigating its central region for a putative starburst in the past and its connection to the AGN and the outflows. Observations were performed with the adaptive optics near infrared integral field spectrograph SINFONI on the VLT, complemented by longslit observations with ISAAC on the VLT, as well as N- and Q-band data from the Spitzer archive. The spatial and spectral resolutions of the SINFONI data are 50 pc and 83 km/s, respectively. The field of view of 3" x 3" corresponds to 450 pc x 450 pc. Br_gamma equivalent width and line fluxes from PAHs were compared to stellar population models to constrain the age of the putative recent star formation. A simple geometric model of two mutually inclined disks and an additional cone to describe an outflow was developed to explain the observed complex velocity field in H_2 1-0S(1). The morphologies of the Br_gamma and the stellar continuum are different suggesting that at least part of the Br_gamma emission comes from the AGN. This is confirmed by PAH emission lines at 6.2 micron and 11.2 micron and the strength of the silicon absorption feature at 9.7 micron, which point to dominant AGN activity with a relatively minor starburst contribution. We find a starburst age of 40 Myr - 50 Myr from Br_gamma line diagnostics and the radio continuum; ongoing star formation can be excluded. Both the energetics and the timescales indicate that the outflows are driven by the AGN rather than the starburst. The complex velocity field observed in H_2 1-0S(1) in the central 450 pc can be explained by the superposition of the galaxy rotation and an outflow.

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

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

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.

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

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

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.

This paper is part of a multi-species survey of line emission from the molecular gas in the circum-nuclear disk (CND) of the Seyfert 2 galaxy NGC1068. Single-dish observations have provided evidence that the abundance of silicon monoxide(SiO) in the CND of NGC1068 is enhanced by 3-4 orders of magnitude with respect to the values typically measured in quiescent molecular gas in the Galaxy. We aim at unveiling the mechanism(s) underlying the SiO enhancement. We have imaged with the IRAM Plateau de Bure interferometer the emission of the SiO(2-1) and CN(2--1) lines in NGC1068 at 150pc and 60pc spatial resolution, respectively. We have also obtained complementary IRAM 30m observations of HNCO and methanol (CH3OH) lines. SiO is detected in a disk of 400pc size around the AGN. SiO abundances in the CND of (1-5)xE-09 are about 1-2 orders of magnitude above those measured in the starburst ring. The overall abundance of CN in the CND is high: (0.2-1)xE-07. The abundances of SiO and CN are enhanced at the extreme velocities of gas associated with non-circular motions close to the AGN (r<70pc). Abundances measured for CN and SiO, and the correlation of CN/CO and SiO/CO ratios with hard X-ray irradiation, suggest that the CND of NGC1068 has become a giant X-ray dominated region (XDR). The extreme properties of molecular gas in the circum-nuclear molecular disk of NGC1068 result from the interplay between different processes directly linked to nuclear activity. Whereas XDR chemistry offers a simple explanation for CN and SiO in NGC1068, the relevance of shocks deserves further scrutiny. The inclusion of dust grain chemistry would help solve the controversy regarding the abundances of other molecular species, like HCN, which are under-predicted by XDR models.

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

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

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

We present CO observations of 9 ULIRGs at z~2 with S(24\mu m)>1mJy, previously confirmed with the mid-IR spectra in the Spitzer First Look Survey. All targets are required to have accurate redshifts from Keck/GEMINI near-IR spectra. Using the Plateau de Bure millimeter-wave Interferometer (PdBI) at IRAM, we detect CO J(3-2) [7 objects] or J(2-1) [1 object] line emission from 8 sources with integrated intensities Ic ~(5-9)sigma. The CO detected sources have a variety of mid-IR spectra, including strong PAH, deep silicate absorption and power-law continuum, implying that these molecular gas rich objects at z~2 could be either starbursts or dust obscured AGNs. The measured line luminosity L'[CO] is (1.28-3.77)e+10[K km/s pc^2]. The averaged molecular gas mass M(H2) is 1.7e+10Msun, assuming CO-to-H2 conversion factor of 0.8Msun/[K km/s pc^2]. Three sources (33%) -- MIPS506, MIPS16144 & MIPS8342 -- have double peak velocity profiles. The CO double peaks in MIPS506 and MIPS16144 show spatial separations of 45kpc and 10.9kpc, allowing the estimates of the dynamical masses of 3.2e+11*sin^(-2)(i)Msun and 5.4e+11*sin^{-2}(i)Msun respectively. The implied gas fraction, M(gas)/M(dyn), is 3% and 4%, assuming an average inclination angle. Finally, the analysis of the HST/NIC2 images, mid-IR spectra and IR SED revealed that most of our sources are mergers, containing dust obscured AGNs dominating the luminosities at (3-6)um. Together, these results provide some evidence suggesting SMGs, bright 24um z~2 ULIRGs and QSOs could represent three different stages of a single evolutionary sequence, however, a complete physical model would require much more data, especially high spatial resolution spectroscopy.

Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.

We present high spatial resolution (0.4", ~3.5 kpc) PdBI interferometric data on three ultra-luminous infrared galaxies (ULIRGs) at z~2: two submillimetre galaxies and one submillimetre faint star forming radio galaxy. The three galaxies have been ro- bustly detected in CO rotational transitions, either 12CO(J=4-3) or 12CO(J=3-2), allowing their sizes and gas masses to be accurately constrained. These are the highest spatial resolution observations observed to date (by a factor of ~2) for intermediate-excitation CO emission in z~2 ULIRGs. The galaxies appear extended over several resolution elements, having a mean radius of 3.7 kpc. High-resolution (0.3") combined MERLIN-VLA observations of their radio continua allow an analysis of the star formation behaviour of these galaxies, on comparable spatial scales to that of the CO observations. This 'matched beam' approach sheds light on the spatial distribution of both molecular gas and star formation, and we can therefore calculate accurate star formation rates and gas surface densities: this allows us to place the three systems in the context of a Kennicutt-Schmidt (KS)-style star formation law. We find a difference in size between the CO and radio emission regions, and as such we suggest that using the spatial extent of the CO emission region to estimate the surface density of star formation may lead to error. This size difference also causes the star formation efficiencies within systems to vary by up to a factor of 5. We also find, with our new accurate sizes, that SMGs lie significantly above the KS relation, indicating that stars are formed more efficiently in these extreme systems than in other high-z star forming galaxies.

We present results from infrared spectroscopic projects that aim to test the relation between the mass of a black hole, M_BH, and the velocity dispersion of the stars in its host-galaxy bulge. We demonstrate that near-infrared, high-resolution spectroscopy assisted by adaptive optics is key in populating the high-luminosity end of the relation. We show that the velocity dispersions of mid-infrared, high-ionization lines originating from gas in the narrow-line region of the active galactic nucleus follow the same relation. This result provides a way of inferring M_BH estimates for the cosmologically significant population of obscured, type-2 AGN that can be applicable to data from spectrographs on the next generation infrared telescopes.

We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large program that is gathering high quality HI-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05, selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5-5%) is reached. This paper presents the first Data Release, consisting of ~20% of the final GASS sample. We use this data set to explore the main scaling relations of HI gas fraction with galaxy structure and NUV-r colour. A large fraction (~60%) of the galaxies in our sample are detected in HI. We find that the atomic gas fraction decreases strongly with stellar mass, stellar surface mass density and NUV-r colour, but is only weakly correlated with galaxy bulge-to-disk ratio (as measured by the concentration index of the r-band light). We also find that the fraction of galaxies with significant (more than a few percent) HI decreases sharply above a characteristic stellar surface mass density of 10^8.5 Msun kpc^-2. The fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of the key goals of GASS is to identify and quantify the incidence of galaxies that are transitioning between the blue, star-forming cloud and the red sequence of passively-evolving galaxies. Likely transition candidates can be identified as outliers from the mean scaling relations between gas fraction and other galaxy properties. [abridged]

Our adaptive optics observations of nearby AGN at spatial resolutions as small as 0.085arcsec show strong evidence for recent, but no longer active, nuclear star formation. We begin by describing observations that highlight two contrasting methods by which gas can flow into the central tens of parsecs. Gas accumulation in this region will inevitably lead to a starburst, and we discuss the evidence for such events. We then turn to the impact of stellar evolution on the further inflow of gas by combining a phenomenological approach with analytical modelling and hydrodynamic simulations. These complementary perspectives paint a picture in which all the processes are ultimately regulated by the mass accretion rate into the central hundred parsecs, and the ensuing starburst that occurs there. The resulting supernovae delay accretion by generating a starburst wind, which leaves behind a clumpy interstellar medium. This provides an ideal environment for slower stellar outflows to accrete inwards and form a dense turbulent disk on scales of a few parsecs. Such a scenario may resolve the discrepancy between the larger scale structure seen with adaptive optics and the small scale structure seen with VLTI.

We present the SINS survey with SINFONI of high redshift galaxies. With 80 objects observed and 63 detected, SINS is the largest survey of spatially resolved gas kinematics, morphologies, and physical properties of star-forming galaxies at z~1-3. We describe the selection of the targets, the observations, and the data reduction. We then focus on the "SINS Halpha sample" of 62 rest-UV/optically-selected sources at 1.3<z<2.6 for which we targeted primarily the Halpha and [NII] emission lines. Only 30% of this sample had previous near-IR spectroscopic observations. As a whole, the SINS Halpha sample covers a reasonable representation of massive log(M*/Msun)>~10 star-forming galaxies at z~1.5-2.5, with some bias towards bluer systems compared to pure K-selected samples due to the requirement of secure optical redshift. The sample spans two orders of magnitude in stellar mass and in absolute and specific star formation rates, with median values of approximately log(M*/Msun) = 10.5, 70 Msun/yr, and 3/Gyr. The ionized gas distribution and kinematics are spatially resolved on scales ranging from 1.5 kpc for adaptive optics assisted observations to typically 4-5 kpc for seeing-limited data. The Halpha morphologies tend to be irregular and/or clumpy. About one-third are rotation-dominated yet turbulent disks, another third comprises compact and velocity dispersion-dominated objects, and the remaining galaxies are clear interacting/merging systems; the fraction of rotation-dominated systems increases among the more massive part of the sample. The Halpha luminosities and equivalent widths suggest on average roughly twice higher dust attenuation towards the HII regions relative to the bulk of the stars, and comparable current and past-averaged star formation rates. [Abridged]

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

High signal-to-noise, representative spectra of star-forming galaxies at z~2, obtained via stacking, reveal a high-velocity component underneath the narrow H-alpha and [NII] emission lines. When modeled as a single Gaussian, this broad component has FWHM > 1500 km/s; when modeled as broad wings on the H-alpha and [NII] features, it has FWHM > 500 km/s. This feature is preferentially found in the more massive and more rapidly star-forming systems, which also tend to be older and larger galaxies. We interpret this emission as evidence of either powerful starburst-driven galactic winds or active supermassive black holes. If galactic winds are responsible for the broad emission, the observed luminosity and velocity of this gas imply mass outflow rates comparable to the star formation rate. On the other hand, if the broad line regions of active black holes account for the broad feature, the corresponding black holes masses are estimated to be an order of magnitude lower than those predicted by local scaling relations, suggesting a delayed assembly of supermassive black holes with respect to their host bulges.

Using the IRAM Plateau-de-Bure Interferometer (PdBI) we have searched for the upper fine structure line of neutral carbon (CI(2->1), nu_rest = 809 GHz) and 12CO(J=7->6) (nu_rest=806GHz) towards the submillimetre galaxies (SMGs) GN20 (SMM J123711.9+622212, z = 4.055) and GN20.2 (SMM J123708.8+622202, z = 4.051). The far-infrared (FIR) continuum is detected at 8sigma significance in GN20, with a flux density of S_1.8mm = 1.9+-0.2mJy, while no continuum is detected in GN20.2. Both sources are statistically undetected in both CI and CO(7-6) lines; we derive line luminosity limits for both CI and CO of L'<~2x10^10K km s^-1 pc^2. Assuming carbon excitation temperatures of T_ex = 30K (the galaxies' measured dust temperatures), we infer CI mass limits of M_CI < 5.4x10^6 Msun (GN20) and M_CI < 6.8x10^6 Msun (GN20.2). The derived CI abundance limits are <1.8x10^-5 for GN20 and <3.8x10^-5 for GN20.2 implying that the systems have Milky Way level neutral carbon enrichment (X[CI]/X[H2]) or lower, similar to high-redshift carbon-detected systems (at 5x10^-5) but about 50 times less than the neutral carbon enrichment of local starburst galaxies. Observations of GN20 and GN20.2 in high-resolution MERLIN+VLA radio maps of GOODS-N are used to further constrain the sizes and locations of active regions. We conclude that the physical gas properties of young rapidly evolving systems like GN20 and GN20.2 are likely significantly different than starburst/ULIRG environments in the local Universe yet similar to z~2 SMGs. Unless gravitationally amplified examples can be found, observations of galaxies like GN20 will require the order of magnitude increase in sensitivity of the Atacama Large Millimetre Array (ALMA) to constrain their CI and high-J CO content, despite the fact that they are the brightest systems at z~4.

We present spatially resolved distributions and kinematics of the stars and molecular gas in the central 320pc of NGC1097. The stellar continuum confirms the previously reported 3-arm spiral pattern extending into the central 100pc. The stellar kinematics and the gas distribution imply this is a shadowing effect due to extinction by gas and dust in the molecular spiral arms. The molecular gas kinematics show a strong residual (i.e. non-circular) velocity, which is manifested as a 2-arm kinematic spiral. Linear models indicate that this is the line-of-sight velocity pattern expected for a density wave in gas that generates a 3-arm spiral morphology. We estimate the inflow rate along the arms. Using hydrodynamical models of nuclear spirals, we show that when deriving the accretion rate into the central region, outflow in the disk plane between the arms has to be taken into account. For NGC1097, despite the inflow rate along the arms being ~1.2Msun/yr, the net gas accretion rate to the central few tens of parsecs is much smaller. The numerical models indicate that the inflow rate could be as little as ~0.06Msun/yr. This is sufficient to generate recurring starbursts, similar in scale to that observed, every 20-150Myr. The nuclear spiral represents a mechanism that can feed gas into the central parsecs of the galaxy, with the gas flow sustainable for timescales of a Gigayear.

We report the results from a deep HST NICMOS H-band imaging survey of 28 z < 0.3 QSOs from the Palomar-Green (PG) sample. This program is part of QUEST (Quasar / ULIRG Evolution STudy) and complements a similar set of data on 26 highly-nucleated ULIRGs presented in Paper I. Our analysis indicates that the fraction of QSOs with elliptical hosts is higher among QSOs with undetected far-infrared (FIR) emission, small infrared excess, and luminous hosts. The hosts of FIR-faint QSOs show a tendency to have less pronounced merger-induced morphological anomalies and larger QSO-to-host luminosity ratios on average than the hosts of FIR-bright QSOs, consistent with late-merger evolution from FIR-bright to FIR-faint QSOs. The spheroid sizes and total host luminosities of the radio-quiet PG QSOs in our sample are statistically indistinguishable from the ULIRG hosts presented in Paper I, while those of radio-loud PG QSOs are systematically larger and more luminous. ULIRGs and PG QSOs with elliptical hosts fall near, but not exactly on, the fundamental plane of inactive spheroids. We confirm the systematic trend noted in Paper I for objects with small (< 2 kpc) spheroids to be up to ~1 mag. brighter than inactive spheroids. The host colors and wavelength dependence of their sizes support the idea that these deviations are due at least in part to non-nuclear star formation. However, the amplitudes of these deviations does not depend on host R-H colors. Taken at face value (i.e., no correction for extinction or the presence of a young stellar population), the H-band spheroid-host luminosities imply BH masses ~5 -- 200 x 10^7 M_sun and sub-Eddington mass accretion rates for both QSOs and ULIRGs. These results are compared with published BH mass estimates derived from other methods. (abridged)

We present high-resolution CO(1-0) and CO(2-1) maps of the LINER 2 galaxy NGC1961. This galaxy is unusual among late-type (Sc) disk galaxies in having a very large radial extent and inferred dynamical mass. We propose a head-on collision scenario to explain the perturbed morphology of this galaxy- both the off-centered rings and the inflated radius. This scenario is supported by the detection of a steep velocity gradient in the CO(1-0) map at the position of a southwest peak in radio continuum and near-infrared emission. This peak would represent the remnant of the disrupting companion. We use numerical models to demonstrate the plausibility of the scenario. While ram pressure stripping could in principle be important for shocking the atomic gas and produce the striking head-tail morphology, the non detection of this small galaxy group in X-ray emission suggests that any hot intragroup medium has too low a density. The collision with the companion of mass ratio 1:4 produces a lopsided wave, which accounts for the sharp boundary observed in the atomic gas on the southern side. We argue that NGC1961 has a strongly warped disk, which gives the false impression of a nearly face-on system; the main disk is actually more edge-on, and this error in the true inclination has led to the surprisingly high dynamical mass for a morphologically late-type galaxy. In addition, the outward propagating ring artificially enlarges the disk. The collision de-stabilizes the inner disk and can provide gas inflow to the active nucleus.

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)

In a sample of local active galactic nuclei studied at a spatial resolution on the order of 10 pc we show that the interstellar medium traced by the molecular hydrogen v=1-0 S(1) 2.1um line forms a geometrically thick, clumpy disk. The kinematics of the molecular gas reveals general rotation, although an additional significant component of random bulk motion is required by the high local velocity dispersion. The size scale of the typical gas disk is found to have a radius of ~30 pc with a comparable vertical height. Within this radius the average gas mass is estimated to be ~10^7 Msun based on a typical gas mass fraction of 10%, which suggests column densities of Nh ~ 5x10^23 cm^-2. Extinction of the stellar continuum within this same region suggest lower column densities of Nh ~ 2x10^22 cm^-2, indicating that the gas distribution on these scales is dominated by dense clumps. In half of the observed Seyfert galaxies this lower column density is still great enough to obscure the AGN at optical/infrared wavelengths. We conclude, based on the spatial distribution, kinematics, and column densities that the molecular gas observed is spatially mixed with the nuclear stellar population and is likely to be associated with the outer extent of any smaller scale nuclear obscuring structure. Furthermore, we find that the velocity dispersion of the molecular gas is correlated with the star formation rate per unit area, suggesting a link between the two phenomena, and that the gas surface density follows known "Schmidt-Kennicutt" relations. The molecular/dusty structure on these scales may be dynamic since it is possible that the velocity dispersion of the gas, and hence the vertical disk height, is maintained by a short, massive inflow of material into the nuclear region and/or by intense, short-lived nuclear star formation.

We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 from Halpha line emission. The galaxies are selected from the larger sample of the SINS survey, based on the prominence of ordered rotational motions with respect to more complex merger induced dynamics. The quality of the data allows us to carefully select systems with kinematics dominated by rotation, and to model the gas dynamics across the whole galaxy using suitable exponential disk models. We obtain a good correlation between the dynamical mass and the stellar mass, finding that large gas fractions Mgas~M*) are required to explain the difference between the two quantities. We use the derived stellar mass and maximum rotational velocity Vmax from the modeling to construct for the first time the stellar mass Tully-Fisher relation at z ~ 2.2. The relation obtained shows a slope similar to what is observed at lower redshift, but we detect an evolution of the zero point. We find that at z ~ 2.2 there is an offset in log(M*) for a given rotational velocity of 0.41+-0.11 with respect to the local Universe. This result is consistent with the predictions of the latest N-body/hydrodynamical simulations of disk formation and evolution, which invoke gas accretion onto the forming disk in filaments and cooling flows. This scenario is in agreement with other dynamical evidence from SINS, where gas accretion from the halo is required to reproduce the observed properties of a large fraction of the z ~ 2 galaxies.

In these proceedings, we summarize recent results from our "SINS" VLT/SINFONI integral-field survey, focusing on the 52 detected UV/optically-selected star-forming galaxies at z~2. Our H-alpha emission-line imaging and kinematic data of these systems illustrates that a substantial fraction (> 1/3) of these galaxies are large, rotating disks and that these disks are clumpy, thick, and forming stars rapidly. We compare these systems to local disk scaling relations and find that the backbones of these relations are already in place at z~2. Detailed analysis of the large disks in our sample provides strong evidence that this population cannot result from a merger-dominated formation history and instead must be assembled by the smooth but rapid inflow of gas along filaments. These systems will then secularly evolve from clump-dominated disks to bulge-dominated disks on short timescales, a phenomenon that is observed in our SINS observations and is consistent with predictions from numerical simulations. These results provide new and exciting insights into the formation of bulge-dominated galaxies in the local Universe.

We present a review on bar pattern speeds from preliminary results in the context of the Nuclei of Galaxies (NUGA) project. The large variety of molecular circumnuclear morphologies found in NUGA is a challenging result that implies the refinement of current dynamical models of galaxies.

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

We report the first direct observations of neutral, molecular gas streaming in the nucleus of NGC1068 on scales of <30 pc using SINFONI near-infrared integral field spectroscopy. At a resolution of 0.075", the flux map of 2.12 $\mu$m 1-0 S(1) molecular hydrogen emission around the nucleus in the central arcsec reveals two prominent linear structures leading to the AGN from the north and south. The kinematics of the gas in these features are dominated by non-circular motions and indicate that material is streaming towards the nucleus on highly elliptical or parabolic trajectories whose orientations are compatible with that of the disk plane of the galaxy. We interpret the data as evidence for fueling of gas to the central region. The radial transport rate from ~30 pc to a few parsec from the nucleus is ~15 M$_\sun$ yr$^{-1}$. One of the infalling clouds lies directly in front of the central engine. We interpret it as a tidally disrupted streamer that forms the optically thick outerpart of an amorphous clumpy molecular/dusty structure which contributes to the nuclear obscuration.

Recent observations of UV-/optically selected, massive star forming galaxies at z~2 indicate that the baryonic mass assembly and star formation history is dominated by continuous rapid accretion of gas and internal secular evolution, rather than by major mergers. We use the Millennium Simulation to build new halo merger trees, and extract halo merger fractions and mass accretion rates. We find that even for halos not undergoing major mergers the mass accretion rates are plausibly sufficient to account for the high star formation rates observed in z~2 disks. On the other hand, the fraction of major mergers in the Millennium Simulation is sufficient to account for the number counts of submillimeter galaxies (SMGs), in support of observational evidence that these are major mergers. When following the fate of these two populations in the Millennium Simulation to z=0, we find that subsequent mergers are not frequent enough to convert all z~2 turbulent disks into elliptical galaxies at z=0. Similarly, mergers cannot transform the compact SMGs/red sequence galaxies at z~2 into observed massive cluster ellipticals at z=0. We argue therefore, that secular and internal evolution must play an important role in the evolution of a significant fraction of z~2 UV-/optically and submillimeter selected galaxy populations.

We present Ha integral field spectroscopy of well resolved, UV/optically selected z~2 star-forming galaxies as part of the SINS survey with SINFONI on the ESO VLT. Our laser guide star adaptive optics and good seeing data show the presence of turbulent rotating star forming rings/disks, plus central bulge/inner disk components, whose mass fractions relative to total dynamical mass appears to scale with [NII]/Ha flux ratio and star formation age. We propose that the buildup of the central disks and bulges of massive galaxies at z~2 can be driven by the early secular evolution of gas-rich proto-disks. High redshift disks exhibit large random motions. This turbulence may in part be stirred up by the release of gravitational energy in the rapid cold accretion flows along the filaments of the cosmic web. As a result dynamical friction and viscous processes proceed on a time scale of <1 Gyr, at least an order of magnitude faster than in z~0 disk galaxies. Early secular evolution thus drives gas and stars into the central regions and can build up exponential disks and massive bulges, even without major mergers. Secular evolution along with increased efficiency of star formation at high surface densities may also help to account for the short time scales of the stellar buildup observed in massive galaxies at z~2.

We present the first use of the Gemini North laser guide star adaptive optics (LGS AO) system and an integral field unit (IFU) to measure the stellar velocity dispersion of the host of a luminous quasar. The quasar PG1426+015 (z=0.086) was observed with the Near-Infrared Integral Field Spectrometer (NIFS) on the 8m Gemini North telescope in the H-band as part of the Science Verification phase of the new ALTAIR LGS AO system. The NIFS IFU and LGS AO are well suited for host studies of luminous quasars because one can achieve a large ratio of host to quasar light. We have measured the stellar velocity dispersion of PG1426+015 from 0.1'' to 1'' (0.16 kpc to 1.6 kpc) to be 217+/-15 km/s based on high signal-to-noise ratio measurements of Si I, Mg I, and several CO bandheads. This new measurement is a factor of four more precise than a previous measurement obtained with long-slit spectroscopy and good, natural seeing, yet was obtained with a shorter net integration time. We find that PG1426+015 has a velocity dispersion that places it significantly above the M-sigma relation of quiescent galaxies and lower-luminosity active galactic nuclei with black hole masses estimated from reverberation mapping. We discuss several possible explanations for this discrepancy that could be addressed with similar observations of a larger sample of luminous quasars.

We present Spitzer rest-frame mid-infrared spectroscopy of twelve z~2 mm-bright type 1 QSOs, selected from unlensed and lensed QSO samples and covering a range of AGN optical luminosities L_5100=10^45 to 10^47 erg/s. On top of the AGN continuum, we detect PAH emission from luminous star formation in nine objects individually as well as in the composite spectrum for the full sample. PAH luminosity and rest frame far-infrared luminosity correlate and extend the similar correlation for lower luminosity local QSOs. This provides strong evidence for intense star formation in the hosts of these mm-bright QSOs, sometimes exceeding 1000 Msun/yr and dominating their rest frame far-infrared emission. The PAH-based limit on star formation rates is lower for luminous z~2 QSOs that are not preselected for their mm emission. Partly dependent on systematic changes of the AGN dust covering factor and the dust spectral energy distribution of the AGN proper, the spectral energy distributions of mm-faint high-z QSOs may be AGN dominated out to rest frame far-infrared wavelengths. Towards the most luminous high-z QSOs, there is a flattening of the relation between star formation and AGN luminosity that is observed for lower redshift QSOs. No QSO in our sample has a PAH-measured star formation rate in excess of 3000 Msun/yr.

We present the far-IR, millimeter, and radio photometry as well as optical and near-IR spectroscopy of a sample of 48 z~1-3 Spitzer-selected ULIRGs with IRS mid-IR spectra. Our goals are to compute their bolometric emission, and to determine both the presence and relative strength of their AGN and starburst components. We find that strong-PAH sources tend to have higher 160um and 1.2mm fluxes than weak-PAH sources. The depth of the 9.7um silicate feature does not affect MAMBO detectability. We fit the far-IR SEDs of our sample and find an average <L_{IR}>~7x10^{12}Lsun for our z>1.5 sources. Spectral decomposition suggests that strong-PAH sources typically have ~20-30% AGN fractions. Weak-PAH sources by contrast tend to have >~70% AGN fractions, with a few sources having comparable contributions of AGN and starbursts. The optical line diagnostics support the presence of AGN in the bulk of the weak-PAH sources. With one exception, our sources are narrow-line sources, show no obvious correspondence between the optical extinction and the silicate feature depth, and, in two cases, show evidence for outflows. Radio AGN are present in both strong-PAH and weak-PAH sources. This is supported by our sample's far-IR-to-radio ratios (q) being consistently below the average value of 2.34 for local star-forming galaxies. We use survival analysis to include the lower-limits given by the radio-undetected sources, arriving at <q>=2.07+/-0.01 for our z>1.5 sample. In total, radio and, where available, optical line diagnostics support the presence of AGN in 57% of the z>1.5 sources, independent of IR-based diagnostics. For higher-z sources, the AGN luminosities alone are estimated to be >10^{12}Lsun, which, supported by the [OIII] luminosities, implies that the bulk of our sources host obscured quasars.

We present CO(1-0) and CO(2-1) maps of the starburst/Seyfert 1 galaxy NGC 2782, obtained with the IRAM interferometer. The CO emission is aligned along the stellar nuclear bar of radius 1 kpc, configured in an elongated structure with two spiral arms at high pitch angle. At the extremity of the nuclear bar, the CO changes direction to trace two more extended spiral features at a lower pitch angle. These are the beginning of two straight dust lanes, which are aligned parallel to an oval distortion, reminiscent of a primary bar, almost perpendicular to the nuclear one. The two embedded bars appear in Spitzer IRAC near-infrared images, and HST color images, although highly obscured by dust in the latter. We compute the torques exerted by the stellar bars on the gas, and find systematically negative average torques down to the resolution limit of the images, providing evidence of gas inflow tantalizingly close to the nucleus of NGC 2782. The observations are well reproduced by numerical simulations, including gas dissipation, which predict the secondary bar decoupling, the formation of an elongated ring at the 1 kpc-radius Inner Lindblad Resonance (ILR) of the primary bar, and the gas inflow to the ILR of the nuclear bar. The presence of molecular gas inside the ILR of the primary bar, transported by a second nuclear bar, is a potential ``smoking gun''; the gas there is certainly fueling the central starburst, and in a second step could fuel directly the AGN.

We present a simple set of kinematic criteria that can distinguish between galaxies dominated by ordered rotational motion and those involved in major merger events. Our criteria are based on the dynamics of the warm ionized gas (as traced by H-alpha) within galaxies, making this analysis accessible to high-redshift systems, whose kinematics are primarily traceable through emission features. Using the method of kinemetry (developed by Krajnovic and co-workers), we quantify asymmetries in both the velocity and velocity dispersion maps of the warm gas, and the resulting criteria enable us to empirically differentiate between non-merging and merging systems at high redshift. We apply these criteria to 11 of our best-studied rest-frame UV/optical-selected z~2 galaxies for which we have near infrared integral field spectroscopic data from SINFONI on the VLT. Of these 11 systems, we find that >50% have kinematics consistent with a single rotating disk interpretation, while the remaining systems are more likely undergoing major mergers. This result, combined with the short formation timescales of these systems, provides evidence that rapid, smooth accretion of gas plays a significant role in galaxy formation at high redshift.

This paper addresses the origin of the silicate emission observed in PG QSOs, based on observations with the Spitzer Space Telescope. Scenarios based on the unified model suggest that silicate emission in AGN arises mainly from the illuminated faces of the clouds in the torus at temperatures near sublimation. However, detections of silicate emission in Type 2 QSOs, and the estimated cool dust temperatures, argue for a more extended emission region.To investigate this issue we present the mid-infrared spectra of 23 QSOs. These spectra, and especially the silicate emission features at ~10 and ~18 mu can be fitted using dusty narrow line region (NLR) models and a combination of black bodies. The bolometric luminosities of the QSOs allow us to derive the radial distances and covering factors for the silicate-emitting dust. The inferred radii are 100-200 times larger than the dust sublimation radius, much larger than the expected dimensions of the inner torus. Our QSO mid-IR spectra are consistent with the bulk of the silicate dust emission arising from the dust in the innermost parts of the NLR.

We report sub-arcsecond resolution IRAM PdBI millimeter CO interferometry of four z~2 submillimeter galaxies (SMGs), and sensitive CO (3-2) flux limits toward three z~2 UV-/optically selected star forming galaxies. The new data reveal for the first time spatially resolved CO gas kinematics in the observed SMGs. Two of the SMGs show double or multiple morphologies, with complex, disturbed gas motions. The other two SMGs exhibit CO velocity gradients of ~500 km/s across 0.2 arcsec (1.6 kpc) diameter regions, suggesting that the star forming gas is in compact, rotating disks. Our data provide compelling evidence that these SMGs represent extreme, short-lived 'maximum' star forming events in highly dissipative mergers of gas rich galaxies. The resulting high mass surface and volume densities of SMGs are similar to those of compact quiescent galaxies in the same redshift range, and much higher than those in local spheroids. From the ratio of the comoving volume densities of SMGs and quiescent galaxies in the same mass and redshift ranges, and from the comparison of gas exhaustion time scales and stellar ages, we estimate that the SMG phase duration is about 100 Myrs. Our analysis of SMGs and optically/UV selected high redshift star forming galaxies supports a 'universal' Chabrier IMF as being valid over the star forming history of these galaxies. We find that the 12CO luminosity to total gas mass conversion factors at z~2-3 are probably similar to those assumed at z~0. The implied gas fractions in our sample galaxies range from 20 to 50%.