The Lya line in the UV and the [CII] line in the FIR are widely used tools to identify galaxies and to obtain insights into ISM properties in the early Universe. By combining data obtained with ALMA in band 7 at ~ 320 GHz as part of the ALMA Large Program to INvestigate [CII] at Early Times (ALPINE) with spectroscopic data from DEIMOS at Keck, VIMOS and FORS2 at the VLT, we assembled a unique sample of 53 main-sequence star-forming galaxies at 4.4 < z < 6 in which we detect both the Lya line and the [CII]. We used [CII], observed with ALMA, as a tracer of the systemic velocity of the galaxies, and we find that 90% of the selected objects have Lya-[CII] velocity offsets in the range 0 < Dv_Lya-[CII] < 400 km/s, in line with the few measurements available so far in the early Universe, and significantly smaller than those observed at lower z. We observe ISM-[CII] offsets in the range -500 < Dv_ISM-[CII] < 0 km/s, in line with values at all redshifts. We find significant anticorrelations between Dv_Lya-[CII] and the Lya rest-frame equivalent width EW0(Lya) (or equivalently, the Lya escape fraction f_esc(Lya)). According to available models for the radiative transfer of Lya photons, the escape of Lya photons would be favored in galaxies with high outflow velocities, in agreement with our observations. The uniform shell model would also predict that the Lya escape in galaxies with slow outflows (0 < v_out < 300 km/s) is mainly determined by the neutral hydrogen column density (NHI), while the alternative model by Steidel+10 would favor a combination of NHI and covering fraction as driver of the Lya escape. We suggest that the observed increase in Lya escape that is observed between z~2 and z~6 is not due to a higher incidence of fast outflows at high redshift, but rather to a decrease in average NHI along the line of sight, or alternatively, a decrease in HI covering fraction. [abridged]

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) observations towards 27 low-redshift ($0.02< z<0.2$) star-forming galaxies taken from the Valpara\'iso ALMA/APEX Line Emission Survey (VALES). We perform stacking analyses of the $^{12}$CO($1-0$), $^{13}$CO($1-0$) and C$^{18}$O($1-0$) emission lines to explore the $L'$ ($^{12}$CO($1-0$))/$L'$($^{13}$CO($1-0$))) (hereafter $L'$($^{12}$CO)/$L'$($^{13}$CO)) and $L'$($^{13}$CO($1-0$))/$L'$(C$^{18}$O($1-0$)) (hereafter $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratio dependence as a function of different global galaxy parameters related to the star formation activity. The sample has far-IR luminosities $10^{10.1-11.9}$L$_{\odot}$ and stellar masses of $10^{9.8-10.9}$M$_{\odot}$ corresponding to typical star-forming and starburst galaxies at these redshifts. On average we find a $L'$($^{12}$CO)/$L'$($^{13}$CO) line luminosity ratio value of 16.1$\pm$2.5. Galaxies with evidences of possible merging activity tend to show higher $L'$($^{12}$CO)/$L'$($^{13}$CO) ratios by a factor of two, while variations of this order are also found in galaxy samples with higher star formation rates or star formation efficiencies. We also find an average $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratio of 2.5$\pm$0.6, which is in good agreement with those previously reported for starburst galaxies. We find that galaxy samples with high $L_{\text{IR}}$, SFR and SFE show low $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratios with high $L'$($^{12}$CO)/$L'$($^{13}$CO) line luminosity ratios, suggesting that these trends are produced by selective enrichment of massive stars in young starbursts.

The Lyman-$\alpha$ (Ly$\alpha$) emission line has been ubiquitously used to confirm and study high redshift galaxies. We report on the line morphology as seen in the 2D spectra from the VIMOS Ultra Deep Survey in a sample of 914 Ly$\alpha$ emitters from a parent sample of 4192 star-forming galaxies at $2<z_\mathrm{spec}\lesssim6$. The study of the spatial extent of Ly$\alpha$ emission provides insight into the escape of Ly$\alpha$ photons from galaxies. We classify the line emission as either non-existent, coincident, projected spatial offset, or extended with respect to the observed 2D UV continuum emission. The line emitters in our sample are classified as ~45% coincident, ~24% extended and ~11% offset emitters. For galaxies with detected UV continuum, we show that extended Ly$\alpha$ emitters (LAEs) correspond to the highest equivalent width galaxies (with an average $W_\mathrm{Ly\alpha}\sim-22${\AA}). This means that this class of objects is the most common in narrow-band selected samples, which usually select high equivalent width LAEs, $<-20${\AA}. Extended Ly$\alpha$ emitters are found to be less massive, less star-forming, with lower dust content, and smaller UV continuum sizes ($r_{50}\sim0.9$kpc) of all the classes considered here. We also find that galaxies with larger UV-sizes have lower fractions of Ly$\alpha$ emitters. By stacking the spectra per emitter class we find that the weaker Ly$\alpha$ emitters have stronger low ionization inter-stellar medium (ISM) absorption lines. Interestingly, we find that galaxies with Ly$\alpha$ offset emission (median separation of $1.1_{-0.8}^{+1.3}$kpc from UV continuum) show similar velocity offsets in the ISM as those with no visible emission (and different from other Ly$\alpha$ emitting classes). This class of objects may hint at episodes of gas accretion, bright offset clumps, or on-going merging activity into the larger galaxies.

In the second work of this series, we analyse the connection between the availability of gas and the position of a region with respect to the spatially resolved main sequence (MS) relation. Following the procedure presented in Paper I we obtain 500pc scales estimates of stellar mass and star formation rate surface densities ($\Sigma_{\star}$ and $\Sigma_{\rm{SFR}}$). Our sample consists of five face-on, grand design spiral galaxies located on the MS. Thanks to HI 21cm and $^{12}$CO(2-1) maps, we connect the gas surface densities and gas fractions to the observed star formation properties of each region. We find that the spatially resolved MS ($\sigma=0.23$ dex) is the combination of two relations: the Kennicutt-Schmidt law ($\sigma=0.19$ dex) and the molecular gas MS (MGMS, $\sigma=0.22$ dex); $\Sigma_{\star}$, $\Sigma_{\rm{SFR}}$ and the surface density of the molecular gas, $\Sigma_{\rm{H_2}}$, define a 3D relation as proposed by \citet{2019ApJ...884L..33L}. We find that $\Sigma_{\rm{H_2}}$ steadily increases along the MS relation, varies little towards higher $\Sigma_{\rm{SFR}}$ at fixed stellar surface densities (not enough to sustain the change in SFR), and it is almost constant perpendicular to the relation. The surface density of neutral gas ($\Sigma_{\rm{HI}}$) is constant along the MS, and increases in its upper envelop. $\Sigma_{\rm{SFR}}$ can be expressed as a function of $\Sigma_{\star}$ and $\Sigma_{\rm{HI}}$, following the Equation: $\log\Sigma_{\rm{SFR}}$ = 0.97$\log\Sigma_{\star}$ + 1.99$\log\Sigma_{\rm{HI}}$ - 11.11. Finally, we show that f$_{\rm{gas}}$ increases significantly towards the starburst region in the $\log\Sigma_{\star}$ - $\log\Sigma_{\rm{SFR}}$ plane, accompanied by a slight increase in SFE.

The ALMA-ALPINE [CII] survey is aimed at characterizing the properties of a sample of normal star-forming galaxies (SFGs). The ALMA Large Program to INvestigate (ALPINE) features 118 galaxies observed in the [CII]-158$\mu$m line and far infrared (FIR) continuum emission during the period of rapid mass assembly, right after the end of the HI reionization, at redshifts of 4<z<6. We present the survey science goals, the observational strategy, and the sample selection of the 118 galaxies observed with ALMA, with an average beam minor axis of about 0.85 arcsec, or $\sim$5 kpc at the median redshift of the survey. The properties of the sample are described, including spectroscopic redshifts derived from the UV-rest frame, stellar masses, and star-formation rates obtained from a spectral energy distribution (SED) fitting. The observed properties derived from the ALMA data are presented and discussed in terms of the overall detection rate in [CII] and FIR continuum, with the observed signal-to-noise distribution. The sample is representative of the SFG population in the main sequence at these redshifts. The overall detection rate in [CII] is 64% for a signal-to-noise ratio (S/N) threshold larger than 3.5 corresponding to a 95% purity (40% detection rate for S/N>5). Based on a visual inspection of the [CII] data cubes together with the large wealth of ancillary data, we find a surprisingly wide range of galaxy types, including 40% that are mergers, 20% extended and dispersion-dominated, 13% compact, and 11% rotating discs, with the remaining 16% too faint to be classified. This diversity indicates that a wide array of physical processes must be at work at this epoch, first and foremost, those of galaxy mergers. This paper sets a reference sample for the gas distribution in normal SFGs at 4<z<6.

We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA Band-7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 16. per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (7 show a counterpart only in the mid-IR and no HST or near-IR identification, while 2 are detected as [CII] emitters at z=5). The 7 HST and near-IR dark galaxies with mid-IR counterpart contribute for about 15 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3.

We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim 4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A} < -20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [C II] detections and 43 are upper limits. With this UV-selected sample, we set stringent lower bounds to the volume density of [C II] emitters at $z\sim4-6$. Our derived LFs are consistent with the $z\sim0$ [C II] LF at $\leq 10^9L_\odot$, but show a hint of an excess of [C II] emitters at the luminous end ($>10^9L_\odot$) compared to $z\sim0$. This excess is consistent with a potential population of [C II]-bright UV-faint sources that are largely excluded from the ALPINE target selection. To further constrain such sources, in a companion analysis, Loiacono et al. (2020) estimate the [C II] LF based on serendipitous [C II] sources in the ALPINE maps. Their resulting single data point is consistent with the excess we observe at high-L. Combined, these studies suggest possible redshift evolution in the [C II] LFs at high luminosity between $z\sim 5$ and $0$. This conclusion is further supported by the inferred [C II] LFs from the far-IR and CO LFs at $z\sim4-6$. Integrating our estimate on the best-fit [C II] LF, we set a constraint on molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$ Mpc$^{-3}$. This is broadly consistent with previous studies, although the uncertainties span more than an order of magnitude. Finally, we find that available model predictions tend to underestimate the number densities of [C II] emitters at $z\sim5$.

We present the first [CII] 158 $\mu$m luminosity function (LF) at $z\sim 5$ from a sample of serendipitous lines detected in the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). A search performed over the 118 ALPINE pointings revealed several serendipitous lines. Based on their fidelity, we selected 14 lines for the final catalog. According to the redshift of their counterparts, we identified 8 out of 14 detections as [CII] lines at $z\sim 5$, and two as CO transitions at lower redshifts. The remaining 4 lines have an elusive identification in the available catalogs and we considered them as [CII] candidates. We used the 8 confirmed [CII] and the 4 [CII] candidates to build one of the first [CII] LFs at $z\sim 5$. We found that 11 out of these 12 sources have a redshift very similar to that of the ALPINE target in the same pointing, suggesting the presence of overdensities around the targets. Therefore, we split the sample in two (a "clustered" and "field" sub-sample) according to their redshift separation and built two separate LFs. Our estimates suggest that there could be an evolution of the [CII] LF between $z \sim 5$ and $z \sim 0$. By converting the [CII] luminosity to star formation rate we evaluated the cosmic star formation rate density (SFRD) at $z\sim 5$. The clustered sample results in a SFRD $\sim 10$ times higher than previous measurements from UV-selected galaxies. On the other hand, from the field sample (likely representing the average galaxy population) we derived a SFRD $\sim 1.6$ higher compared to current estimates from UV surveys but compatible within the errors. Because of the large uncertainties, observations of larger samples are necessary to better constrain the SFRD at $z\sim 5$. This study represents one of the first efforts aimed at characterizing the demography of [CII] emitters at $z\sim 5$ using a mm-selection of galaxies.

The ALPINE-ALMA large program targets the [CII] 158 $\mu$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates, we measured the conversion factor from the ALMA 158 $\mu$m rest-frame dust continuum luminosity to the total infrared luminosity (L$_{\rm IR}$) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L$_{\rm IR}$ of 4.4$\times 10^{11}$ L$_\odot$. We also detected 57 additional continuum sources in our ALMA pointings. They are at lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5$\pm$0.2. We measured the 850 $\mu$m number counts between 0.35 and 3.5 mJy, improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40 % of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8$\times$10$^8$ L$_\odot$ and their median full width at half maximum is 252 km/s. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR-L$_{\rm [CII]}$ relations of De Looze et al. (2014) and Lagache et al. (2018).

The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature we examine the relation between the [CII] luminosity and the SFR over the entire redshift range from $z \sim 4-8$. ALPINE galaxies, which are both detected in [CII] and dust continuum, show a good agreement with the local L([CII])-SFR relation. Galaxies undetected in the continuum with ALMA are found to be over-luminous in [CII], when the UV SFR is used. After accounting for dust-obscured star formation, by an amount SFR(IR)$\approx$SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L([CII])-SFR relation comparable to the local one. When [CII] non-detections are taken into account, the slope may be marginally steeper at high-z, although this is still somewhat uncertain. When compared in a homogeneous manner, the $z>6 $ [CII] measurements (detections and upper limits) do not behave very differently from the $z \sim 4-6$ data. We find a weak dependence of L([CII])/SFR on the Lyman-alpha equivalent width. Finally, we find that the ratio L([CII])/LIR $\sim (1-3) \times 10^{-3}$ for the ALPINE sources, comparable to that of "normal" galaxies at lower redshift. Our analysis, which includes the largest sample ($\sim 150$ galaxies) of [CII] measurements at $z>4$ available so far, suggests no or little evolution of the L([CII])-SFR relation over the last 13 Gyr of cosmic time.

We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope images. The effective radius r_e,[CII] exceeds r_e,UV by factors of ~2-3 and the ratio of r_e,[CII]/r_e,UV increases as a function of M_star. We do not find strong evidence that [CII] line, the rest-frame UV, and FIR continuum are always displaced over ~ 1-kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [CII] component over 10-kpc scales detected at 4.1$\sigma$-10.9$\sigma$ beyond the size of rest-frame UV and far-infrared (FIR) continuum. One object has tentative rotating features up to ~10-kpc, where the 3D model fit shows the rotating [CII]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [CII] line structure have high star-formation rate (SFR), stellar mass (M_star), low Lya equivalent-width, and more blue-shifted (red-shifted) rest-frame UV metal absorption (Lya line), as compared to galaxies without such extended [CII] structures. Although we cannot rule out the possibility that a selection bias towards luminous objects may be responsible for such trends, the star-formation driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [CII] line structures are ubiquitous to high-z star-forming galaxies.

We present ALMA observations of a merging system at z ~ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII] 158 micron and far-infrared continuum data with multi-wavelength ancillary data we find that the system is composed of two massive (Mstar >~ 10^10 Msun) star-forming galaxies experiencing a major merger (stellar mass ratio r_mass ~ 0.9) at close spatial (~13 kpc; projected) and velocity (delta_v < 300 km/s) separations, and two additional faint narrow [CII]-emitting satellites. The overall system belongs to a larger-scale protocluster environment and is coincident to one of its overdensity peaks. ALMA reveals also the presence of [CII] emission arising from a circumgalactic gas structure, extending out to ~30 kpc. Our morpho-spectral decomposition analysis shows that about 50% of the total flux resides between the individual galaxy components, in a metal-enriched gaseous envelope characterized by a disturbed morphology and complex kinematics. Similarly to observations of shock-excited [CII] emitted from tidal tails in local groups, our results can be interpreted as a possible signature of interstellar gas stripped by strong gravitational interactions, with a possible contribution from material ejected by galactic outflows and emission triggered by star formation in small faint satellites. Our findings suggest that mergers could be an efficient mechanism of gas mixing in the circumgalactic medium around high-z galaxies, and thus play a key role in the galaxy baryon cycle at early epochs.

