We explore X-ray emission from a sample of 18 He II 1640 emitting star-forming galaxies at z ~ 2.3-3.6 from the VANDELS survey in the Chandra Deep Field South, to set constraints on the role of X-ray sources in powering the He II emission. We find that 4 He II emitters have tentative detections with S/N ~ 2 and have X-ray luminosities, L_X = 1.5-4.9 x 10^41 erg/s. The stacked luminosity of all 18 He II emitters is 2.6 x 10^41 erg/s, and that of a subset of 13 narrow He II emitters (FHWM(He II) < 1000 km/s) is 3.1 x 10^41 erg/s. We also measure stacked L_X for non-He II emitters through bootstrapping of matched samples, and find L_X = 2.5 x 10^41 erg/s, which is not significantly different from L_X measured for He II emitters. The L_X per star-formation rate for He II emitters (log (L_X/SFR) ~ 40.0) and non-emitters (log (L_X/SFR) ~ 39.9) are also comparable and in line with the redshift evolution and metallicity dependence predicted by models. Due to the non-significant difference between the X-ray emission from galaxies with and without He II, we conclude that X-ray binaries or weak or obscured AGNs are unlikely to be the dominant producers of He II ionising photons in VANDELS star-forming galaxies at z ~ 3. Given the comparable physical properties of both He II emitters and non-emitters reported previously, alternative He II ionising mechanisms such as localised low-metallicity stellar populations, Pop-III stars, etc. may need to be explored.

We study the projected spatial offset between the ultraviolet continuum and Ly$\alpha$ emission for a sample of 65 lensed and unlensed galaxies in the Epoch of Reionization (5<z<7), the first such study at these redshifts, in order to understand the potential for these offsets to confuse estimates of the Ly$\alpha$ properties of galaxies observed in slit spectroscopy. While we find that ~40% of galaxies in our sample show significant projected spatial offsets ($|\Delta_{Ly\alpha-UV}|$), we find a modest average offset of 0.61$\pm$0.08 kpc. A small fraction of our sample, ~10%, exhibits offsets of 2-4 kpc, sizes that are larger than the effective radii of typical galaxies at these redshifts. An internal comparison and a comparison to studies at lower redshift yielded no significant evidence of evolution of $|\Delta_{Ly\alpha-UV}|$ with redshift. In our own sample, UV-bright galaxies showed offsets a factor of three greater than their fainter counterparts, 0.89$\pm$0.18 vs. 0.27$\pm$0.05 kpc, respectively. We argue that offsets are likely not the result of merging processes, but are rather due to anisotropic internal processes internal to the galaxies as a result of stellar feedback, which facilitates Ly$\alpha$ fluorescence and/or backscattering from nearby or outflowing gas. The reduction in the Ly$\alpha$ flux due to offset effects for various observational setups was quantified through mock observations of simple simulations. It was found that the loss of Ly$\alpha$ photons for galaxies with average offsets is not, if corrected for, a limiting factor for all but the narrowest slit widths (<0.4''). However, for the largest offsets, if such offsets are mostly perpendicular to the slit major axis, slit losses were found to be extremely severe in cases where slit widths of $\le$1'' were employed, such as those planned for James Webb Space Telescope/NIRSpec observations. (abridged)

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

ALMA observations have revealed that [CII] 158$\mu$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z>6 galaxies observed in both [CII] and [OIII]. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect [CII] emission in the whole sample, with an extent systematically larger than the [CII] emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40% of the extended [CII] component might be missed at an angular resolution of 0.8$^{\prime\prime}$, implying that $L_{\rm [CII]}$ is underestimated by a factor $\approx2$ in data at low (<7) signal-to-noise ratio . By combining these results, we conclude that $L_{\rm [CII]}$ of z>6 galaxies lies, on average, slightly below the local $L_{\rm [CII]}-SFR$ relation ($\Delta^{z=6-9}=-0.07\pm0.3$), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm [OIII]}/L_{\rm [CII]}<10$, i.e. more in line with current hydrodynamical simulations.

We present a sample of 151 massive ($M_* > 10^{10}\mathrm{M_\odot}$) quiescent galaxies at $2 < z < 5$, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at $z>3$, of which 2 are at $z>4$. We report formation redshifts, demonstrating that the oldest objects formed at $z > 6$, however individual ages from our photometric data have significant uncertainties, typically $\sim0.5$ Gyr. We demonstrate that the UVJ colours of the quiescent population evolve with redshift at $z>3$, becoming bluer and more similar to post-starburst galaxies at lower redshift. Based upon this we construct a model for the time-evolution of quiescent galaxy UVJ colours, concluding that the oldest objects are consistent with forming the bulk of their stellar mass at $z\sim6-7$ and quenching at $z\sim5$. We report spectroscopic redshifts for two of our objects at $z=3.440$ and $3.396$, which exhibit extremely weak Ly$\alpha$ emission in ultra-deep VANDELS spectra. We calculate star-formation rates based on these line fluxes, finding that these galaxies are consistent with our quiescent selection criteria, provided their Ly$\alpha$ escape fractions are $>3$ and $>10$ per cent respectively. We finally report that our highest-redshift robust object exhibits a continuum break at $\lambda\sim7000$A in a spectrum from VUDS, consistent with our photometric redshift of $z_\mathrm{phot}=4.72^{+0.06}_{-0.04}$. If confirmed as quiescent this object would be the highest-redshift known quiescent galaxy. To obtain stronger constraints on the times of the earliest quenching events, high-SNR spectroscopy must be extended to $z\gtrsim3$ quiescent objects.

The reionization of hydrogen is closely linked to the first structures in the universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly$\alpha$) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly$\alpha$ observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity-halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly$\alpha$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$\alpha$ visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{HI}}$, at $z \sim 7$ using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{HI}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{HI}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies.

Next generation telescopes such as the James Webb Space Telescope (JWST) and the Wide Field Infrared Survey Telescope (WFIRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work we use the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, that couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the Stellar Mass Function (SMF) in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above 100% for M$_{UV} \sim$ -17.5 (-20) at z = 12 (6) for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin$^2$); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at z > 9, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields) of Lyman Break Galaxies during the EoR. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas.

