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 study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843~MHz radio AGN population in galaxy clusters from two large catalogs created using the Dark Energy Survey (DES): $\sim$11,800 optically selected RM-Y3 and $\sim$1,000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGN are highly concentrated around cluster centers to $z\sim1$. We measure the halo occupation number for cluster radio AGN above a threshold luminosity, finding that the number of radio AGN per cluster increases with cluster halo mass as $N\propto M^{1.2\pm0.1}$ ($N\propto M^{0.68\pm0.34}$) for the RM-Y3 (MARD-Y3) sample. Together, these results indicate that radio mode feedback is favoured in more massive galaxy clusters. Using optical counterparts for these sources, we demonstrate weak redshift evolution in the host broad band colors and the radio luminosity at fixed host galaxy stellar mass. We use the redshift evolution in radio luminosity to break the degeneracy between density and luminosity evolution scenarios in the redshift trend of the radio AGN luminosity function (LF). The LF exhibits a redshift trend of the form $(1+z)^\gamma$ in density and luminosity, respectively, of $\gamma_{\rm D}=3.0\pm0.4$ and $\gamma_{\rm P}=0.21\pm0.15$ in the RM-Y3 sample, and $\gamma_{\rm D}=2.6\pm0.7$ and $\gamma_{\rm P}=0.31\pm0.15$ in MARD-Y3. We discuss the physical drivers of radio mode feedback in cluster AGN, and we use the cluster radio galaxy LF to estimate the average radio-mode feedback energy as a function of cluster mass and redshift and compare it to the core ($<0.1R_{500}$) X-ray radiative losses for clusters at $z<1$.

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.