The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z<2. We discuss three possible effects, namely outflows, halt of gas supplying, and over-efficient star formation, which may jointly contribute to the gas fraction plateau of the latter massive galaxies.

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at redshift ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame $158\,\rm{\mu m}$ obtained with the ALMA large program ALPINE. We study the connection between the UV spectral slope ($\beta$), stellar mass ($M_{\star}$), and infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$). Twenty-three galaxies are individually detected in the continuum at $>3.5\,\sigma$ significance. We perform a stacking analysis using both detections and non-detections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$\beta$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$\beta$ relation as a function of redshift suggests an evolution of dust attenuation properties at $z>4$. Similarly, we find that our galaxies have lower IRX values up to 1 dex on average at fixed mass compared to previously studied IRX-$M_{\star}$ relations at $z\lesssim4$, albeit with significant scatter. This implies a lower obscured fraction of star-formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at $z\gtrsim4$ are characterised by (i) a steeper attenuation curve than at $z\lesssim4$, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies ($\log M_{\star}/M_\odot > 10$) at z~5-6 already exhibit an obscured fraction of star formation of $\sim45\%$, indicating a rapid build-up of dust during the epoch of reionization.

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$\sigma$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($\beta$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($\Delta\beta$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$\lambda$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows.

We analyse the spatially resolved relation between stellar mass (M$_{\star}$) and star formation rate (SFR) in disk galaxies (i.e. the Main Sequence, MS). The studied sample includes eight nearby face-on grand-design spirals, e.g. the descendant of high-redshift, rotationally-supported star-forming galaxies. We exploit photometric information over 23 bands, from the UV to the far-IR, from the publicly available DustPedia database to build spatially resolved maps of stellar mass and star formation rates on sub-galactic scales of 0.5-1.5 kpc, by performing a spectral energy distribution fitting procedure that accounts for both the observed and the obscured star formation processes, over a wide range of internal galaxy environments (bulges, spiral arms, outskirts). With more than 30 thousands physical cells, we have derived a definition of the local spatially resolved MS per unit area for disks, $\log(\Sigma_{SFR})$=0.82log$(\Sigma_{*})$-8.69. This is consistent with the bulk of recent results based on optical IFU, using the H$\alpha$ line emission as a SFR tracer. Our work extends the analysis at lower sensitivities in both M$_{\star}$ and SFR surface densities, up to a factor $\sim$ 10. The self consistency of the MS relation over different spatial scales, from sub-galactic to galactic, as well as with a rescaled correlation obtained for high redshift galaxies, clearly proves its universality.

We present the detection of CO(5-4) with S/N> 7 - 13 and a lower CO transition with S/N > 3 (CO(4-3) for 4 galaxies, and CO(3-2) for one) with ALMA in band 3 and 4 in five main sequence star-forming galaxies with stellar masses 3-6x10^10 M/M_sun at 3 < z < 3.5. We find a good correlation between the total far-infrared luminosity LFIR and the luminosity of the CO(5-4) transition L'CO(5-4), where L'CO(5-4) increases with SFR, indicating that CO(5-4) is a good tracer of the obscured SFR in these galaxies. The two galaxies that lie closer to the star-forming main sequence have CO SLED slopes that are comparable to other star-forming populations, such as local SMGs and BzK star-forming galaxies; the three objects with higher specific star formation rates (sSFR) have far steeper CO SLEDs, which possibly indicates a more concentrated episode of star formation. By exploiting the CO SLED slopes to extrapolate the luminosity of the CO(1-0) transition, and using a classical conversion factor for main sequence galaxies of alpha_CO = 3.8 M_sun(K km s^-1 pc^-2)^-1, we find that these galaxies are very gas rich, with molecular gas fractions between 60 and 80%, and quite long depletion times, between 0.2 and 1 Gyr. Finally, we obtain dynamical masses that are comparable with the sum of stellar and gas mass (at least for four out of five galaxies), allowing us to put a first constraint on the alpha_CO parameter for main sequence galaxies at an unprecedented redshift.

We report the serendipitous discovery of a bright galaxy (Gal-A) observed as part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). While this galaxy is detected both in line and continuum emission in ALMA Band 7, it is completely dark in UV/optical filters and only presents a marginal detection in the UltraVISTA Ks band. We discuss the nature of the observed ALMA line, i.e. whether the emission comes from [CII] at z~4.6, or from high-J CO transitions at z~2.2. In the first case we find a [CII]-to-FIR luminosity ratio of log(L_[CII]/L_FIR)=-2.5, consistent with the average value for local star-forming galaxies (SFGs); in the second case, instead, the source would lie outside of the empirical relations between L_CO and L_FIR found in the literature. At both redshifts, we derive the star-formation rate (SFR) from the ALMA continuum, and the stellar mass (M*) by using stellar population synthesis models as input for LePHARE spectral energy distribution (SED) fitting. Exploiting our results, we believe that Gal-A is a Main-Sequence (MS), dusty SFG at z=4.6 (i.e. [CII] emitter) with log(SFR/[M/yr])~1.4 and log(M*/M)~9.7. This work underlines the crucial role of the ALPINE survey in making a census of this class of objects, in order to unveil their contribution to the global star-formation rate density (SFRD) of the Universe at the end of the Reionisation epoch.

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

Aims. We seek is to identify old and massive galaxies at 0.5<z<2.1 on the basis of the magnesium index MgUV and then study their physical properties. We computed the MgUV index based on the best spectral fitting template of $\sim$3700 galaxies using data from the VLT VIMOS Deep Survey (VVDS) and VIMOS Ultra Deep Survey (VUDS) galaxy redshift surveys. Based on galaxies with the largest signal to noise and the best fit spectra we selected 103 objects with the highest spectral MgUV signature. We performed an independent fit of the photometric data of these galaxies and computed their stellar masses, star formation rates, extinction by dust and age, and we related these quantities to the MgUV index. We find that the MgUV index is a suitable tracer of early-type galaxies at an advanced stage of evolution. Selecting galaxies with the highest MgUV index allows us to choose the most massive, passive, and oldest galaxies at any epoch. The formation epoch t_f computed from the fitted age as a function of the total mass in stars supports the downsizing formation paradigm in which galaxies with the highest mass formed most of their stars at an earlier epoch.

We report the detection of [CII]158um emission from a system of three closely-separated sources in the COSMOS field at z~4.56, as part of the ALMA Large Program to INvestigate CII at Early times (ALPINE). The two dominant sources are closely associated, both spatially (1.6"~11kpc) and in velocity (~100km/s), while the third source is slightly more distant (2.8"~18kpc, ~300km/s). The second strongest source features a slight velocity gradient, while no significant velocity gradient is seen in the other two sources. Using the observed [CII] luminosities, we derive a total log(SFR_[CII]/[Msol/year])=2.8+/-0.2, which may be split into contributions of 59%, 31%, and 10% from the central, east, and west sources, respectively. Comparison of these [CII] detections to recent zoom-in cosmological simulations suggests an ongoing major merger. We are thus witnessing a system in a major phase of mass build-up by merging, including an on-going major merger and an upcoming minor merger, which is expected to end up in a single massive galaxy by z~2.5.

Based on HST and MUSE data, we probe the stellar and gas properties (i.e. kinematics, stellar mass, star formation rate) of the radio-loud brightest cluster galaxy (BCG) located at the centre of the X-ray luminous cool core cluster Abell 2667 (z = 0.2343). The bi-dimensional modelling of the BCG surface brightness profile reveals the presence of a complex system of substructures extending all around the galaxy. Clumps of different size and shape plunged into a more diffuse component constitute these substructures, whose intense 'blue' optical colour hints to the presence of a young stellar population. Our results depict the BCG as a massive (M_star ~ 1.38 x 10^11 M_sun) dispersion-supported spheroid (v_star < 150 km/s, sigma_0 ~ 216 km/s) hosting an active supermassive black hole (M_SMBH ~ 3.8 x 10^9 M_sun) whose optical features are typical of low ionisation nuclear emission line regions. Although the velocity pattern of the stars in the BCG is irregular, the stellar kinematics in the regions of the clumps show a positive velocity of ~ 100 km/s, similarly to the gas component. An analysis of the mechanism giving rise to the observed lines in the clumps through empirical diagnostic diagrams points out that the emission is composite, suggesting the contribution from both star formation and AGN. We conclude our analysis describing how scenarios of both chaotic cold accretion and merging with a gas-rich disc galaxy can efficaciously explain the phenomena the BCG is undergoing.

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

We establish a robust statistical description of the star-forming galaxy population at the end of cosmic HI reionization ($5.0\le{}z\le6.6$) from a large sample of 52 galaxies with spectroscopically confirmed redshifts from the VIMOS UltraDeep Survey. We identify galaxies with Ly$\alpha$ either in absorption or in emission, at variance with most spectroscopic samples in the literature where Ly$\alpha$ emitters dominate. We find that star-forming galaxies at these redshifts are distributed along a main sequence in the stellar mass vs. SFR plane. We report a flat evolution of the sSFR(z) in 3<z<6 compared to lower redshift measurements. UV-continuum slopes vary with luminosity, with a large dispersion. We determing UV and Ly$\alpha$ luminosity functions using V$_{max}$ method and use them to derive star formation rate densities (SFRD). We find that both UV-derived and Ly$\alpha$-derived SFRDs are in excellent agreement after correcting Ly$\alpha$ luminosity density for IGM absorption. Our new SFRD measurements at a mean redshift z=5.6 confirm the steep decline of the SFRD at z>2. The bright end of the Ly$\alpha$ luminosity function has a high number density, indicating a significant star formation activity concentrated in the brightest Ly$\alpha$ emitters (LAE) at these redshifts. LAE with EW>25\AA ~contribute to about 75\% of the total UV-derived SFRD. While our analysis favors a low dust content in 5.0<z<6.6, uncertainties on the dust extinction correction and associated degeneracies in spectral fitting will remain an issue to estimate the total SFRD until future survey extending spectroscopy to the NIR rest-frame spectral domain, e.g. with JWST.

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

Thanks to deep optical to near-IR imaging and spectroscopy, significant progress is made in characterizing the rest-frame UV to optical properties of galaxies in the early universe (z > 4). Surveys with Hubble, Spitzer, and ground-based facilities (Keck, Subaru, and VLT) provide spectroscopic and photometric redshifts, measurements of the spatial structure, stellar masses, and optical emission lines for large samples of galaxies. Recently, the Atacama Large (Sub) Millimeter Array (ALMA) has become a major player in pushing studies of high redshift galaxies to far-infrared wavelengths, hence making panchromatic surveys over many orders of frequencies possible. While past studies focused mostly on bright sub-millimeter galaxies, the sensitivity of ALMA now enables surveys like ALPINE, which focuses on measuring the gas and dust properties of a large sample of normal main-sequence galaxies at z > 4. Combining observations across different wavelengths into a single, panchromatic picture of galaxy formation and evolution is currently and in the future an important focus of the astronomical community.

[Abridged] We unveil the complex shape of a proto-supercluster at z~2.45 in the COSMOS field using the spectroscopic redshifts of the VIMOS Ultra-Deep Survey (VUDS), complemented by the zCOSMOS-Deep sample and high-quality photometric redshifts. The method relies on a 2D Voronoi tessellation in overlapping redshift slices that is converted into a 3D density field. The galaxy distribution is constructed using a statistical treatment of spectroscopic and photometric redshifts. We identify a proto-supercluster, dubbed "Hyperion" for its immense size and mass, which extends over a volume of ~60x60x150 comoving Mpc^3 and has an estimated total mass of ~4.8x 10^15 M_sun. This immensely complex structure contains at least 7 density peaks within 2.4 < z < 2.5, connected by filaments. Based on the peaks average matter density, we estimate their total mass, M_tot, and find a range of masses from ~0.1x10^14 M_sun to ~2.7x10^14 M_sun. By using spectroscopic members of each peak, we obtain the velocity dispersion of the galaxies in the peaks, and then their virial mass M_vir (under the strong assumption that they are virialised). The agreement between M_vir and M_tot is surprisingly good, considering that almost all the peaks are probably not yet virialised. According to the spherical collapse model, the peaks have already started or are about to start collapsing, and they are all predicted to be virialised by redshift z~0.8-1.6. We finally perform a careful comparison with the literature, given that smaller components of this proto-supercluster had previously been identified using heterogeneous galaxy samples. With VUDS, we obtain a panoramic view of this large structure, that encompasses, connects, and considerably expands in a homogeneous way on all previous detections of the various sub-components. This provides us the unique possibility to study a rich supercluster in formation.

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

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.

Using spectroscopic observations taken for the VIMOS Ultra-Deep Survey (VUDS) we report here on the discovery of PCl J1001+0220, a massive proto-cluster located at $z_{spec}\sim4.57$ in the COSMOS field. The proto-cluster was initially detected as a $\sim12\sigma$ overdensity of typical star-forming galaxies in the blind spectroscopic survey of the early universe ($2<z<6$) performed by VUDS. It was further mapped using a new technique developed that statistically combines spectroscopic and photometric redshifts, the latter derived from a recent compilation of deep multi-band imaging. Through various methods, the descendant halo mass of PCl J1001+0220 is estimated to be $\log(M_{h}/M_{\odot})_{z=0}\sim14.5-15$ with a large amount of mass apparently already in place at $z\sim4.57$. Tentative evidence is found for a fractional excess of older and more massive galaxies within the proto-cluster, an observation which suggests the pervasive early onset of vigorous star formation. No evidence is found for the differences in the star formation rates of member and a matched sample of coeval field galaxies either through rest-frame ultraviolet methods or through stacking extremely deep Very Large Array 3 GHz imaging. Additionally, no evidence for pervasive strong active galactic nuclei (AGN) activity is observed. Analysis of Hubble Space Telescope images provides weak evidence for for an elevated incidence of galaxy-galaxy interaction within the proto-cluster. The spectral properties of the two samples are compared, with a definite suppression of Ly$\alpha$ seen in the average member galaxy relative to the coeval field ($f_{esc,Ly\alpha}=1.8^{+0.3}_{-1.7}$% and $4.0^{+1.0}_{-0.8}$%, respectively). This observation along with other lines of evidence leads us to infer the possible presence of a large, cool diffuse medium within the proto-cluster environment evocative of a nascent intracluster medium.

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

To identify the galaxies responsible for the reionization of the Universe, we must rely on the investigation of the Lyman Continuum (LyC) properties of z<5 star-forming galaxies, where we can still directly observe their ionizing radiation. We selected a sample of 201 star-forming galaxies from the Vimos Ultra Deep Survey at 3.5<z<4.3 to explore the validity of some of the proposed indirect indicators of LyC radiation. We created subsamples of galaxies with EWLya>70{\AA}, Lya(ext)<5.7kpc, rUV<0.3kpc and |Dv Lya|<200km/s, stacked all the galaxies in each subsample and measured the flux density ratio fnu(895)/fnu(1470), that we consider to be a proxy for LyC emission. We compared these ratios to those obtained for the complementary samples. We find that the stacks of galaxies which are UV compact (rUV<0.3kpc) and have bright Lya emission (EWLya>70{\AA}), have much higher LyC fluxes compared to the rest of the galaxy population in agreement with theoretical studies and previous observational works. We also find that galaxies with low Lya(ext) have the highest LyC flux: this new correlation seems even stronger than the correlations with high EWLya and small rUV. These results assume that the stacks from all the subsamples present the same statistical contamination from lower redshift interlopers. If we subtract a statistical contamination from low redshift interlopers obtained with dedicated Monte Carlo simulations, from the flux density ratios (fnu(895)/fnu(1470)) of the significant subsamples we find that these samples contain real LyC leaking flux with a very high probability, but the true average escape fractions remain uncertain. Our work indicates that galaxies with very high EWLya, small Lya(ext) and small rUV are very likely the best candidates to show LyC radiation at z=4 and could therefore be the galaxies that contributed more to reionization.

We analyze the CIII]-1908\AA emission properties in a sample of 3899 star-forming galaxies (SFGs) at 2<z<3.8 from the VIMOS Ultra-Deep Survey (VUDS). We find a median EW(CIII])=2.0$\pm$0.2 to 2.2$\pm$0.2\AA ~for the SFG population at 2<z<3 and 3<z<4, resp. About 24% of SFGs are showing EW(CIII])>3\AA, including 20% with 3<EW(CIII])<10\AA ~and 4% with strong emission EW(CIII])>10\AA. A significant fraction of 1.2% of SFGs presents strong CIII] emission 20<EW(CIII])<40\AA. This makes CIII] the second most-frequent emission line in the UV-rest spectra of SFGs after Lyman-$\alpha$. We find a large dispersion in the weak EW(CIII]) - EW(Ly\alpha) correlation, with galaxies showing strong CIII] and no Ly$\alpha$, and vice-versa. SFGs with 10<EW(CIII])<0\AA ~present strong emission lines including CIV-1549, HeII-1640, OIII-1664, but also weaker emission of highly ionized elements SiIV-1403, NIV-1485, NIII-1750, or SiIII-1888, indicating a hard radiation field. We present a broad range of observational evidence supporting the presence of AGN in the strong CIII] emitting population. As EW(CIII]) is rising, we identify powerful outflows with velocities up to 1014 km/s. The strongest CIII] emitters are preferentially located below the SFGs main sequence, with the SFR reduced by x2. In addition, the median stellar age of the strongest emitters is ~0.8 Gyr, three times that of galaxies with EW(CIII])<10\AA. Spectral line analysis presented in a joint paper by Nakajima et al. (2017) confirms that the strongest emitters require an AGN. We conclude that these properties are indicative of AGN feedback in SFGs at 2<z<3.8, contributing to star-formation quenching. We find that quenching timescales of 0.25-0.5x10^9 years are necessary for such AGN feedback to turn part of the star-forming galaxy population with Mstar>10^10 MSun at z~3 into the population of quiescent galaxies observed at redshift z~1-2.

(Abridged) The properties of stellar clumps in star forming galaxies and their evolution over the redshift range $2\lesssim z \lesssim 6$ are presented and discussed in the context of the build-up of massive galaxies at early cosmic times. We use HST/ACS images of galaxies with spectroscopic redshifts from the VIMOS Ultra Deep Survey (VUDS) to identify clumps within a 20 kpc radius. We find that the population of galaxies with more than one clump is dominated by galaxies with two clumps, representing $\sim21-25$\% of the population, while the fraction of galaxies with 3, or 4 and more, clumps is 8-11 and 7-9\%, respectively. The fraction of clumpy galaxies is in the range $\sim35-55\%$ over $2<z<6$, increasing at higher redshifts, indicating that the fraction of irregular galaxies remains high up to the highest redshifts. The large and bright clumps (M$_{\star}\sim10^9$ up to $\sim10^{10}$M$_\odot$) are found to reside predominantly in galaxies with two clumps. Smaller and lower luminosity clumps ($\log_{10}\left(M_{\star}/\mathrm{M_\odot}\right)<9$) are found in galaxies with three clumps or more. We interpret these results as evidence for two different modes of clump formation working in parallel. The small low luminosity clumps are likely the result of disc fragmentation, with violent disc instabilities (VDI) forming several long-lived clumps in-situ, as suggested from simulations. A fraction of these clumps is also likely coming from minor mergers. The clumps in the dominating population of galaxies with two clumps are significantly more massive and have properties akin to those in merging pairs observed at similar redshifts; they appear as more massive than the most massive clumps observed in VDI numerical simulations.

We selected a sample of 76 Lya emitting galaxies from the VIMOS Ultra Deep Survey (VUDS) at 2<z<4. We estimated the velocity of the neutral gas flowing out of the interstellar medium as the velocity offset, Deltav, between the systemic redshift (zsys) and the center of low-ionization absorption line systems (LIS). To increase the SN of VUDS spectra, we stacked subsamples. We measured the systemic redshift from the rest-frame UV spectroscopic data using the CIII]1908 nebular emission line, and we considered SiII1526 as the highest signal-to-noise LIS line. We calculated the Lya peak shift with respect to the zsys, the EW(Lya), and the Lya spatial extension, Ext(Lya-C), from the profiles in the 2D stacked spectra. The galaxies that are faint in the rest-frame UV continuum, strong in Lya and CIII], with compact UV morphology, and localized in an underdense environment are characterized by outflow velocities of the order of a few hundreds of km/sec. The subsamples with smaller Deltav are characterized by larger Lya peak shifts, larger Ext(Lya-C), and smaller EW(Lya). In general we find that EW(Lya) anti-correlates with Ext(Lya-C) and Lya peak shift. We interpret these trends using a radiative-transfer shell model. The model predicts that an HI gas with a column density larger than 10^20/cm^2 is able to produce Lya peak shifts larger than >300km/sec. An ISM with this value of NHI would favour a large amount of scattering events, especially when the medium is static, so it can explain large values of Ext(Lya-C) and small EW(Lya). On the contrary, an ISM with a lower NHI, but large velocity outflows would lead to a Lya spatial profile peaked at the galaxy center (i.e. low values of Ext(Lya-C)) and to a large EW(Lya), as we see in our data. Our results and their interpretation via radiative-transfer models tell us that it is possible to use Lya to study the properties of the HI gas.

We present an extragalactic survey using observations from the Atacama Large Millimeter/submillimeter Array (ALMA) to characterise galaxy populations up to $z=0.35$: the Valpara\'iso ALMA Line Emission Survey (VALES). We use ALMA Band-3 CO(1--0) observations to study the molecular gas content in a sample of 67 dusty normal star-forming galaxies selected from the $Herschel$ Astrophysical Terahertz Large Area Survey ($H$-ATLAS). We have spectrally detected 49 galaxies at $>5\sigma$ significance and 12 others are seen at low significance in stacked spectra. CO luminosities are in the range of $(0.03-1.31)\times10^{10}$ K km s$^{-1}$ pc$^2$, equivalent to $\log({\rm M_{gas}/M_{\odot}}) =8.9-10.9$ assuming an $\alpha_{\rm CO}$=4.6(K km s$^{-1}$ pc$^{2}$)$^{-1}$, which perfectly complements the parameter space previously explored with local and high-z normal galaxies. We compute the optical to CO size ratio for 21 galaxies resolved by ALMA at $\sim 3$."$5$ resolution (6.5 kpc), finding that the molecular gas is on average $\sim$ 0.6 times more compact than the stellar component. We obtain a global Schmidt-Kennicutt relation, given by $\log [\Sigma_{\rm SFR}/({\rm M_{\odot} yr^{-1}kpc^{-2}})]=(1.26 \pm 0.02) \times \log [\Sigma_{\rm M_{H2}}/({\rm M_{\odot}\,pc^{-2}})]-(3.6 \pm 0.2)$. We find a significant fraction of galaxies lying at `intermediate efficiencies' between a long-standing mode of star-formation activity and a starburst, specially at $\rm L_{IR}=10^{11-12} L_{\odot}$. Combining our observations with data taken from the literature, we propose that star formation efficiencies can be parameterised by $\log [{\rm SFR/M_{H2}}]=0.19 \times {\rm (\log {L_{IR}}-11.45)}-8.26-0.41 \times \arctan[-4.84 (\log {\rm L_{IR}}-11.45) ]$. Within the redshift range we explore ($z<0.35$), we identify a rapid increase of the gas content as a function of redshift.

We present results on the dust attenuation of galaxies at redshift $\sim3-6$ by studying the relationship between the UV spectral slope ($\beta_{\rm UV}$) and the infrared excess (IRX; $L_{\rm IR}$/$L_{\rm UV}$) using ALMA far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of $M_{\ast}\sim 10^{10.7}\,M_{\rm \odot}$ spanning a redshift range $z=2.6-3.7$ (median $z=3.2$) that were observed with ALMA at $\lambda_{rest}=300\,{\rm \mu m}$. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX-$\beta_{\rm UV}$ relationship at $z\sim3$ is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to $\pm0.5$ dex. Nevertheless, the locally calibrated dust correction factors based on the IRX-$\beta_{\rm UV}$ relation are on average applicable to main-sequence $z\sim3$ galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of $z\sim4-6$ galaxies we find evidence for a significant evolution in the IRX-$\beta_{\rm UV}$ and the IRX-$M_{\ast}$ relations beyond $z\sim3$ toward lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger $z>3$ galaxy samples will be required to confirm this intriguing redshift evolution.

We report the detection of morphology dependent stellar age in massive quenched galaxies (QGs) at z~1.2. The sense of the dependence is that compact QGs are 0.5-2 Gyr older than normal-sized ones. The evidence comes from three different age indicators, Dn4000, Hdelta and fits to spectral synthesis models, applied to their stacked optical spectra. All age indicators consistently show that the stellar populations of compact QGs are older than their normally-sized counterparts. We detect weak [OII] emission in a fraction of QGs, and the strength of the line, when present, is similar between the two samples; however, compact galaxies exhibit significantly lower frequency of [OII] emission than normal ones. A fraction of both samples are individually detected in 7 Ms Chandra X-ray images (luminosities$\sim10^{40}$-$10^{41}$ erg/sec). 7 Ms stacks of non-detected galaxies show similarly low luminosities in the soft band only, consistent with a hot gas origin for the X-ray emission. While both [OII] emitters and non-emitters are also X-ray sources among normal galaxies, no compact galaxy with [OII] emission is an X-ray source, arguing against an AGN powering the line in compact galaxies. We interpret the [OII] properties as further evidence that compact galaxies are older and further along into the process of quenching star-formation and suppressing gas accretion. Finally, we argue that the older age of compact QGs is evidence of progenitor bias: compact QGs simply reflect the smaller sizes of galaxies at their earlier quenching epoch, with stellar density most likely having nothing directly to do with cessation of star-formation.

Determining the average fraction of Lyman continuum (LyC) photons escaping high redshift galaxies is essential for understanding how reionization proceeded in the z>6 Universe. We want to measure the LyC signal from a sample of sources in the Chandra Deep Field South (CDFS) and COSMOS fields for which ultra-deep VIMOS spectroscopy as well as multi-wavelength Hubble Space Telescope (HST) imaging are available. We select a sample of 46 galaxies at $z\sim 4$ from the VIMOS Ultra Deep Survey (VUDS) database, such that the VUDS spectra contain the LyC part of the spectra, that is, the rest-frame range $880-910\AA$. Taking advantage of the HST imaging, we apply a careful cleaning procedure and reject all the sources showing nearby clumps with different colours, that could potentially be lower-redshift interlopers. After this procedure, the sample is reduced to 33 galaxies. We measure the ratio between ionizing flux (LyC at $895\AA$) and non-ionizing emission (at $\sim 1500 \AA$) for all individual sources. We also produce a normalized stacked spectrum of all sources. Assuming an intrinsic average $L_{\nu}(1470)/L_{\nu}(895)$ of 3, we estimate the individual and average relative escape fraction. We do not detect ionizing radiation from any individual source, although we identify a possible LyC emitter with very high Ly$\alpha$ equivalent width (EW). From the stacked spectrum and assuming a mean transmissivity for the sample, we measure a relative escape fraction $f_{esc}^{rel}=0.09\pm0.04$. We also look for correlations between the limits in the LyC flux and source properties and find a tentative correlation between LyC flux and the EW of the Ly$\alpha$ emission line. Our results imply that the LyC flux emitted by $V=25-26$ star-forming galaxies at z$\sim$4 is at most very modest, in agreement with previous upper limits from studies based on broad and narrow band imaging.

