We present the detection of multiple CO line transitions with ALMA in a few tens of infrared-selected galaxies on and above the main sequence at z=1.1-1.7. We reliably detected the emission of CO(5-4), CO(2-1), and CO(7-6)+[CI](2-1) in 50, 33, and 13 galaxies, respectively, and we complemented this information with available CO(4-3) and [CI](1-0) fluxes for part of the sample, and modeling of the optical-to-mm SEDs. We retrieve a quasi-linear relation between LIR and CO(5-4) or CO(7-6) for main-sequence galaxies and starbursts, corroborating the hypothesis that these transitions can be used as SFR tracers. We find the CO excitation to steadily increase as a function of the star formation efficiency, the mean intensity of the radiation field warming the dust, the surface density of SFR, and, less distinctly, with the distance from the main sequence. This adds to the tentative evidence for higher excitation of the CO+[CI] SLED of starbursts relative to that for main-sequence objects, where the dust opacities play a minor role in shaping the high-J CO transitions in our sample. However, the distinction between the average SLED of upper main-sequence and starburst galaxies is blurred, driven by a wide variety of intrinsic shapes. LVG radiative transfer modeling demonstrates the existence of a highly excited component that elevates the CO SLED of high-redshift main-sequence and starbursting galaxies above the typical values observed in the disk of the Milky Way. This excited component is dense and it encloses ~50% of the total molecular gas mass in main-sequence objects. We interpret the observed trends involving the CO excitation as mainly driven by a combination of large SFRs and compact sizes, naturally connected with enhanced dense molecular gas fractions and higher dust and gas temperatures, due to increasing UV radiation fields, cosmic ray rates, and dust/gas coupling. [Abridged]

We have discovered a 300kpc-wide giant Lya nebula centered on the massive galaxy group RO-1001 at z=2.91 in the COSMOS field. Keck Cosmic Web Imager observations reveal three cold gas filaments converging into the center of the potential well of its ~4x10^{13}Msun dark matter halo, hosting 1200Msun/yr of star formation as probed by ALMA and NOEMA observations. The inwards deceleration of the filament and the prevalence of blueshifted components in the Lya spectra support a scenario of gas accretion. The upper limits on AGN activity and overall energetics favor gravity as the primary Lya powering source and infall as the main source of gas flows to the system. Although interpretational difficulties remain, with outflows and likely also recombinations still playing a role, this finding provides arguably an ideal environment to quantitatively test models of cold gas accretion and galaxy feeding inside a massive halo at high redshift.

We present a characterization of the physical properties of a sample of 35 securely-detected, dusty galaxies in the deep ALMA 1.2-mm image obtained as part of the ALMA Spectroscopic Survey in the {\it Hubble} Ultra Deep Field (ASPECS) Large Program. This sample is complemented by 26 additional sources identified via an optical/infrared source positional prior. Using their well-characterized spectral energy distributions, we derive median stellar masses and star formation rates (SFR) of $4.8\times10^{10}~M_\odot$ and 30 $M_\odot$ yr$^{-1}$, and interquartile ranges of $(2.4-11.7)\times10^{10}~M_\odot$ and $20-50~M_\odot$ yr$^{-1}$, respectively. We derive a median spectroscopic redshift of 1.8 with an interquartile range $1.1-2.6$, significantly lower than sub-millimeter galaxies detected in shallower, wide-field surveys. We find that 59\%$\pm$13\%, 6\%$\pm$4\%, and 34\%$\pm$9\% of our sources are within, above and below $\pm0.4$ dex from the SFR- stellar mass relation or main-sequence (MS), respectively. The ASPECS galaxies closely follow the SFR-molecular gas mass relation and other previously established scaling relations, confirming a factor of five increase of the gas-to-stellar mass ratio from $z=0.5$ to $z=2.5$, and a mild evolution of the gas depletion timescales with a typical value of 0.7 Gyr at $z=1-3$. ASPECS galaxies located significantly below the MS, a poorly exploited parameter space, have low gas-to-stellar mass ratios of $\sim0.1-0.2$ and long depletion timescales $>1$ Gyr. Galaxies along the MS dominate the cosmic density of molecular gas at all redshifts. Systems above the main sequence have an increasing contribution to the total gas reservoirs from $z<1$ to $z=2.5$, while the opposite is found for galaxies below the MS.

We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2.

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues.

We report the detection of CO($J$=2$\to$1) emission from three massive dusty starburst galaxies at $z$$>$5 through molecular line scans in the NSF's Karl G. Jansky Very Large Array (VLA) CO Luminosity Density at High Redshift (COLDz) survey. Redshifts for two of the sources, HDF 850.1 ($z$=5.183) and AzTEC-3 ($z$=5.298), were previously known. We revise a previous redshift estimate for the third source GN10 ($z$=5.303), which we have independently confirmed through detections of CO $J$=1$\to$0, 5$\to$4, 6$\to$5, and [CII] 158 $\mu$m emission with the VLA and the NOrthern Extended Milllimeter Array (NOEMA). We find that two currently independently confirmed CO sources in COLDz are "optically dark", and that three of them are dust-obscured galaxies at $z$$>$5. Given our survey area of $\sim$60 arcmin$^2$, our results appear to imply a $\sim$6-55 times higher space density of such distant dusty systems within the first billion years after the Big Bang than previously thought. At least two of these $z$$>$5 galaxies show star-formation rate surface densities consistent with so-called "maximum" starbursts, but we find significant differences in CO excitation between them. This result may suggest that different fractions of the massive gas reservoirs are located in the dense, star-forming nuclear regions - consistent with the more extended sizes of the [CII] emission compared to the dust continuum and higher [CII]-to-far-infrared luminosity ratios in those galaxies with lower gas excitation. We thus find substantial variations in the conditions for star formation between $z$$>$5 dusty starbursts, which typically have dust temperatures $\sim$57%$\pm$25% warmer than starbursts at $z$=2-3 due to their enhanced star formation activity.

We have obtained spectroscopic confirmation with Hubble Space Telescope WFC3/G141 of a first sizeable sample of nine quiescent galaxies at 2.4<z<3.3. Their average near-UV/optical rest-frame spectrum is characterized by low attenuation (Av$\sim$0.6 mag) and a strong Balmer break, larger than the 4000 A break, corresponding to a fairly young age of $\sim$300 Myr. This formally classifies a substantial fraction of classically selected quiescent galaxies at z $\sim$ 3 as post-starbursts, marking their convergence to the quenching epoch. The rapid spectral evolution with respect to z $\sim$ 1.5 quiescent galaxies is not matched by an increase of residual star-formation, as judged from the weak detection of [O II]3727A emission, pointing to a flattening of the steep increase in gas fractions previously seen from z $\sim$ 0 to 1.8. However, radio 3GHz stacked emission implies either much stronger dust-obscured star formation or substantial further evolution in radio-mode AGN activity with respect to z $\sim$ 1.5.

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$. At all redshifts, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($\lambda_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $\lambda_{\rm EDD}$ distribution, and with a positive shift of the break $\lambda^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($\lambda_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $\lambda_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rm\mu Jy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rm\mu Jy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_\nu \sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

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

[Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value.

We report ALMA observations of the neutral atomic carbon transitions [CI] and multiple CO lines in a sample of $\sim30$ main sequence galaxies at $z\sim1$, including novel information on [CI](2-1) and CO(7-6) for 7 of such normal objects. We complement our observations with a collection of $>200$ galaxies with coverage of similar transitions, spanning the $z=0-4$ redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [CI]/IR and [CI]/high-$J$ ($J=7$) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main sequence galaxies. We further report constant $L'_{\rm [CI]2-1}$/$L'_{\rm [CI]1-0}$ ratios across the galaxy populations and redshifts, suggesting that gas temperatures $T_{\rm exc}$ traced by [CI] do not strongly vary. We find only a mild correlation with $T_{\rm dust}$ and that, generally, $T_{\rm exc} \lesssim T_{\rm dust}$. We fit the line ratios with classical PDR models, retrieving consistently larger densities and intensities of the UV radiation fields in submm galaxies than in main sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submm emission. While we confirm the success in reproducing the CO lines, we find systematically larger [CI] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [CI] emission on galactic scales. We release our data compilation to the community.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

We present hitherto the largest sample of gas-phase metallicity radial gradients measured at sub-kiloparsec resolution in star-forming galaxies in the redshift range of $z\in[1.2, 2.3]$. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 79 galaxies with stellar mass ranging from 10$^7$ to 10$^{10}$ $M_\odot$, instantaneous star-formation rate in the range of [1, 100] $M_\odot$/yr, and global metallicity [$\frac{1}{12}$, 2] solar. At 2-$\sigma$ confidence level, 14/79 galaxies in our sample show negative radial gradients, whereas 7/79 show inverted gradients. Combining ours and all other metallicity gradients obtained at similar resolution currently available in the literature, we measure a negative mass dependence of $\Delta\log({\rm O/H})/\Delta r~ [\mathrm{dex~kpc^{-1}}] = (-0.020\pm0.007) + (-0.014\pm0.008) \log(M_\ast/10^{9.4}M_\odot)$ with the intrinsic scatter being $\sigma=0.060\pm0.006$ over four orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemo-structural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star-formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of $M_\ast\lesssim10^9 M_\odot$.

We report Atacama Large Millimetre Array (ALMA) observations of four high-redshift dusty star-forming galaxy candidates selected from far-Infrared (FIR)/submm observations in the COSMOS field. We securely detect all galaxies in the continuum and spectroscopically confirm them at z=3.62--5.85 using ALMA 3mm line scans, detecting multiple CO and/or [CI] transitions. This includes the most distant dusty galaxy currently known in the COSMOS field, ID85001929 at z=5.847. These redshifts are lower than we had expected as these galaxies have substantially colder dust temperatures (i.e., their SEDs peak at longer rest frame wavelengths) than most literature sources at z>4. The observed cold dust temperatures are best understood as evidence for optically thick dust continuum in the FIR, rather than the result of low star formation efficiency with rapid metal enrichment. We provide direct evidence that, given their cold spectral energy distributions, CMB plays a significant role biasing their observed Rayleigh-Jeans (RJ) slopes to unlikely steep values and, possibly, reducing their CO fluxes by a factor of two. We recover standard RJ slopes when the CMB contribution is taken into account. High resolution ALMA imaging shows compact morphology and evidence for mergers. This work reveals a population of cold dusty star-forming galaxies that were under-represented in current surveys, and are even colder than typical Main Sequence galaxies at the same redshift. High FIR dust optical depth might be a widespread feature of compact starbursts at any redshift.

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

We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time.

We present ALMA observations of the CO(1-0) line and 3-mm continuum emission in eight ultraluminous infrared (IR) quasi-stellar objects (QSOs) at z = 0.06-0.19. All eight IR QSO hosts are clearly resolved in their CO molecular gas emission with a median source size of 3.2 kpc, and seven out of eight sources are detected in 3-mm continuum, which is found to be more centrally concentrated with respect to molecular gas with sizes of 0.4-1.0 kpc. Our observations reveal a diversity of CO morphology and kinematics for the IR QSO systems which can be roughly classified into three categories, rotating gas disk with ordered velocity gradient, compact CO peak with disturbed velocity, and multiple CO distinct sources undergoing a merger between luminous QSO and a companion galaxy separated by a few kpc. The molecular gas in three of IR QSO hosts are found to be rotation-dominated with the ratio of the maximum rotation velocity to the local velocity dispersion of $V_{\rm rot}/\sigma=4-6$. Basic estimates of the dynamical masses within the CO-emitting regions give masses between $7.4\times10^9$ and $6.9\times10^{10}$ $M_\odot$. We find an increasing trend between BH mass accretion rate and star formation rate (SFR) over three orders of magnitude in far-IR luminosity/SFR, in line with the correlation between QSO bolometric luminosity and SF activity, indicative of a likely direct connection between AGN and SF activity over galaxy evolution timescales.

Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $\epsilon \sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$).

While the formation of stellar clumps in distant galaxies is usually attributed to gravitational violent disk instabilities, we show here that major mergers also represent a competitive mechanism to form bright clumps. Using ~0.1'' resolution ACS F814W images in the entire COSMOS field, we measure the fraction of clumpy emission in 109 main sequence (MS) and 79 Herschel-detected starbursts (off-MS) galaxies at 0.5 < z < 0.9, representative of normal versus merger induced star-forming activity, respectively. We additionally identify merger samples from visual inspection and from Gini-M20 morphological parameters. Regardless of the merger criteria adopted, the clumpiness distribution of merging systems is different from that of normal isolated disks at > 99.5 % confidence level, with the former reaching higher clumpiness values, up to 20 % of the total galaxy emission. We confirm the merger induced clumpiness enhancement with novel hydrodynamical simulations of colliding galaxies with gas fractions typical of z ~ 0.7. Multi-wavelength images of three starbursts in the CANDELS field support the young nature of clumps, which are likely merger products rather than older pre-existing structures. Finally, for a subset of 19 starbursts with existing near-IR rest frame spectroscopy, we find that the clumpiness is mildly anti-correlated with the merger phase, decreasing towards final coalescence. Our result can explain recent ALMA detections of clumps in hyperluminous high-z starbursts, while normal objects are smooth. This work raises a question on the role of mergers on the origin of clumps in high redshift galaxies in general.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

We report the first sub-kiloparsec spatial resolution measurements of strongly inverted gas-phase metallicity gradients in two dwarf galaxies at $z$$\sim$2. The galaxies have stellar masses $\sim$$10^9M_\odot$, specific star-formation rate $\sim$20 Gyr$^{-1}$, and global metallicity $12+\log({\rm O/H})\sim8.1$ (1/4 solar), assuming the Maiolino et al. (2008) strong line calibrations of OIII/Hb and OII/Hb. Their metallicity radial gradients are measured to be highly inverted, i.e., 0.122$\pm$0.008 and 0.111$\pm$0.017 dex/kpc, which is hitherto unseen at such small masses in similar redshift ranges. From the Hubble Space Telescope observations of the source nebular emission and stellar continuum, we present the 2-dimensional spatial maps of star-formation rate surface density, stellar population age, and gas fraction, which show that our galaxies are currently undergoing rapid mass assembly via disk inside-out growth. More importantly, using a simple chemical evolution model, we find that the gas fractions for different metallicity regions cannot be explained by pure gas accretion. Our spatially resolved analysis based on a more advanced gas regulator model results in a spatial map of net gaseous outflows, triggered by active central starbursts, that potentially play a significant role in shaping the spatial distribution of metallicity by effectively transporting stellar nucleosynthesis yields outwards. The relation between wind mass loading factors and stellar surface densities measured in different regions of our galaxies shows that a single type of wind mechanism, driven by either energy or momentum conservation, cannot explain the entire galaxy. These sources present a unique constraint on the effects of gas flows on the early phase of disk growth from the perspective of spatially resolved chemical evolution within individual systems.

We present high-fidelity, 30 milliarcsecond (200-pc) resolution ALMA rest-frame 240 $\mu$m observations of cold dust emission in three typical main-sequence star-forming galaxies (SFGs) at $z \sim 3$ in the Hubble Ultra-Deep Field (HUDF). The cold dust is distributed within the smooth disk-like central regions of star formation $1 - 3$ kpc in diameter, despite their complex and disturbed rest-frame UV and optical morphologies. No dust substructures or clumps are seen down to $\simeq 1- 3$ $M_\odot$yr$^{-1}$ (1$\sigma$) per 200-pc beam. No dust emission is observed at the locations of UV-emitting clumps, which lie $\simeq 2-10$ kpc from the bulk of star formation. Clumpy substructures can contribute no more than $1-7$% of the total star formation in these galaxies (3$\sigma$ upper limits). The lack of star-forming substructures in our HUDF galaxies is to be contrasted with the multiple substructures characteristic of submillimeter-selected galaxies (SMGs) at the same cosmic epoch, particularly the far-IR-bright SMGs with similarly high-fidelity ALMA observations of Hodge et al. (2019). Individual star-forming substructures in these SMGs contain $\sim10-30$% of their total star formation. A substructure in these SMGs is often comparably bright in the far-infrared as (or in some cases brighter than) our typical SFGs, suggesting that these SMGs originate from a class of disruptive event involving multiple objects at the scale of our HUDF galaxies. The scale of the disruptive event found in our main-sequence SFGs, characterized by the lack of star-forming substructures at our resolution and sensitivity, could be less violent, e.g., gas-rich disk instability or minor mergers.

We present an enhanced version of the multiwavelength spectral modeling code MAGPHYS that allows the estimation of galaxy photometric redshift and physical properties (e.g., stellar mass, star formation rate, dust attenuation) simultaneously, together with robust characterization of their uncertainties. The self-consistent modeling over ultraviolet to radio wavelengths in MAGPHYS+photo-z is unique compared to standard photometric redshift codes. The broader wavelength consideration is particularly useful for breaking certain degeneracies in color vs. redshift for dusty galaxies with limited observer-frame ultraviolet and optical data (or upper limits). We demonstrate the success of the code in estimating redshifts and physical properties for over 4,000 infrared-detected galaxies at 0.4<z<6.0 in the COSMOS field with robust spectroscopic redshifts. We achieve high photo-z precision ($\sigma_{\Delta z/(1+z_{spec})}\lesssim0.04$), high accuracy (i.e., minimal offset biases; median$(\Delta z/(1+z_{spec}))\lesssim0.02$), and low catastrophic failure rates ($\eta\simeq4\%$) over all redshifts. Interestingly, we find that a weak 2175A absorption feature in the attenuation curve models is required to remove a subtle systematic photo-z offset ($z_{phot}-z_{spec}\simeq-0.03$) that occurs when this feature is not included. As expected, the accuracy of derived physical properties in MAGPHYS+photo-z decreases strongly as redshift uncertainty increases. The all-in-one treatment of uncertainties afforded with this code is beneficial for accurately interpreting physical properties of galaxies in large photometric datasets. Finally, we emphasize that MAGPHYS+photo-z is not intended to replace existing photo-z codes, but rather offer flexibility to robustly interpret physical properties when spectroscopic redshifts are unavailable. The MAGPHYS+photo-z code is publicly available online.

We investigate the contribution of clumps and satellites to the galaxy mass assembly. We analyzed spatially-resolved Hubble Space Telescope observations (imaging and slitless spectroscopy) of 53 star-forming galaxies at z ~ 1 - 3. We created continuum and emission line maps and pinpointed residual "blobs" detected after subtracting the galaxy disk. Those were separated into compact (unresolved) and extended (resolved) components. Extended components have sizes ~ 2 kpc and comparable stellar mass and age as the galaxy disks, whereas the compact components are 1.5 dex less massive and 0.4 dex younger than the disks. Furthermore the extended blobs are typically found at larger distances from the galaxy barycenter than the compact ones. Prompted by these observations and by the comparison with simulations, we suggest that compact blobs are in-situ formed clumps, whereas the extended ones are accreting satellites. Clumps and satellites enclose respectively ~ 20% and ~< 80% of the galaxy stellar mass, ~ 30% and ~ 20% of its star formation rate. Considering the compact blobs, we statistically estimated that massive clumps (Mstar >~ 10^9 Msun) have lifetimes of ~ 650 Myr, and the less massive ones (10^8 < Mstar < 10^9 Msun) of ~ 145 Myr. This supports simulations predicting long-lived clumps (lifetime > 100 Myr). Finally, ~< 30% (13%) of our sample galaxies are undergoing single (multiple) merger(s), they have a projected separation ~< 10 kpc, and the typical mass ratio of our satellites is 1:5 (but ranges between 1:10 and 1:1), in agreement with literature results for close pair galaxies.

We present Atacama Large Millimetre Array and Atacama Compact Array observations of the Sunyaev-Zel'dovich effect in the z = 2 galaxy cluster Cl J1449+0856, an X-ray-detected progenitor of typical massive clusters in the present day Universe. While in a cleaned but otherwise untouched 92 GHz map of this cluster, little to no negative signal is visible, careful subtraction of known sub-millimetre emitters in the uv plane reveals a decrement at 5$\sigma$ significance. The total signal is -190$\pm$36 $\mu$Jy, with a peak offset by 5"-9" ($\sim$50 kpc) from both the X-ray centroid and the still-forming brightest cluster galaxy. A comparison of the recovered uv-amplitude profile of the decrement with different pressure models allows us to derive total mass constraints consistent with the $\sim$6$\times$10$^{13}$ M$_{\odot}$ estimated from X-ray data. Moreover, we find no strong evidence for a deviation of the pressure profile with respect to local galaxy clusters, although a slight tension at small-to-intermediate spatial scales suggests a flattened central profile, opposite to what seen in a cool core and possibly an AGN-related effect. This analysis of the lowest mass single SZ detection so far illustrates the importance of interferometers when observing the SZ effect in high-redshift clusters, the cores of which cannot be considered quiescent, such that careful subtraction of galaxy emission is necessary.

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

We used the Plateau De Bure Interferometer to observe multiple CO and neutral carbon transitions in a z=2.2 main sequence disk galaxy, BX610. Our observation of CO(7-6), CO(4-3), and both far-infrared(FIR) [CI] lines complements previous observations of H$\alpha$ and low-J CO, and reveals a galaxy that is vigorously forming stars with UV fields (Log($G$ G$_0^{-1}) \lesssim3.25);$ although less vigorously than local ultra-luminous infrared galaxies or most starbursting submillimeter galaxies in the early universe. Our observations allow new independent estimates of the cold gas mass which indicate $M_\textrm{gas}\sim2\times10^{11}$M$_\odot$, and suggest a modestly larger $\alpha_{\textrm{CO}}$ value of $\sim$8.2. The corresponding gas depletion timescale is $\sim$1.5 Gyr. In addition to gas of modest density (Log($n$ cm$^3)\lesssim3$ ) heated by star formation, BX610 shows evidence for a significant second gas component responsible for the strong high-J CO emission. This second component might either be a high-density molecular gas component heated by star formation in a typical photodissociation region, or could be molecular gas excited by low-velocity C shocks. The CO(7-6)-to-FIR luminosity ratio we observe is significantly higher than typical star-forming galaxies and suggests that CO(7-6) is not a reliable star-formation tracer in this galaxy.

We investigate the fraction of close pairs and morphologically identified mergers on and above the star-forming main sequence (MS) at 0.2$\leq z\leq$2.0. The novelty of our work lies in the use of a non-parametric morphological classification performed on resolved stellar mass maps, reducing the contamination by non-interacting, high-redshift clumpy galaxies. We find that the merger fraction rapidly rises to $\geq$70% above the MS, implying that -- already at $z{\gtrsim}1$ -- starburst (SB) events ($\Delta_{\rm MS}\geq$0.6) are almost always associated with a major merger (1:1 to 1:6 mass ratio). The majority of interacting galaxies in the SB region are morphologically disturbed, late-stage mergers. Pair fractions show little dependence on MS-offset and pairs are more prevalent than late-stage mergers only in the lower half of the MS. In our sample, major mergers on the MS occur with a roughly equal frequency of $\sim$5-10% at all masses ${\gtrsim} 10^{10}M_{\odot}$. The MS major merger fraction roughly doubles between $z=0.2$ and $z=2$, with morphological mergers driving the overall increase at $z{\gtrsim}1$. The differential redshift evolution of interacting pairs and morphologically classified mergers on the MS can be reconciled by evolving observability timescales for both pairs and morphological disturbances. The observed variation of the late-stage merger fraction with $\Delta_{\rm MS}$ follows the perturbative 2-Star Formation Mode model, where any MS galaxy can experience a continuum of different SFR enhancements. This points to a starburst-merger connection not only for extreme events, but also more moderate bursts which merely scatter galaxies upward within the MS, rather than fully elevating them above it.

We present ALMA Band~7 observations at $850\mu$m of 20 luminous ($\log\, L_{\rm bol}>46.9$ [erg s$^{-1}$]) unobscured quasars at $z\sim2$. We detect continuum emission for 19/20 quasars. After subtracting an AGN contribution, we measure the total far-IR luminosity for 18 quasars, assuming a modified blackbody model, and attribute the emission as indicative of the star formation rate (SFR). Our sample can be characterized with a log-normal SFR distribution having a mean of 140 $M_\odot$ yr$^{-1}$ and a dispersion of 0.5 dex. Based on an inference of their stellar masses, the SFRs are similar, in both the mean and dispersion, with star-forming main-sequence galaxies at the equivalent epoch. Thus, there is no evidence for a systematic enhancement or suppression (i.e., regulation or quenching) of star formation in the hosts of the most luminous quasars at $z\sim2$. These results are consistent with the Magneticum cosmological simulation, while in disagreement with a widely recognized phenomenological model that predicts higher SFRs than observed here based on the high bolometric luminosities of this sample. Furthermore, there is only a weak relation between SFR and accretion rate onto their supermassive black holes both for average and individual measurements. We interpret these results as indicative of star formation and quasar accretion being fed from the available gas reservoir(s) in their host with a disconnect due to their different physical sizes, temporal scales, and means of gas processing.

We investigate a sample of 152 dusty sources at 1:5 < z < 2:5 to understand the connection of enhanced Star-Formation-Rate (SFR) and Black-Hole-Accretion-Rate (BHAR). The sources are Herschel-selected, having stellar masses M*>10^10Msun and SFR (100-1000Msun/yr) elevated(>4?) above the star-forming "main sequence", classifying them as Starbursts (SB). Through a multiwavelength fitting approach (including a dusty torus component), we divided the sample into active SBs (dominated by an AGN emission, SBs-AGN, ? 23% of the sample) and purely star-forming SBs (SBs-SFR). We visually inspected their HST/UV-restframe maps: SBs-SFR are generally irregular and composite systems; ? 50% of SBs-AGN are instead dominated by regular compact morphologies. We then found archival ALMA continuum counterparts for 33 galaxies (12 SBs-AGN and 21 SBs-SFR). For these sources we computed dust masses, and, with standard assumptions, we also guessed total molecular gas-masses. SBs turn to be gas rich systems (fgas = Mgas=Mgas/(Mgas+M*) \sim 20%-70%), and the gas fractions of the two SB classes are very similar (fgas = 43 +/-4% and fgas = 42+/- 2%). Our results show that SBs are consistent with a mixture of: 1) highly star-forming merging systems (dominating the SBs-SFR), and 2) primordial galaxies, rapidly growing their M* together with their Black Hole (mainly the more compact SBs-AGN). Anyway, feedback effects have not reduced their fgas yet. Indeed, SBs at z = 2, with relatively low bolometric AGN luminosities in the range 10^44 < Lbol(AGN) < 10^46 erg/s (compared to bright optical and X-ray quasars), are still relatively far from the epoch when the AGN feedback will quench the SFR in the host and will substantially depress the gas fractions.

ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having $10^{8} \leqslant M_{\star} \leqslant 10^{11.5} \ M_{\odot}$ and $20 \leqslant SFR \leqslant 680 \ M_{\odot} \rm yr^{-1}$. The 32 size measurements and 61 upper limits are measured on ALMA images that combine observations of CO(5-4), CO(4-3), CO(2-1) and $\lambda_{\rm obs} \sim 1.1-1.3 \ \rm mm$ continuum, all tracing the star-forming molecular gas. These compact galaxies have instead normally extended $K_{band}$ sizes, suggesting strong specific $SFR$ gradients. Compact galaxies comprise the $50\pm18 \%$ of MS galaxies at $M_{\star} > 10^{11} M_{\odot}$. This is not expected in standard bi-modal scenarios where MS galaxies are mostly steadily-growing extended disks. We suggest that compact MS objects are early post-starburst galaxies in which the merger-driven boost of SF has subsided. They retain their compact SF size until either further gas accretion restores pre-merger galaxy-wide SF, or until becoming quenched. The fraction of merger-affected SF inside the MS seems thus larger than anticipated and might reach $\sim 50$\% at the highest $M_{\star}$. The presence of large galaxies above the MS demonstrates an overall poor correlation between galaxy SF size and specific $SFR$.