We present the results of a new study investigating the relationship between observed Ly$\alpha$ equivalent width ($W_{\lambda}$(Ly$\alpha$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_{\lambda}$(Ly$\alpha$) across the range $-58 \unicode{xC5} \lesssim$ $W_{\lambda}$(Ly$\alpha$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_{\lambda}$(Ly$\alpha$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=-18\unicode{xC5}$) and the highest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_{\lambda}$(Ly$\alpha$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_{\lambda}$(Ly$\alpha$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_{\lambda}$(Ly$\alpha$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.

Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $\lambda 1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to study their physical properties. After identifying seven potential AGNs in our sample and discarding them from further analysis, we divided the sample of \emph{Bright} emitters into 20 \emph{Narrow} (FWHM < 1000 km s$^{-1}$) and 6 \emph{Broad} (FWHM > 1000 km s$^{-1}$) He II emitters. We created stacks of Faint, Narrow, and Broad emitters and measured other rest-frame UV lines such as O III] and C III] in both individual galaxies and stacks. We then compared the UV line ratios with the output of stellar population-synthesis models to study the ionising properties of He II emitters. We do not see a significant difference between the stellar masses, star-formation rates, and rest-frame UV magnitudes of galaxies with He II and no He II emission. The stellar population models reproduce the observed UV line ratios from metals in a consistent manner, however they under-predict the total number of \heii ionising photons, confirming earlier studies and suggesting that additional mechanisms capable of producing He II are needed, such as X-ray binaries or stripped stars. The models favour subsolar metallicities ($\sim0.1Z_\odot$) and young stellar ages ($10^6 - 10^7$ years) for the He II emitters. However, the metallicity measured for He II emitters is comparable to that of non-He II emitters at similar redshifts. We argue that galaxies with He II emission may have undergone a recent star-formation event, or may be powered by additional sources of He II ionisation.

We present the star formation histories of 39 galaxies with high quality rest-frame optical spectra at 0.5<z<1.3 selected to have strong Balmer absorption lines and/or Balmer break, and compare to a sample of spectroscopically selected quiescent galaxies at the same redshift. Photometric selection identifies a majority of objects that have clear evidence for a recent short-lived burst of star formation within the last 1.5 Gyr, i.e. "post-starburst" galaxies, however we show that good quality continuum spectra are required to obtain physical parameters such as burst mass fraction and burst age. Dust attenuation appears to be the primary cause for misidentification of post-starburst galaxies, leading to contamination in spectroscopic samples where only the [OII] emission line is available, as well as a small fraction of objects lost from photometric samples. The 31 confirmed post-starburst galaxies have formed 40-90% of their stellar mass in the last 1-1.5 Gyr. We use the derived star formation histories to find that the post-starburst galaxies are visible photometrically for 0.5-1 Gyr. This allows us to update a previous analysis to suggest that 25-50% of the growth of the red sequence at z~1 could be caused by a starburst followed by rapid quenching. We use the inferred maximum historical star formation rates of several 100-1000 Msun/yr and updated visibility times to confirm that sub-mm galaxies are likely progenitors of post-starburst galaxies. The short quenching timescales of 100-200 Myr are consistent with cosmological hydrodynamic models in which rapid quenching is caused by the mechanical expulsion of gas due to an AGN.

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 36 spectroscopically confirmed intrinsically UV-faint Ly$\alpha$ emitting galaxies from follow-up observations with Keck/DEIMOS of gravitationally lensed high-redshift candidates. Candidates were selected to be between $5\lesssim z \lesssim 7$ from photometric data using \textit{HST} and \textit{Spitzer} imaging surveys. We used photometric redshift information to perform an integrated photometric redshift probability cut $>1\%$ between $5<z<7$ to construct a sample of 198 high-redshift objects. Our high-redshift sample spans intrinsic UV luminosities from a few $L^*$ down to $0.001L^*$. We identified 19 high-confidence detections of Ly$\alpha$ and an additional 17 likely detections. We divided our sample into lower-redshift ($z\sim5.5$) and higher-redshift ($z\sim6.5$) bins and ran Monte Carlo trials, incorporating the strength of the Ly$\alpha$ emission and the photometric redshift of the non-detections. Considering only objects where Ly$\alpha$ could be detected at EW(Ly$\alpha$)$>$25{\AA} at $3\sigma$ at the fiducial depth of our survey, and only those galaxies with EW(Ly$\alpha$)$>$25{\AA} as true LAEs, and finally, only objects with $m_{AB}<26.8$, we found the LAE fraction to be flat, or modestly increase from 0.26$\pm0.04$ to 0.30$\pm0.04$. These values relative to those for lower-redshift samples are consistent with a rising LAE fraction with redshift out to $z\sim6$, but at $z\sim6.5$ there is some tension between our results and results from surveys at intrinsically brighter luminosities. We conclude intrinsically fainter galaxies have Ly$\alpha$ emission, and there is a steep drop in the LAE fraction from our high-redshift sample at $z\sim6.5$ and from similar galaxies at $z\sim7.5$. This likely indicates we are witnessing the tail end of the epoch of reionization, as such a drop is not expected due to changes of intrinsic galaxy properties between these redshifts.

We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies.

We report on the HST detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z=3.794, dubbed Ion1 (Vanzella et al. 2012). The LyC from Ion1 is detected at rest-frame wavelength 820$\sim$890 \AA with HST WFC3/UVIS in the F410M band ($m_{410}=27.60\pm0.36$ magnitude (AB), peak SNR = 4.17 in a circular aperture with radius r = 0.12'') and at 700$\sim$830 \AA with the VLT/VIMOS in the U-band ($m_U = 27.84\pm0.19$ magnitude (AB), peak SNR = 6.7 with a r = 0.6'' aperture). A 20-hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines, the P-Cygni profile of CIV and Ly$\alpha$ in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly$\alpha$ emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest, presumably because of lower dust obscuration. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the non-ionizing UV emission by 0.12''$\pm$0.03'', corresponding to 0.85$\pm$0.21 kpc (physical), and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low HI column-density "cavities" in the ISM, possibly carved by stellar winds and/or supernova. We also collect the VIMOS U-band images of a sample of 107 Lyman-break galaxies with spectroscopic redshifts at $3.40<z<3.95$, i.e. sampling the LyC, and stack them with inverse-variance weights. No LyC emission is detected in the stacked image, resulting in a 32.5 magnitude (AB) flux limit (1$\sigma$) and an upper limit of absolute LyC escape fraction $f_{esc}^{abs} < 0.63\%$. LyC emitters like Ion1 are very likely at the bright-end of the LyC luminosity function.

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We firstly constrain the stellar mass vs stellar age relationship, finding a steep trend towards earlier average formation with increasing stellar mass of $1.48^{+0.34}_{-0.39}$ Gyr per decade in mass, although this shows signs of flattening at $M_* > 10^{11} M_\odot$. We show that this is consistent with other spectroscopic studies from $0 < z < 2$. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by Simba and IllustrisTNG agree well with observations at $z=0.1$, they are too shallow at $z=1$, predicting an evolution of $<0.5$ Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies ($M_* \sim 10^{10.5} M_\odot$) are consistent with formation in recent ($z<2$), intense starburst events, with timescales of $\lesssim500$ Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by $z=1$, they formed at $z\sim5$ and quenched by $z=3$. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity.

We present spectra and near-infrared images of a sample of faint radio sources initially selected as promising high-redshift radio galaxy (HzRG) candidates. We have determined redshifts for a total of 13 radio galaxies with redshifts ranging from $0.52\le z \le 5.72$. Our sample probes radio luminosities that are almost an order of magnitude fainter than previous large samples at the highest redshifts. We use near-infrared photometry for a subsample of these galaxies to calculate stellar masses using simple stellar population models, and find stellar masses to be in the range $10^{10.8} - 10^{11.7} M_\odot$. We then compare our faint radio galaxies with brighter radio galaxies at $z\ge2$ from the literature. We find that fainter radio galaxies have lower Ly$\alpha$ luminosities and narrower full width at half maxima compared to the bright ones, implying photoionisation by weaker AGN. We also rule out the presence of strong shocks in faint HzRGs. The stellar masses determined for faint HzRGs are lower than what has been observed for brighter ones. We find that faint HzRG population in the redshift range $2-4$ forms a bridge between star-forming and narrow-line AGN, whereas the ones at $z>4$ are dominated by star-formation, and may be building up their stellar mass through cold accretion of gas. Finally, we show that the overall redshift evolution of radio sizes at $z>2$ is fully compatible with increased inverse Compton scattering losses at high redshifts.

We search the five CANDELS fields (COSMOS, EGS, GOODS-North, GOODS-South and UDS) for passively evolving a.k.a. "red and dead" massive galaxies in the first 2 Gyr after the Big Bang, integrating and updating the work on GOODS-South presented in our previous paper. We perform SED-fitting on photometric data, with top-hat star-formation histories to model an early and abrupt quenching, and using a probabilistic approach to select only robust candidates. Using libraries without (with) spectral lines emission, starting from a total of more than 20,000 $z>3$ sources we end up with 102 (40) candidates, including one at $z=6.7$. This implies a minimal number density of $1.73 \pm 0.17 \times 10^{-5}$ ($6.69 \pm 1.08 \times 10^{-6}$) Mpc$^{-3}$ for $3<z<5$; applying a correction factor to account for incompleteness yields $2.30 \pm 0.20 \times 10^{-5}$. We compare these values with those from five recent hydrodynamical cosmological simulations, finding a reasonable agreement at $z<4$; tensions arise at earlier epochs. Finally, we use the star-formation histories from the best-fit models to estimate the contribution of the high-redshift passive galaxies to the global Star Formation Rate Density during their phase of activity, finding that they account for $\sim5-10\%$ of the total star formation at $3<z<8$, despite being only $\sim0.5\%$ of the total in number. The resulting picture is that early and strong star formation activity, building massive galaxies on short timescales and followed by a quick and abrupt quenching, is a rare but crucial phenomenon in the early Universe: the evolution of the cosmos must be heavily influenced by the short but powerful activity of these pristine monsters.

The study of the space density of bright AGNs at $z>4$ has been subject to extensive effort given its importance for the estimate of the cosmological ionizing emissivity and growth of supermassive black holes. In this context we have recently derived high space densities of AGNs at $z\sim 4$ and $-25<M_{1450}<-23$ in the COSMOS field from a spectroscopically complete sample. In the present paper we attempt to extend the knowledge of the AGN space density at fainter magnitudes ($-22.5<M_{1450}<-18.5$) in the $4<z<6.1$ redshift interval by means of a multiwavelength sample of galaxies in the CANDELS GOODS-South, GOODS-North and EGS fields. We use an updated criterion to extract faint AGNs from a population of NIR (rest-frame UV) selected galaxies at photometric $z>4$ showing X-ray detection in deep Chandra images available for the three CANDELS fields. We have collected a photometric sample of 32 AGN candidates in the selected redshift interval, six of which having spectroscopic redshifts. Including our COSMOS sample as well as other bright QSO samples allows a first guess on the shape of the UV luminosity function at $z\sim 4.5$. The resulting emissivity and photoionization rate appear consistent with that derived from the photoionization level of the intergalactic medium at $z\sim 4.5$. An extrapolation to $z\sim 5.6$ suggests an important AGN contribution to the IGM ionization if there are no significant changes in the shape of the UV luminosity function.

We investigate the prevalence of galactic-scale outflows in post-starburst (PSB) galaxies at high redshift ($1 < z < 1.4$), using the deep optical spectra available in the UKIDSS Ultra Deep Survey (UDS). We use a sample of $\sim40$ spectroscopically confirmed PSBs, recently identified in the UDS field, and perform a stacking analysis in order to analyse the structure of strong interstellar absorption features such as Mg ii ($\lambda2800$ Ang.). We find that for massive ($M_* > 10^{10}\rm\,M_{\odot}$) PSBs at $z > 1$, there is clear evidence for a strong blue-shifted component to the Mg ii absorption feature, indicative of high-velocity outflows ($v_{\rm out}\sim1150\pm160\rm\,km\,s^{-1}$) in the interstellar medium. We conclude that such outflows are typical in massive PSBs at this epoch, and potentially represent the residual signature of a feedback process that quenched these galaxies. Using full spectral fitting, we also obtain a typical stellar velocity dispersion $\sigma_*$ for these PSBs of $\sim200\rm\,km\,s^{-1}$, which confirms they are intrinsically massive in nature (dynamical mass $M_{\rm d}\sim10^{11}\rm\,M_{\odot}$). Given that these high-$z$ PSBs are also exceptionally compact ($r_{\rm e}\sim1$--$2\rm\,kpc$) and spheroidal (Sersic index $n\sim3$), we propose that the outflowing winds may have been launched during a recent compaction event (e.g. major merger or disc collapse) that triggered either a centralised starburst or active galactic nuclei (AGN) activity. Finally, we find no evidence for AGN signatures in the optical spectra of these PSBs, suggesting they were either quenched by stellar feedback from the starburst itself, or that if AGN feedback is responsible, the AGN episode that triggered quenching does not linger into the post-starburst phase.

We recommend a conceptual design study for a spectroscopic facility in the southern hemisphere comprising a large diameter telescope, fiber system, and spectrographs collectively optimized for massively-multiplexed spectroscopy. As a baseline, we propose an 11.4-meter aperture, optical spectroscopic survey telescope with a five square degree field of view. Using current technologies, the facility could be equipped with 15,000 robotically-controlled fibers feeding spectrographs over 360<lambda<1330 nm with options for fiber-fed spectrographs at high resolution and a panoramic IFU at a separate focus. This would enable transformational progress via its ability to access a larger fraction of objects from Gaia, LSST, Euclid, and WFIRST than any currently funded or planned spectroscopic facility. An ESO-sponsored study (arXiv:1701.01976) discussed the scientific potential in ambitious new spectroscopic surveys in Galactic astronomy, extragalactic astronomy, and cosmology. The US community should establish links with European and other international communities to plan for such a powerful facility and maximize the potential of large aperture multi-object spectroscopy given the considerable investment in deep imaging surveys.

We constrain the distribution of spatially offset Lyman-alpha emission (Ly$\alpha$) relative to rest-frame ultraviolet emission in $\sim300$ high redshift ($3<z<5.5$) Lyman-break galaxies (LBGs) exhibiting Ly$\alpha$ emission from VANDELS, a VLT/VIMOS slit-spectroscopic survey of the CANDELS Ultra Deep Survey and Chandra Deep Field South fields (${\simeq0.2}~\mathrm{deg}^2$ total). Because slit spectroscopy compresses two-dimensional spatial information into one spatial dimension, we use Bayesian inference to recover the underlying Ly$\alpha$ spatial offset distribution. We model the distribution using a 2D circular Gaussian, defined by a single parameter $\sigma_{r,\mathrm{Ly}\alpha}$, the standard deviation expressed in polar coordinates. Over the entire redshift range of our sample ($3<z<5.5$), we find $\sigma_{r,\mathrm{Ly}\alpha}=1.70^{+0.09}_{-0.08}$ kpc ($68\%$ conf.), corresponding to $\sim0.25$ arcsec at $\langle z\rangle=4.5$. We also find that $\sigma_{r,\mathrm{Ly}\alpha}$ decreases significantly with redshift. Because Ly$\alpha$ spatial offsets can cause slit-losses, the decrease in $\sigma_{r,\mathrm{Ly}\alpha}$ with redshift can partially explain the increase in the fraction of Ly$\alpha$ emitters observed in the literature over this same interval, although uncertainties are still too large to reach a strong conclusion. If $\sigma_{r,\mathrm{Ly}\alpha}$ continues to decrease into the reionization epoch, then the decrease in Ly$\alpha$ transmission from galaxies observed during this epoch might require an even higher neutral hydrogen fraction than what is currently inferred. Conversely, if spatial offsets increase with the increasing opacity of the IGM, slit losses may explain some of the drop in Ly$\alpha$ transmission observed at $z>6$. Spatially resolved observations of Ly$\alpha$ and UV continuum at $6<z<8$ are needed to settle the issue.

We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($\eta=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $\eta \propto M_{\ast}^{\beta}$ with $\beta \simeq -0.4$.

We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.

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

Detections and non-detections of Lyman alpha (Ly$\alpha$) emission from $z>6$ galaxies ($<1$ Gyr after the Big Bang) can be used to measure the timeline of cosmic reionization. Of key interest to measuring reionization's mid-stages, but also increasing observational challenge, are observations at z > 7, where Ly$\alpha$ redshifts to near infra-red wavelengths. Here we present a search for z > 7.2 Ly$\alpha$ emission in 53 intrinsically faint Lyman Break Galaxy candidates, gravitationally lensed by massive galaxy clusters, in the KMOS Lens-Amplified Spectroscopic Survey (KLASS). With integration times of ~7-10 hours, we detect no Ly$\alpha$ emission with S/N>5 in our sample. We determine our observations to be 80% complete for 5$\sigma$ spatially and spectrally unresolved emission lines with integrated line flux $>5.7\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We define a photometrically selected sub-sample of 29 targets at $z=7.9\pm0.6$, with a median 5$\sigma$ Ly$\alpha$ EW limit of 58A. We perform a Bayesian inference of the average intergalactic medium (IGM) neutral hydrogen fraction using their spectra. Our inference accounts for the wavelength sensitivity and incomplete redshift coverage of our observations, and the photometric redshift probability distribution of each target. These observations, combined with samples from the literature, enable us to place a lower limit on the average IGM neutral hydrogen fraction of $> 0.76 \; (68\%), \; > 0.46 \; (95\%)$ at z ~ 8, providing further evidence of rapid reionization at z~6-8. We show that this is consistent with reionization history models extending the galaxy luminosity function to $M_\textrm{UV} \lesssim -12$, with low ionizing photon escape fractions, $f_\textrm{esc} \lesssim 15\%$.

We analyse a unique sample of 517 bright ($L>L^{*}$) LBGs at redshift z$\sim$3 in order to characterise the distribution of their UV slopes $\beta$ and infer their dust extinction under standard assumptions. We exploited multi-band observations over 750 arcmin$^2$ of the COSMOS field that were acquired with three different ground-based facilities: the Large Binocular Camera (LBC) on the Large Binocular Telescope (LBT), the Suprime-Cam on the SUBARU telescope, and the VIRCAM on the VISTA telescope (ULTRAVISTA DR2). Our multi-band photometric catalogue is based on a new method that is designed to maximise the signal-to-noise ratio in the estimate of accurate galaxy colours from images with different point spread functions (PSF). We adopted an improved selection criterion based on deep Y-band data to isolate a sample of galaxies at $z\sim 3 $ to minimise selection biases. We measured the UV slopes ($\beta$) of the objects in our sample and then recovered the intrinsic probability density function of $\beta$ values (PDF($\beta$)), taking into account the effect of observational uncertainties through detailed simulations. The galaxies in our sample are characterised by mildly red UV slopes with $<\beta>\simeq -1.70$ throughout the enitre luminosity range that is probed by our data ($-24\lesssim M_{1600}\lesssim -21$). The resulting dust-corrected star formation rate density (SFRD) is $log(SFRD)\simeq-1.6 M_{\odot}/yr/Mpc^{3}$, corresponding to a contribution of about 25% to the total SFRD at z$\sim$3 under standard assumptions. Ultra-bright LBGs at $z \sim 3$ match the known trends, with UV slopes being redder at decreasing redshifts, and brighter galaxies being more highly dust extinct and more frequently star-forming than fainter galaxies. [abridged]

The selection of red, passive galaxies in the early Universe is very challenging, especially beyond z~3, and it is crucial to constrain theoretical modelling of the processes responsible for their rapid assembly and abrupt shut-down of the star formation. We present here the analysis of ALMA archival observations of 26 out of the 30 galaxies in the deep CANDELS GOODS-South field that we identified as passive at z~3-5 by means of a careful and conservative SED fitting analysis. ALMA data are used to verify the potential contamination from red, dusty but star--forming sources that could enter the sample due to similar optical--nearIR colours. With the exception of a few marginal detections at <3sigma, we could only infer upper limits, both on individual sources and on the stacks. We translated the ALMA continuum measurements into corresponding SFRs, using a variety of far-IR models. These SFRs are compared with those predicted by secondary star-forming solutions of the optical fits and with the expected position of the star formation Main Sequence. This analysis confirms the passive nature of 9 candidates with high confidence and suggests that the classification is correct for at least half of the sample in a statistical sense. For the remaining sources the analysis remain inconclusive because available ALMA data is not deep enough, although the stacking results corroborate their passive nature. Despite the uncertainties, this work provides decisive support to the existence of passive galaxies beyond z~3.

Realizing the utility of Ly$\alpha$ emission to trace the evolution of the intergalactic medium (IGM) during the epoch of reionization requires deep spectroscopy across the boundary of optical and near-infrared (NIR) spectrographs at $z\sim7.2$ when Ly$\alpha$ emission is at $\sim$1$\mu$m. Our Texas Spectroscopic Search for Ly$\alpha$ Emission at the End of Reionization includes 18 nights of deep spectroscopic observations using the Keck DEIMOS (optical) and MOSFIRE (NIR) spectrographs. Within this dataset we observe Ly$\alpha$ emission from 183 photometric-redshift selected galaxies at $z =$ 5.5 - 8.3 from the Cosmic Assembly Near infrared Deep Extragalactic Legacy Survey (CANDELS). Our overlapping MOSFIRE observations, over 84 galaxies total, provide the deepest NIR spectroscopic data yet obtained for Ly$\alpha$ from galaxies $z > 7$, with $>16$ hr integration time for four observed galaxies. Here we analyze these four targets, and we report the discovery of a new $z = 7.60$ Ly$\alpha$ detection as well as provide an updated observation of the previously confirmed $z=7.51$ Ly$\alpha$ emission from Finkelstein et al. (2013) with a $\sim$3$\times$ longer exposure time. Our analysis of these Ly$\alpha$ emission line profiles reveal a significant asymmetric shape. The two detected Ly$\alpha$ emission lines from bright sources ($M_{\text{UV}}<-20.25$) could imply that these bright galaxies inhabit ionized bubbles in a partially neutral IGM, although deeper exposures may yet reveal Ly$\alpha$ emission in the fainter sources.

We provide access to the current VANDELS data release (DR2) through a dedicated database, accessible at the link http://vandels.inaf.it/db. The data release includes the spectra for all galaxies for which the scheduled integration time was completed during the first two observing seasons and the spectra for those galaxies for which the scheduled integration time was 50\% complete at the end of season two (i.e. 20/40 hours and 40/80 hours). The total number of spectra released is 1362 (586 in CDFS and 776 in UDS). Spectra are stored as multi-extension FITS files containing the 1D extracted spectrum, the 2D linearly re-sampled spectrum, the 1D sky spectrum, the 1D noise estimate and the image thumbnail of the object.

We report the discovery of the most distant radio galaxy to date, TGSS1530 at a redshift of $z=5.72$ close to the presumed end of the Epoch of Reionisation. The radio galaxy was selected from the TGSS ADR1 survey at 150 MHz for having to an ultra-steep spectral index, $\alpha^{\textrm{150 MHz}}_{\textrm{1.4 GHz}} = -1.4$ and a compact morphology obtained using VLA imaging at 1.4 GHz. No optical or infrared counterparts for it were found in publicly available sky surveys. Follow-up optical spectroscopy at the radio position using GMOS on Gemini North revealed the presence of a single emission line. We identify this line as Lyman alpha at $z=5.72$, because of its asymmetric line profile, the absence of other optical/UV lines in the spectrum and a high equivalent width. With a Ly$\alpha$ luminosity of $5.7 \times 10^{42}$ erg s$^{-1}$ and a FWHM of $370$ km s$^{-1}$, TGSS1530 is comparable to `non-radio' Lyman alpha emitters (LAEs) at a similar redshift. However, with a radio luminosity of $\log L_{\textrm{150 MHz}} = 29.1$ W Hz$^{-1}$ and a deconvolved physical size $3.5$ kpc, its radio properties are similar to other known radio galaxies at $z>4$. Subsequent $J$ and $K$ band imaging using LUCI on the Large Binocular Telescope resulted in non-detection of the host galaxy down to $3\sigma$ limits of $J>24.4$ and $K>22.4$ (Vega). The $K$ band limit is consistent with $z>5$ from the $K-z$ relation for radio galaxies, suggesting stellar mass limits using simple stellar population models of $M_{\textrm{stars}}< 10^{10.5}$ $M_\odot$. Its high redshift coupled with relatively small radio and Ly$\alpha$ sizes suggest that TGSS1530 may be a radio galaxy in an early phase of its evolution.

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

We present the results of CANDELSz7, an ESO large program aimed at confirming spectroscopically a homogeneous sample of z~6 and z~7 star forming galaxies. The candidates were selected in the GOODS-South, UDS and COSMOS fields using the official CANDELS catalogs based on H160-band detections. Standard color criteria, which were tailored depending on the ancillary multi-wavelength data available for each field, were applied to select more than 160 candidate galaxies at z~6 and z~7. Deep medium resolution FORS2 spectroscopic observations were then conducted with integration times ranging from 12 to 20 hours, to reach a Lyalpha flux limit of approximately 1-3x 10-18 erg/s/cm^2 at 3sigma. For about 40% of the galaxies we could determine a spectroscopic redshift, mainly through the detection of a single emission line that we interpret as Lyalpha emission, or for some of the brightest objects (H160< 25.5) from the presence of faint continuum and sharp drop that we interpret as a Lyman break. In this paper we present the redshifts and main properties of 65 newly confirmed high redshift galaxies. Adding previous proprietary and archival data we assemble a sample of ~260 galaxies that we use to explore the evolution of the Lyalpha fraction in Lyman break galaxies and the change in the shape of the emission line between z~6 and z~7. We also discuss the accuracy of the CANDELS photometric redshifts in this redshift range.

We present deep spectroscopic follow-up observations of the Bremer Deep Field (BDF) where the two $z\sim$7 bright Ly$\alpha$ emitters (LAE) BDF521 and BDF3299 were previously discovered by Vanzella et al. (2011) and where a factor of $\sim$3-4 overdensity of faint LBGs has been found by Castellano et al. (2016). We confirm a new bright Ly$\alpha$ emitter, BDF2195, at the same redshift of BDF521, $z=7.008$, and at only $\sim$90 kpc physical distance from it, confirming that the BDF area is likely an overdense, reionized region. A quantitative assessment of the Ly$\alpha$ fraction shows that the number of detected bright emitters is much higher than the average found at z$\sim$7, suggesting a high Ly$\alpha$ transmission through the inter-galactic medium (IGM). However, the line visibility from fainter galaxies is at odds with this finding, since no Ly$\alpha$ emission is found in any of the observed candidates with $M_{UV}>$-20.25. This discrepancy can be understood either if some mechanism prevents Ly$\alpha$ emission from fainter galaxies within the ionized bubbles from reaching the observer, or if faint galaxies are located outside the reionized area and bright LAEs are solely responsible for the creation of their own HII regions. A thorough assessment of the nature of the BDF region and of its sources of re-ionizing radiation will be made possible by JWST spectroscopic capabilities.

When combined with the huge collecting area of the ELT, MOSAIC will be the most effective and flexible Multi-Object Spectrograph (MOS) facility in the world, having both a high multiplex and a multi-Integral Field Unit (Multi-IFU) capability. It will be the fastest way to spectroscopically follow-up the faintest sources, probing the reionisation epoch, as well as evaluating the evolution of the dwarf mass function over most of the age of the Universe. MOSAIC will be world-leading in generating an inventory of both the dark matter (from realistic rotation curves with MOAO fed NIR IFUs) and the cool to warm-hot gas phases in z=3.5 galactic haloes (with visible wavelenth IFUs). Galactic archaeology and the first massive black holes are additional targets for which MOSAIC will also be revolutionary. MOAO and accurate sky subtraction with fibres have now been demonstrated on sky, removing all low Technical Readiness Level (TRL) items from the instrument. A prompt implementation of MOSAIC is feasible, and indeed could increase the robustness and reduce risk on the ELT, since it does not require diffraction limited adaptive optics performance. Science programmes and survey strategies are currently being investigated by the Consortium, which is also hoping to welcome a few new partners in the next two years.

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.

Recent detections of Lyman alpha (Ly$\alpha$) emission from $z>7.5$ galaxies were somewhat unexpected given a dearth of previous non-detections in this era when the intergalactic medium (IGM) is still highly neutral. But these detections were from UV bright galaxies, which preferentially live in overdensities which reionize early, and have significantly Doppler-shifted Ly$\alpha$ line profiles emerging from their interstellar media (ISM), making them less affected by the global IGM state. Using a combination of reionization simulations and empirical ISM models we show, as a result of these two effects, UV bright galaxies in overdensities have $>2\times$ higher transmission through the $z\sim7$ IGM than typical field galaxies, and this boosted transmission is enhanced as the neutral fraction increases. The boosted transmission is not sufficient to explain the observed high Ly$\alpha$ fraction of $M_\mathrm{UV} \lesssim -22$ galaxies (Stark et al. 2017), suggesting Ly$\alpha$ emitted by these galaxies must be stronger than expected due to enhanced production and/or selection effects. Despite the bias of UV bright galaxies to reside in overdensities we show Ly$\alpha$ observations of such galaxies can accurately measure the global neutral hydrogen fraction, particularly when Ly$\alpha$ from UV faint galaxies is extinguished, making them ideal candidates for spectroscopic follow-up into the cosmic Dark Ages.

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

We present a new flexible Bayesian framework for directly inferring the fraction of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR, z~6-10) from detections and non-detections of Lyman Alpha (Ly$\alpha$) emission from Lyman break galaxies (LBGs). Our framework combines sophisticated reionization simulations with empirical models of the interstellar medium (ISM) radiative transfer effects on Ly$\alpha$. We assert that the Ly$\alpha$ line profile emerging from the ISM has an important impact on the resulting transmission of photons through the IGM, and that these line profiles depend on galaxy properties. We model this effect by considering the peak velocity offset of Ly$\alpha$ lines from host galaxies' systemic redshifts, which are empirically correlated with UV luminosity and redshift (or halo mass at fixed redshift). We use our framework on the sample of LBGs presented in Pentericci et al. (2014) and infer a global neutral fraction at z~7 of $\overline{x}_\mathrm{HI} = 0.59_{-0.15}^{+0.11}$, consistent with other robust probes of the EoR and confirming reionization is on-going ~700 Myr after the Big Bang. We show that using the full distribution of Ly$\alpha$ equivalent width detections and upper limits from LBGs places tighter constraints on the evolving IGM than the standard Ly$\alpha$ emitter fraction, and that larger samples are within reach of deep spectroscopic surveys of gravitationally lensed fields and JWST NIRSpec.

We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2<z<3$ in the \emph{Hubble} \emph{Ultra} \emph{Deep} \emph{Field} (HUDF), as well as empirical $A_V - M_{\star}$ relations predicted by a Calzetti-like law. Our results, combined with other literature data, suggest that the $A_V - M_{\star}$ relation does not evolve over the redshift range $0<z<5$, at least for galaxies with log$(M_{\star}/M_{\odot}) \gtrsim 9.5$. Finally, we present tentative evidence which suggests that the attenuation curve may become steeper at log$(M_{\star}/M_{\odot}) \lesssim 9.0$.

Finding the sources responsible for the hydrogen reionization is one of the most pressing issues in cosmology. Bright QSOs are known to ionize their surrounding neighborhood, but they are too few to ensure the required HI ionizing background. A significant contribution by faint AGNs, however, could solve the problem, as recently advocated on the basis of a relatively large space density of faint active nuclei at z>4. We have carried out an exploratory spectroscopic program to measure the HI ionizing emission of 16 faint AGNs spanning a broad U-I color interval, with I~21-23 and 3.6<z<4.2. These AGNs are three magnitudes fainter than the typical SDSS QSOs (M1450<~-26) which are known to ionize their surrounding IGM at z>~4. The LyC escape fraction has been detected with S/N ratio of ~10-120 and is between 44 and 100% for all the observed faint AGNs, with a mean value of 74% at 3.6<z<4.2 and -25.1<M1450<-23.3, in agreement with the value found in the literature for much brighter QSOs (M1450<~-26) at the same redshifts. The LyC escape fraction of our faint AGNs does not show any dependence on the absolute luminosities or on the observed U-I colors. Assuming that the LyC escape fraction remains close to ~75% down to M1450~-18, we find that the AGN population can provide between 16 and 73% (depending on the adopted luminosity function) of the whole ionizing UV background at z~4, measured through the Lyman forest. This contribution increases to 25-100% if other determinations of the ionizing UV background are adopted. Extrapolating these results to z~5-7, there are possible indications that bright QSOs and faint AGNs can provide a significant contribution to the reionization of the Universe, if their space density is high at M1450~-23.

We present deep spectroscopic observations of a Lyman-break galaxy candidate (hereafter MACS1149-JD) at $z\sim9.5$ with the $\textit{Hubble}$ Space Telescope ($\textit{HST}$) WFC3/IR grisms. The grism observations were taken at 4 distinct position angles, totaling 34 orbits with the G141 grism, although only 19 of the orbits are relatively uncontaminated along the trace of MACS1149-JD. We fit a 3-parameter ($z$, F160W mag, and Ly$\alpha$ equivalent width) Lyman-break galaxy template to the three least contaminated grism position angles using an MCMC approach. The grism data alone are best fit with a redshift of $z_{\mathrm{grism}}=9.53^{+0.39}_{-0.60}$ ($68\%$ confidence), in good agreement with our photometric estimate of $z_{\mathrm{phot}}=9.51^{+0.06}_{-0.12}$ ($68\%$ confidence). Our analysis rules out Lyman-alpha emission from MACS1149-JD above a $3\sigma$ equivalent width of 21 \AA{}, consistent with a highly neutral IGM. We explore a scenario where the red $\textit{Spitzer}$/IRAC $[3.6] - [4.5]$ color of the galaxy previously pointed out in the literature is due to strong rest-frame optical emission lines from a very young stellar population rather than a 4000 \AA{} break. We find that while this can provide an explanation for the observed IRAC color, it requires a lower redshift ($z\lesssim9.1$), which is less preferred by the $\textit{HST}$ imaging data. The grism data are consistent with both scenarios, indicating that the red IRAC color can still be explained by a 4000 \AA{} break, characteristic of a relatively evolved stellar population. In this interpretation, the photometry indicate that a $340^{+29}_{-35}$ Myr stellar population is already present in this galaxy only $\sim500~\mathrm{Myr}$ after the Big Bang.

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

We present a large-scale galaxy structure Cl J021734-0513 at z~0.65 discovered in the UKIDSS UDS field, made of ~20 galaxy groups and clusters, spreading over 10 Mpc. We report on a VLT/VIMOS spectroscopic follow-up program that, combined with past spectroscopy, allowed us to confirm four galaxy clusters (M200 ~ 10^14 Msol) and a dozen associated groups and star-forming galaxy overdensities. Two additional filamentary structures at z ~ 0.62 and z ~ 0.69 and foreground and background clusters at 0.6 < z < 0.7 were also confirmed along the line-of-sight. The structure subcomponents are at different formation stages. The clusters have a core dominated by passive galaxies and an established red sequence. The remaining structures are a mix of star-forming galaxy overdensities and forming groups. The presence of quiescent galaxies in the core of the latter shows that `preprocessing' has already happened before the groups fall into their more massive neighbours. Our spectroscopy allows us to derive spectral index measurements e.g. emission/absorption line equivalent widths, strength of the 4000A break, valuable to investigate the star formation history of structure members. Based on these line measurements, we select a population of `poststarburst' galaxies. These galaxies are preferentially found within the virial radius of clusters, supporting a scenario in which their recent quenching could be prompted by gas stripping by the dense intracluster medium. We derive stellar age estimates using MCMC-based spectral fitting for quiescent galaxies and find a correlation between ages and colours/stellar masses which favours a top-down formation scenario of the red sequence. A catalogue of ~650 redshifts in UDS will be released alongside the paper.

Context: Ultraviolet (UV) emission-line spectra are used to spectroscopically confirm high-z galaxies and increasingly also to determine their physical properties. Aims: To interpret the observed UV spectra of distant galaxies in terms of the dominant radiation field and the physical condition of the interstellar medium. Methods: We construct a large grid of photoionization models and derive new spectral UV line diagnostics using equivalent widths (EWs) of CIII]1909, CIV1549 and the line ratios of CIII], CIV, and HeII1640 recombination lines. We apply these diagnostics to a sample of 450 CIII]-emitting galaxies at z=2-4 previously identified in the VIMOS Ultra Deep Survey. Results: We show that the average star-forming galaxy (EW(CIII])~2A) is well described by stellar photoionization from single and binary stars. The inferred metallicity and ionization parameter is typically Z=0.3-0.5Zsun and logU=-2.7 to -3, in agreement with earlier works at similar redshifts. The models also indicate an average age of 50-200Myr since the beginning of the current star-formation, and an ionizing photon production rate, xi_ion, of log(xi_ion/[Hz/erg])~25.3-25.4. Among the sources with EW(CIII])=10-20A, ~30% are likely dominated by AGNs. Their derived metallicity is low, Z=0.02-0.2Zsun, and the ionization parameter higher (logU=-1.7). To explain the average UV observations of the strongest but rarest CIII] emitters (EW(CIII])>20A), we find that stellar photoionization is clearly insufficient. A radiation field consisting of a mix of a young stellar population (log(xi_ion/[Hz/erg])~25.7) plus an AGN component is required. Furthermore an enhanced C/O abundance ratio is needed for metallicities Z=0.1-0.2Zsun and logU=-1.7 to -1.5. Conclusions: The UV diagnostics we propose should serve as an important basis for the interpretation of observations of high-redshift galaxies. [abridged]