Deep observations are revealing a growing number of young galaxies in the first billion year of cosmic time. Compared to typical galaxies at later times, they show more extreme emission-line properties, higher star formation rates, lower masses, and smaller sizes. However, their faintness precludes studies of their chemical abundances and ionization conditions, strongly limiting our understanding of the physics driving early galaxy build-up and metal enrichment. Here we study a rare population of UV-selected, sub$-L^{*}$(z=3) galaxies at redshift 2.4$<z<$3.5 that exhibit all the rest-frame properties expected from primeval galaxies. These low-mass, highly-compact systems are rapidly-forming galaxies able to double their stellar mass in only few tens million years. They are characterized by very blue UV spectra with weak absorption features and bright nebular emission lines, which imply hard radiation fields from young hot massive stars. Their highly-ionized gas phase has strongly sub-solar carbon and oxygen abundances, with metallicities more than a factor of two lower than that found in typical galaxies of similar mass and star formation rate at $z\lesssim$2.5. These young galaxies reveal an early and short stage in the assembly of their galactic structures and their chemical evolution, a vigorous phase which is likely to be dominated by the effects of gas-rich mergers, accretion of metal-poor gas and strong outflows.

We investigate of the properties of $\sim$2000 Herschel/SPIRE-selected galaxies from $0<z<4$ using a combination of extensive spectroscopy, deep imaging from CFHT, VLA, Spitzer, XMM-Newton, and Herschel, and well-calibrated SED fitting. Herschel galaxies are observed to span a range of stellar masses, colors, and absolute magnitudes equivalent to galaxies undetected in SPIRE. Though many Herschel galaxies appear to be in transition, such galaxies are largely consistent with normal star-forming galaxies when rest-frame colors are utilized. The nature of the star-forming "main sequence" is studied and we warn against adopting this framework unless the main sequence is determined precisely. Herschel galaxies at different total infrared luminosities ($L_{TIR}$) are compared. Bluer colors, larger nebular extinctions, and larger contributions from younger stellar populations are observed for galaxies with larger $L_{TIR}$, suggesting that low-$L_{TIR}$ galaxies are undergoing rejuvenated starbursts while galaxies with higher $L_{TIR}$ are forming a larger percentage of their stellar mass. A variety of methods are used to select powerful active galactic nuclei (AGN). Galaxies hosting AGN are observed to be undergoing starbursts more commonly and vigorously than a matched sample of galaxies without powerful AGN and, additionally, the fraction of galaxies with an AGN increases with increasing star formation rate at all redshifts. At all redshifts ($0<z<4$) the most prodigious star-forming galaxies are found to contain the highest fraction of powerful AGN. For redshift bins that allow a comparison ($z>0.5$), the highest $L_{TIR}$ galaxies in a given redshift bin are unobserved by SPIRE at subsequently lower redshifts, a trend linked to downsizing. In conjunction with other results, this evidence is used to argue for prevalent AGN-driven quenching in starburst galaxies across cosmic time.

The observed evolution of the gas fraction and its associated depletion time in main sequence (MS) galaxies provides insights on how star formation proceeds over cosmic time. We report ALMA detections of the rest-frame $\sim$300$\mu$m continuum observed at 240 GHz for 45 massive ($\rm \langle log(M_{\star}(M_{\odot}))\rangle=10.7$), normal star forming ($\rm \langle log(sSFR(yr^{-1}))\rangle=-8.6$), i.e. MS, galaxies at $\rm z\approx3.2$ in the COSMOS field. From an empirical calibration between cold neutral, i.e. molecular and atomic, gas mass $\rm M_{gas}$ and monochromatic (rest-frame) infrared luminosity, the gas mass for this sample is derived. Combined with stellar mass $\rm M_{\star}$ and star formation rate (SFR) estimates (from {\sc MagPhys} fits) we obtain a median gas fraction of $\rm \mu_{gas}=M_{gas}/M_{\star}=1.65_{-0.19}^{+0.18}$ and a median gas depletion time $\rm t_{depl.}(Gyr)=M_{gas}/SFR=0.68_{-0.08}^{+0.07}$; correction for the location on the MS will only slightly change the values. The reported uncertainties are the $\rm 1\sigma$ error on the median. Our results are fully consistent with the expected flattening of the redshift evolution from the 2-SFM (2 star formation mode) framework that empirically prescribes the evolution assuming a universal, log-linear relation between SFR and gas mass coupled to the redshift evolution of the specific star formation rate (sSFR) of main sequence galaxies. While $\rm t_{dep.}$ shows only a mild dependence on location within the MS, a clear trend of increasing $\rm \mu_{gas}$ across the MS is observed (as known from previous studies). Further we comment on trends within the MS and (in)consistencies with other studies.

We present a study of a large filamentary structure at z~0.73 in the field of the COSMOS survey, the so-called COSMOS Wall. This structure encompasses a comprehensive range of environments from a dense cluster and a number of galaxy groups to filaments, less dense regions, and adjacent voids. It thus provides a valuable laboratory for the accurate mapping of environmental effects on galaxy evolution at a look-back time of ~6.5 Gyr, when the Universe was roughly half its present age. We performed deep spectroscopic observations with VIMOS at VLT of a K-band selected sample of galaxies in this complex structure, building a sample of galaxies complete in galaxy stellar mass down to a lower limit of log(M/M_sun)~ 9.8, which is significantly deeper than previously available data. Thanks to its location within the COSMOS survey, each galaxy benefits from a wealth of ancillary information. In this paper we detail the survey strategy and weighting scheme adopted to account for the biases introduced by the photometric pre-selection of our targets. We present our galaxy stellar mass and rest-frame magnitudes estimates together with a group catalog obtained with our new data and their member galaxies color/mass distribution. Owing to to our new sample we can perform a detailed, high definition mapping of the complex COSMOS Wall structure. The sharp environmental information, coupled with high quality spectroscopic information and rich ancillary data available in the COSMOS field, enables a detailed study of galaxy properties as a function of local environment in a redshift slice where environmental effects are important, and in a stellar mass range where mass and environment driven effects are both at work.

In this paper we investigate the impact of different star formation histories (SFHs) on the relation between stellar mass M$_{*}$ and star formation rate (SFR) using a sample of galaxies with reliable spectroscopic redshift zspec>2 drawn from the VIMOS Ultra-Deep Survey (VUDS). We produce an extensive database of dusty model galaxies, calculated starting from the new library of single stellar population (SSPs) models presented in Cassara' et al. 2013 and weighted by a set of 28 different SFHs based on the Schmidt function, and characterized by different ratios of the gas infall time scale $\tau_{infall}$ to the star formation efficiency $\nu$. The treatment of dust extinction and re-emission has been carried out by means of the radiative transfer calculation. The spectral energy distribution (SED) fitting technique is performed by using GOSSIP+, a tool able to combine both photometric and spectroscopic information to extract the best value of the physical quantities of interest, and to consider the Intergalactic Medium (IGM) attenuation as a free parameter. We find that the main contribution to the scatter observed in the $SFR-M_{*}$ plane is the possibility of choosing between different families of SFHs in the SED fitting procedure, while the redshift range plays a minor role. The majority of the galaxies, at all cosmic times, are best-fit by models with SFHs characterized by a high $\tau_{\rm infall}/\nu$ ratio. We discuss the reliability of the presence of a small percentage of dusty and highly star forming galaxies, in the light of their detection in the FIR.

We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3<z<6$, using the unique database of spectroscopic redshifts in COSMOS, we find $\sigma_{\Delta z/(1+z_s)}$ = 0.021 and $\eta=13.2\% $. The deepest regions reach a 90\% completeness limit of 10$^{10}M_\odot$ to $z=4$. Detailed comparisons of the color distributions, number counts, and clustering show excellent agreement with the literature in the same mass ranges. COSMOS2015 represents a unique, publicly available, valuable resource with which to investigate the evolution of galaxies within their environment back to the earliest stages of the history of the Universe. The COSMOS2015 catalog is distributed via anonymous ftp (ftp://ftp.iap.fr/pub/from_users/hjmcc/COSMOS2015/) and through the usual astronomical archive systems (CDS, ESO Phase 3, IRSA).

This paper describes the first data release (DR1) of the VIMOS Ultra Deep Survey (VUDS). The DR1 includes all low-resolution spectroscopic data obtained in 276.9 arcmin2 of the CANDELS-COSMOS and CANDELS-ECFDS survey areas, including accurate spectroscopic redshifts z_spec and individual spectra obtained with VIMOS on the ESO-VLT. A total of 698 objects have a measured redshift, with 677 galaxies, two type-I AGN and a small number of 19 contaminating stars. The targets of the spectroscopic survey are selected primarily on the basis of their photometric redshifts to ensure a broad population coverage. About 500 galaxies have z_spec>2, 48 with z_spec>4, and the highest reliable redshifts reach beyond z_spec=6. This dataset approximately doubles the number of galaxies with spectroscopic redshifts at z>3 in these fields. We discuss the general properties of the sample in terms of the spectroscopic redshift distribution, the distribution of Lyman-alpha equivalent widths, and physical properties including stellar masses M_star and star formation rates (SFR) derived from spectral energy distribution fitting with the knowledge of z_spec. We highlight the properties of the most massive star-forming galaxies, noting the large range in spectral properties, with Lyman-alpha in emission or in absorption, and in imaging properties with compact, multi-component or pair morphologies. We present the catalogue database and data products. All data are publicly available and can be retrieved from a dedicated query-based database available at http://cesam.lam.fr/vuds.

We aim at improving constraints on the epoch of galaxy formation by measuring the ages of 3597 galaxies with spectroscopic redshifts 2<z<6.5 in the VIMOS Ultra Deep Survey (VUDS). We derive ages and other physical parameters from the simultaneous fitting with the GOSSIP+ software of observed UV rest-frame spectra and photometric data from the u-band up to 4.5 microns using composite stellar population models. We conclude from extensive simulations that at z>2 the joint analysis of spectroscopy and photometry combined with restricted age possibilities when taking into account the age of the Universe substantially reduces systematic uncertainties and degeneracies in the age derivation. We find galaxy ages ranging from very young with a few tens of million years to substantially evolved with ages up to ~1.5-2 Gyr. The formation redshifts z_f derived from the measured ages indicate that galaxies may have started forming stars as early as z_f~15. We produce the formation redshift function (FzF), the number of galaxies per unit volume formed at a redshift z_f, and compare the FzF in increasing redshift bins finding a remarkably constant 'universal' FzF. The FzF is parametrized with (1+z)^\zeta, with \zeta~0.58+/-0.06, indicating a smooth 2 dex increase from z~15 to z~2. Remarkably this observed increase is of the same order as the observed rise in the star formation rate density (SFRD). The ratio of the SFRD with the FzF gives an average SFR per galaxy of ~7-17Msun/yr at z~4-6, in agreement with the measured SFR for galaxies at these redshifts. From the smooth rise in the FzF we infer that the period of galaxy formation extends from the highest possible redshifts that we can probe at z~15 down to redshifts z~2. This indicates that galaxy formation is a continuous process over cosmic time, with a higher number of galaxies forming at the peak in SFRD at z~2 than at earlier epochs. (Abridged)

We measure galaxy sizes on a sample of $\sim1200$ galaxies with confirmed spectroscopic redshifts $2 \leq z_{spec} \leq 4.5$ in the VIMOS Ultra Deep Survey (VUDS), representative of star-forming galaxies with $i_\mathrm{AB} \leq 25$. We first derive galaxy sizes applying a classical parametric profile fitting method using GALFIT. We then measure the total pixel area covered by a galaxy above a given surface brightness threshold, which overcomes the difficulty of measuring sizes of galaxies with irregular shapes. We then compare the results obtained for the equivalent circularized radius enclosing 100\% of the measured galaxy light $r_T^{100}$ to those obtained with the effective radius $r_{e,\mathrm{circ}}$ measured with GALFIT. We find that the sizes of galaxies computed with our non-parametric approach span a large range but remain roughly constant on average with a median value $r_T^{100}\sim2.2$ kpc for galaxies with $2<z<4.5$. This is in stark contrast with the strong downward evolution of $r_e$ with increasing redshift, down to sizes of $<1$ kpc at $z\sim4.5$. We analyze the difference and find that parametric fitting of complex, asymmetric, multi-component galaxies is severely underestimating their sizes. By comparing $r_T^{100}$ with physical parameters obtained through SED fitting we find that the star-forming galaxies that are the largest at any redshift are, on average, more massive and more star-forming. We discover that galaxies present more concentrated light profiles as we move towards higher redshifts. We interpret these results as the signature of several, possibly different, evolutionary paths of galaxies in their early stages of assembly, including major and minor merging or star-formation in multiple bright regions. (abridged)

We measure the H{\alpha} and [OIII] emission line properties as well as specific star-formation rates (sSFR) of spectroscopically confirmed 3<z<6 galaxies in COSMOS from their observed colors vs. redshift evolution. Our model describes consistently the ensemble of galaxies including intrinsic properties (age, metallicity, star-formation history), dust-attenuation, and optical emission lines. We forward-model the measured H{\alpha} equivalent-widths (EW) to obtain the sSFR out to z~6 without stellar mass fitting. We find a strongly increasing rest-frame H{\alpha} EW that is flattening off above z~2.5 with average EWs of 300-600A at z~6. The sSFR is increasing proportional to (1+z)^2.4 at z<2.2 and (1+z)^1.5 at higher redshifts, indicative of a fast mass build-up in high-z galaxies within e-folding times of 100-200Myr at z~6. The redshift evolution at z>3 cannot be fully explained in a picture of cold accretion driven growth. We find a progressively increasing [OIII]{\lambda}5007/H{\beta} ratio out to z~6, consistent with the ratios in local galaxies selected by increasing H{\alpha} EW (i.e., sSFR). This demonstrates the potential of using "local high-z analogs" to investigate the spectroscopic properties and relations of galaxies in the re-ionization epoch.

Aim: We aim to measure the LyC signal from a sample of sources in the Chandra deep field south. We collect star-forming galaxies (SFGs) and active galactic nuclei (AGN) with accurate spectroscopic redshifts, for which Hubble Space Telescope (HST) coverage and multi-wavelength photometry are available. Method: We selected a sample of about 200 sources at z~3. Taking advantage of HST resolution, we applied a careful cleaning procedure and rejected sources showing nearby clumps with different colours, which could be lower-z interlopers. Our clean sample consisted of 86 SFGs (including 19 narrow-band selected Lya emitters) and 8 AGN (including 6 detected in X-rays). We measured the LyC flux from aperture photometry in four narrow-band filters covering wavelengths below a 912 A rest frame (3.11<z<3.53). We estimated the ratio between ionizing (LyC flux) and 1400 A non-ionizing emissions for AGN and galaxies. Results: By running population synthesis models, we assume an average intrinsic L(1400 A)/L(900 A) ratio of 5 as the representative value for our sample. With this value and an average treatment of the lines of sight of the inter-galactic medium, we estimate the LyC escape fraction relative to the intrinsic value (fesc_rel(LyC)). We do not directly detect ionizing radiation from any individual SFG, but we are able to set a 1(2)sigma upper limit of fesc_rel(LyC)<12(24)%. This result is consistent with other non-detections published in the literature. No meaningful limits can be calculated for the sub-sample of Lya emitters. We obtain one significant direct detection for an AGN at z=3.46, with fesc_rel(LyC) = (72+/-18)%. Conclusions: Our upper limit on fescrel(LyC) implies that the SFGs studied here do not present either the physical properties or the geometric conditions suitable for efficient LyC-photon escape.

The aim of this paper is to investigate spectral and photometric properties of 854 faint ($i_{AB}$<~25 mag) star-forming galaxies (SFGs) at 2<z<2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy based on their photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes ($\beta$) as well as Ly$\alpha$ equivalent widths (EW). On average, the spectroscopically measured $\beta$ (-1.36$\pm$0.02), is comparable to the photometrically measured $\beta$ (-1.32$\pm$0.02), and has smaller measurement uncertainties. The positive correlation of $\beta$ with the Spectral Energy Distribution (SED)-based measurement of dust extinction, E$_{\rm s}$(B-V), emphasizes the importance of $\beta$ as an alternative dust indicator at high redshifts. To make a proper comparison, we divide these SFGs into three subgroups based on their rest-frame Ly$\alpha$ EW: SFGs with no Ly$\alpha$ emission (SFG$_{\rm N}$; EW$\le$0\AA), SFGs with Ly$\alpha$ emission (SFG$_{\rm L}$; EW$>$0\AA), and Ly$\alpha$ emitters (LAEs; EW$\ge$20\AA). The fraction of LAEs at these redshifts is $\sim$10%, which is consistent with previous observations. We compared best-fit SED-estimated stellar parameters of the SFG$_{\rm N}$, SFG$_{\rm L}$ and LAE samples. For the luminosities probed here ($\sim$L$^*$), we find that galaxies with and without Ly$\alpha$ in emission have small but significant differences in their SED-based properties. We find that LAEs have less dust, and lower star-formation rates (SFR) compared to non-LAEs. We also find that LAEs are less massive compared to non-LAEs, though the difference is smaller and less significant compared to the SFR and E$_{\rm s}$(B-V). [abridged]

We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wide range of CANDELS science goals. This scheme includes the basic Hubble sequence types, but also includes a detailed look at mergers and interactions, the clumpiness of galaxies, $k$-corrections, and a variety of other structural properties. In this paper, we focus on the first field to be completed -- GOODS-S. The wide area coverage spanning the full field includes 7634 galaxies that have been classified by at least three different people. In the deep area of the field, 2534 galaxies have been classified by at least five different people at three different depths. With this paper, we release to the public all of the visual classifications in GOODS-S along with the GUI that we developed to classify galaxies. We find that the level of agreement among classifiers is good and depends on both the galaxy magnitude and the galaxy type, with disks showing the highest level of agreement and irregulars the lowest. A comparison of our classifications with the Sersic index and rest-frame colors shows a clear separation between disk and spheroid populations. Finally, we explore morphological k-corrections between the V-band and H-band observations and find that a small fraction (84 galaxies in total) are classified as being very different between these two bands. These galaxies typically have very clumpy and extended morphology or are very faint in the V-band.

The Lyman continuum (LyC) flux escaping from high-z galaxies into the IGM is a fundamental quantity to understand the physical processes involved in the reionization epoch. We have investigated a sample of star-forming galaxies at z~3.3 in order to search for possible detections of LyC photons escaping from galaxy halos. UV deep imaging in the COSMOS field obtained with the prime focus camera LBC at the LBT telescope was used together with a catalog of spectroscopic redshifts obtained by the VIMOS Ultra Deep Survey (VUDS) to build a sample of 45 galaxies at z~3.3 with L>0.5L*. We obtained deep LBC images of galaxies with spectroscopic redshifts in the interval 3.27<z<3.40 both in the R and deep U bands. A sub-sample of 10 galaxies apparently shows escape fractions>28% but a detailed analysis of their properties reveals that, with the exception of two marginal detections (S/N~2) in the U band, all the other 8 galaxies are most likely contaminated by the UV flux of low-z interlopers located close to the high-z targets. The average escape fraction derived from the stacking of the cleaned sample was constrained to fesc_rel<2%. The implied HI photo-ionization rate is a factor two lower than that needed to keep the IGM ionized at z~3, as observed in the Lyman forest of high-z QSO spectra or by the proximity effect. These results support a scenario where high redshift, relatively bright (L>0.5L*) star-forming galaxies alone are unable to sustain the level of ionization observed in the cosmic IGM at z~3. Star-forming galaxies at higher redshift and at fainter luminosities (L<<L*) can be the major contributors to the reionization of the Universe only if their physical properties are subject to rapid changes from z~3 to z~6-10. Alternatively, ionizing sources could be discovered looking for fainter sources among the AGN population at high-z.

We report the discovery of 31 low-luminosity (-14.5 > M_{AB}(B) > -18.8), extreme emission line galaxies (EELGs) at 0.2 < z < 0.9 identified by their unusually high rest-frame equivalent widths (100 < EW[OIII] < 1700 A) as part of the VIMOS Ultra Deep Survey (VUDS). VIMOS optical spectra of unprecedented sensitivity ($I_{AB}$ ~ 25 mag) along with multiwavelength photometry and HST imaging are used to investigate spectrophotometric properties of this unique sample and explore, for the first time, the very low stellar mass end (M* < 10^8 M$_{\odot}$) of the luminosity-metallicity (LZR) and mass-metallicity (MZR) relations at z < 1. Characterized by their extreme compactness (R50 < 1 kpc), low stellar mass and enhanced specific star formation rates (SFR/M* ~ 10^{-9} - 10^{-7} yr^{-1}), the VUDS EELGs are blue dwarf galaxies likely experiencing the first stages of a vigorous galaxy-wide starburst. Using T_e-sensitive direct and strong-line methods, we find that VUDS EELGs are low-metallicity (7.5 < 12+log(O/H) < 8.3) galaxies with high ionization conditions, including at least three EELGs showing HeII 4686A emission and four EELGs of extremely metal-poor (<10% solar) galaxies. The LZR and MZR followed by EELGs show relatively large scatter, being broadly consistent with the extrapolation toward low luminosity and mass from previous studies at similar redshift. However, we find evidences that galaxies with younger and more vigorous star formation -- as characterized by their larger EWs, ionization and sSFR -- tend to be more metal-poor at a given stellar mass.

Context. Lyman Break Galaxies (LBG) and Narrow Band (NB) surveys have been successful at detecting large samples of high-redshift galaxies. Both methods are subject to contamination from low-redshift interlopers. Aims. In this paper, our aim is to investigate the nature of low-redshift interlopers in NB Lyman-$\alpha$ emitters (LAE) searches. Methods. From previous HAWK-I NB imaging at z $\sim$7.7 we identify three objects that would have been selected as high-redshift LAEs had our optical data been one magnitude shallower (but still one to two magnitudes fainter than the near infrared data). We follow-up these objects in spectroscopy with XSHOOTER at the VLT. Results. Despite low quality data due to bad weather conditions, for each of the three objects we identify one, and only one emission line, in the spectra of the objects, that we identify as the O[III]5007A line. This result combined to spectral energy density fitting and tests based on line ratios of several populations of galaxies we infer that the 3 objects are ultrastrong line emitters at redshifts $\sim$1.1. Conclusions. From this work and the literature we remark that the O[III] line appears to be a common source of contamination in high-redshift LBG and LAE samples and we suggest that efforts be put to characterize with high accuracy the O[III] luminosity function out to redshift $\sim$3 or higher.

The relation between the galaxy stellar mass M_star and the dark matter halo mass M_h gives important information on the efficiency in forming stars and assembling stellar mass in galaxies. We present the stellar mass to halo mass ratio (SMHR) measurements at redshifts 2<z<5, obtained from the VIMOS Ultra Deep Survey. We use halo occupation distribution (HOD) modelling of clustering measurements on ~3000 galaxies with spectroscopic redshifts to derive the dark matter halo mass M_h, and SED fitting over a large set of multi-wavelength data to derive the stellar mass M_star and compute the SMHR=M_star/M_h. We find that the SMHR ranges from 1% to 2.5% for galaxies with M_star=1.3x10^9 M_sun to M_star=7.4x10^9 M_sun in DM halos with M_h=1.3x10^{11} M_sun} to M_h=3x10^{11} M_sun. We derive the integrated star formation efficiency (ISFE) of these galaxies and find that the star formation efficiency is a moderate 6-9% for lower mass galaxies while it is relatively high at 16% for galaxies with the median stellar mass of the sample ~7x10^9 M_sun. The lower ISFE at lower masses may indicate that some efficient means of suppressing star formation is at work (like SNe feedback), while the high ISFE for the average galaxy at z~3 is indicating that these galaxies are efficiently building-up their stellar mass at a key epoch in the mass assembly process. We further infer that the average mass galaxy at z~3 will start experiencing star formation quenching within a few hundred millions years.