We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1\lambda_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) provides new constraints for galaxy formation models on the molecular gas properties of galaxies. We compare results from ASPECS to predictions from two cosmological galaxy formation models: the IllustrisTNG hydrodynamical simulations and the Santa Cruz semi-analytic model (SC SAM). We explore several recipes to model the H$_2$ content of galaxies, finding them to be consistent with one another, and take into account the sensitivity limits and survey area of ASPECS. For a canonical CO-to-H$_2$ conversion factor of $\alpha_{\rm CO} = 3.6\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$ the results of our work include: (1) the H$_2$ mass of $z>1$ galaxies predicted by the models as a function of their stellar mass is a factor of 2-3 lower than observed; (2) the models do not reproduce the number of H$_2$-rich ($M_{\rm H2} > 3\times 10^{10}\,\rm{M}_\odot$) galaxies observed by ASPECS; (3) the H$_2$ cosmic density evolution predicted by IllustrisTNG (the SC SAM) is in tension (only just agrees) with the observed cosmic density, even after accounting for the ASPECS selection function and field-to-field variance effects. The tension between models and observations at $z>1$ can be alleviated by adopting a CO-to-H$_2$ conversion factor in the range $\alpha_{\rm CO} = 2.0 - 0.8\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$. Additional work on constraining the CO-to-H$_2$ conversion factor and CO excitation conditions of galaxies through observations and theory will be necessary to more robustly test the success of galaxy formation models.

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $\beta=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

We report NOrthern Extended Millimetre Array (NOEMA) observations of warm molecular gas traced by CO($5-4$) in a $z \sim 3.2$ gas-rich main-sequence galaxy (MS), initially serendipitously detected in CO($3-2$) emission in `blind' deep NOEMA observations. Our target shows a gas excitation consistent with that seen in $z \sim 1.5$ MS galaxies ($L'_{\rm CO( 5 - 4)}/L'_{\rm CO (3 - 2)} = 0.41 \pm 0.14$), albeit toward the low end, as well as a similar star formation efficiency based on the CO($3-2$) line luminosity and the $L_{\rm IR}$. However, it shows a high molecular gas fraction ($f_{\rm gas} = 0.9\pm 0.2$) as compared to $z\sim 1.5$ MS galaxies ($f_{\rm gas} \sim 0.4$), consistent with a cosmologically increasing gas fraction beyond $z\gtrsim3$ and our current understanding of scaling relations between $z$, $f_{\rm gas}$, the stellar mass $M_*$, and the specific star formation rate sSFR. Our results are consistent with recent findings by the COLDz and ASPECS molecular line scan surveys and suggest that deep searches for CO emission are a powerful means to identify gas-rich, star-forming galaxies at high redshift.

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers.

The physical origin of enhanced star formation activity in interacting galaxies remains an open question. Knowing whether starbursts are triggered by an increase of the quantity of dense gas or an increase of the star formation efficiency would improve our understanding of galaxy evolution and allow to transpose the results obtained in the local Universe to high redshift galaxies. In this paper, we analyze a parsec-resolution simulation of an Antennae-like model of interacting galaxies. We find that the interplay of physical processes has complex and important variations in time and space, through different combinations of mechanisms like tides, shear and turbulence. These can have similar imprints on observables like depletion time and CO emission. The densest gas only constitutes the tail of the density distribution of some clouds, but exists in large excess in others. The super-linearity of the star formation rate dependence on gas density implies that this excess translates into a reduction of depletion times, and thus leads to a deviation from the classical star formation regime, visible up to galactic scales. These clouds are found in all parts of the galaxies, but their number density varies from one region to the next, due to different cloud assembly mechanisms. Therefore, the dependence of cloud and star formation-related quantities (like CO flux and depletion time) on the scale at which they are measured also varies across the galaxies. We find that the $\alpha_{\rm CO}$ conversion factor between the CO luminosity and molecular gas mass has even stronger spatial than temporal variations in a system like the Antennae. These results raise a number of cautionary notes for the interpretation of observations of unresolved star-forming regions, but also predict that the diversity of environments for star formation will be better captured by the future generations of instruments.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

We study a population of significantly sub-solar enrichment galaxies at z=1.99, to investigate how molecular gas, dust and star-formation relate in low-metallicity galaxies at the peak epoch of star-formation. We target our sample with several deep ALMA and VLA datasets, and find no individual detections of CO[4-3], CO[1-0] or dust, in stark contrast to the $>$60% detection rate expected for solar-enrichment galaxies with these MS Halpha SFRs. We find that both low and high density molecular gas (traced by CO[1-0] and CO[4-3] respectively) are affected by the low enrichment, showing sample average (stacked) luminosity deficits $>$0.5-0.7 dex below expectations. This is particularly pertinent for the use of high-J CO emission as a proxy of instantaneous star-formation rate. Our individual galaxy data and stacked constraints point to a strong inverse dependence $\propto Z^{\gamma}$ of gas-to-dust ratios (G/D) and CO-to-H$_{2}$ conversion factors (aco) on metallicity at z$\sim$2, with $\gamma_{\rm G/D}<$-2.2 and $\gamma_{\alpha_{\rm CO}}<$-0.8, respectively. We quantify the importance of comparing G/D and aco vs. metallicity trends from the literature on a common, suitably normalised metallicity scale. When accounting for systematic offsets between different metallicity scales, our z$\sim$2 constraints on these scaling relations are consistent with the corresponding relations for local galaxies. However, among those local relations, we favour those with a steep/double power-law dependence of G/D on metallicity. Finally, we discuss the implications of these findings for (a) gas mass measurements for sub-M$^{*}$ galaxies, and (b) efforts to identify the characteristic galaxy mass scale contributing most to the comoving molecular gas density at z=2.

Based on optical/near-IR Magellan FIRE spectra of 25 starburst galaxies at 0.5 < z < 0.9, Calabr\`o et al.(2018) showed that their attenuation properties can be explained by a single-parameter sequence of total obscurations ranging from A(V)=2 to A(V)=30 towards the starburst core centers in a mixed stars and dust configuration. We investigate here the origin of this sequence for the same sample. We show that total attenuations anti-correlate with the starburst sizes in radio (3 GHz) with a significance larger than 5sigma and a scatter of 0.26 dex. More obscured and compact starbursts also show enhanced N2 (=[NII]/Halpha) ratios and larger line velocity widths that we attribute to an increasing shock contribution toward later merger phases, driven by deeper gravitational potential wells at the coalescence. Additionally, the attenuation is also linked to the equivalent width (EW) of hydrogen recombination lines, which is sensitive to the luminosity weighted age of the relatively unobscured stellar populations. Overall, the correlations among A(V), radio size, line width, N2 and EW of Balmer/Paschen lines converge towards suggesting an evolutionary sequence of merger stages: all of these quantities are likely to be good time-tracers of the merger phenomenon, and their large spanned range appears to be characteristic of the different merger phases. Half of our sample at higher obscurations have radio sizes approximately 3 times smaller than early type galaxies at the same redshift, suggesting that, in analogy with local Ultraluminous Infrared galaxies (ULIRGs), these cores cannot be directly forming elliptical galaxies. Finally, we detect mid-IR AGN torus for half of our sample and additional X-ray emission for 6 starbursts; intriguingly, the latter have systematically more compact sizes, suggestive of emerging AGNs towards later merger stages, possibly precursors of a later QSO phase.

We report the first detailed measurement of the shape of the CO luminosity function at high redshift, based on $>$320 hr of the NSF's Karl G. Jansky Very Large Array (VLA) observations over an area of $\sim$60 arcmin$^2$ taken as part of the CO Luminosity Density at High Redshift (COLDz) survey. COLDz "blindly" selects galaxies based on their cold gas content through CO($J$=1$\to$0) emission at $z$$\sim$2-3 and CO($J$=2$\to$1) at $z$$\sim$5-7 down to a CO luminosity limit of log($L'_{\rm CO}$/K km s$^{-1}$ pc$^2$)$\simeq$9.5. We find that the characteristic luminosity and bright end of the CO luminosity function are substantially higher than predicted by semi-analytical models, but consistent with empirical estimates based on the infrared luminosity function at $z$$\sim$2. We also present the currently most reliable measurement of the cosmic density of cold gas in galaxies at early epochs, i.e., the cold gas history of the universe, as determined over a large cosmic volume of $\sim$375,000 Mpc$^3$. Our measurements are in agreement with an increase of the cold gas density from $z$$\sim$0 to $z$$\sim$2-3, followed by a possible decline towards $z$$\sim$5-7. These findings are consistent with recent surveys based on higher-$J$ CO line measurements, upon which COLDz improves in terms of statistical uncertainties by probing $\sim$50-100 times larger areas and in the reliability of total gas mass estimates by probing the low-$J$ CO lines accessible to the VLA. Our results thus appear to suggest that the cosmic star-formation rate density follows an increased cold molecular gas content in galaxies towards its peak about 10 billion years ago, and that its decline towards the earliest epochs is likely related to a lower overall amount of cold molecular gas (as traced by CO) bound in galaxies towards the first billion years after the Big Bang.

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

We present a new data release from the Fiber Multi-Object Spectrograph (FMOS)-COSMOS survey, which contains the measurements of spectroscopic redshift and flux of rest-frame optical emission lines (H$\alpha$, [NII], [SII], H$\beta$, [OIII]) for 1931 galaxies out of a total of 5484 objects observed over the 1.7 deg$^2$ COSMOS field. We obtained $H$-band and $J$-band medium-resolution ($R\sim3000$) spectra with FMOS mounted on the Subaru telescope, which offers an in-fiber line flux sensitivity limit of $\sim 1 \times 10^{-17}~\mathrm{erg~s^{-1}~cm^{-2}}$ for an on-source exposure time of five hours. The full sample contains the main population of star-forming galaxies at $z\sim1.6$ over the stellar mass range $10^{9.5}\lesssim M_\ast/M_\odot \lesssim 10^{11.5}$, as well as other subsamples of infrared-luminous galaxies detected by Spitzer and Herschel at the same and lower ($z\sim0.9$) redshifts and X-ray emitting galaxies detected by Chandra. This paper presents an overview of our spectral analyses, a description of the sample characteristics, and a summary of the basic properties of emission-line galaxies. We use the larger sample to re-define the stellar mass--star formation rate relation based on the dust-corrected H$\alpha$ luminosity, and find that the individual galaxies are better fit with a parametrization including a bending feature at $M_\ast\approx10^{10.2}~M_\odot$, and that the intrinsic scatter increases with $M_\ast$ from 0.19 to $0.37$ dex. We also confirm with higher confidence that the massive ($M_\ast\gtrsim10^{10.5}~M_\odot$) galaxies are chemically mature as much as local galaxies with the same stellar masses, and that the massive galaxies have lower [SII]/H$\alpha$ ratios for their [OIII]/H$\beta$, as compared to local galaxies, which is indicative of enhancement in ionization parameter.

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex.

The conversion factor $\alpha_{\rm CO}$ from the observable CO(1-0) luminosity to the mass of molecular gas is known to vary between isolated galaxies and some mergers, but the underlying reasons are not clearly understood. Thus, the value(s) of $\alpha_{\rm CO}$ to be adopted remain highly uncertain. To provide better constraints, we apply the large velocity gradient method to a series of hydrodynamical simulations of galaxies and derive the evolution of $\alpha_{\rm CO}$. We report significant variations of $\alpha_{\rm CO}$, and identify three distinct regimes: disk galaxies, starbursts and post-burst phases. We show that estimating the star formation rate over 20 Myr smooths out some of these differences, but still maintains a distinction between disks and starbursts. We find a tighter correlation of $\alpha_{\rm CO}$ with the gas depletion time than with star formation rate, yet with deviations induced by the transitions to and from the starburst episodes. We conclude that $\alpha_{\rm CO}$ fluctuates because of both the feedback energy and the velocity dispersion. Identifying the phase of an interaction by classical means (e.g. morphology, luminosity) could then help selecting the relevant conversion factor to be used and get more accurate estimates of the molecular masses of galaxies.

We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a log-normal distribution (i.e. the same shape as that of MS galaxies), galaxies hosting more powerful X-ray AGNs (i.e. Lx>2x10^43erg s-1) display a narrower RMS distribution that is shifted to higher values compared to their lower Lx counterparts. In addition, we find that more powerful X-ray AGNs have SFRs that are more consistent with that of MS galaxies compared to lower Lx AGNs. Despite this, the mean SFRs (as opposed to RMS) measured from these distributions are consistent with the previously observed flat relationship between SFR and Lx. Our results suggest that the typical star-forming properties of AGN hosts change with Lx , and that more powerful AGNs typically reside in more MS-like star-forming galaxies compared to lower Lx AGNs.

We present the first results of an ALMA survey of the lower fine structure line of atomic carbon [C I]$(^3P_1\,-\,^{3}P_0)$ in far infrared-selected galaxies on the main sequence at $z\sim1.2$ in the COSMOS field. We compare our sample with a comprehensive compilation of data available in the literature for local and high-redshift starbursting systems and quasars. We show that the [C I]($^3P_1$$\rightarrow$$^3P_0$) luminosity correlates on global scales with the infrared luminosity $L_{\rm IR}$ similarly to low-$J$ CO transitions. We report a systematic variation of $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L_{\rm IR}$ as a function of the galaxy type, with the ratio being larger for main-sequence galaxies than for starbursts and sub-millimeter galaxies at fixed $L_{\rm IR}$. The $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L'_{\rm CO(2-1)}$ and $M_{\rm{[C I]}}$/$M_{\rm dust}$ mass ratios are similar for main-sequence galaxies and for local and high-redshift starbursts within a 0.2 dex intrinsic scatter, suggesting that [C I] is a good tracer of molecular gas mass as CO and dust. We derive a fraction of $f_{\rm{[C\,I]}} = M_{\rm{[C\,I]}} / M_{\rm{C}}\sim3-13$% of the total carbon mass in the atomic neutral phase. Moreover, we estimate the neutral atomic carbon abundance, the fundamental ingredient to calibrate [C I] as a gas tracer, by comparing $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$ and available gas masses from CO lines and dust emission. We find lower [C I] abundances in main-sequence galaxies than in starbursting systems and sub-millimeter galaxies, as a consequence of the canonical $\alpha_{\rm CO}$ and gas-to-dust conversion factors. This argues against the application to different galaxy populations of a universal standard [C I] abundance derived from highly biased samples.

We present a census of the molecular gas properties of galaxies in the most distant known X-ray cluster, CLJ1001, at z=2.51, using deep observations of CO(1-0) with JVLA. In total 14 cluster members with $M_{*} > 10^{10.5} M_{\odot}$ are detected, including all the massive star-forming members within the virial radius, providing the largest galaxy sample in a single cluster at $z > 2$ with CO(1-0) measurements. We find a large variety in the gas content of these cluster galaxies, which is correlated with their relative positions (or accretion states), with those closer to the cluster core being increasingly gas-poor. Moreover, despite their low gas content, the galaxies in the cluster center exhibit an elevated star formation efficiency (SFE=SFR/$M_{\rm gas}$) compared to field galaxies, suggesting that the suppression on the SFR is likely delayed compared to that on the gas content. Their gas depletion time is around $t_{\rm dep} \sim 400$ Myrs, comparable to the cluster dynamical time. This implies that they will likely consume all their gas within a single orbit around the cluster center, and form a passive cluster core by $z\sim2$. This result is one of the first direct pieces of evidence for the influence of environment on the gas reservoirs and SFE of $z > 2$ cluster galaxies, thereby providing new insights into the rapid formation and quenching of the most massive galaxies in the early universe.

We present near-IR spectroscopy in J- and H-band for a large sample of 243 X-ray selected moderate-luminosity type-1 AGN in the COSMOS, SXDS and E-CDF-S survey fields using the multi-object spectrograph Subaru/FMOS. Our sample covers the redshift range 0.5<z<3.0 and an X-ray luminosity range of $10^{43}\lesssim L_X \lesssim 10^{45}$~erg s$^{-1}$. We provide emission-line properties and derived virial black hole mass estimates, bolometric luminosities and Eddington ratios, based on H$\alpha$ (211), H$\beta$ (63) and MgII (4). We compare line widths, luminosities and black hole mass estimates from H$\alpha$ and H$\beta$ and augment these with commensurate measurements of MgII and CIV detected in optical spectra. We demonstrate the robustness of using H$\alpha$, H$\beta$ and MgII as reliable black hole mass estimators for high-z moderate-luminosity AGN, while the use of CIV is prone to large uncertainties (>0.4 dex). We extend a recently proposed correction based on the CIV blueshift to lower luminosities and black hole masses. While our sample shows an improvement in their CIV black hole mass estimates, the deficit of high blueshift sources reduces its overall importance for moderate-luminosity AGN, compared to the most luminous quasars. In addition, we revisit luminosity correlations between $L_{\rm{bol}}$, $L_X$, $L_{\rm{[OIII]}}$, $L_{5100}$ and $L_{\rm{H}\alpha}$ and find them to be consistent with a simple empirical model, based on a small number of well-established scaling relations. We finally highlight our highest redshift AGN, CID 781, at z=4.6 which shows the lowest black hole mass ($\sim10^8$M$_\odot$) among current near-IR samples at this redshift, and is in a state of fast growth.

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

We present an analysis of the molecular gas properties, based on CO(2 - 1) emission, of twelve starburst galaxies at z~1.6 selected by having a boost (>~4x) in their star formation rate (SFR) above the average star-forming galaxy at an equivalent stellar mass. ALMA observations are acquired of six additional galaxies than previously reported through our effort. As a result of the larger statistical sample, we significantly detect, for the first time at high-z, a systematically lower L'_CO/L_IR ratio in galaxies lying above the star-forming `main sequence' (MS). Based on an estimate of alpha_CO (i.e., the ratio of molecular gas mass to L'_CO(1-0)), we convert the observational quantities (e.g., L'_CO/L_IR) to physical units (M_gas/SFR) that represent the gas depletion time or its inverse, the star formation efficiency. We interpret the results as indicative of the star formation efficiency increasing in a continuous fashion from the MS to the starburst regime, whereas the gas fractions remain comparable to those of MS galaxies. Although, the balance between an increase in star-formation efficiency or gas fraction depends on the adopted value of alpha_CO as discussed.

We report on the discovery of a merger-driven starburst at z = 1.52, PACS-787, based on high signal-to-noise ALMA observations. CO(5-4) and continuum emission (850um) at a spatial resolution of 0.3" reveal two compact (r_1/2 ~ 1 kpc) and interacting molecular gas disks at a separation of 8.6 kpc thus indicative of an early stage in a merger. With a SFR of 991 Msun/yr, this starburst event should occur closer to final coalescence, as usually seen in hydrodynamical simulations. From the CO size, inclination, and velocity profile for both disks, the dynamical mass is calculated through a novel method that incorporates a calibration using simulations of galaxy mergers. Based on the dynamical mass, we measure (1) the molecular gas mass, independent from the CO luminosity, (2) the ratio of the total gas mass and the CO(1 - 0) luminosity (alpha_CO = M_gas/L'_CO(1-0)), and (3) the gas-to-dust ratio, with the latter two being lower than typically assumed. We find that the high star formation, triggered in both galaxies, is caused by a set of optimal conditions: a high gas mass/fraction, a short depletion time (t_depl=85 and 67 Myrs) to convert gas into stars, and the interaction of likely counter-rotating molecular disks that may accelerate the loss of angular momentum. The state of interaction is further established by the detection of diffuse CO and continuum emission, tidal debris that bridges the two nuclei and is associated with stellar emission seen by HST/WFC3. This observation demonstrates the power of ALMA to study the dynamics of galaxy mergers at high redshift.

We present a multiwavelength study of an atypical submillimeter galaxy in the GOODS-North field, with the aim to understand its physical properties of stellar and dust emission, as well as the central AGN activity. Although it is shown that the source is likely an extremely dusty galaxy at high redshift, its exact position of submillimeter emission is unknown. With the new NOEMA interferometric imaging, we confirm that the source is a unique dusty galaxy. It has no obvious counterpart in the optical and even NIR images observed with HST at lambda~<1.4um. Photometric-redshift analyses from both stellar and dust SED suggest it to likely be at z~>4, though a lower redshift at z~>3.1 cannot be fully ruled out (at 90% confidence interval). Explaining its unusual optical-to-NIR properties requires an old stellar population (~0.67 Gyr), coexisting with a very dusty ongoing starburst component. The latter is contributing to the FIR emission, with its rest-frame UV and optical light being largely obscured along our line of sight. If the observed fluxes at the rest-frame optical/NIR wavelengths were mainly contributed by old stars, a total stellar mass of ~3.5x10^11Msun would be obtained. An X-ray spectral analysis suggests that this galaxy harbors a heavily obscured AGN with N_H=3.3x10^23 cm^-2 and an intrinsic 2-10 keV luminosity of L_X~2.6x10^44 erg/s, which places this object among distant type 2 quasars. The radio emission of the source is extremely bright, which is an order of magnitude higher than the star-formation-powered emission, making it one of the most distant radio-luminous dusty galaxies. The combined characteristics of the galaxy suggest that the source appears to have been caught in a rare but critical transition stage in the evolution of submillimeter galaxies, where we are witnessing the birth of a young AGN and possibly the earliest stage of its jet formation and feedback.

We present ALMA Band 9 observations of the [C II]158um emission for a sample of 10 main-sequence galaxies at redshift z ~ 2, with typical stellar masses (log M*/Msun ~ 10.0 - 10.9) and star formation rates (~ 35 - 115 Msun/yr). Given the strong and well understood evolution of the interstellar medium from the present to z = 2, we investigate the behaviour of the [C II] emission and empirically identify its primary driver. We detect [C II] from six galaxies (four secure, two tentative) and estimate ensemble averages including non detections. The [C II]-to-infrared luminosity ratio (L[C II]/LIR) of our sample is similar to that of local main-sequence galaxies (~ 2 x 10^-3), and ~ 10 times higher than that of starbursts. The [C II] emission has an average spatial extent of 4 - 7 kpc, consistent with the optical size. Complementing our sample with literature data, we find that the [C II] luminosity correlates with galaxies' molecular gas mass, with a mean absolute deviation of 0.2 dex and without evident systematics: the [C II]-to-H2 conversion factor (alpha_[C II] ~ 30 Msun/Lsun) is largely independent of galaxies' depletion time, metallicity, and redshift. [C II] seems therefore a convenient tracer to estimate galaxies' molecular gas content regardless of their starburst or main-sequence nature, and extending to metal-poor galaxies at low and high redshifts. The dearth of [C II] emission reported for z > 6 - 7 galaxies might suggest either a high star formation efficiency or a small fraction of UV light from star formation reprocessed by dust.

We describe the CO Luminosity Density at High-z (COLDz) survey, the first spectral line deep field targeting CO(1-0) emission from galaxies at $z=1.95-2.85$ and CO(2-1) at $z=4.91-6.70$. The main goal of COLDz is to constrain the cosmic density of molecular gas at the peak epoch of cosmic star formation. By targeting both a wide ($\sim$51 arcmin$^2$) and a deep area ($\sim$9 arcmin$^2$), the survey is designed to robustly constrain the bright end and the characteristic luminosity of the CO(1-0) luminosity function. An extensive analysis of the reliability of our line candidates, and new techniques provide detailed completeness and statistical corrections as necessary to determine the best constraints to date on the CO luminosity function. Our blind search for CO(1-0) uniformly selects starbursts and massive Main Sequence galaxies based on their cold molecular gas masses. Our search also detects CO(2-1) line emission from optically dark, dusty star-forming galaxies at $z>5$. We find a range of spatial sizes for the CO-traced gas reservoirs up to $\sim40$ kpc, suggesting that spatially extended cold molecular gas reservoirs may be common in massive, gas-rich galaxies at $z\sim2$. Through CO line stacking, we constrain the gas mass fraction in previously known typical star-forming galaxies at $z=2$-3. The stacked CO detection suggests lower molecular gas mass fractions than expected for massive Main Sequence galaxies by a factor of $\sim3-6$. We find total CO line brightness at $\sim34\,$GHz of $0.45\pm0.2\,\mu$K, which constrains future line intensity mapping and CMB experiments.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

We present a "super-deblended" far-infrared to (sub)millimeter photometric catalog in the Cosmic Evolution Survey (COSMOS), prepared with the method recently developed by Liu et al. 2018, with key adaptations. We obtain point spread function (PSF) fitting photometry at fixed prior positions including 88,008 galaxies detected in either VLA 1.4~GHz, 3~GHz and/or MIPS 24~$\mu$m images. By adding a specifically carved mass-selected sample (with an evolving stellar mass limit), a highly complete prior sample of 194,428 galaxies is achieved for deblending FIR/(sub)mm images. We performed ``active' removal of non relevant priors at FIR/(sub)mm bands using spectral energy distribution (SED) fitting and redshift information. In order to cope with the shallower COSMOS data we subtract from the maps the flux of faint non-fitted priors and explicitly account for the uncertainty of this step. The resulting photometry (including data from Spitzer, Herschel, SCUBA2, AzTEC, MAMBO and NSF's Karl G. Jansky Very Large Array at 3~GHz and 1.4~GHz) displays well behaved quasi-Gaussian uncertainties, calibrated from Monte Carlo simulations and tailored to observables (crowding, residual maps). Comparison to ALMA photometry for hundreds of sources provide a remarkable validation of the technique. We detect 11,220 galaxies over the 100--1200~$\mu$m range, extending to $z_{\rm phot}\sim7$. We conservatively selected a sample of 85 $z>4$ high redshift candidates, significantly detected in the FIR/(sub)mm, often with secure radio and/or Spitzer/IRAC counterparts. This provides a chance to investigate the first generation of vigorous starburst galaxies (SFRs$\sim1000\mathrm{M}_\odot$~yr$^{-1}$). The photometric and value added catalogs are publicly released.