We investigate the morphology of the [CII] emission in a sample of "normal" star-forming galaxies at $5<z<7.2$ in relation to their UV (rest-frame) counterpart. We use new ALMA observations of galaxies at $z\sim6-7$ as well as a careful re-analysis of archival ALMA data. In total 29 galaxies were analysed, 21 of which are detected in [CII]. For several of the latter the [CII] emission breaks into multiple components. Only a fraction of these [CII] components, if any, is associated with the primary UV systems, while the bulk of the [CII] emission is associated either with fainter UV components, or not associated with any UV counterpart at the current limits. By taking into account the presence of all these components, we find that the L$_{\rm [CII]}$-SFR relation at early epochs is fully consistent with the local relation, but it has a dispersion of 0.48$\pm$0.07 dex, which is about two times larger than observed locally. We also find that the deviation from the local L$_{\rm [CII]}$-SFR relation has a weak anti-correlation with the EW(Ly$\alpha$). The morphological analysis also reveals that [CII] emission is generally much more extended than the UV emission. As a consequence, these primordial galaxies are characterised by a [CII] surface brightness generally much lower than expected from the local ${\rm \Sigma_{[CII}-\Sigma_{SFR}}$ relation. These properties are likely a consequence of a combination of different effects, namely: gas metallicity, [CII] emission from obscured star-forming regions, strong variations of the ionisation parameter, and circumgalactic gas in accretion or ejected by these primeval galaxies.

We search for passive galaxies at $z>3$ in the GOODS-South field, using different techniques based on photometric data, and paying attention to develop methods that are sensitive to objects that have become passive shortly before the epoch of observation. We use CANDELS HST catalogues, ultra-deep $Ks$ data and new IRAC photometry, performing spectral energy distribution fitting using models with abruptly quenched star formation histories. We then single out galaxies which are best fitted by a passively evolving model, and having only low probability $( < 5 \%)$ star-forming solutions. We verify the effects of including nebular lines emission, and we consider possible solutions at different redshifts. The number of selected sources dramatically depends on the models used in the SED fitting. Without including emission lines and with photometric redshifts fixed at the CANDELS estimate, we single out 30 candidates; the inclusion of nebular lines emission reduces the sample to 10 objects; allowing for solutions at different redshifts, only 2 galaxies survive as robust candidates. Most of the candidates are not far-infrared emitters, corroborating their association with passive galaxies. Our results translate into an upper limit in the number density of $\sim 0.173$ arcmin$^2$ above the detection limit. However, we conclude that the selection of passive galaxies at $z>3$ is still subject to significant uncertainties, being sensitive to assumptions in the SED modeling adopted and to the relatively low S/N of the objects. By means of dedicated simulations, we show that JWST will greatly enhance the accuracy, allowing for a much more robust classification.

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

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

The fraction of Lyman-$\alpha$ emitters among the galaxy population has been found to increase from $z\sim0$ to $z\sim6$ and drop dramatically at $z>6$. This drop has been interpreted as an effect of an increasingly neutral intergalactic medium with increasing redshift, while a LyC escape fraction evolving with redshift. We report the result of a large VLT/FORS2 program aiming to confirm spectroscopically a large galaxy sample at $z\geq6$ that has been selected in several independent fields through the Lyman Break technique. Combining those data with archival data, we create a large and homogeneous sample of $z\sim6$ galaxies ($N=127$), complete in terms of Ly$\alpha$ detection at $>95\%$ for EW(Ly$\alpha)\geq25\AA$. We use this sample to derive a new measurement of the LAE fraction at $z\sim6$ and derive the physical properties of these galaxies through spectral energy distribution fitting. We find a median LAE fraction at $z\sim6$ lower than in previous studies. The main difference between LAEs and non-LAEs is that the latter are significantly dustier. Using predictions of our SED fitting code accounting for nebular emission, we find an effective Ly$\alpha$ escape fraction $f^{eff}_{esc}(Ly\alpha)=0.23^{+0.36}_{-0.17}$ remarkably consistent with the value derived by comparing UV luminosity function with Ly$\alpha$ luminosity function. We conclude that the drop in the LAE fraction from $z\sim6$ to $z>6$ is less dramatic than previously found and the effect of an increasing IGM neutral fraction is possibly observed at $5<z<6$. Based on our derived $f^{eff}_{esc}(Ly\alpha)$, we find that the IGM has a relatively small impact on Ly$\alpha$ photon visibility at $z\sim6$, with a lower limit for the IGM transmission to \lya\ photons, $T_{IGM}\gtrsim0.20$, likely due to the presence of outflows. [abdridged]

Fast and energetic winds are invoked by galaxy formation models as essential processes in the evolution of galaxies. These outflows can be powered either by star-formation and/or AGN activity, but the relative dominance of the two mechanisms is still under debate. We use spectroscopic stacking analysis to study the properties of the low-ionization phase of the outflow in a sample of 1330 star-forming galaxies (SFGs) and 79 X-ray detected (42<log(L_X)<45 erg/s) Type 2 AGN at 1.7<z<4.6 selected from a compilation of deep optical spectroscopic surveys, mostly zCOSMOS-Deep and VUDS. We measure mean velocity offsets of -150 km/s in the SFGs while in the AGN sample the velocity is much higher (-950 km/s), suggesting that the AGN is boosting the outflow up to velocities that could not be reached only with the star- formation contribution. The sample of X-ray AGN has on average a lower SFR than non-AGN SFGs of similar mass: this, combined with the enhanced outflow velocity in AGN hosts, is consistent with AGN feedback in action. We further divide our sample of AGN into two X-ray luminosity bins: we measure the same velocity offsets in both stacked spectra, at odds with results reported for the highly ionized phase in local AGN, suggesting that the two phases of the outflow may be mixed only up to relatively low velocities, while the highest velocities can be reached only by the highly ionized phase.

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

We report the discovery of 11 very faint (r< 23), low surface brightness ({\mu}_r< 27 mag/arcsec^2) dwarf galaxies in one deep field in the Virgo cluster, obtained by the prime focus cameras (LBC) at the Large Binocular Telescope (LBT). These extend our previous sample to reach a total number of 27 galaxies in a field of just of 0.17 deg^2 located at a median distance of 390 kpc from the cluster center. Their association with the Virgo cluster is supported by their separate position in the central surface brightness - total magnitude plane with respect to the background galaxies of similar total magnitude. For a significant fraction (26\%) of the sample the association to the cluster is confirmed by spectroscopic follow-up. We show that the mere abundance of satellite galaxies corresponding to our observed number in the target field provides extremely tight constraints on Dark Matter models with suppressed power spectrum compared to the Cold Dark Matter case, independently of the galaxy luminosity distribution. In particular, requiring the observed number of satellite galaxies not to exceed the predicted abundance of Dark Matter sub-halos yields a limit m_X >3 keV at 1-{\sigma} and m_X > 2.3 keV at 2-{\sigma} confidence level for the mass of thermal Warm Dark Matter particles. Such a limit is competitive with other limits set by the abundance of ultra-faint satellite galaxies in the Milky Way, is completely independent of baryon physics involved in galaxy formation, and has the potentiality for appreciable improvements with next observations. We extend our analysis to Dark Matter models based on sterile neutrinos, showing that our observations set tight constraints on the combination of sterile neutrino mass m_{\nu} and mixing parameter sin^2(2{\theta}). We discuss the robustness of our results with respect to systematics.

We present the multiwavelength photometry of two Frontier Fields massive galaxy clusters MACS-J0717 and MACS-J1149 and their parallel fields, ranging from HST to ground based K and Spitzer IRAC bands, and the public release of photometric redshifts and rest frame properties of galaxies found in cluster and parallel pointings. This work was done within ASTRODEEP project and aims to provide a reference for future investigations of the extragalactic populations. To fully exploit the depth of the images and detect faint sources we used an accurate procedure which carefully removes the foreground light of bright cluster sources and the intra-cluster light thus enabling detection and measurement of accurate fluxes in crowded cluster regions. This same procedure has been successfully used to derive the photometric catalogue of MACS-J0416 and Abell-2744. The obtained multi-band photometry was used to derive photometric redshifts, magnification and physical properties of sources. In line with the first two FF catalogues released by ASTRODEEP, the photometric redshifts reach $\sim$4$\%$ accuracy. Moreover we extend the presently available samples to galaxies intrinsically as faint as H160$\sim$32-34 mag thanks the magnification factors induced to strong gravitational lensing. Our analysis allows us to probe galaxy masses larger then 10$^{7}$ M$\odot$ and/or SFR=0.1-1M$\odot$/yr out to redshift z$>6$.

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 present new ALMA observations of the [OIII]88$\mu$m line and high angular resolution observations of the [CII]158$\mu$m line in a normal star forming galaxy at z$=$7.1. Previous [CII] observations of this galaxy had detected [CII] emission consistent with the Ly$\alpha$ redshift but spatially slightly offset relative to the optical (UV-rest frame) emission. The new [CII] observations reveal that the [CII] emission is partly clumpy and partly diffuse on scales larger than about 1kpc. [OIII] emission is also detected at high significance, offset relative to the optical counterpart in the same direction as the [CII] clumps, but mostly not overlapping with the bulk of the [CII] emission. The offset between different emission components (optical/UV and different far-IR tracers) is similar to what observed in much more powerful starbursts at high redshift. We show that the [OIII] emitting clump cannot be explained in terms of diffuse gas excited by the UV radiation emitted by the optical galaxy, but it requires excitation by in-situ (slightly dust obscured) star formation, at a rate of about 7 M$_{\odot}$/yr. Within 20 kpc from the optical galaxy the ALMA data reveal two additional [OIII] emitting systems, which must be star forming companions. We discuss that the complex properties revealed by ALMA in the z$\sim$7.1 galaxy are consistent with expectations by recent models and cosmological simulations, in which differential dust extinction, differential excitation and different metal enrichment levels, associated with different subsystems assembling a galaxy, are responsible for the different appearance of the system when observed with different tracers.

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.

Within one billion years of the Big Bang, intergalactic hydrogen was ionized by sources emitting ultraviolet and higher energy photons. This was the final phenomenon to globally affect all the baryons (visible matter) in the Universe. It is referred to as cosmic reionization and is an integral component of cosmology. It is broadly expected that intrinsically faint galaxies were the primary ionizing sources due to their abundance in this epoch. However, at the highest redshifts ($z>7.5$; lookback time 13.1 Gyr), all galaxies with spectroscopic confirmations to date are intrinsically bright and, therefore, not necessarily representative of the general population. Here, we report the unequivocal spectroscopic detection of a low luminosity galaxy at $z>7.5$. We detected the Lyman-$\alpha$ emission line at $\sim 10504$ {\AA} in two separate observations with MOSFIRE on the Keck I Telescope and independently with the Hubble Space Telescope's slit-less grism spectrograph, implying a source redshift of $z = 7.640 \pm 0.001$. The galaxy is gravitationally magnified by the massive galaxy cluster MACS J1423.8+2404 ($z = 0.545$), with an estimated intrinsic luminosity of $M_{AB} = -19.6 \pm 0.2$ mag and a stellar mass of $M_{\star} = 3.0^{+1.5}_{-0.8} \times 10^8$ solar masses. Both are an order of magnitude lower than the four other Lyman-$\alpha$ emitters currently known at $z > 7.5$, making it probably the most distant representative source of reionization found to date.

The CIII] and CIV rest-frame UV emission lines are powerful probes of the ionizations states of galaxies. They have furthermore been suggested as alternatives for spectroscopic redshift confirmation of objects at the epoch of reionization ($z>6$), where the most frequently used redshift indicator, Ly$\alpha$, is attenuated by the high fraction of neutral hydrogen in the inter-galactic medium. However, currently only very few confirmations of carbon UV lines at these high redshifts exist, making it challenging to quantify these claims. Here, we present the detection of CIV$\lambda\lambda$1548,1551\AA\ in \HST\ slitless grism spectroscopy obtained by GLASS of a Ly$\alpha$ emitter at $z=6.11$ multiply imaged by the massive foreground galaxy cluster RXJ2248. The CIV emission is detected at the 3--5$\sigma$ level in two images of the source, with marginal detection in two other images. We do not detect significant CIII]$\lambda\lambda$1907,1909\AA\ emission implying an equivalent width EW$_\textrm{CIII]}<20$\AA\ (1$\sigma$) and $\textrm{CIV/CIII}>0.7$ (2$\sigma$). Combined with limits on the rest-frame UV flux from the HeII$\lambda$1640\AA\ emission line and the OIII]$\lambda\lambda$1661,1666\AA\ doublet, we put constraints on the metallicity and the ionization state of the galaxy. The estimated line ratios and equivalent widths do not support a scenario where an AGN is responsible for ionizing the carbon atoms. SED fits including nebular emission lines imply a source with a mass of log(M/M$_\odot)\sim9$, SFR of around 10M$_\odot$/yr, and a young stellar population $<50$Myr old. The source shows a stronger ionizing radiation field than objects with detected CIV emission at $z<2$ and adds to the growing sample of low-mass (log(M/M$_\odot)\lesssim9$) galaxies at the epoch of reionization with strong radiation fields from star formation.

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27<z<3.40 to faint magnitude limits (L=0.2L*, or equivalently M1500~-19). We have measured through stacks a stringent upper limit (<1.7% at 1sigma) for the relative escape fraction of HI ionizing photons from bright galaxies (L>L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields.

Exploiting the data of the Grism Lens-Amplified Survey from Space (GLASS), we characterize the spatial distribution of star formation in 76 high star forming galaxies in 10 clusters at 0.3< z <0.7. All these galaxies are likely restricted to first infall. In a companion paper we contrast the properties of field and cluster galaxies, whereas here we correlate the properties of H{\alpha} emitters to a number of tracers of the cluster environment to investigate its role in driving galaxy transformations. H{\alpha} emitters are found in the clusters out to 0.5 virial radii, the maximum radius covered by GLASS. The peak of the H{\alpha} emission is offset with respect to the peak of the UV-continuum. We decompose this offsets into a radial and tangential component. The radial compo- nent points away from the cluster center in 60% of the cases, with 95% confidence. The decompositions agree with cosmological simulations, i.e. the H{\alpha} emission offset correlates with galaxy velocity and ram-pressure stripping signatures. Trends between H{\alpha} emitter properties and surface mass density distributions and X-ray emissions emerge only for unrelaxed clusters. The lack of strong correlations with the global environment does not allow us to identify a unique environmental effect originating from the cluster center. In contrast, correla- tions between H{\alpha} morphology and local number density emerge. We conclude that local effects, uncorrelated to the cluster-centric radius, play a more important role in shaping galaxy properties.

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.

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

We present the discovery and spectrophotometric characterization of a large sample of 164 faint ($i_{AB}$ $\sim$ $23$-$25$ mag) star-forming dwarf galaxies (SFDGs) at redshift $0.13$ $\leq z \leq$ $0.88$ selected by the presence of bright optical emission lines in the VIMOS Ultra Deep Survey (VUDS). We investigate their integrated physical properties and ionization conditions, which are used to discuss the low-mass end of the mass-metallicity relation (MZR) and other key scaling relations. We use optical VUDS spectra in the COSMOS, VVDS-02h, and ECDF-S fields, as well as deep multiwavelength photometry, to derive stellar masses, star formation rates (SFR) and gas-phase metallicities. The VUDS SFDGs are compact (median $r_{e}$ $\sim$ $1.2$ kpc), low-mass ($M_{*}$ $\sim$ $10^7-10^9$ $M_{\odot}$) galaxies with a wide range of star formation rates (SFR($H\alpha$) $\sim 10^{-3}-10^{1}$ $M_{\odot}/yr$) and morphologies. Overall, they show a broad range of subsolar metallicities (12+log(O/H)=$7.26$-$8.7$; $0.04$ $\lesssim Z/Z_{\odot} \lesssim$ $1$). The MZR of SFDGs shows a flatter slope compared to previous studies of galaxies in the same mass range and redshift. We find the scatter of the MZR partly explained in the low mass range by varying specific SFRs and gas fractions amongst the galaxies in our sample. Compared with simple chemical evolution models we find that most SFDGs do not follow the predictions of a "closed-box" model, but those from a gas regulating model in which gas flows are considered. While strong stellar feedback may produce large-scale outflows favoring the cessation of vigorous star formation and promoting the removal of metals, younger and more metal-poor dwarfs may have recently accreted large amounts of fresh, very metal-poor gas, that is used to fuel current star formation.

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

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

The largest galaxies in the Universe reside in galaxy clusters. Using sensitive observations of carbon-monoxide, we show that the Spiderweb Galaxy -a massive galaxy in a distant protocluster- is forming from a large reservoir of molecular gas. Most of this molecular gas lies between the protocluster galaxies and has low velocity dispersion, indicating that it is part of an enriched inter-galactic medium. This may constitute the reservoir of gas that fuels the widespread star formation seen in earlier ultraviolet observations of the Spiderweb Galaxy. Our results support the notion that giant galaxies in clusters formed from extended regions of recycled gas at high redshift.

Exploiting the slitless spectroscopy taken as part of the Grism Lens-Amplified Survey from Space (GLASS), we present an extended analysis of the spatial distribution of star formation in 76 galaxies in 10 clusters at 0.3< z <0.7. We use 85 foreground and background galaxies in the same redshift range as a field sample. The samples are well matched in stellar mass (10^8-10^11 M_sun) and star formation rate (0.5-50 M_sun/yr). We visually classify galaxies in terms of broad-band morphology, Halpha morphology and likely physical process acting on the galaxy. Most Halpha emitters have a spiral morphology (41+/-8% in clusters, 51+/-8% in the field), followed by mergers/interactions (28+/-8%, 31+/-7%, respectively) and early-type galaxies (remarkably as high as 29+/-8% in clusters and 15+/-6% in the field). A diversity of Halpha morphologies is detected, suggesting a diversity of physical processes. In clusters, 30+/-8% of the galaxies present a regular morphology, mostly consistent with star formation diffused uniformly across the stellar population (mostly in the disk component, when present). The second most common morphology (28+/-8%) is asymmetric/jellyfish, consistent with ram pressure stripping or other non-gravitational processes in 18+/-8% of the cases. Ram pressure stripping appears significantly less prominent in the field (2+/-2%), where the most common morphology/mechanism appears to be consistent with minor gas rich mergers or clump accretion. This work demonstrates that while environment specific mechanisms affect galaxy evolution at this redshift, they are diverse and their effects subtle. A full understanding of this complexity requires larger samples and detailed and spatially resolved physical models.

We present new results on [CII]158$\mu$ m emission from four galaxies in the reionization epoch. These galaxies were previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly$\alpha$ emission line in their spectra. The Ly$\alpha$ emission line is redshifted by 100-200 km/s compared to the systemic redshift given by the [CII] line. These velocity offsets are smaller than what is observed in z~3 Lyman break galaxies with similar UV luminosities and emission line properties. Smaller velocity shifts reduce the visibility of Ly$\alpha$ and hence somewhat alleviate the need for a very neutral IGM at z~7 to explain the drop in the fraction of Ly$\alpha$ emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0.6-1.6 x10$^8 L_\odot$: these luminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming and metal poor galaxies and also below z =5.5 Lyman break galaxies with similar star formation rates. We argue that previous undetections of [CII] in z~7 galaxies with similar or smaller star formation rates are due to selection effects: previous targets were mostly strong Ly$\alpha$ emitters and therefore probably metal poor systems, while our galaxies are more representative of the general high redshift star forming population .

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

We report the detection of Ly$\alpha$ emission at $\sim9538$\AA{} in the Keck/DEIMOS and \HST WFC3 G102 grism data from a triply-imaged galaxy at $z=6.846\pm0.001$ behind galaxy cluster MACS J2129.4$-$0741. Combining the emission line wavelength with broadband photometry, line ratio upper limits, and lens modeling, we rule out the scenario that this emission line is \oii at $z=1.57$. After accounting for magnification, we calculate the weighted average of the intrinsic Ly$\alpha$ luminosity to be $\sim1.3\times10^{42}~\mathrm{erg}~\mathrm{s}^{-1}$ and Ly$\alpha$ equivalent width to be $74\pm15$\AA{}. Its intrinsic UV absolute magnitude at 1600\AA{} is $-18.6\pm0.2$ mag and stellar mass $(1.5\pm0.3)\times10^{7}~M_{\odot}$, making it one of the faintest (intrinsic $L_{UV}\sim0.14~L_{UV}^*$) galaxies with Ly$\alpha$ detection at $z\sim7$ to date. Its stellar mass is in the typical range for the galaxies thought to dominate the reionization photon budget at $z\gtrsim7$; the inferred Ly$\alpha$ escape fraction is high ($\gtrsim 10$\%), which could be common for sub-$L^*$ $z\gtrsim7$ galaxies with Ly$\alpha$ emission. This galaxy offers a glimpse of the galaxy population that is thought to drive reionization, and it shows that gravitational lensing is an important avenue to probe the sub-$L^*$ galaxy population.

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 multiwavelength photometric catalogues (HST, Spitzer and Hawk-I K band) for the first two of the Frontier Fields, Abell2744 and MACSJ0416 (plus their parallel fields). To detect faint sources even in the central regions of the clusters, we develop a robust and repeatable procedure that uses the public codes Galapagos and Galfit to model and remove most of the light contribution from both the brightest cluster members as well as the ICL. We perform the detection on the HST H160 processed image to obtain a pure H-selected sample. We also add a sample of sources which are undetected in the H160 image but appear on a stacked infrared image. Photometry in the other HST bands is obtained using SExtractor, performed again on residual images after the Galfit procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands has been obtained using our PSF-matching deconfusion code T-PHOT. A similar procedure, but without the need for the foreground light removal, is adopted for the Parallel fields. The procedure allows for the detection and the photometric measurements of ~2500 sources per field. We deliver and release complete photometric H-detected catalogues, with the addition of a complementary sample of infrared-detected sources. All objects have multiwavelength coverage including B to H HST bands, plus K band from Hawk-I, and 3.6 - 4.5 {\mu}m from Spitzer. Full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. The multiwavelength catalogues are publicly available and are ready to be used for scientific purposes. Our procedures allows for the detection of outshined objects near the bright galaxies, which, coupled with the magnification effect of the clusters, can reveal extremely faint high redshift sources. Full analysis on photometric redshifts is presented in a companion Paper II. [abridged]

We present the first public release of photometric redshifts, galaxy rest-frame properties and associated magnification values in the cluster and parallel pointings of the first two Frontier Fields, Abell-2744 and MACS-J0416. We exploit a multi-wavelength catalogue ranging from HST to ground-based K and Spitzer IRAC which is specifically designed to enable detection and measurement of accurate fluxes in crowded cluster regions. The multi-band information is used to derive photometric redshifts and physical properties of sources detected either in the H-band image alone or from a stack of four WFC3 bands. To minimize systematics median photometric redshifts are assembled from six different approaches to photo-z estimates. Their reliability is assessed through a comparison with available spectroscopic samples. State of the art lensing models are used to derive magnification values on an object-by-object basis by taking into account sources positions and redshifts. We show that photometric redshifts reach a remarkable ~3-5% accuracy. After accounting for magnification the H band number counts are found in agreement at bright magnitudes with number counts from the CANDELS fields, while extending the presently available samples to galaxies intrinsically as faint as H160~32-33 thanks to strong gravitational lensing. The Frontier Fields allow to probe the galaxy stellar mass distribution at 0.5-1.5 dex lower masses, depending on magnification, with respect to extragalactic wide fields, including sources at Mstar~ 10^7-10^8 Msun at z>5. Similarly, they allow the detection of objects with intrinsic SFRs>1dex lower than in the CANDELS fields reaching 0.1-1 Msun/yr at z~6-10. [abridged]

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

Star-forming galaxies are considered to be the leading candidate sources that dominate the cosmic reionization at z>7, and the search for analogs at moderate redshift showing Lyman continuum (LyC) leakage is currently a active line of research. We have observed a star-forming galaxy at z=3.2 with Hubble/WFC3 in the F336W filter, corresponding to the 730-890A rest-frame, and detect LyC emission. This galaxy is very compact and also has large Oxygen ratio [OIII]5007/[OII]3727 (>=10). No nuclear activity is revealed from optical/near-infrared spectroscopy and deep multi-band photometry (including the 6Ms X-ray, Chandra). The measured escape fraction of ionizing radiation spans the range 50-100\%, depending on the IGM attenuation. The LyC emission is detected at S/N=10 with m(F336W)=27.57+/-0.11 and it is spatially unresolved, with effective radius R_e<200pc. Predictions from photoionization and radiative transfer models are in line with the properties reported here, indicating that stellar winds and supernova explosions in a nucleated star-forming region can blow cavities generating density-bounded conditions compatible with optically thin media. Irrespective to the nature of the ionizing radiation, spectral signatures of these sources over the entire electromagnetic spectrum are of central importance for their identification during the epoch of reionization, when the LyC is unobservable. Intriguingly, the Spitzer/IRAC photometric signature of intense rest-frame optical emissions ([OIII]+Hbeta) observed recently at z~7.5-8.5 is similar to what is observed in this galaxy. Only the James Webb Space Telescope will measure optical line ratios at z>7 allowing a direct comparison with lower redshift LyC emitters, as reported here.

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.

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

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

[abbreviated] We present a census of Ly\alpha\ emission at $z\gtrsim7$ utilizing deep near infrared HST grism spectroscopy from the first six completed clusters of the Grism Lens-Amplified Survey from Space (GLASS). In 24/159 photometrically selected galaxies we detect emission lines consistent with Ly\alpha\ in the GLASS spectra. Based on the distribution of signal-to-noise ratios and on simulations we expect the completeness and the purity of the sample to be 40-100% and 60-90%, respectively. For the objects without detected emission lines we show that the observed (not corrected for lensing magnification) 1$\sigma$ flux limits reaches $5\times10^{-18}$erg/s/cm$^{2}$ per position angle over the full wavelength range of GLASS (0.8-1.7$\mu$m). Based on the conditional probability of Ly\alpha\ emission measured from the ground at $z\sim7$ we would have expected 12-18 Ly\alpha\ emitters. This is consistent with the number of detections, within the uncertainties, confirming the drop in Ly\alpha\ emission with respect to $z\sim6$. These candidates include a promising source at $z=8.1$. The spatial extent of Ly\alpha\ in a deep stack of the most convincing Ly\alpha\ emitters with $\langle z\rangle=7.2$ is consistent with that of the rest-frame UV continuum. Extended Ly$\alpha$ emission, if present, has a surface brightness below our detection limit, consistent with the properties of lower redshift comparison samples. From the stack we estimate upper limits on rest-frame UV emission line ratios and find $f_\textrm{CIV} / f_\textrm{Ly${\alpha}$} \lesssim 0.32$ and $f_\textrm{CIII]} / f_\textrm{Ly$\alpha$} \lesssim 0.23$ in good agreement with other values published in the literature.

We present the first study of the spatial distribution of star formation in z~0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS0717.5+3745 and MACS1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10^8-10^11 M_sun, and star formation rates in the range 1-20 M_sun/yr. Both in clusters and in the field, H{\alpha} is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside out growth. In ~20% of the cases, the H{\alpha} emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the H{\alpha} emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models and find no conclusive results. The diversity of morphologies and sizes observed in H_alpha illustrates the complexity of the environmental process that regulate star formation. Upcoming analysis of the full GLASS dataset will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial dataset.

We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$\mu$Jy, to derive differential number counts down to 60$\mu$Jy and 100$\mu$Jy at $\lambda=$1.3 mm and $\lambda=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$\mu$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{\lambda = 1.1~mm} = 130~\mu$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{\lambda = 1.3~mm} = 60~\mu$Jy. At the faintest flux limits, i.e. 30$\mu$Jy and 40$\mu$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.

We have discovered 11 ultra-faint ($r\lesssim 22.1$) low surface brightness (LSB, central surface brightness $23\lesssim \mu_r\lesssim 26$) dwarf galaxy candidates in one deep Virgo field of just $576$ arcmin$^2$ obtained by the Large Binocular Camera (LBC) at the Large Binocular Telescope (LBT). Their association with the Virgo cluster is supported by their distinct position in the central surface brightness - total magnitude plane with respect to the background galaxies of similar total magnitude. They have typical absolute magnitudes and scale sizes, if at the distance of Virgo, in the range $-13\lesssim M_r\lesssim -9$ and $250\lesssim r_s\lesssim 850$ pc, respectively. Their colors are consistent with a gradually declining star formation history with a specific star formation rate of the order of $10^{-11}$ yr$^{-1}$, i.e. 10 times lower than that of main sequence star forming galaxies. They are older than the cluster formation age and appear regular in morphology. They represent the faintest extremes of the population of low luminosity LSB dwarfs that has been recently detected in wider surveys of the Virgo cluster. Thanks to the depth of our observations we are able to extend the Virgo luminosity function down to $M_r\sim -9.3$ (corresponding to total masses $M\sim 10^7$ M$_{\odot}$), finding an average faint-end slope $\alpha\simeq -1.4$. This relatively steep slope puts interesting constraints on the nature of the Dark Matter and in particular on warm Dark Matter (WDM) often invoked to solve the overprediction of the dwarf number density by the standard CDM scenario. We derive a lower limit on the WDM particle mass $>1.5$ keV.