Quenched galaxies at z>2 are nearly all very compact relative to z~0, suggesting a physical connection between high stellar density and efficient, rapid cessation of star-formation. We present restframe UV spectra of Lyman-break galaxies (LBGs) at z~3 selected to be candidate progenitors of quenched galaxies at z~2 based on their compact restframe optical sizes and high surface density of star-formation. We compare their UV properties to those of more extended LBGs of similar mass and star formation rate (non-candidates). We find that candidate progenitors have faster ISM gas velocities and higher equivalent widths of interstellar absorption lines, implying larger velocity spread among absorbing clouds. Candidates deviate from the relationship between equivalent widths of Lyman-alpha and interstellar absorption lines in that their Lyman-alpha emission remains strong despite high interstellar absorption, possibly indicating that the neutral HI fraction is patchy such that Lyman-alpha photons can escape. We detect stronger CIV P-Cygni features (emission and absorption) and HeII emission in candidates, indicative of larger populations of metal rich Wolf-Rayet stars compared to non-candidates. The faster bulk motions, broader spread of gas velocity, and Lyman-alpha properties of candidates are consistent with their ISM being subject to more energetic feedback than non-candidates. Together with their larger metallicity (implying more evolved star-formation activity) this leads us to propose, if speculatively, that they are likely to quench sooner than non-candidates, supporting the validity of selection criteria used to identify them as progenitors of z~2 passive galaxies. We propose that massive, compact galaxies undergo more rapid growth of stellar mass content, perhaps because the gas accretion mechanisms are different, and quench sooner than normally-sized LBGs at these early epochs.

Using new spectroscopic observations obtained as part of the VIMOS Ultra-Deep Survey (VUDS), we perform a systematic search for overdense environments in the early universe ($z>2$) and report here on the discovery of Cl J0227-0421, a massive protocluster at $z=3.29$. This protocluster is characterized by both the large overdensity of spectroscopically confirmed members, $\delta_{gal}=10.5\pm2.8$, and a significant overdensity in photometric redshift members. The halo mass of this protocluster is estimated, by a variety of methods, to be roughly $3\times10^{14}$ $\mathcal{M}_{\odot}$ at $z\sim3.3$, which, evolved to $z=0$ results in a halo mass rivaling or exceeding that of the Coma cluster. The properties of 19 spectroscopically confirmed member galaxies are compared with a large sample of VUDS/VVDS galaxies in lower density field environments at similar redshifts. We find tentative evidence for an excess of redder, brighter, and more massive galaxies within the confines of the protocluster relative to the field population, which suggests that we may be observing the beginning of environmentally-induced quenching. The properties of these galaxies are investigated, including a discussion of the brightest protocluster galaxy which appears to be undergoing vigorous coeval nuclear and starburst activity. The remaining member galaxies appear to have characteristics which are largely similar to the field population. Though we find weaker evidence of the suppression of the median star formation rates amongst and differences in stacked spectra of member galaxies with respect to the field, we defer any conclusions of these trends to future work with the ensemble of protostructures that are found in the full VUDS sample.

We investigate the evolution of galaxy clustering for galaxies in the redshift range 2.0<$z$<5.0 using the VIMOS Ultra Deep Survey (VUDS). We present the projected (real-space) two-point correlation function $w_p(r_p)$ measured by using 3022 galaxies with robust spectroscopic redshifts in two independent fields (COSMOS and VVDS-02h) covering in total 0.8 deg$^2$. We quantify how the scale dependent clustering amplitude $r_0$ changes with redshift making use of mock samples to evaluate and correct the survey selection function. Using a power-law model $\xi(r) = (r/r_0)^{-\gamma}$ we find that the correlation function for the general population is best fit by a model with a clustering length $r_0$=3.95$^{+0.48}_{-0.54}$ h$^{-1}$Mpc and slope $\gamma$=1.8$^{+0.02}_{-0.06}$ at $z$~2.5, $r_0$=4.35$\pm$0.60 h$^{-1}$Mpc and $\gamma$=1.6$^{+0.12}_{-0.13}$ at $z$~3.5. We use these clustering parameters to derive the large-scale linear galaxy bias $b_L^{PL}$, between galaxies and dark matter. We find $b_L^{PL}$ = 2.68$\pm$0.22 at redshift $z$~3 (assuming $\sigma_8$ = 0.8), significantly higher than found at intermediate and low redshifts. We fit an HOD model to the data and we obtain that the average halo mass at redshift $z$~3 is $M_h$=10$^{11.75\pm0.23}$ h$^{-1}$M$_{\odot}$. From this fit we confirm that the large-scale linear galaxy bias is relatively high at $b_L^{HOD}$ = 2.82$\pm$0.27. Comparing these measurements with similar measurements at lower redshifts we infer that the star-forming population of galaxies at $z$~3 should evolve into the massive and bright ($M_r$<-21.5) galaxy population which typically occupy haloes of mass $\langle M_h\rangle$ = 10$^{13.9}$ h$^{-1}$ $M_{\odot}$ at redshift $z$=0.

We study the evolution of the star formation rate (SFR) - stellar mass (M_star) relation and specific star formation rate (sSFR) of star forming galaxies (SFGs) since a redshift z~5.5 using 2435 (4531) galaxies with highly reliable (reliable) spectroscopic redshifts in the VIMOS Ultra-Deep Survey (VUDS). It is the first time that these relations can be followed over such a large redshift range from a single homogeneously selected sample of galaxies with spectroscopic redshifts. The log(SFR) - log(M_star) relation for SFGs remains roughly linear all the way up to z=5 but the SFR steadily increases at fixed mass with increasing redshift. We find that for stellar masses M_star>3.2 x 10^9 M_sun the SFR increases by a factor ~13 between z=0.4 and z=2.3. We extend this relation up to z=5, finding an additional increase in SFR by a factor 1.7 from z=2.3 to z=4.8 for masses M_star > 10^10 M_sun. We observe a turn-off in the SFR-M_star relation at the highest mass end up to a redshift z~3.5. We interpret this turn-off as the signature of a strong on-going quenching mechanism and rapid mass growth. The sSFR increases strongly up to z~2 but it grows much less rapidly in 2<z<5. We find that the shape of the sSFR evolution is not well reproduced by cold gas accretion-driven models or the latest hydrodynamical models. Below z~2 these models have a flatter evolution (1+z)^{Phi} with Phi=2-2.25 compared to the data which evolves more rapidly with Phi=2.8+-0.2. Above z~2, the reverse is happening with the data evolving more slowly with Phi=1.2+-0.1. The observed sSFR evolution over a large redshift range 0<z<5 and our finding of a non linear main sequence at high mass both indicate that the evolution of SFR and M_star is not solely driven by gas accretion. The results presented in this paper emphasize the need to invoke a more complex mix of physical processes {abridge}

Aims. The aim of this work is to constrain the evolution of the fraction of Lya emitters among UV selected star forming galaxies at 2<z<6, and to measure the stellar escape fraction of Lya photons over the same redshift range. Methods. We exploit the ultradeep spectroscopic observations collected by the VIMOS Ultra Deep Survey (VUDS) to build an unique, complete and unbiased sample of 4000 spectroscopically confirmed star forming galaxies at 2<z<6. Our galaxy sample UV luminosities brighter than M* at 2<z<6, and luminosities down to one magnitude fainter than M* at 2<z<3.5. Results. We find that 80% of the star forming galaxies in our sample have EW0(Lya)<10A, and correspondingly fesc(Lya)<1%. By comparing these results with literature, we conclude that the bulk of the Lya luminosity at 2<z<6 comes from galaxies that are fainter in the UV than those we sample in this work. The strong Lya emitters constitute, at each redshift, the tail of the distribution of the galaxies with extreme EW0(Lya) and fesc(Lya) . This tail of large EW0 and fesc(Lya) becomes more important as the redshift increases, and causes the fraction of Lya with EW0> 25A to increase from 5% at z=2 to 30% at z=6, with the increase being relatively stronger beyond z=4. We observe no difference, for the narrow range of UV luminosities explored in this work, between the fraction of strong Lya emitters among galaxies fainter or brighter than M*, although the fraction for the FUV faint galaxies evolves faster, at 2<z<3.5, than for the bright ones. We do observe an anticorrelation between E(B-V) and fesc(Lya): generally galaxies with high fesc(Lya) have also small amounts of dust (and viceversa). However, when the dust content is low (E(B-V)<0.05) we observe a very broad range of fesc(Lya), ranging from 10^-3 to 1. This implies that the dust alone is not the only regulator of the amount of escaping Lya photons.

We present the first observations of foreground Lyman-$\alpha$ forest absorption from high-redshift galaxies, targeting 24 star-forming galaxies (SFGs) with $z\sim 2.3-2.8$ within a $5' \times 15'$ region of the COSMOS field. The transverse sightline separation is $\sim 2\,h^{-1}\mathrm{Mpc}$ comoving, allowing us to create a tomographic reconstruction of the 3D Ly$\alpha$ forest absorption field over the redshift range $2.20\leq z\leq 2.45$. The resulting map covers $6\,h^{-1}\mathrm{Mpc} \times 14\,h^{-1}\mathrm{Mpc}$ in the transverse plane and $230\,h^{-1}\mathrm{Mpc}$ along the line-of-sight with a spatial resolution of $\approx 3.5\,h^{-1}\mathrm{Mpc}$, and is the first high-fidelity map of large-scale structure on $\sim\mathrm{Mpc}$ scales at $z>2$. Our map reveals significant structures with $\gtrsim 10\,h^{-1}\mathrm{Mpc}$ extent, including several spanning the entire transverse breadth, providing qualitative evidence for the filamentary structures predicted to exist in the high-redshift cosmic web. Simulated reconstructions with the same sightline sampling, spectral resolution, and signal-to-noise ratio recover the salient structures present in the underlying 3D absorption fields. Using data from other surveys, we identified 18 galaxies with known redshifts coeval with our map volume enabling a direct comparison to our tomographic map. This shows that galaxies preferentially occupy high-density regions, in qualitative agreement with the same comparison applied to simulations. Our results establishes the feasibility of the CLAMATO survey, which aims to obtain Ly$\alpha$ forest spectra for $\sim 1000$ SFGs over $\sim 1 \,\mathrm{deg}^2$ of the COSMOS field, in order to map out IGM large-scale structure at $\langle z \rangle \sim 2.3$ over a large volume $(100\,h^{-1}\mathrm{Mpc})^3$.

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

[Abridged] We characterise a massive proto-cluster at z=2.895 that we found in the COSMOS field using the spectroscopic sample of the VIMOS Ultra-Deep Survey (VUDS). This is one of the rare structures at z~3 not identified around AGNs or radio galaxies, so it is an ideal laboratory to study galaxy formation in dense environments. The structure comprises 12 galaxies with secure spectroscopic redshift in an area of 7'x8', in a z bin of Dz=0.016. The measured galaxy number overdensity is delta_g=12+/-2. This overdensity has total mass of M~8.1x10^(14)M_sun in a volume of 13x15x17 Mpc^3. Simulations indicate that such an overdensity at z~2.9 is a proto-cluster that will collapse in a cluster of total mass M~2.5x10^(15)M_sun at z=0. We compare the properties of the galaxies within the overdensity with a control sample at the same z but outside the overdensity. We did not find any statistically significant difference between the properties (stellar mass, SFR, sSFR, NUV-r, r-K) of the galaxies inside and outside the overdensity. The stacked spectrum of galaxies in the overdensity background shows a significant absorption feature at the wavelength of Lya redshifted at z=2.895 (lambda=4736 A), with a rest frame EW = 4+/- 1.4 A. Stacking only background galaxies without intervening sources at z~2.9 along their line of sight, we find that this absorption feature has a rest frame EW of 10.8+/-3.7 A, with a detection S/N of ~4. These EW values imply a high column density (N(HI)~3-20x10^(19)cm^(-2)), consistent with a scenario where such absorption is due to intervening cold gas streams, falling into the halo potential wells of the proto-cluster galaxies. However, we cannot exclude the hypothesis that this absorption is due to the diffuse gas within the overdensity.

We present the VIMOS Ultra Deep Survey (VUDS), a spectroscopic redshift survey of ~10.000 very faint galaxies to study the major phase of galaxy assembly 2<z<~6. The survey covers 1 deg^2 in 3 separate fields: COSMOS, ECDFS and VVDS-02h, with targets selection based on an inclusive combination of photometric redshifts and color properties. Spectra covering 3650<lambda<9350 A are obtained with VIMOS on the ESO-VLT with integration times of 14h. Here we present the survey strategy, the target selection, the data processing, as well as the redshift measurement process, emphasizing the specific methods adapted to this high redshift range. The spectra quality and redshift reliability are discussed, and we derive a completeness in redshift measurement of 91%, or 74% for the most reliable measurements, down to i_AB=25, and measurements are performed all the way down to i_AB=27. The redshift distribution of the main sample peaks at z=3-4 and extends over a large redshift range mainly in 2 < z < 6. At 3<z<5, the galaxies cover a large range of luminosities -23< M_U < -20.5, stellar mass 10^9 M_sun< M_star < 10^{11} M_sun, and star formation rates 1 M_sun/yr< SFR < 10^3 M_sun/yr. We discuss the spectral properties of galaxies using individual as well as stacked spectra. The comparison between spectroscopic and photometric redshifts as well as color selection demonstrate the effectiveness of our selection scheme. With ~6000 galaxies with reliable spectroscopic redshifts in 2<z<6 expected when complete, this survey is the largest at these redshifts and offers the opportunity for unprecedented studies of the star-forming galaxy population and its distribution in large scale structures during the major phase of galaxy assembly.

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

We describe the completed VIMOS VLT Deep Survey, and the final data release of 35016 galaxies and type-I AGN with measured spectroscopic redshifts up to redshift z~6.7, in areas 0.142 to 8.7 square degrees, and volumes from 0.5x10^6 to 2x10^7h^-3Mpc^3. We have selected samples of galaxies based solely on their i-band magnitude reaching i_{AB}=24.75. Spectra have been obtained with VIMOS on the ESO-VLT, integrating 0.75h, 4.5h and 18h for the Wide, Deep, and Ultra-Deep nested surveys. A total of 1263 galaxies have been re-observed independently within the VVDS, and from the VIPERS and MASSIV surveys. They are used to establish the redshift measurements reliability, to assess completeness, and to provide a weighting scheme taking into account the survey selection function. We describe the main properties of the VVDS samples, and the VVDS is compared to other spectroscopic surveys. In total we have obtained spectroscopic redshifts for 34594 galaxies, 422 type-I AGN, and 12430 Galactic stars. The survey has enabled to identify galaxies up to very high redshifts with 4669 redshifts in 1<=z_{spec}<=2, 561 in 2<=z_{spec}<=3 and 468 with z_{spec}>3, and specific populations like LAE have been identified out to z=6.62. We show that the VVDS occupies a unique place in the parameter space defined by area, depth, redshift coverage, and number of spectra. The VVDS provides a comprehensive survey of the distant universe, covering all epochs since z, or more than 12 Gyr of cosmic time, with a uniform selection, the largest such sample to date. A wealth of science results derived from the VVDS have shed new light on the evolution of galaxies and AGN, and their distribution in space, over this large cosmic time. A final public release of the complete VVDS spectroscopic redshift sample is available at http://cesam.lam.fr/vvds.

We use GOODS and CANDELS images to identify progenitors of massive (log M > 10 Msun) compact "early-type" galaxies (ETGs) at z~1.6. Since merging and accretion increase the size of the stellar component of galaxies, if the progenitors are among known star-forming galaxies, these must be compact themselves. We select candidate progenitors among compact Lyman-break galaxies at z~3 based on their mass, SFR and central stellar density and find that these account for a large fraction of, and possibly all, compact ETGs at z~1.6. We find that the average far-UV SED of the candidates is redder than that of the non-candidates, but the optical and mid-IR SED are the same, implying that the redder UV of the candidates is inconsistent with larger dust obscuration, and consistent with more evolved (aging) star-formation. This is in line with other evidence that compactness is a sensitive predictor of passivity among high-redshift massive galaxies. We also find that the light distribution of both the compact ETGs and their candidate progenitors does not show any extended "halos" surrounding the compact "core", both in individual images and in stacks. We argue that this is generally inconsistent with the morphology of merger remnants, even if gas-rich, as predicted by N-body simulations. This suggests that the compact ETGs formed via highly dissipative, mostly gaseous accretion of units whose stellar components are very small and undetected in the HST images, with their stellar mass assembling in-situ, and that they have not experienced any major merging until the epoch of observations at z~1.6.

We measure and analyse the redshift distribution N(z) of magnitude-selected samples using spectroscopic redshift measurement from the magnitude-selected VIMOS VLT Deep Survey (VVDS) with 17<iAB<24.75. We compute the N(z) and provide reference parametric fits for i band, J, H and Ks band magnitud limited samples. The N(z) of a sample with iAB<24 has a mean redshift z=0.92, with 8.2% of the galaxies with z>2. Down to iAB<24.75 the sample has a mean redshift z=1.15 and 17.1% of the galaxies are beyond z=2. The projected sky density is 2.07+/-0.12 gal/arcmin2 at 1.4<z<2.5 and KsAB<22.5, 1.72+/-0.15 gal/arcmin2 at 2.7<z<3.4 and 0.59+/-0.09 gal/arcmin2 at 3.4<z<4.5 brighter than iAB=24.75. Galaxies at z~3 identified from magnitude-selected samples are 1.5 to 3 times more numerous than when they are colour-colour selected. We demonstrate that colour-colour selected samples over 1.4<z<4.5 are strongly contaminated by galaxies at other redshifts. Semi-analytic models on the Millennium simulations under-predict the number of luminous star-forming galaxies at z\sim1.8-2, as well as over-predict the number of low-luminosity galaxies at z<0.8. Our study provides comprehensive galaxy number counts N(z) from galaxies with spectroscopic redshifts over a large redshift domain 0<z<5, a solid basis for the measurement of volume-complete quantities. Magnitude-selected surveys identify a higher number of galaxies at z>2 than in colour-colour selected samples, and we use the magnitude-selected VVDS to emphasize the large uncertainties associated to other surveys using colour or colour-colour selected samples. Our results further demonstrate that semi-analytical models on dark matter simulations have yet to find the right balance of physical processes and time-scales to properly reproduce a fundamental galaxy population property like the observed N(z).