We obtained optical/near-IR rest-frame Magellan FIRE spectra (including Pa$\beta$ and Pa$\gamma$) of 25 starburst galaxies at 0.5<z<0.9, with average star formation rates (SFR) x7 above the Main Sequence (MS). We find that Paschen-to-Balmer line ratios saturate around a constant value corresponding to $A_{\rm V}\sim$2-3 mag, while line to IR luminosity ratios suggest a large range of more extreme obscurations and appear to be uncorrelated to the former. This behavior is not consistent with standard attenuation laws derived for local and distant galaxies, while being remarkably consistent with observations of starburst cores in which young stars and dust are homogeneously mixed. This model implies $A_{\rm V}=$2-30 mag attenuation to the center of starburst cores, with a median of ~9 mag (a factor of 4000). X-ray hardness ratios for 6 AGNs in our sample and column densities derived from observed dust masses and radio sizes independently confirm this level of attenuation. In these conditions observed optical/near-IR emission comes from surface regions, while inner starburst cores are invisible. We thus attribute the high [NII]/H$\alpha$ ratios to widespread shocks from accretion, turbulence and dynamic disturbances rather than to AGNs. The large range of optical depths demonstrates that substantial diversity is present within the starburst population, possibly connected to different merger phases or progenitor properties. The majority of our targets are, in fact, morphologically classified as mergers. We argue that the extreme obscuration provides in itself smoking gun evidence of their merger origin, and a powerful tool for identifying mergers at even higher redshifts.

The standard AGN-galaxy co-evolutionary scenario predicts a phase of deeply buried supermassive black hole growth coexisting with a starburst (SB) before feedback phenomena deplete the cold molecular gas reservoir of the galaxy and an optically luminous QSO is revealed ('SB-QSO evolutionary sequence'). The aim of this work is to measure the cold gas reservoir of three highly obscured QSOs to test if their gas fraction is similar to that of sub-millimeter galaxies (SMGs), as expected by some models, and place these measurements in the context of the SB-QSO framework. We target CO(1-0) transition in BzK4892, a Compton Thick (CT) QSO at z=2.6, CO(1-0) in BzK8608 and CO(2-1) in CDF153, two highly obscured QSOs at z=2.5 and z=1.5, respectively. For all these targets, we place 3$\sigma$ upper limits on the CO, with $L'_{CO} < (1.5\div 2.8)\times 10^{10}$ K km/s pc$^2$. We also compare the molecular gas conditions of our targets with those of other systems at z>1, considering normal star forming galaxies and SMGs, unobscured and obscured AGN from the literature. For the AGN samples, we provide an updated and (almost) complete collection of targets with CO follow-up. BzK4892 displays a high star formation efficiency (SFE$=L_{IR}/L'_{CO}>410$ L$_{\odot}$/(K km s$^{-1}$ pc$^2$)) and a gas fraction $f_{gas}<0.1$. Less stringent constraints are derived for the other two targets ($f_{gas}<0.5$ and SFE$>10$). From the comparison with literature data, we found that a) obscured AGN at z>1 are associated with higher SFE and lower $f_{gas}$ with respect to star forming galaxies; b) mildly and highly obscured active galaxies have comparable gas fractions; c) the SFE of CT and obscured AGN are similar to those of unobscured AGN. Within the SB-QSO framework, these findings could be consistent with a scenario where feedback can impact the host galaxy already from the early phases of the SB-QSO sequence.

We present ALMA observations of the 870$\mu$m continuum and CO(4-3) line emission in the core of the galaxy cluster ClJ1449+0856 at z=2, a NIR-selected, X-ray detected system in the mass range of typical progenitors of today's massive clusters. The 870$\mu$m map reveals six F$_{870\mu m}$ > 0.5 mJy sources spread over an area of 0.07 arcmin$^2$, giving an overdensity of a factor ~10 (6) with respect to blank field counts down to F$_{870\mu m}$ > 1 (0.5) mJy. On the other hand, deep CO(4-3) follow-up confirms membership of three of these sources, but suggests that the remaining three, including the brightest 870$\mu$m sources in the field (F$_{870\mu m}\gtrsim$2 mJy), are likely interlopers. The measurement of 870$\mu$m continuum and CO(4-3) line fluxes at the positions of previously-known cluster members provides a deep probe of dusty star formation occurring in the core of this high-redshift structure, adding up to a total SFR~700$\pm$100 M$_{\odot}$/yr and yielding an integrated star formation rate density of ~10$^4$ M$_{\odot}$/yr/Mpc$^3$, five orders of magnitude larger than in the field at the same epoch, due to the concentration of star-forming galaxies in the small volume of the dense cluster core. The combination of these observations with previously available HST imaging highlights the presence in this same volume of a population of galaxies with already suppressed star formation. This diverse composition of galaxy populations in ClJ1449+0856 is especially highlighted at the very cluster center, where a complex assembly of quiescent and star-forming sources is likely forming the future Brightest Cluster Galaxy.

We use ALMA and JVLA observations of the galaxy cluster Cl J1449+0856 at z=1.99, in order to study how dust-obscured star-formation, ISM content and AGN activity are linked to environment and galaxy interactions during the crucial phase of high-z cluster assembly. We present detections of multiple transitions of $^{12}$CO, as well as dust continuum emission detections from 11 galaxies in the core of Cl J1449+0856. We measure the gas excitation properties, star-formation rates, gas consumption timescales and gas-to-stellar mass ratios for the galaxies. We find evidence for a large fraction of galaxies with highly-excited molecular gas, contributing $>$50% to the total SFR in the cluster core. We compare these results with expectations for field galaxies, and conclude that environmental influences have strongly enhanced the fraction of excited galaxies in this cluster. We find a dearth of molecular gas in the galaxies' gas reservoirs, implying a high star-formation efficiency (SFE) in the cluster core, and find short gas depletion timescales $\tau$<0.1-0.4 Gyrs for all galaxies. Interestingly, we do not see evidence for increased specific star-formation rates (sSFRs) in the cluster galaxies, despite their high SFEs and gas excitations. We find evidence for a large number of mergers in the cluster core, contributing a large fraction of the core's total star-formation compared with expectations in the field. We conclude that the environmental impact on the galaxy excitations is linked to the high rate of galaxy mergers, interactions and active galactic nuclei in the cluster core.

We have undertaken a comprehensive search for both Lyman Alpha Emitters (LAEs) and Lyman Break Galaxies (LBGs) in the SHARDS Survey of the GOODS-N field. SHARDS is a deep imaging survey, made with the 10.4 m Gran Telescopio Canarias (GTC), employing 25 medium band filters in the range from 500 to 941 nm. This is the first time that both LAEs and LBGs are surveyed simultaneously in a systematic way in a large field. We draw a sample of 1558 sources; 528 of them are LAEs. Most of the sources (1434) show rest-frame UV continua. A minority of them (124) are pure LAEs with virtually no continuum detected in SHARDS. We study these sources from $z\sim3.35$ up to $z\sim6.8$, well into the epoch of reionization. Note that surveys done with just one or two narrow band filters lack the possibility to spot the rest-frame UV continuum present in most of our LAEs. We derive redshifts, Star Formation Rates (SFRs), Ly$\alpha$ Equivalent Widths (EWs) and Luminosity Functions (LFs). Grouping within our sample is also studied, finding 92 pairs or small groups of galaxies at the same redshift separated by less than 60 comoving kpc. In addition, we relate 87 and 55 UV-selected objects with two known overdensities at $z=4.05$ and $z=5.198$, respectively. Finally, we show that surveys made with broad band filters are prone to introduce many unwanted sources ($\sim20$% interlopers), which means that previous studies may be overestimating the calculated LFs, specially at the faint end.

We use high-resolution continuum images obtained at 870microns with the Atacama Large Millimeter Array (ALMA) to probe the surface density of star-formation in z~2 galaxies and study the different physical properties between galaxies within and above the star-formation main sequence of galaxies. This sample of eight star-forming galaxies at z~2 selected among the most massive Herschel galaxies in the GOODS-South field is supplemented with eleven galaxies from the public data of the 1.3 mm survey of the Hubble Ultra-Deep Field. ALMA reveals systematically dense concentrations of dusty star-formation close to the center of the stellar component of the galaxies. We identify two different starburst regimes: (i) the classical population of starbursts located above the SFR-M* main sequence, with enhanced gas fractions and short depletion times and (ii) a sub-population of galaxies located within the scatter of the main sequence that experience compact star formation with depletion timescales typical of starbursts of ~150 Myr. In both starburst populations, the far infrared and UV are distributed in distinct regions and dust-corrected star formation rates estimated using UV-optical-NIR data alone underestimate the total star formation rate. Starbursts hidden in the main sequence show instead the lowest gas fractions of our sample and could represent the last stage of star-formation before they become passive. Being Herschel-selected, these main sequence galaxies are located in the high-mass end of the main sequence, hence we do not know whether these "starbursts hidden in the main sequence" also exist below 10^11 Msun. Active galactic nuclei are found to be ubiquitous in these compact starbursts, suggesting that the triggering mechanism also feeds the central black hole or that the active nucleus triggers star formation.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

[abridged] Building upon a previous study, we define a method to blindly identify galaxies that underwent, and may still be undergoing, a fast downfall of their star-formation activity, that is, a more than 80% drop in star-formation rate (SFR) occurring in less than 500 Myr. Modeling galaxies' SED with a delayed-$\tau$ star formation history (SFH), with and without allowing an instantaneous SFR drop within the last hundreds Myr, we isolate 102 candidates out of a subsample of 6,680 galaxies classified as star-forming from the UVJ criterion in the ZFOURGE catalogues. These galaxies are mostly located in the lower part of the SFR-M$_*$ main sequence (MS) and extend up to a factor 100 below it. They also lie close to the limit between the passive and active regions on the UVJ diagram, indicating that they are in a transition phase. We show that the selected candidates have different physical properties compared to galaxies with similar UVJ colors, namely, lower star-formation rates and different stellar masses. The morphology of the candidates show no preference for a particular type. Among the 102 candidates, only 4 show signs of an AGN activity (from X-ray luminosity or UV-IR SED fitting decomposition). This low fraction of AGNs among the candidates implies that AGN activity may not be the main driver of the recent downfall, although timescale differences and duty cycle must be taken into account. We finally attempt to recover the past position of these galaxies on the SFR-M$_*$ plane, before the downfall of their star-formation and show that some of them were in the starburst region before and are back on the MS. These candidates constitute a promising sample that need more investigation in order to understand the different mechanisms at the origin of the star formation decrease of the Universe since $z$$\sim$2.

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]

The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific X-ray luminosities (i.e. the ratio of $L_X$ to host stellar mass; a common proxy for Eddington ratio). However, recent studies have demonstrated that this universal Eddington ratio distribution fails to reproduce the observed X-ray luminosity functions beyond z$\sim$1.2. Furthermore, empirical studies have recently shown that the Eddington ratio distribution may instead depend upon host galaxy properties, such as SFR and/or stellar mass. To investigate this further we develop a population synthesis model in which the Eddington ratio distribution is different for star-forming and quiescent host galaxies. We show that, although this model is able to reproduce the observed X-ray luminosity functions out to z$\sim$2, it fails to simultaneously reproduce the observed flat relationship between SFR and X-ray luminosity. We can solve this, however, by incorporating a mass dependency in the AGN Eddington ratio distribution for star-forming host galaxies. Overall, our models indicate that a relative suppression of low Eddington ratios ($\lambda_{\rm Edd}\lesssim$0.1) in lower mass galaxies (M<$10^{10-11}$Msun) is required to reproduce both the observed X-ray luminosity functions and the observed flat SFR/X-ray relationship.

We present ALMA observations at 107.291 GHz (band 3) and 214.532 GHz (band 6) of GMASS 0953, a star-forming galaxy at z=2.226 hosting an obscured AGN that has been proposed as a progenitor of compact quiescent galaxies (QG). We measure for the first time the size of the dust and molecular gas emission of GMASS 0953 that we find to be extremely compact ($\sim$1 kpc). This result, coupled with a very high ISM density (n$\sim$10$^{5.5}$ cm$^{-3}$), a low gas mass fraction ($\sim$0.2) and a short gas depletion timescale ($\sim$150 Myr) imply that GMASS 0953 is experiencing an episode of intense star-formation in its central region that will rapidly exhaust its gas reservoirs, likely aided by AGN-induced feedback, confirming its fate as a compact QG. Kinematic analysis of the CO(6-5) line shows evidence of rapidly-rotating gas ($V_{rot}$=320$^{+92}_{-53}$ km s$^{-1}$), as observed also in a handful of similar sources at the same redshift. On-going quenching mechanisms could either destroy the rotation or leave it intact leading the galaxy to evolve into a rotating QG.

We present a new technique to measure multi-wavelength "Super-deblended" photometry from highly confused images, which we apply to Herschel and ground-based far-infrared (FIR) and (sub-)millimeter (mm) data in the northern field of the Great Observatories Origins Deep Survey (GOODS). There are two key novelties. First, starting with a large database of deep Spitzer 24{\mu}m and VLA 20cm detections that are used to define prior positions for fitting the FIR/submm data, we perform an active selection of useful priors independently at each frequency band, moving from less to more confused bands. Exploiting knowledge of redshift and all available photometry, we identify hopelessly faint priors that we remove from the fitting pool. This approach significantly reduces blending degeneracies and allows reliable photometry to be obtained for galaxies in FIR+mm bands. Second, we obtain well-behaved, nearly Gaussian flux density uncertainties, individually tailored to all fitted priors in each band. This is done by exploiting extensive simulations that allow us to calibrate the conversion of formal fitting uncertainties to realistic uncertainties depending on quantities directly measurable. We achieve deeper detection limits with high fidelity measurements and uncertainties at FIR+mm bands. As an illustration of the utility of these measurements, we identify 70 galaxies with z>3 and reliable FIR+mm detections. We present new constraints on the cosmic star formation rate density at 3<z<6, finding a significant contribution from z>3 dusty galaxies that are missed by optical-to-near-infrared color selection. Photometric measurements for 3306 priors, including over 1000 FIR+mm detections are released publicly with our catalog.

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

We report the detection of the CO(3-2) emission line from the host galaxy of Gamma Ray Burst (GRB) 080207 at $z$ = 2.086. This is the first detection of molecular gas in emission from a GRB host galaxy beyond redshift 1. We find this galaxy to be rich in molecular gas with a mass of $1.1 \times 10^{11}\,\rm M_{\odot}$ assuming $\alpha_{\rm CO}=$ 4.36 $\rm M_{\odot}(\rm K\,km\,s^{-1}\,pc^2)^{-1}$. The molecular gas mass fraction of the galaxy is $\sim$ 0.5, typical of star forming galaxies (SFGs) with similar stellar masses and redshifts. With a $\rm SFR_{FIR}$ of 260 $\rm M_{\odot}\,yr^{-1}$, we measure a molecular-gas-depletion timescale of 0.43 Gyr, near the peak of the depletion timescale distribution of SFGs at similar redshifts. Our findings are therefore in contradiction with the proposed molecular gas deficiency in GRB host galaxies. We argue that the reported molecular gas deficiency for GRB hosts could be the artifact of improper comparisons or neglecting the effect of the typical low metallicities of GRB hosts on the CO-to-molecular-gas conversion factor. We also compare the kinematics of the CO(3-2) emission line to that of the H$\alpha$ emission line from the host galaxy. We find the H$\alpha$ emission to have contributions from two separate components, a narrow and a broad one. The narrow component matches the CO emission well in velocity space. The broad component, with a FWHM of $\sim$ 1100 $\rm km\,s^{-1}$, is separated by $+390$ $\rm km\,s^{-1}$ in velocity space from the narrow component. We speculate this broad component to be associated with a powerful outflow in the host galaxy or in an interacting system.

Early type galaxies (ETG) contain most of the stars present in the local Universe and, above a stellar mass of ~5e10 Msun, vastly outnumber spiral galaxies like the Milky Way. These massive spheroidal galaxies have, in the present day, very little gas or dust, and their stellar populations have been evolving passively for over 10 billion years. The physical mechanisms that led to the termination of star formation in these galaxies and depletion of their interstellar medium remain largely conjectural. In particular, there are currently no direct measurements of the amount of residual gas that might be still present in newly quiescent spheroids at high redshift. Here we show that quiescent ETGs at z~1.8, close to their epoch of quenching, contained 2-3 orders of magnitude more dust at fixed stellar mass than local ETGs. This implies the presence of substantial amounts of gas (5-10%), which was however consumed less efficiently than in more active galaxies, probably due to their spheroidal morphology, and consistently with our simulations. This lower star formation efficiency, and an extended hot gas halo possibly maintained by persistent feedback from an active galactic nucleus (AGN), combine to keep ETGs mostly passive throughout cosmic time.

The most heavily-obscured, luminous quasars might represent a specific phase of the evolution of actively accreting supermassive black holes and their host galaxies, possibly related to mergers. We investigated a sample of the most luminous quasars at $z\approx 1-3$ in the GOODS fields, selected in the mid-infrared band through detailed spectral energy distribution (SED) decomposition. The vast majority of these quasars (~80%) are obscured in the X-ray band and ~30% of them to such an extent, that they are undetected in some of the deepest (2 and 4 Ms) Chandra X-ray data. Although no clear relation is found between the star-formation rate of the host galaxies and the X-ray obscuration, we find a higher incidence of heavily-obscured quasars in disturbed/merging galaxies compared to the unobscured ones, thus possibly representing an earlier stage of evolution, after which the system is relaxing and becoming unobscured.

We investigate the relation between the star formation rate (SFR) and the stellar mass (M), i.e. the Main Sequence (MS) relation of star-forming galaxies, at 1.3<=z<6 in the first four HST Frontier Fields, based on rest-frame UV observations. Gravitational lensing combined with deep HST observations allows us to extend the analysis of the MS down to stellar masses as low as logM/Msun~7.5 at z<~4 and logM/Msun~8 at higher redshifts, a factor of ~10 below most previous results. We perform an accurate simulation to take into account the effect of observational uncertainties on the MS and correct for the Eddington bias. This step allows us to reliably measure the MS and in particular its slope. While the normalization increases with redshift, we fit an unevolving and approximately linear slope. We nicely extend to lower masses the results of brighter surveys. Thanks to the large dynamic range in mass and by making use of the simulation, we analyzed any possible mass dependence of the dispersion around the MS. We find tentative evidence that the scatter decreases with increasing stellar masses, suggesting a larger variety of star formation histories in low mass galaxies. This trend agrees with theoretical predictions, and is explained as either a consequence of the smaller number of progenitors of low mass galaxies in a hierarchical scenario and/or of the efficient but intermittent stellar feedback processes in low mass halos. Finally, we observe an increase in the SFR per unit stellar mass with redshift milder than predicted by theoretical models, implying a still incomplete understanding of the processes responsible for galaxy growth.