[Abridged] We investigate the physical properties of a Lyman continuum emitter candidate at $z=3.212$ with photometric coverage from $U$ to MIPS 24$\mu$m band and VIMOS/VLT and MOSFIRE/Keck spectroscopy. Investigation of the UV spectrum confirms a direct spectroscopic detection of the Lyman continuum emission with $S/N>5$. Non-zero Ly$\alpha$ flux at the systemic redshift and high Lyman-$\alpha$ escape fraction suggest a low HI column density. The weak C and Si low-ionization absorption lines are also consistent with a low covering fraction along the line of sight. The [OIII]$\lambda\lambda4959,5007+\mathrm{H}\beta$ equivalent width is one of the largest reported for a galaxy at $z>3$ ($\mathrm{EW}([\mathrm{OIII}]\lambda\lambda4959,5007+\mathrm{H}\beta) \simeq 1600\AA$, rest-frame) and the NIR spectrum shows that this is mainly due to an extremely strong [OIII] emission. The large observed [OIII]/[OII] ratio ($>10$) and high ionization parameter are consistent with prediction from photoionization models in case of a density-bounded nebula scenario. Furthermore, the $\mathrm{EW}([\mathrm{OIII}]\lambda\lambda4959,5007+\mathrm{H}\beta)$ is comparable to recent measurements reported at $z\sim7-9$, in the reionization epoch. We also investigate the possibility of an AGN contribution to explain the ionizing emission but most of the AGN identification diagnostics suggest that stellar emission dominates instead. This source is currently the first high-$z$ example of a Lyman continuum emitter exhibiting indirect and direct evidences of a Lyman continuum leakage and having physical properties consistent with theoretical expectation from Lyman continuum emission from a density-bounded nebula.

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

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters of galaxies at redshift z=0.308-0.686, including the Hubble Frontier Fields (HFF). The Wide Field Camera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81-1.69mum through grisms G102 and G141, while the Advanced Camera for Surveys in parallel mode provides G800L spectra of the infall regions of the clusters. The WFC3 spectra are taken at two almost orthogonal position angles in order to minimize the effects of confusion. After summarizing the scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data processing pipeline. We then utilize MACSJ0717.5+3745, a HFF cluster and the first one observed by GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than H_AB=23 is visually inspected by at least two co-authors and a redshift is measured when sufficient information is present in the spectra. Furthermore, we conducted a thorough search for emission lines through all the GLASS WFC3 spectra with the aim of measuring redshifts for sources with continuum fainter than H_AB=23. We provide a catalog of 139 emission-line based spectroscopic redshifts for extragalactic sources, including three new redshifts of multiple image systems (one probable, two tentative). In addition to the data itself we also release software tools that are helpful to navigate the data.

We report new deep ALMA observations aimed at investigating the [CII]158um line and continuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 5-15 Msun/yr these systems are much more representative of the high-z galaxy population than other systems targeted in the past by millimeter observations. For the galaxy with the deepest observation we detect [CII] emission at redshift z=7.107, fully consistent with the Lyalpha redshift, but spatially offset by 0.7" (4 kpc) from the optical emission. At the location of the optical emission, tracing both the Lyalpha line and the far-UV continuum, no [CII] emission is detected in any of the three galaxies, with 3sigma upper limits significantly lower than the [CII] emission observed in lower reshift galaxies. These results suggest that molecular clouds in the central parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emission mostly arises from more external accreting/satellite clumps of neutral gas. These findings are in agreement with recent models of galaxy formation. Thermal far-infrared continuum is not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV emission ratio do not exceed those derived in metal- and dust-poor galaxies.

In the low-redshift Universe, the most powerful radio sources are often associated with gas-rich galaxy mergers or interactions. We here present evidence for an advanced, gas-rich (`wet') merger associated with a powerful radio galaxy at a redshift of z~2. This radio galaxy, MRC 0152-209, is the most infrared-luminous high-redshift radio galaxy known in the southern hemisphere. Using the Australia Telescope Compact Array, we obtained high-resolution CO(1-0) data of cold molecular gas, which we complement with HST/WFPC2 imaging and WHT long-slit spectroscopy. We find that, while roughly M(H2) ~ 2 x 10$^{10}$ M$_{\odot}$ of molecular gas coincides with the central host galaxy, another M(H2) ~ 3 x 10$^{10}$ M$_{\odot}$ is spread across a total extent of ~60 kpc. Most of this widespread CO(1-0) appears to follow prominent tidal features visible in the rest-frame near-UV HST/WFPC2 imaging. Ly$\alpha$ emission shows an excess over HeII, but a deficiency over L(IR), which is likely the result of photo-ionisation by enhanced but very obscured star formation that was triggered by the merger. In terms of feedback, the radio source is aligned with widespread CO(1-0) emission, which suggests that there is a physical link between the propagating radio jets and the presence of cold molecular gas on scales of the galaxy's halo. Its optical appearance, combined with the transformational stage at which we witness the evolution of MRC 0152-209, leads us to adopt the name `Dragonfly Galaxy'.

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

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

We aim to investigate the effect of the escaping ionizing radiation on the color selection of high redshift galaxies and identify candidate Lyman continuum (LyC) emitters. The intergalactic medium prescription of Inoue et al.(2014) and galaxy synthesis models of Bruzual&Charlot (2003) have been used to properly treat the ultraviolet stellar emission, the stochasticity of the intergalactic transmission and mean free path in the ionizing regime. Color tracks are computed by turning on/off the escape fraction of ionizing radiation. At variance with recent studies, a careful treatment of IGM transmission leads to no significant effects on the high-redshift broad-band color selection. The decreasing mean free path of ionizing photons with increasing redshift further diminishes the contribution of the LyC to broad-band colors. We also demonstrate that prominent LyC sources can be selected under suitable conditions by calculating the probability of a null escaping ionizing radiation. The method is applied to a sample of galaxies extracted from the GOODS-S field. A known LyC source at z=3.795 is successfully recovered as a LyC emitter candidate and another convincing candidate at z=3.212 is reported. A detailed analysis of the two sources (including their variability and morphology) suggests a possible mixture of stellar and non-stellar (AGN) contribution in the ultraviolet. Conclusions: Classical broad-band color selection of 2.5<z<4.5 galaxies does not prevent the inclusion of LyC emitters in the selected samples. Large fesc in relatively bright galaxies (L>0.1L*) could be favored by the presence of a faint AGN not easily detected at any wavelength. A hybrid stellar and non-stellar (AGN) ionizing emission could coexist in these systems and explain the tensions found among the UV excess and the stellar population synthesis models reported in literature.

We present the public release of the stellar mass catalogs for the GOODS-S and UDS fields obtained using some of the deepest near-IR images available, achieved as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) project. We combine the effort from ten different teams, who computed the stellar masses using the same photometry and the same redshifts. Each team adopted their preferred fitting code, assumptions, priors, and parameter grid. The combination of results using the same underlying stellar isochrones reduces the systematics associated with the fitting code and other choices. Thanks to the availability of different estimates, we can test the effect of some specific parameters and assumptions on the stellar mass estimate. The choice of the stellar isochrone library turns out to have the largest effect on the galaxy stellar mass estimates, resulting in the largest distributions around the median value (with a semi interquartile range larger than 0.1 dex). On the other hand, for most galaxies, the stellar mass estimates are relatively insensitive to the different parameterizations of the star formation history. The inclusion of nebular emission in the model spectra does not have a significant impact for the majority of galaxies (less than a factor of 2 for ~80% of the sample). Nevertheless, the stellar mass for the subsample of young galaxies (age < 100 Myr), especially in particular redshift ranges (e.g., 2.2 < z < 2.4, 3.2 < z < 3.6, and 5.5 < z < 6.5), can be seriously overestimated (by up to a factor of 10 for < 20 Myr sources) if nebular contribution is ignored.

In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN candidates at $z>4$ in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwavelength coverage from Chandra, HST, Spitzer and various groundbased telescopes. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR $H$ band down to very faint levels ($H\leq27$) using reliable photometric redshifts. This corresponds at $z>4$ to a selection criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they show X-ray fluxes above a threshold of $F_X\sim 1.5\times 10^{-17}$ cgs (0.5-2 keV). We have found 22 AGN candidates at $z>4$ and we have derived the first estimate of the UV luminosity function in the redshift interval $4<z<6.5$ and absolute magnitude interval $-22.5\lesssim M_{1450} \lesssim -18.5$ typical of local Seyfert galaxies. The faint end of the derived luminosity function is about two/four magnitudes fainter at $z\sim 4-6$ than that derived from previous UV surveys. We have then estimated ionizing emissivities and hydrogen photoionization rates in the same redshift interval under reasonable assumptions and after discussion of possible caveats, the most important being the large uncertainties involved in the estimate of photometric redshift for sources with featureless, almost power-law SEDs and/or low average escape fraction of ionizing photons from the AGN host galaxies. We argue that, under reasonable evaluations of possible biases, the probed AGN population can produce at $z=4-6.5$ photoionization rates consistent with that required to keep highly ionized the intergalactic medium observed in the Lyman-$\alpha$ forest of high redshift QSO spectra, providing an important contribution to the cosmic reionization.

We present spatially resolved gas-phase metallicity for a system of three galaxies at z=1.85 detected in the Grism Lens-Amplified Survey from Space (GLASS). The combination of HST's diffraction limit and strong gravitational lensing by the cluster MACS J0717+3745 results in a spatial resolution of ~200-300 pc, enabling good spatial sampling despite the intrinsically small galaxy sizes. The galaxies in this system are separated by 50-200 kpc in projection and are likely in an early stage of interaction, evidenced by relatively high specific star formation rates. Their gas-phase metallicities are consistent with larger samples at similar redshift, star formation rate, and stellar mass. We obtain a precise measurement of the metallicity gradient for one galaxy and find a shallow slope compared to isolated galaxies at high redshift, consistent with a flattening of the gradient due to gravitational interaction. An alternative explanation for the shallow metallicity gradient and elevated star formation rate is rapid recycling of metal-enriched gas, but we find no evidence for enhanced gas-phase metallicities which should result from this effect. Notably, the measured stellar masses log(M/Msun) = 7.2-9.1 probe to an order of magnitude below previous mass-metallicity studies at this redshift. The lowest mass galaxy has properties similar to those expected for Fornax at this redshift, indicating that GLASS is able to directly study the progenitors of local group dwarf galaxies on spatially resolved scales. Larger samples from the full GLASS survey will be ideal for studying the effects of feedback, and the time evolution of metallicity gradients. These initial results demonstrate the utility of HST spectroscopy combined with gravitational lensing for characterizing resolved physical properties of galaxies at high redshift.

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

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

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.

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

(arXiv abridged abstract) The observed UV rest-frame spectra of distant galaxies are the result of their intrinsic emission combined with absorption along the line of sight produced by the inter-galactic medium (IGM). Here we analyse the evolution of the mean IGM transmission Tr(Ly_alpha) and its dispersion along the line of sight for 2127 galaxies with 2.5<z<5.5 in the VIMOS Ultra Deep Survey (VUDS). We fit model spectra combined with a range of IGM transmission to the galaxy spectra using the spectral fitting algorithm GOSSIP+. We use these fits to derive the mean IGM transmission towards each galaxy for several redshift slices from z=2.5 to z=5.5. We find that the mean IGM transmission defined as Tr(Ly_alpha)=e^{-tau} (with tau the HI optical depth) is 79%, 69%, 59%, 55% and 46% at redshifts 2.75, 3,22, 3.70, 4.23, 4.77, respectively. We compare these results to measurements obtained from quasars lines of sight and find that the IGM transmission towards galaxies is in excellent agreement with quasar values up to redshift z~4. We find tentative evidence for a higher IGM transmission at z>= 4 compared to results from QSOs, but a degeneracy between dust extinction and IGM prevents to draw firm conclusions if the internal dust extinction for star-forming galaxies at z>4 takes a mean value significantly in excess of E(B-V)>0.15. Most importantly, we find a large dispersion of IGM transmission along the lines of sight towards distant galaxies with 68% of the distribution within 10 to 17% of the median value in delta z=0.5 bins, similar to what is found on the LOS towards QSOs. We demonstrate the importance of taking into account this large range of IGM transmission when selecting high redshift galaxies based on their colour properties (e.g. LBG or photometric redshift selection) or otherwise face a significant incompleteness in selecting high redshift galaxy populations.

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

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

We used near-infrared integral field spectroscopic observations from the AMAZE and LSD programs to constrain the metallicity in a sample of 40 star forming galaxies at 3<z<5 (most of which at z~3.4). We measure metallicities by exploiting strong emission line diagnostics. We found that a significant fraction of star-forming galaxies at z~3.4 deviate from the Fundamental Metallicity Relation (FMR), with a metallicity up to a factor of ten lower than expected according to the FMR. This deviation does not correlate with the dynamical properties of the galaxy or with the presence of interactions. To investigate the origin of the metallicity deviations in more detail, we also infer information on the gas content, by inverting the Schmidt-Kennicutt relation. In agreement with recent CO observational data, we found that, in contrast with the steeply rising trend at 0<z<2, the gas fraction in massive galaxies remains constant, with indication of a marginal decline, at 2<z<3.5. When combined with the metallicity information, we infer that to explain both the low metallicity and gas content in z~3.4 galaxies, both prominent outflows and massive pristine gas inflows are needed. In ten galaxies we can also spatially resolve the metallicity distribution. We found that the metallicity generally anticorrelates with the distribution of star formation and with the gas surface density. We discuss these findings in terms of pristine gas inflows towards the center, and outflows of metal rich gas from the center toward the external regions. (Abridged)

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

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

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

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

We aim to determine the redshift of GDS-1408, the most solid z~7 galaxy candidate lying in the Hubble Ultra Deep Field. We have used all the VLT spectra of GDS-1408 collected by us and two other groups with FORS2 at VLT in the last five years, for a total integration time of 52hr. The combined spectrum is the deepest ever obtained of a galaxy in the Reionization epoch. We do not detect any emission line or continuum over the whole wavelength range, up to 10100A. Based on an accurate set of simulations, we are able to put a stringent upper limit of f(Lya)<3x10^(-18) erg/s/cm2 at 3-9 sigma in the explored wavelength range, corresponding to a rest-frame equivalent width EW<9A. Combining this limit with the SED modelling we refine the redshift to be z=6.82+/- 0.1 (1-sigma). The same SED fitting indicates that GDS-1408 is relatively extinct (A1600~1) with a dust corrected star formation rate of ~ 20 Msol/yr. The comparison between the un-attenuated equivalent width predicted by the case-B recombination theory and the observed upper limit, provides a limit on the effective Lya escape fraction of f_(esc)^(eff)(Lya)<8%. Even though we cannot rule out a major contribution of the inter/circum galactic medium in damping the line, a plausible interpretation is that GDS-1408 is moderately evolved and contains sufficient gas and dust to attenuate the Lya emission, before it reaches the intergalactic medium. The redshift confirmation of even the best z~7 candidates is very hard to achieve (unless the Lya or unusually strong rest-UV nebular emission lines are present) with the current generation of 8-10m class telescopes. We show that both JWST and E-ELT will be necessary to make decisive progresses. Currently, the increased redshift accuracy obtained with this kind of analysis makes ALMA an interesting option for the redshift confirmation.

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

The large cross-section of the Lyman alpha (Lya) line makes it a sensitive probe of the ionization state of the intergalactic medium (IGM). Here we present the most complete study to date of the IGM Lya opacity, and its application to the redshift evolution of the 'Lya fraction', i.e. the fraction of color-selected galaxies with a detectable Lya emission line. We use a tiered approach, which combines large-scale semi-numeric simulations of reionization with moderate-scale hydrodynamic simulations of the ionized IGM. This allows us to simultaneously account for evolution in both: (i) the opacity from an incomplete (patchy) reionization, parameterized by the filling factor of ionized regions, Q_HII; and (ii) the opacity from self-shielded systems in the ionized IGM, parameterized by the average photo-ionization rate inside HII regions, \Gamma. In contrast to recent empirical models, attenuation from patchy reionization has a unimodal distribution along different sightlines, while attenuation from self-shielded systems is more bimodal. We quantify the average IGM transmission in our (Q_HII, \Gamma) parameter space, which can easily be used to interpret new data sets. Using current observations, we predict that the Lya fraction cannot drop by more than a factor of ~2 with IGM attenuation alone, even for HII filling factors as low as Q_HII>0.1. Larger changes in the Lya fraction could result from a co-evolution with galaxy properties. Marginalizing over \Gamma, we find that current observations constrain Q_HII < 0.6 at z=7 [68% confidence level (C.L.)]. However, all of our parameter space is consistent with observations at 95% C.L., highlighting the need for larger observational samples at z >= 6.

We investigate various galaxy population properties of the massive X-ray luminous galaxy cluster XDCP J0044.0-2033 at z=1.58, which constitutes the most extreme matter density peak at this redshift currently known. We analyze deep VLT/HAWK-I NIR data in the J- and Ks-bands, complemented by Subaru imaging in i and V, Spitzer observations at 4.5 micron, and new spectroscopic observations with VLT/FORS2. We detect a cluster-associated excess population of about 90 galaxies, which follows a centrally peaked, compact NFW galaxy surface density profile with a concentration of c200~10. Based on the Spitzer 4.5 micron imaging data, we measure a stellar mass fraction of fstar,500=(3.3+-1.4)% consistent with local values. The total J- and Ks-band galaxy luminosity functions of the core region yield characteristic magnitudes J* and Ks* consistent with expectations from simple z_f=3 burst models. However, a detailed look at the morphologies and color distributions of the spectroscopically confirmed members reveals that the most massive galaxies are undergoing a very active mass assembly epoch through merging processes. Consequently, the bright end of the cluster red-sequence is not in place, while at intermediate magnitudes [Ks*,Ks*+1.6] a red-locus population is present, which is then sharply truncated at magnitudes fainter than Ks*+1.6. The dominant cluster core population comprises post-quenched galaxies transitioning towards the red-sequence at intermediate magnitudes, while additionally a significant blue cloud population of faint star-forming galaxies is present even in the densest central regions. Our observations lend support to the scenario in which the dominant effect of the dense z~1.6 cluster environment is an accelerated mass assembly timescale through merging activity that is responsible for driving core galaxies across the mass quenching threshold of log(Mstar/Msun)~10.4.

We analyse 14 LBGs at z~2.8-3.8 constituting the only sample where both a spectroscopic measurement of their metallicity and deep IR observations (CANDELS+HUGS survey) are available. Fixing the metallicity of population synthesis models to the observed values, we determine best-fit physical parameters under different assumptions about the star-formation history and also consider the effect of nebular emission. For comparison we determine the UV slope of the objects, and use it to estimate their SFR_UV99 by correcting the UV luminosity following Meurer et al. (1999). A comparison between SFR obtained through SED-fitting (SFR_fit) and the SFR_UV99 shows that the latter are underestimated by a factor 2-10, regardless of the assumed SFH. Other SFR indicators (radio, far-IR, X-ray, recombination lines) coherently indicate SFRs a factor of 2-4 larger than SFR_UV99 and in closer agreement with SFR_fit. This discrepancy is due to the solar metallicity implied by the usual beta-A1600 conversion factor. We propose a refined relation, appropriate for sub-solar metallicity LBGs: A1600 = 5.32+1.99beta. This relation reconciles the dust-corrected UV with the SED-fitting and the other SFR indicators. We show that the fact that z~3 galaxies have sub-solar metallicity implies an upward revision by a factor of ~1.5-2 of the global SFRD, depending on the assumptions about the age of the stellar populations. We find very young best-fit ages (10-500 Myrs) for all our objects. From a careful examination of the uncertainties in the fit and the amplitude of the Balmer break we conclude that there is little evidence of the presence of old stellar population in at least half of the LBGs in our sample, suggesting that these objects are probably caught during a huge star-formation burst, rather than being the result of a smooth evolution.

We investigate the gas-phase metallicity and Lyman Continuum (LyC) escape fraction of a strongly gravitationally lensed, extreme emission-line galaxy at z=3.417, J1000+0221S, recently discovered by the CANDELS team. We derive ionization and metallicity sensitive emission-line ratios from H+K band LBT/LUCI medium resolution spectroscopy. J1000+0221S shows high ionization conditions, as evidenced by its enhanced [OIII]/[OII] and [OIII]/Hbeta ratios. Consistently, strong-line methods based on the available line ratios suggest that J1000+0221S is an extremely metal-poor galaxy, with a metallicity of 12+log(O/H) < 7.44 (< 5% solar), placing it among the most metal-poor star-forming galaxies at z > 3 discovered so far. In combination with its low stellar mass (2x10^8 Msun) and high star formation rate (5 Msun/yr), the metallicity of J1000+0221S is consistent with the extrapolation to low masses of the mass-metallicity relation traced by Lyman-break galaxies at z > 3, but it is 0.55 dex lower than predicted by the fundamental metallicity relation at z < 2.5. These observations suggest the picture of a rapidly growing galaxy, possibly fed by the massive accretion of pristine gas. Additionally, deep LBT/LBC in the UGR bands are used to derive a limit to the LyC escape fraction, thus allowing us to explore for the first time the regime of sub-L* galaxies at z > 3. We find a 1sigma upper limit to the escape fraction of 23%, which adds a new observational constraint to recent theoretical models predicting that sub-L* galaxies at high-z have high escape fractions and thus are the responsible for the reioization of the Universe.

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

The rapid decline in the number of strong Lyman Alpha (Lya) emitting galaxies at z > 6 provides evidence for neutral hydrogen in the IGM, but is difficult to explain with plausible models for reionization. We demonstrate that the observed reduction in Lya flux from galaxies at z > 6 can be explained by evolution in the escape fraction of ionizing photons, f_esc. We find that the median observed drop in the fraction of galaxies showing strong Lya emission, as well as the observed evolution of the Lya luminosity function both follow from a small increase in f_esc of Delta f_esc ~ 0.1 from f_esc ~ 0.6 at z ~ 6. This high escape fraction may be at odds with current constraints on the ionising photon escape fraction, which favor smaller values of f_esc < 20%. However, models that invoke a redshift evolution of f_ esc that is consistent with these constraints can suppress the z~7 Lya flux to the observed level, if they also include a small evolution in global neutral fraction of Delta x_HI ~ 0.2. Thus, an evolving escape fraction of ionising photons can be a plausible part of the explanation for evolution in the Lya emission of high redshift galaxies. More generally, our analysis also shows that the drop in the Lya fraction is quantitatively consistent with the observed evolution in the Lya luminosity functions of Lya Emitters.

Aims. We present a deep multiwavelength imaging survey (UGR) in 3 different fields, Q0933, Q1623, and COSMOS, for a total area of ~1500arcmin^2. The data were obtained with the Large Binocular Camera on the Large Binocular Telescope. Methods. To select our Lyman break galaxy (LBG) candidates, we adopted the well established and widely used color-selection criterion (U-G vs. G-R). One of the main advantages of our survey is that it has a wider dynamic color range for U-dropout selection than in previous studies. This allows us to fully exploit the depth of our R-band images, obtaining a robust sample with few interlopers. In addition, for 2 of our fields we have spectroscopic redshift information that is needed to better estimate the completeness of our sample and interloper fraction. Results. Our limiting magnitudes reach 27.0(AB) in the R band (5\sigma) and 28.6(AB) in the U band (1\sigma). This dataset was used to derive LBG candidates at z~3. We obtained a catalog with a total of 12264 sources down to the 50% completeness magnitude limit in the R band for each field. We find a surface density of ~3 LBG candidates arcmin^2 down to R=25.5, where completeness is >=95% for all 3 fields. This number is higher than the original studies, but consistent with more recent samples.