The aim of this work is to identify HeII emitters at 2<z<4.6 and to constrain the source of the hard ionizing continuum that powers the HeII emission. We have assembled a sample of 277 galaxies with a high quality spectroscopic redshift at 2<z<4.6 from the VVDS survey, and we have identified 39 HeII1640A emitters. We study their spectral properties, measuring the fluxes, equivalent widths (EW) and FWHM for most relevant lines. About 10% of galaxies at z~3 show HeII in emission, with rest frame equivalent widths EW0~1-7A, equally distributed between galaxies with Lya in emission or in absorption. We find 11 high-quality HeII emitters with unresolved HeII line (FWHM_0<1200km/s), 13 high-quality emitters with broad He II emission (FWHM_0>1200km/s), 3 AGN, and an additional 12 possible HeII emitters. The properties of the individual broad emitters are in agreement with expectations from a W-R model. On the contrary, the properties of the narrow emitters are not compatible with such model, neither with predictions of gravitational cooling radiation produced by gas accretion. Rather, we find that the EW of the narrow HeII line emitters are in agreement with expectations for a PopIII star formation, if the episode of star formation is continuous, and we calculate that a PopIII SFR of 0.1-10 Mo yr-1 only is enough to sustain the observed HeII flux. We conclude that narrow HeII emitters are either powered by the ionizing flux from a stellar population rare at z~0 but much more common at z~3, or by PopIII star formation. As proposed by Tornatore et al. (2007), incomplete ISM mixing may leave some small pockets of pristine gas at the periphery of galaxies from which PopIII may form, even down to z~2 or lower. If this interpretation is correct, we measure at z~3 a SFRD in PopIII stars of 10^6Mo yr^-1 Mpc^-3 qualitatively comparable to the value predicted by Tornatore et al. (2007).

We have made a serendipitous discovery of a massive cD galaxy at z=1.096 in a candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy is the most distant cD galaxy confirmed to date. Ultra-deep HST/WFC3 images reveal an extended envelope starting from ~10 kpc and reaching ~70 kpc in radius along the semi-major axis. The spectral energy distributions indicate that both its inner component and outer envelope are composed of an old, passively-evolving stellar population. The cD galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The cD galaxy has a higher stellar mass surface density but a similar velocity dispersion to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of today's more massive cDs, its size and stellar mass have had to increase on average by factors of $3.4\pm1.1$ and $3.3\pm1.3$ over the past ~8 Gyrs, respectively. Such increases in size and stellar mass without being accompanied by significant increases in velocity dispersion are consistent with evolutionary scenarios driven by both major and minor dry mergers. If such cD envelopes originate from dry mergers, our discovery of even one example proves that some BCGs entered the dry merger phase at epochs earlier than z=1. Our data match theoretical models which predict that the continuance of dry mergers at z<1 can result in structures similar to those of massive cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the HUDF is essential to reveal such an extended cD envelope at z>1 and, yet, the HUDF covers only a minuscule region of sky. Adding that cDs are rare, Our serendipitous discovery hints that such cDs may be more common than expected. [Abridged]

The mass assembly of galaxies can proceed through different physical processes. Here we report on the spectroscopic identification of close physical pairs of galaxies at redshifts 2 <z < 4 and discuss the impact of major mergers in shaping galaxies at these early cosmological times. We aim to identify and characterize close physical pairs of galaxies destined to merge and use their properties to infer the contribution of merging processes to the early mass assembly of galaxies. We search for galaxy pairs with a transverse separation r_p < 25 h^{-1} kpc and a velocity difference Delta v < 500 km s^{-1} using the VIMOS VLT Deep Survey (VVDS) and early data from the VIMOS Ultra Deep Survey (VUDS) that comprise a sample of 1111 galaxies with spectroscopic redshifts measurements at redshifts 1.8 < z < 4 in the COSMOS, ECDFS, and VVDS-02h fields. We analyse their spectra and associated visible and near-infrared photometry to assess the main properties of merging galaxies that have an average stellar mass M_{star}=2.3 x 10^{10} M_{sun} at these redshifts. Using the 12 physical pairs found in our sample we obtain a first estimate of the merger fraction at these rdshifts, fmerg ~(15-20)%. The pair separations indicate that these pairs will merge within 1 Gyr, on average, with each producing a more massive galaxy by the time the cosmic star formation peaks at z~1-2. From the average mass ratio between galaxies in the pairs, the stellar mass of the resulting galaxy after merging will be ~60% larger than the most massive galaxy in the pair before merging. We conclude that major merging of galaxy pairs is on-going at 2<z<4 and significantly contributing to the major assembly phase of galaxies at this early epoch.

We study the evolution of the number density, as a function of the size, of passive early-type galaxies with a wide range of stellar masses 10^10<M*/Msun<10^11.5) from z~3 to z~1, exploiting the unique dataset available in the GOODS-South field, including the recently obtained WFC3 images as a part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). In particular, we select a sample of 107 massive (M*>10^10 M_sun), passive (SSFR<10^-2 Gyr^-1) and morphologically spheroidal galaxies at 1.2<z<3, taking advantage of the panchromatic dataset available for GOODS, including VLT, CFHT, Spitzer, Chandra and HST ACS+WFC3 data. We find that at 1<z<3 the passively evolving early-type galaxies are the reddest and most massive objects in the Universe, and we prove that a correlation between mass, morphology, color and star-formation activity is already in place at that epoch. We measure a significant evolution in the mass-size relation of passive early-type galaxies (ETGs) from z~3 to z~1, with galaxies growing on average by a factor of 2 in size in a 3 Gyr timescale only. We witness also an increase in the number density of passive ETGs of 50 times over the same time interval. We find that the first ETGs to form at z>2 are all compact or ultra-compact, while normal sized ETGs (meaning ETGs with sizes comparable to those of local counterparts of the same mass) are the most common ETGs only at z<1. The increase of the average size of ETGs at 0<z<1 is primarily driven by the appearance of new large ETGs rather than by the size increase of individual galaxies.

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

This paper studies the properties of kiloparsec-scale clumps in star-forming galaxies at z~2 through multi-wavelength broad band photometry. A sample of 40 clumps is identified through auto-detection and visual inspection from 10 galaxies with 1.5<z<2.5 in the Hubble Ultra Deep Field, where deep and high-resolution HST/WFC3 and ACS images enable us to resolve structures of z~2 galaxies down to kpc scale in the rest-frame UV and optical bands as well as to detect clumps toward the faint end. The physical properties of clumps are measured through fitting spatially resolved seven-band (BVizYJH) spectral energy distribution to models. On average, the clumps are blue and have similar median rest-frame UV--optical color as the diffuse components of their host galaxies, but the clumps have large scatter in their colors. Although the star formation rate (SFR)--stellar mass relation of galaxies is dominated by the diffuse components, clumps emerge as regions with enhanced specific SFRs, contributing individually ~10% and together ~50% of the SFR of the host galaxies. However, the contributions of clumps to the rest-frame UV/optical luminosity and stellar mass are smaller, typically a few percent individually and ~20% together. On average, clumps are younger by 0.2 dex and denser by a factor of 8 than diffuse components. Clump properties have obvious radial variations in the sense that central clumps are redder, older, more extincted, denser, and less active on forming stars than outskirts clumps. Our results are broadly consistent with a widely held view that clumps are formed through gravitational instability in gas-rich turbulent disks and would eventually migrate toward galactic centers and coalesce into bulges. Roughly 40% of the galaxies in our sample contain a massive clump that could be identified as a proto-bulge, which seems qualitatively consistent with such a bulge-formation scenario.

A new set of color selection criteria (VJL) analogous with the BzK method is designed to select both star-forming galaxies (SFGs) and passively-evolving galaxies (PEGs) at 2.3<z<3.5 by using rest-frame UV--optical (V-J vs. J-L) colors. The criteria are thoroughly tested with theoretical stellar population synthesis models and real galaxies with spectroscopic redshifts to evaluate their efficiency and contamination. We apply the well-tested VJL criteria to the HST/WFC3 Early Release Science field and study the physical properties of selected galaxies. The redshift distribution of selected SFGs peaks at z~2.7, slightly lower than that of Lyman Break Galaxies at z~3. Comparing the observed mid-infrared fluxes of selected galaxies with the prediction of pure stellar emission, we find that our VJL method is effective at selecting massive dusty SFGs that are missed by the Lyman Break Technique. About half of the star formation in massive (M_{star}>10^{10}M_{Sun}) galaxies at 2.3<z<3.5 is contributed by dusty (extinction E(B-V)>0.4) SFGs, which however, only account for ~20% of the number density of massive SFGs. We also use the mid-infrared fluxes to clean our PEG sample, and find that galaxy size can be used as a secondary criterion to effectively eliminate the contamination of dusty SFGs. The redshift distribution of the cleaned PEG sample peaks at z~2.5. We find 6 PEG candidates at z>3 and discuss possible methods to distinguish them from dusty contamination. We conclude that at least part of our candidates are real PEGs at z~3, implying that this type of galaxies began to form their stars at z>5. We measure the integrated stellar mass density of PEGs at z~2.5 and set constraints on it at z>3. We find that the integrated stellar mass density grows by at least about factor of 10 in 1 Gyr at 3<z<5 and by another factor of 10 in next 3.5 Gyr (1<z<3).

We present the analysis of a large sample of early-type galaxies (ETGs) at 0<z<3 aimed at tracing the cosmic evolution of their size and compare it with a model of pure dissipationless (dry) merging in the LambdaCDM framework. The effective radius R_e depends on stellar mass M as R_e(M) \propto M}^{alpha} with alpha ~ 0.5 at all redshifts. The redshift evolution of the mass- or SDSS-normalized size can be reproduced as \propto (1+z)^beta with beta ~ -1, with the most massive ETGs possibly showing the fastest evolutionary rate (beta ~ -1.4). This size evolution slows down significantly to beta ~ -0.6 if the ETGs at z>2 are removed from the sample, suggesting an accelerated increase of the typical sizes at z>2, especially for the ETGs with the largest masses. A pure dry merging LambdaCDM model is marginally consistent with the average size evolution at 0<z<1.7, but predicts descendants too compact for z>2 progenitor ETGs. This opens the crucial question on what physical mechanism can explain the accelerated evolution at z>2, or whether an unclear observational bias is partly responsible for that.

[Abridged] We investigate the global galaxy evolution over 12 Gyr (0.05<z<4.5), from the star formation rate density (SFRD), combining the VVDS Deep (17.5<=I<=24.0) and Ultra-Deep (23.00<=i<=24.75) surveys. We obtain a single homogeneous spectroscopic redshift sample, totalizing about 11000 galaxies. We estimate the rest-frame FUV luminosity function (LF) and luminosity density (LD), extract the dust attenuation of the FUV radiation using SED fitting, and derive the dust-corrected SFRD. We find a constant and flat faint-end slope alpha in the FUV LF at z<1.7. At z>1.7, we set alpha steepening with (1+z). The absolute magnitude M*_FUV brightens in the entire range 0<z<4.5, and at z>2 it is on average brighter than in the literature, while phi* is smaller. Our total LD shows a peak at z=2, present also when considering all sources of uncertainty. The SFRD history peaks as well at z=2. It rises by a factor of 6 during 2 Gyr (from z=4.5 to z=2), and then decreases by a factor of 12 during 10 Gyr down to z=0.05. This peak is mainly produced by a similar peak within the population of galaxies with -21.5<=M_FUV<=-19.5 mag. As times goes by, the total SFRD is dominated by fainter and fainter galaxies. The presence of a clear peak at z=2 and a fast rise at z>2 of the SFRD is compelling for models of galaxy formation. The mean dust attenuation A_FUV of the global galaxy population rises by 1 mag during 2 Gyr from z=4.5 to z=2, reaches its maximum at z=1 (A_FUV=2.2 mag), and then decreases by 1.1 mag during 7 Gyr down to z=0. The dust attenuation maximum is reached 2 Gyr after the SFRD peak, implying a contribution from the intermediate-mass stars to the dust production at z<2.

(Abridged) We analyze the environments and galactic properties (morphologies and star-formation histories) of a sample of 153 close kinematic pairs in the redshift range 0.2 < z < 1 identified in the zCOSMOS-bright 10k spectroscopic sample of galaxies. Correcting for projection effects, the fraction of close kinematic pairs is three times higher in the top density quartile than in the lowest one. This translates to a three times higher merger rate because the merger timescales are shown, from mock catalogues based on the Millennium simulation, to be largely independent of environment once the same corrections for projection is applied. We then examine the morphologies and stellar populations of galaxies in the pairs, comparing them to control samples that are carefully matched in environment so as to remove as much as possible the well-known effects of environment on the properties of the parent population of galaxies. Once the environment is properly taken into account in this way, we find that the early-late morphology mix is the same as for the parent population, but that the fraction of irregular galaxies is boosted by 50-75%, with a disproportionate increase in the number of irregular-irregular pairs (factor of 4-8 times), due to the disturbance of disk galaxies. Future dry-mergers, involving elliptical galaxies comprise less than 5% of all close kinematic pairs. In the closest pairs, there is a boost in the specific star-formation rates of star-forming galaxies of a factor of 2-4, and there is also evidence for an increased incidence of post star-burst galaxies. Although significant for the galaxies involved, the "excess" star-formation associated with pairs represents only about 5% of the integrated star-formation activity in the parent sample. Although most pair galaxies are in dense environments, the effects of interaction appear to be largest in the lower density environments.

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

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

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

Our goal is to develop a new and reliable statistical method to classify galaxies from large surveys. We probe the reliability of the method by comparing it with a three-dimensional classification cube, using the same set of spectral, photometric and morphological parameters.We applied two different methods of classification to a sample of galaxies extracted from the zCOSMOS redshift survey, in the redshift range 0.5 < z < 1.3. The first method is the combination of three independent classification schemes, while the second method exploits an entirely new approach based on statistical analyses like Principal Component Analysis (PCA) and Unsupervised Fuzzy Partition (UFP) clustering method. The PCA+UFP method has been applied also to a lower redshift sample (z < 0.5), exploiting the same set of data but the spectral ones, replaced by the equivalent width of H$\alpha$. The comparison between the two methods shows fairly good agreement on the definition on the two main clusters, the early-type and the late-type galaxies ones. Our PCA-UFP method of classification is robust, flexible and capable of identifying the two main populations of galaxies as well as the intermediate population. The intermediate galaxy population shows many of the properties of the green valley galaxies, and constitutes a more coherent and homogeneous population. The fairly large redshift range of the studied sample allows us to behold the downsizing effect: galaxies with masses of the order of $3\cdot 10^{10}$ Msun mainly are found in transition from the late type to the early type group at $z>0.5$, while galaxies with lower masses - of the order of $10^{10}$ Msun - are in transition at later epochs; galaxies with $M <10^{10}$ Msun did not begin their transition yet, while galaxies with very large masses ($M > 5\cdot 10^{10}$ Msun) mostly completed their transition before $z\sim 1$.

We report on the evolution of the number density and size of early-type galaxies from z~2 to z~0. We select a sample of 563 massive (M>10^{10} Msun), passively evolving (SSFR<10^{-2} Gyr^{-1}) and morphologically spheroidal galaxies at 0<z<2.5, using the panchromatic photometry and spectroscopic redshifts available in the GOODS fields. We combine ACS and WFC3 HST images to study the morphology of our galaxies in their optical rest-frame in the whole 0<z<2.5 range. We find that throughout the explored redshift range the passive galaxies selected with our criteria have weak morphological K-correction, with size being slightly smaller in the optical than in the UV rest-frame (by ~20 and ~10 at z>1.2 and z<1.2, respectively). We measure a significant evolution of the mass-size relation of early-type galaxies, with the fractional increment that is almost independent on the stellar mass. Early-type galaxies (ETGs) formed at z>1 appear to be preferentially small, and the evolution of the mass-size relation at z<1 is driven by both the continuous size growth of the compact galaxies and the appearance of new ETGs with large sizes. We also find that the number density of all passive early-type galaxies increases rapidly, by a factor of 5, from z~2 to z~1, and then more mildly by another factor of 1.5 from z~1 to z~0. We interpret these results as the evidence that the bulk of the ETGs are formed at 1<z<3 through a mechanism that leaves very compact remnants. At z<1 the compact ETGs grow gradually in size, becoming normal size galaxies, and at the same time new ETGs with normal-large sizes are formed.

We report the discovery of large amounts of cold (T ~ 10^4 K), chemically young gas in an overdensity of galaxies at redshift z ~ 1.6 in the Great Observatories Origins Deep Survey southern field (GOODS-S). The gas is identified thanks to the ultra-strong Mg II absorption features it imprints in the rest-frame UV spectra of galaxies in the background of the overdensity. There is no evidence that the optically-thick gas is part of any massive galaxy (i.e. M_star > 4x10^9 M_sun), but rather is associated with the overdensity; less massive and fainter galaxies (25.5 < z_850 < 27.5 mag) have too large an impact parameter to be causing ultra-strong absorption systems, based on our knowledge of such systems. The lack of corresponding Fe II absorption features, not detected even in co-added spectra, suggests that the gas is chemically more pristine than the ISM and outflows of star-forming galaxies at similar redshift, including those in the overdensity itself, and comparable to the most metal-poor stars in the Milky Way halo. A crude estimate of the projected covering factor of the high-column density gas (N_H >~ 10^20 cm-2) based on the observed fraction of galaxies with ultra-strong absorbers is C_F ~ 0.04. A broad, continuum absorption profile extending to the red of the interstellar Mg II absorption line by <~ 2000 km/s is possibly detected in two independent co-added spectra of galaxies of the overdensity, consistent with a large-scale infall motion of the gas onto the overdensity and its galaxies. Overall, these findings provides the first tentative evidence of accretion of cold, chemically young gas onto galaxies at high redshift, possibly feeding their star formation activity. The fact that the galaxies are members of a large structure, as opposed to field galaxies, might play a significant role in our ability to detect the accreting gas.

We report the detection of color gradients in six massive (stellar mass > 10^{10} M_{sun}) and passively evolving (specific SFR < 10^{-11}/yr) galaxies at redshift 1.3<z<2.5 identified in the HUDF using HST ACS and WFC3/IR images. After matching different PSFs, we obtain color maps and multi-band optical/near-IR photometry (BVizYJH) in concentric annuli, from the smallest resolved radial (~1.7 kpc) up to several times the H-band effective radius. We find that the inner regions of these galaxies have redder rest-frame UV--optical colors than the outer parts. The slopes of the color gradients mildly depend on the overall dust obscuration and rest-frame (U-V) color, with more obscured or redder galaxies having steeper color gradients. The z~2 color gradients are also steeper than those of local early-types. The gradient of a single parameter (age, extinction or metallicity) cannot fully explain the observed color gradients. Fitting spatially resolved HST seven-band photometry to stellar population synthesis models, we find that, regardless of assumptions for metallicity gradient, the redder inner regions of the galaxies have slightly higher dust obscuration than the bluer outer regions, although the magnitude depends on the assumed extinction law. The derived age gradient depends on the assumptions for metallicity gradient. We discuss the implications of a number of assumptions for metallicity gradient on the formation and evolution of these galaxies. We find that the evolution of the mass--size relationship from z~2 to z~0 cannot be driven by in--situ extended star formation, implying that accretion or merger is mostly responsible for the evolution. The lack of a correlation between color gradient and stellar mass argues against the metallicity gradient predicted by the monolithic collapse, which would require significant major mergers to evolve into the one observed at z~0. (Abridged)

The contribution of major mergers to galaxy mass assembly along cosmic time is an important ingredient to the galaxy evolution scenario. We aim to measure the evolution of the merger rate for both luminosity/mass selected galaxy samples and investigate its dependence with the local environment. We use a sample of 10644 spectroscopically observed galaxies from the zCOSMOS redshift survey to identify pairs of galaxies destined to merge, using only pairs for which the velocity difference and projected separation of both components with a confirmed spectroscopic redshift indicate a high probability of merging. We have identified 263 spectroscopically confirmed pairs with r_p^{max} = 100 h^{-1} kpc. We find that the density of mergers depends on luminosity/mass, being higher for fainter/less massive galaxies, while the number of mergers a galaxy will experience does not depends significantly on its intrinsic luminosity but rather on its stellar mass. We find that the pair fraction and merger rate increase with local galaxy density, a property observed up to redshift z=1. We find that the dependence of the merger rate on the luminosity or mass of galaxies is already present up to redshifts z=1, and that the evolution of the volumetric merger rate of bright (massive) galaxies is relatively flat with redshift with a mean value of 3*10^{-4} (8*10^{-5} respectively) mergers h^3 Mpc^{-3} Gyr^{-1}. The dependence of the merger rate with environment indicates that dense environments favors major merger events as can be expected from the hierarchical scenario. The environment therefore has a direct impact in shapping-up the mass function and its evolution therefore plays an important role on the mass growth of galaxies along cosmic time.

