We present a new prospective analysis of deep multi-band imaging with the James Webb Space Telescope (JWST). In this work, we investigate the recovery of high-redshift $5<z<12$ galaxies through extensive image simulations of accepted JWST programs such as CEERS in the EGS field and HUDF GTO. We introduce complete samples of $\sim300,000$ galaxies with stellar masses $\log(M_*/M_\odot)>6$ and redshifts $0<z<15$, as well as galactic stars, into realistic mock NIRCam, MIRI and HST images to properly describe the impact of source blending. We extract the photometry of the detected sources as in real images and estimate the physical properties of galaxies through spectral energy distribution fitting. We find that the photometric redshifts are primarily limited by the availability of blue-band and near-infrared medium-band imaging. The stellar masses and star-formation rates are recovered within $0.25$ and $0.3$ dex respectively, for galaxies with accurate photometric redshifts. Brown dwarfs contaminating the $z>5$ galaxy samples can be reduced to $<0.01$ arcmin$^{-2}$ with a limited impact on galaxy completeness. We investigate multiple high-redshift galaxy selection techniques and find the best compromise between completeness and purity at $5<z<10$ using the full redshift posterior probability distributions. In the EGS field, the galaxy completeness remains higher than $50\%$ for $m_\text{UV}<27.5$ sources at all redshifts, and the purity is maintained above $80$ and $60\%$ at $z\leq7$ and $10$ respectively. The faint-end slope of the galaxy UV luminosity function is recovered with a precision of $0.1-0.25$, and the cosmic star-formation rate density within $0.1$ dex. We argue in favor of additional observing programs covering larger areas to better constrain the bright end.

We discovered a strongly lensed (\mu >40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.

We present the results from the search for Lyman Alpha emitters (LAEs) in the proximity of 11 CIV absorption systems at $z > 4.7$ in the spectrum of the QSO J1030+0524, using data from MUSE. We have found multiple LAE candidates close to four CIV systems at $z_{CIV} = 4.94$-5.74 with log$_{10}$(N CIV[cm$^{-2}$])$ > 13.5.$ At z = 5-6, CIV systems with W0(CIV) $> 0.2$ {\AA} seem more likely to have galaxies with Ly${\alpha}$ emission within $\rho < 200$ proper kpc (4/5 cases), than the CIV systems with W0(CIV)$ < 0.2$ {\AA} (0/6 cases). The impact parameter of LAE-CIV systems with equivalent widths W0(CIV)$ > 0.5 $ {\AA} is in the range 11 $ < \rho <$ 200 proper kpc (pkpc). Furthermore, all candidates are in the luminosity range 0.18-1.15 L*$(Ly{ \alpha }, z = 5.7)$ indicating that the environment of CIV systems within 200 pkpc is populated by the faint end of the Ly${\alpha}$ luminosity function. We report a 0.28 L*$(Ly{\alpha})$ galaxy at a separation of $\rho = 10$ pkpc from a strong CIV absorption (log10(NCIV[cm$^{-2}$]) = 14.52) at $z_{CIV} = 5.72419$. The prevalence of sub-L*(Ly$\alpha$) galaxies in the proximity of z $> 4.9$ CIV systems suggest that the absorbing material is rather young, likely ejected in the recent past of the galaxies at these redshifts. The connection between faint LAEs and high-ionization absorption systems reported in this work, is potentially a consequence of the role of low mass galaxies in the early evolution of the circum-galactic and intergalactic media.

We carry out a detailed strong lensing analysis of a sub-sample of eight galaxy clusters of the CLASH survey, in the redshift range of $ z_{\rm cluster} = [0.23-0.59]$, using extensive spectroscopic information, primarily MUSE archival data complemented with CLASH-VLT redshift measurements. Different models are tested in each cluster depending on the complexity of its mass distribution and on the number of detected multiple images. Four clusters show more than five spectroscopically confirmed multiple image families. In this sample, we do not make use of families that are only photometrically identified, in order to reduce model degeneracies and systematics due to the potential misidentifications of some multiple images. We present spectroscopic confirmation of 27 multiply lensed sources, with no previous spectroscopic measurements, spanning over the redshift range of $z_{\rm src}=[0.7-6.1]$. Moreover, we confirm an average of $48$ galaxy members in the core of each cluster, thanks to the high efficiency and large field of view of MUSE. Despite having different properties (i.e., number of mass components, total mass, redshift, etc), the projected total mass and mass density profiles of all clusters have very similar shapes, when rescaled by independent measurements of $M_{200c}$ and $R_{200c}$. Specifically, we measure the mean value of the projected total mass of our cluster sample within 10 (20)% of $R_{200c}$ to be 0.13 (0.32) of $M_{200c}$, with a remarkably small scatter of 5 (6)%. Furthermore, the large number of high-z sources and the precise magnification maps derived in this work for four clusters add up to the sample of high-quality gravitational telescopes to be used to study the faint and distant Universe. The strong lensing models and the full redshift catalogues from MUSE are publicly available.

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

We search for bright Ly$\rm \alpha$ emitters among Spitzer SMUVS galaxies at z > 2.9 with homogeneous MUSE data. Although it only covers a small region of COSMOS, MUSE has the unique advantage of providing spectral information over the entire field, without the need of target pre-selection. This gives an unbiased detection of all the brightest Ly$\rm \alpha$ emitters among SMUVS sources, which by design are stellar-mass selected galaxies. Within the studied area, ~14% of the SMUVS galaxies at z > 2.9 have Ly$\rm \alpha$ fluxes F$\rm _\lambda$ > 7 x 10$^{-18}$ erg s$^{-1}$ cm$^{-2}$. These Ly$\rm \alpha$ emitters are characterized by three types of emission, 47% show a single line profile, 19% present a double peak or a blue bump and 31% show a red tail. One object (3%) shows both a blue bump and a red tail. We also investigate the spectral energy distribution (SED) properties of the SMUVS MUSE-detected galaxies and MUSE non-detections. After stellar-mass matching both populations, we find that MUSE detected galaxies have generally lower extinction than SMUVS-only objects, while there is no clear intrinsic difference in the mass and age distributions. For the MUSE-detected SMUVS galaxies, we compare the instantaneous SFR lower limit given by Ly$\rm \alpha$ flux with its past average derived from SED fitting, and find evidence for rejuvenation in some of our oldest objects. We also study the spectra of those Ly$\rm \alpha$ emitters which are not detected in SMUVS in the same field. We find different distributions of the emission line profiles, which could be ascribed to the fainter Ly$\rm \alpha$ luminosities of the MUSE-only sources and an intrinsically different mass distribution. Finally, we search for the presence of galaxy associations. MUSE's integral coverage is 20 times more likely to find associations than all other existing spectral data in COSMOS, biased by target pre-selection.

We present 16 new ultrabright $H_{AB}\lesssim25$ galaxy candidates at z~8 identified over the COSMOS/UltraVISTA field. The new search takes advantage of the deepest-available ground-based optical and near-infrared observations, including the DR3 release of UltraVISTA and full-depth Spitzer/IRAC observations from the SMUVS and SPLASH programs. Candidates are selected using Lyman-break criteria, combined with strict optical non-detection and SED-fitting criteria, minimizing contamination by low-redshift galaxies and low-mass stars. HST/WFC3 coverage from the DASH program reveals that one source evident in our ground-based near-IR data has significant substructure and may actually correspond to 3 separate z~8 objects, resulting in a sample of 18 galaxies, 10 of which seem to be fairly robust (with a >97% probability of being at z>7). The UV-continuum slope $\beta$ for the bright z~8 sample is $\beta=-2.2\pm0.6$, bluer but still consistent with that of similarly bright galaxies at z~6 ($\beta=-1.55\pm0.17$) and z~7 ($\beta=-1.75\pm0.18$). Their typical stellar masses are 10$^{9.1^{+0.5}_{-0.4}}M_{\odot}$, with the SFRs of $32^{+44}_{-32}M_{\odot}$/year, specific SFR of $4^{+8}_{-4}$ Gyr$^{-1}$, stellar ages of $\sim22^{+69}_{-22}$\,Myr, and low dust content A$_V=0.15^{+0.30}_{-0.15}$ mag. Using this sample we constrain the bright end of the z~8 UV luminosity function (LF). When combined with recent empty field LF estimates at z~8-9, the resulting z~8 LF can be equally well represented by either a Schechter or a double power-law (DPL) form. Assuming a Schechter parameterization, the best-fit characteristic magnitude is $M^*= -20.95^{+0.30}_{-0.35}$ mag with a very steep faint end slope $\alpha=-2.15^{+0.20}_{-0.19}$. These new candidates include amongst the brightest yet found at these redshifts, 0.5-1.0 mag brighter than found over CANDELS, providing excellent targets for follow-up studies.

We present a VLT/X-Shooter spectroscopy of the Lyman continuum (LyC) emitting galaxy 'Ion2' at z=3.2121 and compare it to that of the recently discovered strongly lensed LyC-emitter at z=2.37, known as the 'Sunburst' arc. Three main results emerge from the X-Shooter spectrum: (a) the Lya has three distinct peaks with the central one at the systemic redshift, indicating a ionised tunnel through which both Lya and LyC radiation escape; (b) the large O32 oxygen index ([OIII]4959-5007 / [OII]3727-3729) of 9.18(-1.32/+1.82) is compatible to those measured in local (z~0.4) LyC leakers; (c) there are narrow nebular high-ionisation metal lines with \sigma_v < 20 km/s, which confirms the presence of young hot, massive stars. The HeII1640 appears broad, consistent with a young stellar component including Wolf-Rayet stars. Similarly, the Sunburst LyC-emitter shows a triple-peaked Lya profile and from VLT/MUSE spectroscopy the presence of spectral features arising from young hot and massive stars. The strong lensing magnification, (\mu > 20), suggests that this exceptional object is a gravitationally-bound star cluster observed at a cosmological distance, with a stellar mass M <~ 10^7 Msun and an effective radius smaller than 20 pc. Intriguingly, sources like Sunburst but without lensing magnification might appear as Ion2-like galaxies, in which unresolved massive star clusters dominate the ultraviolet emission. This work supports the idea that dense young star clusters can contribute to the ionisation of the IGM through holes created by stellar feedback.

The James Webb Space Telescope will provide deep imaging and spectroscopy for sources at redshifts above 6, covering the Epoch of Reionization (EoR, 6 < z < 10). The Mid-IR instrument (MIRI) integral field spectrograph (MRS) will be the only instrument on board JWST able to observe the brightest optical emission lines H$\alpha$ and [OIII]0.5007$\mu$m at redshifts above 7 and 9, respectively. This paper presents a study of the H$\alpha$ fluxes predicted by FIRSTLIGHT cosmological simulations for galaxies at redshifts of 6.5 to 10.5, and its detectability with MIRI. Deep (40 ks) spectroscopic integrations with MRS will be able to detect (S/N > 5) EoR sources at redshifts above 7 with intrinsic star formation rates of more than 2 M$_{\odot}$ yr$^{-1}$, and stellar masses above 4-9 $\times$ 10$^7$ M$_{\odot}$. In addition, the paper presents realistic MRS simulated observations of the expected (rest-frame) optical and near-infrared spectra for some spectroscopically confirmed EoR sources detected by ALMA as [OIII]88$\mu$m emitters. The MRS simulated spectra cover a wide range of low metallicities from about 0.2 to 0.02Z$_{\odot}$, and different [OIII]88$\mu$m/[OIII]0.5007$\mu$m line ratios. The simulated 10ks MRS spectra show S/N in the range of 5 to 90 for H$\beta$, [OIII]0.4959,0.5007$\mu$m, H$\alpha$ and HeI1.083$\mu$m emission lines of MACS1149-JD1 at z = 9.11, independent of metallicity. In addition, deep 40 ks simulated spectra of the luminous merger candidate B14-65666 at z=7.15 shows the MRS capabilities of detecting, or putting strong upper limits, on the [NII]0.6584$\mu$m, [SII]0.6717,0.6731$\mu$m, and [SIII]0.9069,0.9532$\mu$m emission lines. In summary, MRS will enable the detailed study of key physical properties like internal extinction, instantaneous star formation, hardness of the ionising continuum, and metallicity, in bright (intrinsic or lensed) EoR sources.

Stellar masses in future James Webb Space Telescope (JWST) deep blank-fields will be mainly derived fitting the spectral energy distribution with theoretical galaxy templates. We investigate the uncertainties and biases of the stellar masses derived by using the LePhare code for spectral energy distribution fitting and the Yggdrasil theoretical templates. We consider a sample of mock galaxies at z$=$7-10 with mock JWST observations with S/N$_{F150W}\geqslant$ 10. Our goal is to provide a list of statistical stellar mass corrections to include on the stellar mass derivation for different output galaxy properties and JWST filter combinations to correct for template degeneracies. Median statistical stellar mass corrections vary from -0.83 dex to 0.87 dex, while 25$\%$ (75$\%$) quartiles range from -0.83 dex (-0.67 dex) to 0.51 dex (0.88 dex), depending on filter combinations and galaxy models. The most challenging cases are galaxies with nebular emission lines, especially the ones that are wrongly identified as galaxies without, relative dust-free galaxies and galaxies with small metallicities (i.e. Z$=1/50Z_{\odot}$). The stellar mass estimation of galaxies correctly identified without emission lines is generally fine, except at z$=$10 when considering only the 8 NIRCam bands, which make the MIRI bands very valuable. We have tested our stellar mass corrections using the public JAGUAR galaxy catalogue, deriving that the average discrepancy in the recovered stellar mass distribution decreases by 20-50$\%$ at z$>$7 after the correction. We found that without the stellar-mass corrections the number of low-mass galaxies (M$^{*}<10^{7}M_{\odot}$) is overestimated, which can potentially lead to systematic errors in the calculation of the galaxy stellar mass function faint-end slope at high z.

We select a sample of 90 obscured (type2) AGN with 1.45<z<3.05 from the zCOSMOS-deep galaxy sample by 5 sigma-detection of the high-ionization CIV {\lambda}1549 narrow emission line. The presence of this feature in a galaxy spectrum is often associated with nuclear activity, and the selection effectiveness has been also confirmed by ultraviolet (UV) emission line ratio diagnostic diagrams. Applying the same selection technique, a sample of 102 unobscured (type 1) AGN was collected. Taking advantage of the large amount of multi-band data available in the COSMOS field, we investigate the properties of the CIV-selected type 2 AGN, focusing on their host galaxies, X-ray emission and UV emission lines. Finally, we investigate the physical properties of the ionized gas in the Narrow Line Region (NLR) of this type 2 AGN sample, combining the analysis of strong UV emission lines with predictions from photo-ionization models. We find that, in order to successfully reproduce the relative intensity of UV emission lines of the selected high-z type 2 AGN, two new ingredients in the photo-ionization models are fundamental,i.e. small inner radii of the NLR (~90pc for LAGN = 10^45erg/s) and the internal dissipative micro-turbulence of the gas emitting clouds (with vmicr~100km/s). With these modified models, we compute the gas-phase metallicity of the NLR, and our measurements indicate a statistically significant evolution of the metal content with redshift. Finally, we do not observe, in our CIV-selected type 2 AGN sample, a strong relationship between the NLR gas metallicity and the stellar mass of the host galaxy.

In this work, we use measurements of galaxy stellar mass and two-point angular correlation functions to constrain the stellar-to-halo mass ratios (SHMRs) of passive and \np\ galaxies at $z\sim2-3$, as identified in the \emph{Spitzer} Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS). We adopt a sophisticated halo modeling approach to statistically divide our two populations into central and satellite galaxies. For central galaxies, we find that the normalization of the SHMR is greater for our passive population. Through the modeling of $\Lambda$ cold dark matter halo mass accretion histories, we show that this can only arise if the conversion of baryons into stars was more efficient at higher redshifts and additionally that passive galaxies can be plausibly explained as residing in halos with the highest formation redshifts (i.e., those with the lowest accretion rates) at a given halo mass. At a fixed stellar mass, satellite galaxies occupy host halos with a greater mass than central galaxies, and we find further that the fraction of passive galaxies that are satellites is higher than for the combined population. This, and our derived satellite quenching timescales, combined with earlier estimates from the literature, support dynamical/environmental mechanisms as the dominant process for satellite quenching at $z\lesssim3$.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

We report detections of two 1.2 mm continuum sources ($S_\mathrm{1.2mm}$ ~ 0.6 mJy) without any counterparts in the deep $H$- and/or $K$-band image (i.e., $K$-band magnitude $\gtrsim$ 26 mag). These near-infrared-dark faint millimeter sources are uncovered by ASAGAO, a deep and wide-field ($\simeq$ 26 arcmin$^2$) Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm survey. One has a red IRAC (3.6 and 4.5 $\mu$m) counterpart, and the other has been independently detected at 850 and 870 $\mu$m using SCUBA2 and ALMA Band 7, respectively. Their optical to radio spectral energy distributions indicate that they can lie at $z \gtrsim$ 3-5 and can be in the early phase of massive galaxy formation. Their contribution to the cosmic star formation rate density is estimated to be ~ 1 $\times$ 10$^{-3}$ $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ if they lie somewhere in the redshift range of $z$ ~ 3-5. This value can be consistent with, or greater than that of bright submillimeter galaxies ($S_\mathrm{870\mu m}>$ 4.2 mJy) at $z$ ~ 3-5. We also uncover 3 more candidates near-infrared-dark faint ALMA sources without any counterparts ($S_\mathrm{1.2mm}$ ~ 0.45-0.86 mJy). These results show that an unbiased ALMA survey can reveal the dust-obscured star formation activities, which were missed in previous deep optical/near-infrared surveys.