We report the LBT/MODS1 spectroscopic confirmation of two images of faint Lyman alpha emitters at $z=6.4$ behind the Frontier Fields galaxy cluster MACSJ0717.5+3745. A wide range of lens models suggests that the two images are highly magnified, with a strong lower limit of mu>5. These are the faintest z>6 candidates spectroscopically confirmed to date. These may be also multiple images of the same z=6.4 source as supported by their similar intrinsic properties, but the lens models are inconclusive regarding this interpretation. To be cautious, we derive the physical properties of each image individually. Thanks to the high magnification, the observed near-infrared (restframe ultraviolet) part of the spectral energy distributions and Ly-alpha lines are well detected with S/N(m_1500)>~10 and S/N(Ly-alpha)~10-15. Adopting mu>5, the absolute magnitudes, M_1500, and Ly-alpha fluxes, are fainter than -18.7 and 2.8x10^(-18)erg/s/cm2, respectively. We find a very steep ultraviolet spectral slope beta=-3.0+/-0.5 (F_lambda=lambda^(beta)), implying that these are very young, dust-free and low metallicity objects, made of standard stellar populations or even extremely metal poor stars (age<~30Myr, E(B-V)=0 and metallicity 0.0-0.2 Z/Zsolar). The objects are compact (< 1 kpc^(2)), and with a stellar mass M* < 10^(8) M_solar. The very steep beta, the presence of the Ly-alpha line and the intrinsic FWHM (<300 km/s) of these newborn objects do not exclude a possible leakage of ionizing radiation. We discuss the possibility that such faint galaxies may resemble those responsible for cosmic reionization.

The Grism Lens-Amplified Survey from Space (GLASS) is a Hubble Space Telescope (HST) Large Program, which will obtain 140 orbits of grism spectroscopy of the core and infall regions of 10 galaxy clusters, selected to be among the very best cosmic telescopes. Extensive HST imaging is available from many sources including the CLASH and Frontier Fields programs. We introduce the survey by analyzing spectra of faint multiply-imaged galaxies and $z\gtrsim6$ galaxy candidates obtained from the first seven orbits out of fourteen targeting the core of the Frontier Fields cluster MACS0717.5+3745. Using the G102 and G141 grisms to cover the wavelength range 0.8-1.7$\mu$m, we confirm 4 strongly lensed systems by detecting emission lines in each of the images. For the 9 $z\gtrsim6$ galaxy candidates clear from contamination, we do not detect any emission lines down to a seven-orbit 1$\sigma$ noise level of $\sim$5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. Taking lensing magnification into account, our flux sensitivity reaches $\sim$0.2-5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. These limits over an uninterrupted wavelength range rule out the possibility that the high-$z$ galaxy candidates are instead strong line emitters at lower redshift. These results show that by means of careful modeling of the background - and with the assistance of lensing magnification - interesting flux limits can be reached for large numbers of objects, avoiding pre-selection and the wavelength restrictions inherent to ground-based multi-slit spectroscopy. These observations confirm the power of slitless HST spectroscopy even in fields as crowded as a cluster core.

We have evaluated the diffuse intracluster light (ICL) in the central core of the galaxy cluster CL0024+17 at $z\sim 0.4$ observed with the prime focus camera (LBC) at LBT. The measure required an accurate removal of the galaxies light within $\sim 200$ kpc from the center. The residual background intensity has then been integrated in circular apertures to derive the average ICL intensity profile. The latter shows an approximate exponential decline as expected from theoretical cold dark matter models. The radial profile of the ICL over the galaxies intensity ratio (ICL fraction) is increasing with decreasing radius but near the cluster center it starts to bend and then decreases where the overlap of the halos of the brightest cluster galaxies becomes dominant. Theoretical expectations in a simplified CDM scenario show that the ICL fraction profile can be estimated from the stripped over galaxy stellar mass ratio in the cluster. It is possible to show that the latter quantity is almost independent of the properties of the individual host galaxies but mainly depends on the average cluster properties. The predicted ICL fraction profile is thus very sensitive to the assumed CDM profile, total mass and concentration parameter of the cluster. Adopting values very similar to those derived from the most recent lensing analysis in CL0024+17 we find a good agreement with the observed ICL fraction profile. The galaxy counts in the cluster core have then been compared with that derived from composite cluster samples in larger volumes, up to the clusters virial radius. The galaxy counts in the CL0024+17 core appear flatter and the amount of bending respect to the average cluster galaxy counts imply a loss of total emissivity in broad agreement with the measured ICL fraction.

We present a CO(1-0) survey for cold molecular gas in a representative sample of 13 high-z radio galaxies (HzRGs) at 1.4<z<2.8, using the Australia Telescope Compact Array. We detect CO(1-0) emission associated with five sources: MRC 0114-211, MRC 0152-209, MRC 0156-252, MRC 1138-262 and MRC 2048-272. The CO(1-0) luminosities are in the range $L'_{\rm CO} \sim (5 - 9) \times 10^{10}$ K km/s pc$^{2}$. For MRC 0152-209 and MRC 1138-262 part of the CO(1-0) emission coincides with the radio galaxy, while part is spread on scales of tens of kpc and likely associated with galaxy mergers. The molecular gas mass derived for these two systems is M$_{\rm H2} \sim 6 \times 10^{10}\, {\rm M}_{\odot}$ (M$_{\rm H2}$/$L'_{\rm CO}$=0.8). For the remaining three CO-detected sources, the CO(1-0) emission is located in the halo (~50-kpc) environment. These three HzRGs are among the fainter far-IR emitters in our sample, suggesting that similar reservoirs of cold molecular halo gas may have been missed in earlier studies due to pre-selection of IR-bright sources. In all three cases the CO(1-0) is aligned along the radio axis and found beyond the brightest radio hot-spot, in a region devoid of 4.5$\mu$m emission in Spitzer imaging. The CO(1-0) profiles are broad, with velocity widths of ~ 1000 - 3600 km/s. We discuss several possible scenarios to explain these halo reservoirs of CO(1-0). Following these results, we complement our CO(1-0) study with detections of extended CO from the literature and find at marginal statistical significance (95% level) that CO in HzRGs is preferentially aligned towards the radio jet axis. For the eight sources in which we do not detect CO(1-0), we set realistic upper limits of $L'_{\rm CO} \sim 3-4 \times 10^{10}$ K km/s pc$^{2}$. Our survey reveals a CO(1-0) detection rate of 38%, allowing us to compare the CO(1-0) content of HzRGs with that of other types of high-z galaxies.

We present the spectroscopic confirmation of a structure of galaxies surrounding the radio galaxy MRC0156-252 at z = 2.02. The structure was initially discovered as an overdensity of both near-infrared selected z > 1.6 and mid-infrared selected z > 1.2 galaxy candidates. We used the VLT/FORS2 multi-object spectrograph to target ~80 high-redshift galaxy candidates, and obtain robust spectroscopic redshifts for more than half the targets. The majority of the confirmed sources are star-forming galaxies at z > 1.5. In addition to the radio galaxy, two of its close-by companions (< 6'') also show AGN signatures. Ten sources, including the radio galaxy, lie within |z - 2.020 | < 0.015 (i.e., velocity offsets < 1500 km/s) and within projected 2 Mpc comoving of the radio galaxy. Additional evidence suggests not only that the galaxy structure associated with MRC0156-252 is a forming galaxy cluster but also that this structure is most probably embedded in a larger scale structure.

We present multicolour Hubble Space Telescope images of the powerful z=2.4 radio galaxy MRC 0406-244 and model its complex morphology with several components including a host galaxy, a point source, and extended nebular and continuum emission. We suggest that the main progenitor of this radio galaxy was a normal, albeit massive (M ~10^{11} solar masses), star-forming galaxy. The optical stellar disc of the host galaxy is smooth and well described by a S\'ersic profile, which argues against a recent major merger, however there is also a point-source component which may be the remnant of a minor merger. The half-light radius of the optical disc is constrained to lie in the range 3.5 to 8.2kpc, which is of similar size to coeval star forming galaxies. Biconical shells of nebular emission and UV-bright continuum extend out from the host galaxy along the radio jet axis, which is also the minor axis of the host galaxy. The origin of the continuum emission is uncertain, but it is most likely to be young stars or dust-scattered light from the AGN, and it is possible that stars are forming from this material at a rate of 200^{+1420}_{-110} solar masses per year.

We present the results of observations taken with the X-shooter spectrograph devoted to the study of quasars at z ~ 6. This paper focuses on the properties of metals at high redshift traced, in particular, by the C IV doublet absorption systems. Six objects were observed with resolutions ~27 and 34 km/s in the visual, and 37.5 and 53.5 km/s in the near infrared. We detected 102 C IV lines in the range: 4.35 < z < 6.2 of which 27 are above z ~ 5. Thanks to the characteristics of resolution and spectral coverage of X-shooter, we could also detect 25 Si IV doublets associated with the C IV at z>5. The column density distribution function of the C IV line sample is observed to evolve in redshift for z>~ 5.3, with respect to the normalization defined by low redshift (1.5 < z <4) C IV lines. This behaviour is reflected in the redshift evolution of the C IV cosmic mass density, Omega_CIV, of lines with column density in the range 13.4 < log N(C IV) < 15, which is consistent with a drop of a factor of ~2 for z>~ 5.3. Considering only the stronger C IV lines (13.8 < log N(CIV) < 15), Omega_CIV gently rises by a factor of ~10 between z ~ 6.2 and z~ 1.5 with a possible flattening toward z ~ 0. The increase is well fitted by a power law: Omega_CIV = (2+-1)x10^{-8} [(1+z)/4]^{-3.1+-0.1}. An insight into the properties of the C IV absorbers and their evolution with redshift is obtained by comparing the observed column densities of associated C IV, Si IV and C II absorptions with the output of a set of CLOUDY photo-ionization models. As already claimed by cosmological simulations, we find that C IV is a good tracer of the metallicity in the low-density intergalacitc medium gas at z ~ 5-6 while at z ~ 3 it arises in gas with over density delta ~100.

We present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources.

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

The workhorse instruments of the 8-10m class observatories have become their multi-object spectrographs (MOS), providing comprehensive follow-up to both ground-based and space-borne imaging. With the advent of deeper imaging surveys from, e.g., the HST and VISTA, there are a plethora of spectroscopic targets which are already beyond the sensitivity limits of current facilities. This wealth of targets will grow even more rapidly in the coming years, e.g., after the completion of ALMA, the launch of the JWST and Euclid, and the advent of the LSST. Thus, one of the key requirements underlying plans for the next generation of ground-based telescopes, the Extremely Large Telescopes (ELTs), is for even greater sensitivity for optical and infrared spectroscopy. Here we revisit the scientific motivation for a MOS capability on the European ELT, combining updated elements of science cases advanced from the Phase A instrument studies with new science cases which draw on the latest results and discoveries. These science cases address key questions related to galaxy evolution over cosmic time, from studies of resolved stellar populations in nearby galaxies out to observations of the most distant galaxies, and are used to identify the top-level requirements on an 'E-ELT/MOS'. We argue that several of the most compelling ELT science cases demand MOS observations, in highly competitive areas of modern astronomy. Recent technical studies have demonstrated that important issues related to e.g. sky subtraction and multi-object AO can be solved, making fast- track development of a MOS instrument feasible. To ensure that ESO retains world leadership in exploring the most distant objects in the Universe, galaxy evolution and stellar populations, we are convinced that a MOS should have high priority in the instrumentation plan for the E-ELT.

Determining the AGN content in structures of different mass/velocity dispersion and comparing them to higher mass/lower redshift analogs is important to understand how the AGN formation process is related to environmental properties. We use our well-tested cluster finding algorithm to identify structures in the GOODS North and South fields, exploiting the available spectroscopic redshifts and accurate photometric redshifts. We identify 9 structures in GOODS-south (presented in a previous paper) and 8 new structures in GOODS-north. We only consider structures where at least 2/3 of the members brighter than M_R=-20 have a spectroscopic redshift. For those group members that coincide with X-ray sources in the 4 and 2 Msec Chandra source catalogs respectively, we determine if the X-ray emission originates from AGN activity or it is related to the galaxies' star-formation activity. We find that the fraction of AGN with Log L_H > 42 erg/s in galaxies with M_R < -20 is on average 6.3+-1.3%, much higher than in lower redshift groups of similar mass and more than double the fraction found in massive clusters at a similarly high redshift. We then explore the spatial distribution of AGN in the structures and find that they preferentially populate the outer regions. The colors of AGN host galaxies in structures tend to be confined to the green valley, thus avoiding the blue cloud and, partially, also the red-sequence, contrary to what happens in the field. We finally compare our results to the predictions of two sets of semi analytic models to investigate the evolution of AGN and evaluate potential triggering and fueling mechanisms. The outcome of this comparison attests the importance of galaxy encounters, not necessarily leading to mergers, as an efficient AGN triggering mechanism. (abridged)

The high-redshift radio galaxy MRC 1138-262 (`Spiderweb Galaxy'; z = 2.16), is one of the most massive systems in the early Universe and surrounded by a dense `web' of proto-cluster galaxies. Using the Australia Telescope Compact Array, we detected CO(1-0) emission from cold molecular gas -- the raw ingredient for star formation -- across the Spiderweb Galaxy. We infer a molecular gas mass of M(H2) = 6x10^10 M(sun) (for M(H2)/L'(CO)=0.8). While the bulk of the molecular gas coincides with the central radio galaxy, there are indications that a substantial fraction of this gas is associated with satellite galaxies or spread across the inter-galactic medium on scales of tens of kpc. In addition, we tentatively detect CO(1-0) in the star-forming proto-cluster galaxy HAE 229, 250 kpc to the west. Our observations are consistent with the fact that the Spiderweb Galaxy is building up its stellar mass through a massive burst of widespread star formation. At maximum star formation efficiency, the molecular gas will be able to sustain the current star formation rate (SFR ~ 1400 M(sun)/yr, as traced by Seymour et al.) for about 40 Myr. This is similar to the estimated typical lifetime of a major starburst event in infra-red luminous merger systems.

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

The exploration of the relation between galaxy sizes and other physical parameters has provided important clues for understanding galaxy formation. We use the CANDELS Deep+Wide surveys in the GOODS-South, UDS and EGS fields, complemented by data from the HUDF09 program, to address the relation between size and luminosity at z\sim7. We select 153 z-band drop-out galaxies in six different fields characterized by a wide combination of depth and areal coverage, ideally suited to sample without biases the observed size-magnitude plane. Detailed simulations allow us to derive the completeness as a function of size and magnitude and to quantify measurements errors/biases. We find that the half light radius distribution function of z\sim7 galaxies fainter than J=26.6 is peaked at <0.1 arcsec (or equivalently 0.5 kpc proper), while at brighter magnitudes high-z galaxies are typically larger than ~0.15 arcsec. We also find a well defined size-luminosity relation, Rh\simL^1/2. We compute the Luminosity Function in the HUDF and P12HUDF fields, finding large spatial variation on the number density of faint galaxies. Adopting the size distribution and the size-luminosity relation found for faint galaxies at z=7, we derive a mean slope of -1.7\pm0.1 for the luminosity function of LBGs at this redshift. Using this LF, we find that the amount of ionizing photons cannot keep the Universe re-ionized if the IGM is clumpy (C_HII>3) and the Lyman continuum escape fraction of high-z LBGs is relatively low (f_esc<0.3). With future CANDELS data, we can put severe limits to the role of galaxies in the reionization of the Universe.

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.

We have evaluated the contribution of the AGN population to the ionization history of the Universe based on a semi-analytic model of galaxy formation and evolution in the CDM cosmological scenario. The model connects the growth of black holes and of the ensuing AGN activity to galaxy interactions. In the model we have included a self consistent physical description of the escape of ionizing UV photons; this is based on the blast-wave model for the AGN feedback we developed in a previous paper to explain the distribution of hydrogen column densities in AGNs of various redshifts and luminosities, due to absorption by the host galaxy gas. The model predicts UV luminosity functions for AGNs which are in good agreement with those derived from the observations especially at low and intermediate redshifts (z=3). At higher redshifts (z>5) the model tends to overestimate the data at faint luminosities. Critical biases both in the data and in the model are discussed to explain such apparent discrepancies. The predicted hydrogen photoionization rate as a function of redshift is found to be consistent with that derived from the observations. All that suggests to reconsider the role of the AGNs as the main driver of the ionization history of the Universe.

We plan to analyse dust extinction in Lyman Break Galaxies (LBGs) by introducing a new and more reliable approach to their selection and to the characterization of their distribution of UV slopes beta, using deep IR images from HST. We exploit deep WFC3 IR observations of the ERS and HUDF fields over GOODS-South, combined with HST-ACS optical data, to select z~4 LBGs through a new (B-V) vs. (V-H) colour diagram. The UV slope of the selected galaxies is robustly determined by a linear fit over their observed I, Z, Y, J magnitudes, coherently with the original definition of beta. The same fit is used to determine their rest-frame UV magnitudes M1600 through a simple interpolation. We estimate the effect of observational uncertainties with detailed simulations that we also exploit, under a parametric maximum-likelihood approach, to constrain the probability density function of UV slopes PDF(beta) as a function of rest-frame magnitude. We find 142 and 25 robust LBGs in the ERS and HUDF fields respectively, limiting our sample to S/N(H)>10 objects. Our newly defined criteria improve the selection of z~4 LBGs and allow us to exclude red interlopers at lower redshift, especially z~3-3.5 objects. We find that z~4 LBGs are characterized by blue UV slopes, suggesting a low dust extinction: all L<L* galaxies have an average UV slope <beta>\simeq-2.1, while brighter objects only are slightly redder (<beta> -1.9). We find an intrinsic dispersion ~ 0.3 for PDF(beta) at all magnitudes. The SFRD at z~4 corrected according to these estimates turns out to be lower than previously found: log(SFRD)\simeq-1.09 M_sun/yr/Mpc^3. Finally, we discuss how the UV slope of z~4 galaxies changes as a function of the dust-corrected UV magnitude (i.e. SFR) showing that most galaxies with a high SFR (> 80 M_sun/yr) are highly extincted objects. [Abridged]