We estimate the number counts of line emitters at high redshift and their evolution with cosmic time based on a combination of photometry and spectroscopy. We predict the H$\alpha$, H$\beta$, [OII], and [OIII] line fluxes for more than $35,000$ galaxies down to stellar masses of $\sim10^9$ $M_{\odot}$ in the COSMOS and GOODS-S fields, applying standard conversions and exploiting the spectroscopic coverage of the FMOS-COSMOS survey at $z\sim1.55$ to calibrate the predictions. We calculate the number counts of H$\alpha$, [OII], and [OIII] emitters down to fluxes of $1\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the range $1.4 < z < 1.8$ covered by the FMOS-COSMOS survey. We model the time evolution of the differential and cumulative H$\alpha$ counts, steeply declining at the brightest fluxes. We expect $\sim9,300-9,700$ and $\sim2,300-2,900$ galaxies deg$^{-2}$ for fluxes $\geq1\times10^{-16}$ and $\geq2\times10^{-16}$ erg cm$^{-2}$ s$^{-1}$ over the range $0.9<z<1.8$. We show that the observed evolution of the Main Sequence of galaxies with redshift is enough to reproduce the observed counts variation at $0.2<z<2.5$. We characterize the physical properties of the H$\alpha$ emitters with fluxes $\geq2\times10^{-16}$ erg cm$^{-2}$ s$^{-1}$, including their stellar masses, UV sizes, [NII]/H$\alpha$ ratios, and H$\alpha$ equivalent widths. An aperture of $R\sim R_{\rm e}\sim0.5$" maximizes the signal-to-noise ratio for a detection, while causing a factor of $\sim2\times$ flux losses, influencing the recoverable number counts, if neglected. Our approach, based on deep and large photometric datasets, reduces the uncertainties on the number counts due to the selection and spectroscopic samplings, while exploring low fluxes. We publicly release the line flux predictions for the explored photometric samples.

Follow-up observations at high-angular resolution of submillimeter galaxies showed that the single-dish sources are comprised of a blend of several galaxies. Consequently, number counts derived from low and high angular resolution observations are in disagreement. This demonstrates the importance of resolution effects and the need to have realistic simulations to explore them. We built a new 2deg^2 simulation of the extragalactic sky from the far-infrared to the submillimeter. It is based on an updated version of the two star-formation mode galaxy evolution model. Using global galaxy properties, we use the abundance matching technique to populate a dark-matter lightcone and thus simulate the clustering. We produce maps and extract the sources, and show that the limited angular resolution of single-dish instruments have a strong impact on (sub)millimeter continuum observations. Taking into account these resolution effects, we are reproducing a large set of observables, including number counts, redshift distributions, and cosmic infrared background power spectra. Our simulation describes consistently the number counts from single-dish telescopes and interferometers. In particular, at 350 and 500 um, we find that number counts measured by Herschel between 5 and 50 mJy are biased towards high values by a factor 2, and that redshift distributions are biased towards low z. We also show that the clustering has an important impact on the Herschel pixel histogram used to derive number counts from P(D) analysis. Finally, we demonstrate that the large number density of red Herschel sources found in observations but not in models could be an observational artifact caused by the combination of noise, resolution effects, and the steepness of color and flux density distributions. Our simulation (SIDES) is available at http://cesam.lam.fr/sides

We show that the most distant X-ray detected cluster known to date, ClJ1001 at z=2.506, hosts a strong overdensity of radio sources. Six of them are individually detected (within 10") in deep 0.75" resolution VLA 3GHz imaging, with S(3GHz)>8uJy. Of the six, AGN likely affects the radio emission in two galaxies while star formation is the dominant source powering the remaining four. We searched for cluster candidates over the full COSMOS 2-square degree field using radio-detected 3GHz sources and looking for peaks in Sigma5 density maps. ClJ1001 is the strongest overdensity by far with >10sigma, with a simple z_phot>1.5 preselection. A cruder photometric rejection of z<1 radio foregrounds leaves ClJ1001 as the second strongest overdensity, while even using all radio sources ClJ1001 remains among the four strongest projected overdensities. We conclude that there are great prospects for future, deep and wide-area radio surveys to discover large samples of the first generation of forming galaxy clusters. In these remarkable structures widespread star formation and AGN activity of massive galaxy cluster members, residing within the inner cluster core, will ultimately lead to radio continuum as one of the most effective means for their identification, with detection rates expected in the ballpark of 0.1-1 per square degree at z>2.5. Samples of hundreds such high-redshift clusters could potentially constrain cosmological parameters and test cluster and galaxy formation models.

We present deep ALMA CO(5-4) observations of a main sequence, clumpy galaxy at z=1.5 in the HUDF. Thanks to the ~0.5" resolution of the ALMA data, we can link stellar population properties to the CO(5-4) emission on scales of a few kpc. We detect strong CO(5-4) emission from the nuclear region of the galaxy, consistent with the observed $L_{\rm IR}$-$L^{\prime}_{\rm CO(5-4)}$ correlation and indicating on-going nuclear star formation. The CO(5-4) gas component appears more concentrated than other star formation tracers or the dust distribution in this galaxy. We discuss possible implications of this difference in terms of star formation efficiency and mass build-up at the galaxy centre. Conversely, we do not detect any CO(5-4) emission from the UV-bright clumps. This might imply that clumps have a high star formation efficiency (although they do not display unusually high specific star formation rates) and are not entirely gas dominated, with gas fractions no larger than that of their host galaxy (~50%). Stellar feedback and disk instability torques funnelling gas towards the galaxy centre could contribute to the relatively low gas content. Alternatively, clumps could fall in a more standard star formation efficiency regime if their actual star-formation rates are lower than generally assumed. We find that clump star-formation rates derived with several different, plausible methods can vary by up to an order of magnitude. The lowest estimates would be compatible with a CO(5-4) non-detection even for main-sequence like values of star formation efficiency and gas content.

We aim to characterize the relationship between dust properties. We also aim to provide equations to estimate accurate dust properties from limited observational datasets. We assemble a sample of 1,630 nearby (z<0.1) galaxies-over a large range of Mstar, SFR - with multi-wavelength observations available from wise, iras, planck and/or SCUBA. The characterization of dust emission comes from SED fitting using Draine & Li dust models, which we parametrize using two components (warm and cold ). The subsample of these galaxies with global measurements of CO and/or HI are used to explore the molecular and/or atomic gas content of the galaxies. The total Lir, Mdust and dust temperature of the cold component (Tc) form a plane that we refer to as the dust plane. A galaxy's sSFR drives its position on the dust plane: starburst galaxies show higher Lir, Mdust and Tc compared to Main Sequence and passive galaxies. Starburst galaxies also show higher specific Mdust (Mdust/Mstar) and specific Mgas (Mgas/Mstar). The Mdust is more closely correlated with the total Mgas (atomic plus molecular) than with the individual components. Our multi wavelength data allows us to define several equations to estimate Lir, Mdust and Tc from one or two monochromatic luminosities in the infrared and/or sub-millimeter. We estimate the dust mass and infrared luminosity from a single monochromatic luminosity within the R-J tail of the dust emission, with errors of 0.12 and 0.20dex, respectively. These errors are reduced to 0.05 and 0.10 dex, respectively, if the Tc is used. The Mdust is correlated with the total Mism (Mism \propto Mdust^0.7). For galaxies with Mstar 8.5<log(Mstar/Msun) < 11.9, the conversion factor \alpha_850mum shows a large scatter (rms=0.29dex). The SF mode of a galaxy shows a correlation with both the Mgass and Mdust: high Mdust/Mstar galaxies are gas-rich and show the highest SFRs.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

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

We study the properties of dark matter halos that contain star-forming galaxies at $1.43 \le z \le 1.74$ using the FMOS-COSMOS survey. The sample consists of 516 objects with a detection of the H$\alpha$ emission line, that represent the star-forming population at this epoch having a stellar mass range of $10^{9.57}\le M_\ast/M_\odot \lesssim 10^{11.4}$ and a star formation rate range of $15\lesssim \mathrm{SFR}/(M_\odot \mathrm{yr^{-1}}) \lesssim 600$. We measure the projected two-point correlation function while carefully taking into account observational biases, and find a significant clustering amplitude at scales of $0.04$-$10~h^{-1}~\mathrm{cMpc}$, with a correlation length $r_0 = 5.21^{+0.70}_{-0.67}~h^{-1}~\mathrm{cMpc}$ and a bias $b=2.59^{+0.41}_{-0.34}$. We interpret our clustering measurement using a halo occupation distribution model. The sample galaxies appear to reside in halos with mass $M_\mathrm{h} = 4.6^{+1.1}_{-1.6}\times10^{12}~h^{-1}M_\odot$ on average that will likely become present-day halos of mass $M_\mathrm{h} (z=0) \sim2\times10^{13}~h^{-1}M_\odot$, equivalent to the typical halo mass scale of galaxy groups. We then confirm the decline of the stellar-to-halo mass ratio at $M_\mathrm{h}<10^{12}~M_\odot$, finding $M_\ast/M_\mathrm{h} \approx 5\times10^{-3}$ at $M_\mathrm{h}=10^{11.86}~M_\odot$, which is lower by a factor of 2-4 than those measured at higher masses. Finally, we use our results to illustrate the future capabilities of Subaru's Prime-Focus Spectrograph, a next-generation instrument that will provide strong constraints on the galaxy-formation scenario by obtaining precise measurements of galaxy clustering at $z>1$.

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

We report observations of dense molecular gas in the star-forming galaxy EGS 13004291 (z=1.197) using the Plateau de Bure Interferometer. We tentatively detect HCN and HNC (J=2-1) emission when stacked together at ~4sigma significance, yielding line luminosities of L_HCN (J=2-1) =(9 +/- 3) x 10^9 K km s^-1 pc^2 and L_HNC (J=2-1)= (5 +/-2) x 10^9 K km s^-1 pc^2 respectively. We also set 3sigma upper limits of < 7-8 x 10^9 K km s^-1 pc^2 on the HCO+, H2O (3_13-2_20) and HC3N (J=20-19) line luminosities. We serendipitously detect CO emission from two sources at z~1.8 and z~3.2 in the same field of view. We also detect CO(J=2-1) emission in EGS 13004291, showing that the excitation in the previously detected CO(J=3-2) line is subthermal (r_32=0.65 +/- 0.15). We find a line luminosity ratio of L_HCN/L_CO=0.17 +/- 0.07 , as an indicator of the dense gas fraction. This is consistent with the median ratio observed in z>1 galaxies (L_HCN/L_CO=0.16 +/- 0.07) and nearby ULIRGs (L_HCN/L_CO=0.13 +/- 0.03), but higher than in local spirals (L_HCN/L_CO=0.04 +/- 0.02). Although EGS 13004291 lies significantly above the galaxy main sequence at z~1, we do not find an elevated star formation efficiency (traced by L_FIR/L_CO) as in local starbursts, but a value consistent with main-sequence galaxies. The enhanced dense gas fraction, the subthermal gas excitation, and the lower than expected star formation efficiency of the dense molecular gas in EGS 13004291 suggest that different star formation properties may prevail in high-z starbursts. Thus, using L_FIR/L_CO as a simple recipe to measure the star formation efficiency may be insufficient to describe the underlying mechanisms in dense star-forming environments inside the large gas reservoirs.

We explore the dependence of the incidence of moderate-luminosity ($L_{X} = 10^{41.9-43.7}$ erg s$^{-1}$) AGNs and the distribution of their accretion rates on host color at 0.5 < z < 2.5, using deep X-ray data in GOODS fields. We use extinction-corrected rest-frame U-V colors to divide both AGN hosts and non-AGN galaxies into red sequence (quiescent), green valley (transition), and blue cloud (star-forming) populations. We find that both the AGN fraction at fixed stellar mass and its evolution with redshift are dependent on host colors. Most notably, red galaxies have the lowest AGN fraction (~5\%) at z~1 yet with most rapid redshift evolution, increasing by a factor of 5 (~24\%) at z~2. Green galaxies exhibit the highest AGN fraction across all redshifts, which is most pronounced at z~2 with more than half of them hosting an AGN at $M_{*} > 10^{10.6} M_{\odot}$. Together with the high AGN fraction in red galaxies at z~2, this indicates that X-ray AGNs could be important in both transforming blue galaxies into red ones and subsequently maintaining their quiescence at high redshift. Furthermore, consistent with low-redshift studies, we find that the probability of hosting an AGN in the total galaxy population can be characterized by a universal Eddington ratio ($p(\lambda_{Edd}) \sim \lambda_{Edd}^{-0.4}$) and a moderate redshift evolution. Yet consistent with their different AGN fractions, different populations appear to also have different $p(\lambda_{Edd})$ with red galaxies exhibiting more rapid redshift evolution than green and blue ones. Evidence for a steeper power law of $p(\lambda_{Edd})$ in red galaxies is also presented, though larger samples are needed to confirm. These results suggest that the AGN accretion or the growth of supermassive black holes is related to their host properties, and may also influence their hosts in a different mode dependent on the host color.

We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages <1 Gyr is at least 3 times higher in the outer regions of the pBzK galaxies than in their cores. This results in a mild age gradient of ~<0.4 Gyr over ~6 kpc and suggests either the occurrence of widespread rejuvenation episodes or that inside-out quenching played a role in the passivization of this galaxy population. We also report on low-level star formation rates derived from the [OII]3727A emission line, with SFR_OII = 3.7-4.5 Msun/yr. This estimate is confirmed by an independent measurement on a separate sample of similarly-selected quiescent galaxies in the COSMOS field, using stacked far-infrared data (SFR_FIR = 2-4 Msun/yr). This second, photometric sample also displays significant excess at 1.4 GHz, suggestive of the presence of radio-mode AGN activity.

We investigate the physical conditions of ionized gas in high-z star-forming galaxies using diagnostic diagrams based on the rest-frame optical emission lines. The sample consists of 701 galaxies with an Ha detection at $1.4\lesssim z\lesssim1.7$, from the FMOS-COSMOS survey, that represent the normal star-forming population over the stellar mass range $10^{9.6} \lesssim M_\ast/M_\odot \lesssim 10^{11.6}$ with those at $M_\ast>10^{11}~M_\odot$ being well sampled. We confirm an offset of the average location of star-forming galaxies in the BPT diagram ([OIII]/Hb vs. [NII]/Ha), primarily towards higher [OIII]/Hb, compared with local galaxies. Based on the [SII] ratio, we measure an electron density ($n_e=220^{+170}_{-130}~\mathrm{cm^{-3}}$), that is higher than that of local galaxies. Based on comparisons to theoretical models, we argue that changes in emission-line ratios, including the offset in the BPT diagram, are caused by a higher ionization parameter both at fixed stellar mass and at fixed metallicity with additional contributions from a higher gas density and possibly a hardening of the ionizing radiation field. Ionization due to AGNs is ruled out as assessed with Chandra. As a consequence, we revisit the mass-metallicity relation using [NII]/Ha and a new calibration including [NII]/[SII] as recently introduced by Dopita et al. Consistent with our previous results, the most massive galaxies ($M_\ast\gtrsim10^{11}~M_\odot$) are fully enriched, while those at lower masses have metallicities lower than local galaxies. Finally, we demonstrate that the stellar masses, metallicities and star formation rates of the FMOS sample are well fit with a physically-motivated model for the chemical evolution of star-forming galaxies.

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

We use HST/WFC3 imaging to study the red population in the IR-selected, X-ray detected, low-mass cluster Cl J1449+0856 at z=2, one of the few bona-fide established clusters discovered at this redshift, and likely a typical progenitor of an average massive cluster today. This study explores the presence and significance of an early red sequence in the core of this structure, investigating the nature of red sequence galaxies, highlighting environmental effects on cluster galaxy populations at high redshift, and at the same time underlining similarities and differences with other distant dense environments. Our results suggest that the red population in the core of Cl J1449+0856 is made of a mixture of quiescent and dusty star-forming galaxies, with a seedling of the future red sequence already growing in the very central cluster region, and already characterising the inner cluster core with respect to lower density environments. On the other hand, the color-magnitude diagram of this cluster is definitely different from that of lower-redshift (z<1) clusters, as well as of some rare particularly evolved massive clusters at similar redshift, and it is suggestive of a transition phase between active star formation and passive evolution occurring in the proto-cluster and established lower-redshift cluster regimes.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

We report the discovery of a remarkable concentration of massive galaxies with extended X-ray emission at $z_{spec} = 2.506$, which contains 11 massive ($M_{*} \gtrsim 10^{11} M_{\odot}$) galaxies in the central 80kpc region (11.6$\sigma$ overdensity). We have spectroscopically confirmed 17 member galaxies with 11 from CO and the remaining ones from $H\alpha$. The X-ray luminosity, stellar mass content and velocity dispersion all point to a collapsed, cluster-sized dark matter halo with mass $M_{200c} = 10^{13.9\pm0.2} M_{\odot}$, making it the most distant X-ray-detected cluster known to date. Unlike other clusters discovered so far, this structure is dominated by star-forming galaxies (SFGs) in the core with only 2 out of the 11 massive galaxies classified as quiescent. The star formation rate (SFR) in the 80kpc core reaches $\sim$3400 $M_{\odot}$ yr$^{-1}$ with a gas depletion time of $\sim 200$ Myr, suggesting that we caught this cluster in rapid build-up of a dense core. The high SFR is driven by both a high abundance of SFGs and a higher starburst fraction ($\sim25\%$, compared to 3\%-5\% in the field). The presence of both a collapsed, cluster-sized halo and a predominant population of massive SFGs suggests that this structure could represent an important transition phase between protoclusters and mature clusters. It provides evidence that the main phase of massive galaxy passivization will take place after galaxies accrete onto the cluster, providing new insights into massive cluster formation at early epochs. The large integrated stellar mass at such high redshift challenges our understanding of massive cluster formation.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

We present the discovery of a giant $\gtrsim$100~kpc Ly$\alpha$ nebula detected in the core of the X-ray emitting cluster CL~J1449+0856 at $z=1.99$ through Keck/LRIS narrow-band imaging. This detection extends the known relation between Ly$\alpha$ nebulae and overdense regions of the Universe to the dense core of a $5-7\times10^{13}$ M$_{\odot}$ cluster. The most plausible candidates to power the nebula are two Chandra-detected AGN host cluster members, while cooling from the X-ray phase and cosmological cold flows are disfavored primarily because of the high Ly$\alpha$ to X-ray luminosity ratio ($L_{\mathrm{Ly\alpha}}/L_{\mathrm{X}} \approx0.3$, $\gtrsim10-1000\times$ higher than in local cool-core clusters) and by current modeling. Given the physical conditions of the Ly$\alpha$-emitting gas and the possible interplay with the X-ray phase, we argue that the Ly$\alpha$ nebula would be short-lived ($\lesssim10$ Myr) if not continuously replenished with cold gas at a rate of $\gtrsim1000$ M$_{\odot}$ yr$^{-1}$. We investigate the possibility that cluster galaxies supply the required gas through outflows and we show that their total mass outflow rate matches the replenishment necessary to sustain the nebula. This scenario directly implies the extraction of energy from galaxies and its deposition in the surrounding intracluster medium, as required to explain the thermodynamic properties of local clusters. We estimate an energy injection of the order of $\thickapprox2$ keV per particle in the intracluster medium over a $2$ Gyr interval. In our baseline calculation AGN provide up to $85$% of the injected energy and 2/3 of the mass, while the rest is supplied by supernovae-driven winds.

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$.

We use $\sim$83,000 star-forming galaxies at $0.04<z<0.3$ from the Sloan Digital Sky Survey to study the so-called fundamental metallicity relation (FMR) and report on the disappearance of its anti-correlation between metallicity and star formation rate (SFR) when using the new metallicity indicator recently proposed by Dopita et al. In this calibration, metallicity is primarily sensitive to the emission line ratio [NII]$\lambda$6584 / [SII]$\lambda\lambda$6717, 6731 that is insensitive to dilution by pristine infalling gas that may drive the FMR anti-correlation with SFR. Therefore, we conclude that the apparent disappearance of the FMR (using this new metallicity indicator) does not rule out its existence.

We investigate the star-forming properties of 1620 X-ray selected AGN host galaxies as a function of their specific X-ray luminosity (i.e., X-ray luminosity per unit host stellar mass) -- a proxy of the Eddington ratio. Our motivation is to determine whether there is any evidence of a suppression of star-formation at high Eddington ratios, which may hint toward "AGN feedback" effects. Star-formation rates (SFRs) are derived from fits to Herschel-measured far-infrared spectral energy distributions, taking into account any contamination from the AGN. Herschel-undetected AGNs are included via stacking analyses to provide average SFRs in bins of redshift and specific X-ray luminosity (spanning $0.01 \lesssim L_{\rm X}/M_{\ast} \lesssim 100~L_{\odot} ~M_{\odot}^{-1}$). After normalising for the effects of mass and redshift arising from the evolving galaxy main sequence, we find that the SFRs of high specific luminosity AGNs are slightly enhanced compared to their lower specific luminosity counterparts. This suggests that the SFR distribution of AGN hosts changes with specific X-ray luminosity, a result reinforced by our finding of a significantly higher fraction of starbursting hosts among high specific luminosity AGNs compared to that of the general star-forming galaxy population (i.e., 8-10 per cent vs. 3 per cent). Contrary to our original motivation, our findings suggest that high specific luminosity AGNs are more likely to reside in galaxies with enhanced levels of star-formation.

We study the relationship between stellar mass, star formation rate (SFR),ionization state, and gas-phase metallicity for a sample of 41 normal star-forming galaxies at $3 \lesssim z \lesssim 3.7$. The gas-phase oxygen abundance, ionization parameter, and electron density of ionized gas are derived from rest-frame optical strong emission lines measured on near-infrared spectra obtained with Keck/MOSFIRE. We remove the effect of these strong emission lines in the broad-band fluxes to compute stellar masses via spectral energy distribution fitting, while the SFR is derived from the dust-corrected ultraviolet luminosity. The ionization parameter is weakly correlated with the specific SFR, but otherwise the ionization parameter and electron density do not correlate with other global galaxy properties such as stellar mass, SFR, and metallicity. The mass-metallicity relation (MZR) at $z\simeq3.3$ shows lower metallicity by $\simeq 0.7$ dex than that at $z=0$ at the same stellar mass. Our sample shows an offset by $\simeq 0.3$ dex from the locally defined mass-metallicity-SFR relation, indicating that simply extrapolating such relation to higher redshift may predict an incorrect evolution of MZR. Furthermore, within the uncertainties we find no SFR-metallicity correlation, suggesting a less important role of SFR in controlling the metallicity at high redshift. We finally investigate the redshift evolution of the MZR by using the model by Lilly et al. (2013), finding that the observed evolution from $z=0$ to $z\simeq3.3$ can be accounted for by the model assuming a weak redshift evolution of the star formation efficiency.

In this paper we study the causes of the reported mass-dependence of the slope of SFR-M* relation, the so-called "Main Sequence" of star-forming galaxies, and discuss its implication on the physical processes that shaped the star formation history of massive galaxies over cosmic time. We use the CANDELS near-IR imaging from the Hubble Space Telescope to perform the bulge-to-disk decomposition of distant galaxies and measure for the first time the slope of the SFR-Mdisk relation at z=1. We find that this relation follows very closely the shape of the SFR-M* correlation, still with a pronounced flattening at the high-mass end. This is clearly excluding, at least at z=1, the secular growth of quiescent bulges in star-forming galaxies as the main driver for the change of slope of the Main Sequence. Then, by stacking the Herschel data available in the CANDELS field, we estimate the total gas mass and the star formation efficiency at different positions on the SFR-M* relation. We find that the relatively low SFRs observed in massive galaxies (M* > 5e10 Msun) are caused by a decreased star formation efficiency, by up to a factor of 3 as compared to lower stellar mass galaxies, and not by a reduced gas content. The trend at the lowest masses is likely linked to the dominance of atomic over molecular gas. We argue that this stellar-mass-dependent SFE can explain the varying slope of the Main Sequence since z=1.5, hence over 70% of the Hubble time. The drop of SFE occurs at lower masses in the local Universe (M* > 2e10 Msun) and is not present at z=2. Altogether this provides evidence for a slow downfall of the star formation efficiency in massive Main Sequence galaxies. The resulting loss of star formation is found to be rising starting from z=2 to reach a level comparable to the mass growth of the quiescent population by z=1. We finally discuss the possible physical origin of this phenomenon.

We present high resolution (0.3") Atacama Large Millimeter Array (ALMA) 870um imaging of five z~1.5-4.5 X-ray detected AGN with luminosities of L(2-8keV)>10^42 erg/s. These data provide a >~20x improvement in spatial resolution over single-dish rest-frame FIR measurements. The sub-millimetre emission is extended on scales of FWHM~0.2"-0.5", corresponding to physical sizes of 1-3 kpc (median value of 1.8 kpc). These sizes are comparable to the majority of z=1-5 sub-millimetre galaxies (SMGs) with equivalent ALMA measurements. In combination with spectral energy distribution analyses, we attribute this rest-frame far-infrared (FIR) emission to dust heated by star formation. The implied star-formation rate surface densities are ~20-200 Mo/yr/kpc^2, which are consistent with SMGs of comparable FIR luminosities (i.e., L(IR)~[1-5]x10^(12)Lo). Although limited by a small sample of AGN, which all have high FIR luminosities, our study suggests that the kpc-scale spatial distribution and surface density of star formation in high-redshift star-forming galaxies is the same irrespective of the presence of X-ray detected AGN.

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

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

We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars ($\nu$L(6 micron)>6x10$^{44}$ erg/s) at redshift z~1-3, identified through detailed spectral energy distribution analyses of distant star-forming galaxies, using the deepest IR data from Spitzer and Herschel in the GOODS-Herschel fields. The aim is to constrain the fraction of obscured, and Compton-thick (CT, N$_H$>1.5x10$^{24}$ cm$^{-2}$) quasars at the peak era of nuclear and star-formation activities. Despite being very bright in the mid-IR band, ~30% of these quasars are not detected in the extremely deep 2 Ms and 4 Ms Chandra X-ray data available in these fields. X-ray spectral analysis of the detected sources reveals that the majority (~67%) are obscured by column densities N$_H$>10$^{22}$ cm$^{-2}$; this fraction reaches ~80% when including the X-ray undetected sources (9 out of 33), which are likely to be the most heavily-obscured, CT quasars. We constrain the fraction of CT quasars in our sample to be ~24-48%, and their space density to be $\Phi$=(6.7$\pm$2.2)x10$^{-6}$ Mpc$^{-3}$. From the investigation of the quasar host galaxies in terms of star-formation rates (SFRs) and morphological distortions, as a sign of galaxy mergers/interactions, we do not find any direct relation between SFRs and quasar luminosity or X-ray obscuration. On the other hand, there is tentative evidence that the most heavily-obscured quasars have, on average, more disturbed morphologies than the unobscured/moderately-obscured quasar hosts, which preferentially live in undisturbed systems. However, the fraction of quasars with disturbed morphology amongst the whole sample is ~40%, suggesting that galaxy mergers are not the main fuelling mechanism of quasars at z~2.

We study the debated contribution from thermally pulsing asymptotic giant branch (TP-AGB) stars in evolutionary population synthesis models. We investigate the Spectral Energy Distributions (SEDs) of a sample of 51 spectroscopically confirmed, high-z ($1.3<z_{\rm spec}<2.7$), galaxies using three evolutionary population synthesis models with strong, mild and light TP-AGB. Our sample is the largest of spectroscopically confirmed galaxies on which such models are tested so far. Galaxies were selected as passive, but we model them using a variety of star formation histories in order not to be dependent on this pre-selection. We find that the observed SEDs are best fitted with a significant contribution of TP-AGB stars or with substantial dust attenuation. Without including reddening, TP-AGB-strong models perform better and deliver solutions consistent within $1\sigma$ from the best-fit ones in the vast majority of cases. Including reddening, all models perform similarly. Using independent constraints from observations in the mid- and far-IR, we show that low/negligible dust attenuation, i.e. $E(B-V)\lesssim 0.05$ , should be preferred for the SEDs of passively-selected galaxies. Given that TP-AGB-light models give systematically older ages for passive galaxies, we suggest number counts of passive galaxies at higher redshifts as a further test to discriminate among stellar population models.

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

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermediate redshifts (z~1). We combine the ultra-deep optical spectro-photometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) with deep UV-to-FIR observations in the GOODS-N field. Exploiting two of the 25 SHARDS medium-band filters, F687W17 and F823W17, we select [OII] emission line galaxies at z~0.84 and z~1.23 and characterize their physical properties. Their rest-frame equivalent widths (EW$_{\mathrm{rf}}$([OII])), line fluxes, luminosities, star formation rates (SFRs) and dust attenuation properties are investigated. The evolution of the EW$_{\mathrm{rf}}$([OII]) closely follows the SFR density evolution of the universe, with a trend of EW$_{\mathrm{rf}}$([OII])$\propto$(1+z)$^3$ up to redshift z~1, followed by a possible flattening. The SF properties of the galaxies selected on the basis of their [OII] emission are compared with complementary samples of SFGs selected by their MIR and FIR emission, and also with a general mass-selected sample of galaxies at the same redshifts. We demonstrate observationally that the UVJ diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent galaxies from the SFGs. The SFR-M$_*$ relation is investigated for the different samples, yelding a logarithmic slope ~1, in good agreement with previous results. The dust attenuations derived from different SFR indicators (UV(1600), UV(2800), [OII], IR) are compared and show clear trends with respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies being affected by stronger dust attenuation.

We present correlations between 9 CO transition ($J=4-3$ to $12-11$) and beam-matched far-infrared (Far-IR) luminosities ($L_{\mathrm{FIR},\,b}$) among 167 local galaxies, using {\it{Herschel}} Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE; FTS) spectroscopic data and Photoconductor Array Camera and Spectrometer (PACS) photometry data. We adopt entire-galaxy FIR luminosities ($L_{\mathrm{FIR},\,e}$) from the {\it{IRAS}} Revised Bright Galaxy Sample and correct to $L_{\mathrm{FIR},\,b}$ using PACS images to match the varying FTS beam sizes. All 9 correlations between $L'_{\mathrm{CO}}$ and $L_{\mathrm{FIR},\,b}$ are essentially linear and tight ($\sigma$=0.2-0.3 dex dispersion), even for the highest transition, $J=12-11$. This supports the notion that the star formation rate (SFR) is linearly correlated with the dense molecular gas ($n_{\mathrm{H}_2}\gtrsim10^{4-6}\,cm^{-3}$). We divide the entire sample into three subsamples and find that smaller sample sizes can induce large differences in the correlation slopes. We also derive an average CO spectral line energy distribution (SLED) for the entire sample and discuss the implied average molecular gas properties for these local galaxies. We further extend our sample to high-{\it{z}} galaxies with CO $J=5-4$ data from the literature as an example, including submillimeter galaxies (SMGs) and "normal" star-forming BzKs. BzKs have similar FIR/CO(5-4) ratios as that of local galaxies, an follow well the locally-determined correlation, whereas SMG ratios fall around or slightly above the local correlation with large uncertainties. Finally, by including Galactic CO($J=10-9$) data as well as very limited high-{\it{z}} CO $J=10-9$ data, we verify that the CO(10-9) -- FIR correlation successfully extends to Galactic young stellar objects, suggesting that linear correlations are valid over 15 orders of magnitude.