We present the survey design, data reduction, construction of images, and source catalog of the Atacama Large Millimeter/submillimeter Array (ALMA) twenty-six arcmin^2 survey of GOODS-S at one-millimeter (ASAGAO). ASAGAO is a deep (1sigma ~ 61 uJy/beam for a 250 klambda-tapered map with a synthesized beam size of 0.51" x 0.45") and wide area (26 arcmin^2) survey on a contiguous field at 1.2 mm. By combining with ALMA archival data in the GOODS-South field, we obtained a deeper map in the same region (1sigma ~ 30 uJy/beam for a deep region with a 250 klambda-taper, and a synthesized beam size of 0.59" x 0.53"), providing the largest sample of sources (25 sources at >=5.0sigma, 45 sources at >=4.5sigma) among ALMA blank-field surveys to date. The number counts shows that 52(+11 -8)% of the extragalactic background light at 1.2 mm is resolved into discrete sources at S1.2m > 135 uJy. We create infrared (IR) luminosity functions (LFs) in the redshift range of z = 1-3 from the ASAGAO sources with KS-band counterparts, and constrain the faintest luminosity of the LF at 2.0 < z < 3.0. The LFs are consistent with previous results based on other ALMA and SCUBA-2 observations, which suggest a positive luminosity evolution and negative density evolution with increasing redshift. We find that obscured star-formation of sources with IR luminosities of log(L(IR)/Lsun)} ~> 11.8 account for ~~60%-90% of the z ~ 2 cosmic star-formation rate density.