We use the spatially-resolved, multi-band photometry in the GOODS South field acquired by the CANDELS project to constrain the nature of candidate Lyman continuum (LyC) emitters at redshift z~3.7 identified using ultra-deep imaging below the Lyman limit (1-sigma limit of ~30 AB in a 2" diameter aperture). In 18 candidates, out of a sample of 19 with flux detected at >3-sigma level, the light centroid of the candidate LyC emission is offset from that of the LBG by up to 1.5". We fit the SED of the LyC candidates to spectral population synthesis models to measure photometric redshifts and the stellar population parameters. We also discuss the differences in the UV colors between the LBG and the LyC candidates, and how to estimate the escape fraction of ionizing radiation (f_esc) in cases, like in most of our galaxies, where the LyC emission is spatially offset from the host galaxy. In all but one case we conclude that the candidate LyC emission is most likely due to lower redshift interlopers. Based on these findings, we argue that the majority of similar measurements reported in the literature need further investigation before it can be firmly concluded that LyC emission is detected. Our only surviving LyC candidate is a LBG at z=3.795, which shows the bluest (B-V) color among LBGs at similar redshift, a stellar mass of M~2 x 10^9 Msun, weak interstellar absorption lines and a flat UV spectral slope with no Lya in emission. We estimate its f_esc to be in the range 25%-100%, depending on the dust and intergalactic attenuation.

We used Early Release Science (ERS) observations taken with the Wide Field Camera 3 (WFC3) in the GOODS-S field to study the galaxy stellar mass function (GSMF) at 0.6<=z<4.5. Deep WFC3 near-IR data (for Y as faint as 27.3, J and H as faint as 27.4 AB mag at 5 sigma), as well as deep Ks (as faint as 25.5 at 5 sigma) Hawk-I band data, provide an exquisite data set with which determine in an unprecedented way the low-mass end of the GSMF, allowing an accurate probe of masses as low as M~7.6 10^9 Msun at z~3. Although the area used is relatively small (~33 arcmin^2), we found generally good agreement with previous studies on the entire mass range. Our results show that the slope of the faint-end increases with redshift, from alpha=-1.44+/-0.03 at z~0.8 to alpha=-1.86+/-0.16 at z~3, although indications exist that it does not steepen further between z~3 and z~4. This result is insensitive to any uncertainty in the M* parameter. The steepness of the GSMF faint-end solves the well-known disagreement between the stellar mass density (SMD) and the integrated star formation history at z>2. However, we confirm the that there appears to be an excess of integrated star formation with respect to the SMD at z<2, by a factor of ~2-3. Our comparison of the observations with theoretical predictions shows that the models forecast a greater abundance of low mass galaxies, at least up to z~3, as well as a dearth of massive galaxies at z~4 with respect to the data, and that the predicted SMD is generally overestimated at z<~2.

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.

Aims. Ly-alpha emitters (LAEs) can be detected out to very high redshifts during the epoch of reionization. The evolution of the LAE luminosity function with redshift is a direct probe of the Ly-alpha transmission of the intergalactic medium (IGM), and therefore of the IGM neutral-hydrogen fraction. Measuring the Ly-alpha luminosity function (LF) of LAEs at redshift z = 7.7 therefore allows us to constrain the ionizing state of the Universe at this redshift. Methods. We observed three 7.5'x7.5' fields with the HAWK-I instrument at the VLT with a narrow band filter centred at 1.06 $\mu$m and targeting LAEs at redshift z ~ 7.7. The fields were chosen for the availability of multiwavelength data. One field is a galaxy cluster, the Bullet Cluster, which allowed us to use gravitational amplification to probe luminosities that are fainter than in the field. The two other fields are subareas of the GOODS Chandra Deep Field South and CFHTLS-D4 deep field. We selected z=7.7 LAE candidates from a variety of colour criteria, in particular from the absence of detection in the optical bands. Results. We do not find any LAE candidates at z = 7.7 in ~2.4 x 10^4 Mpc^3 down to a narrow band AB magnitude of ~ 26, which allows us to infer robust constraints on the Ly-alpha LAE luminosity function at this redshift. Conclusions. The predicted mean number of objects at z = 6.5, derived from somewhat different LFs of Hu et al. (2010), Ouchi et al. (2010), and Kashikawa et al. (2011) are 2.5, 13.7, and 11.6, respectively. Depending on which of these LFs we refer to, we exclude a scenario with no evolution from z = 6.5 to z = 7.7 at 85% confidence without requiring a strong change in the IGM Ly-alpha transmission, or at 99% confidence with a significant quenching of the IGM Ly-alpha transmission, possibly from a strong increase in the high neutral-hydrogen fraction between these two redshifts.

We present the final results from our ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 "z--band dropout" candidates selected from our VLT/Hawk-I imaging survey over three independent fields. In particular we report on two newly discovered galaxies at redshift ~6.7 in the NTT deep field: both galaxies show a Ly-alpha emission line with rest-frame EWs of the order 15-20 A and luminosities of 2-4 X 10^{42} erg/s. We also present the results of ultra-deep observations of a sample of i-dropout galaxies, from which we set a solid upper limit on the fraction of interlopers. Out of the 20 z-dropouts observed we confirm 5 galaxies at 6.7 < z < 7.1. This is systematically below the expectations drawn on the basis of lower redshift observations: in particular there is a significant lack of objects with intermediate Ly-alpha EWs (between 20 and 55 A). We conclude that the trend for the fraction of Ly-alpha emission in LBGs that is constantly increasing from z~3 to z~6 is most probably reversed from z~6 to z~7. Explaining the observed rapid change in the LAE fraction among the drop-out population with reionization requires a fast evolution of the neutral fraction of hydrogen in the Universe. Assuming that the Universe is completely ionized at z=6 and adopting the semi-analytical models of Dijkstra et al. (2011), we find that our data require a change of the neutral hydrogen fraction of the order Delta chi_{HI} ~ 0.6 in a time Delta z ~ 1, provided that the escape fraction does not increase dramatically over the same redshift interval.

We present the first results obtained using a tunable narrowband filter in the search for high-z protoclusters. Using the recently commissioned red tunable filter on the Gran Telescopio Canarias we have searched for Lya emitters in a 75 arcmin^2 field centered on the z=4.413 radio galaxy 6C0140+326. With three different wavelength tunings we find a total of 27 unique candidate Lya emitters. The availability of three different wavelength tunings allows us to make estimates of the redshifts for each of the objects. It also allows us to separate a possible protocluster from structure in the immediate foreground. This division shows that the foreground region contains significantly fewer Lya emitters. Also, the spatial distribution of the objects in the protocluster field deviates from a random distribution at the 2.5 sigma level. The observed redshift distribution of the emitters is different from the expected distribution of a blank field at the ~3 sigma level, with the Lya emitters concentrated near the radio galaxy at z>4.38. The 6C0140+326 field is denser by a factor of 9+/-5 than a blank field, and the number density of Lya emitters close to the radio galaxy is similar to that of the z~4.1 protocluster around TNJ1338-1942. We thus conclude that there is an overdensity of Lya emitters around the radio galaxy 6C0140+326. This is one of few known overdensities at such a high redshift.

[Abridged] We use the large public spectroscopic database available in the GOODS-South field to estimate the dynamical mass and the virialization status of cluster ClG 0332-2747 at z=0.734. Cluster members selected from their photometric redshift are used with spectroscopic ones to analyse the galaxy population of the cluster. In the newly released Chandra 4Ms observations we detect a faint extended X-ray emission associated to the cluster. Finally, we compare the optical and X-ray properties of ClG 0332-2747 with the predictions of a well tested semianalytic model. We estimate the velocity dispersion and the virial mass considering all 44 spectroscopic members, or 20 red-sequence members only. We obtain sigma_v=634 +/- 105 Km/s, M_200=3.07 ^{+1.57}_{-1.16}~10^{14} M_sun in the former case, and slightly lower values in the latter case. The cluster appears to have reached the virial equilibrium: it shows a perfectly Gaussian velocity distribution and no evidence for substructures. ClG 0332-2747 contains a high fraction of bright red galaxies, and is dominated by a very massive (1.1 x 10^{12} M_sun) old brightest cluster galaxy (BCG), suggesting that it formed at an early epoch. We detect a faint extended X-ray source centered on the BCG, with a total X-ray luminosity of L_X ~ 2 x 10^{42} erg s^-1 (0.1-2.4 keV). This L_X is lower by a factor of ~10-20 than expected according to the M-L_X relation. We provide a possible explanation of this discrepancy as due to the effects of AGN feedback on the ICM: the semianalytic model reproduces the M-L_X relation measured from "X-ray bright" clusters, and it predicts a high scatter at low masses due to heating and expulsion of the cluster gas. Interestingly, in the model clusters with an evolved galaxy population like ClG 0332-2747 present the largest scatter in X-ray luminosity.

We present an upper limit for the relative escape fraction (f_{esc}^{rel}) of ionizing radiation at z~3.3 using a sample of 11 Lyman Break Galaxies (LBGs) with deep imaging in the U band obtained with the Large Binocular Camera, mounted on the prime focus of the Large Binocular Telescope. We selected 11 LBGs with secure redshift in the range 3.27<z<3.35, from 3 independent fields. We stacked the images of our sources in the R and U band, which correspond to an effective rest-frame wavelength of 1500\AA and 900\AA respectively, obtaining a limit in the U band image of >=30.7(AB)mag at 1 sigma. We derive a 1 sigma upper limit of f_{esc}^{rel}~5%, which is one of the lowest values found in the literature so far at z~3.3. Assuming that the upper limit for the escape fraction that we derived from our sample holds for all galaxies at this redshift, the hydrogen ionization rate that we obtain (Gamma_{-12}<0.3 s^{-1}) is not enough to keep the IGM ionized and a substantial contribution to the UV background by faint AGNs is required. Since our sample is clearly still limited in size, larger z~3 LBG samples, at similar or even greater depths are necessary to confirm these results on a more firm statistical basis.

The study of the Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at z=7 is important for ascertaining their role in the reionization of the Universe. We perform a detailed and critical analysis of the statistical and systematic errors in the z~7 LF determination: we have assembled a large sample of candidate LBGs at z~7 from different surveys, spanning a large variety of areas and depths. In particular, we have combined data from the deep (J<27.4) and ultradeep (J<29.2) surveys recently acquired with the new WFC3 NIR camera on HST, over the GOODS-ERS and the HUDF fields, with ground based surveys in wide and shallow areas from VLT and Subaru. We have used public ACS images in the z-band to select z-dropout galaxies, and other public data both in the blue (BVI) and in the red bands to reject possible low-redshift interlopers. We have compared our results with extensive simulations to quantify the observational effects of our selection criteria as well as the effects of photometric scatter, color selections or the morphology of the candidates. We have found that the number density of faint LBGs at z~7 is only marginally sensitive to the color selection adopted, but it is strongly dependent from the assumption made on the half light distributions of the simulated galaxies, used to correct the observed sample for incompleteness. The slope of the faint end of the LBGs LF has thus a rather large uncertainty, due to the unknown distribution of physical sizes of the z~7 LBGs. We conclude that galaxies at z~7 are unable to reionize the Universe unless there is a significant evolution in the clumpiness of the IGM or in the escape fraction of ionising photons or, alternatively, there is a large population of z~7 LBGs with large physical dimensions but still not detected by the present observations.

The environment of the high-z radio galaxy PKS 1138-262 at z~2.2 is a prime example of a forming galaxy cluster. We use deep SINFONI data to perform a detailed study of the kinematics of the galaxies within 60 kpc of the radio core and we link this to the kinematics of the protocluster on the megaparsec scale. Identification of optical emission lines shows that 11 galaxies are at the redshift of the protocluster. The density of line emitters is more than an order of magnitude higher in the core of the protocluster than the larger scale environment. This implies a matter overdensity in the core of delta_m~70 which is similar to the outskirts of local galaxy clusters. The velocity distribution of the confirmed satellite galaxies shows a broad, double-peaked velocity structure with sigma=1360+/-206 km/s. A similar broad, double-peaked distribution was found in a previous study targeting the large scale protocluster structure, indicating that a common process is acting on both small and large scales. Including all spectroscopically confirmed protocluster galaxies, a velocity dispersion of 1013+/-87 km/s is found. We show that the protocluster has likely decoupled from the Hubble flow and is a dynamically evolved structure. Comparison to the Millenium simulation indicates that the protocluster velocity distribution is consistent with that of the most massive haloes at z~2, but we rule out that the protocluster is a fully virialized structure based on dynamical arguments and its X-ray luminosity. Comparison to merging haloes in the Millennium simulation shows that the structure as observed in and around the Spiderweb galaxy is best interpreted as being the result of a merger between two massive haloes. We propose that this merger can result in an increase in star formation and AGN activity in the protocluster core and is possibly an important stage in the evolution of massive cD galaxies.

We have studied the large scale distribution of matter in the Chandra Deep Field South on the basis of photometric redshifts and we have identifed several over-densities between redshift 0.6 and 2.3. We analyse two of these structures using the deepest X-ray observations ever obtained: 4 Ms with the Chandra satellite and 2.5 Ms with XMM-Newton. We set a very faint upper limit on the X-ray luminosity of a structure at redshift 1.6, and we find an extended X-ray emission from a structure at redshift 0.96 of which we can estimate the gas temperature and make a comparison with the scaling relations between the X-ray luminosity and mass or temperature of high redshift galaxy clusters.

This is a study of H-alpha emitters in two dense galaxy proto-clusters surrounding radio galaxies at z~2. We show that the proto-cluster surrounding MRC 1138-262 contains 14+/-2 times more H-alpha candidates than the average field (9 sigma significance), and the z=2.35 radio galaxy 4C+10.48 is surrounded by 12+/-2 times more emitters than the field (5 sigma), so it is also likely to reside in a dense proto-cluster environment. We compared these H-alpha emitters, situated in dense environments, to a control field sample selected from 3 separate fields forming a total area of 172 arcmin^2. We constructed and compared H-alpha and rest-frame R continuum luminosity functions of the emitters in both environments. The star formation density is on average 13 times greater in the proto-clusters than the field at z~2, so the total star formation rate within the central 1.5Mpc of the proto-clusters exceeds 3000Msun/yr. However, we found no significant difference in the shape of the H-alpha luminosity functions, implying that environment does not substantially affect the strength of the H-alpha line from strongly star forming galaxies. The proto-cluster emitters are typically 0.8mag brighter in rest-frame R continuum than field emitters, implying they are twice as massive as their field counterparts at the same redshift. We also show the proto-cluster galaxies have lower specific star formation rates than field galaxies, meaning the emitters in the dense environments formed more of their stars earlier than the field galaxies. We conclude that galaxy growth in the early Universe was accelerated in dense environments, and that cluster galaxies differed from field galaxies even before the cluster had fully formed.

We report the spectroscopic confirmation of two Lyman break galaxies at redshift > 7. The galaxies were observed as part of an ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 ``z--band dropout'' candidates selected from our VLT/Hawk-I imaging survey. Both galaxies show a prominent emission line at 9735A and 9858A respectively: the lines have fluxes of ~ 1.6-1.2 x 10^(-17) erg/s/cm2 and exhibit a sharp decline on the blue side and a tail on the red side. The asymmetry is quantitatively comparable to the observed asymmetry in z ~ 6 Ly-alpha lines, where absorption by neutral hydrogen in the IGM truncates the blue side of the emission line profile. We carefully evaluate the possibility that the galaxies are instead at lower redshift and we are observing either [OII], [OIII] or H-alpha emission: however from the spectroscopic and the photometric data we conclude that there are no other plausible identifications, except for Ly-alpha at redshift >7, implying that these are two of the most robust redshift determination for galaxies in the reionization epoch. Based on their redshifts and broad--band photometry, we derive limits on the star formation rate and on the ultraviolet spectral slopes of the two galaxies. We argue that these two galaxies alone are unlikely to have ionized the IGM in their surroundings.

X-shooter, with its characteristics of resolution, spectral coverage and efficiency, provides a unique opportunity to obtain spectra of the highest-redshift quasars (z ~ 6) that will allow us to carry out successful investigations on key cosmological issues, from the details of the re-ionization process, to the evolution of the first galaxies and AGNs. In this paper, we present the spectra of three z ~ 6 quasars: one obtained during the commissioning of X-shooter and two in the context of our ongoing GTO programme. Combining this sample with data in the literature, we update the value of the C IV cosmic mass density in the range 4.5 < z < 5, confirming the constant trend with redshift between 2.5 and 5.

We present ultradeep optical spectroscopy obtained with FORS2 on VLT of seven Lyman-break galaxy (LBG) candidates at z>6.5 selected in the GOODS-S field from Hawk-I/VLT and WFC3/HST imaging. For one galaxy we detect a low significance emission line (S/N< 7), located at 9691.5 +/- 0.5A and with flux 3.4 x 10^(-18)erg/cm^2/s. If identified as Lyman alpha, it places the LBG at redshift z=6.972+/- 0.002, with a rest-frame equivalent width EW}=13A. Using Monte Carlo simulations and conservative EW distribution functions at 2<z<6, we estimate that the probability of observing no galaxies in our data with S/N>10 is ~ 2%, and that of observing only one galaxy out of seven with S/N=5 is ~4%, but these can be as small as ~1E-3, depending on the details of the EW distribution. We conclude that either a significant fraction of the candidates is not at high redshift or that some physical mechanism quenches the Lyman alpha emission emerging from the galaxies at z>6.5, abruptly reversing the trend of the increasing fraction of strong emitters with increasing redshift observed up to z~ 6.5. We discuss the possibility that an increasingly neutral intergalactic medium is responsible for such quenching.

(Abridged) We present here the second half of an ESO Large Programme, which exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained - 30 observing hours with Hawk-I in the Y-band of two high galactic latitude fields. We combined the Y-band data with deep J and K Hawk-I observations, and with FORS1/FORS2 U, B, V, R, I, and Z observations to select z-drop galaxies having Z - Y > 1, no optical detection and flat Y - J and Y - K colour terms. We detect 8 high-quality candidates in the magnitude range Y = 25.5 - 26.5 that we add to the z-drop candidates selected in two Hawk-I pointings over the GOODS-South field. We use this full sample of 15 objects found in -161 arcmin^2 of our survey to constrain the average physical properties and the evolution of the number density of z ~ 7 LBGs. A stacking analysis yields a best-fit SED with photometric redshift z= 6.85 +0.20 -0.15 and an E(B-V)=0.05 +0.15 -0.05. We compute a binned estimate of the z ~ 7 LF and explore the effects of photometric scatter and model uncertainties on the statistical constraints. After accounting for the expected incompleteness through MonteCarlo simulations, we strengthen our previous finding that a Schechter luminosity function constant from z=6 to z=7 is ruled out at a >99% confidence level, even including the effects of cosmic variance. For galaxies brighter than M_1500= -19.0, we derive a luminosity density rho_UV = 1.5^{+2.1}{-0.8} x 10^25 erg/s/Hz/Mpc^3, implying a decrease by a factor 3.5 from z=6 to z=6.8. We find that, under standard assumptions, the emission rate of ionizing photons coming from UV bright galaxies is lower by at least a factor of two than the value required for reionization. Finally, we exploit deep Hawk-I J and K band observations to derive an upper limit on the number density of M1500<~ -22.0 LBGs at z-8 (Y-dropouts).

Aims: We investigate the nature of the source GDS J033218.92-275302.7at redshift ~ 5.56. Methods: The SED of the source is well sampled by 16 bands photometry, from UV-optical, near infrared and mid-infrared (MID-IR).The detection of signal in the MID-IR Spitzer/IRAC bands 5.8, 8.0 um -- where the nebular emission contribution is less effective -- suggests the presence of a Balmer break, signature of an underlying stellar population formed at earlier epochs. The optical spectrum shows a clear Lya emission line together with semi-forbidden NIV] 1483.3-1486.5 also in emission. Results: From the SED fitting and the Lya modelling it turns out that the source seems to have an evolved component with stellar mass of ~5 x10^(10) Msolar and age ~ 0.4 Gyrs, and a young component with an age of ~ 0.01 Gyrs and SFR in the range of 30-200 Msolar yr^(-1). The limits on the effective radius derived from the ACS/z850 and VLT/Ks bands indicate that this galaxy is denser than the local ones with similar mass. A relatively high nebular gas column density is favored from the Lya line modelling (NHI>=10^(21) cm^(-2)). A vigorous outflow (~ 450 km/s) has been measured from the optical spectrum,consistent with the Lya modelling. From ACS observations it turns out that the region emitting Lya photons is spatially compact and of the same order of the effective radius estimated at the ~1400A rest-frame wavelength, whose emission is dominated by the stellar continuum and/or AGN. The gas is blown out from the central region,but given the mass of the galaxy it is uncertain whether it will pollute the IGM to large distances. We argue that a burst of star formation in a dense gas environment is active (possibly containing hot and massive stars and/or a low luminosity AGN), superimposed to an already formed fraction of stellar mass (abridged).

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

We present a population study of several types of galaxies within the protocluster surrounding the radio galaxy MRC0316-257 at z~3.1. In addition to the known population of Ly_alpha emitters (LAEs) and [OIII] emitters, we use colour selection techniques to identify protocluster candidates that are Lyman break galaxies (LBG) and Balmer break galaxies (BBGs). The radio galaxy field contains an excess of LBG candidates, with a surface density 1.6\pm0.3 times larger than found for comparable blank fields. This surface overdensity corresponds to an LBG volume overdensity of ~8\pm4. The BBG photometric redshift distribution peaks at the protocluster's redshift, but we detect no significant surface overdensity of BBG. This is not surprising because a volume overdensity similar to the LBGs would have resulted in a surface density of ~1.2 that found in the blank field. This could not have been detected in our sample. Masses and star formation rates of the candidate protocluster galaxies are determined using SED fitting. These properties are not significantly different from those of field galaxies. The galaxies with the highest masses and star formation rates are located near the radio galaxy, indicating that the protocluster environment influences galaxy evolution at z~3. We conclude that the protocluster around MRC0316-257 is still in the early stages of formation.

We investigate the physical and morphological properties of LBGs at z ~2.5 to ~3.5, to determine if and how they depend on the nature and strength of the Lyalpha emission. We selected U-dropout galaxies from the z-detected GOODS MUSIC catalog, by adapting the classical Lyman Break criteria on the GOODS filter set. We kept only those galaxies with spectroscopic confirmation, mainly from VIMOS and FORS public observations. Using the full multi-wavelength 14-bands photometry, we determined the physical properties of the galaxies, through a standard spectral energy distribution fitting with the updated Charlot & Bruzual (2009) templates. We also added other relevant observations, i.e. the 24mu m observations from Spitzer/MIPS and the 2 MSec Chandra X-ray observations. Finally, using non parametric diagnostics (Gini, Concentration, Asymmetry, M_20 and ellipticity), we characterized the rest-frame UV morphology of the galaxies. We then analyzed how these physical and morphological properties correlate with the presence of the Lyalpha line in the optical spectra. We find that, unlike at higher redshift, the dependence of physical properties on the Lyalpha line is milder: galaxies without Lyalpha in emission tend to be more massive and dustier than the rest of the sample, but all other parameters, ages, SFRs, X-ray emission as well as UV morphology do not depend strongly on the presence of the line emission. A simple scenario where all LBGs have intrinsically high Lyalpha emission, but where dust and neutral hydrogen content (which shape the final appearance of the Lyalpha) depend on the mass of the galaxies, is able to reproduce the majority of the observed properties at z~3. Some modification might be needed to account for the observed evolution of these properties with cosmic epoch, which is also discussed.

We perform a deep search for galaxies in the redshift range 6.5<z<7.5, to measure the evolution of the number density of luminous galaxies in this redshift range and derive useful constraints on the evolution of their Luminosity Function. We present here the first results of an ESO Large Program, that exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained two Hawk-I pointings on the GOODS South field for a total of 32 observing hours, covering ~90 arcmin2. The images reach Y=26.7 mags for the two fields. We have used public ACS images in the z band to select z-dropout galaxies with the colour criteria Z-Y>1, Y-J<1.5 and Y-K<2. The other public data in the UBVRIJHK bands are used to reject possible low redshift interlopers. The output has been compared with extensive Monte Carlo simulations to quantify the observational effects of our selection criteria as well as the effects of photometric errors. We detect 7 high quality candidates in the magnitude range Y=25.5-26.7. This interval samples the critical range for M* at z>6 (M_1500 ~- 19.5 to -21.5). After accounting for the expected incompleteness, we rule out at a 99% confidence level a Luminosity Function constant from z=6 to z=7, even including the effects of cosmic variance. For galaxies brighter than M_1500=-19.0 we derive a luminosity density rho_UV=1.5^{+2.0}_{-0.9} 10^25 erg/s/Hz/Mpc3, implying a decrease by a factor 3.5 from z=6 to z~6.8. On the basis of our findings, we make predictions for the surface densities expected in future surveys surveys, based on ULTRA-VISTA, HST-WFC3 or JWST-NIRCam, evaluating the best observational strategy to maximise their impact.

We study the star formation and the mass assembly process of 0.3<=z<2.5 galaxies using their IR emission from MIPS 24um band. We used an updated version of the GOODS-MUSIC catalog, extended by the addition of mid-IR fluxes. We compared two different estimators of the Star Formation Rate: the total infrared emission derived from 24um, estimated using both synthetic and empirical IR templates, and the multiwavelength fit to the full galaxy SED. For both estimates, we computed the SFR Density and the Specific SFR. The two SFR tracers are roughly consistent, given the uncertainties involved. However, they show a systematic trend, IR-based estimates exceeding the fit-based ones as the SFR increases. We show that: a) at z>0.3, the SFR is well correlated with stellar mass, and this relationship seems to steepen with redshift (using IR-based SFRs); b) the contribution to the global SFRD by massive galaxies increases with redshift up to ~2.5, more rapidly than for galaxies of lower mass, but appears to flatten at higher z; c) despite this increase, the most important contributors to the SFRD at any z are galaxies of about, or immediately lower than, the characteristic stellar mass; d) at z~2, massive galaxies are actively star-forming, with a median SFR 300 Msun/yr. During this epoch, they assemble a substantial part of their final stellar mass; e) the SSFR shows a clear bimodal distribution. The analysis of the SFRD and the SSFR seems to support the downsizing scenario, according to which high mass galaxies have formed their stars earlier and faster than their low mass counterparts. A comparison with theoretical models indicates that they follow the global increase in the SSFR with redshift and predict the existence of quiescent galaxies even at z>1.5, but they systematically underpredict the average SSFR.

Aims: The aim of this work is to collect a complete, mass--selected sample of galaxies with very low specific star formation rate, for a comparison with the prediction of recent theoretical models. Method: We use the 24/K flux ratio, complemented by the SED fitting to the full 0.35-8.0 mum spectral distribution, to select quiescent galaxies from z~0.4 to z~4 in the GOODS--MUSIC sample. Our observational selection can be translated into thresholds on the specific star formation rate SFR/M_*, that can be used to compare with the theoretical predictions. Results: We find that, in the framework of the well known global decline of the quiescent fraction with redshift, a non-negligible fraction ~15-20% of massive galaxies with very low specific star formation rate exists up to z~4, including a tail of "Red&Dead" galaxies with SFR/M_*<10^{-11}/yr. Recent theoretical models vary to a large extent in the prediction of the fraction of galaxies with very low specific star formation rates, and are unable to provide a global match to our data.

The aim of this paper is to identify and study the properties and galactic content of groups and clusters in the GOODS-South field up to z\sim2.5, and to analyse the physical properties of galaxies as a continuous function of environmental density up to high redshift. We use the deep (z850\sim26), multi-wavelength GOODS-MUSIC catalogue, which has a 15% of spectroscopic redshifts and accurate photometric redshifts for the remaining fraction. On these data, we apply a (2+1)D algorithm, previously developed by our group, that provides an adaptive estimate of the 3D density field. We support our analysis with simulations to evaluate the purity and the completeness of the cluster catalogue produced by our algorithm. We find several high density peaks embedded in larger structures in the redshift range 0.4-2.5. From the analysis of their physical properties (mass profile, M200, \sigmav, LX, U-B vs. B diagram), we derive that most of them are groups of galaxies, while two are poor clusters with masses of few times 10^14 Mo. For these two clusters we find, from the Chandra 2Ms data, an X-ray emission significantly lower than expected from their optical properties, suggesting that the two clusters are either not virialised or gas poor. We also analyse the dependance on environment of galaxy colours, luminosities, stellar masses, ages and star formations. We find that galaxies in high density regions are, on average, more luminous and massive than field galaxies up to z\sim 2. The fraction of red galaxies increases with luminosity and with density up to z\sim 1.2. At higher z this dependance on density disappears. The variation of galaxy properties as a function of redshift and density suggests that a significant change occurs at z\sim 1.5-2. (abridged)

We discuss the uncertainties and the systematic effects that exist in the estimates of the stellar masses of high redshift galaxies, using broad band photometry, and how they affect the deduced galaxy stellar mass function. We use at this purpose the latest version of the GOODS-MUSIC catalog. In particular, we discuss the impact of different synthetic models, of the assumed initial mass function and of the selection band. Using Charlot & Bruzual 2007 and Maraston 2005 models we find masses lower than those obtained from Bruzual & Charlot 2003 models. In addition, we find a slight trend as a function of the mass itself comparing these two mass determinations with that from Bruzual & Charlot 2003 models. As consequence, the derived galaxy stellar mass functions show diverse shapes, and their slope depends on the assumed models. Despite these differences, the overall results and scenario remains unchanged. The masses obtained with the assumption of the Chabrier initial mass function are in average 0.24 dex lower than those from the Salpeter assumption, at all redshifts, causing a shift of galaxy stellar mass function of the same amount. Finally, using a 4.5 um-selected sample instead of a Ks-selected one, we add a new population of highly absorbed, dusty galaxies at z\simeq 2-3 of relatively low masses, yielding stronger constraints on the slope of the galaxy stellar mass function at lower masses.

A significant fraction of high redshift starburst galaxies presents strong Ly alpha emission. Understanding the nature of these galaxies is important to assess the role they played in the early Universe and to shed light on the relation between the narrow band selected Lyalpha emitters and the Lyman break galaxies: are the Lyalpha emitters a subset of the general LBG population? or do they represent the youngest galaxies in their early phases of formation? We studied a sample of UV continuum selected galaxies from z~2.5 to z~6 (U, B, V and i-dropouts) from the GOODS-South survey, that have been observed spectroscopically. Using the GOODS-MUSIC catalog we investigated their physical properties, such as total masses, ages, SFRs, extinction etc as determined from a spectrophotometric fit to the multi-wavelength (U band to mid-IR) SEDs, and their dependence on the emission line characteristics. In particular we determined the nature of the LBGs with Lyalpha in emission and compared them to the properties of narrow band selected Lyalpha emitters. For U and B-dropouts we also compared the properties of LBGs with and without the Lyalpha emission line.