We present the first estimate of age, stellar metallicity and chemical abundance ratios, for an individual early-type galaxy at high-redshift (z = 1.426) in the COSMOS field. Our analysis is based on observations obtained with the X-Shooter instrument at the VLT, which cover the visual and near infrared spectrum at high (R >5000) spectral resolution. We measure the values of several spectral absorptions tracing chemical species, in particular Magnesium and Iron, besides determining the age-sensitive D4000 break. We compare the measured indices to stellar population models, finding good agreement. We find that our target is an old (t > 3 Gyr), high-metallicity ([Z/H] > 0.5) galaxy which formed its stars at z_{form} > 5 within a short time scale ~0.1 Gyr, as testified by the strong [\alpha/Fe] ratio ( > 0.4), and has passively evolved in the first > 3-4 Gyr of its life. We have verified that this result is robust against the choice and number of fitted spectral features, and stellar population model. The result of an old age and high-metallicity has important implications for galaxy formation and evolution confirming an early and rapid formation of the most massive galaxies in the Universe.

We use a spectroscopic sample of 286 star-forming galaxies (SFGs) at 1<z<3 from the GMASS survey to study different star formation rate (SFR) estimators. Infrared (IR) data are used to derive empirical calibrations to correct ultraviolet (UV) and [OII]{\lambda}3727 luminosities for dust extinction and dust-corrected estimates of SFR. In the selection procedure we fully exploit the available spectroscopic information. On the basis of three continuum indices, we are able to identify and exclude from the sample galaxies in which old stellar populations might bring a non-negligible contribution to IR luminosity (LIR) and continuum reddening. Using Spitzer-MIPS and Herschel-PACS data we derive LIR for two-thirds of our sample. The LIR/LUV ratio is used as a probe of effective attenuation (AIRX) to search for correlations with continuum and spectroscopic features. The relation between AIRX and UV continuum slope ({\beta}) was tested for our sample and found to be broadly consistent with the literature results at the same redshift, though with a larger dispersion with respect to UV-selected samples. We find a correlation between the rest-frame equivalent width (EW) of the [OII]{\lambda}3727 line and {\beta}, which is the main result of this work. We therefore propose the [OII]{\lambda}3727 line EW as a dust attenuation probe and calibrate it through AIRX, though the assumption of a reddening curve is still needed to derive the actual attenuation towards the [OII]{\lambda}3727 line. We tested the issue of differential attenuation towards stellar continuum and nebular emission: our results are in line with the traditional prescription of extra attenuation towards nebular lines. A set of relations is provided that allows the recovery of the total unattenuated SFR from UV and [OII]{\lambda}3727 luminosities. (Abridged)

Using deep Herschel and ALMA observations, we investigate the star formation rate (SFR) distributions of X-ray selected AGN host galaxies at 0.5<z<1.5 and 1.5<z<4, comparing them to that of normal, star-forming (i.e., "main-sequence", or MS) galaxies. We find 34--55 per cent of AGNs in our sample have SFRs at least a factor of two below that of the average MS galaxy, compared to ~15 per cent of all MS galaxies, suggesting significantly different SFR distributions. Indeed, when both are modelled as log-normal distributions, the mass and redshift-normalised SFR distributions of X-ray AGNs are roughly twice as broad, and peak ~0.4 dex lower, than that of MS galaxies. However, like MS galaxies, the normalised SFR distribution of AGNs in our sample appears not to evolve with redshift. Despite X-ray AGNs and MS galaxies having different SFR distributions, the linear-mean SFR of AGNs derived from our distributions is remarkably consistent with that of MS galaxies, and thus with previous results derived from stacked Herschel data. This apparent contradiction is due to the linear-mean SFR being biased by bright outliers, and thus does not necessarily represent a true characterisation of the typical SFR of X-ray AGNs.

We present a spectroscopic survey of galaxies in the COSMOS field using the Fiber Multi-Object Spectrograph (FMOS), a near-infrared instrument on the Subaru Telescope. Our survey is specifically designed to detect the Halpha emission line that falls within the H-band (1.6-1.8 um) spectroscopic window from star-forming galaxies with 1.4 < z < 1.7 and M_stellar>~10^10 Msolar. With the high multiplex capability of FMOS, it is now feasible to construct samples of over one thousand galaxies having spectroscopic redshifts at epochs that were previously challenging. The high-resolution mode (R~2600) effectively separates Halpha and [NII]6585 thus enabling studies of the gas-phase metallicity and photoionization state of the interstellar medium. The primary aim of our program is to establish how star formation depends on stellar mass and environment, both recognized as drivers of galaxy evolution at lower redshifts. In addition to the main galaxy sample, our target selection places priority on those detected in the far-infrared by Herschel/PACS to assess the level of obscured star formation and investigate, in detail, outliers from the star formation rate - stellar mass relation. Galaxies with Halpha detections are followed up with FMOS observations at shorter wavelengths using the J-long (1.11-1.35 um) grating to detect Hbeta and [OIII]5008 that provides an assessment of extinction required to measure star formation rates not hampered by dust, and an indication of embedded Active Galactic Nuclei. With 460 redshifts measured from 1153 spectra, we assess the performance of the instrument with respect to achieving our goals, discuss inherent biases in the sample, and detail the emission-line properties. Our higher-level data products, including catalogs and spectra, are available to the community.

We study the properties of satellites in the environment of massive star-forming galaxies at z~1.8 in the COSMOS field, using a sample of 215 galaxies on the main sequence of star formation with an average mass of 10^11 Msun. At z>1.5, these galaxies typically trace halos of mass >10^13 Msun. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of centrals and satellites down to ~6*10^9 Msun. We stack data around 215 central galaxies to statistically detect their satellite halos, finding an average of ~3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple beta-models, and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to 68% of the central galaxy, while SED modeling and far-infrared photometry consistently show their total SFR to be 25-35% of the central's rate. We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at <100 kpc. After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut down star formation in satellites on relatively short timescales of ~0.3 Gyr.

We investigate the stellar population properties of a sample of 24 massive quenched galaxies at $1.25<z_\mathrm{spec}<2.09$ identified in the COSMOS field with our Subaru/MOIRCS near-IR spectroscopic observations. Tracing the stellar population properties as close to their major formation epoch as possible, we try to put constraints on the star formation history, post-quenching evolution, and possible progenitor star-forming populations for such massive quenched galaxies. By using a set of Lick absorption line indices on a rest-frame optical composite spectrum, the average age, metallicity [Z/H], and $\alpha$-to-iron element abundance ratio [$\alpha$/Fe] are derived as $\log(\mathrm{age}/\mathrm{Gyr})=0.04_{-0.08}^{+0.10}$, $\mathrm{[Z/H]}=0.24_{-0.14}^{+0.20}$, and $[\alpha/\mathrm{Fe}]=0.31_{-0.12}^{+0.12}$, respectively. If our sample of quenched galaxies at $\langle z \rangle = 1.6$ is evolved passively to $z=0$, their stellar population properties will align in excellent agreement with local counterparts at similar stellar velocity dispersions, which qualifies them as progenitors of local massive early-type galaxies. Redshift evolution of stellar population ages in quenched galaxies combined with low redshift measurements from the literature suggests a formation redshift of $z_\mathrm{f} \sim 2.3$ around which the bulk of stars in these galaxies have been formed. The measured [$\alpha$/Fe] value indicates a star formation timescale of $\lesssim 1$ Gyr, which can be translated into a specific star formation rate of $\simeq 1\,\mathrm{Gyr}^{-1}$ prior to quenching. Based on these findings, we discuss identifying possible progenitor star-forming galaxies at $z \simeq 2.3$. We identify normal star-forming galaxies, i.e, those on the star-forming main sequence, followed by a rapid quenching event, as likely precursors of the quenched galaxies at $\langle z \rangle = 1.6$ presented here.

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

Gas and dust in star-forming galaxies at the peak epoch of galaxy assembly are presently the topic of intense study, but little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular (H$_2$) gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$\sigma$ upper limit of 0.32 Jy km/s on the flux of the CO($J$=2$\rightarrow$1) line or $L'_{\rm CO}$$<$8.8$\times$10$^{9}$ K km/s pc$^2$, assuming a disk-like CO-morphology and a circular velocity scaling with the stellar velocity dispersion as in local early-type galaxies (ETGs). This translates to an H$_2$ mass of $<$3.9$\times$10$^{10}$($\alpha_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$6% assuming a Salpeter initial mass function (IMF) and an ISM dominated by H$_2$, as observed in many local, high-mass ellipticals. This low value approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a lesser role for inhibitive mechanisms like morphological quenching.

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.

At z>1, the distinction between merging and 'normal' star-forming galaxies based on single band morphology is often hampered by the presence of large clumps which result in a disturbed, merger-like appearance even in rotationally supported disks. In this paper we discuss how a classification based on canonical, non-parametric structural indices measured on resolved stellar mass maps, rather than on single-band images, reduces the misclassification of clumpy but not merging galaxies. We calibrate the mass-based selection of mergers using the MIRAGE hydrodynamical numerical simulations of isolated and merging galaxies which span a stellar mass range of $10^{9.8}$-$10^{10.6}M_{sun}$ and merger ratios between 1:1-1:6.3. These simulations are processed to reproduce the typical depth and spatial resolution of observed HUDF data. We test our approach on a sample of real z~2 galaxies with kinematic classification into disks or mergers and on ~100 galaxies in the HUDF field with photometric/spectroscopic redshift between 1.5$\leqslant z \leqslant$3 and $M>10^{9.4}M_{sun}$. We find that a combination of the asymmetry $A_{MASS}$ and $M_{20, MASS}$ indices measured on the stellar mass maps can efficiently identify real (major) mergers with $\lesssim$20% contamination from clumpy disks in the merger sample. This mass-based classification cannot be reproduced in star-forming galaxies by $H-$band measurements alone, which instead result in a contamination from clumpy galaxies that can be as high as 50%. Moreover, we find that the mass-based classification always results in a lower contamination from clumpy galaxies than an $H-$band classification, regardless of the depth of the imaging used (e.g., CANDELS vs. HUDF).

We have used FMOS on Subaru to obtain near-infrared spectroscopy of 123 far-infrared selected galaxies in COSMOS and obtain the key rest-frame optical emission lines. This is the largest sample of infrared galaxies with near-infrared spectroscopy at these redshifts. The far-infrared selection results in a sample of galaxies that are massive systems that span a range of metallicities in comparison with previous optically selected surveys, and thus has a higher AGN fraction and better samples the AGN branch. We establish the presence of AGN and starbursts in this sample of (U)LIRGs selected as Herschel-PACS and Spitzer-MIPS detections in two redshift bins (z~0.7 and z~1.5) and test the redshift dependence of diagnostics used to separate AGN from star-formation dominated galaxies. In addition, we construct a low redshift (z~0.1) comparison sample of infrared selected galaxies and find that the evolution from z~1.5 to today is consistent with an evolving AGN selection line and a range of ISM conditions and metallicities from the models of Kewley et al. (2013b). We find that a large fraction of (U)LIRGs are BPT-selected AGN using their new, redshift-dependent classification line. We compare the position of known X-ray detected AGN (67 in total) with the BPT selection and find that the new classification line accurately selects most of these objects (> 70%). Furthermore, we identify 35 new (likely obscured) AGN not selected as such by their X-ray emission. Our results have direct implications for AGN selection at higher redshift with either current (MOSFIRE, KMOS) or future (PFS, MOONS) spectroscopic efforts with near-infrared spectral coverage.

When the cosmic star formation history peaks (z ~ 2), galaxies vigorously fed by cosmic reservoirs are gas dominated and contain massive star-forming clumps, thought to form by violent gravitational instabilities in highly turbulent gas-rich disks. However, a clump formation event has not been witnessed yet, and it is debated whether clumps survive energetic feedback from young stars, thus migrating inwards to form galaxy bulges. Here we report spatially resolved spectroscopy of a bright off-nuclear emission line region in a galaxy at z = 1.987. Although this region dominates the star formation in the galaxy disk, its stellar continuum remains undetected in deep imaging, revealing an extremely young (age < 10 Myr) massive clump, forming through the gravitational collapse of > 10$^9$ M$_{\odot}$ of gas. Gas consumption in this young clump is > 10 times faster than in the host galaxy, displaying high star formation efficiency during this phase, in agreement with our hydrodynamic simulations. The frequency of older clumps with similar masses coupled with our initial estimate of their formation rate (~ 2.5 Gyr$^{-1}$) supports long lifetimes (~ 500 Myr), favouring scenarios where clumps survive feedback and grow the bulges of present-day galaxies.

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

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

We interpret the large variety of redshift distributions of galaxies found by far-infrared and (sub-)millimeter deep surveys depending on their depth and wavelength using the B\'ethermin et al. (2012) phenomenological model of galaxy evolution. This model reproduces without any new parameter tuning the observed redshift distributions from 100 $\mu$m to 1.4 mm, and especially the increase of the median redshift with survey wavelength. This median redshift varies also significantly with the depth of the surveys, and deeper surveys do necessarily not probe higher redshifts. Paradoxically, at fixed wavelength and flux limit, the lensed sources are not always at higher redshift. We found that the higher redshift of 1.4 mm-selected south pole telescope (SPT) sources compared to other SMG surveys is not only caused by the lensing selection, but also by the longer wavelength. This SPT sample is expected to be dominated by a population of lensed main-sequence galaxies and a minor contribution ($\sim$10\%) of unlensed extreme starbursts.

Even 10 billion years ago, the cores of the first galaxy clusters are often found to host a characteristic population of massive galaxies with already suppressed star formation. Here we search for distant cluster candidates at z~2 using massive passive galaxies as tracers. With a sample of ~40 spectroscopically confirmed passive galaxies at 1.3<z<2.1, we tune photometric redshifts of several thousands passive sources in the full 2 sq.deg. COSMOS field. This allows us to map their density in redshift slices, probing the large scale structure in the COSMOS field as traced by passive sources. We report here on the three strongest passive galaxy overdensities that we identify in the redshift range 1.5<z<2.5. While the actual nature of these concentrations is still to be confirmed, we discuss their identification procedure, and the arguments supporting them as candidate galaxy clusters (likely mid-10^13 M_sun range). Although this search approach is likely biased towards more evolved structures, it has the potential to select still rare, cluster-like environments close to their epoch of first appearance, enabling new investigations of the evolution of galaxies in the context of structure growth.

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

We investigate the CO excitation of normal star forming disk galaxies at z=1.5 using IRAM PdBI observations of the CO[2-1], CO[3-2] and CO[5-4] transitions for 4 galaxies, including VLA observations of CO[1-0] for 3 of them, with the aim of constraining the average state of H2 gas. Exploiting prior knowledge of the velocity range, spatial extent and size of the CO emission we measure reliable line fluxes with S/N>4-7 for individual transitions. While the average CO Spectral Line Energy Distribution (SLED) has a sub-thermal excitation similar to the Milky Way (MW) up to CO[3-2], we show that the average CO[5-4] emission is four times stronger than assuming MW excitation. This demonstrates the presence of an additional component of more excited, denser and possibly warmer molecular gas. The ratio of CO[5-4] to lower-J CO emission is lower than in local (U)LIRGs and high-redshift SMGs, correlating closely with the average intensity of the radiation field <U> and with the star formation surface density, but not with the SF efficiency (SFE). The CO[5-4] luminosity correlates linearly with LIR over 4 orders of magnitudes, with z=1.5 BzK galaxies following the same trend as local spirals and (U)LIRGs and high redshift star bursting SMGs. The CO[5-4] luminosity is thus empirically related to the dense gas, and might be a more convenient way to probe it than standard high--density tracers that are much fainter than CO. We see excitation variations among our sample galaxies, that can be linked to their evolutionary state and clumpiness in optical rest frame images. In one galaxy we see spatially resolved excitation variations, where the more highly excited part of the galaxy corresponds to the location of massive SF clumps. This supports to models that suggest that giant clumps are the main source of the high excitation CO emission in high redshift disk-like galaxies.

We investigate the environmental effect on the metal enrichment of star-forming galaxies (SFGs) in the farthest spectroscopically confirmed and X-ray detected cluster, CL~J1449+0856 at $z=1.99$. We combined HST/WFC3~G141 slitless spectroscopic data, our 13-band photometry, and a recent Subaru/MOIRCS near infrared spectroscopic follow-up to constrain the physical properties of SFGs in CL~J1449+0856 and in a mass-matched field sample. After a conservative active galactic nuclei (AGN) removal, stacking individual MOIRCS spectra of $6$ ($31$) sources in the cluster (field) in the mass range $10 \leq \mathrm{log(M/M_{\odot}}) \leq 11$, we find a $\sim4\sigma$ significant lower [N II]/H$\alpha$ ratio in the cluster than in the field. Stacking a subsample of $16$ field galaxies with H$\beta$ and [O III] in the observed range, we measure a [O III]/H$\beta$ ratio fully compatible with the cluster value. Converting these ratios into metallicities, we find that the cluster SFGs are up to $0.25\,\mathrm{dex}$ poorer in metals than their field counterparts, depending on adopted calibrations. The low metallicity in cluster sources is confirmed using alternative indicators. Furthermore, we observe a significantly higher H$\alpha$ luminosity and equivalent width in the average cluster spectrum than in the field. This is likely due to enhanced specific star formation rate, even if lower dust reddening and/or an uncertain environmental dependence of the continuum-to-nebular emission differential reddening may play a role. Our findings might be explained by the accretion of pristine gas around galaxies at $z=2$ and from cluster-scale reservoirs, possibly connected with a phase of rapid halo mass assembly at $z>2$ and of high galaxy merging rate.