The Spitzer Matching Survey of the UltraVISTA Ultra-deep Stripes (SMUVS) has obtained the largest ultra-deep Spitzer maps to date in a single field of the sky. We considered the sample of about 66,000 SMUVS sources at $z=2-6$ to investigate the evolution of dusty and non-dusty galaxies with stellar mass through the analysis of the galaxy stellar mass function (GSMF). We further divide our non-dusty galaxy sample with rest-frame optical colours to isolate red quiescent (`passive') galaxies. At each redshift, we identify a characteristic stellar mass in the GSMF above which dusty galaxies dominate, or are at least as important as non-dusty galaxies. Below that stellar mass, non-dusty galaxies comprise about 80% of all sources, at all redshifts except at $z=4-5$. The percentage of dusty galaxies at $z=4-5$ is unusually high: 30-40% for $M_{*}=10^9 - 10^{10.5} \, \rm M_\odot$ and $>80\%$ at $M_*>10^{11} \, \rm M_\odot$, which indicates that dust obscuration is of major importance in this cosmic period. The overall percentage of massive ($\log_{10} (M_*/M_\odot)>10.6$) galaxies that are quiescent increases with decreasing redshift, reaching $>30\%$ at $z\sim2$. Instead, the quiescent percentage among intermediate-mass galaxies (with $\log_{10} (M_*/M_\odot)=9.7-10.6$) stays roughly constant at a $\sim 10\%$ level. Our results indicate that massive and intermediate-mass galaxies clearly have different evolutionary paths in the young Universe, and are consistent with the scenario of galaxy downsizing.

This paper describes new deep 3.6 and 4.5 micron imaging of three UltraVISTA near-infrared survey stripes within the COSMOS field. The observations were carried out with Spitzer's Infrared Array Camera (IRAC) for the Spitzer Matching Survey of the Ultra-VISTA Deep Stripes (SMUVS). In this work we present our data reduction techniques, and document the resulting mosaics, coverage maps, and catalogs in both IRAC passbands for the three easternmost UltraVISTA survey stripes, covering a combined area of about 0.66 square degrees, of which 0.45 square degrees have at least 20 hr integration time. SMUVS reaches point-source sensitivities of about 25.0 AB mag at both 3.6 and 4.5 microns with a significance of 4-sigma accounting for both survey sensitivity and source confusion. To this limit the SMUVS catalogs contain a total of about 350,000 sources, each of which is detected significantly in at least one IRAC band. Because of its uniform and high sensitivity, relatively large area coverage, and the wide array of ancillary data available in COSMOS, the SMUVS survey will be useful for a large number of cosmological investigations. We will make all images and catalogues described herein publicly available via the Spitzer Science Center.

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

We report on the serendipitous discovery of a z=4.0, M1500=-22.20 star-forming galaxy (Ion3) showing copious Lyman continuum (LyC) leakage (~60% escaping), a remarkable multiple peaked Lya emission, and significant Lya radiation directly emerging at the resonance frequency. This is the highest redshift confirmed LyC emitter in which the ionising and Lya radiation possibly share a common ionised cavity (with N_HI<10^17.2 cm^-2). Ion3 is spatially resolved, it shows clear stellar winds signatures like the P-Cygni NV1240 profile, and has blue ultraviolet continuum (\beta = -2.5 +/- 0.25, F_\lambda~ \lambda^\beta) with weak low-ionisation interstellar metal lines. Deep VLT/HAWKI Ks and Spitzer/IRAC 3.6um and 4.5um imaging show a clear photometric signature of the Halpha line with equivalent width of 1000A rest-frame emerging over a flat continuum (Ks-4.5um ~ 0). From the SED fitting we derive a stellar mass of 1.5x10^9 Msun, SFR of 140 Msun/yr and age of ~10 Myr, with a low dust extinction, E(B-V)< 0.1, placing the source in the starburst region of the SFR-M^* plane. Ion3 shows similar properties of another LyC emitter previously discovered (z=3.21, Ion2, Vanzella et al. 2016). Ion3 (and Ion2) represents ideal high-redshift reference cases to guide the search for reionising sources at z>6.5 with JWST.

We use Multi Unit Spectroscopic Explorer (MUSE) observations of the galaxy cluster MACS J1149.5+2223 to explore the kinematics of the grand-design spiral galaxy Sp1149 hosting the SN Refsdal. Sp1149 lies at $z\simeq1.49$, has a stellar mass $M_*\simeq5\times10^9 \, \mathrm{M_\odot}$, a star-formation rate $\mathrm{SFR} \simeq1-6 \, \mathrm{M_\odot/yr}$ and represents a likely progenitor of a Milky-Way-like galaxy. All the four multiple images of Sp1149 in our data show strong OII-line emissions pointing to a clear rotation pattern. We take advantage of the gravitational lensing magnification effect ($\simeq 4 \times$) on the OII emission of the least distorted image to fit 3D kinematic models to the MUSE data-cube and derive the rotation curve and the velocity dispersion profile of Sp1149. We find that the rotation curve steeply rises, peaks at $R\simeq1$ kpc and then (initially) declines and flattens to an average $V_\mathrm{flat} = 128^{+29}_{-19}$ km/s. The shape of the rotation curve is well determined but the actual value of $V_\mathrm{flat}$ is quite uncertain because of the nearly face-on configuration of the galaxy. The intrinsic velocity dispersion due to gas turbulence is almost constant across the entire disc with an average of $27\pm5$ km/s. This value is consistent with $z=0$ measurements in the ionized gas component and a factor 2-4 lower than other estimates in different galaxies at similar redshifts. The average stellar-to-total mass fraction is of the order of one fifth. Our kinematic analysis returns the picture of a regular star-forming, mildly turbulent, rotation-dominated ($V / \sigma\simeq5$) spiral galaxy in a 4 Gyr old Universe.

The \emph{Spitzer} Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS) provides unparalleled depth at $3.6$ and $4.5$~$\mu$m over $\sim0.66$~deg$^2$ of the COSMOS field, allowing precise photometric determinations of redshift and stellar mass. From this unique dataset we can connect galaxy samples, selected by stellar mass, to their host dark matter halos for $1.5<z<5.0$, filling in a large hitherto unexplored region of the parameter space. To interpret the observed galaxy clustering we utilize a phenomenological halo model, combined with a novel method to account for uncertainties arising from the use of photometric redshifts. We find that the satellite fraction decreases with increasing redshift and that the clustering amplitude (e.g., comoving correlation length / large-scale bias) displays monotonic trends with redshift and stellar mass. Applying $\Lambda$CDM halo mass accretion histories and cumulative abundance arguments for the evolution of stellar mass content we propose pathways for the coevolution of dark matter and stellar mass assembly. Additionally, we are able to estimate that the halo mass at which the ratio of stellar to halo mass is maximized is $10^{12.5_{-0.08}^{+0.10}}$~M$_{\odot}$ at $z\sim2.5$. This peak halo mass is here inferred for the first time from stellar mass-selected clustering measurements at $z\gtrsim2$, and implies mild evolution of this quantity for $z\lesssim3$, consistent with constraints from abundance-matching techniques.

We investigate the X-ray active galactic nucleus (AGN) properties of millimeter galaxies in the Great Observatories Origins Deep Survey South (GOODS-S) field detected with the Atacama Large Millimeter/submillimeter Array (ALMA), by utilizing the Chandra 7-Ms data, the deepest X-ray survey to date. Our millimeter galaxy sample comes from the ASAGAO survey covering 26 arcmin$^2$ (12 sources at a 1.2-mm flux-density limit of $\approx$0.6 mJy), supplemented by the deeper but narrower 1.3-mm survey of a part of the ASAGAO field by Dunlop et al.\ (2017). Fourteen out of the total 25 millimeter galaxies have Chandra counterparts. The observed AGN fractions at $z=1.5-3$ is found to be 90$^{+8}_{-19}$\% and $57^{+23}_{-25}$\% for the ultra/luminous infrared galaxies with $\log L_{\rm IR}/L_{\odot} = 12-12.8$ and $\log L_{\rm IR}/L_{\odot} = 11.5-12$, respectively. The majority ($\sim$2/3) of the ALMA and/or Herschel detected X-ray AGNs at $z=1.5-3$ appear to be star-formation dominant populations, having $L_{\rm X}$/ $L_{\rm IR}$ ratios smaller than the "simultaneous evolution" value expected from the local black-hole mass to stellar mass ($M_{\rm BH}$-$M_*$) relation. On the basis of the $L_{\rm X}$ and stellar mass relation, we infer that a large fraction of star-forming galaxies at $z=1.5-3$ have black hole masses smaller than those expected from the local $M_{\rm BH}$-$M_*$ relation. This is opposite to previous reports on luminous AGN at same redshifts detected in wider and shallower surveys, which are subject to selection biases against lower luminosity AGN. Our results are consistent with an evolutionary scenario that star formation occurs first, and an AGN-dominant phase follows later, in objects finally evolving into galaxies with classical bulges.

We report the study of far-IR sizes of submillimeter galaxies (SMGs) in relation to their dust-obscured star formation rate (SFR) and active galactic nuclei (AGN) presence, determined using mid-IR photometry. We determined the millimeter-wave ($\lambda_{\rm obs}=1100 \mu$m) sizes of 69 ALMA-identified SMGs, selected with $\geq10$$\sigma$ confidence on ALMA images ($F_{\rm 1100 \mu m}=1.7$--7.4 mJy). We found that all the SMGs are located above an avoidance region in the millimeter size-flux plane, as expected by the Eddington limit for star formation. In order to understand what drives the different millimeter-wave sizes in SMGs, we investigated the relation between millimeter-wave size and AGN fraction for 25 of our SMGs at $z=1$--3. We found that the SMGs for which the mid-IR emission is dominated by star formation or AGN have extended millimeter-sizes, with respective median $R_{\rm c,e} = 1.6^{+0.34}_{-0.21}$ and 1.5$^{+0.93}_{-0.24}$ kpc. Instead, the SMGs for which the mid-IR emission corresponds to star-forming/AGN composites have more compact millimeter-wave sizes, with median $R_{\rm c,e}=1.0^{+0.20}_{-0.20}$ kpc. The relation between millimeter-wave size and AGN fraction suggests that this size may be related to the evolutionary stage of the SMG. The very compact sizes for composite star-forming/AGN systems could be explained by supermassive black holes growing rapidly during the SMG coalescing, star-formation phase.

The analysis of galaxies on the star formation rate - stellar mass (SFR-M*) plane is a powerful diagnostic for galaxy evolution at different cosmic times. We consider a sample of 24463 galaxies from the CANDELS/GOODS-S survey to conduct a detailed analysis of the SFR-M* relation at redshifts 0.5$\leqslant z<$3 over more than three dex in stellar mass. To obtain SFR estimates, we utilise mid- and far-IR photometry when available, and rest-UV fluxes for all the other galaxies. We perform our analysis in different redshift bins, with two different methods: 1) a linear regression fitting of all star-forming galaxies, defined as those with specific star formation rates $\rm log_{10}(sSFR/yr^{-1}) > -9.8$, similarly to what is typically done in the literature; 2) a multi-Gaussian decomposition to identify the galaxy main sequence (MS), the starburst sequence and the quenched galaxy cloud. We find that the MS slope becomes flatter when higher stellar mass cuts are adopted, and that the apparent slope change observed at high masses depends on the SFR estimation method. In addition, the multi-Gaussian decomposition reveals the presence of a starburst population which increases towards low stellar masses and high redshifts. We find that starbursts make up ~5% of all galaxies at z=0.5-1.0, while they account for ~16% of galaxies at 2$<z<$3 with log$_{10}(M^{*})=$8.25-11.25. We conclude that the dissection of the SFR-M* in multiple components over a wide range of stellar masses is necessary to understand the importance of the different modes of star formation through cosmic time.

We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-$K_{S} \sim 25.3-24.8$~mag, $5\sigma$) NIR data from the UltraVISTA DR3 release, deep ground based optical imaging from the CFHTLS and Subaru Suprime Cam programs and Spitzer/IRAC mosaics combining observations from the SMUVS and SPLASH programs. Through the refined spectral energy distributions, which now also include new HyperSuprime Cam g, r, i, z and Y band data, we confirm that 3/5 galaxies have robust $z_{phot}\sim8.0-8.7$, consistent with the initial selection. The remaining 2/5 galaxies have a nominal $z_{phot}\sim2$. However, if we use the HST data alone, these objects have increased probability of being at $z\sim9$. Furthermore, we measure mean UV continuum slopes $\beta=-1.91\pm0.26$ for the three $z\sim8-9$ galaxies, marginally bluer than similarly luminous $z\sim4-6$ in CANDELS but consistent with previous measurements of similarly luminous galaxies at $z\sim7$. The circularized effective radius for our brightest source is $0.9\pm0.2$ kpc, similar to previous measurements for a bright $z\sim11$ galaxy and bright $z\sim7$ galaxies. Finally, enlarging our sample to include the six brightest $z\sim8$ LBGs identified over UltraVISTA (i.e., including three other sources from Labbe et al. 2017, in prep.) we estimate for the first time the volume density of galaxies at the extreme bright ($M_{UV}\sim-22$~mag) end of the $z\sim8$ UV LF. Despite this exceptional result, the still large statistical uncertainties do not allow us to discriminate between a Schechter and a double power-law form.

We present a new strong lensing analysis of the galaxy cluster MACS J1206.2-0847 (MACS 1206), at z=0.44, using deep spectroscopy from CLASH-VLT and VLT/MUSE archival data in combination with imaging from the Cluster Lensing and Supernova survey with Hubble. MUSE observations enable the spectroscopic identification of 23 new multiply imaged sources, extending the previous compilations by a factor of approximately five. In total, we use the positional measurements of 82 spectroscopic multiple images belonging to 27 families at z=1.0-6.1 to reconstruct the projected total mass distribution of MACS 1206. Remarkably, 11 multiple images are found within 50 kpc of the brightest cluster galaxy, making this an unprecedented set of constraints for the innermost projected mass distribution of a galaxy cluster. We thus find that, although dynamically relaxed, the smooth matter component (dark matter plus hot gas) of MACS 1206 shows a significant asymmetry, which closely follows the asymmetric distribution of the stellar component (galaxy members and intracluster light). We determine the value of the innermost logarithmic slope of the projected total mass density profile and find it to be close to the canonical Navarro-Frenk-White value. We demonstrate that this quantity is very robust against different parametrizations of the diffuse mass component; however, this is not the case when only one central image is used in the mass reconstruction. We also show that the mass density profile from our new strong lensing model is in very good agreement with dynamical and X-ray measurements at larger radii, where they overlap.

We report the spectroscopic confirmation of 22 new multiply lensed sources behind the Hubble Frontier Field (HFF) galaxy cluster MACS~J0416.1$-$2403 (MACS 0416), using archival data from the Multi Unit Spectroscopic Explorer (MUSE) on the VLT. Combining with previous spectroscopic measurements of 15 other multiply imaged sources, we obtain a sample of 102 secure multiple images with measured redshifts, the largest to date in a single strong lensing system. The newly confirmed sources are largely low-luminosity Lyman-$\alpha$ emitters with redshift in the range [3.08-6.15]. With such a large number of secure constraints, and a significantly improved sample of galaxy members in the cluster core, we have improved our previous strong lensing model and obtained a robust determination of the projected total mass distribution of MACS 0416. We find evidence of three cored dark-matter halos, adding to the known complexity of this merging system. The total mass density profile, as well as the sub-halo population, are found in good agreement with previous works. We update and make public the redshift catalog of MACS 0416 from our previous spectroscopic campaign with the new MUSE redshifts. We also release lensing maps (convergence, shear, magnification) in the standard HFF format.

We present physical properties [redshifts (z), star-formation rates (SFRs) and stellar masses (Mstar)] of bright (S850>4mJy) submm galaxies in the ~2deg2 COSMOS and UDS fields selected with SCUBA-2/JCMT. We complete the galaxy identification process for all (~2000) S/N>3.5 850-um sources, but focus our scientific analysis on a high-quality sub-sample of 651 S/N>4 sources with complete multi-wavelength coverage including 1.1-mm imaging. We check the reliability of our identifications, and the robustness of the SCUBA-2 fluxes by revisiting the recent ALMA follow-up of 29 sources in our sample. Considering >4mJy ALMA sources, our identification method has a completeness of ~86 per cent with a reliability of ~92 per cent, and only ~15-20 per cent of sources are significantly affected by multiplicity (when a secondary component contributes >1/3 of the primary source flux). The impact of source blending on the 850-um source counts as determined with SCUBA-2 is modest; scaling the single-dish fluxes by ~0.9 reproduces the ALMA source counts. For our final SCUBA-2 sample we find median z=2.40+0.10-0.04, SFR=287+-6Moyr-1, and log(Mstar/Mo)=11.12+-0.02 (the latter for 349/651 sources with optical identifications). These properties clearly locate bright submm galaxies on the high-mass end of the 'main sequence' of star-forming galaxies out to z~6, suggesting that major mergers are not a dominant driver of the high-redshift submm-selected population. Their number densities are also consistent with the evolving galaxy stellar mass function. Hence, the submm galaxy population is as expected, albeit reproducing the evolution of the main sequence of star-forming galaxies remains a challenge for theoretical models/simulations.

We study a large galaxy sample from the Spitzer Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS) to search for sources with enhanced 3.6 micron fluxes indicative of strong Halpha emission at z=3.9-4.9. We find that the percentage of "Halpha excess" sources reaches 37-40% for galaxies with stellar masses log10(M*/Msun) ~ 9-10, and decreases to <20% at log10(M*/Msun) ~ 10.7. At higher stellar masses, however, the trend reverses, although this is likely due to AGN contamination. We derive star formation rates (SFR) and specific SFR (sSFR) from the inferred Halpha equivalent widths (EW) of our "Halpha excess" galaxies. We show, for the first time, that the "Halpha excess" galaxies clearly have a bimodal distribution on the SFR-M* plane: they lie on the main sequence of star formation (with log10(sSFR/yr^{-1})<-8.05) or in a starburst cloud (with log10(sSFR/yr^{-1}) >-7.60). The latter contains ~15% of all the objects in our sample and accounts for >50% of the cosmic SFR density at z=3.9-4.9, for which we derive a robust lower limit of 0.066 Msun yr^{-1} Mpc^{-3}. Finally, we identify an unusual >50sigma overdensity of z=3.9-4.9 galaxies within a 0.20 x 0.20 sq. arcmin region. We conclude that the SMUVS unique combination of area and depth at mid-IR wavelengths provides an unprecedented level of statistics and dynamic range which are fundamental to reveal new aspects of galaxy evolution in the young Universe.

Imaging with the James Webb Space Telescope (JWST) will allow for observing the bulk of distant galaxies at the epoch of reionisation. The recovery of their properties, such as age, color excess E(B-V), specific star formation rate (sSFR) and stellar mass, will mostly rely on spectral energy distribution fitting, based on the data provided by JWST's two imager cameras, namely the Near Infrared Camera (NIRCam) and the Mid Infrared Imager (MIRI). In this work we analyze the effect of choosing different combinations of NIRCam and MIRI broad-band filters, from 0.6 {\mu}m to 7.7 {\mu}m, on the recovery of these galaxy properties. We performed our tests on a sample of 1542 simulated galaxies, with known input properties, at z=7-10. We found that, with only 8 NIRCam broad-bands, we can recover the galaxy age within 0.1 Gyr and the color excess within 0.06 mag for 70% of the galaxies. Besides, the stellar masses and sSFR are recovered within 0.2 and 0.3 dex, respectively, at z=7-9. Instead, at z=10, no NIRCam band traces purely the {\lambda}> 4000 {\AA} regime and the percentage of outliers in stellar mass (sSFR) increases by > 20% (> 90%), in comparison to z=9. The MIRI F560W and F770W bands are crucial to improve the stellar mass and the sSFR estimation at z=10. When nebular emission lines are present, deriving correct galaxy properties is challenging, at any redshift and with any band combination. In particular, the stellar mass is systematically overestimated in up to 0.3 dex on average with NIRCam data alone and including MIRI observations improves only marginally the estimation.

We study the spectrophotometric properties of a highly magnified (\mu~40-70) pair of stellar systems identified at z=3.2222 behind the Hubble Frontier Field galaxy cluster MACS~J0416. Five multiple images (out of six) have been spectroscopically confirmed by means of VLT/MUSE and VLT/X-Shooter observations. Each image includes two faint (m_uv~30.6), young (<100 Myr), low-mass (<10^7 Msun), low-metallicity (12+Log(O/H)~7.7, or 1/10 solar) and compact (30 pc effective radius) stellar systems separated by ~300pc, after correcting for lensing amplification. We measured several rest-frame ultraviolet and optical narrow (\sigma_v <~ 25 km/s) high-ionization lines. These features may be the signature of very hot (T>50000 K) stars within dense stellar clusters, whose dynamical mass is likely dominated by the stellar component. Remarkably, the ultraviolet metal lines are not accompanied by Lya emission (e.g., CIV / Lya > 15), despite the fact that the Lya line flux is expected to be 150 times brighter (inferred from the Hbeta flux). A spatially-offset, strongly-magnified (\mu>50) Lya emission with a spatial extent <~7.6 kpc^2 is instead identified 2 kpc away from the system. The origin of such a faint emission can be the result of fluorescent Lya induced by a transverse leakage of ionizing radiation emerging from the stellar systems and/or can be associated to an underlying and barely detected object (with m_uv > 34 de-lensed). This is the first confirmed metal-line emitter at such low-luminosity and redshift without Lya emission, suggesting that, at least in some cases, a non-uniform covering factor of the neutral gas might hamper the Lya detection.

We analyze a sample of $z$-dropout galaxies in the CANDELS GOODS South and UDS fields that have been targeted by a dedicated spectroscopic campaign aimed at detecting their Ly$\alpha$ line. Deep IRAC observations at 3.6 and 4.5 $\mu$m are used to determine the strength of optical emission lines affecting these bands at z$\sim$6.5-6.9 in order to i) investigate possible physical differences between Ly$\alpha$ emitting and non-emitting sources; ii) constrain the escape fraction of ionizing photons; iii) provide an estimate of the specific star-formation rate at high redshifts. We find evidence of strong [OIII]+H$\beta$ emission in the average (stacked) SEDs of galaxies both with and without Ly$\alpha$ emission. The blue IRAC [3.6]-[4.5] color of the stack with detected Ly$\alpha$ line can be converted into a rest-frame equivalent width EW([OIII]+H$\beta$)=1500$^{+530}_{-440}\AA$ assuming a flat intrinsic stellar continuum. This strong optical line emission enables a first estimate of f$_{esc}\lesssim$20% on the escape fraction of ionizing photons from Ly$\alpha$ detected objects. The objects with no Ly$\alpha$ line show less extreme EW([OIII]+H$\beta$)=520$^{+170}_{-150}\AA$ suggesting different physical conditions of the HII regions with respect to Ly$\alpha$-emitting ones, or a larger f$_{esc}$. The latter case is consistent with a combined evolution of f$_{esc}$ and the neutral hydrogen fraction as an explanation of the lack of bright Ly$\alpha$ emission at z$>$6. A lower limit on the specific star formation rate, SSFR$>$9.1$Gyr^{-1}$ for $M_{star}=2 \times 10^9 M_{\odot}$ galaxies at these redshifts can be derived from the spectroscopically confirmed sample.

In this paper, we use stacking analysis to trace the mass-growth, colour evolution, and structural evolution of present-day massive galaxies ($\log(M_{*}/M_{\odot})=11.5$) out to $z=5$. We utilize the exceptional depth and area of the latest UltraVISTA data release, combined with the depth and unparalleled seeing of CANDELS to gather a large, mass-selected sample of galaxies in the NIR (rest-frame optical to UV). Progenitors of present-day massive galaxies are identified via an evolving cumulative number density selection, which accounts for the effects of merging to correct for the systematic biases introduced using a fixed cumulative number density selection, and find progenitors grow in stellar mass by $\approx1.5~\mathrm{dex}$ since $z=5$. Using stacking, we analyze the structural parameters of the progenitors and find that most of the stellar mass content in the central regions was in place by $z\sim2$, and while galaxies continue to assemble mass at all radii, the outskirts experience the largest fractional increase in stellar mass. However, we find evidence of significant stellar mass build up at $r<3~\mathrm{kpc}$ beyond $z>4$ probing an era of significant mass assembly in the interiors of present day massive galaxies. We also compare mass assembly from progenitors in this study to the EAGLE simulation and find qualitatively similar assembly with $z$ at $r<3~\mathrm{kpc}$. We identify $z\sim1.5$ as a distinct epoch in the evolution of massive galaxies where progenitors transitioned from growing in mass and size primarily through in-situ star formation in disks to a period of efficient growth in $r_{e}$ consistent with the minor merger scenario.

We report on five compact, extremely young (<10Myr) and blue (\beta_UV<-2.5, F_\lambda =\lambda^\beta) objects observed with VLT/MUSE at redshift 3.1169, 3.235, in addition to three objects at z=6.145. These sources are magnified by the Hubble Frontier Field galaxy clusters MACS~J0416 and AS1063. Their de-lensed half light radii (Re) are between 16 to 140pc, the stellar masses are ~1-20 X 10^6 Msun, the magnitudes are m_uv=28.8 - 31.4 (-17<Muv<-15) and specific star formation rates can be as large as ~800Gyr^-1. Multiple images of these systems are widely separated in the sky (up to 50'') and individually magnified by factors 3-40. Remarkably, the inferred physical properties of two objects are similar to those expected in some globular cluster formation scenarios, representing the best candidate proto-globular clusters (proto-GC) discovered so far. Rest-frame optical high dispersion spectroscopy of one of them at z=3.1169 yields a velocity dispersion \sigma_v~20km/s, implying a dynamical mass dominated by the stellar mass. Another object at z=6.145, with de-lensed Muv ~ -15.3 (m_uv ~ 31.4), shows a stellar mass and a star-formation rate surface density consistent with the values expected from popular GC formation scenarios. An additional star-forming region at z=6.145, with de-lensed m_uv ~ 32, a stellar mass of 0.5 X 10^6 Msun and a star formation rate of 0.06 Msun/yr is also identified. These objects currently represent the faintest spectroscopically confirmed star-forming systems at z>3, elusive even in the deepest blank fields. We discuss how proto-GCs might contribute to the ionization budget of the universe and augment Lya visibility during reionization. This work underlines the crucial role of JWST in characterizing the rest-frame optical and near-infrared properties of such low-luminosity high-z objects.

We present the detailed characterization of two extremely red submillimeter galaxies (SMGs), ASXDF1100.053.1 and 231.1, with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). These SMGs were selected originally using AzTEC at 1100 micron, and are observed by Herschel to be faint at 100--500 micron. Their (sub)millimeter colors are as red as -- or redder -- than known z>~5 SMGs; indeed, ASXDF1100.053.1 is redder than HFLS 3, which lies at z=6.3. They are also faint and red in the near-/mid-infrared: ~1 microJy at IRAC 4.5 micron and <0.2 microJy in the Ks filter. These SMGs are also faint in the radio waveband, where F_6GHz=4.5 microJy for ASXDF1100.053.1 and F_1.4GHz=28 microJy for ASXDF1100.231.1, suggestive of z=6.5^{+1.4}_{-1.1} and z=4.1^{+0.6}_{-0.7} for ASXDF1100.053.1 and 231.1, respectively. ASXDF1100.231.1 has a flux excess in the 3.6-micron filter, probably due to H$\alpha$ emission at z=4--5. Derived properties of ASXDF1100.053.1 for z=5.5--7.5 and 231.1 for z=3.5--5.5 are as follows: their infrared luminosities are [6.5-7.4]x10^{12} and [4.2-4.5]x10^{12} L_sun; their stellar masses are [0.9-2]x10^{11} and [0.4-3]x10^{10} M_sun; their circularized half-light radii in the ALMA maps are ~1 and <~0.2 kpc (~2--3 kpc for 90% of the total flux). Lastly, their surface infrared luminosity densities, Sigma_IR, are ~1x10^{12} and >~1.5x10^{13} L_sun kpc^{-2}, similar to values seen for local (U)LIRGs. These data suggest that ASXDF1100.053.1 and 231.1 are compact SMGs at z>~4 and can plausibly evolve into z>~3 compact quiescent galaxies.

The determination of galaxy redshifts in James Webb Space Telescope (JWST)'s blank-field surveys will mostly rely on photometric estimates, based on the data provided by JWST's Near-Infrared Camera (NIRCam) at 0.6-5.0 {\mu}m and Mid Infrared Instrument (MIRI) at {\lambda}>5.0 {\mu}m. In this work we analyse the impact of choosing different combinations of NIRCam and MIRI broad-band filters (F070W to F770W), as well as having ancillary data at {\lambda}<0.6 {\mu}m, on the derived photometric redshifts (zphot) of a total of 5921 real and simulated galaxies, with known input redshifts z=0-10. We found that observations at {\lambda}<0.6 {\mu}m are necessary to control the contamination of high-z samples by low-z interlopers. Adding MIRI (F560W and F770W) photometry to the NIRCam data mitigates the absence of ancillary observations at {\lambda}<0.6 {\mu}m and improves the redshift estimation. At z=7-10, accurate zphot can be obtained with the NIRCam broad bands alone when S/N>=10, but the zphot quality significantly degrades at S/N<=5. Adding MIRI photometry with one magnitude brighter depth than the NIRCam depth allows for a redshift recovery of 83-99%, depending on SED type, and its effect is particularly noteworthy for galaxies with nebular emission. The vast majority of NIRCam galaxies with [F150W]=29 AB mag at z=7-10 will be detected with MIRI at [F560W, F770W]<28 mag if these sources are at least mildly evolved or have spectra with emission lines boosting the mid-infrared fluxes.

We present stacking analyses on our ALMA deep 1.1 mm imaging in the SXDF using 1.6 {\mu}m and 3.6 {\mu}m selected galaxies in the CANDELS WFC3 catalog. We detect a stacked flux of ~0.03-0.05 mJy, corresponding to LIR < 10^11 Lsun and a star formation rate (SFR) of ~ 15 Msun/yr at z = 2. We find that galaxies brighter in the rest-frame near-infrared tend to be also brighter at 1.1 mm, and galaxies fainter than m[3.6um] = 23 do not produce detectable 1.1 mm emission. This suggests a correlation between stellar mass and SFR, but outliers to this correlation are also observed, suggesting strongly boosted star formation or extremely large extinction. We also find tendencies that redder galaxies and galaxies at higher redshifts are brighter at 1.1 mm. Our field contains z ~ 2.5 H-alpha emitters and a bright single-dish source. However, we do not find evidence of bias in our results caused by the bright source. By combining the fluxes of sources detected by ALMA and fluxes of faint sources detected with stacking, we recover a 1.1 mm surface brightness of up to 20.3 +/- 1.2 Jy/deg, comparable to the extragalactic background light measured by COBE. Based on the fractions of optically faint sources in our and previous ALMA studies and the COBE measurements, we find that approximately half of the cosmic star formation may be obscured by dust and missed by deep optical surveys, Much deeper and wider ALMA imaging is therefore needed to better constrain the obscured cosmic star formation history.

We report the discovery of a multiply lensed Lyman-$\alpha$ blob (LAB) behind the galaxy cluster AS1063 using the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT). The background source is at $z=$ 3.117 and is intrinsically faint compared to almost all previously reported LABs. We used our highly precise strong lensing model to reconstruct the source properties, and we find an intrinsic luminosity of $L_{\rm Ly\alpha}$=$1.9\times10^{42}$ erg s$^{-1}$, extending to 33 kpc. We find that the LAB is associated with a group of galaxies, and possibly a protocluster, in agreement with previous studies that find LABs in overdensities. In addition to Lyman-$\alpha$ (Ly$\alpha$) emission, we find \ion{C}{IV}, \ion{He}{II}, and \ion{O}{III}] ultraviolet (UV) emission lines arising from the centre of the nebula. We used the compactness of these lines in combination with the line ratios to conclude that the \Lya nebula is likely powered by embedded star formation. Resonant scattering of the \Lya photons then produces the extended shape of the emission. Thanks to the combined power of MUSE and strong gravitational lensing, we are now able to probe the circumgalatic medium of sub-$L_{*}$ galaxies at $z\approx 3$.

In spite of their conjectured importance for the Epoch of Reionization, the properties of low-mass galaxies are currently still under large debate. In this article, we study the stellar and gaseous properties of faint, low-mass galaxies at z>3. We observed the Frontier Fields cluster Abell S1063 with MUSE over a 2 arcmin^2 field, and combined integral-field spectroscopy with gravitational lensing to perform a blind search for intrinsically faint Lya emitters (LAEs). We determined in total the redshift of 172 galaxies of which 14 are lensed LAEs at z=3-6.1. We increased the number of spectroscopically-confirmed multiple-image families from 6 to 17 and updated our gravitational-lensing model accordingly. The lensing-corrected Lya luminosities are with L(Lya) <= 10^41.5 erg/s among the lowest for spectroscopically confirmed LAEs at any redshift. We used expanding gaseous shell models to fit the Lya line profile, and find low column densities and expansion velocities. This is to our knowledge the first time that gaseous properties of such faint galaxies at z>=3 are reported. We performed SED modelling to broadband photometry from the U-band through the infrared to determine the stellar properties of these LAEs. The stellar masses are very low (10^{6-8} Msun), and are accompanied by very young ages of 1-100 Myr. The very high specific star formation rates (~100 Gyr^-1) are characteristic of starburst galaxies. The UV-continuum slopes beta are low in our sample, with beta<-2 for all galaxies with M_star<10^8 Msun. We conclude that our low-mass galaxies at 3<z<6 are forming stars at higher rates when correcting for stellar mass effects than seen locally or in more massive galaxies. The young stellar populations with high star-formation rates and low HI column densities lead to continuum slopes and LyC-escape fractions expected for a scenario where low mass galaxies reionise the Universe.

We report on the discovery of extended Ly-alpha nebulae at z~3.3 in the Hubble Ultra Deep Field (HUDF, ~ 40 kpc X 80 kpc) and behind the Hubble Frontier Fields galaxy cluster MACSJ0416 (~ 40kpc), spatially associated with groups of star-forming galaxies. VLT/MUSE integral field spectroscopy reveals a complex structure with a spatially-varying double peaked Ly-alpha emission. Overall, the spectral profiles of the two Ly-alpha nebulae are remarkably similar, both showing a prominent blue emission, more intense and slightly broader than the red peak. From the first nebula, located in the HUDF, no X-ray emission has been detected, disfavoring the possible presence of AGNs. Spectroscopic redshifts have been derived for 11 galaxies within two arcsec from the nebula and spanning the redshift range 1.037<z<5.97. The second nebula, behind MACSJ0416, shows three aligned star-forming galaxies plausibly associated to the emitting gas. In both systems, the associated galaxies reveal possible intense rest-frame-optical nebular emissions lines [OIII]4959-5007+Hbeta with equivalent widths as high as 1500A rest-frame and star formation rates ranging from a few to tens of solar masses per year. A possible scenario is that of a group of young, star-forming galaxies sources of escaping ionising radiation that induce Ly-alpha fluorescence, therefore revealing the kinematics of the surrounding gas. Also Ly-alpha powered by star-formation and/or cooling radiation may resemble the double peaked spectral properties and the morphology observed here. If the intense blue emission is associated with inflowing gas, then we may be witnessing an early phase of galaxy or a proto-cluster (or group) formation.

We present VLT/X-Shooter and MUSE spectroscopy of an faint F814W=28.60+/-0.33 (Muv=-17.0), low mass (~<10^7 Msun) and compact (Reff=62pc) freshly star-forming galaxy at z=3.1169 magnified (16x) by the Hubble Frontier Fields galaxy cluster Abell S1063. Gravitational lensing allows for a significant jump toward low-luminosity regimes, in moderately high resolution spectroscopy (R=lambda/dlambda ~ 3000-7400). We measured CIV1548,1550, HeII1640, OIII]1661,1666, CIII]1907,1909, Hbeta, [OIII]4959,5007, emission lines with FWHM< 50 km/s and (de-lensed) fluxes spanning the interval 1.0x10^-19 - 2.0x10^-18 erg/s/cm2 at S/N=4-30. The double peaked Lya emission with Delta_v(red-blue) = 280(+/-7)km/s and de-lensed fluxes 2.4_(blue)|8.5_(red)x10^-18 erg/s/cm2 (S/N=38_(blue)|110_(red)) indicate a low column density of neutral hydrogen gas consistent with a highly ionized interstellar medium as also inferred from the large [OIII]5007/[OII]3727>10 ratio. We detect CIV1548,1550 resonant doublet in emission, each component with FWHM ~< 45 km/s, and redshifted by +51(+/-10)km/s relative to the systemic redshift. We interpret this as nebular emission tracing an expanding optically-thin interstellar medium. Both CIV1548,1550 and HeII1640 suggest the presence of hot and massive stars (with a possible faint AGN). The ultraviolet slope is remarkably blue, beta =-2.95 +/- 0.20 (F_lambda=lambda^beta), consistent with a dust-free and young ~<20 Myr galaxy. Line ratios suggest an oxygen abundance 12+log(O/H)<7.8. We are witnessing an early episode of star-formation in which a relatively low NHI and negligible dust attenuation might favor a leakage of ionizing radiation. This galaxy currently represents a unique low-luminosity reference object for future studies of the reionization epoch with JWST.

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

We report 1.1 mm number counts revealed with the Atacama Large Millimeter/submillimeter Array (ALMA) in the Subaru/XMM-Newton Deep Survey Field (SXDF). The advent of ALMA enables us to reveal millimeter-wavelength number counts down to the faint end without source confusion. However, previous studies are based on the ensemble of serendipitously-detected sources in fields originally targeting different sources and could be biased due to the clustering of sources around the targets. We derive number counts in the flux range of 0.2-2 mJy by using 23 (>=4sigma) sources detected in a continuous 2.0 arcmin$^2$ area of the SXDF. The number counts are consistent with previous results within errors, suggesting that the counts derived from serendipitously-detected sources are not significantly biased, although there could be field-to-field variation due to the small survey area. By using the best-fit function of the number counts, we find that ~40% of the extragalactic background light at 1.1 mm is resolved at S(1.1mm) > 0.2 mJy.

We present MUSE observations in the core of the HFF galaxy cluster MACS J1149.5+2223, where the first magnified and spatially-resolved multiple images of SN 'Refsdal' at redshift 1.489 were detected. Thanks to a DDT program with the VLT and the extraordinary efficiency of MUSE, we measure 117 secure redshifts with just 4.8 hours of total integration time on a single target pointing. We spectroscopically confirm 68 galaxy cluster members, with redshift values ranging from 0.5272 to 0.5660, and 18 multiple images belonging to 7 background, lensed sources distributed in redshifts between 1.240 and 3.703. Starting from the combination of our catalog with those obtained from extensive spectroscopic and photometric campaigns using the HST, we select a sample of 300 (164 spectroscopic and 136 photometric) cluster members, within approximately 500 kpc from the BCG, and a set of 88 reliable multiple images associated to 10 different background source galaxies and 18 distinct knots in the spiral galaxy hosting SN 'Refsdal'. We exploit this valuable information to build 6 detailed strong lensing models, the best of which reproduces the observed positions of the multiple images with a rms offset of only 0.26". We use these models to quantify the statistical and systematic errors on the predicted values of magnification and time delay of the next emerging image of SN 'Refsdal'. We find that its peak luminosity should should occur between March and June 2016, and should be approximately 20% fainter than the dimmest (S4) of the previously detected images but above the detection limit of the planned HST/WFC3 follow-up. We present our two-dimensional reconstruction of the cluster mass density distribution and of the SN 'Refsdal' host galaxy surface brightness distribution. We outline the roadmap towards even better strong lensing models with a synergetic MUSE and HST effort.

We have conducted 1.1 mm ALMA observations of a contiguous $105'' \times 50''$ or 1.5 arcmin$^2$ window in the SXDF-UDS-CANDELS. We achieved a 5$\sigma$ sensitivity of 0.28 mJy, providing a flat sensus of dusty star-forming galaxies with $L_{\rm IR} \sim6\times10^{11}$ $L_\odot$ (for $T_{\rm dust}$ =40K) up to $z\sim10$ thanks to the negative K-correction at this wavelength. We detected 5 brightest sources (S/N$>$6) and 18 low-significance sources (5$>$S/N$>$4; these may contain spurious detections, though). One of the 5 brightest ALMA sources ($S_{\rm 1.1mm} = 0.84 \pm 0.09$ mJy) is extremely faint in the WFC3 and VLT/HAWK-I images, demonstrating that a contiguous ALMA imaging survey is able to uncover a faint dust-obscured population that is invisible in deep optical/near-infrared surveys. We found a possible [CII]-line emitter at $z=5.955$ or a low-$z$ CO emitting galaxy within the field, which may allow us to constrain the [CII] and/or the CO luminosity functions across the history of the universe.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

Supernova (SN) Refsdal is the first multiply-imaged, highly-magnified, and spatially-resolved SN ever observed. The SN exploded in a highly-magnified spiral galaxy at z=1.49 behind the Frontier Fields Cluster MACS1149, and provides a unique opportunity to study the environment of SNe at high z. We exploit the time delay between multiple images to determine the properties of the SN and its environment, before, during, and after the SN exploded. We use the integral-field spectrograph MUSE on the VLT to simultaneously target all observed and model-predicted positions of SN Refsdal. We find MgII emission at all positions of SN Refsdal, accompanied by weak FeII* emission at two positions. The measured ratios of [OII] to MgII emission of 10-20 indicate a high degree of ionization with low metallicity. Because the same high degree of ionization is found in all images, and our spatial resolution is too coarse to resolve the region of influence of SN Refsdal, we conclude that this high degree of ionization has been produced by previous SNe or a young and hot stellar population. We find no variability of the [OII] line over a period of 57 days. This suggests that there is no variation in the [OII] luminosity of the SN over this period, or that the SN has a small contribution to the integrated [OII] emission over the scale resolved by our observations.

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.

We have analysed a sample of 574 Spitzer 4.5 micron-selected galaxies with [4.5]<23 and Ks_auto>24 (AB) over the UltraVISTA ultra-deep COSMOS field. Our aim is to investigate whether these mid-IR bright, near-IR faint sources contribute significantly to the overall population of massive galaxies at redshifts z>=3. By performing a spectral energy distribution (SED) analysis using up to 30 photometric bands, we have determined that the redshift distribution of our sample peaks at redshifts z~2.5-3.0, and ~32% of the galaxies lie at z>=3. We have studied the contribution of these sources to the galaxy stellar mass function (GSMF) at high redshifts. We found that the [4.5]<23, Ks_auto>24 galaxies produce a negligible change to the GSMF previously determined for Ks_auto<24 sources at 3=<z<4, but their contribution is more important at 4=<z<5, accounting for >~50% of the galaxies with stellar masses Mst>~6 x 10^10 Msun. We also constrained the GSMF at the highest-mass end (Mst>~2 x 10^11 Msun) at z>=5. From their presence at 5=<z<6, and virtual absence at higher redshifts, we can pinpoint quite precisely the moment of appearance of the first most massive galaxies as taking place in the ~0.2 Gyr of elapsed time between z~6 and z~5. Alternatively, if very massive galaxies existed earlier in cosmic time, they should have been significantly dust-obscured to lie beyond the detection limits of current, large-area, deep near-IR surveys.

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

We report the source size distribution, as measured by ALMA millimetric continuum imaging, of a sample of 13 AzTEC-selected submillimeter galaxies (SMGs) at z_photo ~ 3-6. Their infrared luminosities and star-formation rates (SFR) are L_IR ~ 2-6 x 10^12 L_sun and ~ 200-600 M_sun yr-1, respectively. The size of z ~ 3-6 SMGs ranges from 0".10 to 0".38 with a median of 0".20+0".03-0".05 (FWHM), corresponding to a median circularized effective radius (Rc,e) of 0.67+0.13-0.14 kpc, comparable to the typical size of the stellar component measured in compact quiescent galaxies at z ~ 2 (cQGs) --- R ~ 1 kpc. The median surface SFR density of our z ~ 3-6 SMGs is 100+42-26 M_sun yr-1 kpc-2, comparable to that seen in local merger-driven (U)LIRGsrather than in extended disk galaxies at low and high redshifts. The discovery of compact starbursts in z >~ 3 SMGs strongly supports a massive galaxy formation scenario wherein z ~ 3-6 SMGs evolve into the compact stellar components of z ~ 2 cQGs. These cQGs are then thought to evolve into the most massive ellipticals in the local Universe, mostly via dry mergers. Our results thus suggest that z >~ 3 SMGs are the likely progenitors of massive local ellipticals, via cQGs, meaning that we can now trace the evolutionary path of the most massive galaxies over a period encompassing ~ 90% of the age of the Universe.

We present a thorough characterization of a large sample of 183 extreme emission-line galaxies (EELGs) at redshift 0.11 < z < 0.93 selected from the 20k zCOSMOS Bright Survey because of their unusually large emission line equivalent widths. We use multiwavelength COSMOS photometry, HST-ACS I-band imaging and optical zCOSMOS spectroscopy to derive the main global properties of EELGs, such as sizes, masses, SFRs, reliable metallicities from both "direct" and "strong-line" methods. The EELGs are compact (R_50 ~ 1.3 kpc), low-mass (log(M*/Msol)~7-10) galaxies forming stars at unusually high specific SFR (log(sSFR/yr) up to ~ -7) compared to main sequence SFGs of the same stellar mass and redshift. At UV wavelengths, the EELGs are luminous and show high surface brightness and include strong Ly$\alpha$ emitters, as revealed by GALEX spectroscopy. We show that zCOSMOS EELGs are high-ionization, low-metallicity systems, with median 12+log(O/H)=8.16, including a handful of extremely metal-deficient galaxies (<10% solar). While ~80% of the EELGs show non-axisymmetric morphologies, including clumpy and tadpole galaxies, we find that ~29% of them show additional low surface-brightness features, which strongly suggest recent or ongoing interactions. As star-forming dwarfs in the local Universe, EELGs are most often found in relative isolation. While only very few EELGs belong to compact groups, almost one third of them are found in spectroscopically confirmed loose pairs or triplets. We conclude that EELGs are galaxies caught in a transient and probably early period of their evolution, where they are efficiently building-up a significant fraction of their present-day stellar mass in an ongoing galaxy-wide starburst. Therefore, the EELGs constitute an ideal benchmark for comparison studies between low- and high-redshift low-mass star-forming galaxies.

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

We use percent-level precision photometric redshifts in the UltraVISTA-DR1 near-infrared survey to investigate the changing relationship be- tween galaxy stellar mass and the dark matter haloes hosting them to $z\sim2$. We achieve this by measuring the clustering properties and abundances of a series of volume-limited galaxy samples selected by stellar mass and star-formation activity. We interpret these results in the framework of a phenomenological halo model and numerical simulations. Our measurements span a uniquely large range in stellar mass and redshift and reach below the characteristic stellar mass to $z\sim2$. Our results are: 1. At fixed redshift and scale, clustering amplitude depends monotonically on sample stellar mass threshold; 2. At fixed angular scale, the projected clustering amplitude decreases with redshift but the co-moving correlation length remains constant; 3. Characteristic halo masses and galaxy bias increase with increasing median stellar mass of the sample; 4. The slope of these relationships is modified in lower mass haloes; 5. Concerning the passive galaxy population, characteristic halo masses are consistent with a simply less-abundant version of the full galaxy sample, but at lower redshifts the fraction of satellite galaxies in the passive population is very different from the full galaxy sample; 6. Finally we find that the ratio between the characteristic halo mass and median stellar mass at each redshift bin reaches a peak at ${\log(M_h/M_\odot)}\sim12.2$ and the position of this peak remains constant out to $z\sim2$. The behaviour of the full and passively evolving galaxy samples can be understood qualitatively by considering the slow evolution of the characteristic stellar mass in the redshift range probed by our survey.

We present the first observations of the Frontier Fields Cluster Abell S1063 taken with the newly commissioned Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph. Because of the relatively large field of view (1 arcmin^2), MUSE is ideal to simultaneously target multiple galaxies in blank and cluster fields over the full optical spectrum. We analysed the four hours of data obtained in the Science Verification phase on this cluster and measured redshifts for 53 galaxies. We confirm the redshift of five cluster galaxies, and determine the redshift of 29 other cluster members. Behind the cluster, we find 17 galaxies at higher redshift, including three previously unknown Lyman-alpha emitters at z>3, and five multiply-lensed galaxies. We report the detection of a new z=4.113 multiply lensed galaxy, with images that are consistent with lensing model predictions derived for the Frontier Fields. We detect C III], C IV, and He II emission in a multiply lensed galaxy at z=3.116, suggesting the likely presence of an active galactic nucleus. We also created narrow-band images from the MUSE datacube to automatically search for additional line emitters corresponding to high-redshift candidates, but we could not identify any significant detections other than those found by visual inspection. With the new redshifts, it will become possible to obtain an accurate mass reconstruction in the core of Abell S1063 through refined strong lensing modelling. Overall, our results illustrate the breadth of scientific topics that can be addressed with a single MUSE pointing. We conclude that MUSE is a very efficient instrument to observe galaxy clusters, enabling their mass modelling, and to perform a blind search for high-redshift galaxies.

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

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

The remarkable progress made in infrared (IR) astronomical instruments over the last 10-15 years has radically changed our vision of the extragalactic IR sky, and overall understanding of galaxy evolution. In particular, this has been the case for the study of active galactic nuclei (AGN), for which IR observations provide a wealth of complementary information that cannot be derived from data in other wavelength regimes. In this review, I summarize the unique contribution that IR astronomy has recently made to our understanding of AGN and their role in galaxy evolution, including both physical studies of AGN at IR wavelengths, and the search for AGN among IR galaxies in general. Finally, I identify and discuss key open issues that it should be possible to address with forthcoming IR telescopes.

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

Galaxy formation models invoke the presence of strong feedback mechanisms that regulate the growth of massive galaxies at high redshifts. In this paper we aim to: (1) confirm spectroscopically the redshifts of a sample of massive galaxies selected with photometric redshifts z > 2.5; (2) investigate the properties of their stellar and interstellar media; (3) detect the presence of outflows, and measure their velocities. To achieve this, we analysed deep, high-resolution (R~2000) FORS2 rest-frame UV spectra for 11 targets. We confirmed that 9 out of 11 have spectroscopic redshifts z > 2.5. We also serendipitously found two mask fillers at redshift z > 2.5, which originally were assigned photometric redshifts 2.0 < z < 2.5. In the four highest-quality spectra we derived outflow velocities by fitting the absorption line profiles with models including multiple dynamical components. We found strongly asymmetric, high-ionisation lines, from which we derived outflow velocities ranging from 480 to 1518 km/s. The two galaxies with highest velocity show signs of AGN. We revised the spectral energy distribution fitting U-band through 8 micron photometry, including the analysis of a power-law component subtraction to identify the possible presence of active galactic nuclei (AGN). The revised stellar masses of all but one of our targets are >1e10 Msun, with four having stellar masses > 5e10 Msun. Three galaxies have a significant power-law component in their spectral energy distributions, which indicates that they host AGN. We conclude that massive galaxies are characterised by significantly higher velocity outflows than the typical Lyman break galaxies at z ~ 3. The incidence of high-velocity outflows (~40% within our sample) is also much higher than among massive galaxies at z < 1, which is consistent with the powerful star formation and nuclear activity that most massive galaxies display at z > 2.

We studied the chronology of galactic bulge and disc formation by analysing the relative contributions of these components to the B-band rest-frame luminosity density at different epochs. We present the first estimate of the evolution of the fraction of rest-frame B-band light in galactic bulges and discs since redshift z~0.8. We performed a bulge-to-disc decomposition of HST/ACS images of 3266 galaxies in the zCOSMOS-bright survey with spectroscopic redshifts in the range 0.7 < z < 0.9. We find that the fraction of B-band light in bulges and discs is $(26 \pm 4)%$ and $(74 \pm 4)%$, respectively. When compared with rest-frame B-band measurements of galaxies in the local Universe in the same mass range ($10^{9} M_{\odot}\lessapprox M \lessapprox 10^{11.5} M_{\odot}$), we find that the B-band light in discs decreases by ~30% from z~0.7-0.9 to z~0, while the light from the bulge increases by ~30% over the same period of time. We interpret this evolution as the consequence of star formation and mass assembly processes, as well as morphological transformation, which gradually shift stars formed at half the age of the Universe from star-forming late-type/irregular galaxies toearlier types and ultimately into spheroids.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2<z<3 to explore whether they are progenitors of compact quiescent galaxies at z~2. The optical/NIR and far-IR Spitzer/Herschel colors indicate that most compact SFGs are heavily obscured. Nearly half (47%) host an X-ray bright AGN. In contrast, only about 10% of other massive galaxies at that time host AGNs. Compact SFGs have centrally-concentrated light profiles and spheroidal morphologies similar to quiescent galaxies, and are thus strikingly different from other SFGs. Most compact SFGs lie either within the SFR-M main sequence (65%) or below (30%), on the expected evolutionary path towards quiescent galaxies. These results show conclusively that galaxies become more compact before they lose their gas and dust, quenching star formation. Using extensive HST photometry from CANDELS and grism spectroscopy from the 3D-HST survey, we model their stellar populations with either exponentially declining (tau) star formation histories (SFHs) or physically-motivated SFHs drawn from semi-analytic models (SAMs). SAMs predict longer formation timescales and older ages ~2 Gyr, which are nearly twice as old as the estimates of the tau models. While both models yield good SED fits, SAM SFHs better match the observed slope and zero point of the SFR-M main sequence. Some low-mass compact SFGs (log(M/Msun)=10-10.6) have younger ages but lower sSFRs than that of more massive galaxies, suggesting that the low-mass galaxies reach the red sequence faster. If the progenitors of compact SFGs are extended SFGs, state-of-the-art SAMs show that mergers and disk instabilities are both able to shrink galaxies, but disk instabilities are more frequent (60% versus 40%) and form more concentrated galaxies. We confirm this result via high-resolution hydrodynamic simulations.

We identify 42 candidate groups lying between 1.8<z<3.0 from a sample of 3502 galaxies with spectroscopic redshifts in the zCOSMOS-deep redshift survey within the same redshift interval. These systems contain three to five spectroscopic galaxies that lie within 500kpc in projected distance (in physical space) and within 700km/s in velocity. Based on extensive analysis of mock catalogues that have been generated from the Millennium simulation, we examine the likely nature of these systems at the time of observation, and what they will evolve into down to the present epoch. Although few of the "member" galaxies are likely to reside in the same halo at the epoch we observe them, 50% of the systems will eventually bring them all into the same halo, and almost all (93%) will have at least part of the member galaxies in the same halo by the present epoch. Most of the candidate groups can therefore be described as "proto-groups". An estimate of the overdensities is also consistent with the idea that these systems are being seen at the start of the assembly process. We also examine present-day haloes and ask whether their progenitors would have been seen amongst our candidate groups. For present-day haloes between 10^14-10^15Msun/h, 35% should have appeared amongst our candidate groups, and this would have risen to 70% if our survey had been fully-sampled, so we can conclude that our sample can be taken as representative of a large fraction of such systems. There is a clear excess of massive galaxies above 10^10Msun around the locations of the candidate groups in a large independent COSMOS photo-z sample, but we see no evidence in this latter data for any colour differentiation with respect to the field. This is however consistent with the idea that such differentiation arises in satellite galaxies, as indicated at z<1, if the candidate groups are indeed only starting to be assembled.

We present the analysis of photometric, spectroscopic, and morphological properties for differently selected samples of passive galaxies up to z=1 extracted from the zCOSMOS-20k spectroscopic survey. This analysis intends to explore the dependence of galaxy properties on the selection criterion adopted, study the degree of contamination due to star-forming outliers, and provide a comparison between different commonly used selection criteria. We extracted from the zCOSMOS-20k catalog six different samples of passive galaxies, based on morphology, optical colors, specific star-formation rate, a best fit to the observed spectral energy distribution, and a criterion that combines morphological, spectroscopic, and photometric information. The morphological sample has the higher percentage of contamination in colors, specific star formation rate and presence of emission lines, while the red & passive ETGs sample is the purest, with properties mostly compatible with no star formation activity; however, it is also the less economic criterion in terms of information used. The best performing among the other criteria are the red SED and the quiescent ones, providing a percentage of contamination only slightly higher than the red & passive ETGs criterion (on average of a factor of ~2) but with absolute values of the properties of contaminants still compatible with a red, passively evolving population. We also provided two revised definitions of early type galaxies based on restframe color-color and color-mass criteria, that better reproduce the observed bimodalities. The analysis of the number densities shows evidences of mass-assembly downsizing, with galaxies at 10.25<log(M/Msun)<10.75 increasing their number by a factor ~2-4 from z=0.6 to z=0.2, by a factor ~2-3 from z=1 to z=0.2 at 10.75<log(M/Msun)<11, and by only ~10-50% from z=1 to z=0.2 at 11<log(M/Msun)<11.5.

We present an analysis of cool outflowing gas around galaxies, traced by MgII absorption lines in the co-added spectra of a sample of 486 zCOSMOS galaxies at 1 < z < 1.5. These galaxies span a range of stellar masses (9.45< log[M*/Msun]<10.7) and star formation rates (0.14 < log [SFR/Msun/yr] < 2.35). We identify the cool outflowing component in the MgII absorption and find that the equivalent width of the outflowing component increases with stellar mass. The outflow equivalent width also increases steadily with the increasing star formation rate of the galaxies. At similar stellar masses the blue galaxies exhibit a significantly higher outflow equivalent width as compared to red galaxies. The outflow equivalent width shows strong effect with star formation surface density ({\Sigma}SFR) of the sample. For the disk galaxies, the outflow equivalent width is higher for the face-on systems as compared to the edge-on ones, indicating that for the disk galaxies, the outflowing gas is primarily bipolar in geometry. Galaxies typically exhibit outflow velocities ranging from -200 km/s to -300 km/s and on average the face-on galaxies exhibit higher outflow velocity as compared to the edge-on ones. Galaxies with irregular morphologies exhibit outflow equivalent width as well as outflow velocities comparable to face on disk galaxies. These galaxies exhibit minimum mass outflow rates > 5-7 Msun/yr and a mass loading factor ({\eta} = dMout/dt /SFR) comparable to the star formation rates of the galaxies.

We explore the role of environment in the evolution of galaxies over 0.1<z<0.7 using the final zCOSMOS-bright data set. Using the red fraction of galaxies as a proxy for the quenched population, we find that the fraction of red galaxies increases with the environmental overdensity and with the stellar mass, consistent with previous works. As at lower redshift, the red fraction appears to be separable in mass and environment, suggesting the action of two processes: mass and environmental quenching. The parameters describing these appear to be essentially the same at z~0.7 as locally. We explore the relation between red fraction, mass and environment also for the central and satellite galaxies separately, paying close attention to the effects of impurities in the central-satellite classification and using carefully constructed samples matched in stellar mass. There is little evidence for a dependence of the red fraction of centrals on overdensity. Satellites are consistently redder at all overdensities, and the satellite quenching efficiency increases with overdensity at 0.1<z<0.4. This is less marked at higher redshift, but both are nevertheless consistent with the equivalent local measurements. At a given stellar mass, the fraction of galaxies that are satellites also increases with the overdensity. At a given overdensity and mass, the obtained relation between the environmental quenching and the satellite fraction agrees well with the satellite quenching efficiency, demonstrating that the environmental quenching in the overall population is consistent with being entirely produced through the satellite quenching process at least up to z=0.7. However, despite the unprecedented size of our high redshift samples, the associated statistical uncertainties are still significant and our statements should be understood as approximations to physical reality, rather than physically exact formulae.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

A sample of 94 narrow line AGN with 0.65<z<1.20 has been selected from the 20k-Bright zCOSMOS galaxy sample by detection of the high-ionization [NeV]3426 line. Taking advantage of the large amount of data available in the COSMOS field, the properties of the [NeV]-selected Type-2 AGN have been investigated, focusing on their host galaxies, X-ray emission, and optical line flux ratios. Finally, the diagnostic developed by Gilli et al. (2010), based on the X-ray to [NeV] luminosity ratio, has been exploited to search for the more heavily obscured AGN. We found that [Ne v]-selected narrow line AGN have Seyfert 2-like optical spectra, although with emission line ratios diluted by a star-forming component. The ACS morphologies and stellar component in the optical spectra indicate a preference for our Type-2 AGN to be hosted in early-spirals with stellar masses greater than 10^(9.5-10)Msun, on average higher than those of the galaxy parent sample. The fraction of galaxies hosting [NeV]-selected obscured AGN increases with the stellar mass, reaching a maximum of about 3% at 2x10^11 Msun. A comparison with other selection techniques at z~1 shows that the detection of the [Ne v] line is an effective method to select AGN in the optical band, in particular the most heavily obscured ones, but can not provide by itself a complete census of AGN2. Finally, the high fraction of [NeV]-selected Type-2 AGN not detected in medium-deep Chandra observations (67%) is suggestive of the inclusion of Compton-thick sources in our sample. The presence of a population of heavily obscured AGN is corroborated by the X-ray to [NeV] ratio; we estimated, by mean of X-ray stacking technique and simulations, that the Compton-thick fraction in our sample of Type-2 AGN is 43+-4%, in good agreement with standard assumptions by the XRB synthesis models.

We examine the red fraction of central and satellite galaxies in the large zCOSMOS group catalog out to z ~ 0.8 correcting for both the incompleteness in stellar mass and for the less than perfect purities of the central and satellite samples. We show that, at all masses and at all redshifts, the fraction of satellite galaxies that have been quenched, i.e., are red, is systematically higher than that of centrals, as seen locally in the Sloan Digital Sky Survey (SDSS). The satellite quenching efficiency, which is the probability that a satellite is quenched because it is a satellite rather than a central, is, as locally, independent of stellar mass. Furthermore, the average value is about 0.5, which is also very similar to that seen in the SDSS. We also construct the mass functions of blue and red centrals and satellites and show that these broadly follow the predictions of the Peng et al. analysis of the SDSS groups. Together, these results indicate that the effect of the group environment in quenching satellite galaxies was very similar when the universe was about a half its present age, as it is today.

I present a generalized power-law diagnostic that allows to identify the presence of active galactic nuclei (AGN) in infrared (IR) galaxies at z>1, down to flux densities at which the extragalactic IR background is mostly resolved. I derive this diagnostic from the analysis of 174 galaxies with Snu(24)>80 microJy and spectroscopic redshifts zspec>1 in the Chandra Deep Field South, for which I study the rest-frame UV/optical/near-IR spectral energy distributions (SEDs), after subtracting a hot-dust, power-law component with three possible spectral indices alpha=1.3, 2.0 and 3.0. I obtain that 35% of these 24micron sources are power-law composite galaxies (PLCGs), which I define as those galaxies for which the SED fitting with stellar templates, without any previous power-law subtraction, can be rejected with >2sigma confidence. Subtracting the power-law component from the PLCG SEDs produces stellar-mass correction factors <1.5 in >80% of cases. The PLCG incidence is especially high (47%) at z=1.0-1.5. To unveil which PLCGs host AGN, I conduct a combined analysis of 4Ms X-ray data, galaxy morphologies, and a greybody modelling of the hot dust. I find that: 1) 77% of all the X-ray AGN in my 24micron sample at z=1.0-1.5 are recognised by the PLCG criterion; 2) PLCGs with alpha=1.3 or 2.0 have regular morphologies and T_dust >~1000 K, indicating nuclear activity. Instead, PLCGs with alpha=3.0 are characterised by disturbed galaxy dynamics, and a hot interstellar medium can explain their dust temperatures T_dust ~700-800 K. Overall, my results indicate that the fraction of AGN among 24 micron sources is between ~30% and 52% at z=1.0-1.5.

We investigate the link between AGN activity, star-formation and stellar mass of the host galaxy at 0<z<1, looking for spectroscopic traces of AGN and aging of the host. This work provides an extension of the existing studies at z<0.1 and contributes to shed light on galaxy evolution at intermediate redshifts. We used the zCOSMOS 20k data to create a sample of galaxies at z<1. We divided the sample in several mass-redshift bins to obtain stacked galaxy spectra with an improved S/N. We exploited emission-line diagnostic diagrams to separate AGN from star-forming galaxies. We found indication of a role for the total galaxy stellar mass in leading galaxy classification. Stacked spectra show AGN signatures above the log M_*/M_sun>10.2 threshold. Moreover, the stellar populations of AGN hosts are found to be older with respect to star-forming and composites galaxies. This could be due to the the tendency of AGN to reside in massive hosts. The dependence of the AGN classification on the stellar mass is in agreement with what has been already found in previous studies. It is consistent with, together with the evidence of older stellar populations inhabiting the AGN-like galaxies, the downsizing scenario. In particular, our evidence points to an evolutionary scenario where the AGN-feedback is capable of quenching the star formation in the most massive galaxies. Therefore, the AGN-feedback is the best candidate for initiating the passive evolutionary phase of galaxies.

In this paper we describe the first data release of the UltraVISTA near-infrared imaging survey of the COSMOS field. We summarise the key goals and design of the survey and provide a detailed description of our data reduction techniques . We provide stacked, sky-subtracted images in $YJHK_{\rm s}$ and narrow-band filters constructed from data collected during the first year of UltraVISTA observations. Our stacked images reach $5\sigma$ $AB$ depths in an aperture of $2\arcsec$ diameter of $\sim 25$ in $Y$ and $\sim 24$ in $JHK_{\rm s}$ bands and all have sub-arcsecond seeing. To this $5\sigma$ limit, our $K_{\rm s}$ catalogue contains 216,268 sources. We carry out a series of quality assessment tests on our images and catalogues, comparing our stacks with existing catalogues. The $1\sigma$ astrometric RMS in both directions for stars selected with $17.0<K_{\rm s}\rm {(AB)} <19.5$ is $\sim 0.08\arcsec$ in comparison to the publicly-available COSMOS ACS catalogues. Our images are resampled to the same pixel scale and tangent point as the publicly available COSMOS data and so may be easily used to generate multi-colour catalogues using this data. All images and catalogues presented in this paper are publicly available through ESO's "phase 3" archiving and distribution system and from the UltraVISTA web site.

We present new improved constraints on the Hubble parameter H(z) in the redshift range 0.15 < z < 1.1, obtained from the differential spectroscopic evolution of early-type galaxies as a function of redshift. We extract a large sample of early-type galaxies (\sim11000) from several spectroscopic surveys, spanning almost 8 billion years of cosmic lookback time (0.15 < z < 1.42). We select the most massive, red elliptical galaxies, passively evolving and without signature of ongoing star formation. Those galaxies can be used as standard cosmic chronometers, as firstly proposed by Jimenez & Loeb (2002), whose differential age evolution as a function of cosmic time directly probes H(z). We analyze the 4000 {\AA} break (D4000) as a function of redshift, use stellar population synthesis models to theoretically calibrate the dependence of the differential age evolution on the differential D4000, and estimate the Hubble parameter taking into account both statistical and systematical errors. We provide 8 new measurements of H(z) (see Tab. 4), and determine its change in H(z) to a precision of 5-12% mapping homogeneously the redshift range up to z \sim 1.1; for the first time, we place a constraint on H(z) at z \neq 0 with a precision comparable with the one achieved for the Hubble constant (about 5-6% at z \sim 0.2), and covered a redshift range (0.5 < z < 0.8) which is crucial to distinguish many different quintessence cosmologies. These measurements have been tested to best match a \Lambda CDM model, clearly providing a statistically robust indication that the Universe is undergoing an accelerated expansion. This method shows the potentiality to open a new avenue in constrain a variety of alternative cosmologies, especially when future surveys (e.g. Euclid) will open the possibility to extend it up to z \sim 2.

The chemical evolution of galaxies on a cosmological timescale is still a matter of debate despite the increasing number of available data provided by spectroscopic surveys of star-forming galaxies at different redshifts. The fundamental relations involving metallicity, such as the mass-metallicity relation (MZR) or the fundamental-metallicity relation, give controversial results about the reality of evolution of the chemical content of galaxies at a given stellar mass. In this work we shed some light on this issue using the completeness reached by the 20k bright sample of the zCOSMOS survey and using for the first time the nitrogen-to-oxygen ratio (N/O) as a tracer of the gas phase chemical evolution of galaxies that is independent of the star formation rate. Emission-line galaxies both in the SDSS and 20k zCOSMOS bright survey were used to study the evolution from the local Universe of the $MZR up to a redshift of 1.32 and the relation between stellar mass and nitrogen-to-oxygen ratio (MNOR) up to a redshift of 0.42 using the N2S2 parameter. All the physical properties derived from stellar continuum and gas emission-lines, including stellar mass, star formation rates, metallicity and N/O, were calculated in a self-consistent way over the full redshift range. We confirm the trend to find lower metallicities in galaxies of a given stellar mass in a younger Universe. This trend is even observed when taking possible selection effects into account that are due to the observed larger median star formation rates for galaxies at higher redshifts. We also find a significant evolution of the MNOR up to z = 0.4. Taking the slope of the O/H vs. N/O relation into account for the secondary-nitrogen production regime, the observed evolution of the MNOR is consistent with the trends found for both the MZR and its equivalent relation using new expressions to reduce its dependence on SFR.

We report the results of a comprehensive study of the relationship between galaxy size, stellar mass and specific star-formation rate (sSFR) at redshifts 1.3<z<1.5. Based on a mass complete (M_star >= 6x10^10 Msun), spectroscopic sample from the UKIDSS Ultra-deep Survey (UDS), with accurate stellar-mass measurements derived from spectro photometric fitting, we find that at z~1.4 the location of massive galaxies on the size-mass plane is determined primarily by their sSFR. At this epoch we find that massive galaxies which are passive (sSFR <= 0.1 Gyr^-1) follow a tight size-mass relation, with half-light radii a factor f=2.4+/-0.2 smaller than their local counterparts. Moreover, amongst the passive sub-sample we find no evidence that the off-set from the local size-mass relation is a function of stellar population age. Based on a sub-sample with dynamical mass estimates we also derive an independent estimate of f=2.3+/-0.3 for the typical growth in half-light radius between z~1.4 and the present day. Focusing on the passive sub-sample, we conclude that to produce the necessary evolution predominantly via major mergers would require an unfeasible number of merger events and over populate the high-mass end of the local stellar mass function. In contrast, we find that a scenario in which mass accretion is dominated by minor mergers can produce the necessary evolution, whereby an increase in stellar mass by a factor of ~2, accompanied by an increase in size by a factor of ~3.5, is sufficient to reconcile the size-mass relation at z~1.4 with that observed locally. Finally, we note that a significant fraction (44+/-12%) of the passive galaxies in our sample have a disk-like morphology, providing additional evidence that separate physical processes are responsible for the quenching of star-formation and the morphological transformation of massive galaxies (abridged).

In this paper we measure the merger fraction and rate, both minor and major, of massive early-type galaxies (M_star >= 10^11 M_Sun) in the COSMOS field, and study their role in mass and size evolution. We use the 30-band photometric catalogue in COSMOS, complemented with the spectroscopy of the zCOSMOS survey, to define close pairs with a separation 10h^-1 kpc <= r_p <= 30h-1 kpc and a relative velocity Delta v <= 500 km s^-1. We measure both major (stellar mass ratio mu = M_star,2/M_star,1 >= 1/4) and minor (1/10 <= mu < 1/4) merger fractions of massive galaxies, and study their dependence on redshift and on morphology. The merger fraction and rate of massive galaxies evolves as a power-law (1+z)^n, with major mergers increasing with redshift, n_MM = 1.4, and minor mergers showing little evolution, n_mm ~ 0. When split by their morphology, the minor merger fraction for early types is higher by a factor of three than that for spirals, and both are nearly constant with redshift. Our results show that massive early-type galaxies have undergone 0.89 mergers (0.43 major and 0.46 minor) since z ~ 1, leading to a mass growth of ~30%. We find that mu >= 1/10 mergers can explain ~55% of the observed size evolution of these galaxies since z ~ 1. Another ~20% is due to the progenitor bias (younger galaxies are more extended) and we estimate that very minor mergers (mu < 1/10) could contribute with an extra ~20%. The remaining ~5% should come from other processes (e.g., adiabatic expansion or observational effects). This picture also reproduces the mass growth and velocity dispersion evolution of these galaxies. We conclude from these results that merging is the main contributor to the size evolution of massive ETGs at z <= 1, accounting for ~50-75% of that evolution in the last 8 Gyr. Nearly half of the evolution due to mergers is related to minor (mu < 1/4) events.

We report the discovery of a massive ultra-compact quiescent galaxy that has been strongly-lensed into multiple images by a foreground galaxy at z = 0.960. This system was serendipitously discovered as a set of extremely Ks-bright high-redshift galaxies with red J - Ks colors using new data from the UltraVISTA YJHKs near-infrared survey. The system was also previously identified as an optically-faint lens/source system using the COSMOS ACS imaging by Faure et al. (2008, 2011). Photometric redshifts for the three brightest images of the source galaxy determined from twenty-seven band photometry place the source at z = 2.4 +/- 0.1. We provide an updated lens model for the system which is a good fit to the positions and morphologies of the galaxies in the ACS image. The lens model implies that the magnification of the three brightest images is a factor of 4 - 5. We use the lens model, combined with the Ks-band image to constrain the size and Sersic profile of the galaxy. The best-fit model is an ultra-compact galaxy (Re = 0.64^{+0.08}_{-0.18} kpc, lensing-corrected), with a Sersic profile that is intermediate between a disk and bulge profile (n = 2.2^{+2.3}_{-0.9}). We present aperture photometry for the source galaxy images which have been corrected for flux contamination from the central lens. The best-fit stellar population model is a massive galaxy (Log(M_{star}/M_{sol}) = 10.8^{+0.1}_{-0.1}, lensing-corrected) with an age of 1.0^{+1.0}_{-0.4} Gyr, moderate dust extinction (Av = 0.8^{+0.5}_{-0.6}), and a low specific star formation rate (Log(SSFR) < -11.0 yr^{-1}). This is typical of massive "red-and-dead" galaxies at this redshift and confirms that this source is the first bona fide strongly-lensed massive ultra-compact quiescent galaxy to be discovered. We conclude with a discussion of the prospects of finding a larger sample of these galaxies.

We present an optical group catalog between 0.1 < z < 1 based on 16,500 high-quality spectroscopic redshifts in the completed zCOSMOS-bright survey. The catalog published herein contains 1498 groups in total and 192 groups with more than five observed members. The catalog includes both group properties and the identification of the member galaxies. Based on mock catalogs, the completeness and purity of groups with three and more members should be both about 83% with respect to all groups that should have been detectable within the survey, and more than 75% of the groups should exhibit a one-to-one correspondence to the "real" groups. Particularly at high redshift, there are apparently more galaxies in groups in the COSMOS field than expected from mock catalogs. We detect clear evidence for the growth of cosmic structure over the last seven billion years in the sense that the fraction of galaxies that are found in groups (in volume-limited samples) increases significantly with cosmic time. In the second part of the paper, we develop a method for associating galaxies that only have photo-z to our spectroscopically identified groups. We show that this leads to improved definition of group centers, improved identification of the most massive galaxies in the groups, and improved identification of central and satellite galaxies, where we define the former to be galaxies at the minimum of the gravitational potential wells. Subsamples of centrals and satellites in the groups can be defined with purities up to 80%, while a straight binary classification of all group and non-group galaxies into centrals and satellites achieves purities of 85% and 75%, respectively, for the spectroscopic sample.

We present a group-galaxy cross-correlation analysis using a group catalog produced from the 16,500 spectra from the optical zCOSMOS galaxy survey. Our aim is to perform a consistency test in the redshift range 0.2 < z < 0.8 between the clustering strength of the groups and mass estimates that are based on the richness of the groups. We measure the linear bias of the groups by means of a group-galaxy cross-correlation analysis and convert it into mass using the bias-mass relation for a given cosmology, checking the systematic errors using realistic group and galaxy mock catalogs. The measured bias for the zCOSMOS groups increases with group richness as expected by the theory of cosmic structure formation and yields masses that are reasonably consistent with the masses estimated from the richness directly, considering the scatter that is obtained from the 24 mock catalogs. An exception are the richest groups at high redshift (estimated to be more massive than 10^13.5 M_sun), for which the measured bias is significantly larger than for any of the 24 mock catalogs (corresponding to a 3-sigma effect), which is attributed to the extremely large structure that is present in the COSMOS field at z ~ 0.7. Our results are in general agreement with previous studies that reported unusually strong clustering in the COSMOS field.

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

The AzTEC 1.1 mm survey of the SCUBA HAlf Degree Extragalactic Survey (SHADES) fields is the largest (0.7 deg2) blank-field millimetre-wavelength survey undertaken to date at a resolution of ~18" and a depth of ~1 mJy. We have used the deep optical-to-radio multi-wavelength data in the SHADES Lockman Hole East and SXDF/UDS fields to obtain galaxy identifications for ~64% (~80% with tentative identifications) of the 148 AzTEC-SHADES 1.1 mm sources, exploiting deep radio and 24 um data complemented by methods based on 8 um flux-density and red optical-infrared (i-K) colour. This unusually high identification rate can be attributed to the relatively bright millimetre-wavelength flux-density threshold, combined with the relatively deep supporting multi-frequency data. We have further exploited the optical-mid-infrared-radio data to derive a ~60% (~75% with tentative identifications) complete redshift distribution for the AzTEC-SHADES sources, yielding a median redshift of z~2.2, with a high-redshift tail extending to at least z~4. Despite the larger area probed by the AzTEC survey relative to the original SCUBA SHADES imaging, the redshift distribution of the AzTEC sources is consistent with that displayed by the SCUBA sources, and reinforces tentative evidence that the redshift distribution of mm/sub-mm sources in the Lockman Hole field is significantly different from that found in the SXDF/UDS field. Comparison with simulated surveys of similar scale extracted from semi-analytic models based on the Millennium simulation indicates that this is as expected if the mm/sub-mm sources are massive (M>10^11 Mo) star-forming galaxies tracing large-scale structures over scales of 10-20 Mpc. This confirms the importance of surveys covering several square degrees (as now underway with SCUBA2) to obtain representative samples of bright (sub)mm-selected galaxies.

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

The COSMOS field has been the subject of a wide range of observations, with a number of studies focusing on reconstructing the 3D dark matter density field. Typically, these studies have focused on one given method or tracer. In this paper, we reconstruct the distribution of mass in the COSMOS field out to a redshift z=1 by combining Hubble Space Telescope weak lensing measurements with zCOSMOS spectroscopic measurements of galaxy clustering. The distribution of galaxies traces the distribution of mass with high resolution (particularly in redshift, which is not possible with lensing), and the lensing data empirically calibrates the mass normalisation (bypassing the need for theoretical models). Two steps are needed to convert a galaxy survey into a density field. The first step is to create a smooth field from the galaxy positions, which is a point field. We investigate four possible methods for this: (i) Gaussian smoothing, (ii) convolution with truncated isothermal sphere, (iii) fifth nearest neighbour smoothing and (iv) a muliti-scale entropy method. The second step is to rescale this density field using a bias prescription. We calculate the optimal bias scaling for each method by comparing predictions from the smoothed density field with the measured weak lensing data, on a galaxy-by-galaxy basis. In general, we find scale-independent bias for all the smoothing schemes, to a precision of 10%. For the nearest neighbour smoothing case, we find the bias to be 2.51\pm 0.25. We also find evidence for a strongly evolving bias, increasing by a factor of ~3.5 between redshifts 0<z<0.8. We believe this strong evolution can be explained by the fact that we use a flux limited sample to build the density field.

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

One key piece of information missing from high redshift galaxy surveys is the galaxies' cold gas contents. We present a new method to indirectly determine cold gas surface densities and integrated gas masses from galaxy star formation rates and to separate the atomic and molecular gas components. Our predicted molecular and total gas surface densities and integrated masses are in very good agreement with direct measurements quoted in the literature for low and high-z galaxies. We apply this method to predict the gas content for a sample of $\sim 57000$ galaxies in the COSMOS field at $0.5 \leq z \leq 2.0$, selected to have $I_{AB} < 24$ mag. This approach allows us to investigate in detail the redshift evolution of galaxy cold and molecular gas content versus stellar mass and to provide fitting formulae for galaxy gas fractions. We find a clear trend between galaxy gas fraction, molecular gas fraction and stellar mass with redshift, suggesting that massive galaxies consume and/or expel their gas at higher redshift than less massive objects and have lower fractions of their gas in molecular form. The characteristic stellar mass separating gas- from stellar-dominated galaxies decreases with time. This indicates that massive galaxies reach a gas-poor state earlier than less massive objects. These trends can be considered to be another manifestation of downsizing in star formation activity.

Using the group catalog obtained from zCOSMOS spectroscopic data and the complementary photometric data from the COSMOS survey, we explore segregation effects occurring in groups of galaxies at intermediate/high redshifts. We built two composite groups at intermediate (0.2 <= z <= 0.45) and high (0.45 < z <= 0.8) redshifts, and we divided the corresponding composite group galaxies into three samples according to their distance from the group center. We explored how galaxy stellar masses and colors - working in narrow bins of stellar masses - vary as a function of the galaxy distance from the group center. We found that the most massive galaxies in our sample (Log(M_gal/M_sun) >= 10.6) do not display any strong group-centric dependence of the fractions of red/blue objects. For galaxies of lower masses (9.8 <= Log(M_gal/M_sun) <= 10.6) there is a radial dependence in the changing mix of red and blue galaxies. This dependence is most evident in poor groups, whereas richer groups do not display any obvious trend of the blue fraction. Interestingly, mass segregation shows the opposite behavior: it is visible only in rich groups, while poorer groups have a a constant mix of galaxy stellar masses as a function of radius. We suggest a simple scenario where color- and mass-segregation originate from different physical processes. While dynamical friction is the obvious cause for establishing mass segregation, both starvation and galaxy-galaxy collisions are plausible mechanisms to quench star formation in groups at a faster rate than in the field. In poorer groups the environmental effects are caught in action superimposed to secular galaxy evolution. Their member galaxies display increasing blue fractions when moving from the group center to more external regions, presumably reflecting the recent accretion history of these groups.

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

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.

We present a photometric and spectroscopic study of galaxies at 0.5<z<1 as a function of environment based on data from the zCOSMOS survey. There is a fair amount of evidence that galaxy properties depend on mass of groups and clusters, in the sense that quiescent galaxies prefer more massive systems. We base our analysis on a mass-selected environment using X-ray groups of galaxies and define the group membership using a large number of spectroscopic redshifts from zCOSMOS. We show that the fraction of red galaxies is higher in groups than in the field at all redshifts probed in our study. Interestingly, the fraction of [OII] emitters on the red sequence increases at higher redshifts in groups, while the fraction does not strongly evolve in the field. This is due to increased dusty star formation activities and/or increased activities of active galactic nuclei (AGNs) in high redshift groups. We study these possibilities using the 30-band photometry and X-ray data. We find that the stellar population of the red [OII] emitters in groups is old and there is no clear hint of dusty star formation activities in those galaxies. The observed increase of red [OII] emitters in groups is likely due to increased AGN activities. However, our overall statistics is poor and any firm conclusions need to be drawn from a larger statistical sample of z~1 groups.

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 explore the connection between black hole growth at the center of obscured quasars selected from the XMM-COSMOS survey and the physical properties of their host galaxies. We study a bolometric regime (<Lbol > 8 x 10^45 erg/s) where several theoretical models invoke major galaxy mergers as the main fueling channel for black hole accretion. We confirm that obscured quasars mainly reside in massive galaxies (Mstar>10^10 Msun) and that the fraction of galaxies hosting such powerful quasars monotonically increases with the stellar mass. We stress the limitation of the use of rest-frame color-magnitude diagrams as a diagnostic tool for studying galaxy evolution and inferring the influence that AGN activity can have on such a process. We instead use the correlation between star-formation rate and stellar mass found for star-forming galaxies to discuss the physical properties of the hosts. We find that at z ~1, ~62% of Type-2 QSOs hosts are actively forming stars and that their rates are comparable to those measured for normal star-forming galaxies. The fraction of star-forming hosts increases with redshift: ~71% at z ~2, and 100% at z ~3. We also find that the the evolution from z ~1 to z ~3 of the specific star-formation rate of the Type-2 QSO hosts is in excellent agreement with that measured for star-forming galaxies. From the morphological analysis, we conclude that most of the objects are bulge-dominated galaxies, and that only a few of them exhibit signs of recent mergers or disks. Finally, bulge-dominated galaxies tend to host Type-2 QSOs with low Eddington ratios (lambda<0.1), while disk-dominated or merging galaxies have at their centers BHs accreting at high Eddington ratios (lambda > 0.1).

Close encounters between galaxies are expected to be a viable mechanism, as predicted by numerical simulations, by which accretion onto supermassive black holes can be initiated. To test this scenario, we construct a sample of 562 galaxies (M*>2.5x10^10 M_sun) in kinematic pairs over the redshift range 0.25 < z < 1.05 that are more likely to be interacting than a well-matched control sample of 2726 galaxies not identified as being in a pair, both from the zCOSMOS 20k spectroscopic catalog. Galaxies that harbor an active galactic nucleus (AGN) are identified on the basis of their X-ray emission (L_x>2x10^42 erg s^-1) detected by Chandra. We find a higher fraction of AGN in galaxies in pairs relative to isolated galaxies of similar stellar mass. Our result is primarily due to an enhancement of AGN activity, by a factor of 1.9 (observed) and 2.6 (intrinsic), for galaxies in pairs of projected separation less than 75 kpc and line-of-sight velocity offset less than 500 km s^-1. This study demonstrates that close kinematic pairs are conducive environments for black hole growth either indicating a causal physical connection or an inherent relation, such as, to enhanced star formation. In the Appendix, we describe a method to estimate the intrinsic fractions of galaxies (either in pairs or the field) hosting an AGN with confidence intervals, and an excess fraction in pairs. We estimate that 17.8_{-7.4}^{+8.4}% of all moderate-luminosity AGN activity takes place within galaxies undergoing early stages of interaction that leaves open the question as to what physical processes are responsible for fueling the remaining ~80% that may include late-stage mergers.

With this paper, we release accurate photometric redshifts for 1692 counterparts to Chandra sources in the central square degree of the COSMOS field. The availability of a large training set of spectroscopic redshifts that extends to faint magnitudes enabled photometric redshifts comparable to the highest quality results presently available for normal galaxies. We demonstrate that morphologically extended, faint X-ray sources without optical variability are more accurately described by a library of normal galaxies (corrected for emission lines) than by AGN-dominated templates, even if these sources have AGN-like X-ray luminosities. Preselecting the library on the bases of the source properties allowed us to reach an accuracy sigma_(Delta z/(1+z_spec)) \sim0.015 with a fraction of outliers of 5.8% for the entire Chandra-COSMOS sample. In addition, we release revised photometric redshifts for the 1735 optical counterparts of the XMM-detected sources over the entire 2 sq. deg.of COSMOS. For 248 sources, our updated photometric redshift differs from the previous release by Delta z>0.2. These changes are predominantly due to the inclusion of newly available deep H-band photometry H_AB=24 mag. We illustrate once again the importance of a spectroscopic training sample and how an assumption about the nature of a source together with the number and the depth of the available bands influence the accuracy of the photometric redshifts determined for AGN. These considerations should be kept in mind when defining the observational strategies of upcoming large surveys targeting AGN, such as eROSITA at X-ray energies and ASKAP/EMU in the radio band.

We map the radial and azimuthal distribution of Mg II gas within 200 kpc (physical) of 4000 galaxies at redshifts 0.5 < z < 0.9 using co-added spectra of more than 5000 background galaxies at z > 1. We investigate the variation of Mg II rest frame equivalent width as a function of the radial impact parameter for different subsets of foreground galaxies selected in terms of their rest-frame colors and masses. Blue galaxies have a significantly higher average Mg II equivalent width at close galactocentric radii as compared to the red galaxies. Amongst the blue galaxies, there is a correlation between Mg II equivalent width and galactic stellar mass of the host galaxy. We also find that the distribution of Mg II absorption around group galaxies is more extended than that for non-group galaxies, and that groups as a whole have more extended radial profiles than individual galaxies. Interestingly, these effects can be satisfactorily modeled by a simple superposition of the absorption profiles of individual member galaxies, assuming that these are the same as those of non-group galaxies, suggesting that the group environment may not significantly enhance or diminish the Mg II absorption of individual galaxies. We show that there is a strong azimuthal dependence of the Mg II absorption within 50 kpc of inclined disk-dominated galaxies, indicating the presence of a strongly bipolar outflow aligned along the disk rotation axis. There is no significant dependence of Mg II absorption on the apparent inclination angle of disk-dominated galaxies.

We have studied the evolution of high redshift quiescent galaxies over an effective area of ~1.7 deg^2 in the COSMOS field. Galaxies have been divided according to their star-formation activity and the evolution of the different populations has been investigated in detail. We have studied an IRAC (mag_3.6 < 22.0) selected sample of ~18000 galaxies at z > 1.4 with multi-wavelength coverage. We have derived accurate photometric redshifts (sigma=0.06) and other important physical parameters through a SED-fitting procedure. We have divided our sample into actively star-forming, intermediate and quiescent galaxies depending on their specific star formation rate. We have computed the galaxy stellar mass function of the total sample and the different populations at z=1.4-3.0. We have studied the properties of high redshift quiescent galaxies finding that they are old (1-4 Gyr), massive (log(M/M_sun)~10.65), weakly star forming stellar populations with low dust extinction (E(B-V) < 0.15) and small e-folding time scales (tau ~ 0.1-0.3 Gyr). We observe a significant evolution of the quiescent stellar mass function from 2.5 < z < 3.0 to 1.4 < z < 1.6, increasing by ~ 1 dex in this redshift interval. We find that z ~ 1.5 is an epoch of transition of the GSMF. The fraction of star-forming galaxies decreases from 60% to 20% from z ~ 2.5-3.0 to z ~ 1.4-1.6 for log(M/M_sun) > 11, while the quiescent population increases from 10% to 50% at the same redshift and mass intervals. We compare the fraction of quiescent galaxies derived with that predicted by theoretical models and find that the Kitzbichler & White (2007) model is the one that better reproduces the data. Finally, we calculate the stellar mass density of the star-forming and quiescent populations finding that there is already a significant number of quiescent galaxies at z > 2.5 (rho~6.0 MsunMpc^-3).

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 have analysed a sample of 1292 4.5 micron-selected galaxies at z>=3, over 0.6 square degrees of the UKIRT Infrared Deep Survey (UKIDSS) Ultra Deep Survey (UDS). Using photometry from the U band through 4.5 microns, we have obtained photometric redshifts and derived stellar masses for our sources. Only two of our galaxies potentially lie at z>5. We have studied the galaxy stellar mass function at 3<=z<5, based on the 1213 galaxies in our catalogue with [4.5]<= 24.0. We find that: i) the number density of M > 10^11 Msun galaxies increased by a factor > 10 between z=5 and 3, indicating that the assembly rate of these galaxies proceeded > 20 times faster at these redshifts than at 0<z<2; ii) the Schechter function slope alpha is significantly steeper than that displayed by the local stellar mass function, which is both a consequence of the steeper faint end and the absence of a pure exponential decline at the high-mass end; iii) the evolution of the comoving stellar mass density from z=0 to 5 can be modelled as log10 (rho_M) =-(0.05 +/- 0.09) z^2 - (0.22 -/+ 0.32) z + 8.69. At 3<=z<4, more than 30% of the M > 10^11 Msun galaxies would be missed by optical surveys with R<27 or z<26. Thus, our study demonstrates the importance of deep mid-IR surveys over large areas to perform a complete census of massive galaxies at high z and trace the early stages of massive galaxy assembly.

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.

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.

The Cosmic Far-Infrared Background (CIB) at wavelengths around 160 {\mu}m corresponds to the peak intensity of the whole Extragalactic Background Light, which is being measured with increasing accuracy. However, the build up of the CIB emission as a function of redshift, is still not well known. Our goal is to measure the CIB history at 70 {\mu}m and 160 {\mu}m at different redshifts, and provide constraints for infrared galaxy evolution models. We use complete deep Spitzer 24 {\mu}m catalogs down to about 80 {\mu}Jy, with spectroscopic and photometric redshifts identifications, from the GOODS and COSMOS deep infrared surveys covering 2 square degrees total. After cleaning the Spitzer/MIPS 70 {\mu}m and 160 {\mu}m maps from detected sources, we stacked the far-IR images at the positions of the 24 {\mu}m sources in different redshift bins. We measured the contribution of each stacked source to the total 70 and 160 {\mu}m light, and compare with model predictions and recent far-IR measurements made with Herschel/PACS on smaller fields. We have detected components of the 70 and 160 {\mu}m backgrounds in different redshift bins up to z ~ 2. The contribution to the CIB is maximum at 0.3 <= z <= 0.9 at 160{\mu}m (and z <= 0.5 at 70 {\mu}m). A total of 81% (74%) of the 70 (160) {\mu}m background was emitted at z < 1. We estimate that the AGN relative contribution to the far-IR CIB is less than about 10% at z < 1.5. We provide a comprehensive view of the CIB buildup at 24, 70, 100, 160 {\mu}m. IR galaxy models predicting a major contribution to the CIB at z < 1 are in agreement with our measurements, while our results discard other models that predict a peak of the background at higher redshifts. Our results are available online http://www.ias.u-psud.fr/irgalaxies/ .

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

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for ~50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (~95%) have IRAC counterparts, and 1394 (~78%) have MIPS 24micron detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ~150 high redshift (z>1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L_X>10^{44} erg s^{-1}) X-ray luminosity is ~15-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

Using Chandra observations we have identified a sample of seven off-nuclear X-ray sources, in the redshift range z=0.072-0.283, located within optically bright galaxies in the COSMOS Survey. Using the multi-wavelength coverage available in the COSMOS field, we study the properties of the host galaxies of these ULXs. In detail, we derived their star formation rate from H_alpha measurements and their stellar masses using SED fitting techniques with the aim to compute the probability to have an off-nuclear source based on the host galaxy properties. We divide the host galaxies in different morphological classes using the available ACS/HST imaging. We find that our ULXs candidates are located in regions of the SFR versus M$_star$ plane where one or more off-nuclear detectable sources are expected. From a morphological analysis of the ACS imaging and the use of rest-frame colours, we find that our ULXs are hosted both in late and early type galaxies. Finally, we find that the fraction of galaxies hosting a ULX ranges from ~0.5% to ~0.2% going from L[0.5-2 keV]=3 x 10^39 erg s^-1 to L[0.5-2 keV]= 2 x 10^40 erg s^-1.

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 present a large robust sample of 1503 reliable and unconfused 70microm selected sources from the multiwavelength data set of the Cosmic Evolution Survey (COSMOS). Using the Spitzer IRAC and MIPS photometry, we estimate the total infrared luminosity, L_IR (8--1000 microns), by finding the best fit template from several different template libraries. The long wavelength 70 and 160 micron data allow us to obtain a reliable estimate of L_IR, accurate to within 0.2 and 0.05 dex, respectively. The 70 micron data point enables a significant improvement over the luminosity estimates possible with only a 24 micron detection. The full sample spans a wide range in L_IR, L_IR ~ 10^8-10^14 L_sun, with a median luminosity of 10^11.4 L_sun. We identify a total of 687 luminous, 303 ultraluminous, and 31 hyperluminous infrared galaxies (LIRGs, ULIRGs, and HyLIRGs) over the redshift range 0.01<z<3.5 with a median redshift of 0.5. Presented here are the full spectral energy distributions for each of the sources compiled from the extensive multiwavelength data set from the ultraviolet (UV) to the far-infrared (FIR). Using SED fits we find possible evidence for a subset of cooler ultraluminous objects than observed locally. However, until direct observations at longer wavelengths are obtained, the peak of emission and the dust temperature cannot be well constrained. We use these SEDs, along with the deep radio and X-ray coverage of the field, to identify a large sample of candidate active galactic nuclei (AGN). We find that the fraction of AGN increases strongly with L_IR, as it does in the local universe, and that nearly 70% of ULIRGs and all HyLIRGs likely host a powerful AGN.

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.

(Abriged) We report on the measurement of the rest frame K-band luminosity and total stellar mass of the hosts of 89 broad line Active Galactic Nuclei detected in the zCOSMOS survey in the redshift range 1<z<2.2. The unprecedented multiwavelength coverage of the survey field allows us to disentangle the emission of the host galaxy from that of the nuclear black hole in their Spectral Energy Distributions. We derive an estimate of black hole masses through the analysis of the broad Mg II emission lines observed in the medium-resolution spectra taken with VIMOS/VLT as part of the zCOSMOS project. We found that, as compared to the local value, the average black hole to host galaxy mass ratio appears to evolve positively with redshift, with a best fit evolution of the form (1+z)^{0.68 \pm0.12 +0.6 -0.3}, where the large asymmetric systematic errors stem from the uncertainties in the choice of IMF, in the calibration of the virial relation used to estimate BH masses and in the mean QSO SED adopted. A thorough analysis of observational biases induced by intrinsic scatter in the scaling relations reinforces the conclusion that an evolution of the MBH-M* relation must ensue for actively growing black holes at early times: either its overall normalization, or its intrinsic scatter (or both) appear to increase with redshift. This can be interpreted as signature of either a more rapid growth of supermassive black holes at high redshift, a change of structural properties of AGN hosts at earlier times, or a significant mismatch between the typical growth times of nuclear black holes and host galaxies.

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.

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

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.

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.

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

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.

We present photometric redshifts and spectral energy distribution (SED) classifications for a sample of 1542 optically identified sources detected with XMM in the COSMOS field. Our template fitting classifies 46 sources as stars and 464 as non-active galaxies, while the remaining 1032 require templates with an AGN contribution. High accuracy in the derived photometric redshifts was accomplished as the result of 1) photometry in up to 30 bands with high significance detections, 2) a new set of SED templates including 18 hybrids covering the far-UV to mid-infrared, which have been constructed by the combination of AGN and non-active galaxies templates, and 3) multi-epoch observations that have been used to correct for variability (most important for type 1 AGN). The reliability of the photometric redshifts is evaluated using the sub-sample of 442 sources with measured spectroscopic redshifts. We achieved an accuracy of $\sigma_{\Delta z/(1+z_{spec})} = 0.014$ for i$_{AB}^*<$22.5 ($\sigma_{\Delta z/(1+z_{spec})} \sim0.015$ for i$_{AB}^*<$24.5). The high accuracies were accomplished for both type 2 (where the SED is often dominated by the host galaxy) and type 1 AGN and QSOs out to $z=4.5$. The number of outliers is a large improvement over previous photometric redshift estimates for X-ray selected sources (4.0% and 4.8% outliers for i$_{AB}^*<$22.5 and i$_{AB}^*<$24.5, respectively). We show that the intermediate band photometry is vital to achieving accurate photometric redshifts for AGN, whereas the broad SED coverage provided by mid infrared (Spitzer/IRAC) bands is important to reduce the number of outliers for normal galaxies.

We present accurate photometric redshifts in the 2-deg2 COSMOS field. The redshifts are computed with 30 broad, intermediate, and narrow bands covering the UV (GALEX), Visible-NIR (Subaru, CFHT, UKIRT and NOAO) and mid-IR (Spitzer/IRAC). A chi2 template-fitting method (Le Phare) was used and calibrated with large spectroscopic samples from VLT-VIMOS and Keck-DEIMOS. We develop and implement a new method which accounts for the contributions from emission lines (OII, Hbeta, Halpha and Ly) to the spectral energy distributions (SEDs). The treatment of emission lines improves the photo-z accuracy by a factor of 2.5. Comparison of the derived photo-z with 4148 spectroscopic redshifts (i.e. Delta z = zs - zp) indicates a dispersion of sigma_{Delta z/(1+zs)}=0.007 at i<22.5, a factor of 2-6 times more accurate than earlier photo-z in the COSMOS, CFHTLS and COMBO-17 survey fields. At fainter magnitudes i<24 and z<1.25, the accuracy is sigma_{Delta z/(1+zs)}=0.012. The deep NIR and IRAC coverage enables the photo-z to be extended to z~2 albeit with a lower accuracy (sigma_{Delta z/(1+zs)}=0.06 at i~24). The redshift distribution of large magnitude-selected samples is derived and the median redshift is found to range from z=0.66 at 22<i<22.5 to z=1.06 at 24.5<i<25. At i<26.0, the multi-wavelength COSMOS catalog includes approximately 607,617 objects. The COSMOS-30 photo-z enable the full exploitation of this survey for studies of galaxy and large scale structure evolution at high redshift.

We study zCOSMOS-bright optical spectra for 609 Spitzer/MIPS 24 micron-selected galaxies with S(24um)> 0.30 mJy and I<22.5 (AB mag) over 1.5 sq. deg. of the COSMOS field. From emission-line diagnostics we find that: 1) star-formation rates (SFR) derived from the observed Halpha and Hbeta lines underestimate, on average, the total SFR by factors ~5 and 10, respectively; 2) both the Calzetti et al. and the Milky Way reddening laws are suitable to describe the extinction observed in infrared (IR) sources in most cases; 3) some IR galaxies at z<0.3 have low abundances, but many others with similar IR luminosities and redshifts are chemically enriched; 4) The average [OIII]/Hbeta ratios of nuLnu(24um)>10^11 Lsun galaxies at 0.6<z<0.7 are ~0.6 dex higher than the average ratio of all zCOSMOS galaxies at similar redshifts. Massive star formation and active galactic nuclei (AGN) could simultaneously be present in those galaxies with the highest ionising fluxes; 5) ~1/3 of the galaxies with metallicity measurements at 0.5<z<0.7 lie below the general mass-metallicity relation at the corresponding redshifts. The strengths of the 4000 Angstrom break and the Hdelta EW of our galaxies show that secondary bursts of star formation are needed to explain the spectral properties of most IR sources. The LIRG and ULIRG phases occur, on average, between 10^7 and 10^8 years after the onset of a starburst on top of underlying older stellar populations. These results are valid for galaxies of different IR luminosities at 0.6<z<1.0 and seem independent of the mechanisms triggering star formation.

We present here some of the first results we have obtained on the study of the optical spectra of Spitzer/MIPS 24 micron-selected galaxies in the COSMOS field. This is part of a series of studies we are conducting to analyse the optical spectral properties of mid-infrared (mid-IR) galaxies with different IR luminosities up to high redshifts. The results shown here correspond to the brightest S(24 micron)>2 mJy IR galaxy population at z<1.

We present the observed correlations between rest-frame 8, 24, 70 and 160 um monochromatic luminosities and measured total infrared luminosities L_IR of galaxies detected by Spitzer. Our sample consists of 372 star-forming galaxies with individual detections and flux measurements at 8, 24, 70 and 160 um. We have spectroscopic redshifts for 93% of these sources, and accurate photometric redshifts for the remainder. We also used a stacking analysis to measure the IR fluxes of fainter sources at higher redshifts. We show that the monochromatic mid and far-infrared luminosities are strongly correlated with the total infrared luminosity and our stacking analysis confirms that these correlations also hold at higher redshifts. We provide relations between monochromatic luminosities and total infrared luminosities L_IR that should be reliable up to z~2 (z~1.1) for ULIRGs (LIRGs). In particular, we can predict L_IR with accuracies of 37% and 54% from the 8 and 24 um fluxes, while the best tracer is the 70 um flux. Combining bands leads to slightly more accurate estimates. For example, combining the 8 and 24 um luminosities predicts L_IR with an accuracy of 34%. Our results are generally compatible with previous studies, and the small changes are probably due to differences in the sample selection criteria. We can rule out strong evolution in dust properties with redshift up to z~1. Finally, we show that infrared and sub-millimeter observations are complementary means of building complete samples of star-forming galaxies, with the former being more sensitive for z<~2 and the latter at higher z>~2.

Because of their unique quality, Hubble Space Telescope (HST) data have played an important complementary role in studies of infrared (IR) galaxies conducted with major facilities, as VLT or Spitzer, and will be as well very valuable for future telescopes as Herschel and ALMA. I review here some of the most recent works led by European astronomers on IR galaxies, and discuss the role that HST has had in the study of different IR galaxy populations. I particularly focus the analysis on the GOODS fields, where the multiwavelength data and unique HST coverage have enabled to jointly put constraints on the evolution of star formation activity and stellar-mass growth with cosmic time.

We present the rest-frame 8 micron luminosity function (LF) at redshifts z=1 and ~2, computed from Spitzer 24 micron-selected galaxies in the GOODS fields over an area of 291 sq. arcmin. Using classification criteria based on X-ray data and IRAC colours, we identify the AGN in our sample. The rest-frame 8 micron LF for star-forming galaxies at redshifts z=1 and ~2 have the same shape as at z~0, but with a strong positive luminosity evolution. The number density of star-forming galaxies with log_{10}(nu L_nu(8 micron))>11 increases by a factor >250 from redshift z~0 to 1, and is basically the same at z=1 and ~2. The resulting rest-frame 8 micron luminosity densities associated with star formation at z=1 and ~2 are more than four and two times larger than at z~0, respectively. We also compute the total rest-frame 8 micron LF for star-forming galaxies and AGN at z~2 and show that AGN dominate its bright end, which is well-described by a power-law. Using a new calibration based on Spitzer star-forming galaxies at 0<z<0.6 and validated at higher redshifts through stacking analysis, we compute the bolometric infrared (IR) LF for star-forming galaxies at z=1 and ~2. We find that the respective bolometric IR luminosity densities are (1.2+/-0.2) x 10^9 and (6.6^{+1.2}_{-1.0}) x 10^8 L_sun Mpc^{-3}, in agreement with previous studies within the error bars. At z~2, around 90% of the IR luminosity density associated with star formation is produced by luminous and ultraluminous IR galaxies (LIRG and ULIRG), with the two populations contributing in roughly similar amounts. Finally, we discuss the consistency of our findings with other existing observational results on galaxy evolution.

zCOSMOS is a large redshift survey that is being undertaken in the COSMOS field using 600 hours of observation with the VIMOS spectrograph on the 8-m VLT. The survey is designed to characterise the environments of COSMOS galaxies from the 100 kpc scales of galaxy groups up to the 100 Mpc scale of the cosmic web and to produce diagnostic information on galaxies and active galactic nuclei. The zCOSMOS survey consists of two parts: (a) zCOSMOS-bright, a magnitude-limited I-band IAB < 22.5 sample of about 20,000 galaxies with 0.1 < z < 1.2 covering the whole 1.7 deg2 COSMOS ACS field and designed to mimic the parameters of the 2dfGRS; and (b) zCOSMOS-deep, a survey of approximately 10,000 galaxies selected through colour-selection criteria to have 1.4 < z < 3.0, within the central 1 deg2. This paper describes the survey design and the construction of the target catalogues, and briefly outlines the observational program and the data pipeline. In the first observing season, spectra of 1303 zCOSMOS-bright targets and of 977 zCOSMOS-deep targets have been obtained. These are briefly analysed to demonstrate the characteristics that may be expected from zCOSMOS, and particularly zCOSMOS-bright, when it is finally completed between 2008-2009. The power of combining spectroscopic and photometric redshifts is demonstrated, especially in correctly identifying the emission line in single-line spectra and in determining which of the less reliable spectroscopic redshifts are correct and which are incorrect. Our zCOSMOS-deep spectra demonstrate the effectiveness of our selection techniques to isolate high redshift galaxies at 1.4 < z < 3.0 and of VIMOS to measure their redshifts using ultraviolet absorption lines.

We have performed a deep spectroscopic survey of extremely red galaxies on the GOODS-South field, using GMOS on Gemini South. We present here spectra and redshifts for 16 ERGs at 0.87<z<2.02, to a limit K=20.2. The ERHs are a mixture of spheroidals, mergers and spirals, with one AGN. For at least 10 of these galaxies we observe [OII] emission lines. We perform an age-dating analysis by fitting the spectra and 9-band photometry of the ERGs with models of passively evolving stellar populations combined with a younger star-forming component. The best-fitting ages for the old stellar components range from 0.6 to 4.5 Gyr, with a mean 2.1 Gyr. Masses range from 3 to 20 times 10^10 solar masses. The star-forming component typically forms a few per cent of the total mass, with dust reddening averaging E(B-V)=0.35. Its timescale tends to be short for mergers (<50 Myr) and longer (200-800 Myr) for spiral ERGs.

We investigate the role of the luminous infrared galaxy (LIRG) and ultra-luminous infrared galaxy (ULIRG) phases in the evolution of Ks-selected galaxies and, in particular, Extremely Red Galaxies (ERGs). With this aim, we compare the properties of a sample of 2905 Ks<21.5 (Vega mag) galaxies in the GOODS/CDFS with the sub-sample of those 696 sources which are detected at 24 microns. We find that LIRGs constitute 30% of the galaxies with stellar mass M>1x10^{11} Msun assembled at redshift z=0.5. A minimum of 65% of the galaxies with M>2.5x10^{11} Msun at z~2-3 are ULIRGs at those redshifts. 60% of the ULIRGs in our sample have the characteristic colours of ERGs. Conversely, 40% of the ERGs with stellar mass M>1.3x10^{11} Msun at 1.5<z<2.0 and a minimum of 52% of those with the same mass cut at 2.0<z<3.0 are ULIRGs. The average optical/near-IR properties of the massive ERGs at similar redshifts that are identified with ULIRGs and that are not have basically no difference, suggesting that both populations contain the same kind of objects in different phases of their lives. LIRGs and ULIRGs have an important role in galaxy evolution and mass assembly, and, although they are only able to trace a fraction of the massive (M>1x10^{11} Msun) galaxies present in the Universe at a given time, this fraction becomes very significant (>50%) at redshifts z>~2.

We present the results of our most recent works on Spitzer/MIPS 24 micron galaxies. Through a multiwavelength analysis, we study different properties (redshifts, luminosities, stellar masses) characterising the sources which produce the bulk of the mid-IR background. From a comparative study with the total population of Ks-selected galaxies, we determine that 24 micron sources account for an important fraction of the most massive galaxies present at different redshifts. On the other hand, we determine that 24 micron galaxies also produce most of the energy contained in the far-IR cosmic background at 70 and 160 microns. Furthermore, we are able to set tight constraints on the Cosmic Infrared Background (CIB) spectral energy distribution (SED). Our results help to clarify the links between these presumably different IR galaxy populations.

Aims. We quantify the contributions of 24um galaxies to the Far-Infrared (FIR) Background at 70 and 160um. We provide new estimates of the Cosmic Infrared Background (CIB), and compare it with the Cosmic Optical Background (COB). Methods. Using Spitzer data at 24, 70 and 160um in three deep fields, we stacked more than 19000 MIPS 24um sources with S24>60uJy at 70 and 160um, and measured the resulting FIR flux densities. Results. This method allows a gain up to one order of magnitude in depth in the FIR. We find that the Mid-Infrared (MIR) 24um selected sources contribute to more than 70% of the CIB at 70 and 160um. This is the first direct measurement of the contribution of MIR-selected galaxies to the FIR CIB. Galaxies contributing the most to the total CIB are thus z~1 luminous infrared galaxies, which have intermediate stellar masses. We estimate that the CIB will be resolved at 0.9 mJy at 70 and 3 mJy at 160um. By combining the extrapolation of the 24um source counts below 60uJy, with 160/24 and 70/24 colors as measured with the stacking analysis, we obtain lower limits of 7.1+/-1.0 and 13.4+/-1.7 nW/m2/sr for the CIB at 70 and 160um, respectively. Conclusions. The MIPS surveys have resolved more than three quarters of the MIR and FIR CIB. By carefully integrating the Extragalactic Background Light (EBL) SED, we also find that the CIB has the same brightness as the COB, around 24 nW/m2/sr. The EBL is produced on average by 115 infrared photons for one visible photon. Finally, the galaxy formation and evolution processes emitted a brightness equivalent to 5% of the primordial electromagnetic background (CMB).

We have selected and analysed the properties of a sample of 2905 Ks<21.5 galaxies in ~ 131 sq.arcmin of the Great Observatories Origins Deep Survey (GOODS) Chandra Deep Field South (CDFS), to obtain further constraints on the evolution of Ks-selected galaxies with respect to the results already obtained in previous studies. We made use of the public deep multiwavelength imaging from the optical B through the infrared (IR) 4.5 micron bands, in conjunction with available spectroscopic and COMBO17 data in the CDFS, to construct an optimised redshift catalogue for our galaxy sample. We computed the Ks-band LF and determined that its characteristic magnitude has a substantial brightening and a decreasing total density from z=0 to <z>=2.5. We also analysed the colours and number density evolution of galaxies with different stellar masses. Within our sample, and in contrast to what is observed for less massive systems, the vast majority (~ 85-90%) of the most massive (M>2.5x10^11 Msun) local galaxies appear to be in place before redshift z ~1. Around 65-70% of the total assemble between redshifts z=1 and z=3 and most of them display extremely red colours, suggesting that plausible star formation in these very massive systems should mainly proceed in obscured, short-timescale bursts. The remaining fraction (up to ~ 20%) could be in place at even higher redshifts z=3-4, pushing the first epoch of formation of massive galaxies beyond the limits of current near-IR surveys.

We present deep Ks<21.5 (Vega) identifications, redshifts and stellar masses for most of the sources composing the bulk of the 24 micron background in the GOODS/CDFS. Our identified sample consists of 747 Spitzer/MIPS 24 micron objects, and includes ~94% of all the 24 micron sources in the GOODS-South field which have fluxes Snu(24)>83 microJy (the 80% completeness limit of the Spitzer/GTO 24 micron catalog). 36% of our galaxies have spectroscopic redshifts (mostly at z<1.5) and the remaining ones have photometric redshifts of very good quality, with a median of |dz|=|zspec-zphot|/(1+zspec)=0.02. We find that MIPS 24 micron galaxies span the redshift range z~0-4, and that a substantial fraction (28%) lie at high redshifts z>1.5. We determine the existence of a bump in the redshift distribution at z~1.9, indicating the presence of a significant population of galaxies with PAH emission at these redshifts. Massive (M>10^11 Msun) star-forming galaxies at redshifts 2<z<3 are characterized by very high star-formation rates (SFR>500 Msun/yr), and some of them are able to construct a mass of 10^10-10^11 Msun in a single burst lifetime (~0.01-0.1 Gyr). At lower redshifts z<2, massive star-forming galaxies are also present, but appear to be building their stars on long timescales, either quiescently or in multiple modest burst-like episodes. At redshifts z~1-2, the ability of the burst-like mode to produce entire galaxies in a single event is limited to some lower (M<7x10^10 Msun) mass systems, and it is basically negligible at z<1. Our results support a scenario where star-formation activity is differential with assembled stellar mass and redshift, and where the relative importance of the burst-like mode proceeds in a down-sizing way from high to low redshifts. (abridged)

We present estimated redshifts and derived properties for a sample of 1663 galaxies with Ks <= 22 (Vega), selected from 50.4 sq.arcmin of the GOODS/CDFS field with deep ISAAC imaging, and make an extensive comparison of their properties with those of the extremely red galaxies (ERGs) selected in the same field. We study in detail the evolution of Ks-selected galaxies up to redshifts z ~ 4 and clarify the role of ERGs within the total Ks-band galaxy population. We compute the total Ks-band luminosity function (LF) and compare its evolution with the ERG LF. Up to <z_phot>=2.5, the bright end of the Ks-band LF shows no sign of decline, and is progressively well reproduced by the ERGs with increasing redshift. We also explore the evolution of massive systems present in our sample: up to 20%-25% of the population of local galaxies with assembled stellar mass M>1x10^11 Msun has been formed before redshift z ~ 4, and contains ~ 45% to 70% of the stellar mass density of the Universe at that redshift. Within our sample, the comoving number density of these massive systems is then essentially constant down to redshift z ~ 1.5, by which point most of them have apparently evolved into (I-Ks)-selected ERGs. The remaining massive systems observed in the local Universe are assembled later, at redshifts z <= 1.5. Our results therefore suggest a two-fold assembly history for massive galaxies, in which galaxy/star formation proceeds very efficiently in high mass haloes at very high redshift.

We have analysed 5-epoch GOODS HST-ACS B, V, I_775 and z datasets (V1.0 release), in conjunction with existing VLT-ISAAC imaging in the J, H and Ks bands, to derive estimated redshifts for the sample of 198 Extremely Red Galaxies (ERGs) with Ks<22 (Vega) and I_775-Ks>3.92 selected by Roche et al.(2003) from 50.4 sq.arcmin of the GOODS/CDFS field. We find that, at this depth, the ERG population spans the redshift range 0.5<z_phot<4.75 and over two decades in mass (~ 3x10^9 M_sun to ~3x10^11 M_sun). We explore the evolution of the ERG luminosity function (LF) from redshifts <z_phot>=1.0 to <z_phot>=2.5 and compare it with the global Ks-band LF at redshifts 1<z_phot<2. We find that the bright end of the ERG LF does not decrease from redshifts <z_phot>=2.0 to <z_phot>=2.5 and we connect this fact with the presence of progenitors of the local L>L* population at redshifts z_phot>2. We determine lower limits of rho_c=(6.1 +/- 1.9) x 10^{-5} Mpc^{-3} and rho_c=(2.1 +/- 1.1) x 10^{-5}Mpc^{-3} on the comoving densities of progenitors of local massive galaxies already assembled at redshifts <z_phot>=2.5 and <z_phot>=3.5, respectively. We have investigated the existence of high-redshift Lyman break galaxies massive enough to be included in this ERG sample. Out of an initial list of 12 potential very high redshift candidates, we have identified 2 ERGs which have a high probability of lying at z_phot>4. We discuss the advantages of multi-colour to single-colour selection techniques in obtaining reliable lists of very high-redshift candidate sources, and present revised lower redshift estimates for sources previously claimed as potential z> 5 dropouts in recent studies.(abridged)

We present the photometric redshift distribution of a sample of 198 Extremely Red Galaxies (ERGs) with Ks<22 and I-Ks>3.92 (Vega), selected by Roche et al. in 50.4 sq. arcmin of the Chandra Deep Field South (CDFS). The sample has been obtained using ISAAC-VLT and ACS-HST GOODS public data. We also show the results of a morphological study of the 72 brightest ERGs in the z band (z<25, AB).