We have analyzed a sample of LBGs from z =3.5 to z=6 selected from the GOODS-S field as B,V and i-dropouts, and with spectroscopic observations showing that they have the Lyalpha line in emission. Our main aim is to investigate their physical properties and their dependence on the emission line characteristics, to shed light on the relation between galaxies with Lyalpha emission and the general LBG population.The objects were selected from their continuum colors and then spectroscopically confirmed by the GOODS collaboration and other campaigns. From the spectra we derived the line flux and EW. We then used U-band to mid-IR photometry from GOODS-MUSIC to derive the physical properties of the galaxies, such as total stellar mass, age and so on, through standard SED fitting techniques.Although most galaxies are fit by young stellar populations, a small but non negligible fraction has SEDs that require considerably older stellar component, up to 1 Gyr. There is no apparent relation between age and EW: some of the oldest galaxies have large EW, and should be also selected in narrow band surveys. Therefore not all Lyalpha emitters are primeval galaxies in the very early stages of formation,as is commonly assumed. We also find a large range of stellar populations, with masses from 5x10^8 Msol to 5x10^10 Msol and SFR from few to 60 Msol/yr. Although there is no correlation between mass and EW, we find a significant lack of massive galaxies with high EW, which could be explained if the most massive galaxies were more dusty and/or contained more neutral gas than less massive objects. Finally we find that more than half of the galaxies contain small but non negligible amounts of dust: the mean E(B-V) and the EW are well correlated, although with a large scatter, as already found at lower redshift

We present initial results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z>3 by means of deep near-IR spectroscopy. Gas metallicities are measured, for an initial sample of nine star forming galaxies at z~3.5, by means of optical nebular lines redshifted into the near-IR. Stellar masses are accurately determined by using Spitzer-IRAC data, which sample the rest-frame near-IR stellar light in these distant galaxies. When compared with previous surveys, the mass-metallicity relation inferred at z~3.5 shows an evolution much stronger than observed at lower redshifts. The evolution is prominent even in massive galaxies, indicating that z~3 is an epoch of major action in terms of star formation and metal enrichment also for massive systems. There are also indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.5 is difficult to reconcile with the predictions of some hierarchical evolutionary models. Such discrepancies suggest that at z>3 galaxies are assembled mostly with relatively un-evolved sub-units, i.e. small galaxies with low star formation efficiency. The bulk of the star formation and metallicity evolution probably occurs once small galaxies are already assembled into bigger systems.

Venemans et al. (2005) found evidence for an overdensity of Ly-alpha emission line galaxies associated with the radio galaxy MRC 0316-257 at z=3.13 indicating the presence of a massive protocluster. Here, we present the results of a search for additional star-forming galaxies and AGN within the protocluster. Narrow-band infrared imaging was used to select candidate [O III] emitters in a 1.1 x 1.1 Mpc^2 region around the radio galaxy. Thirteen candidates have been detected. Four of these are among the previously confirmed sample of Ly-alpha galaxies, and an additional three have been confirmed through follow-up infrared spectroscopy. The three newly confirmed objects lie within a few hundred km/s of each other, but are blueshifted with respect to the radio galaxy and Ly-alpha emitters by ~2100 km/s. Although the sample is currently small, our results indicate that the radio-selected protocluster is forming at the centre of a larger, ~60 co-moving Mpc super-structure. On the basis of an HST/ACS imaging study we calculate dust-corrected star-formation rates and investigate morphologies and sizes of the [O III] candidate emitters. From a comparison of the star formation rate derived from UV-continuum and [O III] emission, we conclude that at least two of the [O III] galaxies harbour an AGN which ionized the O+ gas.

The X-ray spectrum of NGC7213 is known to present no evidence for Compton reflection, a unique result among bright Seyfert 1s. The observed neutral iron K$\alpha$ line, therefore, cannot be associated with a Compton-thick material, like the disc or the torus, but is due to Compton-thin gas, with the Broad Line Region (BLR) as the most likely candidate. To check this hypothesis, a long Chandra HETG observation, together with a quasi-simultaneous optical spectroscopic observation at the ESO NTT EMMI were performed. We found that the iron line is resolved with a FWHM=$2 400^{+1 100}_{-600}$ km/s, in perfect agreement with the value measured for the broad component of the H$\alpha$, $2640^{+110}_{-90}$ km/s. Therefore, NGC7213 is the only Seyfert 1 galaxy whose iron K$\alpha$ line is unambiguously produced in the BLR. We also confirmed the presence of two ionised iron lines and studied them in greater detail than before. The resonant line is the dominant component in the Fe XXV triplet, therefore suggesting an origin in collisionally ionised gas. If this is the case, the blueshift of around 1000 km/s of the two ionised iron lines could be the first measure of the velocity of a starburst wind from its X-ray emission.

We present new constraints on the evolution of the early-type galaxy color-magnitude relation (CMR) based on deep near-infrared imaging of a galaxy protocluster at z=2.16 obtained using NICMOS on-board the Hubble Space Telescope. This field contains a spectroscopically confirmed space-overdensity of Lyman-alpha and H-alpha emitting galaxies which surrounds the powerful radio galaxy MRC 1138-262. Using these NICMOS data we identify a significant surface-overdensity (= 6.2x) of red J-H galaxies in the color-magnitude diagram (when compared with deep NICMOS imaging from the HDF-N and UDF). The optical-NIR colors of these prospective red-sequence galaxies indicate the presence of on-going dust-obscured star-formation or recently formed (<~ 1.5 Gyr)stellar populations in a majority of the red galaxies. We measure the slope and intrinsic scatter of the CMR for three different red galaxy samples selected by a wide color cut, and using photometric redshifts both with and without restrictions on rest-frame optical morphology. In all three cases both the rest-frame $U-B$ slope and intrinsic color scatter are considerably higher than corresponding values for lower redshift galaxy clusters. These results suggest that while some relatively quiescent galaxies do exist in this protocluster both the majority of the galaxy population and hence the color-magnitude relation are still in the process of forming, as expected.

We present a catalogue of overdensities in the GOODS-South field. We find various high density peaks that are embedded in structures diffused on the entire field, up to z ~ 2.5. The slope of their colour-magnitude relation does not show significative evolution with z. We find evidence that galaxies forming these structures are more massive than galaxies located in low density regions. We also analyse the variation of galaxy properties with the associated environmental density and we find that the segregation of red galaxies with density is stronger at low redshift and at high luminosities while it gets much weaker for increasing z.

We present preliminary results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z~3 by means of deep near-IR spectroscopy. Gas metallicities and stellar masses are measured for an initial sample of nine star forming galaxies at z~3.3. When compared with previous surveys, the mass-metallicity relation inferred at z~3.3 shows an evolution significantly stronger than observed at lower redshifts. There are also some indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered as the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.3 is difficult to reconcile with the predictions of some hierarchical evolutionary models. We shortly discuss the possible implications of such discrepancies.

We study the evolution of the galaxy population up to z\sim3 as a function of its colour properties. In particular, luminosity functions and luminosity densities have been derived as a function of redshift for the blue/late and red/early populations. We use data from the GOODS-MUSIC catalogue which have typical magnitude limits z<26 and Ks<23.5 for most of the sample. About 8% of the galaxies have spectroscopic redshifts; the remaining have well calibrated photometric redshifts derived from the extremely wide multi-wavelength coverage in 14 bands (from the U band to the Spitzer 8 \mu m band). We have derived a catalogue of galaxies complete in rest-frame B-band, which has been divided in two subsamples according to their rest-frame U-V colour (or derived specific star formation rate, SSFR) properties. We confirm a bimodality in the U-V colour and SSFR of the galaxy sample up to z\sim 3. This bimodality is used to compute the LFs of the blue/late and red/early subsamples. The LFs of the blue/late and total samples are well represented by steep Schechter functions evolving in luminosity with increasing redshifts. The volume density of the LFs of the red/early populations decreases with increasing redshift. The shape of the red/early LFs shows an excess of faint red dwarfs with respect to the extrapolation of a flat Schechter function and can be represented by the sum of two Schechter functions. Our model for galaxy formation in the hierarchical clustering scenario, which also includes external feedback due to a diffuse UV background, shows a general broad agreement with the LFs of both populations, the larger discrepancies being present at the faint end for the red population. Hints on the nature of the red dwarf population are given on the basis of their stellar mass and spatial distributions.

We report the discovery of a localized overdensity at z~1.6 in the GOODS-South Field, presumably a poor cluster in the process of formation. The three-dimensional galaxy density has been estimated on the basis of well calibrated photometric redshifts from the multiband photometric GOODS-MUSIC catalog using the (2+1)D technique. The density peak is embedded in the larger scale overdensity of galaxies known to exist at z=1.61 in the area. The properties of the member galaxies are compared to those of the surrounding field and we found that the two populations are significantly different supporting the reality of the structure. The reddest galaxies, once evolved according to their best fit models, have colors consistent with the red sequence of lower redshift clusters. The estimated M_200 total mass of the cluster is in the range 1.3 x 10^14 - 5.7x 10^14 Msun, depending on the assumed bias factor b. An upper limit for the 2-10 keV X-ray luminosity, based on the 1Ms Chandra observations, is L_X=0.5 x 10^43 erg s^-1, suggesting that the cluster has not yet reached the virial equilibrium.

We present sensitive near-infrared spectroscopic observations for a sample of five z ~ 6 quasars. These are amongst the most distant, currently known quasars in the universe. The spectra have been obtained using ISAAC at the VLT and include the CIV, MgII and FeII lines. We measure the FeII/MgII line ratio, as an observational proxy for the Fe/alpha element ratio. We derive a ratio of 2.7+/-0.8 for our sample, which is similar to that found for lower redshift quasars, i.e., we provide additional evidence for the lack of evolution in the FeII/MgII line ratio of quasars up to the highest redshifts. This result demonstrates that the sample quasars must have undergone a major episode of iron enrichment in less than one Gyr and star formation must have commenced at z > 8. The linewidths of the MgII and CIV lines give two estimates for the black hole masses. A third estimate is given by assuming that the quasars emit at their Eddington luminosity. The derived masses using these three methods agree well, implying that the quasars are not likely to be strongly lensed. We derive central black hole masses of 0.3-5.2 10^9 solar masses. We use the difference between the redshift of MgII (a proxy for the systemic redshift of the quasar) and the onset of the Gunn Peterson trough to derive the extent of the ionized Stromgren spheres around our target quasars. The derived physical radii are about five Mpc. Using a simple ionization model, the emission of the central quasars would need of order 10^6-10^8 year to create these cavities in a surrounding intergalactic medium with a neutral fraction between 0.1 and 1.0. As the e-folding time scale for the central accreting black hole is on the order of a few times 10^7 year, it can grow by one e-folding or less within this time span.

We have analysed the physical properties of z~4 Lyman Break Galaxies observed in the GOODS-S survey, in order to investigate the possible differences between galaxies where the Ly-alpha is present in emission, and those where the line is absent or in absorption. The objects have been selected from their optical color and then spectroscopically confirmed by Vanzella et al. (2005). From the public spectra we assessed the nature of the Ly-alpha emission and divided the sample into galaxies with Ly-alpha in emission and objects without Ly-alpha line (i.e. either absent or in absorption). We have then used the complete photometry, from U band to mid infrared from the GOODS-MUSIC database, to study the observational properties of the galaxies, such as UV spectral slopes and optical to mid-infrared colors, and the possible differences between the two samples. Finally through standard spectral fitting tecniques we have determined the physical properties of the galaxies, such as total stellar mass, stellar ages and so on, and again we have studied the possible differences between the two samples. Our results indicate that LBG with Ly-alpha in emission are on average a much younger and less massive population than the LBGs without Ly-alpha emission. Both populations are forming stars very actively and are relatively dust free, although those with line emission seem to be even less dusty on average. We briefly discuss these results in the context of recent models for the evolution of Lyman break galaxies and Ly-alpha emitters.

[Abridged] We present the results of a large program conducted with the Very Large Telescope and Keck telescope to search for forming clusters of galaxies near powerful radio galaxies at 2.0 < z < 5.2. We obtained narrow- and broad-band images of nine radio galaxies and their surroundings. The imaging was used to select candidate Lyman alpha emitting galaxies in ~3x3 Mpc^2 areas near the radio galaxies. A total of 337 candidate emitters were found with a rest-frame Lyman alpha equivalent width of EW_0 > 15 A and Sigma = EW_0/Delta EW_0 > 3. Follow-up spectroscopy confirmed 168 Lyman alpha emitters near eight radio galaxies. The success rate of our selection procedure is 91%. At least six of our eight fields are overdense in Lyman alpha emitters by a factor 3-5. Also, the emitters show significant clustering in velocity space. In the overdense fields, the width of the velocity distributions of the emitters is a factor 2-5 smaller than the width of the narrow-band filters. Taken together, we conclude that we have discovered six forming clusters of galaxies (protoclusters). We estimate that roughly 75% of powerful (L_2.7GHz > 10^33 erg/s/Hz/sr) high redshift radio galaxies reside in a protocluster, with a sizes of at least 1.75 Mpc. We estimate that the protoclusters have masses in the range 2-9 x 10^14 Msun and they are likely to be progenitors of present-day (massive) clusters of galaxies. For the first time, we have been able to estimate the velocity dispersion of cluster progenitors from z~5 to ~2. The velocity dispersion of the emitters increases with cosmic time, in agreement with the dark matter velocity dispersion in numerical simulations of forming massive clusters.

We present a deep image of the radio galaxy MRC 1138-262 taken with the Hubble Space Telescope (HST) at a redshift of z = 2.2. The galaxy is known to have properties of a cD galaxy progenitor and be surrounded by a 3 Mpc-sized structure, identified with a protocluster. The morphology shown on the new deep HST/ACS image is reminiscent of a spider's web. More than 10 individual clumpy features are observed, apparently star-forming satellite galaxies in the process of merging with the progenitor of a dominant cluster galaxy 11 Gyr ago. There is an extended emission component, implying that star formation was occurring over a 50 times 40 kpc region at a rate of more than 100 M_sun/yr. A striking feature of the newly named ``Spiderweb galaxy'' is the presence of several faint linear galaxies within the merging structure. The dense environments and fast galaxy motions at the centres of protoclusters may stimulate the formation of these structures, which dominate the faint resolved galaxy populations in the Hubble Ultra Deep Field. The new image provides a unique testbed for simulations of forming dominant cluster galaxies.

The goal of this work is to measure the evolution of the Galaxy Stellar Mass Function and of the resulting Stellar Mass Density up to redshift ~4, in order to study the assembly of massive galaxies in the high redshift Universe. We have used the GOODS-MUSIC catalog, containing ~3000 Ks-selected galaxies with multi-wavelength coverage extending from the U band to the Spitzer 8 micron band, of which 27% have spectroscopic redshifts and the remaining fraction have accurate photometric redshifts. On this sample we have applied a standard fitting procedure to measure stellar masses. We compute the Galaxy Stellar Mass Function and the resulting Stellar Mass Density up to redshift ~4, taking into proper account the biases and incompleteness effects. Within the well known trend of global decline of the Stellar Mass Density with redshift, we show that the decline of the more massive galaxies may be described by an exponential timescale of ~6 Gyrs up to z~1.5, and proceeds much faster thereafter, with an exponential timescale of ~0.6 Gyrs. We also show that there is some evidence for a differential evolution of the Galaxy Stellar Mass Function, with low mass galaxies evolving faster than more massive ones up to z~1-1.5 and that the Galaxy Stellar Mass Function remains remarkably flat (i.e. with a slope close to the local one) up to z~1-1.3. The observed behaviour of the Galaxy Stellar Mass Function is consistent with a scenario where about 50% of present-day massive galaxies formed at a vigorous rate in the epoch between redshift 4 and 1.5, followed by a milder evolution until the present-day epoch.

We present Spitzer observations of thirteen z~6 quasars using the Infrared Array Camera (IRAC) and Multiband Imaging Photometer for Spitzer (MIPS). All the quasars except SDSS J000552.34-000655.8 (SDSS J0005-0006) were detected with high S/N in the four IRAC channels and the MIPS 24um band, while SDSS J0005-0006 was marginally detected in the IRAC 8.0um band, and not detected in the MIPS 24um band. We find that most of these quasars have prominent emission from hot dust as evidenced by the observed 24um fluxes. Their spectral energy distributions (SEDs) are similar to those of low-redshift quasars at rest-frame 0.15-3.5 um, suggesting that accretion disks and hot-dust structures for these sources already have reached maturity. However, SDSS J0005-0006 has an unusual SED that lies significantly below low-redshift SED templates at rest-frame 1 and 3.5 um, and thus shows a strong near-IR (NIR) deficit and no hot-dust emission. Type I quasars with extremely small NIR-to-optical flux ratios like SDSS J0005-0006 are not found in low-redshift quasar samples, indicating that SDSS J0005-0006 has different dust properties at high redshift. We combine the Spitzer observations with X-ray, UV/optical, mm/submm and radio observations to determine bolometric luminosities for all the quasars. We find that the four quasars with central black-hole mass measurements have Eddington ratios of order unity.

Observations of the radio galaxy MRC 0316-257 at z=3.13 and the surrounding field are presented. Using narrow- and broad-band imaging obtained with the VLT, 77 candidate Ly-alpha emitters with a rest-frame equivalent width of > 15 A were selected in a ~7'x7' field around the radio galaxy. Spectroscopy of 40 candidate emitters resulted in the discovery of 33 emission line galaxies of which 31 are Ly-alpha emitters with redshifts similar to that of the radio galaxy, while the remaining two galaxies turned out to be [OII] emitters. The Ly-alpha profiles have widths (FWHM) in the range of 120-800 km/s, with a median of 260 km/s. Where the signal-to-noise was large enough, the Ly-alpha profiles were found to be asymmetric, with apparent absorption troughs blueward of the profile peaks, indicative of absorption along the line of sight of an HI mass of at least 2x10^2 - 5x10^4 M_sun. The properties of the Ly-alpha galaxies (faint, blue and small) are consistent with young star forming galaxies which are still nearly dust free. The volume density of Ly-alpha emitting galaxies in the field around MRC 0316-257 is a factor of 3.3+0.5-0.4 larger compared with the density of field Ly-alpha emitters at that redshift. The velocity distribution of the spectroscopically confirmed emitters has a FWHM of 1510 km/s, which is substantially smaller than the width of the narrow-band filter (FWHM ~ 3500 km/s). The peak of the velocity distribution is located within 200 km/s of the redshift of the radio galaxy. We conclude that the confirmed Ly-alpha emitters are members of a protocluster of galaxies at z~3.13. The size of the protocluster is larger than 3.3x3.3 Mpc^2. The mass of this structure is estimated to be > 3-6x10^14 M_sun and could be the progenitor of a cluster of galaxies similar to e.g. the Virgo cluster. (Abridged)

PKS 1138-262 is a massive radio galaxy at z = 2.16 surrounded by overdensities of Lya emitters, Ha emitters, EROs and X-ray emitters. Numerous lines of evidence exist that it is located in a forming cluster. We report on Keck spectroscopy of candidate members of this protocluster, including nine of the 18 X-ray sources detected by Pentericci et al. (2002) in this field. Two of these X-ray sources (not counting PKS 1138-262 itself) were previously confirmed to be members of the protocluster; we have discovered that an additional two (both AGN) are members of a filamentary structure, at least 3.5 Mpc in projection, aligned with the radio jet axis, the 150 kpc-sized emission-line halo, and the extended X-ray emission around the radio galaxy. Three of the nine X-ray sources observed are lower redshift AGN, and three are M-dwarf stars.

We present the results of 20 ksec Chandra observations for each of 5 radio galaxies in the redshift range 2.0 < z < 2.6. For 4 of the 5 targets we detect unresolved X-ray components coincident with the radio nuclei. From spectral analysis of one of the cores and comparison to the empirical radio to X-ray luminosity ratio correlation, we find that obscuring material (n(HI)~10^22 cm^-2) may be surrounding the nuclei. We detect X-ray emission coincident with the radio hotspots or lobes in 4 of the 5 targets, which can be explained by Inverse-Compton (IC) scattering of CMB photons. The magnetic field strengths of ~100-200 muG that we derive agree with the equipartition magnetic field strengths. The relative ease with which the lobe X-ray emission is detected is a consequence of the (1+z)^4 increase in the energy density of the CMB. An HST image of one of the sources shows that the X-ray emission could also be produced by a reservoir of hot, shocked gas, as evidenced by a bright, optical bow-shock. By stacking our data we created a deep, 100 ksec exposure to search for diffuse X-ray emission from intra-cluster gas. We detect no diffuse emission and derive upper limits of ~1e+44 erg/s, thereby ruling out a virialized structure of cluster-size scale at z~2. The average number of soft X-ray sources in our fields is consistent with the number density of AGN in the Chandra Deep Fields. Their angular distribution shows no evidence for large-scale structure associated with the radio galaxies.

In this paper we present deep VLT multi wavelength imaging observations of the field around the radio galaxy PKS 1138-262 aimed at detecting and studying a (potential) proto-cluster centered at this radio source. PKS 1138-262 is a massive galaxy at z=2.16, located in a dense environment as indicated by optical, X-ray and radio observations. We had already found an over-density of Ly-alpha emitting galaxies in this field, consistent with a proto-cluster structure associated with the radio galaxy. In addition, we find 40 candidate H-alpha emitters that have nominal rest frame equivalent width > 25 Angstrom within 1.8 Mpc and 2000 km/s of the radio galaxy. Furthermore, we find 44 objects with I - K > 4.3. This number of extremely red objects (EROs) is about twice the number found in blank field ERO surveys, suggesting that some EROs in this field are part of the proto-cluster. The density of H-alpha emitters and extremely red objects increases towards the radio galaxy, indicating a physical association. From comparisons with other K band, ERO, H-alpha and Ly-alpha surveys, we conclude that PKS 1138-262 is located in a density peak which will evolve into a cluster of galaxies.

Radio galaxy PKS 1138-262 is a massive galaxy at z=2.16, located in a dense environment. We have found an overdensity of Ly-alpha emitting galaxies in this field, consistent with a proto-cluster structure associated with the radio galaxy. Recently, we have discovered forty candidate H-alpha emitters by their excess near infrared narrow band flux. Here, we present infrared spectroscopy of nine of the brightest candidate H-alpha emitters. All these candidates show an emission line at the expected wavelength. The identification of three of these lines with H-alpha is confirmed by accompanying [NII] emission. The spectra of the other candidates are consistent with H-alpha emission at z~2.15, one being a QSO as indicated by the broadness of its emission line. The velocity dispersion of the emitters (360 km/s) is significantly smaller than that of the narrow band filter used for their selection (1600 km/s). We therefore conclude that the emitters are associated with the radio galaxy. The star formation rates (SFRs) deduced from the H-alpha flux are in the range 6-44 M_solar/yr and the SFR density observed is 5-10 times higher than in the HDF-N at z=2.23. The properties of the narrow emission lines indicate that the emitters are powered by star formation and contain very young (< 100 Myr) stellar populations with moderately high metallicities.

We present the results of narrow-band and broad-band imaging with the Very Large Telescope of the field surrounding the radio galaxy TN J0924-2201 at z = 5.2. Fourteen candidate Ly alpha emitters with a rest-frame equivalent width of > 20 A were detected. Spectroscopy of 8 of these objects showed that 6 have redshifts similar to that of the radio galaxy. The density of emitters at the redshift of the radio galaxy is estimated to be a factor 1.5-6.2 higher than in the field, and comparable to the density of Ly alpha emitters in radio galaxy protoclusters at z = 4.1, 3.1 and 2.2. The Ly alpha emitters near TN J0924-2201 could therefore be part of a structure that will evolve into a massive cluster. These observations confirm that substantial clustering of Ly alpha emitters occurs at z > 5 and support the idea that radio galaxies pinpoint high density regions in the early Universe.