We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with very high resolutions (3pc and 10^3 Msun) and detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes and galactic winds. The simulations are analyzed with a Large Velocity Gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer and opacity effects, and the intensity emerging from galaxies, to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially-concentrated stellar feedback sources, although CO(1-0) to CO(3-2) emission is overall dominated by low-excitation gas around the densest clumps. These results overall highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn impact the CO emission properties.

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.

We present high-resolution observations of the 880 $\mu$m (rest-frame FIR) continuum emission in the z$=$4.05 submillimeter galaxy GN20 from the IRAM Plateau de Bure Interferometer (PdBI). These data resolve the obscured star formation in this unlensed galaxy on scales of 0.3$^{\prime\prime}$$\times$0.2$^{\prime\prime}$ ($\sim$2.1$\times$1.3 kpc). The observations reveal a bright (16$\pm$1 mJy) dusty starburst centered on the cold molecular gas reservoir and showing a bar-like extension along the major axis. The striking anti-correlation with the HST/WFC3 imaging suggests that the copious dust surrounding the starburst heavily obscures the rest-frame UV/optical emission. A comparison with 1.2 mm PdBI continuum data reveals no evidence for variations in the dust properties across the source within the uncertainties, consistent with extended star formation, and the peak star formation rate surface density (119$\pm$8 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$) implies that the star formation in GN20 remains sub-Eddington on scales down to 3 kpc$^2$. We find that the star formation efficiency is highest in the central regions of GN20, leading to a resolved star formation law with a power law slope of $\Sigma_{\rm SFR}$ $\sim$ $\Sigma_{\rm H_2}^{\rm 2.1\pm1.0}$, and that GN20 lies above the sequence of normal star-forming disks, implying that the dispersion in the star formation law is not due solely to morphology or choice of conversion factor. These data extend previous evidence for a fixed star formation efficiency per free-fall time to include the star-forming medium on $\sim$kpc-scales in a galaxy 12 Gyr ago.

The broad spectral bandwidth at mm and cm-wavelengths provided by the recent upgrades to the Karl G. Jansky Very Large Array (VLA) has made it possible to conduct unbiased searches for molecular CO line emission at redshifts, z > 1.31. We present the discovery of a gas-rich, star-forming galaxy at z = 2.48, through the detection of CO(1-0) line emission in the COLDz survey, through a sensitive, Ka-band (31 to 39 GHz) VLA survey of a 6.5 square arcminute region of the COSMOS field. We argue that the broad line (FWHM ~570 +/- 80 km/s) is most likely to be CO(1-0) at z=2.48, as the integrated emission is spatially coincident with an infrared-detected galaxy with a photometric redshift estimate of z = 3.2 +/- 0.4. The CO(1-0) line luminosity is L'_CO = (2.2 +/- 0.3) x 10^{10} K km/s pc^2, suggesting a cold molecular gas mass of M_gas ~ (2 - 8)x10^{10}M_solar depending on the assumed value of the molecular gas mass to CO luminosity ratio alpha_CO. The estimated infrared luminosity from the (rest-frame) far-infrared spectral energy distribution (SED) is L_IR = 2.5x10^{12} L_solar and the star-formation rate is ~250 M_solar/yr, with the SED shape indicating substantial dust obscuration of the stellar light. The infrared to CO line luminosity ratio is ~114+/-19 L_solar/(K km/s pc^2), similar to galaxies with similar SFRs selected at UV/optical to radio wavelengths. This discovery confirms the potential for molecular emission line surveys as a route to study populations of gas-rich galaxies in the future.

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged).

Star-forming galaxies (SFGs) display a continuous distribution of specific star formation rates (sSFR) which can be approximated by the superposition of two log-normal distributions. The 1st of these encompasses the main sequence (MS) of SFGs, the 2nd one a rarer population of starbursts (SB). We show that the sSFR-distribution of SBs can be regarded as the result of a physical process (plausibly merging) taking the mathematical form of a log-normal boosting kernel and enhancing star formation activity. We explore the utility of splitting the star-forming population into MS and SB galaxies - an approach we term "2-Star Formation Mode" (2-SFM) framework - for understanding their molecular gas properties. Variations of star formation efficiency (SFE) and gas fractions among SFGs take a simple, redshift-independent form, once these quantities are normalized to the value of an average MS galaxy. The change in SFE for galaxies undergoing a starburst event scales supra-linearly with the SFR-increase, as expected for mergers. This implies a continuous distribution of galaxies in the Schmidt-Kennicutt plane that separates more clearly into loci for SBs and normal galaxies than observed in SFR vs. M* space. SBs with the largest deviations (>10-fold) from the MS, like many local ULIRGs, are not average SBs, but even rarer events with progenitors having larger gas fractions than typical MS galaxies. We statistically infer that the gas fractions of typical SBs decrease by a factor of 2 to 3 with respect to their direct MS progenitors, as expected to occur in short-term SFR-boosts during which internal gas reservoirs are drained more quickly than gas is accreted from the cosmic web. We predict variations of the CO-to-H2 conversion factor in the SFR-M* plane and provide evidence that the higher sSFR of distant galaxies is a direct consequence of larger gas fraction in these systems.

We measure the dust and gas content of the three sub-millimeter galaxies (SMGs) in the GN20 proto-cluster at $z=4.05$ using new IRAM Plateau de Bure interferometer (PdBI) CO(4-3) and 1.2-3.3 mm continuum observations. All these three SMGs are heavily dust obscured, with UV-based star formation rate (SFR) estimates significantly smaller than the ones derived from the bolometric infrared (IR), consistent with the spatial offsets revealed by HST and CO imaging. Based also on evaluations of the specific SFR, CO-to-H$_2$ conversion factor and gas depletion timescale, we classify all the three galaxies as starbursts (SBs), although with a lower confidence for GN20.2b that might be a later stage merging event. We place our measurements in the context of the evolutionary properties of main sequence (MS) and SB galaxies. ULIRGs have 3-5 times larger $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios than $z=0$ MS galaxies, but by $z\sim2$ the difference appears to be blurred, probably due to differential metallicity evolution. SB galaxies appear to slowly evolve in their $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios all the way to $z>6$ (consistent with rapid enrichment of SB events), while MS galaxies rapidly increase in $M_{\rm dust}/M_\star$ from $z=0$ to 2 (due to gas fraction increase, compensated by a decrease of metallicities). While no IR/submm continuum detection is available for indisputably normal massive galaxies at $z>2.5$, we show that if metallicity indeed decrease rapidly for these systems at $z>3$ as claimed in the literature, we should expect a strong decrease of their $M_{\rm dust}/M_\star$, consistent with recent PdBI and ALMA upper limits. We conclude that the $M_{\rm dust}/M_\star$ ratio could be a powerful tool for distinguishing starbursts from normal galaxies at $z>4$.

We investigate the relationships between stellar mass, gas-phase oxygen abundance (metallicity), star formation rate, and dust content of star-forming galaxies at z$\sim$1.6 using Subaru/FMOS spectroscopy in the COSMOS field. The mass-metallicity relation at $z\sim1.6$ is steeper than the relation observed in the local Universe. The steeper MZ relation at $z\sim1.6$ is mainly due to evolution in the stellar mass where the MZ relation begins to turnover and flatten. This turnover mass is 1.2 dex larger at $z\sim1.6$. The most massive galaxies at $z\sim1.6$ ($\sim 10^{11}M_\odot$) are enriched to the level observed in massive galaxies in the local Universe. The mass-metallicity relation we measure at $z\sim1.6$ supports the suggestion of an empirical upper metallicity limit that does not significantly evolve with redshift. We find an anti-correlation between metallicity and star formation rate for galaxies at a fixed stellar mass at $z\sim1.6$ which is similar to trends observed in the local Universe. We do not find a relation between stellar mass, metallicity and star formation rate that is independent of redshift; our data suggest that there is redshift evolution in this relation. We examine the relation between stellar mass, metallicity and dust extinction. We find that at a fixed stellar mass dustier galaxies tend to be more metal rich. From examination of the stellar masses, metallicities, SFRs and dust extinctions we conclude that stellar mass is most closely related to dust extinction.

We present detections of the CO(J= 1-0) emission line in a sample of four massive star-forming galaxies at z~1.5-2.2 obtained with the Karl G. Jansky Very Large Array (VLA). Combining these observations with previous CO(2-1) and CO(3-2) detections of these galaxies, we study the excitation properties of the molecular gas in our sample sources. We find an average line brightness temperature ratios of R_{21}=0.70+\-0.16 and R_{31}=0.50+\-0.29, based on measurements for three and two galaxies, respectively. These results provide additional support to previous indications of sub-thermal gas excitation for the CO(3-2) line with a typically assumed line ratio R_{31}~0.5. For one of our targets, BzK-21000, we present spatially resolved CO line maps. At the resolution of 0.18 arcsec (1.5 kpc), most of the emission is resolved out except for some clumpy structure. From this, we attempt to identify molecular gas clumps in the data cube, finding 4 possible candidates. We estimate that <40 % of the molecular gas is confined to giant clumps (~1.5 kpc in size), and thus most of the gas could be distributed in small fainter clouds or in fairly diffuse extended regions of lower brightness temperatures than our sensitivity limit.

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

Emission line diagnostic diagrams probing the ionization sources in galaxies, such as the Baldwin-Phillips-Terlevich (BPT) diagram, have been used extensively to distinguish AGN from purely star-forming galaxies. Yet, they remain poorly understood at higher redshifts. We shed light on this issue with an empirical approach based on a z~0 reference sample built from ~300,000 SDSS galaxies, from which we mimic selection effects due to typical emission line detection limits at higher redshift. We combine this low-redshift reference sample with a simple prescription for luminosity evolution of the global galaxy population to predict the loci of high-redshift galaxies on the BPT and Mass-Excitation (MEx) diagnostic diagrams. The predicted bivariate distributions agree remarkably well with direct observations of galaxies out to z~1.5, including the observed stellar mass-metallicity (MZ) relation evolution. As a result, we infer that high-redshift star-forming galaxies are consistent with having "normal" ISM properties out to z~1.5, after accounting for selection effects and line luminosity evolution. Namely, their optical line ratios and gas-phase metallicities are comparable to that of low-redshift galaxies with equivalent emission-line luminosities. In contrast, AGN narrow-line regions may show a shift toward lower metallicities at higher redshift. While a physical evolution of the ISM conditions is not ruled out for purely star-forming galaxies, and may be more important starting at z>2, we find that reliably quantifying this evolution is hindered by selections effects. The recipes provided here may serve as a basis for future studies toward this goal. Code to predict the loci of galaxies on the BPT and MEx diagnostic diagrams, and the MZ relation as a function of emission line luminosity limits, is made publicly available.

Using a phenomenological approach, we self-consistently model the redshift evolution of the ultraviolet (UV) and infrared (IR) luminosity functions across cosmic time, as well as a range of observed IR properties of UV-selected galaxy population. This model is an extension of the 2SFM (2 star-formation modes) formalism, which is based on the observed "main-sequence" of star-forming galaxies, i.e. a strong correlation between their stellar mass and their star formation rate (SFR), and a secondary population of starbursts with an excess of star formation. The balance between the UV light from young, massive stars and the dust-reprocessed IR emission is modeled following the empirical relation between the attenuation (IRX for IR excess hereafter) and the stellar mass, assuming a scatter of 0.4\,dex around this relation. We obtain a good overall agreement with the measurements of the IR luminosity function up to z~3 and the UV luminosity functions up to z~6, and show that a scatter on the IRX-M relation is mandatory to reproduce these observables. We also naturally reproduce the observed, flat relation between the mean IRX and the UV luminosity at L$_{\rm UV}>$10$^{9.5}$ L$_\odot$. Finally, we perform predictions of the UV properties and detectability of IR-selected samples and the vice versa, and discuss the results in the context of the UV-rest-frame and sub-millimeter surveys of the next decade.

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

We report the detection of 6.2um polycyclic aromatic hydrocarbon (PAH) and rest-frame 4-7um continuum emission in the z=4.055 submillimeter galaxy GN20, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. This represents the first detection of PAH emission at z>4. The strength of the PAH emission feature is consistent with a very high star formation rate of ~1600Msun/yr. We find that this intense starburst powers at least ~1/3 of the faint underlying 6um continuum emission, with an additional, significant (and perhaps dominant) contribution due to a power-law-like hot dust source, which we interpret to likely be a faint, dust-obscured active galactic nucleus (AGN). The inferred 6um AGN continuum luminosity is consistent with a sensitive upper limit on the hard X-ray emission as measured by the Chandra X-Ray Observatory if the previously undetected AGN is Compton-thick. This is in agreement with the finding at optical/infrared wavelengths that the galaxy and its nucleus are heavily dust-obscured. Despite the strong power-law component enhancing the mid-infrared continuum emission, the intense starburst associated with the photon-dominated regions that give rise to the PAH emission appears to dominate the total energy output in the infrared. GN20 is one of the most luminous starburst galaxies known at any redshift, embedded in a rich protocluster of star-forming galaxies. This investigation provides an improved understanding of the energy sources that power such exceptional systems, which represent the extreme end of massive galaxy formation at early cosmic times.

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 present the first results from a near-IR spectroscopic survey of the COSMOS field, using the Fiber Multi-Object Spectrograph on the Subaru telescope, designed to characterize the star-forming galaxy population at $1.4<z<1.7$. The high-resolution mode is implemented to detect H$\alpha$ in emission between $1.6{\rm -}1.8 \mathrm{\mu m}$ with $f_{\rm H\alpha}\gtrsim4\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$. Here, we specifically focus on 271 sBzK-selected galaxies that yield a H$\alpha$ detection thus providing a redshift and emission line luminosity to establish the relation between star formation rate and stellar mass. With further $J$-band spectroscopy for 89 of these, the level of dust extinction is assessed by measuring the Balmer decrement using co-added spectra. We find that the extinction ($0.6\lesssim A_\mathrm{H\alpha} \lesssim 2.5$) rises with stellar mass and is elevated at high masses compared to low-redshift galaxies. Using this subset of the spectroscopic sample, we further find that the differential extinction between stellar and nebular emission \hbox{$E_\mathrm{star}(B-V)/E_\mathrm{neb}(B-V)$} is 0.7--0.8, dissimilar to that typically seen at low redshift. After correcting for extinction, we derive an H$\alpha$-based main sequence with a slope ($0.81\pm0.04$) and normalization similar to previous studies at these redshifts.

Understanding the relationship between the formation and evolution of galaxies and their central super massive black holes (SMBH) is one of the main topics in extragalactic astrophysics. Links and feedback may reciprocally affect both black hole and galaxy growth. Observations of the CO line at redshifts of 2-4 are crucial to investigate the gas mass, star formation activity and accretion onto SMBHs, as well as the effect of AGN feedback. Potential correlations between AGN and host galaxy properties can be highlighted by observing extreme objects. Despite their luminosity, hyper-luminous QSOs at z=2-4 are still little studied at mm wavelengths. We targeted CO(3-2) in ULAS J1539+0557, an hyper-luminos QSO (Lbol> 10^48 erg/s) at z=2.658, selected through its unusual red colors in the UKIDSS Large Area Survey (ULAS). We find a molecular gas mass of 4.1+-0.8 10^10 Msun, and a gas fraction of 0.4-0.1, depending mostly on the assumed source inclination. We also find a robust lower limit to the star-formation rate (SFR=250-1600 Msun/yr) and star-formation efficiency (SFE=25-350 Lsun/(K km s-1 pc2) by comparing the observed optical-near-infrared spectral energy distribution with AGN and galaxy templates. The black hole gas consumption timescale, M(H_2)/dM(accretion)/dt, is ~160 Myr, similar or higher than the gas consumption timescale. The gas content and the star formation efficiency are similar to those of other high-luminosity, highly obscured QSOs, and at the lower end of the star-formation efficiency of unobscured QSOs, in line with predictions from AGN-galaxy co-evolutionary scenarios. Further measurements of the (sub)-mm continuum in this and similar sources are mandatory to obtain a robust observational picture of the AGN evolutionary sequence.

We report ALMA Cycle 0 observations at 1.3mm of LESS J033229.4-275619 (XID403), an Ultraluminous Infrared Galaxy at $z=4.75$ in the Chandra Deep Field South hosting a Compton-thick QSO. The source is not resolved in our data at a resolution of $\sim$0.75 arcsec, placing an upper-limit of 2.5 kpc to the half-light radius of the continuum emission from heated-dust. After deconvolving for the beam size, however, we found a $\sim3\sigma$ indication of an intrinsic source size of $0.27\pm0.08$ arcsec (Gaussian FWHM), which would correspond to $r_{half}\sim0.9\pm0.3$ kpc. We build the far-IR SED of XID403 by combining datapoints from both ALMA and Herschel and fit it with a modified blackbody spectrum. For the first time, we measure the dust temperature $T_d=58.5\pm5.3$ K in this system, which is comparable to what has been observed in other high-z submillimeter galaxies. The measured star formation rate is SFR=$1020\pm150$ $M_{\odot}$ yr$^{-1}$, in agreement with previous estimates at lower S/N. Based on the measured SFR and source size, we constrain the SFR surface density to be $\Sigma_{SFR}>26\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ ($\sim200\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ for $r_{half}\sim0.9$ kpc). The compactness of this starburst is comparable to what has been observed in other local and high-z starburst galaxies. If the gas mass measured from previous [CII] and CO(2-1) observations at low resolution is confined within the same dust region, assuming $r_{half}\sim0.9\pm0.3$ kpc, this would produce a column density of $N_H\sim0.3-1.1\times10^{24}$cm$^{-2}$ towards the central SMBH, similar to the column density of $\approx1.4\times10^{24}$cm$^{-2}$ measured from the X-rays. Then, in principle, if both gas and dust were confined on sub-kpc scales, this would be sufficient to produce the observed X-ray column density without any need of a pc-scale absorber [abridged].

We present direct constraints on the CO luminosity function at high redshift and the resulting cosmic evolution of the molecular gas density, $\rho_{\rm H2}$(z), based on a blind molecular line scan in the Hubble Deep Field North (HDF-N) using the IRAM Plateau de Bure Interferometer. Our line scan of the entire 3mm window (79-115 GHz) covers a cosmic volume of ~7000 Mpc$^3$, and redshift ranges z<0.45, 1.01<z<1.89 and z>2. We use the rich multiwavelength and spectroscopic database of the HDF-N to derive some of the best constraints on CO luminosities in high redshift galaxies to date. We combine the blind CO detections in our molecular line scan (presented in a companion paper) with stacked CO limits from galaxies with available spectroscopic redshifts (slit or mask spectroscopy from Keck and grism spectroscopy from HST) to give first blind constraints on high-z CO luminosity functions and the cosmic evolution of the H2 mass density $\rho_{\rm H2}$(z) out to redshifts z~3. A comparison to empirical predictions of $\rho_{\rm H2}$(z) shows that the securely detected sources in our molecular line scan already provide significant contributions to the predicted $\rho_{\rm H2}$(z) in the redshift bins <z>~1.5 and <z>~2.7. Accounting for galaxies with CO luminosities that are not probed by our observations results in cosmic molecular gas densities $\rho_{\rm H2}$(z) that are higher than current predictions. We note however that the current uncertainties (in particular the luminosity limits, number of detections, as well as cosmic volume probed) are significant, a situation that is about to change with the emerging ALMA observatory.

We present a molecular line scan in the Hubble Deep Field North (HDF-N) that covers the entire 3mm window (79-115 GHz) using the IRAM Plateau de Bure Interferometer. Our CO redshift coverage spans z<0.45, 1<z<1.9 and all z>2. We reach a CO detection limit that is deep enough to detect essentially all z>1 CO lines reported in the literature so far. We have developed and applied different line searching algorithms, resulting in the discovery of 17 line candidates. We estimate that the rate of false positive line detections is ~2/17. We identify optical/NIR counterparts from the deep ancillary database of the HDF-N for seven of these candidates and investigate their available SEDs. Two secure CO detections in our scan are identified with star-forming galaxies at z=1.784 and at z=2.047. These galaxies have colors consistent with the `BzK' color selection and they show relatively bright CO emission compared with galaxies of similar dust continuum luminosity. We also detect two spectral lines in the submillimeter galaxy HDF850.1 at z=5.183. We consider an additional 9 line candidates as high quality. Our observations also provide a deep 3mm continuum map (1-sigma noise level = 8.6 $\mu$Jy/beam). Via a stacking approach, we find that optical/MIR bright galaxies contribute only to <50% of the SFR density at 1<z<3, unless high dust temperatures are invoked. The present study represents a first, fundamental step towards an unbiased census of molecular gas in `normal' galaxies at high-z, a crucial goal of extragalactic astronomy in the ALMA era.

The relationships between galaxies of intermediate stellar mass and moderate luminosity active galactic nuclei (AGNs) at 1<z<3 are investigated with the Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS) sample complemented with public data in the GOODS-South field. Using X-ray data, hidden AGNs are identified in unsuspected star-forming galaxies with no apparent signs of non-stellar activity. In the color-mass plane, two parallel trends emerge during the ~2 Gyr between the average redshifts z~2.2 and z~1.3: while the red sequence becomes significantly more populated by ellipticals, the majority of AGNs with L(2-10 keV)>10^42.3 erg s^-1 disappear from the blue cloud/green valley where they were hosted predominantly by star-forming systems with disk and irregular morphologies. These results are even clearer when the rest-frame colors are corrected for dust reddening. At z~2.2, the ultraviolet spectra of active galaxies (including two Type 1 AGNs) show possible gas outflows with velocities up to about -500 km s^-1 that are not observed neither in inactive systems at the same redshift, nor at lower redshifts. Such outflows indicate the presence of gas that can move faster than the escape velocities of active galaxies. These results suggest that feedback from moderately luminous AGNs (logL_X<44.5 erg s^-1) played a key role at z>~2 by contributing to outflows capable of ejecting part of the interstellar medium and leading to a rapid decrease in the star formation in host galaxies with stellar masses 10<logM<11 M_Sun.

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

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

We present an analysis of the clustering of galaxies from z ~ 2 to the present day using the WIRCam Deep Survey (WIRDS). WIRDS combines deep near-infrared data with the deep optical data from the CFHTLS Deep fields, providing a photometric data-set over an effective area of 2.4 sq. deg., from which accurate photometric redshifts and stellar masses can be estimated. We use the data to calculate the angular correlation function for galaxy samples split by star-formation activity, stellar mass and redshift. We estimate the real-space clustering for each sample, determining clustering lengths and power-law slopes. For galaxies selected by constant mass, we find that the clustering scale shows no evolution up to z ~ 2. Splitting the galaxy sample by mass, we see that higher mass galaxies have larger clustering scales at all redshifts. We use our results to test the GALFORM semi-analytical galaxy formation model and find the two are consistent. We split the galaxy population into passive and star-forming populations and find that the passive galaxy population shows a significantly larger clustering scale at all redshifts than the star-forming population below masses of ~$10^{11}M_\odot/h$, showing that even at z ~ 2 passive galaxies exist in denser environments than the bulk of the star-forming galaxy population. For star-forming galaxies with stellar masses $>10^{11}M_\odot/h$, we find a clustering strength of ~8Mpc/h across all redshifts, comparable to the measurements for the passive population. Also, for star-forming galaxies we see that clustering strength increases for higher stellar mass systems, however there is little sign of a mass dependence in passive galaxies. Finally, we investigate the connection between galaxy stellar mass and dark matter halo mass, showing a clear correlation between the two in both the WIRDS data and the GALFORM predictions.