We measure the spatial clustering of galaxies as a function of their morphological type at z~0.8, for the first time in a deep redshift survey with full morphological information. This is obtained by combining high-resolution HST imaging and VLT spectroscopy for about 8,500 galaxies to I_AB=22.5 with accurate spectroscopic redshifts from the zCOSMOS-Bright redshift survey. At this epoch, early-type galaxies already show a significantly stronger clustering than late-type galaxies on all probed scales. A comparison to the SDSS at z~0.1, shows that the relative clustering strength between early and late morphological classes tends to increase with cosmic time at small separations, while on large scales it shows no significant evolution since z~0.8. This suggests that most early-type galaxies had already formed in intermediate and dense environments at this epoch. Our results are consistent with a picture in which the relative clustering of different morphological types between z~1 and z~0, reflects the evolving role of environment in the morphological transformation of galaxies, on top of the global mass-driven evolution.

Aims. The aim of this work is to study the contribution of the Ly-a emitters (LAE) to the star formation rate density (SFRD) of the Universe in the interval 2<z<6.6. Methods. We assembled a sample of 217 LAE from the Vimos-VLT Deep Survey (VVDS) with secure spectroscopic redshifts in the redshift range 2 < z < 6.62 and fluxes down to F=1.5x10^18 erg/s/cm^2. 133 LAE are serendipitous identifications in the 22 arcmin^2 total slit area surveyed with the VVDS-Deep and the 3.3 arcmin^2 from the VVDS Ultra-Deep survey, and 84 are targeted identifications in the 0.62 deg^2 surveyed with the VVDS-DEEP and 0.16 deg^2 from the Ultra-Deep survey. Among the serendipitous targets we estimate that 90% of the emission lines are most probably Ly-a, while the remaining 10% could be either [OII]3727 or Ly-a. We computed the LF and derived the SFRD from LAE at these redshifts. Results. The VVDS-LAE sample reaches faint line fluxes F(Lya) = 1.5x1^18 erg/s/cm^2 (corresponding to L(Lya)=10^41 erg/s at z~3) enabling the faint end slope of the luminosity function to be constrained to a=-1.6+-0.12 at redshift z~2.5 and to a=-1.78+0.1-0.12 at z=4, placing on firm statistical grounds trends found in previous LAE studies, and indicating that sub-L* LAE contribute significantly to the SFRD. The projected number density and volume density of faint LAE in 2<z<6.6 with F>1.5x10^18 erg/s/cm^2 are 33 galaxies/arcmin^2 and 4x10^-2 Mpc^-3, respectively. We find that the the observed luminosity function of LAE does not evolve from z=2 to z=6. This implies that, after correction for the redshift-dependent IGM absorption, the intrinsic LF must have evolved significantly over 3 Gyr. The SFRD from LAE contributes to about 20% of the SFRD at z =2-3, while the LAE appear to be the dominant source of star formation producing ionizing photons in the early universe z>5-6, becoming equivalent to that of Lyman Break galaxies.

We study the impact of the environment on the evolution of galaxies in the zCOSMOS 10k sample in the redshift range 0.1<z<1.0 over an area of ~1.5 deg2. The considered sample of secure spectroscopic redshifts contains about 8500 galaxies, with their stellar masses estimated by SED fitting of the multiwavelength optical to NIR photometry. The evolution of the galaxy stellar mass function (GSMF) in high and low density regions provides a tool to study the mass assembly evolution in different environments; moreover, the contributions to the GSMF from different galaxy types, as defined by their SEDs and their morphologies, can be quantified. At redshift z~1, the GSMF is only slightly dependent on environment, but at lower redshifts the shapes of the GSMFs in high- and low-density environments become extremely different, with high density regions exhibiting a marked bimodality. As a result, we infer that galaxy evolution depends on both the stellar mass and the environment, the latter setting the probability of a galaxy to have a given mass: all the galaxy properties related to the stellar mass show a dependence on environment, reflecting the difference observed in the mass functions. The shapes of the GSMFs of early- and late-type galaxies are almost identical for the extremes of the density contrast we consider. The evolution toward z=0 of the mass at which the early- and late-type GSMFs match each other is more rapid in high density environments. The comparison of the observed GSMFs to the same quantities derived from a set of mock catalogues shows that blue galaxies in sparse environments are overproduced in the semi-analytical models at intermediate and high masses, because of a deficit of star formation suppression, while at z<0.5 an excess of red galaxies is present in dense environments at intermediate and low masses, because of the overquenching of satellites. ABRIDGED

We explore the inter-relationships between mass, star-formation rate and environment in the SDSS, zCOSMOS and other surveys. The differential effects of mass and environment are completely separable to z ~ 1, indicating that two distinct processes are operating, "mass-quenching" and "environment-quenching". Environment-quenching, at fixed over-density, evidently does not change with epoch to z ~ 1, suggesting that it occurs as large-scale structure develops in the Universe. The observed constancy of the mass-function shape for star-forming galaxies, demands that the mass-quenching of galaxies around and above M*, must be proportional to their star-formation rates at all z < 2. We postulate that this simple mass-quenching law also holds over a much broader range of stellar mass and epoch. These two simple quenching processes, plus some additional quenching due to merging, then naturally produce (a) a quasi-static Schechter mass function for star-forming galaxies with a value of M* that is set by the proportionality between the star-formation and mass-quenching rates, (b) a double Schechter function for passive galaxies with two components: the dominant one is produced by mass-quenching and has exactly the same M* as the star-forming galaxies but an alpha shallower by +1, while the other is produced by environment effects and has the same M* and alpha as the star-forming galaxies, and is larger in high density environments. Subsequent merging of quenched galaxies modifies these predictions somewhat in the denser environments, slightly increasing M* and making alpha more negative. All of these detailed quantitative relationships between the Schechter parameters are indeed seen in the SDSS, lending strong support to our simple empirically-based model. The model naturally produces for passive galaxies the "anti-hierarchical" run of mean ages and alpha-element abundances with mass.

We present the analysis of the U-V rest-frame color distribution and some spectral features as a function of mass and environment for two sample of early-type galaxies up to z=1 extracted from the zCOSMOS spectroscopic survey. The first sample ("red galaxies") is defined with a photometric classification, while the second ("ETGs") by combining morphological, photometric, and spectroscopic properties to obtain a more reliable sample. We find that the color distribution of red galaxies is not strongly dependent on environment for all mass bins, with galaxies in overdense regions redder than galaxies in underdense regions with a difference of 0.027\pm0.008 mag. The dependence on mass is far more significant, with average colors of massive galaxies redder by 0.093\pm0.007 mag than low-mass galaxies throughout the entire redshift range. We study the color-mass relation, finding a mean slope 0.12\pm0.005, while the color-environment relation is flatter, with a slope always smaller than 0.04. The spectral analysis that we perform on our ETGs sample is in good agreement with our photometric results: we find for D4000 a dependence on mass between high and low-mass galaxies, and a much weaker dependence on environment (respectively a difference of of 0.11\pm0.02 and of 0.05\pm0.02); for the equivalent width of H{\delta}we measure a difference of 0.28\pm0.08 {\AA}across the same mass range and no significant dependence on environment.By analyzing the lookback time of early-type galaxies, we support the possibility of a downsizing scenario, in which massive galaxies with a stronger D4000 and an almost constant equivalent width of $H\delta$ formed their mass at higher redshift than lower mass ones. We also conclude that the main driver of galaxy evolution is the galaxy mass, the environment playing a subdominant role.