By using the deepest available mid and far infrared surveys in the CANDELS, GOODS and COSMOS fields we study the evolution of the Main Sequence (MS) of star forming galaxies (SFGs) from z~0 to` ~2.5 at stellar masses larger than 10^{10} M_{\odot}. The MS slope and scatter are consistent with a re-scaled version of the local relation and distribution, shifted at higher values of SFR according to ~(1+z)^{3.2}. The relation exhibits a bending at the high mass end and a slightly increasing scatter as a function of the stellar mass. We show that the previously reported evolution of the MS slope, in the considered mass and redshift range, is due to a selection effect. The distribution of galaxies in the MS region at fixed stellar mass is well represented by a single log-normal distribution at all redshifts and masses, with starburst galaxies (SBs) occupying the tail at high SFR.

Our current knowledge of cosmic star-formation history during the first two billion years (corresponding to redshift z >3) is mainly based on galaxies identified in rest-frame ultraviolet light. However, this population of galaxies is known to under-represent the most massive galaxies, which have rich dust content and/or old stellar populations. This raises the questions of the true abundance of massive galaxies and the star-formation-rate density in the early universe. Although several massive galaxies that are invisible in the ultraviolet have recently been confirmed at early epochs, most of them are extreme starbursts with star-formation rates exceeding 1000 solar masses per year, suggesting that they are unlikely to represent the bulk population of massive galaxies. Here we report submillimeter (wavelength 870um) detections of 39 massive star-forming galaxies at z > 3, which are unseen in the spectral region from the deepest ultraviolet to the near-infrared. With a space density of about $2 \times 10^{-5}$ per cubic megaparsec (two orders of magnitudes higher than extreme starbursts) and star-formation rates of 200 solar masses per year, these galaxies represent the bulk population of massive galaxies that have been missed from previous surveys. They contribute a total star-formation-rate density ten times larger than that of equivalently massive ultraviolet-bright galaxies at z >3. Residing in the most massive dark matter halos at their redshifts, they are probably the progenitors of the largest present-day galaxies in massive groups and clusters. Such a high abundance of massive and dusty galaxies in the early universe challenges our understanding of massive-galaxy formation.

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 first results from a galaxy population study in the highest redshift galaxy clusters identified in the 2500 deg$^2$ South Pole Telescope Sunyaev Zel'dovich effect (SPT-SZ) survey. The cluster selection is to first order independent of galaxy properties, making the SPT-SZ sample particularly well suited for cluster galaxy population studies. We carry out a 4-band imaging campaign with the {\it Hubble} and {\it Spitzer} Space Telescopes of the five $z\gtrsim 1.4$, S/N$_{SZE}>$5 clusters, that are among the rarest most massive clusters known at this redshift. All five show clear overdensities of red galaxies whose colors agree with the initial cluster redshift estimates. The highest redshift cluster in this sample, SPT-CLJ0459-4947 at $z\sim1.72$, is the most distant $M_{500}>10^{14}~M_{\odot}$ ICM-selected cluster discovered thus far, and is one of only three known clusters in this mass range at $z\gtrsim 1.7$, regardless of selection. Based on UVJ-like photometric classification of quiescent and star-forming galaxies, the passive fraction in the cluster central regions ($r/r_{500}<0.7$) is higher than in the field at the same redshift, with corresponding environmental quenching efficiencies typically in the range $\sim0.5-0.8$ for stellar masses $\log(M/M_{\odot})>10.85$. We have explored the impact of emission from star formation on the selection of this sample, concluding that all five clusters studied here would still have been detected with S/N$_{SZE}>$5, even if they had the same passive fraction as measured in the field. Our results thus point towards an efficient suppression of star formation in massive galaxies in the central regions of the most massive clusters, occurring already earlier than $z\sim1.5$. [Abridged]

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.

By using a set of different SFR indicators, including WISE mid-infrared and Halpha emission, we study the slope of the Main Sequence (MS) of local star forming galaxies at stellar masses larger than 10^{10} M_{\odot}. The slope of the relation strongly depends on the SFR indicator used. In all cases, the local MS shows a bending at high stellar masses with respect to the slope obtained in the low mass regime. While the distribution of galaxies in the upper envelope of the MS is consistent with a log-normal distribution, the lower envelope shows an excess of galaxies, which increases as a function of the stellar mass but varies as a function of the SFR indicator used. The scatter of the best log-normal distribution increases with stellar mass from ~0.3 dex at 10^{10} M_{\odot} to ~0.45 at 10^{11} M_{\odot}. The MS high-mass end is dominated by central galaxies of group sized halos with a red bulge and a disk redder than the lower mass counterparts. We argue that the MS bending in this region is due to two processes: i) the formation of a bulge component as a consequence of the increased merger activity in groups, and ii) the cold gas starvation induced by the hot halo environment, which cuts off the gas inflow onto the disk. Similarly, the increase of the MS scatter at high stellar masses would be explained by the larger spread of star formation histories of central group and cluster galaxies with respect to lower mass systems.

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 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 study the molecular gas properties of two star-forming galaxies separated by 6 kpc in the projected space and belonging to a galaxy cluster selected from the Irac Shallow Cluster Survey, at a redshift $z=1.2$, i.e., $\sim2$ Gyr after the cosmic star formation density peak. This work describes the first CO detection from $1<z<1.4$ star forming cluster galaxies with no reported clear evidence of AGN. We exploit observations taken with the NOEMA interferometer at $\sim3$ mm to detect CO(2-1) line emission from the two selected galaxies, unresolved by our observations. Based on the CO(2-1) spectrum we estimate a total molecular gas mass $M({\rm H_2})=(2.2^{+0.5}_{-0.4})\times10^{10}$ $M_\odot$ and dust mass $M_{\rm dust}<4.2\times10^8~M_\odot$ for the two blended sources. The two galaxies have similar stellar masses and a large relative velocity of $\sim$400 km/s estimated from the CO(2-1) line width. These findings tend to privilege a scenario where both sources contribute to the observed CO(2-1). By using the archival Spitzer MIPS flux at 24$\mu$m we estimate an ${\rm SFR(24\mu m)}=(28^{+12}_{-8})~M_\odot$/yr for each of the two galaxies. Assuming that the two sources equally contribute to the observe CO(2-1) our analysis yields a depletion time scale $\tau_{\rm dep}=(3.9^{+1.4}_{-1.8})\times10^8$ yr, and a molecular gas to stellar mass ratio $0.17\pm0.13$ for each of two sources, separately. Our results are in overall agreement with those of other distant cluster galaxies. The two target galaxies have molecular gas mass and depletion time that are marginally compatible with, but smaller than those of main sequence field galaxies, suggesting that the molecular gas has not been refueled enough. Higher resolution and higher frequency observations will enable us to spatially resolve the two sources and possibly distinguish between different gas processing mechanisms.

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

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection.

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.

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

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 obtained ALMA spectroscopy and imaging to investigate the origin of the unexpected sub-mm emission toward the most distant quiescent galaxy known to date, ZF-COSMOS-20115 at z=3.717. We show here that this sub-mm emission is produced by another massive, compact and extremely obscured galaxy, located only 3.1 kpc away from the quiescent galaxy. We dub the quiescent and dusty galaxies Jekyll and Hyde, respectively. No dust emission is detected at the location of Jekyll, implying SFR < 13 Msun/yr, which is the most stringent upper limit ever obtained for a quiescent galaxy at these redshifts. The two sources are confirmed to lie at the same redshift thanks to the detection of [CII]158 in Hyde, which provides one the few robust redshifts for an "H-dropout" galaxy. The line has a rotating-disk velocity profile blueshifted from Jekyll by 549+/-60 km/s, demonstrating that it is produced by another galaxy. Careful de-blending of the Spitzer imaging confirms the existence of Hyde, and its non-detection with Hubble requires extreme attenuation by dust. Modeling the photometry of both galaxies shows that Jekyll has fully quenched >200 Myr prior to observation and still presents a challenge for models, while Hyde only harbors moderate star-formation (SFR<120 Msun/yr) and is located at least a factor 1.4 below the z~4 main sequence. Hyde could also have stopped forming stars <200 Myr before being observed, which would be consistent with its hight compactness similar to z~4 quiescent galaxies and its low [CII]/FIR ratio, but significant SF cannot be ruled out. Finally, we show that Hyde hosts a dense reservoir of gas comparable to that of extreme starbursts, suggesting that its SFR was reduced without expelling the gas outside of the galaxy. We argue that Jekyll and Hyde can be seen as two stages of quenching, and provide a unique laboratory to study this poorly understood phenomenon. [abridged]

We present a new, publicly available library of dust spectral energy distributions (SEDs). These SEDs are characterized by only three parameters: the dust mass (Mdust), the dust temperature (Tdust), and the mid-to-total infrared color (IR8=LIR/L8). The latter measures the relative contribution of PAH molecules to the LIR. We used this library to model star-forming galaxies at 0.5<z<4 in the CANDELS fields, using both individual detections and stacks of Herschel and ALMA imaging, and extending this sample to z=0 using the Herschel Reference Survey. At first order, the dust SED of a galaxy was found independent of stellar mass, but evolving with redshift. We found trends of increasing Tdust and IR8 with redshift and distance from the SFR--Mstar main sequence (MS), and quantified for the first time their intrinsic scatter. Half of the observed variation of these parameters was captured by these empirical relations, with residual scatters of 12% and 0.18 dex, respectively. Second order variations with stellar mass are discussed. Building on these results, we constructed high-fidelity mock galaxy catalogs to predict the accuracy of LIR and Mdust determined from a single flux measurement. Using a single JWST MIRI band, we found that LIR is typically uncertain by 0.15 dex, with a maximum of 0.25 dex when probing the rest-frame 8 um, and this is not significantly impacted by typical redshift uncertainties. On the other hand, we found that ALMA bands 8-to-7 and 6-to-3 measure Mdust at better than 0.2 and 0.15 dex, respectively, and independently of redshift, while bands 9-to-6 only measure LIR at better than 0.2 dex at z>1, 3.2, 3.8, and 5.7, respectively. Starburst galaxies above the MS have LIR significantly underestimated and Mdust overestimated. These results can be used immediately to interpret more accurately the large amount of archival data from Spitzer, Herschel and ALMA. [abridged]

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.

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 test the impact of using variable star forming histories (SFHs) and the use of the IR luminosity (LIR) as a constrain on the physical parameters of high redshift dusty star-forming galaxies. We explore in particular the stellar properties of galaxies in relation with their location on the SFR-M* diagram. We perform SED fitting of the UV-NIR and FIR emissions of a large sample of GOODS-Herschel galaxies, for which rich multi-wavelength observations are available. We test different SFHs and imposing energy conservation in the SED fitting process, to face issues like the age-extinction degeneracy and produce SEDs consistent with observations. Our models work well for the majority of the sample, with the notable exception of the high LIR end, for which we have indications that our simple energy conservation approach cannot hold true. We find trends in the SFHs fitting our sources depending on stellar mass M* and z. Trends also emerge in the characteristic timescales of the SED models depending on the location on the SFR-M* diagram. We show that whilst using the same available observational data, we can produce galaxies less star-forming than usually inferred, if we allow declining SFHs, while properly reproducing their observables. These sources can be post-starbursts undergoing quenching, and their SFRs are potentially overestimated if inferred from their LIR. Fitting without the IR constrain leads to a strong preference for declining SFHs, while its inclusion increases the preference of rising SFHs, more so at high z, in tentative agreement with the cosmic star formation history. Keeping in mind that the sample is biased towards high LIR, the evolution shaped by our model appears as both bursty (initially) and steady-lasting (later on). The global SFH of the sample follows the cosmic SFH with a small scatter, and is compatible with the "downsizing" scenario of galaxy evolution.

This paper introduces EGG, the Empirical Galaxy Generator, a tool designed within the ASTRODEEP collaboration to generate mock galaxy catalogs for deep fields with realistic fluxes and simple morphologies. The simulation procedure is based exclusively on empirical prescriptions -- rather than first principles -- to provide the most accurate match with observations at 0<z<7. In particular, we consider that galaxies can be either quiescent or star-forming, and use their stellar mass (M*) and redshift (z) as the fundamental properties from which all the other observables can be statistically derived. Drawing z and M* from the observed galaxy stellar mass functions, we associate a star formation rate (SFR) to each galaxy from the tight SFR-M* main sequence, while dust attenuation, optical colors and morphologies (including bulge-to-total ratios, sizes and aspect ratios) are obtained from empirical relations that we establish from the high quality Hubble and Herschel observations available in the CANDELS fields. Random scatter is introduced in each step to reproduce the observed distributions of each parameter. Based on these observables, a panchromatic spectral energy distribution (SED) is selected for each galaxy and synthetic photometry is produced by integrating the redshifted SED in common broad-band filters. Finally, the mock galaxies are placed on the sky at random positions with a fixed angular two-point correlation function to implement basic clustering. The resulting flux catalogs reproduce accurately the observed number counts in all broad bands from the ultraviolet up to the sub-millimeter, and can be directly fed to image simulators such as Skymaker. The images can then be used to test source extraction softwares and image-based techniques such as stacking. EGG is open-source, and is made available to the community together with a set of pre-generated catalogs and images.

We introduce the ALMA Redshift 4 Survey (AR4S), a systematic ALMA survey of all the known galaxies with stellar mass (M*) larger than 5e10 Msun at 3.5<z<5 in the GOODS--south, UDS and COSMOS CANDELS fields. The sample we have analyzed in this paper is composed of 96 galaxies observed with ALMA at 890um (180um rest-frame) with an on-source integration time of 1.3 min per galaxy. We detected 32% of the sample at more than 3 sigma significance. Using the stacked ALMA and Herschel photometry, we derived an average dust temperature of 40+/-2 K for the whole sample, and extrapolate the Lir and SFR for all our galaxies based on their ALMA flux. We then used a forward modeling approach to estimate their intrinsic sSFR distribution, deconvolved of measurement errors and selection effects: we find a linear relation between SFR and M*, with a median sSFR=2.8+/-0.8 Gyr and a dispersion around that relation of 0.28+/-0.13 dex. This latter value is consistent with that measured at lower redshifts, which is proof that the main sequence of star-forming galaxies was already in place at z=4, at least among massive galaxies. These new constraints on the properties of the main sequence are in good agreement with the latest predictions from numerical simulations, and suggest that the bulk of star formation in galaxies is driven by the same mechanism from z=4 to the present day, that is, over at least 90% of the cosmic history. We also discuss the consequences of our results on the population of early quiescent galaxies. This paper is part of a series that will employ these new ALMA observations to explore the star formation and dust properties of the massive end of the z=4 galaxy population.

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.

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge.

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

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 present a new method to search for candidate z~>2 Herschel 500{\mu}m sources in the GOODS-North field, using a S500{\mu}m/S24{\mu}m "color deconfusion" technique. Potential high-z sources are selected against low-redshift ones from their large 500{\mu}m to 24{\mu}m flux density ratios. By effectively reducing the contribution from low-redshift populations to the observed 500{\mu}m emission, we are able to identify counterparts to high-z 500{\mu}m sources whose 24{\mu}m fluxes are relatively faint. The recovery of known z~4 starbursts confirms the efficiency of this approach in selecting high-z Herschel sources. The resulting sample consists of 34 dusty star-forming galaxies at z~>2. The inferred infrared luminosities are in the range 1.5x10^12-1.8x10^13 Lsun, corresponding to dust-obscured star formation rates (SFRs) of ~260-3100 Msun/yr for a Salpeter IMF. Comparison with previous SCUBA 850{\mu}m-selected galaxy samples shows that our method is more efficient at selecting high-z dusty galaxies with a median redshift of z=3.07+/-0.83 and 10 of the sources at z~>4. We find that at a fixed luminosity, the dust temperature is ~5K cooler than that expected from the Td-LIR relation at z<1, though different temperature selection effects should be taken into account. The radio-detected subsample (excluding three strong AGN) follows the far-infrared/radio correlation at lower redshifts, and no evolution with redshift is observed out to z~5, suggesting that the far-infrared emission is star formation dominated. The contribution of the high-z Herschel 500{\mu}m sources to the cosmic SFR density is comparable to that of SMG populations at z~2.5 and at least 40% of the extinction-corrected UV samples at z~4 (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 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.

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

We propose a new method of pushing $Herschel$ to its faintest detection limits using universal trends in the redshift evolution of the far infrared over 24$\mu$m colours in the well-sampled GOODS-North field. An extension to other fields with less multi-wavelength information is presented. This method is applied here to raise the contribution of individually detected $Herschel$ sources to the cosmic infrared background (CIRB) by a factor 5 close to its peak at 250$\mu$m and more than 3 in the 350$\mu$m and 500$\mu$m bands. We produce realistic mock $Herschel$ images of the deep PACS and SPIRE images of the GOODS-North field from the GOODS-$Herschel$ Key Program and use them to quantify the confusion noise at the position of individual sources, i.e., estimate a "local confusion noise". Two methods are used to identify sources with reliable photometric accuracy extracted using 24$\mu$m prior positions. The clean index (CI), previously defined but validated here with simulations, which measures the presence of bright 24$\mu$m neighbours and the photometric accuracy index (PAI) directly extracted from the mock $Herschel$ images. After correction for completeness, thanks to our mock $Herschel$ images, individually detected sources make up as much as 54% and 60% of the CIRB in the PACS bands down to 1.1 mJy at 100$\mu$m and 2.2 mJy at 160$\mu$m and 55, 33, and 13% of the CIRB in the SPIRE bands down to 2.5, 5, and 9 mJy at 250$\mu$m, 350$\mu$m, and 500$\mu$m, respectively. The latter depths improve the detection limits of $Herschel$ by factors of 5 at 250$\mu$m, and 3 at 350$\mu$m and 500$\mu$m as compared to the standard confusion limit. Interestingly, the dominant contributors to the CIRB in all $Herschel$ bands appear to be distant siblings of the Milky Way ($z$$\sim$0.96 for $\lambda$$<$300$\mu$m) with a stellar mass of $M_{\star}$$\sim$9$\times$10$^{10}$M$_{\odot}$.

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

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.

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.

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

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

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

We study dust attenuation at UV wavelengths at high redshift, where the UV is redshifted to the observed visible. In particular, we search for a UV bump and related implications for dust attenuation determinations. We use data in the CDFS, obtained in intermediate and broad band filters by the MUSYC project, to sample the UV rest-frame of 751 galaxies with 0.95<z<2.2. When available, Herschel/PACS data (GOODS-Herschel project), and Spitzer/MIPS measurements, are used to estimate the dust emission. The SED of each source is fit using the CIGALE code. The amount of dust attenuation and the dust attenuation curve are obtained as outputs of the SED fitting process, together with other parameters linked to the SFH. The global amount of dust attenuation at UV wavelengths is found to increase with stellar mass and to decrease as UV luminosity increases. A UV bump at 2175A is securely detected in 20% of the galaxies, and the mean amplitude of the bump for the sample is similar to that observed in the LMC supershell region. This amplitude is found to be lower in galaxies with very high SSFRs, and 90% of the galaxies exhibiting a secure bump are at z<1.5. The attenuation curve is confirmed to be steeper than that of local starburst galaxies for 20$% of the galaxies. The large dispersion found for these two parameters describing the attenuation law is likely to reflect a wide diversity of attenuation laws among galaxies. The relations between dust attenuation, IR-to-UV flux ratio, and the slope of the UV continuum are derived for the mean attenuation curve found for our sample. Deviations from the average trends are found to correlate with the age of the young stellar population and the shape of the attenuation curve.(abriged)

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.

We present detection and analysis of faint X-ray sources in the Chandra deep field south (CDFS) using the 4 Msec Chandra observation and adopting a new detection algorithm, based on a targeted search at the position of known high-z galaxies. This optimized technique results in the identification of 54 z>3 AGNs, 29 of which are new detections. Applying stringent completeness criteria, we derive AGN luminosity functions in the redshift bins 3-4, 4-5 and >5.8 and for 42.75<log L(2-10 keV)<44.5. We join this data with the luminous AGN luminosity functions from optical surveys and find that the evolution of the high-z, wide luminosity range luminosity function can be best modeled by pure luminosity evolution with L* decreasing from 6x10^44 ergs/s at z=3 to L*=2x10^44 ergs/s at z=6. We compare the high-z luminosity function with the prediction of theoretical models using galaxy interactions as AGN triggering mechanism. We find that these models are broadly able to reproduce the high-z AGN luminosity functions. A better agreement is found assuming a minimum Dark Matter halo mass for black hole formation and growth. We compare our AGN luminosity functions with galaxy mass functions to derive high-z AGN duty cycle using observed Eddington ratio distributions to derive black hole masses. We find that the duty cycle increases with galaxy stellar mass and redshift by a factor 10-30 from z=0.25 to z=4-5. We also report on the detection of a large fraction of highly obscured, Compton thick AGN at z>3 (18+17-10%). Their optical counterparts are not strongly reddened and we thus conclude that the size of the X-ray absorber is likely smaller than the dust sublimation radius. We finally report the discovery of a highly star-forming galaxy at z=3.47. If confirmed, this would be one of the farthest objects in which stellar sources are detected in X-rays.

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

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

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

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

We report on 2 fossil groups of galaxies at z=0.425 and 0.372 discovered in the Cosmic Evolution Survey (COSMOS) area. Selected as X-ray extended sources, they have total masses (M_200) of 1.9(+/-0.41)E13 and 9.5(+/-0.42)E13 M_sun, respectively, as obtained from a recent X-ray luminosity-mass scaling relation. The lower mass system appears isolated, whereas the other sits in a well-known large-scale structure (LSS) populated by 27 other X-ray emitting groups. The identification as fossil is based on the i-band photometry of all the galaxies with a photo-z consistent with that of the group at the 2-sigma confidence level and within a projected group-centric distance equal to 0.5R_200, and i_AB<=22.5-mag limited spectroscopy. Both fossil groups exhibit high stellar-to-total mass ratios compared to all the X-ray selected groups of similar mass at 0.3<=z<=0.5 in the COSMOS. At variance with the composite galaxy stellar mass functions (GSMFs) of similarly massive systems, both fossil group GSMFs are dominated by passively evolving galaxies down to M^stars~1E10 M_sun (according to the galaxy broad-band spectral energy distributions). The relative lack of star-forming galaxies with 1E10<=M^stars<=1E11 M_sun is confirmed by the galaxy distribution in the b-r vs i color-magnitude diagram. Hence, the 2 fossil groups appear as more mature than the coeval, similarly massive groups. Their overall star formation activity ended rapidly after an accelerated build up of the total stellar mass; no significant infall of galaxies with M^stars>=1E10 M_sun took place in the last 3 to 6 Gyr. This similarity holds although the 2 fossil groups are embedded in two very different density environments of the LSS, which suggests that their galaxy populations were shaped by processes that do not depend on the LSS. However, their progenitors may do so. ...

We explore the evolution of the specific star formation rate (SSFR) for 3.6um-selected galaxies of different M_* in the COSMOS field. The average SFR for sub-sets of these galaxies is estimated with stacked 1.4GHz radio continuum emission. We separately consider the total sample and a subset of galaxies (SF) that shows evidence for substantive recent star formation in the rest-frame optical SED. At 0.2<z<3 both populations show a strong and M_*-independent decrease in their SSFR towards z=0.2, best described by a power- law (1+z)^n, where n~4.3 for all galaxies and n~3.5 for SF sources. The decrease appears to have started at z>2, at least above 4x10^10M_Sun where our conclusions are most robust. We find a tight correlation with power-law dependence, SSFR (M_*)^beta, between SSFR and M_* at all z. It tends to flatten below ~10^10M_Sun if quiescent galaxies are included; if they are excluded a shallow index beta_SFG -0.4 fits the correlation. On average, higher M_* objects always have lower SSFRs, also among SF galaxies. At z>1.5 there is tentative evidence for an upper SSFR-limit that an average galaxy cannot exceed. It is suggested by a flattening of the SSFR-M_* relation (also for SF sources), but affects massive (>10^10M_Sun) galaxies only at the highest z. Below z=1.5 there thus is no direct evidence that galaxies of higher M_* experience a more rapid waning of their SSFR than lower M_* SF systems. In this sense, the data rule out any strong 'downsizing'. We combine our results with recent measurements of the galaxy (stellar) mass function in order to determine the characteristic mass of a SF galaxy (M_*=10^(10.6\pm0.4)M_Sun). In this sense, too, there is no 'downsizing'. Our analysis constitutes the most extensive SFR density determination with a single technique to z=3. Recent Herschel results are consistent with our results, but rely on far smaller samples.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We present a multi-wavelength study of galaxy populations in the core of the massive, X-ray luminous cluster XMMU J2235 at z=1.39, based on VLT and HST optical and near-infrared photometry. Luminosity functions in the z, H, and Ks bands show a faint-end slope consistent with being flat, and a characteristic magnitude M* close to passive evolution predictions of M* of local massive clusters, with a formation redshift z>2. The color-magnitude and color-mass diagrams show evidence of a tight red sequence of massive galaxies, with overall old stellar populations, generally early-type morphology, typically showing early-type spectral features and rest-frame far-UV emission consistent with very low star formation rates (SFR<0.2Msun/yr). Star forming spectroscopic members, with SFRs of up to ~100Msun/yr, are all located at clustercentric distances >~250kpc, with the central cluster region already appearing effectively quenched. Massive galaxies in the core of this cluster appear to be in an advanced evolutionary stage in terms of both star formation and mass assembly. The high-mass end of the galaxy stellar mass function is essentially already in place, and the stellar mass fraction estimated within r500 (~1%, Kroupa IMF) is already similar to that of local massive clusters. On the other hand, morphological analysis of the massive red sequence galaxies suggests that they are smaller than similarly massive local early-types. While possibly affected by systematics and biases, this result might imply that, in spite of the overall early assembly of these sources, their evolution is not complete, and processes like minor (and likely dry) merging might still shape their structural properties to resemble those of their local counterparts, without substantially affecting their stellar mass or host stellar populations.[abridged]

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

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

We present a study of a 20cm selected sample in the Deep SWIRE VLA Field, reaching a limiting flux density of ~13.5 uJy at the image center. In a 0.6x0.6 square degrees field, we are able to assign an optical/IR counterpart to 97% of the radio sources. Up to 11 passbands from the NUV to 4.5um are then used to sample the spectral energy distribution (SED) of these counterparts in order to investigate the nature of the host galaxies. By means of an SED template library and stellar population synthesis models we estimate photometric redshifts, stellar masses, and stellar population properties, dividing the sample in three sub-classes of quiescent, intermediate and star-forming galaxies. We focus on the radio sample in the redshift range 0.3<z<1.3 where we estimate to have a redshift completeness higher than 90%, and study the properties and redshift evolution of these sub-populations. We find that, as expected, the relative contributions of AGN and star-forming galaxies to the uJy population depend on the flux density limit of the sample. At all flux levels a significant population of "green-valley" galaxies is observed. While the actual nature of these sources is not definitely understood, the results of this work may suggest that a significant fraction of faint radio sources might be composite (and possibly transition) objects, thus a simple "AGN vs star-forming" classification might not be appropriate to fully understand what faint radio populations really are.

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

The Cosmic Evolution Survey (COSMOS) allows for the first time a highly significant census of environments and structures up to redshift one, as well as a full morphological description of the galaxy population. In this paper we present a study aimed to constrain the evolution, in the redshift range 0.2 < z < 1.2, of the mass content of different morphological types and its dependence on the environmental density. We use a deep multicolor catalog, covering an area of ~0.7 square degrees inside the COSMOS field, with accurate photometric redshifts (i < 26.5 and dz/(z+1) ~ 0.035). We estimate galaxy stellar masses by fitting the multi-color photometry to a grid of composite stellar population models. We quantitatively describe the galaxy morphology by fitting PSF convolved Sersic profiles to the galaxy surface brightness distributions down to F814 = 24 mag for a sample of 41300 objects. We confirm an evolution of the morphological mix with redshift: the higher the redshift the more disk-dominated galaxies become important. We find that the morphological mix is a function of the local comoving density: the morphology density relation extends up to the highest redshift explored. The stellar mass function of disk-dominated galaxies is consistent with being constant with redshift. Conversely, the stellar mass function of bulge-dominated systems shows a decline in normalization with redshift. Such different behaviors of late-types and early-types stellar mass functions naturally set the redshift evolution of the transition mass. ABRIDGED

We present first results of a study aimed to constrain the star formation rate and dust content of galaxies at z~2. We use a sample of BzK-selected star-forming galaxies, drawn from the COSMOS survey, to perform a stacking analysis of their 1.4 GHz radio continuum as a function of different stellar population properties, after removing AGN contaminants from the sample. Dust unbiased star formation rates are derived from radio fluxes assuming the local radio-IR correlation. The main results of this work are: i) specific star formation rates are constant over about 1 dex in stellar mass and up to the highest stellar mass probed; ii) the dust attenuation is a strong function of galaxy stellar mass with more massive galaxies being more obscured than lower mass objects; iii) a single value of the UV extinction applied to all galaxies would lead to grossly underestimate the SFR in massive galaxies; iv) correcting the observed UV luminosities for dust attenuation based on the Calzetti recipe provide results in very good agreement with the radio derived ones; v) the mean specific star formation rate of our sample steadily decreases by a factor of ~4 with decreasing redshift from z=2.3 to 1.4 and a factor of ~40 down the local Universe. These empirical SFRs would cause galaxies to dramatically overgrow in mass if maintained all the way to low redshifts, we suggest that this does not happen because star formation is progressively quenched, likely starting from the most massive galaxies.

In this article we describe the search for white dwarfs (WDs) in the multi-band photometric data of the Capodimonte deep field survey. The WD candidates were selected through the V-R_C vs B-V color-color diagram. For two bright objects, the WD nature has been confirmed spectroscopically, and the atmospheric parameters (Teff and logg) have been determined. We have computed synthetic stellar population models for the observed field and the expected number of white dwarfs agrees with the observations. The possible contamination by turn-off and horizontal branch halo stars has been estimated. The quasar (QSO) contamination has been determined by comparing the number of WD candidates in different color bins with state-of-the-art models and previous observations. The WD space density is measured at different distances from the Sun. The total contamination (non-degenerate stars + QSOs) in our sample is estimated to be around 30%. This work should be considered a small experiment in view of more ambitious projects to be performed in the coming years in larger survey contexts.

We present a dynamical study of NGC 4697, an almost edge-on, intermediate-luminosity, E4 elliptical galaxy, combining new surface brightness photometry, new as well as published long-slit absorption line kinematic data, and published planetary nebulae (PNe) velocity data. The combined kinematic data set extends out to ~= 5' ~= 4.5 R_e and allows us to probe the galaxy's outer halo. For the first time, we model such a dataset with the new and flexible Chi^2-made-to-measure particle code NMAGIC. We extend NMAGIC to include seeing effects, introduce an efficient scheme to estimate the mass-to-light ratio, and incorporate a maximum likelihood technique to account for discrete velocity measurements. For modelling the PNe kinematics we use line-of-sight velocities and velocity dispersions computed on two different spatial grids, and we also use the individual velocity measurements with the likelihood method, in order to make sure that our results are not biased by the way we treat the PNe measurements. We generate axisymmetric self-consistent models as well models including various dark matter halos. These models fit all the mean velocity and velocity dispersion data with Chi^2/N<1, both in the case with only luminous matter and in potentials including quite massive halos. The likelihood analysis together with the velocity histograms suggest that models with low density halos such that the circular velocity v_c ~< 200 km/s at R_e are not consistent with the data. A range of massive halos with v_c ~= 250 km/s at 5R_e fit the PN data best. To derive stronger results would require PN velocities at even larger radii. The best fitting models are slightly radially anisotropic; the anisotropy parameter beta ~= 0.3 at the center, increasing to beta ~= 0.5 at radii >~ 2R_e.

In this paper we present a deep and homogeneous i-band selected multi-waveband catalogue in the COSMOS field covering an area of about 0.7 square-degree. Our catalogue with a formal 50 percent completeness limit for point sources of i~26.7 comprises about 290.000 galaxies with information in 8 passbands. We combine publicly available u, B, V, r, i, z, and K data with proprietary imaging in H band. We discuss in detail the observations, the data reduction, and the photometric properties of the H-band data. We estimate photometric redshifts for all the galaxies in the catalogue. A comparison with 162 spectroscopic redshifts in the redshift range 0 < z < 3 shows that the achieved accuracy of the photometric redshifts is (Delta_z / (z_spec+1)) ~0.035 with only ~2 percent outliers. We derive absolute UV magnitudes and investigate the evolution of the luminosity function evaluated in the rest-frame UV at 1500 Angstrom. There is a good agreement between the LFs derived here and the LFs derived in the FORS Deep Field. We see a similar brightening of M_star and a decrease of phi_star with redshift. The catalogue including the photometric redshift information is made publicly available.

Fundamental properties of the extinction curve, like the slope in the rest-frame ultraviolet (UV) and the presence/absence of a broad absorption excess centred at 2175 A (the UV bump), are investigated for a sample of 108 massive, star-forming galaxies at 1 < z < 2.5, selected from the FDF Spectroscopic Survey, the K20 survey, and the GDDS. These characteristics are constrained from a parametric description of the UV spectral energy distribution (SED) of a galaxy, as enforced by combined stellar population and radiative transfer models for different geometries, dust/stars configurations and dust properties. In at least one third of the sample, there is a robust evidence for extinction curves with at least a moderate UV bump. The presence of the carriers of the UV bump is more evident in galaxies with UV SEDs suffering from heavy reddening. We interpret these results as follows. The sample objects possess different mixtures of dust grains and molecules producing extinction curves in between the average ones of the Small and Large Magellanic Cloud, where the UV bump is absent or modest, respectively. Most of the dust embeds the UV-emitting stellar populations or is distributed out of the galaxy mid-plane. Alternatively, even dust with a pronounced UV bump, as for the average Milky-Way extinction curve, can be present and distributed in the galaxy mid-plane. In this case, variations of the continuum scattering albedo with wavelength or an age-dependent extinction are not sufficient to explain the previous trend with reddening. Hence, additional extraplanar dust has to be invoked. The data suggest that the carriers of the UV bump are associated with intermediate-age stellar populations, while they survive in the harshest UV-radiation fields owing to dust self-shielding. [abridged]

The near-infrared integral field spectrograph SINFONI at the ESO VLT opens a new window for the study of central supermassive black holes. With a near-IR spatial resolution similar to HST optical and the ability to penetrate dust it provides the possibility to explore the low-mass end of the M-sigma relation (sigma<120km/s) where so far very few black hole masses were measured with stellar dynamics. With SINFONI we observed the central region of the low-luminosity elliptical galaxy NGC4486a at a spatial resolution of ~0.1arcsec in the K band. The stellar kinematics was measured with a maximum penalised likelihood method considering the region around the CO absorption band heads. We determined a black hole mass of M_BH=1.25^{+0.75}_{-0.79} x 10^7 M_sun (90% C.L.) using the Schwarzschild orbit superposition method including the full 2-dimensional spatial information. This mass agrees with the predictions of the M-sigma relation, strengthening its validity at the lower sigma end.

Fundamental properties of the extinction curve, like the slope in the rest-frame UV and the presence/absence of a broad absorption excess centred at 2175A (the UV bump), are investigated for a sample of 108 massive, star-forming galaxies at 1<z<2.5, selected from the FDF Spectroscopic Survey, the K20 survey, and the GDDS. These characteristics are constrained from a parametric description of the UV spectral energy distribution of a galaxy. It turns out that the sample galaxies host dust producing extinction curves with properties in between those of the Small and Large Magellanic Clouds (SMC and LMC, respectively). LMC-like extinction curves, which exhibit a UV bump, are mainly found among highly-reddened, UV-ultraluminous galaxies at z~2.4 and highly-reddened, near-IR-bright, star-forming galaxies at z~1.2. We discuss star-formation rates, total stellar masses, the morphology, and the chemical properties of our sample galaxies with respect to possible explanations for the different extinction curves.

Our goal is to recover the Galactic Halo spatial density by means of field stars. To this aim, we apply a new technique to the Capodimonte Deep Field (OACDF, Alcala' et al. 2004), as a pilot study in view of the VLT Survey Telescope (VST) stellar projects. Considering the unique chance to collect deep and wide-field photometry with the VST, our method may represent a useful tool towards a definitive mapping of the Galactic Halo. In the framework of synthetic stellar populations, turn-off stars are used to reconstruct the spatial density. The determination of the space density is achieved by comparing the data with synthetic color-magnitude diagrams (CMDs). The only assumptions involve the IMF, age and metallicity of the synthetic halo population. Stars are randomly placed in the solid angle. The contributions of the various Monte Carlo distributions (with a step of 4 kpc) along the line of sight are simultaneously varied to reproduce the observed CMD. Our result on the space density is consistent with a power-law exponent n~3 over a range of Galactocentric distances from 8 to 40 kpc.

The multiple images observed in galaxy cluster Abell 1689 provide strong constraints not only on the mass distribution of the cluster but also on the ensemble properties of the cluster galaxies. Using parametric strong lensing models for the cluster, and by assuming well motivated scaling laws between the truncation radius s and the velocity dispersion sigma of a cluster galaxy we are able to derive sizes of the dark matter halos of cluster galaxies. For the scaling law expected for galaxies in the cluster environment (s propto sigma), we obtain s = 64^{+15}_{-14} (sigma / 220 km/s) kpc. For the scaling law used for galaxies in the field with s propto sigma^2 we find s = 66^{+18}_{-16} (sigma / 220 km/s)^2 kpc. Compared to halos of field galaxies, the cluster galaxy halos in Abell 1689 are strongly truncated.

(Abridged) We measure the mass distribution of galaxy cluster Abell 1689 within 0.3 Mpc/h_70 of the cluster centre using its strong lensing effect on 32 background galaxies. The multiple images are based on those of Broadhurst et al. 2005 with some modifications. The cluster profile is explored further out to ~2.5 Mpc/h_70 with weak lensing shear measurements from Broadhurst et al. 2005b. The masses of ~200 cluster galaxies are measured with Fundamental Plane in order to accurately model the small scale mass structure in the cluster. The galaxies are modelled as elliptical truncated isothermal spheres. The dark matter component of the cluster is described by either non-singular isothermal ellipsoids (NSIE) or elliptical versions of the universal dark matter profile (ENFW). We use two dark matter haloes to model the smooth DM in the cluster. The total mass profile is well described by either an NSIS profile with sigma=1514+-18 km/s and core radius of r_c=71+-5kpc/h_70, or an NFW profile with C=6.0+-0.5 and r_200=2.82+-0.11 Mpc/h_70. The errors are assumed to be due to the error in assigning masses to the individual galaxies in the galaxy component. The derived total mass is in good agreement with the mass profile of Broadhurst et al. 05. Using also weak lensing we can constrain the profile further out to r~2.5 Mpc/h_70. The best fit parameters are then sigma=1499+-15 km/s and r_c=66+-5 kpc/h_70 for the NSIS profile and C=7.6+-0.5 and r_200=2.55+-0.07 Mpc/h_70 for the NFW profile. Using the same image configuration as Broadhurst et al. 2005 we obtain a strong lensing model that is superior to that of Broadhurst et al. 2005 (rms of 2.7'' compared to 3.2'').

We study the evolution of the stellar mass density for the separate families of bulge-dominated and disk-dominated galaxies over the redshift range 0.25 < z < 1.15. We derive quantitative morphology for a statistically significant galaxy sample of 1645 objects selected from the FORS Deep and the GOODS-S Fields. We find that the morphological mix evolves monotonically with time: the higher the redshift, the more disk systems dominate the total mass content. At redshift about 1, massive objects (M_stellar > 7E10 M_solar) host about half of the mass contained in objects of similar mass in the local universe. The contribution from early and late type galaxies to the mass budget at z about 1 is nearly equal. We show that in situ star formation is not sufficient to explain the changing mass budget. Moreover we find that the star formation rate per unit stellar mass of massive galaxies increases with redshift only for the intermediate and early morphological types, while it stays nearly constant for late-type objects. This suggests that merging and/or frequent accretion of small mass objects has a key role in the shaping of the Hubble sequence as we observe it now, and also in decreasing the star formation activity of the bulge-dominated descendants of massive disk galaxies.

We present a measurement of the evolution of the stellar mass function (MF) of galaxies and the evolution of the total stellar mass density at 0<z<5. We use deep multicolor data in the Fors Deep Field (FDF; I-selected reaching I_AB=26.8) and the GOODS-S/CDFS region (K-selected reaching K_AB=25.4) to estimate stellar masses based on fits to composite stellar population models for 5557 and 3367 sources, respectively. The MF of objects from the GOODS-S sample is very similar to that of the FDF. Near-IR selected surveys hence detect the more massive objects of the same principal population as do I-selected surveys. We find that the most massive galaxies harbor the oldest stellar populations at all redshifts. At low z, our MF follows the local MF very well, extending the local MF down to 10^8 Msun. The faint end slope is consistent with the local value of alpha~1.1 at least up to z~1.5. Our MF also agrees very well with the MUNICS and K20 results at z<2. The MF seems to evolve in a regular way at least up to z~2 with the normalization decreasing by 50% to z=1 and by 70% to z=2. Objects having M>10^10 Msun which are the likely progenitors of todays L* galaxies are found in much smaller numbers above z=2. However, we note that massive galaxies with M>10^11 Msun are present even to the largest redshift we probe. Beyond z=2 the evolution of the mass function becomes more rapid. We find that the total stellar mass density at z=1 is 50% of the local value. At z=2, 25% of the local mass density is assembled, and at z=3 and z=5 we find that at least 15% and 5% of the mass in stars is in place, respectively. The number density of galaxies with M>10^11 Msun evolves very similarly to the evolution at lower masses. It decreases by 0.4 dex to z=1, by 0.6 dex to z=2, and by 1 dex to z=4.

We measure the star formation rate (SFR) as a function of redshift z up to z \~4.5, based on B, I and (I+B) selected galaxy catalogues from the FORS Deep Field (FDF) and the K-selected catalogue from the GOODS-South field. Distances are computed from spectroscopically calibrated photometric redshifts accurate to (Delta_z / (z_spec+1)) ~0.03 for the FDF and ~0.056 for the GOODS-South field. The SFRs are derived from the luminosities at 1500 Angstroem. We find that the total SFR estimates derived from B, I and I+B catalogues agree very well ($\lsim 0.1$ dex) while the SFR from the K catalogue is lower by ~0.2 dex. We show that the latter is solely due to the lower star-forming activity of K-selected intermediate and low luminosity (L<L_*) galaxies. The SFR of bright (L>L_*) galaxies is independent of the selection band, i.e. the same for B, I, (I+B), and K-selected galaxy samples. At all redshifts, luminous galaxies (L>L_*) contribute only ~1/3 to the total SFR. There is no evidence for significant cosmic variance between the SFRs in the FDF and GOODs-South field, ~0.1 dex, consistent with theoretical expectations. The SFRs derived here are in excellent agreement with previous measurements provided we assume the same faint-end slope of the luminosity function as previous works (alpha ~ -1.6). However, our deep FDF data indicate a shallower slope of alpha=-1.07, implying a SFR lower by ~0.3 dex. We find the SFR to be roughly constant up to z ~4 and then to decline slowly beyond, if dust extinctions are assumed to be constant with redshift.

We present the Capodimonte Deep Field (OACDF), a deep field covering an area of 0.5 square degrees in the B, V, R optical bands plus six medium-band filters in the wavelength range 773--913 nm. The field reaches the following limiting AB magnitudes: B~25.3, V~24.8 and R~25.1 and contains about 50000 extended sources in the AB magnitude range 18 < R < 25 . Hence, it is intermediate between deep pencil beam surveys and very wide but shallow surveys. The main scientific goal of the OACDF is the identification and characterization of early-type field galaxies at different look-back times in order to study different scenarios of galaxy formation. Parallel goals include the search for groups and clusters of galaxies and the search for rare and peculiar objects (gravitational lenses, QSOs, halo White Dwarfs). In this paper we describe the OACDF data reduction, the methods adopted for the extraction of the photometric catalogs, the photometric calibration and the quality assessment of the catalogs by means of galaxy number counts, spectroscopic and photometric redshifts and star colors. We also present the first results of the search for galaxy overdensities. The depth of the OACDF and its relatively large spatial coverage with respect to pencil beam surveys make it a good tool for further studies of galaxy formation and evolution in the redshift range 0--1, as well as for stellar studies.

We use a high resolution collisionless simulation of a Virgo--like cluster in a $\Lambda$CDM cosmology to determine the velocity and clustering properties of the diffuse stellar component in the intracluster region at the present epoch. The simulated cluster builds up hierarchically and tidal interactions between member galaxies and the cluster potential produce a diffuse stellar component free-flying in the intracluster medium. Here we adopt an empirical scheme to identify tracers of the stellar component in the simulation and hence study its properties. We find that at $z=0$ the intracluster stellar light is mostly unrelaxed in velocity space and clustered in structures whose typical clustering radii are about 50 kpc at R=400--500 kpc from the cluster center, and predict the radial velocity distribution expected in spectroscopic follow-up surveys. Finally, we compare the spatial clustering in the simulation with the properties of the Virgo intracluster stellar population, as traced by ongoing intracluster planetary nebulae surveys in Virgo. The preliminary results indicate a substantial agreement with the observed clustering properties of the diffuse stellar population in Virgo.

(Abridged) We report on observations of the optical and NIR afterglow of GRB020405. Ground-based optical observations started about 1 day after the GRB and spanned a period of ~10 days; archival HST data extended the coverage up to 70 days after the GRB. We report the first detection of the afterglow in NIR bands. The detection of emission lines in the optical spectrum indicates that the GRB is located at z = 0.691. Absorptions are also detected at z = 0.691 and at z = 0.472. The latter system is likely caused by clouds in a galaxy located 2 arcsec southwest of the GRB host. Hence, for the first time, the galaxy responsible for an intervening absorption system in the spectrum of a GRB afterglow is identified. Optical and NIR photometry indicates that the decay in all bands follows a single power law of index alpha = 1.54. The late-epoch VLT and HST points lie above the extrapolation of this power law, so that a plateau is apparent in the VRIJ light curves at 10-20 days after the GRB. The light curves at epochs later than day ~20 after the GRB are consistent with a power-law decay with index alphaprime = 1.85. We suggest that this deviation can be modeled with a SN having the same temporal profile as SN2002ap, but 1.3 mag brighter at peak, and located at the GRB redshift. Alternatively, a shock re-energization may be responsible for the rebrightening. A polarimetric R-band measurement shows that the afterglow is polarized, with P = 1.5 % and theta = 172 degrees. Optical-NIR spectral flux distributions show a change of slope across the J band which we interpret as due to the presence of nu_c. The analysis of the multiwavelength spectrum within the fireball model suggests that a population of relativistic electrons produces the optical-NIR emission via synchrotron in an adiabatically expanding blastwave, and the X-rays via IC.

CSL-1 (Capodimonte--Sternberg--Lens Candidate, No.1) is an extragalactic double source detected in the OAC-DF (Osservatorio Astronomico di Capodimonte - Deep Field). It can be interpreted either as the chance alignment of two identical galaxies at z=0.46 or as the first case of gravitational lensing by a cosmic string. Extensive modeling shows in fact that cosmic strings are the only type of lens which (at least at low angular resolution) can produce undistorted double images of a background source. We propose an experimentum crucis to disentangle between these two possible explanations. If the lensing by a cosmic string should be confirmed, it would provide the first measurements of energy scale of symmetry breaking and of the energy scale of Grand Unified Theory (GUT).

We have identified intracluster planetary nebulae in a Virgo cluster core field by imaging with the Subaru Suprime-Cam through two narrow band filters centered at the redshifted wavelengths of the [OIII] $\lambda=5007$ \AA and the H$\alpha$ $\lambda = 6563$ \AA lines; broad-band images in V and R were acquired to check for emission in the adjacent continuum. Emission line objects in Virgo are then selected from the two-color diagram ([OIII] -- Halpha) vs. ([OIII] -- (V+R)), which is calibrated using PNe in M84 (Jacoby et al. 1990). Using both [OIII] and Halpha allows us to distinguish bona-fide planetary nebulae from high redshift emission-line galaxies at the bright end of the [OIII] luminosity function. Spectroscopic observations of a subsample of these objects were made at the TNG and at the VLT, in a region around M84 and in an intracluster field respectively. The observations confirm the efficiency of the combined [OIII]+H$\alpha$ imaging to identify true PNe. We also obtained the first spectrum of an intracluster PN which shows the [OIII] doublet with S/N > 10 and its Halpha emission. From the results based on the spectroscopic follow-up, we derive a lower limit to the fraction of the Virgo cluster light contributed by the intracluster stars at the surveyed position in the cluster core: it amounts to 10%.

We have identified 38 candidates of intracluster planetary nebulae (ICPNe) in a $34'\times27'$ field in the core of the Virgo cluster based on the Suprime-Cam imaging through two narrow-band filters centered at the redshifted wavelengths of the [OIII]$\lambda=5007\AA$ and the Halpha \lambda=6563\AA lines. Broad-band images in V and R bands are used to check for any emissions in the adjacent continuum. We describe the method briefly and present the list of intracluster planetary nebulae candidates, together with their finding charts. The ICPN candidates show a highly inhomogeneous distribution, which may suggest an association with the M86-M84 subcluster. Fraction of diffuse intracluster light with respect to total light in galaxies is estimated to be about 10%, leading to an estimate of about 20% for the baryon fraction. Spectroscopic follow up and a wider survey are critical to reveal the nature of intracluster stellar population.

The Capodimonte Deep Field (OACDF) is a multi-colour imaging survey on two 0.5x0.5 square degree fields performed in the BVRI bands and in six medium-band filters (700 - 900 nm) with the Wide Field Imager (WFI) at the ESO 2.2 m telescope at La Silla, Chile. In this contribution the adopted strategies for the OACDF data reduction are discussed. Preliminary scientific results of the survey are also presented.