We present deep Hubble Space Telescope Wide Field Camera 3 slitless spectroscopic observations of the distant cluster Cl J1449+0856. These cover a single pointing with 18 orbits of G141 spectroscopy and F140W imaging, allowing us to derive secure redshifts down to m_140~25.5 AB and 3sigma line fluxes of 5*10^(-18) erg/s/cm^2. In particular, we were able to spectroscopically confirm 12 early-type galaxies in the field up to z~3, 6 of which in the cluster core, which represents the first direct spectroscopic confirmation of passive galaxies in a z=2 cluster environment. With 140 redshifts in a ~6 arcmin^2 field, we can trace the spatial and redshift galaxy distribution in the cluster core and background field. We find two strong peaks at z=2.00 and z=2.07, where only one was seen in our previously published ground-based data. Thanks to the spectroscopic confirmation of the cluster ETGs, we can now re-evaluate the redshift of Cl J1449+0856 at z=2.00, rather than z=2.07, with the background overdensity being revealed to be sparse and "sheet"-like. This presents an interesting case of chance alignment of two close yet unrelated structures, each one preferentially selected by different observing strategies. With 6 quiescent or early-type spectroscopic members and 20 star-forming ones, Cl J1449+0856 is now reliably confirmed to be at z=2.00. The identified members can now allow for a detailed study of galaxy properties in the densest environment at z=2.

We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA.

Star-forming disk galaxies at high redshift are often subject to violent disk instability, characterized by giant clumps whose fate is yet to be understood. The main question is whether the clumps disrupt within their dynamical timescale (<50Myr), like molecular clouds in today's galaxies, or whether they survive stellar feedback for more than a disk orbital time (~300Myr) in which case they can migrate inward and help building the central bulge. We present 3.5-7pc resolution AMR simulations of high-redshift disks including photo-ionization, radiation pressure, and supernovae feedback (Renaud et al. 2013, and Perret et al., this astro-ph issue). Our modeling of radiation pressure determines the mass loading and initial velocity of winds from basic physical principles. We find that the giant clumps produce steady outflow rates comparable to and sometimes somewhat larger than their star formation rate, with velocities largely sufficient to escape galaxy. The clumps also lose mass, especially old stars, by tidal stripping, and the stellar populations contained in the clumps hence remain relatively young (<=200Myr), as observed. The clumps survive gaseous outflows and stellar loss, because they are wandering in gas-rich turbulent disks from which they can re-accrete gas at high rates compensating for outflows and tidal stripping, overall keeping realistic and self-regulated gaseous and stellar masses. Our simulations produce gaseous outflows with velocities, densities and mass loading consistent with observations, and at the same time suggest that the giant clumps survive for hundreds of Myr and complete their migration to the center of high-redshift galaxies, without rapid dispersion and reformation of clumps.

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.

We present a study of the formation of clustered, massive galaxies at large look-back times via spectroscopic imaging of CO in the unique GN20 proto-cluster at z = 4.05. Existing observations show that this is a dense concentration of gas-rich, very active star forming galaxies, including multiple bright submillimeter galaxies (SMGs). Using deep, high-resolution VLA CO(2-1) observations, we image the molecular gas with a resolution of ~1 kpc just 1.6 Gyr after the Big Bang. The SMGs GN20.2a and GN20.2b have deconvolved sizes of ~5 kpc X 3 kpc and ~8 kpc X 5 kpc (Gaussian FWHM) in CO(2-1), respectively, and we measure gas surface densities up to ~12,700/1,700X(sin i) (\alpha_CO/0.8) M_sun/pc^2 for GN20.2a/GN20.2b in the highest-resolution maps. Dynamical mass estimates allow us to constrain the CO-to-H_2 conversion factor to \alpha_CO = 1.7+/-0.8 M_sun (K km s^{-1} pc^2)^-1 for GN20.2a and \alpha_CO = 1.1+/-^{1.5}_{1.1} M_sun (K km/s pc^2)^-1 for GN20.2b. We measure significant offsets (0.5"-1") between the CO and optical emission, indicating either dust obscuration on scales of tens of kpc or that the emission originates from distinct galaxies. CO spectral line energy distributions imply physical conditions comparable to other SMGs and reveal further evidence that GN20.2a and GN20.2b are in different merging stages. We carry out a targeted search for CO emission from the 14 known B-band Lyman break galaxies (LBGs) in the field, tentatively detecting CO in a previously-undetected LBG and placing 3sigma upper limits on the CO luminosities of those that may lie within our bandpass. A blind search for emission line sources down to a 5sigma limiting CO luminosity of L'_CO(2-1) = 8 X 10^9 K km/s pc^2 and covering \Delta z = 0.0273 (~20 comoving Mpc) produces no other strong contenders associated with the proto-cluster.

[Abridged] Recent studies revealed a strong correlation between the star formation rate (SFR) and stellar mass of star-forming galaxies, the so-called star-forming main sequence. An empirical modeling approach (2-SFM) which distinguishes between the main sequence and rarer starburst galaxies is capable of reproducing most statistical properties of infrared galaxies. In this paper, we extend this approach by establishing a connection between stellar mass and halo mass with the technique of abundance matching. Based on a few, simple assumptions and a physically motivated formalism, our model successfully predicts the (cross-)power spectra of the cosmic infrared background (CIB), the cross-correlation between CIB and cosmic microwave background (CMB) lensing, and the correlation functions of bright, resolved infrared galaxies measured by Herschel, Planck, ACT and SPT. We use this model to infer the redshift distribution these observables, as well as the level of correlation between CIB-anisotropies at different wavelengths. We also predict that more than 90% of cosmic star formation activity occurs in halos with masses between 10^11.5 and 10^13.5 Msun. Taking into account subsequent mass growth of halos, this implies that the majority of stars were initially (at z>3) formed in the progenitors of clusters, then in groups at 0.5<z<3 and finally in Milky-Way-like halos at z<0.5. At all redshifts, the dominant contribution to the star formation rate density stems from halos of mass ~10^12 Msun, in which the instantaneous star formation efficiency is maximal (~70%). The strong redshift-evolution of SFR in the galaxies that dominate the CIB is thus plausibly driven by increased accretion from the cosmic web onto halos of this characteristic mass scale. Material available at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

We explore the relationship between gas, dust and star formation in a sample of 12 ultra-luminous infrared galaxies (ULIRGs) at high redshift compared to a similar sample of local galaxies. We present new CO observations and/or Spitzer mid-IR spectroscopy for 6 70 micron selected galaxies at z~1 in order to quantify the properties of the molecular gas reservoir, the contribution of an active galactic nuclei (AGN) to the mid-IR luminosity and the star formation efficiency (SFE=LIR/L'CO). The mid-IR spectra show strong polycyclic aromatic hydrocarbon (PAH) emission and our spectral decomposition suggests that the AGN makes a minimal contribution (<25%) to the mid-IR luminosity. The 70 micron selected ULIRGs which we find to be spectroscopic close pairs, are observed to have high SFE, similar to local ULIRGs and high redshift submillimeter galaxies, consistent with enhanced IR luminosity due to an ongoing major merger. Combined with existing observations of local and high redshift ULIRGs, we further compare the PAH, IR and CO luminosities. We show that the ratio LPAH6.2/LIR decreases with increasing IR luminosity for both local and high redshift galaxies but the trend for high redshift galaxies is shifted to higher IR luminosities; the average LPAH6.2/LIR ratio at a given LIR is ~3 times higher at high redshift. When we normalize by the molecular gas, we find this trend to be uniform for galaxies at all redshifts and that the molecular gas is correlated with the PAH dust emission.The similar trends seen in the [CII] to molecular gas ratios in other studies suggests that PAH emission, like [CII], continues to be a good tracer of photodissociation regions even at high redshift. Together the CO, PAH and far-IR fine structure lines should be useful for constraining the interstellar medium conditions in high redshift galaxies.

We report on results from the analysis of a stellar mass-selected (log M*>9.0) sample of 1644 galaxies at 0.65<z<1.1 with ultra-deep (m<26.5) optical medium-band (R~50) photometry from the Survey for High-z Absorption Red and Dead Sources (SHARDS). The spectral resolution of SHARDS allows us to consistently measure the strength of the 4000 Angstrom spectral break [Dn(4000), an excellent age indicator for the stellar populations of quiescent galaxies] for all galaxies at z~0.9 down to log M*9. The Dn(4000) index cannot be resolved from broad-band photometry, and measurements from optical spectroscopic surveys are typically limited to galaxies at least x10 more massive. When combined with the rest-frame U-V colour, Dn(4000) provides a powerful diagnostic of the extinction affecting the stellar population that is relatively insensitive to degeneracies with age, metallicity or star formation history. We use this novel approach to estimate de-reddened colours and light-weighted stellar ages for individual sources. We explore the relationships linking stellar mass, (U-V), and Dn(4000) for the sources in the sample, and compare them to those found in local galaxies. The main results are: a) both Dn(4000) and (U-V) correlate with M*. The dispersion in Dn(4000) values at a given M* increases with M*, while the dispersion for (U-V) decreases due to the higher average extinction prevalent in massive star-forming galaxies. b) for massive galaxies, we find a smooth transition between the blue cloud and red sequence in the intrinsic U-V colour, in contrast with other recent results. c) at a fixed stellar age, we find a positive correlation between extinction and stellar mass. d) the fraction of sources with declining or halted star formation increases steeply with the stellar mass, from ~5% at log M*~9.0-9.5 to ~80% at log M*>11, in agreement with downsizing scenarios.

Current X-ray surveys have proved to be essential tools in order to identify and study AGNs across cosmic time. However, there is evidence that the most heavily obscured AGNs are largely missing even in the deepest surveys. The search for these obscured AGNs is one of the most outstanding issues of extragalactic astronomy, since they are expected to make a major contribution to the high energy peak of the X-ray background (XRB) and might constitute a particularly active and dusty phase of black hole and galaxy evolution. Using a newly developed SED fitting technique to decompose the AGN from the star-formation emission in the infrared band we identify a sample of 17 IR bright quasars at z=1-3. For the majority of these sources the X-ray spectrum is well characterised by an absorbed power law model, revealing that ~25% of the sources are unabsorbed; the remainder of these IR quasars have moderate-to-high absorption (N_H>10^23 cm^-2), and ~25% are not detected in the X-rays and are likely to be Compton thick (N_H>10^24 cm^-2). We therefore find a much higher fraction of obscured than unobscured quasars, indicating that a large fraction of the luminous black hole accretion was very heavily obscured at z~2.

We present a study of galaxy populations in the central region of the IRAC-selected, X-ray detected galaxy cluster Cl J1449+0856 at z=2. Based on a sample of spectroscopic and photometric cluster members, we investigate stellar populations and morphological structure of cluster galaxies over an area of ~0.7Mpc^2 around the cluster core. The cluster stands out as a clear overdensity both in redshift space, and in the spatial distribution of galaxies close to the center of the extended X-ray emission. The cluster core region (r<200 kpc) shows a clearly enhanced passive fraction with respect to field levels. However, together with a population of massive passive galaxies mostly with early-type morphologies, it also hosts massive actively star-forming, often highly dust-reddened sources. Close to the cluster center, a multi-component system of passive and star-forming galaxies could be the future BCG still assembling. We observe a clear correlation between passive stellar populations and an early-type morphology, in agreement with field studies at similar redshift. Passive early-type galaxies in this clusters are typically a factor 2-3 smaller than similarly massive early-types at z~0, but also on average larger by a factor ~2 than their field analogs at z~2, lending support to recent claims of an accelerated structural evolution in high-redshift dense environments. These results point towards the early formation of a population of massive galaxies, already evolved both in their structure and stellar populations, coexisting with still-actively forming massive galaxies in the central regions of young clusters 10 billion years ago.

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.

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

Modern (sub-)millimeter interferometers enable the measurement of the cool gas and dust emission of high-redshift galaxies (z>5). However, at these redshifts the cosmic microwave background (CMB) temperature is higher, approaching, and even exceeding, the temperature of cold dust and molecular gas observed in the local Universe. In this paper, we discuss the impact of the warmer CMB on (sub-)millimeter observations of high-redshift galaxies. The CMB affects the observed (sub-)millimeter dust continuum and the line emission (e.g. carbon monoxide, CO) in two ways: (i) it provides an additional source of (both dust and gas) heating; and (ii) it is a non-negligible background against which the line and continuum emission are measured. We show that these two competing processes affect the way we interpret the dust and gas properties of high-redshift galaxies using spectral energy distribution models. We quantify these effects and provide correction factors to compute what fraction of the intrinsic dust (and line) emission can be detected against the CMB as a function of frequency, redshift and temperature. We discuss implications on the derived properties of high-redshift galaxies from (sub-)millimeter data. Specifically, the inferred dust and molecular gas masses can be severely underestimated for cold systems if the impact of the CMB is not properly taken into account.

We characterize the incidence of active galactic nuclei (AGNs) is 0.3 < z < 1 star-forming galaxies by applying multi-wavelength AGN diagnostics (X-ray, optical, mid-infrared, radio) to a sample of galaxies selected at 70-micron from the Far-Infrared Deep Extragalactic Legacy survey (FIDEL). Given the depth of FIDEL, we detect "normal" galaxies on the specific star formation rate (sSFR) sequence as well as starbursting systems with elevated sSFR. We find an overall high occurrence of AGN of 37+/-3%, more than twice as high as in previous studies of galaxies with comparable infrared luminosities and redshifts but in good agreement with the AGN fraction of nearby (0.05 < z < 0.1) galaxies of similar infrared luminosities. The more complete census of AGNs comes from using the recently developed Mass-Excitation (MEx) diagnostic diagram. This optical diagnostic is also sensitive to X-ray weak AGNs and X-ray absorbed AGNs, and reveals that absorbed active nuclei reside almost exclusively in infrared-luminous hosts. The fraction of galaxies hosting an AGN appears to be independent of sSFR and remains elevated both on the sSFR sequence and above. In contrast, the fraction of AGNs that are X-ray absorbed increases substantially with increasing sSFR, possibly due to an increased gas fraction and/or gas density in the host galaxies.

[abridged] Modern (sub-)millimeter/radio interferometers will enable us to measure the dust and molecular gas emission from galaxies that have luminosities lower than the Milky Way, out to high redshifts and with unprecedented spatial resolution and sensitivity. This will provide new constraints on the star formation properties and gas reservoir in galaxies throughout cosmic times through dedicated deep field campaigns targeting the CO/[CII] lines and dust continuum emission. In this paper, we present empirical predictions for such (sub-)millimeter line and continuum deep fields. We base these predictions on the deepest available optical/near-infrared ACS and NICMOS data on the Hubble Ultra Deep Field. Using a physically-motivated spectral energy distribution model, we fit the observed optical/near-infrared emission of 13,099 galaxies with redshifts up to z=5, and obtain median likelihood estimates of their stellar mass, star formation rate, dust attenuation and dust luminosity. We derive statistical constraints on the dust emission in the infrared and (sub-)millimeter which are consistent with the observed optical/near-infrared emission in terms of energy balance. This allows us to estimate, for each galaxy, the (sub-)millimeter continuum flux densities in several ALMA, PdBI/NOEMA and JVLA bands. Using empirical relations between the observed CO/[CII] line luminosities and the infrared luminosity, we infer the flux of the CO(1-0) and [CII] lines from the estimated infrared luminosity of each galaxy in our sample. We then predict the fluxes of higher CO transition lines CO(2-1) to CO(7-6) bracketing two extreme gas excitation scenarios. We use our predictions to discuss possible deep field strategies with ALMA. The predictions presented in this study will serve as a direct benchmark for future deep field campaigns in the (sub-)millimeter regime.

We have compiled a large sample of 151 high redshift (z=0.5-4) galaxies selected at 24 microns (S24>100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-infrared spectrum into contributions from star formation and activity in the galactic nuclei. In addition, we have a wealth of photometric data from Spitzer IRAC/MIPS and Herschel PACS/SPIRE. We explore how effective different infrared color combinations are at separating our mid-IR spectroscopically determined active galactic nuclei from our star forming galaxies. We look in depth at existing IRAC color diagnostics, and we explore new color-color diagnostics combining mid-IR, far-IR, and near-IR photometry, since these combinations provide the most detail about the shape of a source's IR spectrum. An added benefit of using a color that combines far-IR and mid-IR photometry is that it is indicative of the power source driving the IR luminosity. For our data set, the optimal color selections are S250/S24 vs. S8.0/S3.6 and S100/S24 vs. S8.0/S3.6; both diagnostics have ~10% contamination rate in the regions occupied primarily by star forming galaxies and active galactic nuclei, respectively. Based on the low contamination rate, these two new IR color-color diagnostics are ideal for estimating both the mid-IR power source of a galaxy when spectroscopy is unavailable and the dominant power source contributing to the IR luminosity. In the absence of far-IR data, we present color diagnostics using the WISE mid-IR bands which can efficiently select out high z (z~2) star forming galaxies.

Many theoretical models require powerful active galactic nuclei (AGNs) to suppress star formation in distant galaxies and reproduce the observed properties of today's massive galaxies. A recent study based on Herschel-SPIRE submillimetre observations claimed to provide direct support for this picture, reporting a significant decrease in the mean star-formation rates (SFRs) of the most luminous AGNs (Lx>10^44 erg/s) at z=1-3 in the Chandra Deep Field-North (CDF-N). In this letter we extend these results using Herschel-SPIRE 250um data in the COSMOS and CDF-S fields to achieve an order of magnitude improvement in the number of sources at Lx>10^44 erg/s. On the basis of our analysis, we find no strong evidence for suppressed star formation in Lx>10^44 erg/s AGNs at z=1-3. The mean SFRs of the AGNs are constant over the broad X-ray luminosity range of Lx~10^43-10^45 erg/s (with mean SFRs consistent with typical star-forming galaxies at z~2; <SFRs>~100-200 Msol/yr). We suggest that the previous CDF-N results were likely due to low number statistics. We discuss our results in the context of current theoretical models.

We present the discovery of a massive, quiescent galaxy at z=2.99. We have obtained a HST/WFC3 spectrum of this object and measured its redshift from the detection of a deep 4000A break consistent with an old population and a high metallicity. By stellar population modeling of both its grism spectrum and broad-band photometry, we derive an age of ~0.7 Gyr, implying a formation redshift of z>4, and a mass >10^11 Msun. Although this passive galaxy is the most distant confirmed so far, we find that it is slightly less compact than other z>2 early-types of similar mass, being overall more analogous to those z~1.6 field early-type galaxies. The discovery of this object shows that early-type galaxies are detectable to at least z=3 and suggests that the diversity of structural properties found in z=1.4-2 ellipticals to earlier epochs could have its origin in a variety of formation histories among their progenitors.

(Abridged) We present the Survey for High-z Absorption Red and Dead Sources (SHARDS), an ESO/GTC Large Program carried out with GTC/OSIRIS. SHARDS is an ultra-deep optical spectro-photometric survey of the GOODS-N field (130 arcmin^2) at wavelengths 500 to 950 nm and using 24 contiguous medium-band filters (spectral resolution R 50). The data reach 26.5 mag (>3-sigma level) with sub-arcsec seeing in all bands. SHARDS main goal is obtaining accurate physical properties of interm- and high-z galaxies using well-sampled optical SEDs with sufficient spectral resolution to measure absorption and emission features. Among the different populations of high-z galaxies, SHARDS principal targets are massive quiescent galaxies at z>1. In this paper, we outline the observational strategy and include a detailed discussion of the special reduction and calibration procedures applied to the GTC/OSIRIS data. We present science demonstration results about the detection and study of emission-line galaxies (star-forming and AGN) at z=0-5. We also analyze the SEDs for a sample of 27 quiescent massive galaxies at 1.0<z<1.4. We discuss on the improvements introduced by the SHARDS dataset in the analysis of their SFH and stellar properties. We discuss the systematics arising from the use of different stellar population libraries. We find that the UV-to-MIR SEDs of the massive quiescent galaxies at z=1.0-1.5 are well described by an exponential decaying SFH with scale tau=100-200 Myr, age 1.5-2.0 Gyr, solar or slightly sub-solar metallicity, and moderate extinction, A(V)~0.5 mag. We also find that galaxies with masses above M* are typically older than lighter galaxies, as expected in a downsizing scenario of galaxy formation. This trend is, however, model dependent, i.e., it is significantly more evident in the results obtained with some stellar population synthesis libraries and almost absent in others.

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

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

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content.

We present Karl G. Jansky Very Large Array (VLA) observations of the CO(2-1) emission in the z=4.05 submillimeter galaxy (SMG) GN20. These high-resolution data allow us to image the molecular gas at 1.3 kpc resolution just 1.6 Gyr after the Big Bang. The data reveal a clumpy, extended gas reservoir, 14 +/- 4 kpc in diameter, in unprecedented detail. A dynamical analysis shows that the data are consistent with a rotating disk of total dynamical mass 5.4 +/- 2.4 X 10^11 M_sun. We use this dynamical mass estimate to constrain the CO-to-H_2 mass conversion factor (alpha_CO), finding alpha_CO=1.1 +/- 0.6 M_sun (K km s^-1 pc^2)^-1. We identify five distinct molecular gas clumps in the disk of GN20 with masses a few percent of the total gas mass, brightness temperatures of 16-31K, and surface densities of >3,200-4,500 X (alpha_CO/0.8) M_sun pc^-2. Virial mass estimates indicate they could be self-gravitating, and we constrain their CO-to-H_2 mass conversion factor to be <0.2-0.7 M_sun (K km s^-1 pc^2)^-1. A multiwavelength comparison demonstrates that the molecular gas is concentrated in a region of the galaxy that is heavily obscured in the rest-frame UV/optical. We investigate the spatially-resolved gas excitation and find that the CO(6-5)/CO(2-1) ratio is constant with radius, consistent with star formation occuring over a large portion of the disk. We discuss the implications of our results in the context of different fueling scenarios for SMGs.

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

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

We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 {\mu}Jy. Material of the model (predictions, SED library, mock catalogs...) is available online at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

MOONS is a new conceptual design for a Multi-Object Optical and Near-infrared Spectrograph for the Very Large Telescope (VLT), selected by ESO for a Phase A study. The baseline design consists of 1000 fibers deployable over a field of view of 500 square arcmin, the largest patrol field offered by the Nasmyth focus at the VLT. The total wavelength coverage is 0.8um-1.8um and two resolution modes: medium resolution and high resolution. In the medium resolution mode (R=4,000-6,000) the entire wavelength range 0.8um-1.8um is observed simultaneously, while the high resolution mode covers simultaneously three selected spectral regions: one around the CaII triplet (at R=8,000) to measure radial velocities, and two regions at R=20,000 one in the J-band and one in the H-band, for detailed measurements of chemical abundances. The grasp of the 8.2m Very Large Telescope (VLT) combined with the large multiplex and wavelength coverage of MOONS - extending into the near-IR - will provide the observational power necessary to study galaxy formation and evolution over the entire history of the Universe, from our Milky Way, through the redshift desert and up to the epoch of re-ionization at z>8-9. At the same time, the high spectral resolution mode will allow astronomers to study chemical abundances of stars in our Galaxy, in particular in the highly obscured regions of the Bulge, and provide the necessary follow-up of the Gaia mission. Such characteristics and versatility make MOONS the long-awaited workhorse near-IR MOS for the VLT, which will perfectly complement optical spectroscopy performed by FLAMES and VIMOS.

Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)