[Abridged] With the first 10000 spectra of the flux limited zCOSMOS sample (I<=22.5) we study the evolution of environmental effects on galaxy properties since z=1.0, and disentangle the dependence among galaxy colour, stellar mass and local density (3D local density contrast `delta', computed with the 5th nearest neighbour approach). We confirm that within a luminosity-limited sample (M_B<=-20.5-z) the fraction of red (U-B>=1) galaxies 'f_red' depends on delta at least up to z=1, with red galaxies residing mainly in high densities. This trend weakens for increasing z, and it is mirrored by the behaviour of the fraction of galaxies with D4000A break >=1.4. We also find that up to z=1 the fraction of galaxies with log(EW[OII]) >=1.15 is higher for lower delta, and also this dependence weakens for increasing z. Given the triple dependence among galaxy colours, stellar mass and delta, the colour-delta relation found in the luminosity-selected sample can be due to the broad range of stellar masses. Thus, we fix the stellar mass and we find that in this case the colour-delta relation is flat up to z=1 for galaxies with log(M/M_sun)>=10.7. This means that for these masses the colour-delta relation found in a luminosity-selected sample is the result of the combined colour-mass and mass-delta relations. In contrast, we find that for 0.1<=z<=0.5 and log(M/M_sun)<=10.7 'f_red' depends on delta even at fixed mass. In these mass and z ranges, environment affects directly also galaxy colours. We suggest a scenario in which the colour depends primarily on stellar mass, but for relatively low mass galaxies the local density modulates this dependence. These galaxies formed more recently, in an epoch when evolved structures were already in place, and their longer SFH allowed environment-driven physical processes to operate during longer periods of time.

We present a galaxy group catalogue spanning the redshift range 0.1 <~ z <~ 1 in the ~1.7 deg^2 COSMOS field, based on the first ~10,000 zCOSMOS spectra. The performance of both the Friends-of-Friends (FOF) and Voronoi-Delaunay-Method (VDM) approaches to group identification has been extensively explored and compared using realistic mock catalogues. We find that the performance improves substantially if groups are found by progressively optimizing the group-finding parameters for successively smaller groups, and that the highest fidelity catalogue, in terms of completeness and purity, is obtained by combining the independently created FOF and VDM catalogues. The final completeness and purity of this catalogue, both in terms of the groups and of individual members, compares favorably with recent results in the literature. The current group catalogue contains 102 groups with N >= 5 spectroscopically confirmed members, with a further ~700 groups with 2 <= N <= 4. Most of the groups can be assigned a velocity dispersion and a dark-matter mass derived from the mock catalogues, with quantifiable uncertainties. The fraction of zCOSMOS galaxies in groups is about 25% at low redshift and decreases toward ~15% at z ~ 0.8. The zCOSMOS group catalogue is broadly consistent with that expected from the semi-analytic evolution model underlying the mock catalogues. Not least, we show that the number density of groups with a given intrinsic richness increases from redshift z ~ 0.8 to the present, consistent with the hierarchical growth of structure.

We investigate how the shape of the galaxy two-point correlation function as measured in the zCOSMOS survey depends on local environment, quantified in terms of the density contrast on scales of 5 Mpc/h. We show that the flat shape previously observed at redshifts between z=0.6 and z=1 can be explained by this volume being simply 10% over-abundant in high-density environments, with respect to a Universal density probability distribution function. When galaxies corresponding to the top 10% tail of the distribution are excluded, the measured w_p(r_p) steepens and becomes indistinguishable from LCDM predictions on all scales. This is the same effect recognised by Abbas & Sheth in the SDSS data at z~0 and explained as a natural consequence of halo-environment correlations in a hierarchical scenario. Galaxies living in high-density regions trace dark matter halos with typically higher masses, which are more correlated. If the density probability distribution function of the sample is particularly rich in high-density regions because of the variance introduced by its finite size, this produces a distorted two-point correlation function. We argue that this is the dominant effect responsible for the observed "peculiar" clustering in the COSMOS field.

We analyze the morphological properties of a large sample of 1503 70 micron selected galaxies in the COSMOS field spanning the redshift range 0.01<z< 3.5 with a median redshift of 0.5 and an infrared luminosity range of 10^8<L_IR<10^14L_sun with a median luminosity of 10^11.4 L_sun. In general these galaxies are massive, with a stellar mass range of 10^10-10^12 M_sun, and luminous, with -25<M_K<-20. We find a strong correlation between the fraction of major mergers and L_IR, with the fraction at the highest luminosity being up to 50%. We also find that the fraction of spirals drops dramatically with L_IR. Minor mergers likely play a role in boosting the infrared luminosity for sources with low luminosities. The precise fraction of mergers in any given L_IR bin varies by redshift due to sources at z>1 being difficult to classify and subject to the effects of band pass shifting, therefore, these numbers can only be considered lower limits. At z<1, where the morphological classifications are most robust, major mergers clearly dominate the ULIRG population (50-80%) and are important for the LIRG population (25-40%). At z>1 the fraction of major mergers is at least 30-40% for ULIRGs. Although the general morphological trends agree with what has been observed for local (U)LIRGs, the fraction of major mergers is slightly lower than seen locally. This is in part due to the difficulty of identifying merger signatures at high redshift. We argue that given the number of major gas-rich mergers observed and the relatively short timescale that they would be observable in the (U)LIRG phase that it is plausible for the observed red sequence of massive ellipticals (<10^12 M_sun) to have been formed entirely by gas-rich major mergers.

We present the Galaxy Stellar Mass Function (MF) up to z~1 from the zCOSMOS-bright 10k spectroscopic sample. We investigate the total MF and the contribution of ETGs and LTGs, defined by different criteria (SED, morphology or star formation). We unveil a galaxy bimodality in the global MF, better represented by 2 Schechter functions dominated by ETGs and LTGs, respectively. For the global population we confirm that low-mass galaxies number density increases later and faster than for massive galaxies. We find that the MF evolution at intermediate-low values of Mstar (logM<10.6) is mostly explained by the growth in stellar mass driven by smoothly decreasing star formation activities. The low residual evolution is consistent with ~0.16 merger per galaxy per Gyr (of which fewer than 0.1 are major). We find that ETGs increase in number density with cosmic time faster for decreasing Mstar, with a median "building redshift" increasing with mass, in contrast with hierarchical models. For LTGs we find that the number density of blue or spiral galaxies remains almost constant from z~1. Instead, the most extreme population of active star forming galaxies is rapidly decreasing in number density. We suggest a transformation from blue active spirals of intermediate mass into blue quiescent and successively (1-2 Gyr after) into red passive types. The complete morphological transformation into red spheroidals, required longer time-scales or follows after 1-2 Gyr. A continuous replacement of blue galaxies is expected by low-mass active spirals growing in stellar mass. We estimate that on average ~25% of blue galaxies is transforming into red per Gyr for logM<11. We conclude that the build-up of galaxies and ETGs follows the same downsizing trend with mass as the formation of their stars, converse to the trend predicted by current SAMs. We expect a negligible evolution of the global Galaxy Baryonic MF.

We present near-IR images, obtained with the Hubble Space Telescope (HST) and the WFC3/IR camera, of six passive and massive galaxies at redshift 1.3<z<2.4 (SSFR<10^{-2} Gyr^{-1}; stellar mass M~10^{11} M_{sun}), selected from the Great Observatories Origins Deep Survey (GOODS). These images, which have a spatial resolution of ~1.5 kpc, provide the deepest view of the optical rest-frame morphology of such systems to date. We find that the light profile of these galaxies is regular and well described by a Sersic model with index typical of today's spheroids. Their size, however, is generally much smaller than today's early types of similar stellar mass, with four out of six galaxies having r_e ~ 1 kpc or less, in quantitative agreement with previous similar measures made at rest-frame UV wavelengths. The images reach limiting surface brightness mu~26.5 mag arcsec^{-2} in the F160W bandpass; yet, there is no evidence of a faint halo in the galaxies of our sample, even in their stacked image. We also find that these galaxies have very weak "morphological k-correction" between the rest-frame UV (from the ACS z-band), and the rest--frame optical (WFC3 H-band): the Sersic index, physical size and overall morphology are independent or only mildly dependent on the wavelength, within the errors.

The identities of the main processes triggering and quenching star-formation in galaxies remain unclear. A key stage in evolution, however, appears to be represented by post-starburst galaxies. To investigate their impact on galaxy evolution, we initiated a multiwavelength study of galaxies with k+a spectral features in the COSMOS field. We examine a mass-selected sample of k+a galaxies at z=0.48-1.2 using the spectroscopic zCOSMOS sample. K+a galaxies occupy the brightest tail of the luminosity distribution. They are as massive as quiescent galaxies and populate the green valley in the colour versus luminosity (or stellar mass) distribution. A small percentage (<8%) of these galaxies have radio and/or X-ray counterparts (implying an upper limit to the SFR of ~8Msun/yr). Over the entire redshift range explored, the class of k+a galaxies is morphologically a heterogeneous population with a similar incidence of bulge-dominated and disky galaxies. This distribution does not vary with the strength of the Hdelta absorption line but instead with stellar mass in a way reminiscent of the well-known mass-morphology relation. Although k+a galaxies are also found in underdense regions, they appear to reside typically in a similarly rich environment as quiescent galaxies on a physical scale of ~2-8Mpc, and in groups they show a morphological early-to-late type ratio similar to the quiescent galaxy class. With the current data set, we do not find evidence of statistical significant evolution in either the number/mass density of k+a galaxies at intermediate redshift with respect to the local values, or the spectral properties. Those galaxies, which are affected by a sudden quenching of their star-formation activity, may increase the stellar mass of the red-sequence by up to a non-negligible level of ~10%.

We follow the galaxy stellar mass assembly by morphological and spectral type in the COSMOS 2-deg^2 field. We derive the stellar mass functions and stellar mass densities from z=2 to z=0.2 using 196,000 galaxies selected at F(3.6 micron) > 1 microJy with accurate photometric redshifts (sigma_((zp-zs)/(1+zs))=0.008 at i<22.5). Using a spectral classification, we find that z~1 is an epoch of transition in the stellar mass assembly of quiescent galaxies. Their stellar mass density increases by 1.1 dex between z=1.5-2 and z=0.8-1 (Delta t ~2.5 Gyr), but only by 0.3 dex between z=0.8-1 and z~0.1 (Delta t ~ 6 Gyr). Then, we add the morphological information and find that 80-90% of the massive quiescent galaxies (log(M)~11) have an elliptical morphology at z<0.8. Therefore, a dominant mechanism links the shutdown of star formation and the acquisition of an elliptical morphology in massive galaxies. Still, a significant fraction of quiescent galaxies present a Spi/Irr morphology at low mass (40-60% at log(M)~9.5), but this fraction is smaller than predicted by semi-analytical models using a ``halo quenching'' recipe. We also analyze the evolution of star-forming galaxies and split them into ``intermediate activity'' and ``high activity'' galaxies. We find that the most massive ``high activity'' galaxies end their high star formation rate phase first. Finally, the space density of massive star-forming galaxies becomes lower than the space density of massive elliptical galaxies at z<1. As a consequence, the rate of ``wet mergers'' involved in the formation of the most massive ellipticals must decline very rapidly at z<1, which could explain the observed slow down in the assembly of these quiescent and massive sources.

We present a catalog of 213 type-2 AGN selected from the zCOSMOS survey. The selected sample covers a wide redshift range (0.15<z<0.92) and is deeper than any other previous study, encompassing the luminosity range 10^{5.5} < Lsun< L[OIII] < 10^{9.1} Lsun. We explore the intrinsic properties of these AGN and the relation to their X-ray emission (derived from the XMM-COSMOS observations). We study their evolution by computing the [OIII]5007A line luminosity function (LF) and we constrain the fraction of obscured AGN as a function of luminosity and redshift. The sample was selected on the basis of the optical emission line ratios, after applying a cut to the signal-to-noise ratio (S/N) of the relevant lines. We used the standard diagnostic diagrams [OIII]/Hbeta versus [NII]/Halpha and ([OIII]/Hbeta versus [SII]/Halpha) to isolate AGN in the redshift range 0.15<z<0.45 and the diagnostic diagram [OIII]/Hbeta versus [OII]/Hbeta to extend the selection to higher redshift (0.5<z<0.92). Combining our sample with one drawn from SDSS, we found that the best description of the evolution of type-2 AGN is a luminosity-dependent density evolution model. Moreover, using the type-1 AGN LF we were able to constrain the fraction of type-2 AGN to the total (type-1 + type-2) AGN population. We found that the type-2 fraction decreases with luminosity, in agreement with the most recent results, and shows signs of a slight increase with redshift. However, the trend with luminosity is visible only after combining the SDSS+zCOSMOS samples. From the COSMOS data points alone, the type-2 fraction seems to be quite constant with luminosity.

We investigate the properties and the environment of radio sources with optical counterpart from the combined VLA-COSMOS and zCOSMOS samples. The advantage of this sample is the availability of optical spectroscopic information, high quality redshifts, and accurate density determination. By comparing the star formation rates estimated from the optical spectral energy distribution with those based on the radio luminosity, we divide the radio sources in three families, passive AGN, non-passive AGN and star forming galaxies. These families occupy specific regions of the 8.0-4.5 $\mu$m infrared color--specific star formation plane, from which we extract the corresponding control samples. Only the passive AGN have a significantly different environment distribution from their control sample. The fraction of radio-loud passive AGN increases from ~2% in underdense regions to ~15% for overdensities (1+delta) greater than 10. This trend is also present as a function of richness of the groups hosting the radio sources. Passive AGN in overdensities tend to have higher radio luminosities than those in lower density environments. Since the black hole mass distribution is similar in both environments, we speculate that, for low radio luminosities, the radio emission is controlled (through fuel disponibility or confinement of radio jet by local gas pressure) by the interstellar medium of the host galaxy, while in other cases it is determined by the structure (group or cluster) in which the galaxy resides.

We have used the zCOSMOS-bright 10k sample to identify 3244 Spitzer/MIPS 24-micron-selected galaxies with 0.06< S(24um)< 0.50 mJy and I(AB)<22.5, over 1.5 deg^2 of the COSMOS field, and studied different spectral properties, depending on redshift. At 0.2<z<0.3, we found that different reddening laws of common use in the literature explain the dust extinction properties of around 80% of our infrared (IR) sources, within the error bars. For up to 16% of objects, instead, the Halpha/Hbeta ratios are too high for their IR/UV attenuations, which is probably a consequence of inhomogenous dust distributions. In only a few of our galaxies at 0.2<z<0.3 the IR emission could be mainly produced by dust heated by old rather than young stars. Besides, the line ratios of ~22% of our galaxies suggest that they might be star-formation/nuclear-activity composite systems. At 0.5<z<0.7, we estimated galaxy metallicities for 301 galaxies: at least 12% of them are securely below the upper-branch mass-metallicity trend, which is consistent with the local relation. Finally, we performed a combined analysis of the Hdelta equivalent-width versus Dn(4000) diagram for 1722 faint and bright 24um galaxies at 0.6<z<1.0, spanning two decades in mid-IR luminosity. We found that, while secondary bursts of star formation are necessary to explain the position of the most luminous IR galaxies in that diagram, quiescent, exponentially-declining star formation histories can well reproduce the spectral properties of ~40% of the less luminous sources. Our results suggest a transition in the possible modes of star formation at total IR luminosities L(TIR)=(3 +/-2)x10^11 Lsun.

We use the overdensity field reconstructed in the volume of the COSMOS area to study the nonlinear biasing of the zCOSMOS galaxies. The galaxy overdensity field is reconstructed using the current sample of ~8500 accurate zCOSMOS redshifts at I(AB)<22.5 out to z~1 on scales R from 8 to 12 Mpc/h. By comparing the probability distribution function (PDF) of galaxy density contrast delta_g to the lognormal approximation of the PDF of the mass density contrast delta, we obtain the mean biasing function b(delta,z,R) between the galaxy and matter overdensity field and its second moments b(hat) and b(tilde) up to z~1. Over the redshift interval 0.4<z<1 the conditional mean function <delta_g|delta> = b(delta,z,R) delta is of the following characteristic shape. The function vanishes in the most underdense regions and then sharply rises in a nonlinear way towards the mean densities. <delta_g|delta> is almost a linear tracer of the matter in the overdense regions, up to the most overdense regions in which it is nonlinear again and the local effective slope of <delta_g|delta> vs. delta is smaller than unity. The <delta_g|delta> function is evolving only slightly over the redshift interval 0.4<z<1. The linear biasing parameter increases from b(hat)=1.24+/-0.11 at z=0.4 to b(hat)=1.64+/-0.15 at z=1 for the M_B<-20-z sample of galaxies. b(hat) does not show any dependence on the smoothing scale from 8 to 12 Mpc/h, but increases with luminosity. The measured nonlinearity parameter b(tilde)/b(hat) is of the order of a few percent (but it can be consistent with 0) and it does not change with redshift, the smoothing scale or the luminosity. By matching the linear bias of galaxies to the halo bias, we infer that the M_B<-20-z galaxies reside in dark matter haloes with a characteristic mass of about 3-6 x 10^12 Msol, depending on the halo bias fit.

(Abridged) We studied the evolution in the B band luminosity function to z~1 in the zCOSMOS 10k sample, for which both accurate galaxy classifications and a detailed description of the local density field are available. The global LF exhibits a brightening of ~0.7 mag in M* from z~0.2 to z~0.9. At low z, late types dominate at faint magnitudes, while the bright end is populated mainly by early types. At higher z, late-type galaxies evolve significantly and, at z~1, the contribution from the various types to the bright end of the LF is comparable. The evolution for early types is in both luminosity and normalization. A similar behaviour is exhibited by late types, but with an opposite trend for the normalization. Studying the role of the environment, we find that the global LF of galaxies in overdense regions has always a brighter M* and a flatter slope. In low density environments, the main contribution to the LF is from blue galaxies, while for high density environments there is an important contribution from red galaxies to the bright end. The differences between the global LF in the two environments are not due to only a difference in the relative numbers of red and blue galaxies, but also to their relative luminosity distributions: the value of M* for both types in underdense regions is always fainter than in overdense environments. The "specular" evolution of late- and early-type galaxies is consistent with a scenario where a part of blue galaxies is transformed in red galaxies with increasing cosmic time, without significant changes in the fraction of intermediate-type galaxies. The bulk of this tranformation in overdense regions probably happened before z~1, while it is still ongoing at lower z in underdense environments.

We study the evolution of galaxies inside and outside of the group environment since z=1 using a large well defined set of groups and galaxies from the zCOSMOS-bright redshift survey in the COSMOS field. The fraction of galaxies with early-type morphologies increases monotonically with M_B luminosity and stellar mass and with cosmic epoch. It is higher in the groups than elsewhere, especially at later epochs. The emerging environmental effect is superposed on a strong global mass-driven evolution, and at z~0.5 and log(M*/Msol)~10.2, the "effect" of group environment is equivalent to (only) about 0.2 dex in stellar mass or 2 Gyr in time. The stellar mass function of galaxies in groups is enriched in massive galaxies. We directly determine the transformation rates from late to early morphologies, and for transformations involving colour and star formation indicators. The transformation rates are systematically about twice as high in the groups as outside, or up to 3-4 times higher correcting for infall and the appearance of new groups. The rates reach values, for masses around the crossing mass 10^10.5 Msol, as high as (0.3-0.7)/Gyr in the groups, implying transformation timescales of 1.4-3 Gyr, compared with less than 0.2/Gyr, i.e. timescales >5 Gyr, outside of groups. All three transformation rates decrease at higher stellar masses, and must decrease also at the lower masses below 10^10 Msol which we cannot well probe. The rates involving colour and star formation are consistently higher than those for morphology, by a factor of about 50%. Our conclusion is that the transformations which drive the evolution of the overall galaxy population since z~1 must occur at a rate 2-4 times higher in groups than outside of them.

We took advantage of the wealth of information provided by the first ~10000 galaxies of the zCOSMOS-bright survey and its group catalogue to study the complex interplay between group environment and galaxy properties. The classical indicator F_blue (fraction of blue galaxies) proved to be a simple but powerful diagnostic tool. We studied its variation for different luminosity and mass selected galaxy samples. Using rest-frame B-band selected samples, the groups galaxy population exhibits significant blueing as redshift increases, but maintains a lower F_blue with respect both to the global and the isolated galaxy population. However moving to mass selected samples it becomes apparent that such differences are largely due to the biased view imposed by the B-band luminosity selection, being driven by the population of lower mass, bright blue galaxies for which we miss the redder, equally low mass, counterparts. By focusing the analysis on narrow mass bins such that mass segregation becomes negligible we find that only for the lowest mass bin explored (logMass <= 10.6) does a significant residual difference in color remain as a function of environment, while this difference becomes negligible toward higher masses. Our results indicate that red galaxies of logMass >= 10.8 are already in place at z ~ 1 and do not exhibit any strong environmental dependence, possibly originating from so-called 'nature'/internal mechanisms. In contrast, for lower galaxy masses and redshifts lower than z ~ 1, we observe the emergence in groups of a population of 'nurture' red galaxies: slightly deviating from the trend of the downsizing scenario followed by the global galaxy population, and more so with cosmic time. These galaxies exhibit signatures of group-related secular physical mechanisms directly influencing galaxy evolution.

We measured stellar velocity dispersions sigma and dynamical masses of 9 massive (M~10^11 Msun) early-type galaxies (ETG) from the GMASS sample at redshift 1.4<z<2.0. The sigma are based on individual spectra for two galaxies at z~1.4 and on a stacked spectrum for 7 galaxies with 1.6<z<2.0, with 202-h of exposure at the ESO Very Large Telescope. We constructed detailed axisymmetric dynamical models for the objects, based on the Jeans equations, taking the observed surface brightness (from deep HST/ACS observations), PSF and slit effects into account. Our dynamical masses M_Jeans agree within ~30% with virial estimates M_vir=5*Re*sigma^2/G, although the latter tend to be smaller. This suggests that sizes are not underestimated by more than a similar fraction. Our M_Jeans also agrees within a factor <2 with the M_pop previously derived using stellar population models and 11 bands photometry. This confirms that the galaxies are intrinsically massive. The inferred mass-to-light ratios M/L_U in the very age-sensitive rest frame U-band are consistent with passive evolution in the past ~1 Gyr (formation redshift z_f~3). A 'bottom-light' stellar Initial Mass Function (IMF) appears to be required to ensure close agreement between M_Jeans and M_pop at z~2, as it does at z~0. The GMASS ETGs are on average more dense than their local counterpart. However a few percent of local ETGs of similar dynamical masses also have comparable sigma and mass surface density Sigma_50 inside Re.

One of the key unanswered questions in the study of galaxy evolution is what physical processes inside galaxies drive the changes in the SFRs in individual galaxies that, taken together, produce the large decline in the global star-formation rate density (SFRD) to redshifts since z~2. Many studies of the SFR at intermediate redshifts have been made as a function of the integrated stellar mass of galaxies but these did not use information on the internal structural properties of the galaxies. In this paper we present a comparative study of the dependence of SFRs on the average surface mass densities (SigmaM) of galaxies of different morphological types up to z~1 using the zCOSMOS and SDSS surveys. The main findings about the evolution of these relatively massive galaxies are: 1) There is evidence that, for both SDSS ans zCOSMOS galaxies, the mean specific SFR within a given population (either disk-dominated or bulge-dominated) is independent of SigmaM; 2) The observed SSFR - SigmaM step-function relation is due, at all investigated redshifts, to the changing mix of disk-dominated and bulge-dominated galaxies as surface density increases and the strong difference in the average SSFR between disks and bulges. We also find a modest differential evolution in the size-mass relations of disk and spheroid galaxies; 3) The shape of the median SSFR - SigmaM relation is similar, but with median SSFR values that are about 5-6 times higher in zCOSMOS galaxies than for SDSS, across the whole range of SigmaM, and in both spheroid and disk galaxies. This increase matches that of the global SFRD of the Universe as a whole, emphasizing that galaxies of all types are contributing, proportionally, to the global increase in SFRD in the Universe back to these redshifts (abridged).

We present a study of a galaxy overdensity at z=1.6, Cl 0332-2742, uncovered in the Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS). We conservatively estimate the overdensity in redshift space for the spike, containing 42 galaxies at z=1.6, to be 8.3+/-1.5 and find a velocity dispersion of 500 km/s. A map of the surface density of galaxies at z=1.6 in the GMASS field shows that its structure is irregular. The differences in the physical properties of Cl 0332-2742 member and field galaxies agree with the latest hierarchical galaxy formation models: for overdensity members, the star formation rate (SFR), and specific SFR, is approximately 50% lower than for the field galaxies; overdensity galaxies are twice the age, on average, of field galaxies; and there is a higher proportion of both massive (M>10^{10.7} M_sol), and early-type galaxies, inside Cl 0332-2742 than in the field. Among the 42 members, 7 have spectra consistent with being passively evolving, massive galaxies. In a z-J colour-magnitude diagram, the photometric data of these early-type galaxies are in close agreement with a theoretical red sequence of a galaxy cluster at redshift z=1.6, which formed most of its stars in a short burst of star formation at z~3. We conclude that the redshift spike at z=1.6 in the GMASS field represents a sheet-like structure in the cosmic web and the area with the highest surface density within this structure, containing already seven passively evolving galaxies, will evolve into a cluster of galaxies at a later time.

For more than two decades we have known that galaxy morphological segregation is present in the Local Universe. It is important to see how this relation evolves with cosmic time. To investigate how galaxy assembly took place with cosmic time, we explore the evolution of the morphology-density relation up to redshift z~1 using about 10000 galaxies drawn from the zCOSMOS Galaxy Redshift Survey. Taking advantage of accurate HST/ACS morphologies from the COSMOS survey, of the well-characterised zCOSMOS 3D environment, and of a large sample of galaxies with spectroscopic redshift, we want to study here the evolution of the morphology-density relation up to z~1 and its dependence on galaxy luminosity and stellar mass. The multi-wavelength coverage of the field also allows a first study of the galaxy morphological segregation dependence on colour. We further attempt to disentangle between processes that occurred early in the history of the Universe or late in the life of galaxies. The zCOSMOS field benefits of high-resolution imaging in the F814W filter from the Advanced Camera for Survey (ACS). We use standard morphology classifiers, optimised for being robust against band-shifting and surface brightness dimming, and a new, objective, and automated method to convert morphological parameters into early, spiral, and irregular types. We use about 10000 galaxies down to I_AB=22.5 with a spectroscopic sampling rate of 33% to characterise the environment of galaxies up to z~1 from the 100 kpc scales of galaxy groups up to the 100 Mpc scales of the cosmic web. ABRIDGED

We study the dependence of galaxy clustering on luminosity and stellar mass at redshifts z ~ [0.2-1] using the first zCOSMOS 10K sample. We measure the redshift-space correlation functions xi(rp,pi) and its projection wp(rp) for sub-samples covering different luminosity, mass and redshift ranges. We quantify in detail the observational selection biases and we check our covariance and error estimate techniques using ensembles of semi-analytic mock catalogues. We finally compare our measurements to the cosmological model predictions from the mock surveys. At odds with other measurements, we find a weak dependence of galaxy clustering on luminosity in all redshift bins explored. A mild dependence on stellar mass is instead observed. At z~0.7, wp(rp) shows strong excess power on large scales. We interpret this as produced by large-scale structure dominating the survey volume and extending preferentially in direction perpendicular to the line-of-sight. We do not see any significant evolution with redshift of the amplitude of clustering for bright and/or massive galaxies. The clustering measured in the zCOSMOS data at 0.5<z<1 for galaxies with log(M/M_\odot)>=10 is only marginally consistent with predictions from the mock surveys. On scales larger than ~2 h^-1 Mpc, the observed clustering amplitude is compatible only with ~1% of the mocks. Thus, if the power spectrum of matter is LCDM with standard normalization and the bias has no unnatural scale-dependence, this result indicates that COSMOS has picked up a particularly rare, ~2-3 sigma positive fluctuation in a volume of ~10^6 h^-1 Mpc^3. These findings underline the need for larger surveys of the z~1 Universe to appropriately characterize the level of structure at this epoch.

The properties of dust attenuation at rest-frame UV wavelengths are inferred from very high-quality FORS2 spectra of 78 galaxies from the GMASS survey at 1<z<2.5. These objects complement a previously investigated sample of 108 UV-luminous galaxies at similar redshifts, selected from the FDF spectroscopic survey, the K20 survey, and the GDDS. The shape of the UV extinction curve is constrained by a parametric description of the rest-frame UV continuum. The UV bump is further characterised by fitting Lorentzian-like profiles. Spectra exhibit a significant 2175A feature in at least 30% of the cases. If attenuation is dominated by dust ejected from the galaxy main body via superwinds, UV extinction curves in-between those of the SMC and LMC characterise the sample galaxies. The fraction of galaxies with extinction curves differing from the SMC one increases, if more dust resides in the galactic plane or dust attenuation depends on stellar age. On average, the width of the manifested UV bumps is about 60% of the values typical of the LMC and Milky Way. This suggests the presence of dust similar to that found in the LMC2 supershell close to 30Dor. The presence of the carriers of the UV bump at 1<z<2.5 argues for outflows from AGB stars being copious then. Consistent with their higher SFRs, the GMASS galaxies with a manifested UV bump are more luminous at rest-frame 8mum, where the emission is dominated by PAHs (also products of AGB stars). In addition, they exhibit stronger UV absorption features, mostly of interstellar origin, which indicates overall more evolved stellar populations. We conclude that diversification of the small-size dust component has already started in the most evolved star-forming systems at 1<z<2.5.

We use the current sample of ~10,000 zCOSMOS spectra of sources selected with I(AB) < 22.5 to define the density field out to z~1, with much greater resolution in the radial dimension than has been possible with either photometric redshifts or weak lensing. We apply new algorithms that we have developed (ZADE) to incorporate objects not yet observed spectroscopically by modifying their photometric redshift probability distributions using the spectroscopic redshifts of nearby galaxies. This strategy allows us to probe a broader range of galaxy environments and reduce the Poisson noise in the density field. The reconstructed overdensity field of the 10k zCOSMOS galaxies consists of cluster-like patterns surrounded by void-like regions, extending up to z~1. Some of these structures are very large, spanning the ~50 Mpc/h transverse direction of the COSMOS field and extending up to Delta z~0.05 in redshift. We present the three dimensional overdensity maps and compare the reconstructed overdensity field to the independently identified virialised groups of galaxies and clusters detected in the visible and in X-rays. The distribution of the overdense structures is in general well traced by these virialised structures. A comparison of the large scale structures in the zCOSMOS data and in the mock catalogues reveals an excellent agreement between the fractions of the volume enclosed in structures of all sizes above a given overdensity between the data and the mocks in 0.2<z<1.

The impact of environment on AGN activity up to z~1 is assessed by utilizing a mass-selected sample of galaxies from the 10k catalog of the zCOSMOS spectroscopic redshift survey. We identify 147 AGN by their X-ray emission as detected by XMM-Newton from a parent sample of 7234 galaxies. We measure the fraction of galaxies with stellar mass M_*>2.5x10^10 Msun that host an AGN as a function of local overdensity using the 5th, 10th and 20th nearest neighbors that cover a range of physical scales (~1-4 Mpc). Overall, we find that AGNs prefer to reside in environments equivalent to massive galaxies with substantial levels of star formation. Specifically, AGNs with host masses between 0.25-1x10^11 Msun span the full range of environments (i.e., field-to-group) exhibited by galaxies of the same mass and rest-frame color or specific star formation rate. Host galaxies having M_*>10^11 Msun clearly illustrate the association with star formation since they are predominantly bluer than the underlying galaxy population and exhibit a preference for lower density regions analogous to SDSS studies of narrow-line AGN. To probe the environment on smaller physical scales, we determine the fraction of galaxies (M_*>2.5x10^10 Msun) hosting AGNs inside optically-selected groups, and find no significant difference with field galaxies. We interpret our results as evidence that AGN activity requires a sufficient fuel supply; the probability of a massive galaxy to have retained some sufficient amount of gas, as evidence by its ongoing star formation, is higher in underdense regions where disruptive processes (i.e., galaxy harrassment, tidal stripping) are lessened.

We study the spatial clustering of 538 X-ray selected AGN in the 2 deg^2 XMM-COSMOS field that are spectroscopically identified to I_{AB}<23 and span the redshift range z=0.2-3.0. The median redshift and luminosity of the sample are z = 0.98 and L_{0.5-10}=6.3 x 10^{43} erg/s, respectively. A strong clustering signal is detected at ~18sigma level, which is the most significant measurement obtained to date for clustering of X-ray selected AGN. By fitting the projected correlation function w(r_p) with a power law on scales of r_p=0.3-40 Mpc/h, we derive a best fit comoving correlation length of r_0 = 8.6 +- 0.5 Mpc/h and slope of gamma=1.88 +- 0.07 (Poissonian errors; bootstrap errors are about a factor of 2 larger). An excess signal is observed in the range r_p~5-15 Mpc/h, which is due to a large scale structure at z ~ 0.36 containing about 40 AGN. When removing the z ~ 0.36 structure, or computing w(r_p) in a narrower range around the peak of the redshift distribution (e.g. z=0.4-1.6), the correlation length decreases to r_0 ~ 5-6 Mpc/h, which is consistent with that observed for bright optical QSOs at the same redshift. We investigate the clustering properties of obscured and unobscured AGN separately. Within the statistical uncertainties, we do not find evidence that AGN with broad optical lines (BLAGN) cluster differently from AGN without broad optical lines (non-BLAGN). The correlation length measured for XMM-COSMOS AGN at z~1 is similar to that of massive galaxies (stellar mass M_*> 3 x 10^{10} M_sun) at the same redshift. This suggests that AGN at z~1 are preferentially hosted by massive galaxies, as observed both in the local and in the distant (z~2) Universe. (shortened)

We present a study of the host galaxies of AGN selected from the zCOSMOS survey to establish if accretion onto supermassive black holes and star formation are explicitly linked up to z~1. We identify 152 galaxies that harbor AGN, based on XMM observations of 7543 galaxies (i<22.5). Star formation rates (SFRs), including those weighted by stellar mass, are determined using the [OII]3727 emission-line, corrected for an AGN contribution. We find that the majority of AGN hosts have significant levels of star formation with a distribution spanning ~1-100 Msun yr^-1. The close association between AGN activity and star formation is further substantiated by an increase in the AGN fraction with the youthfulness of their stars as indicated by the rest-frame color (U-V) and spectral index Dn(4000); we demonstrate that mass-selection alleviates an artifical peak falling in the transition region when using luminosity-limited samples. We also find that the SFRs of AGN hosts evolve with cosmic time in a manner that closely mirrors the overall galaxy population and explains the low SFRs in AGNs (z<0.3) from the SDSS. We conclude that the conditions most conducive for AGN activity are a massive host galaxy and a large reservoir of gas. Furthermore, a direct correlation between mass accretion rate onto SMBHs and SFR is shown to be weak although the average ratio is constant with redshift, effectively shifting the evidence for a co-evolution scenario in a statistical manner to smaller physical scales. Our findings illustrate an intermittent scenario with an AGN lifetime substantially shorter than that of star formation and underlying complexities regarding fueling over vastly different physical scales yet to be determined [Abridged].

Aims. We investigate the relationships between three main optical galaxy observables (spectral properties, colours, and morphology), exploiting the data set provided by the COSMOS/zCOSMOS survey. The purpose of this paper is to define a simple galaxy classification cube, using a carefully selected sample of around 1000 galaxies. Methods. Using medium resolution spectra of the first 1k zCOSMOS-bright sample, optical photometry from the Subaru/COSMOS observations, and morphological measurements derived from ACS imaging, we analyze the properties of the galaxy population out to z~1. Applying three straightforward classification schemes (spectral, photometric, and morphological), we identify two main galaxy types, which appear to be linked to the bimodality of galaxy population. The three parametric classifications constitute the axes of a "classification cube". Results. A very good agreement exists between the classification from spectral data (quiescent/star-forming galaxies) and that based on colours (red/blue galaxies). The third parameter (morphology) is less well correlated with the first two: in fact a good correlation between the spectral classification and that based on morphological analysis (early-/late-type galaxies) is achieved only after partially complementing the morphological classification with additional colour information. Finally, analyzing the 3D-distribution of all galaxies in the sample, we find that about 85% of the galaxies show a fully concordant classification, being either quiescent, red, bulge-dominated galaxies (~20%) or star-forming, blue, disk-dominated galaxies (~65%). These results imply that the galaxy bimodality is a consistent behaviour both in morphology, colour and dominant stellar population, at least out to z~1.

We investigate the close environment of 203 Spitzer 24 micron-selected sources at 0.6<z<1.0 using zCOSMOS-bright redshifts and spectra of I<22.5 AB mag galaxies, over 1.5 sq. deg. of the COSMOS field. We quantify the degree of passivity of the LIRG and ULIRG environments by analysing the fraction of close neighbours with Dn(4000)>1.4. We find that LIRGs at 0.6<z<0.8 live in more passive environments than those of other optical galaxies that have the same stellar mass distribution. Instead, ULIRGs inhabit more active regions (e.g. LIRGs and ULIRGs at 0.6<z<0.8 have, respectively, (42.0 +/- 4.9)% and (24.5 +/- 5.9)% of neighbours with Dn (4000)>1.4 within 1 Mpc and +/- 500 km/s). The contrast between the activities of the close environments of LIRGs and ULIRGs appears especially enhanced in the COSMOS field density peak at z~0.67, because LIRGs on this peak have a larger fraction of passive neighbours, while ULIRGs have as active close environments as those outside the large-scale structure. The differential environmental activity is related to the differences in the distributions of stellar mass ratios between LIRGs/ULIRGs and their close neighbours, as well as in the general local density fields. At 0.8<z<1.0, instead, we find no differences in the environment densities of ULIRGs and other similarly massive galaxies, in spite of the differential activities. We discuss a possible scenario to explain these findings.

The Galaxy Mass Assembly ultra-deep Spectroscopic Survey samples a part of the CDFS to unprecedented depth. The resulting distribution of 150 z>1.4 redshifts reveals a significant peak at z=1.6, part of a larger overdensity found at this redshift. The 42 spectroscopic members of this structure, called Cl 0332-2742, form an overdensity in redshift of a factor 11+/-3 and have a velocity dispersion of 450 km/s. We derive a total mass for Cl 0332-2742 of ~7x10^14 Msol. The colours of its early-type galaxies are consistent with a theoretical red sequence of galaxies with stars formed at z=3.0. In addition, there are more massive, passive and older, but less star forming galaxies in CL 0332-2742 than in the field. We conclude that this structure is a cluster under assembly at z=1.6.

The aim of this work is to study the evolution of the rest frame color distribution of galaxies with the redshift, in particular in the critical interval 1.4<z<3. We combine ultradeep spectroscopy from the GMASS project (Galaxy Mass Assembly ultradeep Spectroscopic Survey) with GOODS multi-band photometry (from optical to mid-infrared) to study a sample of 1021 galaxies up to m(4.5$\mu$m)=23. We find that the distribution of galaxies in the (U-B) color vs stellar mass plane is bimodal up to at least redshift z=2. We define a mass complete sample of galaxies residing on the red-sequence, selecting objects with log(M/M_{\odot})>10.1, and we study their morphological and spectro-photometric properties. We show that the contribution to this sample of early-type galaxies, defined as galaxies with a spheroidal morphology and no star formation, decreases from 60-70% at z<0.5 down to ~50% at redshift z=2. At z>2 we still find red galaxies in the mass complete sample, even if the bimodality is not seen any more. About 25% of these red galaxies at z>2 are passively evolving, with the bulk of their stars formed at redshift z>`3.

We combine ultradeep optical spectroscopy from the GMASS project ("Galaxy Mass Assembly ultradeep Spectroscopic Survey") with GOODS multi-band photometry and HST imaging to study a sample of passive galaxiesat 1.39<z<1.99 selected at 4.5 microns. A stacked spectrum with an equivalent integration time of ~500 hours was obtained is publicly released. The spectral and photometric SED properties indicate very weak or absent star formation, moderately old stellar ages of ~1 Gyr (for solar metallicity) and stellar masses in the range of 10^{10-11} solar masses, thus implying that the major star formation and assembly processes for these galaxies occurred at z>2. These galaxies have morphologies that are predominantly compact and spheroidal.However, their sizes (R_e <~ 1 kpc) are much smaller than those of spheroids in the present--day Universe. Their stellar mass surface densities are consequently higher by ~1 dex if compared to spheroids at z~0 with the same mass. Their rest-frame B-band surface brightness scales with the effective radius, but the offset with respect to the surface brightness of the local Kormendy relation is too large to be explained by simple passive evolution. At z~1, a larger fraction of passive galaxies follows the z~0 size -- mass relation. Superdense relics with R_e~1 kpc are extremely rare at z~0 with respect to z>1, and absent if R_e<1 kpc. Because of the similar sizes and mass densities, we suggest that the superdense passive galaxies at 1<z<2 are the remnants of the powerful starbursts occurring in submillimeter--selected galaxies at z>2. The results are compared with theoretical models and the main implications discussed in the framework of massive galaxy formation and evolution.

Context: Galaxy metallicities have been measured to redshift z~2 by gas-phase oxygen abundances of the interstellar medium using the R23 and N2 methods. Galaxy stellar metallicities provide crucial data for chemical evolution models but have not been assessed reliably much outside the local Universe. Aims: We determine the iron-abundance, stellar metallicity of star-forming galaxies (SFGs) at redshift z~2, observed as part of the Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS). Methods: We compute the equivalent width of a rest-frame mid-ultraviolet, photospheric absorption-line index, the 1978 index found to vary monotonically with stellar metallicity by Rix and collaborators. We normalise and combine 75 SFG spectra from the GMASS survey to produce a spectrum corresponding to a total integration time 1652.5 hours (and a signal-to-noise ratio ~100 for our 1.5 angstrom binning) of FORS2 spectroscopic observations at the Very Large Telescope. Results: We measure an iron-abundance, stellar metallicity of log (Z/Zsolar) = -0.574+/-0.159 for our spectrum representative of a galaxy of stellar mass 9.4 x 10^9 Msolar assuming a Chabrier IMF. We find that the R04 model SFG spectrum for log (Z/Zsolar) = -0.699 solar metallicity provides the best description of our GMASS coadded spectrum. For similar galaxy stellar mass, our stellar metallicity is ~0.25 dex lower than the oxygen-abundance, gas-phase metallicity quantified by Erb and collaborators for UV-selected star-forming galaxies at z=2. Conclusions: We conclude that we are witnessing the establishment of a light-element overabundance in galaxies as they are being formed at redshift z~2. Our measurements are reminiscent of the alpha-element enhancement seen in low-redshift, galactic bulges and early-type galaxies. (Abridged)

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