The clustering properties of clusters, galaxies and AGN as a function of redshift are briefly discussed. It appears that extremely red objects at z ~ 1, and objects with J-K > 1.7 and photometric redshifts 2 < z_phot < 4 are highly clustered, indicating that a majority of these objects constitutes the progenitors of nearby ellipticals. Similarly clustered seem luminous radio galaxies at z\sim 1, indicating that these objects comprise a short lived phase in the lifetime of these red objects. The high level of clustering furthermore suggests that distant powerful radio galaxies (e.g. z>2) might be residing in the progenitors of nearby clusters -- proto-clusters. A number of observational projects targetting fields with distant radio galaxies, including studies of Lya and Ha emitters, Lyman break galaxies and (sub)mm and X-ray emitters, all confirm that such radio galaxies are located in such proto-clusters. Estimates of the total mass of the proto-clusters are similar to the masses of local clusters. If the total star formation rate which we estimate for the entire proto-clusters is sustained up to z~1, the metals in the hot cluster gas of local clusters can easily be accounted for.

We have carried out narrow band imaging targeted at the Lya line of emitters associated with radio galaxies at 2 < z < 5. Subsequent spectroscopy has confirmed the identity of > 120 Lya emitters and led to the discovery of six proto-clusters up to z = 4.1 with velocity dispersions between 300 and 1000 km/s. The number of emission line galaxies in the observed fields is four to fifteen times higher than in blank fields. These results strongly support the idea that high redshift radio galaxies are the progenitors of central brightest cluster galaxies located in the progenitors of clusters of galaxies.

We present the results of the kinematic study of the extended gas in a sample of 10 high redshift radio galaxies (z~2.5) based on high signal to noise Keck II and VLT long slit spectroscopy. In addition to the typical high surface brightness kinematically perturbed regions (FWHM and velocity shifts >1000 km/s), we find in all objects giant low surface brightness halos which show quieter kinematics with typical emission line FWHM and velocity shifts of \~several hundred km/s. The giant halos often extend for more than 100 kpc and sometimes beyond the radio structures. Their emission line specta are typical of high redshift active galaxies. Continuum is also often detected. The halos are enriched with heavy elements at tens of kpc from the active nucleus. The quasar continuum is the dominant source of ionization of the quiescent halos along the radio axis. The detection of giant quiescent halos in all objects suggests that they could be a common ingredient of high redshift radio galaxies. The radio galaxies seem to be embedded within the halos. The nature, origin of the halos and the cosmological implications are also discussed.

We present J-H-K' photometry for a sample of 45 high redshift quasars found by the Sloan Digital Sky Survey. The sample was originally selected on the basis of optical colors and spans a redshift range from 3.6 to 5.03. Our photometry reflects the rest-frame SED longward of Ly alpha for all redshifts. The results show that the near-IR colors of high redshift quasars are quite uniform. We have modelled the continuum shape of the quasars (from just beyond Ly alpha to ~4000 A) with a power law of the form f_nu \propto nu^alpha, and find <alpha > =-0.57 with a scatter of 0.33. This value is similar to what is found for lower redshift quasars over the same restframe wavelength range, and we conclude that there is hardly any evolution in the continuum properties of optically selected quasars up to redshift 5. The spectral indices found by combining near-IR with optical photometry are in general consistent but slightly flatter than what is found for the same quasars using the optical spectra and photometry alone, showing that the continuum region used to determine the spectral indices can somewhat influence the results.

We present the discovery of three new quasars at z>6 in 1300 deg^2 of SDSS imaging data, J114816.64+525150.3 (z=6.43), J104845.05+463718.3 (z=6.23) and J163033.90+401209.6 (z=6.05). The first two objects have weak Ly alpha emission lines; their redshifts are determined from the positions of the Lyman break. They are only accurate to 0.05 and could be affected by the presence of broad absorption line systems. The last object has a Ly alpha strength more typical of lower redshift quasars. Based on a sample of six quasars at z>5.7 that cover 2870 deg^2 presented in this paper and in Paper I, we estimate the comoving density of luminous quasars at z 6 and M_{1450} < -26.8 to be (8 +/- 3)x10^{-10} Mpc^{-3} (for H_0 = 50 km/s/Mpc, Omega = 1). HST imaging of two z>5.7 quasars and high-resolution ground-based images (seeing 0.4'') of three additional z>5.7 quasars show that none of them is gravitationally lensed. The luminosity distribution of the high-redshfit quasar sample suggests the bright end slope of the quasar luminosity function at z 6 is shallower than Psi L^{-3.5} (2-sigma), consistent with the absence of strongly lensed objects.

We present a method for photometric selection of metal-poor halo giants from the imaging data of the Sloan Digital Sky Survey (SDSS). These stars are offset from the stellar locus in the (g-r) vs. (u-g) color-color diagram. Based on a sample of 29 candidates for which spectra were taken, we derive a selection efficiency of the order of 50%, for stars brighter than $r \sim 17^m$. The candidates selected in 400 deg$^2$ of sky from the SDSS Early Data Release trace the known halo structures (tidal streams from the Sagittarius dwarf galaxy, the Draco dwarf spheroidal galaxy), indicating that such a color-selected sample can be used to study the halo structure even without spectroscopic information. This method, and supplemental techniques for selecting halo stars, such as RR Lyrae stars and other blue horizontal branch stars, can produce an unprecedented three-dimensional map of the Galactic halo based on the SDSS imaging survey.

We present results from a Chandra X-ray Observatory study of the field X-ray source population in the vicinity of the radio galaxy MRC 1138-262. Many serendipitous X-ray sources are detected in an area of 8'x8' around the radio source and 90% are identified in our deep VLT images. The space density of such sources is higher than expected on the basis of the statistics of ROSAT and Chandra deep surveys. The most likely explanation is in terms of a concentration of AGN associated with the protocluster at z=2.16 which was found around the radio galaxy in previous studies. Two sources have a confirmed spectroscopic redshift close to that of the radio galaxy, and for three more sources other observations suggest that they are associated with the protocluster. Four of these five X-ray sources form, together with the radio galaxy, a filament in the plane of the sky. The direction of the filament is similar to that of the radio source axis, the large scale distribution of the other protocluster members, the 150 kpc-sized emission-line halo and the extended X-ray emission associated with the radio galaxy. The majority of optically identified X-ray sources in this field have properties consistent with type I AGN, a few could be soft, low luminosity galaxies, one is probably an obscured (type II) AGN and one is a star. These statistics are consistent with the results of deep X-ray surveys.

We present results on the spectroscopic study of the ionized gas in the high redshift radio galaxy USS0828+193 at z=2.57. Thanks to the high S/N of the emission lines in the Keck spectrum, we have been able to perform a detailed kinematic study by means of the spectral decomposition of the emission line profiles. This study reveals the existence of two types of material in this object: a) a low surface brightness component with apparent quiescent kinematics consistent with gravitational motions and b) a perturbed component with rather extreme kinematics. The quiescent halo extends across the entire object for ~80 kpc. It is enriched with heavy elements and apparently ionized by the continuum from the active nucleus. The properties of the quiescent halo and its origin are discussed in this paper. We propose that it could be part of a structure that surrounds the entire object, although its nature is not clear (a rotating disc? low surface brightness satellites? a cooling flow nebula? material ejected in galactic winds? other?).

Imaging and spectroscopy with the Very Large Telescope have revealed 20 Lyman-alpha emitters within a projected distance of 1.3 Mpc and 600 km/s of the luminous radio galaxy TN J1338-1942 at z=4.1. Compared to the field density of Lyman-alpha emitters, this implies an overdensity on the order of 15. The structure has a projected size of at least 2.7 Mpc x 1.8 Mpc and a velocity dispersion of 325 km/s, which makes it the most distant structure known. Using the galaxy overdensity and assuming a bias parameter b = 3 - 5, the mass is estimated to be ~10^15 M_sun. The radio galaxy itself is surrounded by an uniquely asymmetric Lyman-alpha halo. Taken together with our previous data on PKS 1138-262 at z ~ 2.16, these results suggest that luminous radio sources are excellent tracers of high density regions in the early Universe, which evolve into present-day clusters. The statistics of bright radio sources and of concentrations in the Lyman break galaxy population are consistent with the picture that each of those concentrations harbours an active or passive luminous radio source.

We present new VLT spectroscopic observations of the most distant quasar known, SDSS J1030+0524 at z=6.28 which was recently discovered by the Sloan Digital Sky Survey. We confirm the presence of a complete Gunn-Peterson trough caused by neutral hydrogen in the intergalactic medium. There is no detectable flux over the wavelength range from 8450 to 8710 A. We set an improved limit on the drop of the flux level blueward of the Ly alpha line: a factor of > 200. Below 8450 A the spectrum shows a rise in flux, with a large fraction (> 60 %) of the total emission produced by few narrow features of transmitted flux. We discuss the obvious proximity effect around this quasar, with the presence of transmitted flux with many absorption features in a region of about 23h^{-1} comoving Mpc. If assuming the surrounding medium were completely neutral, the size of this region would imply a quasar lifetime of \~1.3x10^7 years. We also present near-IR spectroscopy of both SDSS J1030+0524 and of SDSS J1306+05, the second most distant quasar known at redshift 6.0. We combine measurements of the CIV line and limits on the HeII emission with the NV line measurements from the optical spectra to derive line ratios, and by implication the abundances of these early quasar environments. The results are indistinguishable from those of lower redshift quasars and indicate little or no evolution in the abundances from z ~ 6 to z ~ 2. The line ratios suggest supersolar metallicities, implying that the first stars around the quasars must have formed at least a few hundreds of Myrs prior to the observation, i.e. at redshift higher than 8.

We study the process of cosmic reionization and estimate the ionizing background in the IGM using the Lyman series absorption in the spectra of the four quasars at 5.7<z<6.3 discovered by the SDSS. We derive the evolution of the ionizing background at high redshifts, using both semi-analytic techniques and cosmological simulations to model the density fluctuations in the IGM. The existence of the complete Ly alpha Gunn-Peterson trough in the spectrum of the z=6.28 quasar SDSS 1030+0524 indicates a photoionization rate Gamma_{-12} at z~6 lower than 0.08, at least a factor of 6 smaller than the value at z~3. The Ly beta and Ly gamma Gunn-Peterson troughs give an even stronger limit Gamma_{-12}<0.02 due to their smaller oscillator strengths, indicating that the ionizing background in the IGM at z~6 is more than 20 times lower than that at z~3. Meanwhile, the volume-averaged neutral hydrogen fraction increases from 10^{-5} at z~3 to >10^{-3} at z~6. At this redshift, the mass-averaged neutral hydrogen fraction is larger than 1%; the mildly overdense regions (delta > 3) are still mostly neutral and the comoving mean free path of ionizing photons is shorter than 8 Mpc. Comparison with simulations of cosmological reionization shows that the observed properties of the IGM at z~6 are typical of those in the era at the end of the overlap stage of reionization when the individual HII regions merge. Thus, z~6 marks the end of the reionization epoch. The redshift of reionization constrains the small scale power of the mass density fluctuations and the star forming efficiency of the first generation of objects.

We present high resolution X-ray observations of the narrow line radio galaxy PKS 1138-262 at z = 2.156 with the ACIS-S detector on the Chandra observatory. These observations show that the X-ray emission from 1138-262 is dominated by emission from the AGN with a (rest frame) 2 to 10 keV luminosity of 4x10^{45} erg/s. The relative X-ray and radio properties of the AGN in 1138-262 are similar to those seen for the AGN in the archetype powerful radio galaxy Cygnus A. Between 10% and 25% (depending on energy) of the X-ray emission from 1138--262 is spatially extended on scales of 10'' to 20''. The extended X-ray emission is elongated, with a major axis aligned with that of the radio source. While the X-ray and radio emissions are elongated on similar scales and position angles, there is no one-to-one correspondence between the radio and X-ray features in the source. The most likely origin for the extended X-ray emission in 1138-262 is thermal emission from shocked gas, although we cannot rule-out a contribution from inverse Compton emission. If the emission is from hot gas, the gas density is 0.05 cm^{-3} and the gas mass is 2.5x10^{12} M_solar. The pressure in this hot gas is adequate to confine the radio emitting plasma and the optical line emitting gas. We set an upper limit of 1.5x10^{44} erg/s to the 2 to 10 keV luminosity of any normal cluster atmosphere associated with 1138-262. No emission was detected from any of the Ly alpha emitting galaxies in the (proto-) cluster around 1138-262, outside of the Ly alpha halo of 1138-262 itself, to a 2 to 10keV luminosity limit of 1.2x10^{43}erg/s.

We have recently discovered a forming cluster around the radio galaxy MRC 1138-262 at redshift 2.2. Besides the population of Ly alpha emitting galaxies that have been confirmed spectroscopically, we have detected many candidate H alpha emitters that seem to have a different spatial distribution from the other galaxies: they are more clustered towards the center of the cluster and seem to be distributed along the same direction as the radio source. We present here the characteristics of the Ly alpha and H alpha emitters and study the nature of these populations.

High redshift radio galaxies are amongst the most massive galaxies in the early Universe and have properties expected from central galaxies in forming clusters. We are carrying out an observational programme on the VLT to find and study galaxy proto clusters around radio galaxies at redshifts 2 < z < 4. First, we use narrow band imaging to select candidate galaxies which show excess Lyman alpha emission at redshifts similar to the central radio galaxy. Then, we use multi object spectroscopy to confirm the redshifts of these candidates and measure the velocity dispersion of the cluster members. Our goal is to observe a sample of about 10 targets and investigate galaxy overdensities as a function of redshift. Here, we report on the current progress of the programme and show some preliminary results which include the discovery of a structure of galaxies at redshift 4.1.

We present the results from a survey of i-dropout objects selected from ~1550 deg^2 of multicolor imaging data from the Sloan Digital Sky Survey, to search for luminous quasars at z>5.8. Objects with i*-z*>2.2 and z*<20.2 are selected, and follow-up J band photometry is used to separate L and T type cool dwarfs from high-redshift quasars. We describe the discovery of three new quasars, at z=5.82, 5.99 and 6.28, respectively. Their spectra show strong and broad Ly alpha+NV emission lines, and very strong Ly alpha absorption, with a mean continuum decrement D_A > 0.90. The ARC 3.5m spectrum of the z=6.28 quasar shows that over a range of 300 A immediately blueward of the Ly alpha emission, the average transmitted flux is only 0.003 +/-0.020 times that of the continuum level, consistent with zero flux, and suggesting a tentative detection of the complete Gunn-Peterson trough. The existence of strong metal lines suggests early chemical enrichment in the quasar enviornment. The three new objects, together with the previously published z=5.8 quasar form a complete color-selected flux-limited sample at z>5.8. We estimate that at $z=6$, the comoving density of luminous quasars at M_1450 < -26.89 (h=0.5, Omega=1)is 1.1x10^-9 Mpc^-3. This is a factor of ~2 lower than that at z~5, and is consistent with an extrapolation of the observed quasar evolution at low-z. We discuss the contribution of quasars to the ionizing background at z~6. The luminous quasars discussed in the paper have central black hole masses of several times 10^9 M_sun by the Eddington argument. Their observed space density provides a sensitive test of models of quasar and galaxy formation at high redshift. (Abridged)

The radio galaxy MRC 1138-262 at z=2.16 is most likely a brightest cluster galaxy in an early stage of evolution. Here we present observations of the luminous emission line halo and the stellar components of this radio galaxy. Optical narrow band imaging shows a very extended (about 160 kpc) and luminous Lyman alpha halo. Infrared narrow band imaging reveals a much smaller H alpha halo with a very different morphology. We advocate a model in which the inner part of the halo is photoionized by direct AGN illumination or by UV photons from young stars. Far from the nucleus (about 25 kpc), there is a region of greatly enhanced Lyman alpha emission. At this location, it is likely that shock ionization is important as indicated by a bend in the radio jet. Spectroscopy of several continuum clumps in the halo shows that, although there are striking differences between the emission and absorption features of the spectra of various regions, they have properties similar to those of Lyman-break galaxies. This is further evidence for a scenario in which massive galaxies form hierarchically from smaller building blocks.

We present the results of a near infrared imaging program of a sample of 19 radio galaxies with redshift between 1.7 and 3.2, using the NICMOS Camera 2 on the Hubble Space Telescope. The galaxies were observed in H-band which, for 18 of the 19 targets, samples the rest-frame optical emission longwards of the 4000 A break. The images show a wide range of morphologies, from unresolved to compact systems, to systems with substructures such as multiple emission peaks, to systems comprised of several components spread over large areas. The morphologies of 5 of the lowest redshift targets are well represented by a de Vaucouleurs profile, consistent with them being elliptical galaxies. Their average effective radius is a factor of 2 smaller than that of z=1 3CR radio galaxies at similar restframe wavelengths. The near infrared continuum light is generally well aligned with the radio axis and the aligned light has typical V-H colors of around 3.5-4. For several galaxies where WFPC2 images were available, we computed a high resolution map of the optical-to-infrared spectral indices: all multi-component systems present net color gradients between the various clumps. We argue that in these galaxies the continuum emission has probably a stellar origin. Finally in many cases we observe nearby faint objects close to the radio sources. The number density of these faint objects is slightly higher than what is found in the deep NICMOS parallel observations of random fields. These objects also tend to be aligned with the radio sources axis, suggesting that they may be related to the presence of the AGN. (shortened astro-ph version)

We summarize observations of the properties of powerful radio galaxies and their environments from z = 0 to 5. These data show that some high redshift radio galaxies inhabit regions of gas and galaxy over-densities indicative of proto-cluster environments.

VLT spectroscopic observations of Ly alpha--excess objects in the field of the clumpy radio galaxy 1138-262 at z=2.16 have led to the discovery of 14 galaxies and one QSO at approximately the same distance as the radio galaxy. All galaxies have redshifts in the range 2.16 \pm 0.02, centered around the redshift of the radio galaxy, and are within a projected physical distance of 1.5 Mpc from it. The velocity distribution suggests that there are two galaxy subgroups having velocity dispersions of ~500 km/s and ~300 km/s and a relative velocity of 1800 km/s. If these are virialized structures, the estimated dynamical masses for the subgroups are ~9 and ~3x10E13 solar masses respectively, implying a total mass for the structure of more than 10E14 solar masses. The new observations, together with previous results, suggest that the structure of galaxies around 1138-262 is likely to be a forming cluster.

We present new deep multi-frequency radio-polarimetric images of a sample of high redshift radio galaxies (HzRGs), having redshift between 1.7 and 4.1. The radio data at 4.7 and 8.2 GHz were taken with the Very Large Array in the A configuration and provide a highest angular resolution of 0.2''. Maps of total intensity, radio spectral index, radio polarization and internal magnetic field are presented for each source. The morphology of most objects is that of standard FRII double radio sources, but several contain multiple hot-spots in one or both lobes. Compared to similar samples of HzRGs previously imaged, there is a higher fraction (29%) of compact steep spectrum sources (i.e. sources with a projected linear size less than 20 kpc). Radio cores are identified in about half of the sample and tend to have relatively steep spectra (alpha < -1). Polarization is detected in all but 4 sources, with typical polarization at 8.2 GHz of around 10-20%. The Faraday rotation can be measured in most of the radio galaxies: the observed rotation measure (RM) of 8 radio sources exceeds 100 rad m^{-2} in at least one of the lobes, with large gradients between the two lobes. We find no dependence of Faraday rotation with other properties of the radio sources. If the origin of the Faraday rotation is local to the sources, as we believe, then the intrinsic RM is more than a 1000 rad m^{-2}. Because low redshift radio galaxies residing at the center of clusters usually show extreme RMs, we suggest that the high-z large RM sources also lie in very dense environments. Finally, we find that the fraction of powerful radio galaxies with extreme Faraday rotation increases with redshift, as would be expected if their average environment tends to become denser with decreasing cosmic epoch.

Radio, optical and X-ray observations of the powerful radio galaxy PKS 1138-262 at z=2.156 have suggested that this galaxy is a massive galaxy in the center of a forming cluster. We have imaged 1138-262 and the surrounding 38 square arcminute field with the Very Large Telescope in a broad band and a narrow band encompassing the redshifted Lya emission. We detect 50 objects with rest equivalent width larger than 20 A and a luminous, highly extended Lya halo around 1138-262. If the radio galaxy is at the center of a forming cluster, these objects are candidate Lya emitting cluster galaxies.

Observations of the powerful radio galaxy 1138-262 at z=2.2 suggest that this galaxy is a massive galaxy in the center of a forming cluster. We have imaged the field of 1138-262 with the Very Large Telescope in a narrow band encompassing the redshifted Lyman alpha emission. We detect 34 probable line emitters in the 6'x6' field. These are candidate Lyman alpha emitting galaxies in the supposed cluster. The observations also reveal that the radio galaxy is enclosed by a very extended Lyman alpha halo (160 kpc).

The results of a deep submillimetre observation using SCUBA of the powerful radio galaxy 3C324, at redshift z=1.206, are presented. At 850 microns, emission from the location of the host radio galaxy is marginally detected at the 4.2 sigma level, 3.01 +/- 0.72 mJy, but there is no detection of emission at 450 microns to a 3 sigma limit of 21 mJy. A new 32 GHz radio observation using the Effelsberg 100m telescope confirms that the sub-millimetre signal is not associated with synchrotron emission. These observations indicate that both the mass of warm dust within 3C324, and the star formation rate, lie up to an order of magnitude below the values recently determined for radio galaxies at z = 3 to 4. The results are compared with dust masses and star formation rates derived in other ways for 3C324.

We present Hubble Space Telescope images of 11 high redshift radio galaxies (between $z=2.3$ and $z=3.6$). The galaxies were observed with the WFPC2 camera in a broad band filter (F606W or F707W, roughly equivalent to V or R-band), for 2 orbits each. We find that on the scale of the HST observations there is a wide variety of morphological structures of the hosting galaxies: most objects have a clumpy, irregular appearance, consisting of a bright nucleus and a number of smaller components, suggestive of merging systems. Some observed structures could be due (at least partly) to the presence of dust distributed through the galaxies. The UV continuum emission is generally elongated and aligned with the axis of the radio sources, however the characteristics of the ``alignment effect'' differ from case to case, suggesting that the phenomenon cannot be explained by a single physical mechanism. We compare the properties of our radio galaxies with those of the UV dropout galaxies and conclude that (i) the most massive radio galaxies may well evolve from an aggregate of UV dropout galaxies and (ii) high redshift radio galaxies probably evolve into present day brightest cluster galaxies.

A third possible explanation for the X-ray emission from the z=2.156 radio galaxy 1138-262 is thermal emission from a very dense, sub-cluster `halo', perhaps associated with the (forming) cD galaxy on a scale < 100 kpc. For instance, increasing the gas density from 0.01 cm^{-3} to 0.1 cm^{-3} would decrease the required hot gas mass by the same factor, and could possibly alleviate constraints on cosmological structure formation models. The pressure in this gas would be very high (10^{-9} dynes cm^{-2}), comparable to the pressure in the optical line emitting nebulosity and to the minimum pressure in the radio source, and the cooling time would be short (< fewx10^{8} years). Circumstantial evidence for such very dense, hot gas enveloping some high z radio sources has been reviewed by Fabian et al. (1986). Fabian (1991) suggests that in some cases the `inferred pressures are close to the maximum that can be obtained by gas cooling in a potential well of a galaxy' (ie. the cooling time = gravitational free-fall time), and he designates such systems as `maximal cooling flows', with implied cooling flow rates up to 2000 M_solar year^{-1}. High resolution X-ray imaging with AXAF should be able to test whether 1138-262 has a normal cluster atmosphere, a `maximal cooling flow', or an unusually X-ray loud AGN.

We present X-ray observations of the narrow line radio galaxy 1138-262 at z = 2.156 with the High Resolution Imager (HRI) on ROSAT. Observations at other wave-bands, and in particular extremely high values of Faraday rotation of the polarized radio emission, suggest that the 1138-262 radio source is in a dense environment, perhaps a hot, cluster-type atmosphere. We detect X-ray emission from the vicinity of 1138-262, and we discuss possible origins for this emission. The X-ray, optical, and radio data all favor thermal emission from a hot cluster atmosphere as the mechanism responsible for the X-rays, although we cannot rule out a contribution from the active nucleus. If this interpretation is correct, then 1138-262 becomes the most distant, by far, of known X-ray emitting clusters. The X-ray luminosity for 1138-262 is 6.7+/-1.3 x 10**44 ergs/sec for emitted energies between 2 keV and 10 keV.