We present a new near-infrared imaging survey in the four CFHTLS deep fields: the WIRCam Deep Survey (WIRDS). WIRDS comprises extremely deep, high quality (FWHM ~0.6") J, H and K imaging covering a total effective area of 2.1 deg^2 and reaching AB 50% completeness limits of ~24.5. We combine our images with the CFHTLS to create a unique eight-band ugrizJHK photometric catalogues in the CFHTLS deep fields; these four separate fields allow us to make a robust estimate of the effect of cosmic variance for all our measurements. We use these catalogues to estimate precise photometric redshifts, galaxy types and stellar masses for a unique sample of ~1.8 million galaxies. Our JHK number counts are consistent with previous studies. We apply the BzK selection to our gzK filter set and find that the star forming BzK selection successfully selects 76% of star-forming galaxies in the redshift range 1.4<z<2.5 in our photometric catalogue. The passive BzK selection returns 52% of the passive 1.4<z<2.5 population identified in the photometric catalogue. We present the galaxy stellar mass function as a function of redshift up to z=2 and present fits using double Schechter functions. A mass-dependent evolution of the mass function is seen with the numbers of galaxies with masses of log(M)<10.75 still evolving at z<1, but galaxies of higher mass reaching their present day numbers by z~0.8-1. This is consistent with the present picture of downsizing in galaxy evolution. We compare our results with the predictions of the GALFORM semi-analytical galaxy formation model and find that the simulations provide a relatively successful fit to the observed mass functions at intermediate masses (i.e. 10<log(M)<11). However, the GALFORM results under-predict the mass function at low masses, whilst the fit as a whole degrades beyond redshifts of z~1.2.

We provide evidence for a correlation between the presence of giant clumps and the occurrence of active galactic nuclei (AGN) in disk galaxies. Giant clumps of 10^8-9 Msun arise from violent gravitational instability in gas-rich galaxies, and it has been proposed that this instability could feed supermassive black holes (BH). We use emission line diagnostics to compare a sample of 14 clumpy (unstable) disks and a sample of 13 smoother (stable) disks at redshift z~0.7. The majority of clumpy disks in our sample have a high probability of containing AGN. Their [OIII] emission line is strongly excited, inconsistent with low-metallicity star formation alone. [NeIII] excitation is also higher. Stable disks rarely have such properties. Stacking ultra sensitive Chandra observations (4 Ms) reveals an X-ray excess in clumpy galaxies, which cannot be solely due to star formation and confirms the presence of AGN. The clumpy galaxies in our intermediate-redshift sample have properties typical of gas-rich disk galaxies rather than mergers, being in particular on the Main Sequence of star formation. This suggests that our findings apply to the physically-similar and numerous gas-rich unstable disks at z>1. Using the observed [OIII] and X-ray luminosities, we conservatively estimate that AGN hosted by clumpy disks have typical bolometric luminosities of the order of a few 10^43 erg/s, BH growth rates ~10^-2 Msun/yr, and that these AGN are substantially obscured in X-rays. This moderate-luminosity mode could be sufficient to provide a large fraction of today's BH mass over a couple of Gyr given that our observations suggest a high duty cycle (>10%), accretion bursts with higher luminosities being possible over shorter phases. The observed evolution of disk instabilities with mass and redshift could explain the simultaneous downsizing of star formation and of BH growth.

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive.

Using a mass-limited sample of 24um-detected, star-forming galaxies at 0.5<z<1.3, we study the mass-star formation rate (SFR) correlation and its tightness. The correlation is well defined (sigma=0.28dex) for disk galaxies (n_sersic<1.5), while more bulge-dominated objects often have lower specific SFRs. For disk galaxies, a much tighter correlation (sigma=0.19dex) is obtained if the rest-frame H-band luminosity is used instead of stellar mass derived from multicolor photometry. The specific SFR (sSFR) correlates strongly with rest-frame optical colors (hence luminosity-weighted stellar age) and also with clumpiness (which likely reflects the molecular gas fraction). This implies that most of the observed scatter is real, despite its low level, and not dominated by random measurement errors. After correcting for these differential effects a remarkably small dispersion remains (sigma=0.14dex), suggesting that measurement errors in mass or SFR are ~0.10dex, excluding systematic uncertainties. Measurement errors in stellar masses, the thickening of the correlation due to real sSFR variations, and varying completeness with stellar mass, can spuriously bias the derived slope to lower values due to the finite range over which observables (mass and SFR) are available. When accounting for these effects, the intrinsic slope for the main sequence for disk galaxies gets closer to unity.

[Abridged] We present the results of new near-IR spectroscopic observations of passive galaxies at z>1.4 in a concentration of BzK-selected galaxies in the COSMOS field. The observations have been conducted with Subaru/MOIRCS, and have resulted in absorption lines and/or continuum detection for 18 out of 34 objects. This allows us to measure spectroscopic redshifts for a sample almost complete to K(AB)=21. COSMOS photometric redshifts are found in fair agreement overall with the spectroscopic redshifts, with a standard deviation of ~0.05; however, ~30% of objects have photometric redshifts systematically underestimated by up to ~25%. We show that these systematic offsets in photometric redshifts can be removed by using these objects as a training set. All galaxies fall in four distinct redshift spikes at z=1.43, 1.53, 1.67 and 1.82, with this latter one including 7 galaxies. SED fits to broad-band fluxes indicate stellar masses in the range of ~4-40x10^10Msun and that star formation was quenched ~1 Gyr before the cosmic epoch at which they are observed. The spectra of several individual galaxies have allowed us to measure their Hdelta_F and Dn4000 indices, which confirms their identification as passive galaxies, as does a composite spectrum resulting from the coaddition of 17 individual spectra. The effective radii of the galaxies have been measured on the HST/ACS F814W image, confirming the coexistence at these redshifts of passive galaxies which are substantially more compact than their local counterparts with others that follow the local size-stellar mass relation. For the galaxy with best S/N spectrum we were able to measure a velocity dispersion of 270+/-105 km/s, indicating that this galaxy lies closely on the virial relation given its stellar mass and effective radius.

Using X-ray stacking analyses we estimate the average amounts of supermassive black hole (SMBH) growth taking place in star-forming galaxies (SFGs) at z~1 and z~2 as a function of galaxy stellar mass (M*). We find the average SMBH growth rate follows remarkably similar trends with M* and redshift as the average star-formation rates (SFRs) of their host galaxies (i.e., (i.e., dM_BH/dt ~ M*^(0.86+/-0.39) for the z~1 sample and dM_BH/dt ~ M*^(1.05+/-0.36) for the z~2 sample). It follows that the ratio of SMBH growth rate to SFR is (a) flat with respect to M* (b) not evolving with redshift and (c) close to the ratio required to maintain/establish a SMBH to M* ratio of ~10^{-3} as also inferred from today's M_BH-M_bulge relationship. We interpret this as evidence that SMBHs have, on average, grown in-step with their host galaxies since at least z~2, irrespective of host galaxy mass and AGN triggering mechanism. As such, we suggest that the same secular processes that drive the bulk of star formation are also responsible for the majority of SMBH growth. From this, we speculate that it is the availability of gas reservoirs that regulate both cosmological SMBH growth and star formation.

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 present a parameter-less approach to predict the shape of the infrared (IR) luminosity function (LF) at redshifts z < 2. It requires no tuning and relies on only three observables: (1) the redshift evolution of the stellar mass function for star-forming galaxies, (2) the evolution of the specific star formation rate (sSFR) of main-sequence galaxies, and (3) the double-Gaussian decomposition of the sSFR-distribution at fixed stellar mass into a contribution (assumed redshift- and mass-invariant) from main-sequence and starburst activity. This self-consistent and simple framework provides a powerful tool for predicting cosmological observables: observed IR LFs are successfully matched at all z < 2, suggesting a constant or only weakly redshift-dependent contribution (8-14%) of starbursts to the star formation rate density. We separate the contributions of main-sequence and starburst activity to the global IR LF at all redshifts. The luminosity threshold above which the starburst component dominates the IR LF rises from log(LIR/Lsun) = 11.4 to 12.8 over 0 < z < 2, reflecting our assumed (1+z)^2.8-evolution of sSFR in main-sequence galaxies.

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

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

We present a study of the formation of clustered, massive galaxies at large look--back times, via high resolution spectroscopic imaging of CO in the unique GN20 proto-cluster. The data reveal evidence for rich structure and gas dynamics in unprecedented detail, allowing us to image the molecular gas with a resolution of only 1 kpc just 1.5 Gyr after the Big Bang. These state-of-the-art data give new insight into the detailed physical processes involved in early massive galaxy formation, and they provide a first glimpse of the morphological studies that will become feasible on a regular basis with ALMA.

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

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

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

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.

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.

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to continuous feeding by intense cold streams leading to violent gravitational instability, transient features and giant clumps. Gravitational torques between the perturbations drive angular momentum out and mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates of the inflow for a self-regulated unstable disk at a Toomre stability parameter Q~1, and isolated galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an average inflow rate ~10 Msun/yr through the disk of a 10^11 Msun galaxy, with conditions representative of z~2 stream-fed disks. The inflow rate scales with disk mass and (1+z)^{3/2}. It includes clump migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows the bulge, while only a fraction ~ 10^-3 of it needs to accrete onto a central black hole (BH), in order to obey the observed BH-bulge relation. A galaxy of 10^11 Msun at z~2 is expected to host a BH of ~10^8 Msun, accreting on average with moderate sub-Eddington luminosity L_X ~ 10^42-43 erg/s, accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies at z~6, the same process may feed 10^9 Msun BH at the Eddington rate. High central gas column densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep X-ray surveys.

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

We study the impact of galaxy-galaxy interactions on the FIR properties of galaxies and its evolution at 0<z<1.2. Using the high-z galaxies in the fields of GOODS observed by Herschel in the framework of the GOODS-Herschel key program and the local IRAS or AKARI-selected galaxies in the field of SDSS DR7, we investigate the dependence of galaxy properties on the morphology of and the distance to the nearest neighbor galaxy. We find that the SFRs and the SSFRs of galaxies, on average, depend on the morphology of and the distance to the nearest neighbor galaxy in this redshift range. When a late-type galaxy has a close neighbor galaxy, the SFR and the SSFR increase as it approaches a late-type neighbor, which is supported by K-S and MC tests with a significance level of >99%. However, the SFR and the SSFR decrease or do not change much as it approaches an early-type neighbor. The bifurcations of SFRs and SSFRs depending on the neighbor's morphology seem to occur at R_n~0.5r_vir,nei, which is supported by K-S and MC tests with a significance level of >98%. For all redshift bins, the SSFRs of late-type galaxies interacting with late-type neighbors are increased by factors of about 1.8+/-0.7 and 4.0+/-1.2 compared to those of non-interacting galaxies when the pair separation is smaller than 0.5r_vir,nei and 0.1r_vir,nei, respectively. The Tdust of both local and high-z late-type galaxies strongly interacting with late-type neighbors (i.e. R_n < 0.1r_vir,nei) appears to be higher than that of non-interacting galaxies with a significance level of 96-99%. However, Tdust of local late-type galaxies strongly interacting with early-type neighbors seems to be lower than or similar to that of non-interacting galaxies. Our results suggest that galaxy-galaxy interactions and mergers have been strongly affecting the SFR and the dust properties of star-forming galaxies over at least 8 billion years.

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

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

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.

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

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

We take advantage of the sensitivity and resolution of Herschel at 100 and 160 micron to directly image the thermal dust emission and investigate the infrared luminosities, L(IR), and dust obscuration of typical star-forming (L*) galaxies at high redshift. Our sample consists of 146 UV-selected galaxies with spectroscopic redshifts 1.5<z<2.6 in the GOODS-North field. Supplemented with deep Very Large Array (VLA) and Spitzer imaging, we construct median stacks at the positions of these galaxies at 24, 100, and 160 micron, and 1.4 GHz. The comparison between these stacked fluxes and a variety of dust templates and calibrations implies that typical star-forming galaxies with UV luminosities L(UV)>1e10 Lsun at z~2 are luminous infrared galaxies (LIRGs) with a median L(IR)=(2.2+/-0.3)e11 Lsun. Typical galaxies at 1.5<z<2.6 have a median dust obscuration L(IR)/L(UV) = 7.1+/-1.1, which corresponds to a dust correction factor, required to recover the bolometric star formation rate (SFR) from the unobscured UV SFR, of 5.2+/-0.6. This result is similar to that inferred from previous investigations of the UV, H-alpha, 24 micron, radio, and X-ray properties of the same galaxies studied here. Stacking in bins of UV slope implies that L* galaxies with redder spectral slopes are also dustier, and that the correlation between UV slope and dustiness is similar to that found for local starburst galaxies. Hence, the rest-frame 30 and 50 micron fluxes validate on average the use of the local UV attenuation curve to recover the dust attenuation of typical star-forming galaxies at high redshift. In the simplest interpretation, the agreement between the local and high redshift UV attenuation curves suggests a similarity in the dust production and stellar and dust geometries of starburst galaxies over the last 10 billion years.

Dust attenuation curves in external galaxies are useful to study their dust properties as well as to interpret their intrinsic spectral energy distributions. In particular the presence or absence of a UV bump at 2175 A remains an open issue which has consequences on the interpretation of broad band colours of distant galaxies. We study the dust attenuation curve in the UV range at z >1. In particular we search for the presence of a UV bump. We use deep photometric data of the CDFS obtained with intermediate and broad band filters by the MUSYC project to sample the UV rest-frame of galaxies with 1<z <2. Herschel/PACS and Spitzer/MIPS data are used to measure the dust emission. 30 galaxies were selected with high S/N in all bands. Their SEDs from the UV to the far-IR are fitted using the CIGALE code and the characteristics of the dust attenuation curve are obtained. The mean dust attenuation curve we derive exhibits a UV bump at 2175A whose amplitude corresponds to 35 % (76%) that of the MW (LMC2 supershell) extinction curve. An analytical expression of the average attenuation curve is given, it is found slightly steeper than the Calzetti et al. one, although at a 1 sigma level. Our galaxy sample is used to study the derivation of the slopes of the UV continuum from broad band colours, including the GALEX FUV-NUV colour. Systematic errors induced by the presence of the bump are quantified. We compare dust attenuation factors measured with CIGALE to the slope of the UV continuum and find that there is a large scatter around the relation valid for local starbursts (0.7 mag). The uncertainties on the determination of the UV slope lead to an extra systematic error of the order of 0.3 to 0.7 mag on dust attenuation when a filter overlaps the UV bump.

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

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

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

We use the 4Ms CDF-S survey to place direct X-ray constraints on the ubiquity of z~2 heavily obscured AGNs in K<22 BzK galaxies. Forty seven of the 222 BzK galaxies in the central region of the CDF-S are detected at X-ray energies, 11 of which have hard X-ray spectral slopes (Gamma<1) indicating the presence of heavily obscured AGN activity. The other 36 X-ray detected BzK galaxies appear to be relatively unobscured AGNs and starburst galaxies; we use X-ray variability analyses over a rest-frame baseline of ~3 years to further confirm the presence of AGN activity in many of these systems. The majority (7 out of 11) of the heavily obscured AGNs have excess IR emission over that expected from star formation (termed "IR-excess galaxies"). However, we find that X-ray detected heavily obscured AGNs only comprise ~25% of the IR-excess galaxy population, which is otherwise composed of relatively unobscured AGNs and starburst galaxies. We find that the typical X-ray spectrum of the heavily obscured AGNs is better characterized by a pure reflection model than an absorbed power-law model, suggesting extreme Compton-thick absorption in some systems. We verify this result by producing a composite rest-frame 2-20 keV spectrum, which has a similar shape as a reflection-dominated X-ray spectrum and reveals an emission feature at rest-frame energy ~6.4 keV, likely to be due to Fe K. These heavily obscured AGNs are likely to be the distant analogs of the reflection-dominated AGNs recently identified at z~0 with >10 keV observatories. On the basis of these analyses we estimate the space density for typical (intrinsic X-ray luminosities of L_X>1E43 erg/s) heavily obscured and Compton-thick AGNs at z~2. Our space-density constraints are conservative lower limits but they are already consistent with the range of predictions from X-ray background models.

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

We present observations of the molecular gas in the GN20 proto-cluster of galaxies at $z =4.05$ using the Expanded Very Large Array (EVLA). This group of galaxies is the ideal laboratory for studying the formation of massive galaxies via luminous, gas-rich starbursts within 1.6 Gyr of the Big Bang. We detect three galaxies in the proto-cluster in CO 2-1 emission, with gas masses (H$_2$) between $10^{10}$ and $10^{11} \times (\alpha/0.8)$ M$_\odot$. The emission from the brightest source, GN20, is resolved with a size $\sim 2"$, and has a clear north-south velocity gradient, possibly indicating ordered rotation. The gas mass in GN20 is comparable to the stellar mass ($1.3\times 10^{11} \times (\alpha/0.8)$ M$_\odot$ and $2.3\times 10^{11}$ M$_\odot$, respectively), and the sum of gas plus stellar mass is comparable to the dynamical mass of the system ($\sim 3.4\times 10^{11} [sin(i)/sin(45^o)]^{-2}$ M$_\odot$), within a 5kpc radius. There is also evidence for a tidal tail extending another $2"$ north of the galaxy with a narrow velocity dispersion. GN20 may be a massive, gas rich disk that is gravitationally disturbed, but not completely disrupted. There is one Lyman-break galaxy (BD29079) in the GN20 proto-cluster with an optical spectroscopic redshift within our search volume, and we set a 3$\sigma$ limit to the molecular gas mass of this galaxy of $1.1\times 10^{10} \times (\alpha/0.8)$ M$_\odot$.

Line and continuum studies at centimeter through submillimeter wavelengths address probe deep into the earliest, most active and dust obscured phases of galaxy formation, and reveal the molecular and cool atomic gas. We summarize the techniques of radio astronomy to perform these studies, then review the progress on radio studies of galaxy formation. The dominant work over the last decade has focused on massive, luminous starburst galaxies (submm galaxies and AGN host galaxies). The far infrared luminosities are ~ 1e13 Lsun, implying star formation rates, SFR > 1e3 Msun/year. Molecular gas reservoirs are found with masses: M(H_2) > 1e10 (alpha/0.8}) Msun. The CO excitation in these luminous systems is much higher than in low redshift spiral galaxies. Imaging of the gas distribution and dynamics suggests strongly interacting and merging galaxies, indicating gravitationally induced, short duration (~ 1e7 year) starbursts. These systems correspond to a major star formation episode in massive galaxies in proto-clusters at intermediate to high redshift. Recently, radio observations have probed the more typical star forming galaxy population (SFR ~ 100 Msun/year), during the peak epoch of Universal star formation (z ~ 1.5 to 2.5). These observations reveal massive gas reservoirs without hyper-starbursts, and show that active star formation occurs over a wide range in galaxy stellar mass. The conditions in this gas are comparable to those found in the Milky Way disk. A key result is that the peak epoch of star formation in the Universe also corresponds to an epoch when the baryon content of star forming galaxies was dominated by molecular gas, not stars. We consider the possibility of tracing out the dense gas history of the Universe, and perform initial, admittedly gross, calculations. ABRIDGED

Tracing the formation and evolution of all supermassive black holes, including the obscured ones, understanding how black holes influence their surroundings and how matter behaves under extreme conditions, are recognized as key science objectives to be addressed by the next generation of instruments. These are the main goals of the COSPIX proposal, made to ESA in December 2010 in the context of its call for selection of the M3 mission. In addition, COSPIX, will also provide key measurements on the non thermal Universe, particularly in relation to the question of the acceleration of particles, as well as on many other fundamental questions as for example the energetic particle content of clusters of galaxies. COSPIX is proposed as an observatory operating from 0.3 to more than 100 keV. The payload features a single long focal length focusing telescope offering an effective area close to ten times larger than any scheduled focusing mission at 30 keV, an angular resolution better than 20 arcseconds in hard X-rays, and polarimetric capabilities within the same focal plane instrumentation. In this paper, we describe the science objectives of the mission, its baseline design, and its performances, as proposed to ESA.

In order to understand galaxy evolution through cosmic times it is critical to derive the properties of the molecular gas content of galaxies, i.e. the material out of which stars ultimately form. The last decade has seen rapid progress in this area, with the detection of massive molecular gas reservoirs at high redshifts in submillimeter-selected galaxies and quasars. In the latter case, molecular gas reservoirs have been quantified out to redshifts z>6, i.e. towards the end of cosmic reionization when the universe was less than one Gyr old. The recent discovery of molecular gas in more normal galaxies have extended these studies from the most extreme objects in the universe (SFR~1000M_sun/yr; quasars and submillimeter galaxies) to more 'normal' starforming systems at redshifts 1.5-2.5 (with SFR~100M_sun/yr). However, detecting the molecular gas reservoirs of high-redshift galaxies that only have moderate star formation rates (~<10 M_sun/yr, similar to the faint galaxies seen in the Hubble Ultra Deep Field) will likely have to await the completion of ALMA.

We report the detection of the 6.4 keV Iron K{\alpha} emission line in two infrared-luminous, massive, star-forming BzK galaxies at z = 2.578 and z = 2.90 in the CDF-S. The Chandra 4 Ms spectra of BzK4892 and BzK8608 show a reflection dominated continuum with strong Iron lines, with rest-frame equivalent widths EW-2.3 keV and 1.2 keV, respectively, demonstrating Compton thick obscuration of the central AGN. For BzK8608 the line identification closely matches the existing photometric redshift derived from the stellar emission. We use the observed luminosities of the Iron K{\alpha} line, of the rest-frame mid-IR continuum and of the UV rest-frame narrow emission lines to infer intrinsic L(2-10 keV)>1e44 erg/s, about 1.0-2.5 dex larger than the observed ones, hence confirming the presence of an absorber with N(H) > 1e24 cm-2. The two BzK galaxies have stellar masses of 5e10 M\odot and, based on VLA 1.4 GHz and submm 870{\mu}m observations, they appear to host vigorous starburst activity with SFR - 300-700 M\odot/yr that is also optically thick. We estimate that the AGN might also conceivably account for an important fraction of the bolometric far-IR emission of the galaxies. The implied volume density of Compton thick (CT) AGN with L(2-10 keV)>1e44 erg/s is in agreement with predictions from X-ray background synthesis models. These sources provide one of the first clearcut observations of the long-sought phase of simultaneous, heavily obscured quasar and star formation activity, predicted by models of massive galaxy evolution at high redshifts.

We report evidence of a fully established galaxy cluster at z=2.07, consisting of a ~20sigma overdensity of red, compact spheroidal galaxies spatially coinciding with extended X-ray emission detected with XMM-Newton. We use VLT VIMOS and FORS2 spectra and deep Subaru, VLT and Spitzer imaging to estimate the redshift of the structure from a prominent z=2.07 spectroscopic redshift spike of emission-line galaxies, concordant with the accurate 12-band photometric redshifts of the red galaxies. Using NICMOS and Keck AO observations, we find that the red galaxies have elliptical morphologies and compact cores. While they do not form a tight red sequence, their colours are consistent with that of a >1.3$~Gyr population observed at z~2.1. From an X-ray luminosity of .2*10^43 erg s^-1 and the stellar mass content of the red galaxy population, we estimate a halo mass of 5.3-8*10^13 Msun, comparable to the nearby Virgo cluster. These properties imply that this structure could be the most distant, mature cluster known to date and that X-ray luminous, elliptical-dominated clusters are already forming at substantially earlier epochs than previously known.

We use CFHTLS deep optical data, WIRCam Deep Survey (WIRDS) NIR data and XMM data to identify z>1.1 clusters in the CFHTLS D1 and D4 fields. Counterparts to such clusters can not be identified without deep NIR data and as such the total of =1deg2 of J , H & Ks band imaging provided by WIRDS is an indispensable tool in such work. Using public XMM X-ray data, we identify extended X-ray sources in the two fields. The resulting catalogue of extended X-ray sources was analyzed for optical/NIR counterparts, using a red-sequence algorithm. Redshifts of candidate groups and clusters were estimated using the median photometric redshifts of detected counterparts and where available spectroscopic data. Additionally, we surveyed X-ray point sources for potential group systems at the limit of our detection range in the X-ray data. A catalogue of z > 1.1 cluster candidates in the two fields has been compiled and cluster masses, radii and temperatures have been estimated using the scaling relations. The catalogue consists of 15 z > 1.1 candidates. Three of the detections are previously published extended X-ray sources. Of note is JKSC 041 for which we identify possible structures at z = 0.8, z = 0.96, z = 1.13 and z = 1.49. We also make an independent detection of the massive cluster, XMMXCS J2215.9-1738. We use the z > 1.1 catalogue to compare the cluster number counts in these fields with models based on WMAP 7-year cosmology and find that the models slightly over-predict the observations, whilst at z>1.5 we do not detect any clusters. We note that cluster number counts at z > 1.1 are highly sensitive to the cosmological model, however a significant reduction in present statistical (due to available survey area) and systematic (due to cluster scaling relations) uncertainties is required in order to confidently constrain cosmological parameters using cluster number counts at high redshift.

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.