We study the enviromental dependence and the morphological composition of the galaxy color-magnitude diagram at z~0.7, using a pilot sub-sample of ~2000 galaxies from the COSMOS surve, with I_AB<24, photometric redshift within 0.61<z<0.85 and galaxy morphologies based on the HST-ACS data. The C-M diagram shows a clear red-sequence dominated by early-type galaxies and also a remarkably well-defined "blue sequence'" described by late-type objects. While the percentage of objects populating the two sequences is a function of environment, with a clear morphology/color-density relation, their normalization and slope are independent of local density. We identify and study a number of objects with "anomalous'' colors, given their morphology. Red late-type galaxies are found to be mostly highly-inclined spiral galaxies, with colors dominated by internal reddening (33% contamination with respect to truly passive spheroidals). Conversely, the population of blue early-type galaxies is composed by objects of moderate luminosity and mass, concurring to only ~5% of the mass in spheroidal galaxies. The majority of them (~70%) occupy a position in the surface-brightness/effective-radius plane not consistent with them being precursors of current epoch elliptical galaxies. In a color-mass diagram, color sequences are even better defined, with red galaxies covering in general a wider range of masses at nearly constant color, and blue galaxies showing a more pronounced dependence of color on mass. While the red sequence is adequately reproduced by models of passive evolution, the blue sequence is better interpreted as a specific star-formation sequence. Its substantial invariance with respect to local density suggests that the overall, "secular'' star formation is driven more by galaxy mass than by environment. [...]

We have identified a large-scale structure at z~0.73 in the COSMOS field, coherently described by the distribution of galaxy photometric redshifts, an ACS weak-lensing convergence map and the distribution of extended X-ray sources in a mosaic of XMM observations. The main peak seen in these maps corresponds to a rich cluster with Tx= 3.51+0.60/-0.46 keV and Lx=(1.56+/-0.04) x 10^{44} erg/s ([0.1-2.4] keV band). We estimate an X-ray mass within $r500$ corresponding to M500~1.6 x 10^{14} Msun and a total lensing mass (extrapolated by fitting a NFW profile) M(NFW)=(6+/-3) x 10^15 Msun. We use an automated morphological classification of all galaxies brighter than I_AB=24 over the structure area to measure the fraction of early-type objects as a function of local projected density Sigma_10, based on photometric redshifts derived from ground-based deep multi-band photometry. We recover a robust morphology-density relation at this redshift, indicating, for comparable local densities, a smaller fraction of early-type galaxies than today. Interestingly, this difference is less strong at the highest densities and becomes more severe in intermediate environments. We also find, however, local "inversions'' of the observed global relation, possibly driven by the large-scale environment. In particular, we find direct correspondence of a large concentration of disk galaxies to (the colder side of) a possible shock region detected in the X-ray temperature map and surface brightness distribution of the dominant cluster. We interpret this as potential evidence of shock-induced star formation in existing galaxy disks, during the ongoing merger between two sub-clusters.

(Abbrev.) This paper prepares a series of papers analysing the Intermediate MAss Galaxy Evolution Sequence (IMAGES) up to z=1. Intermediate mass galaxies (MJ <=-20.3) are selected from the Chandra Deep Field South (CDFS) for which we identify a serious lack of spectroscopically determined redshifts..... We have spectroscopically identified 691 objects including 580 gal., 7 QSOs, and 104 stars. This study provides 531 new redshifts in the CDFS. It confirms the presence of several large scale structures in the CDFS. To test the impact of these structures in the GOODS-South field, we ... compare the evolution of rest-frame U, B, V and K galaxy luminosity densities to that derived from the CFRS. The CDFS field shows a significant excess of luminosity densities in the z=0.5-0.75 range, which increases with the wavelength, reaching up to 0.5 dex at 2.1 um. Stellar mass and specific star formation evolutions might be significantly affected by the presence of the peculiar large scale structures at z= 0.668 and at z= 0.735, that contain a significant excess of evolved, massive galaxies when compared to other fields. This leads to a clear warning to results based on the CDFS/GOODS South fields, especially those related to the evolution of red luminosity densities, i.e. stellar mass density and specific star formation rate. Photometric redshift techniques, when applied to that field, are producing quantities which are apparently less affected by cosmic variance (0.25 dex at 2.1 um), however at the cost of the difficulty in disentangling between evolutionary and cosmic variance effects.

One of the main questions of galaxy formation is how and when galaxies assembled their mass. Observations suggest that most mass assembly and star formation occured at 1 < z < 3. However, in the same redshift range, fully assembled old and massive galaxies co-exist with star-bursting systems still in their assembling stage. In order to explore galaxy evolution in this critical redshift range, we started a project called GMASS (Galaxy Mass Assembly ultra-deep Spectroscopic Survey) based on an 145 hours ESO Large Program and aimed at doing ultradeep spectroscopy with the VLT of a sample of high redshift galaxies slected with IRAC/Spitzer m(4.5mu) < 23, z(phot) > 1.4) from the GOODS-South/HUDF region. The selection at 4.5 micron is crucial to select high redshift galaxies based on the stellar mass or objects affected by strong dust extinction. The spectroscopic deep integration times (15 - 50 hours) are needed to identify unambiguously the redshift and nature of these galaxies, which redshifts are mostly in the "desert" where spectroscopy is traditionally difficult due to the absence of strong emission lines at optical wavelengths. The survey is almost complete and some preliminary results are presented here.

We compare the history of the galaxy mass build-up, as inferred from near-IR observations, and the Star Formation Rate of massive stars in the comoving volume traced by deep extensive far-IR surveys, both possible now with the Spitzer Space Telescope. These two independent and complementary approaches to the history of galaxy formation consistently indicate that a wide interval of cosmic epochs between z~0.7 to z~2 brackets the main evolutionary phases. The rate of the integrated galaxy mass growth indicated by the IR-based comoving SFR appears consistent with the observed decrease of the stellar mass densities with redshift. There are also indications that the evolution with z of the total population depends on galaxy mass, being stronger for moderate-mass, but almost absent up to z=1.4 for high-mass galaxies, thus confirming previous evidence for a "downsizing" effect in galaxy formation. The most massive galaxies appear already mostly in place by z~1. Although a precise matching of this galaxy build-up with the growth of nuclear super-massive black-holes is not possible with the present data (due to difficulties for an accurate census of the obscured AGN phenomenon), some preliminary indications reveal a similar mass/luminosity dependence for AGN evolution as for the hosting galaxies.

We constrain the evolution of the galaxy mass and luminosity functions from the analysis of (public) multi-wavelength data in the Chandra Deep Field South (CDFS) area, obtained from the GOODS and other projects, and including very deep high-resolution imaging by HST/ACS. Our reference catalogue of faint high-redshift galaxies, which we have thoroughly tested for completeness and reliability, comes from a deep (S(3.6micron)>1 microJy) image by IRAC on the Spitzer Observatory. These imaging data in the field are complemented with extensive optical spectroscopy by the ESO VLT/FORS2 and VIMOS spectrographs, while deep K-band VLT/ISAAC imaging is also used to derive further complementary statistical constraints and to assist the source identification and SED analysis. We have selected a highly reliable IRAC 3.6micron sub-sample of 1478 galaxies with S(3.6)>10microJy, 47% of which have spectroscopic redshift, while for the remaining objects both COMBO-17 and Hyperz are used to estimate the photometric redshift. This very extensive dataset is exploited to assess evolutionary effects in the galaxy luminosity and stellar mass functions, while luminosity/density evolution is further constrained with the number counts and redshift distributions. The deep ACS imaging allows us to differentiate these evolutionary paths by morphological type, which our simulations show to be reliable at least up to z=1.5 for the two main early- (E/S0) and late-type (Sp/Irr) classes. These data, as well as our direct estimate of the stellar mass function above M=10^(10)M_sun for the spheroidal subclass, consistently evidence a progressive dearth of such objects to occur starting at z=0.7, paralleled by an increase in luminosity. (abridged)

We present results of a morphological analysis of a small subset of the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) galaxy population. The analysis is based on public ACS data taken inside the SWIRE N1 field, which are the deepest optical high-resolution imaging available within the SWIRE fields as of today. Our reference sample includes 156 galaxies detected by both ACS and SWIRE. Among the various galaxy morphologies, we disentangle two main classes, spheroids (or bulge-dominated galaxies) and disk-dominated ones, for which we compute the number counts as a function of flux. We then limit our sample to objects with IRAC fluxes brighter than 10 microJy, estimated ~90% completeness limit of the SWIRE catalogues, and compare the observed counts to model predictions. We find that the observed counts of the spheroidal population agree with the expectations of a hierarchical model while a monolithic scenario predicts steeper counts. Both scenaria, however, under-predict the number of late-type galaxies. These observations show that the large majority (close to 80 per cent) of the 3.6 and 4.5 micron galaxy population, even at these moderately faint fluxes, is dominated by spiral and irregular galaxies or mergers.

We have performed VLT spectroscopy of an almost complete sample of 18 early-type galaxies with 0.88 < z < 1.3 plus two at z=0.67, selected from the K20 survey, and derived the velocity dispersion for 15+2 of them. By combining these data with HST and VLT images, we study the Fundamental Plane (FP), the Faber-Jackson and the Kormendy relations at z~1, and compare them with the local ones. The FP at z~1 keeps a remarkably small scatter, and shows both an offset and a rotation, which we interpret in terms of evolution of the mass-to-light ratio, and possibly of the size. We show evidence that the evolution rate depends on galaxy mass, being faster for less massive galaxies. We discuss the possible factors driving the evolution of spheroids and compare our results with the predictions of the hierachical models of galaxy formation.

We present a detailed analysis of the evolution of the rest-frame B-band morphology of K-selected galaxies with 0 < z < 2.5. This work is based on the K20 spectroscopic sample (Ks<20) located within the Chandra Deep Field South area, coupled with the public deep GOODS HST+ACSmulti-band optical imaging available in that field. Thanks to the spectroscopic completeness of this catalog reaching 94%, we can compare the morphological and spectroscopic properties of galaxies with unprecedented detail. Our morphological analysis includes visual inspection and automatic procedures using both parametric (e.g. the Sersic indices treated by the GALFIT and GASPHOT packages) and non-parametric (the Concentration, Asymmetry and clumpineSs, CAS) methods. Our analysis classifies the K20 galaxies as: 60/300 (20%, class 1) normal ellipticals/S0; 14/300 (4%, class 2) perturbed or peculiar ellipticals; 80/300 (27%, class 3) normal spirals; 48/300 (16%, class 4) perturbed or actively star-forming spirals; 98/300 (33%, class 5) irregulars. The evolution of the merging fraction is constrained up to z~2, by carefully accounting effects of morphological k-correction: both asymmetry criterion and pair statistic show an increasing merging fraction as a function of the redshift. We finally analyse the redshift-dependence of the effective radii for early- and late-type galaxies and find some mild evidence for a decrease with z of the early-type galaxy sizes, while disks and irregulars remain constant. Altogether, this analysis of the K20 sample indicates the large predominance of spirals and irregulars at 0.5<z<1.5 in K-band selected samples at even moderate depths.

More than half of all stars in the local Universe are found in massive spheroidal galaxies, which are characterized by old stellar populations with little or no current star formation. In present models, such galaxies appear rather late as the culmination of a hierarchical merging process, in which larger galaxies are assembled through mergers of smaller precursor galaxies. But observations have not yet established how, or even when, the massive spheroidals formed, nor if their seemingly sudden appearance when the Universe was about half its present age (at redshift z \approx 1) results from a real evolutionary effect (such as a peak of mergers) or from the observational difficulty of identifying them at earlier epochs. Here we report the spectroscopic and morphological identification of four old, fully assembled, massive (>10^{11} solar masses) spheroidal galaxies at 1.6<z<1.9, the most distant such objects currently known. The existence of such systems when the Universe was only one-quarter of its present age, shows that the build-up of massive early-type galaxies was much faster in the early Universe than has been expected from theoretical simulations.

Deep VLT optical spectroscopy, HST+ACS (GOODS) imaging and VLA observations are used to unveil the nature of a complete sample of 47 EROs with R-Ks>5 and Ks<20. The spectroscopic redshift completeness is 62%. Morphological classification was derived for each ERO through visual inspection and surface brightness profile fitting. Three main ERO morphological types are found: E/S0 galaxies (30-37%), spiral-like (24-46%) and irregular systems (17-39%). The only ERO detected in the radio is likely to host an obscured AGN. The average radio luminosity of the star-forming EROs undetected in the radio implies star formation rates of the order of 33 Msun yr^-1. The colors, redshifts and masses of the E/S0 galaxy subsample imply a minimum formation redshift z_f \sim 2. With this z_f there is enough time to have old and massive stellar spheroids already assembled at z \sim 1. We verify that the R-Ks vs. J-Ks color diagram is efficient in segregating old and dusty-star-forming EROs.