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)

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 discovered Bowen emission arising from a strongly lensed ($\mu>20$) transient stellar object hosted in the Sunburst arc at z=2.37. Unique ultraviolet lines emerge at the location of the transient. In particular, narrow ($\sigma_v$ ~ 40 km/s) ionisation lines of several chemical species (Fe, C, Si) fluoresce after being exposed to Lyman series (e.g., Ly$\alpha$, Ly$\beta$) of H, HeI and HeII that pump selectively their atomic levels. Similarities with the Resonance-Enhanced Two Photon-Ionisation (RETPI) spectral features observed in the circum-stellar Weigelt blobs of Eta~Carinae are observed. Data from VLT/MUSE, X-Shooter and ESPRESSO observations (the latter placed at the focus of the four UTs) at increasing spectral resolution of R=2500, 11400 and R=70000, respectively, confirm the lasing action since at least 3.3 years (~ 1 year rest-frame), and probe the circum-stellar dense gas condensations in radiation-rich conditions. We discuss the physical origin of the transient event, which remains unclear. We expect such transient events (including also supernova or impostors) will be easily recognised with ELTs thanks to high angular resolution provided by adaptive optics and large collecting area, especially in modest ($\mu < 3$) magnification regime.

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

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

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

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

State-of-the-art models of massive black hole formation postulate that quasars at $z>6$ reside in extreme peaks of the cosmic density structure in the early universe. Even so, direct observational evidence of these overdensities is elusive, especially on large scales ($\gg$1 Mpc) as the spectroscopic follow-up of $z>6$ galaxies is observationally expensive. Here we present Keck / DEIMOS optical and IRAM / NOEMA millimeter spectroscopy of a $z\sim6$ Lyman-break galaxy candidate originally discovered via broadband selection, at a projected separation of 4.65 physical Mpc (13.94 arcmin) from the luminous $z$=6.308 quasar J1030+0524. This well-studied field presents the strongest indication to date of a large-scale overdensity around a $z>6$ quasar. The Keck observations suggest a $z\sim6.3$ dropout identification of the galaxy. The NOEMA 1.2mm spectrum shows a 3.5$\sigma$ line that, if interpreted as [CII], would place the galaxy at $z$=6.318 (i.e., at a line-of-sight separation of 3.9 comoving Mpc assuming that relative proper motion is negligible). The measured [CII] luminosity is $3\times10^8$ L$_\odot$, in line with expectations for a galaxy with a star formation rate $\sim15$ M$_\odot$ yr$^{-1}$, as inferred from the rest-frame UV photometry. Our combined observations place the galaxy at the same redshift as the quasar, thus strengthening the overdensity scenario for this $z>6$ quasar. This pilot experiment demonstrates the power of millimeter-wavelength observations in the characterization of the environment of early quasars.

We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $\sigma$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}\sigma$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $\sigma_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $\sigma_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $\sigma_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations.

We report the discovery of a galaxy overdensity around a Compton-thick Fanaroff-Riley type II (FRII) radio galaxy at z=1.7 in the deep multiband survey around the z=6.3 QSO SDSS J1030+0524. Based on a 6hr VLT/MUSE and on a 4hr LBT/LUCI observation, we identify at least eight galaxy members in this structure with spectroscopic redshift z=1.687-1.699, including the FRII galaxy at z=1.699. Most of the identified overdensity members are blue, compact galaxies that are actively forming stars at rates of $\sim$8-60 $M_{\odot}$ yr$^{-1}$. Based on a 500ks Chandra ACIS-I observation we found that the FRII nucleus hosts a luminous QSO ($L_{2-10keV}=1.3\times10^{44}$ erg s$^{-1}$, intrinsic and rest-frame) that is obscured by Compton-thick absorption ($N_H=1.5\pm0.6\times 10^{24}$ cm$^{-2}$). Our Chandra observation, the deepest so far for a distant FRII within a galaxy overdensity, revealed significant diffuse X-ray emission within the region covered by the overdensity. In particular, X-ray emission extending for $\sim$240 kpc is found around the Eastern lobe of the FRII. Four out of the six MUSE star forming galaxies in the overdensity are distributed in an arc-like shape at the edge of this diffuse X-ray emission. The probability of observing by chance four out of the six $z=1.7$ sources at the edge of the diffuse emission is negligible. In addition, these four galaxies have the highest specific star formation rates of the MUSE galaxies in the overdensity and lie above the main sequence of field galaxies of equal stellar mass at z=1.7. We propose that the diffuse X-rays originate from an expanding bubble of gas that is shock-heated by the FRII jet, and that star formation is promoted by the compression of the cold interstellar medium of the galaxies around the bubble, which may be remarkable evidence of positive AGN feedback on cosmological scales. [shortened version]

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

By means of grid-based, 3D hydrodynamical simulations we study the formation of second generation (SG) stars in a young globular cluster (GC) of mass 10^7 Msun, the possible progenitor of an old GC with a present mass ~(1-5) * 10^6 Msun. The cluster accretes external gas as its first generation (FG) asymptotic giant branch (AGB) stars release their ejecta and SG stars form. We consider two models characterised by different densities of the external gas. In both cases, we find that a very compact SG subsystem with central density >10^5M sun/pc^3 forms in the innermost regions of the cluster. The low-density model forms a population of extreme SG stars with high helium enhancement, followed by the formation of another SG group out of a mix of pristine gas and AGB ejecta and characterised by a modest helium enhancement. On the other hand, the high-density model forms in prevalence SG stars with modest helium enhancement. Our simulations illustrate the dynamical processes governing the formation of SG populations in GCs and shed light on the structural properties emerging at the end of this phase. The newly born SG groups have different concentrations, with more extreme SG stars more centrally concentrated than those with less extreme chemical abundances. The very high density of the SG subsystems implies that SG massive stars, if formed, might suffer frequent close encounters, collisions and gas stripping, thus possibly contributing further gas to the SG formation.

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

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 report on a superdense star-forming region with an effective radius (R_e) smaller than 13 pc identified at z=6.143 and showing a star-formation rate density \Sigma_SFR~1000 Msun/yr/kpc2 (or conservatively >300 Msun/yr/kpc2). Such a dense region is detected with S/N>40 hosted by a dwarf extending over 440 pc, dubbed D1 (Vanzella et al. 2017b). D1 is magnified by a factor 17.4+/-5.0 behind the Hubble Frontier Field galaxy cluster MACS~J0416 and elongated tangentially by a factor 13.2+/-4.0 (including the systematic errors). The lens model accurately reproduces the positions of the confirmed multiple images with a r.m.s. of 0.35", and the tangential stretch is well depicted by a giant multiply-imaged Lya arc. D1 is part of an interacting star-forming complex extending over 800 pc. The SED-fitting, the very blue ultraviolet slope (\beta ~ -2.5, F(\lambda) ~ \lambda^\beta) and the prominent Lya emission of the stellar complex imply that very young (< 10-100 Myr), moderately dust-attenuated (E(B-V)<0.15) stellar populations are present and organised in dense subcomponents. We argue that D1 (with a stellar mass of 2 x 10^7 Msun) might contain a young massive star cluster of M < 10^6 Msun and Muv~-15.6 (or m_uv=31.1), confined within a region of 13 pc, and not dissimilar from some local super star clusters (SSCs). The ultraviolet appearance of D1 is also consistent with a simulated local dwarf hosting a SSC placed at z=6 and lensed back to the observer. This compact system fits into some popular globular cluster formation scenarios. We show that future high spatial resolution imaging (e.g., E-ELT/MAORY-MICADO and VLT/MAVIS) will allow us to spatially resolve light profiles of 2-8 pc.

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

ESPRESSO is an extremely stable high-resolution spectrograph which is currently being developed for the ESO VLT. With its groundbreaking characteristics it is aimed to be a "science machine", i.e. a fully-integrated instrument to directly extract science information from the observations. In particular, ESPRESSO will be the first ESO instrument to be equipped with a dedicated tool for the analysis of data, the Data Analysis Software (DAS), consisting in a number of recipes to analyze both stellar and quasar spectra. Through the new ESO Reflex GUI, the DAS (which will implement new algorithms to analyze quasar spectra) is aimed to get over the shortcomings of the existing software providing multiple iteration modes and full interactivity with the data.

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.

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.

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

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

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 the results from a $\sim500$ ks Chandra observation of the $z=6.31$ QSO SDSS J1030+0524. This is the deepest X-ray observation to date of a $z\sim6$ QSO. The QSO is detected with a total of 125 net counts in the full ($0.5-7$ keV) band and its spectrum can be modeled by a single power-law model with photon index of $\Gamma = 1.81 \pm 0.18$ and full band flux of $f=3.95\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$. When compared with the data obtained by XMM-Newton in 2003, our Chandra observation in 2017 shows a harder ($\Delta \Gamma \approx -0.6$) spectrum and a 2.5 times fainter flux. Such a variation, in a timespan of $\sim2$ yrs rest-frame, is unexpected for such a luminous QSO powered by a $> 10^9 \: M_{\odot}$ black hole. The observed source hardening and weakening could be related to an intrinsic variation in the accretion rate. However, the limited photon statistics does not allow us to discriminate between an intrinsic luminosity and spectral change, and an absorption event produced by an intervening gas cloud along the line of sight. We also report the discovery of diffuse X-ray emission that extends for 30"x20" southward the QSO with a signal-to-noise ratio of $\sim$6, hardness ratio of $HR=0.03_{-0.25}^{+0.20}$, and soft band flux of $f_{0.5-2 \: keV}= 1.1_{-0.3}^{+0.3} \times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$, that is not associated to a group or cluster of galaxies. We discuss two possible explanations for the extended emission, which may be either associated with the radio lobe of a nearby, foreground radio galaxy (at $z \approx 1-2$), or ascribed to the feedback from the QSO itself acting on its surrounding environment, as proposed by simulations of early black hole formation.

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.

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

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.

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

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]

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

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.

Many cosmological studies predict that early supermassive black holes (SMBHs) can only form in the most massive dark matter halos embedded within large scale structures marked by galaxy over-densities that may extend up to 10 physical Mpc. This scenario, however, has not been confirmed observationally, as the search for galaxy over-densities around high-z quasars has returned conflicting results. The field around the z=6.28 quasar SDSSJ1030+0524 (J1030) is unique for multi-band coverage and represents an excellent data legacy for studying the environment around a primordial SMBH. In this paper we present wide-area (25x25 arcmin), Y- and J-band imaging of the J1030 field obtained with the near infrared camera WIRCam at the Canada-France-Hawaii Telescope (CFHT). We built source catalogues in the Y- and J-band, and matched those with our photometric catalogue in the r, z, i bands presented in Morselli et al. (2014). We used these new infrared data together with H and K and Spitzer/IRAC data to refine our selection of Lyman Break Galaxies (LBGs), extending our selection criteria to galaxies in the range 25.2<zAB<25.7. We selected 21 robust high-z candidates in the J1030 field with photometric redshift around 6 and colors i-z>=1.3. We found a significant asymmetry in the distribution of the high-z galaxies in J1030, supporting the existence of a coherent large-scale structure around the quasar. We compared our results with those of Bowler et al. (2015), who adopted similar LBGs selection criteria, and estimated an over-density of galaxies in the field of delta = 2.4, which is significant at >4 sigma. The over-density value and its significance are higher than those found in Morselli et al. (2014), and we interpret this as evidence of an improved LBG selection.

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.

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

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.

Star forming galaxies have long been considered the dominant sources of the cosmic ultraviolet background radiation at early epochs. However, observing and characterizing the galaxy population with significant ionizing emission has proven to be challenging. In particular, the fraction of ionizing radiation that escapes the local environment to the intergalactic medium is poorly known. We investigate the relation between the escape fraction and galaxy luminosity. We combine the deep ultraviolet observations of Hubble Ultra Deep Field (UVUDF) with the deep Multi Unit Spectroscopic Explorer (MUSE) observations of the same field, collecting a sample of 165 faint star forming galaxies in the $3 < z < 4$ redshift range with deep rest-frame observations of the Lyman continuum. In our sample, we do not find any galaxy with significant emission of LyC radiation. We bin the galaxies in various redshift and brightness intervals and stack their images. From stacked images we estimate the relative escape fraction upper limits as a function of the luminosity. Thanks to the depth of the sample we measure meaningful 1$\sigma$ upper limits of $f_{esc,rel} < 0.07, 0.2$ and 0.6 at $L \sim L_{\rm z=3}^{*}, 0.5L_{\rm z=3}^{*}$ and $0.1L_{\rm z=3}^{*}$, respectively. We use our estimates and theoretical predictions from the literature to study a possible dependence of the escape fraction on galaxy luminosity by modelling the ionizing background with different prescriptions of $f_{\rm esc} (M_{\rm UV})$. We show that the understanding of the luminosity dependence hinges on the ability to constrain the escape fraction down to $M_{\rm UV} \sim -18$ mag in the future.

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.

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

In this work, we investigate the abundance and distribution of metals in the intergalactic medium (IGM) at $\langle z \rangle \simeq 2.8$ through the analysis of an ultra-high signal-to-noise ratio UVES spectrum of the quasar HE0940-1050. In the CIV forest, our deep spectrum is sensitive at $3\,\sigma$ to lines with column density down to $\log N_{\rm CIV} \simeq 11.4$ and in 60 per cent of the considered redshift range down to $\simeq11.1$. In our sample, all HI lines with $\log N_{\rm HI} \ge 14.8$ show an associated CIV absorption. In the range $14.0 \le \log N_{\rm HI} <14.8$, 43 per cent of HI lines has an associated CIV absorption. At $\log N_{\rm HI} < 14.0$, the detection rates drop to $<10$ per cent, possibly due to our sensitivity limits and not to an actual variation of the gas abundance properties. In the range $\log N_{\rm HI} \ge 14$, we observe a fraction of HI lines with detected CIV a factor of 2 larger than the fraction of HI lines lying in the circum-galactic medium (CGM) of relatively bright Lyman-break galaxies hosted by dark matter haloes with $\langle M\rangle \sim10^{12}$ M$_{\odot}$. The comparison of our results with the output of a grid of photoionization models and of two cosmological simulations implies that the volume filling factor of the IGM gas enriched to a metallicity $\log Z/Z_{\odot} \ge -3$ should be of the order of $\sim 10-13$ percent. In conclusion, our results favour a scenario in which metals are found also outside the CGM of bright star-forming galaxies, possibly due to pollution by lower mass objects and/or to an early enrichment by the first sources.

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

We use a sample of 1669 QSOs ($r<20.15$, $3.6<z<4.0$) from the BOSS survey to study the intrinsic shape of their continuum and the Lyman continuum photon escape fraction (f$_{esc}$), estimated as the ratio between the observed flux and the expected intrinsic flux (corrected for the intergalactic medium absorption) in the wavelength range 865-885 \AA\ rest-frame. Modelling the intrinsic QSO continuum shape with a power-law, $F_{\lambda}\propto\lambda^{-\gamma}$, we find a median $\gamma=1.30$ (with a dispersion of $0.38$, no dependence on the redshift and a mild intrinsic luminosity dependence) and a mean f$_{esc}=0.75$ (independent of the QSO luminosity and/or redshift). The f$_{esc}$ distribution shows a peak around zero and a long tail of higher values, with a resulting dispersion of $0.7$. If we assume for the QSO continuum a double power-law shape (also compatible with the data) with a break located at $\lambda_{\rm br}=1000$ \AA\ and a softening $\Delta\gamma=0.72 $ at wavelengths shorter than $\lambda_{\rm br}$, the mean f$_{esc}$ rises to $=0.82$. Combining our $\gamma$ and f$_{esc}$ estimates with the observed evolution of the AGN luminosity function (LF) we compute the AGN contribution to the UV ionizing background (UVB) as a function of redshift. AGN brighter than one tenth of the characteristic luminosity of the LF are able to produce most of it up $z\sim 3$, if the present sample is representative of their properties. At higher redshifts a contribution of the galaxy population is required. Assuming an escape fraction of Lyman continuum photons from galaxies between $5.5$ and $7.6\%$, independent of the galaxy luminosity and/or redshift, a remarkably good fit to the observational UVB data up to $z\sim 6$ is obtained. At lower redshift the extrapolation of our empirical estimate agrees well with recent UVB observations, dispelling the so-called Photon Underproduction Crisis.

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

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

We present 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 VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

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]

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

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.

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

ESPRESSO is an extremely stable high-resolution spectrograph which is currently being developed for the ESO VLT. With its groundbreaking characteristics it is aimed to be a "science machine", i.e., a fully-integrated instrument to directly extract science information from the observations. In particular, ESPRESSO will be the first ESO instrument to be equipped with a dedicated tool for the analysis of data, the Data Analysis Software (DAS), consisting in a number of recipes to analyze both stellar and quasar spectra. Through the new ESO Reflex GUI, the DAS (which will implement new algorithms to analyze quasar spectra) is aimed to get over the shortcomings of the existing software providing multiple iteration modes and full interactivity with the data.

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

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

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.

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.

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 a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

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.

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.

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

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

We present the 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).

Context. The study of the galaxy stellar mass function (SMF) in relation to the galaxy environment and the stellar mass density profile, rho(r), is a powerful tool to constrain models of galaxy evolution. Aims. We determine the SMF of the z=0.44 cluster of galaxies MACS J1206.2-0847 separately for passive and star-forming (SF) galaxies, in different regions of the cluster, from the center out to approximately 2 virial radii. We also determine rho(r) to compare it to the number density and total mass density profiles. Methods. We use the dataset from the CLASH-VLT survey. Stellar masses are obtained by SED fitting on 5-band photometric data obtained at the Subaru telescope. We identify 1363 cluster members down to a stellar mass of 10^9.5 Msolar. Results. The whole cluster SMF is well fitted by a double Schechter function. The SMFs of cluster SF and passive galaxies are statistically different. The SMF of the SF cluster galaxies does not depend on the environment. The SMF of the passive population has a significantly smaller slope (in absolute value) in the innermost (<0.50 Mpc), highest density cluster region, than in more external, lower density regions. The number ratio of giant/subgiant galaxies is maximum in this innermost region and minimum in the adjacent region, but then gently increases again toward the cluster outskirts. This is also reflected in a decreasing radial trend of the average stellar mass per cluster galaxy. On the other hand, the stellar mass fraction, i.e., the ratio of stellar to total cluster mass, does not show any significant radial trend. Conclusions. Our results appear consistent with a scenario in which SF galaxies evolve into passive galaxies due to density-dependent environmental processes, and eventually get destroyed very near the cluster center to become part of a diffuse intracluster medium.

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

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 investigated the presence of galaxy overdensities around four $z\sim6$ QSOs, namely SDSS J1030+0524 (z = 6.28), SDSS J1148+5251 (z = 6.41), SDSS J1048+4637 (z = 6.20) and SDSS J1411+1217 (z = 5.95), through deep $r$-, $i$- and $z$- band imaging obtained with the wide-field ($\sim23'\times25'$) Large Binocular Camera (LBC) at the Large Binocular Telescope (LBT). We adopted color-color selections within the $i-z$ vs $r-z$ plane to identify samples of $i$-band dropouts at the QSO redshift and measure their relative abundance and spatial distribution in the four LBC fields, each covering $\sim8\times8$ physical Mpc at $z\sim6$. The same selection criteria were then applied to $z$-band selected sources in the $\sim$1 deg$^2$ Subaru-XMM Newton Deep Survey to derive the expected number of dropouts over a blank LBC-sized field ($\sim$0.14 deg$^2$). The four observed QSO fields host a number of candidates larger than what is expected in a blank field. By defining as $i$-band dropouts objects with $z_{AB}<25$, $i-z>1.4$ and undetected in the $r$-band, we found 16, 10, 9, 12 dropouts in SDSS J1030+0524, SDSS J1148+5251, SDSS J1048+4637, and SDSS J1411+1217, respectively, whereas only 4.3 such objects are expected over a 0.14 deg$^2$ blank field. This corresponds to overdensity significances of 3.3, 1.9, 1.7, 2.5$\sigma$, respectively. By considering the total number of dropouts in the four LBC fields and comparing it with what is expected in four blank fields of 0.14 deg$^2$ each, we find that high-z QSOs reside in overdense environments at the $3.7\sigma$ level. This is the first direct and unambiguous measurement of the large scale structures around $z\sim6$ QSOs. [shortened]

By combining cosmological simulations with Frontier Fields project lens models we find that, in the most optimistic case, galaxies as faint as $m \approx 33 - 34$ (AB magnitude at $1.6 \rm \mu m$) can be detected in the Frontier Fields. Such faint galaxies are hosted by dark matter halos of mass $\sim10^9 M_\odot$ and dominate the ionizing photon budget over currently observed bright galaxies, thus allowing for the first time the investigation of the dominant reionization sources. In addition, the observed number of these galaxies can be used to constrain the role of feedback in suppressing star formation in small halos: for example, if galaxy formation is suppressed in halos with circular velocity $v_c < 50$ km s$^{-1}$, galaxies fainter than $m=31$ should not be detected in the FFs.

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

We present a strong lensing system in which a double source is imaged 5 times by 2 early-type galaxies. We take advantage in this target of the multi-band photometry obtained as part of the CLASH program, complemented by the spectroscopic data of the VLT/VIMOS and FORS2 follow-up campaign. We use a photometric redshift of 3.7 for the source and confirm spectroscopically the membership of the 2 lenses to the galaxy cluster MACS J1206.2-0847 at redshift 0.44. We exploit the excellent angular resolution of the HST/ACS images to model the 2 lenses in terms of singular isothermal sphere profiles and derive robust effective velocity dispersions of (97 +/- 3) and (240 +/- 6) km/s. The total mass distribution of the cluster is also well characterized by using only the local information contained in this lensing system, that is located at a projected distance of more than 300 kpc from the cluster luminosity center. According to our best-fitting lensing and composite stellar population models, the source is magnified by a total factor of 50 and has a luminous mass of about (1.0 +/- 0.5) x 10^{9} M_{Sun}. By combining the total and luminous mass estimates of the 2 lenses, we measure luminous over total mass fractions projected within the effective radii of 0.51 +/- 0.21 and 0.80 +/- 0.32. With these lenses we can extend the analysis of the mass properties of lens early-type galaxies by factors that are about 2 and 3 times smaller than previously done with regard to, respectively, velocity dispersion and luminous mass. The comparison of the total and luminous quantities of our lenses with those of astrophysical objects with different physical scales reveals the potential of studies of this kind for investigating the internal structure of galaxies. These studies, made possible thanks to the CLASH survey, will allow us to go beyond the current limits posed by the available lens samples in the field.

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 acronym ESPRESSO stems for Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations; this instrument will be the next VLT high resolution spectrograph. The spectrograph will be installed at the Combined-Coud\'e Laboratory of the VLT and linked to the four 8.2 m Unit Telescopes (UT) through four optical Coud\'e trains. ESPRESSO will combine efficiency and extreme spectroscopic precision. ESPRESSO is foreseen to achieve a gain of two magnitudes with respect to its predecessor HARPS, and to improve the instrumental radial-velocity precision to reach the 10 cm/s level. It can be operated either with a single UT or with up to four UTs, enabling an additional gain in the latter mode. The incoherent combination of four telescopes and the extreme precision requirements called for many innovative design solutions while ensuring the technical heritage of the successful HARPS experience. ESPRESSO will allow to explore new frontiers in most domains of astrophysics that require precision and sensitivity. The main scientific drivers are the search and characterization of rocky exoplanets in the habitable zone of quiet, nearby G to M-dwarfs and the analysis of the variability of fundamental physical constants. The project passed the final design review in May 2013 and entered the manufacturing phase. ESPRESSO will be installed at the Paranal Observatory in 2016 and its operation is planned to start by the end of the same year.

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

We study the observed cosmic ionizing background as a constraint on the nature of the sources responsible for the reionization of the Universe. In earlier work, we showed that extrapolations of the Ultra-Violet Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at fixed Lyman continuum photon escape fraction are not able to reproduce the redshift evolution of this background. Here, we employ extrapolations of the high-z LFs to describe the contribution of LBGs to the ionizing photon rate, taking into account the smoothing of the baryonic perturbations, due to the background itself (i.e. the filtering mass), as well as a possible sharp increase of the escape fraction in dwarf galaxies. Under the hypothesis of a dominant contribution of LBGs to cosmic reionization, our results suggest that sources fainter than the current observational limits should be characterised by escape fractions of the order of ~0.1-0.3 (larger than the current estimates for bright galaxies) to account for a z>6 reionization and the measured evolution of cosmic ionizing background, at the same time. The contribution to the background from quasars turns out to be relevant at z<3. Overall, our results support the case for dedicated observations of faint galaxies in the rest-frame UV, in order to better determine their physical properties. Observed escape fractions outside our proposed range bear relevant consequences on the nature of the astrophysical sources responsible for cosmic reionization and/or its buildup process.

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

A crucial challenge in astrophysics over the coming decades will be to understand the origins of supermassive black holes (SMBHs) that lie at the centres of most, if not all, galaxies. The processes responsible for the initial formation of these SMBHs and their early growth via accretion - when they are seen as Active Galactic Nuclei (AGN) - remain unknown. To address this challenge, we must identify low luminosity and obscured z>6 AGNs, which represent the bulk of early SMBH growth. Sensitive X-ray observations are a unique signpost of accretion activity, uncontaminated by star formation processes, which prevent reliable AGN identification at other wavelengths (e.g. optical, infrared). The Athena+ Wide Field Imager will enable X-ray surveys to be carried out two orders of magnitude faster than with Chandra or XMM-Newton, opening a new discovery space and identifying over 400 z>6 AGN, including obscured sources. Athena+ will also play a fundamental role to enhance the scientific return of future multiwavelength facilities that will probe the physical conditions within the host galaxies of early SMBHs, which is vital for understanding how SMBHs form, what fuels their subsequent growth, and to assess their impact on the early Universe. Follow-up of samples of z>6 galaxies with the Athena+ X-ray Integral Field Unit could also reveal the presence of highly obscured AGNs, thanks to the detection of strong iron lines. Thus, Athena+ will enable the first quantitative measurements of the extent and distribution of SMBH accretion in the early Universe.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

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.

We report the results of a study exploring the stellar populations of 13 luminous (L>L*), spectroscopically confirmed, galaxies in the redshift interval 5.5<z<6.5, all with WFC3/IR and IRAC imaging from the HST/CANDELS and Spitzer/SEDS surveys. Based on fitting the observed photometry with SED templates covering a wide range of different star-formation histories, and a self consistent treatment of Lyman-alpha emission, we find that the derived stellar masses lie within the range 10^9 Msun < M*< 10^10 Msun and are robust to within a factor of two. In contrast, we confirm previous reports that the ages of the stellar populations are poorly constrained. Although the best-fitting models for three objects have ages >= 300 Myr, the degeneracies introduced by dust extinction mean that only two of these objects actually require a >300 Myr old stellar population to reproduce the observed photometry. Moreover, when considering only smoothly-varying star-formation histories, we observe a clear tension between the data and models such that a galaxy SED template with an old age is often chosen in order to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break this tension we explore SED fitting with two-component models (burst plus on-going star-formation) and allow for nebular emission. On average, the inclusion of nebular emission leads to lower stellar-mass estimates (median offset 0.18 dex), moderately higher specific star-formation rates, and allows for a wider range of plausible stellar ages. However, based on our SED modelling, we find no strong evidence for extremely young ages in our sample (<50 Myr). Finally, considering all of the different star-formation histories explored, we find that the median best-fitting ages are of the order 200-300 Myr and that the objects with the tightest constraints indicate ages in the range 50-200 Myr (Abridged).

We present the first results of our spectroscopic follow-up of 6.5 < z < 10 candidate galaxies behind clusters of galaxies. We report the spectroscopic confirmation of an intrinsically faint Lyman break galaxy (LBG) identified as a z 850LP-band dropout behind the Bullet Cluster. We detect an emission line at {\lambda} = 9412 {\AA} at >5{\sigma} significance using a 16 hr long exposure with FORS2 VLT. Based on the absence of flux in bluer broadband filters, the blue color of the source, and the absence of additional lines, we identify the line as Ly{\alpha} at z = 6.740 \pm 0.003. The integrated line flux is f = (0.7 \pm 0.1 \pm 0.3) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (the uncertainties are due to random and flux calibration errors, respectively) making it the faintest Ly{\alpha} flux detected at these redshifts. Given the magnification of {\mu} = 3.0 \pm 0.2 the intrinsic (corrected for lensing) flux is f^int = (0.23 \pm 0.03 \pm 0.10 \pm 0.02) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (additional uncertainty due to magnification), which is ~2-3 times fainter than other such measurements in z ~ 7 galaxies. The intrinsic H 160W-band magnitude of the object is m^int(H_160W)=27.57 \pm 0.17, corresponding to 0.5 L* for LBGs at these redshifts. The galaxy is one of the two sub-L* LBG galaxies spectroscopically confirmed at these high redshifts (the other is also a lensed z = 7.045 galaxy), making it a valuable probe for the neutral hydrogen fraction in the early universe.

In this paper we discuss the contribution of different astrophysical sources to the ionization of neutral hydrogen at different redshifts. We critically revise the arguments in favour/against a substantial contribution of Active Galactic Nuclei (AGNs) and/or Lyman Break Galaxies (LBGs) to the reionization of the Universe at z>5. We consider extrapolations of the high-z QSO and LBG luminosity functions and their redshift evolution as well as indirect constraints on the space density of lower luminosity Active Galactic Nuclei based on the galaxy stellar mass function. Since the hypothesis of a reionization due to LBGs alone requires a significant contribution of faint dwarf galaxies and a LyC photon escape fraction (f_esc) of the order of ~20%, in tension with present observational constraints, we examine under which hypothesis AGNs and LBGs may provide a combined relevant contribution to the reionization. We show that a relatively steep faint-end of the AGN luminosity function, consistent with present constraints, provides a relevant (although sub-dominant) contribution, thus allowing us to recover the required ionizing photon rates with f_esc~5% up to z~7. At higher redshifts, we test the case for a luminosity-dependent f_esc scenario and we conclude that, if the observed LBGs are indeed characterized by very low f_esc, values of the order of f_esc~70% are needed for objects below our detection threshold, for this galaxy population to provide a substantial contribution to reionization. Clearly, the study of the properties of faint sources (both AGNs and LBGs) is crucial.

We show the effectiveness of strong lensing in the characterisation of Lyman continuum emission from faint L<~0.1L* star-forming galaxies at redshift >~ 3. Past observations of L>~L* galaxies at redshift >~3 have provided upper limits of the average escape fraction of ionising radiation of fesc~5%. Galaxies with relatively high fesc (>10%) seem to be particularly rare at these luminosities, there is therefore the need to explore fainter limits. Before the advent of giant ground based telescopes, one viable way to probe fesc down to 0.05-0.15L* is to exploit strong lensing magnification. This is investigated with Monte Carlo simulations that take into account the current observational capabilities. Adopting a lensing cross-section of 10 arcmin^2 within which the magnification is higher than 1 (achievable with about 4-5 galaxy clusters), with a U-band survey depth of 30(30.5) (AB, 1-sigma), it is possible to constrain fesc for z~3 star-forming galaxies down to 15(10)% at 3-sigma for L<0.15L* luminosities. This is particularly interesting if fesc increases at fainter luminosities, as predicted from various HI reionization scenarios and radiation transfer modelling. Ongoing observational programs on galaxy clusters are discussed and offer positive prospects for the future, even though from space the HST/WFC3 instrument represents the only option we have to investigate details of the spatial distribution of the Lyman continuum emission arising from z~2-4 galaxies.

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 present the morphological catalog of galaxies in nearby clusters of the WINGS survey (Fasano et al. 2006). The catalog contains a total number of 39923 galaxies, for which we provide the automatic estimates of the morphological type applying the purposely devised tool MORPHOT to the V-band WINGS imaging. For ~3000 galaxies we also provide visual estimates of the morphological types. A substantial part of the paper is devoted to the description of the MORPHOT tool, whose application is limited, at least for the moment, to the WINGS imaging only. The approach of the tool to the automation of morphological classification is a non parametric and fully empiri- cal one. In particular, MORPHOT exploits 21 morphological diagnostics, directly and easily computable from the galaxy image, to provide two independent classifications: one based on a Maximum Likelihood (ML), semi-analytical technique, the other one on a Neural Network (NN) machine. A suitably selected sample of ~1000 visually clas- sified WINGS galaxies is used to calibrate the diagnostics for the ML estimator and as a training set in the NN machine. The final morphological estimator combines the two techniques and proves to be effective both when applied to an additional test sample of ~1000 visually classified WINGS galaxies and when compared with small samples of SDSS galaxies visually classified by Fukugita et al. (2007) and Nair et al. (2010). Finally, besides the galaxy morphology distribution (corrected for field contamination) in the WINGS clusters, we present the ellipticity ({\epsilon}), color (B-V) and Sersic index (n) distributions for different morphological types, as well as the morphological fractions as a function of the clustercentric distance (in units of R200).

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

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

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.

Observations of quasar pairs represent a powerful tool in the study of the intergalactic medium, providing in particular information on the quasar radiative feedback not only along the line of sight, but also in the transverse direction. In this paper, we present the spectra of the four quasar pairs we observed so far with X-shooter, and we discuss two individual spectra affected by the presence of a close quasar in the foreground.

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.

We use ultra-deep ultraviolet VLT/VIMOS intermediate-band and VLT/FORS1 narrow-band imaging in the GOODS Southern field to derive limits on the distribution of the escape fraction (f_esc) of ionizing radiation for L >~ L*(z=3) Lyman Break Galaxies (LBGs) at redshift 3.4--4.5. Only one LBG, at redshift z=3.795, is detected in its Lyman continuum (LyC; S/N~5.5), the highest redshift galaxy currently known with a direct detection. Its ultraviolet morphology is quite compact (R_eff=0.8, kpc physical). Three out of seven AGN are also detected in their LyC, including one at redshift z=3.951 and z850 = 26.1. From stacked data (LBGs) we set an upper limit to the average f_esc in the range 5%--20%, depending on the how the data are selected (e.g., by magnitude and/or redshift). We undertake extensive Monte Carlo simulations that take into account intergalactic attenuation, stellar population synthesis models, dust extinction and photometric noise in order to explore the moments of the distribution of the escaping radiation. Various distributions (exponential, log-normal and Gaussian) are explored. We find that the median f_esc is lower than ~6% with an 84% percentile limit not larger than 20%. If this result remains valid for fainter LBGs down to current observational limits, then the LBG population might be not sufficient to account for the entire photoionization budget at the redshifts considered here, with the exact details dependent upon the assumed ionizing background and QSO contribution thereto. It is possible that f_esc depends on the UV luminosity of the galaxies, with fainter galaxies having higher f_esc, and estimates of f_esc from a sample of faint LBG from the HUDF (i775<28.5) are in broad quantitative agreement with such a scenario.

(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 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 investigate the possibility of contamination by lower-redshift interlopers in the measure of the ionising radiation escaping from high redshift galaxies. Taking advantage of the new ultra-deep VLT/VIMOS U-band number counts in the GOODS-S field,we calculate the expected probability of contamination by low-z interlopers as a function of the U-mag and the image spatial resolution (PSF). Assuming that ground-based observations can not resolve objects lying within a 0.5" radius of each other, then each z>=3 galaxy has a 2.1 and 3.2% chance of foreground contamination, adopting surface density U-band number counts down to 27.5 and 28.5,respectively. Those probabilities increase to 8.5 and 12.6%, assuming 1.0" radius. If applied to the estimates reported in the literature at z~3 for which a Lyman continuum has been observed directly,the probability that at least 1/3 of them are affected by foreground contamination is larger than 50%.From a Monte-Carlo simulation we estimate the median integrated contribution of foreground sources to the Lyman continuum flux (f900). Current estimations from stacked data are >2 sigma of the median integrated pollution by foreground sources.The spatial cross-correlation between the U-band catalog and a sample of galaxies at z>=3.4 in the GOODS-S field,produces a number of U-band detected systems fully consistent with the expected superposition statistics. Indeed, each of them shows the presence of at least one offset contaminant in the ACS images. An exemplary case of a foreground contamination in the HUDF at redshift 3.797 by a foreground blue compact source (U=28.63 +/- 0.2) is reported; if observed with a seeing larger than 0.5" the polluting source would mimic an observed (f1500/f900)_OBS ~ 38, erroneously ascribed to the source at higher redshift. [Abridged]

Massive quiescent galaxies at z>1 have been found to have small physical sizes, hence to be superdense. Several mechanisms, including minor mergers, have been proposed for increasing galaxy sizes from high- to low-z. We search for superdense massive galaxies in the WIde-field Nearby Galaxy-cluster Survey (WINGS) of X-ray selected galaxy clusters at 0.04<z<0.07. We discover a significant population of superdense massive galaxies with masses and sizes comparable to those observed at high redshift. They approximately represent 22% of all cluster galaxies more massive than 3x10^10Msol, are mostly S0 galaxies, have a median effective radius <Re> =1.61+/-0.29kpc, a median Sersic index <n> = 3.0+/-0.6, and very old stellar populations with a median mass-weighted age of 12.1+/-1.3Gyr. We calculate a number density of 2.9x10^-2Mpc^-3 for superdense galaxies in local clusters, and a hard lower limit of 1.3x10^-5Mpc^-3 in the whole comoving volume between z = 0.04 and z = 0.07. We find a relation between mass, effective radius and luminosity-weighted age in our cluster galaxies, which can mimic the claimed evolution of the radius with redshift, if not properly taken into account. We compare our data with spectroscopic high-z surveys and find that -when stellar masses are considered- there is consistency with the local WINGS galaxy sizes out to z~2, while a discrepancy of a factor of 3 exists with the only spectroscopic z>2 study. In contrast, there is strong evidence for a large evolution in radius for the most massive galaxies with M*>4x10^11Msol compared to similarly massive galaxies in WINGS, i.e. the BCGs.

We present the full data set of the VIMOS spectroscopic campaign of the ESO/GOODS program in the CDFS, which complements the FORS2 ESO/GOODS spectroscopic campaign. The GOODS/VIMOS spectroscopic campaign is structured in two separate surveys using two different VIMOS grisms. The VIMOS Low Resolution Blue (LR-Blue) and Medium Resolution (MR) orange grisms have been used to cover different redshift ranges. The LR-Blue campaign is aimed at observing galaxies mainly at 1.8<z<3.5, while the MR campaign mainly aims at galaxies at z<1 and Lyman Break Galaxies (LBGs) at z>3.5. The full GOODS/VIMOS spectroscopic campaign consists of 20 VIMOS masks. This release adds 8 new masks to the previous release (12 masks, Popesso et al. 2009). In total we obtained 5052 spectra, 3634 from the 10 LR-Blue masks and 1418 from the 10 MR masks. A significant fraction of the extracted spectra comes from serendipitously observed sources: ~21% in the LR-Blue and ~16% in the MR masks. We obtained 2242 redshifts in the LR-Blue campaign and 976 in the MR campaign for a total success rate of 62% and 69% respectively, which increases to 66% and 73% if only primary targets are considered. The typical redshift uncertainty is estimated to be ~0.0012 (~255 km/s) for the LR-Blue grism and ~0.00040 (~120 km/s) for the MR grism. By complementing our VIMOS spectroscopic catalog with all existing spectroscopic redshifts publicly available in the CDFS, we compiled a redshift master catalog with 7332 entries, which we used to investigate large scale structures out to z~3.7. We produced stacked spectra of LBGs in a few bins of equivalent width (EW) of the Ly-alpha and found evidence for a lack of bright LBGs with high EW of the Ly-alpha. Finally, we obtained new redshifts for 12 X-ray sources of the CDFS and extended-CDFS.

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.

We present deep imaging in the {\em U} band covering an area of 630 arcmin$^{2}$ centered on the southern field of the Great Observatories Origins Deep Survey (GOODS). The data were obtained with the VIMOS instrument at the ESO Very Large Telescope. The final images reach a magnitude limit $U_{lim} \approx 29.8$ (AB, 1$\sigma$, in a 1\arcsec radius aperture), and have good image quality, with full width at half maximum $\approx 0.8\arcsec$. They are significantly deeper than previous U--band images available for the GOODS fields, and better match the sensitivity of other multi--wavelength GOODS photometry. The deeper U--band data yield significantly improved photometric redshifts, especially in key redshift ranges such as $2<z<4$, and deeper color--selected galaxy samples, e.g., Lyman--break galaxies at $z\approx 3$. We also present the coaddition of archival ESO VIMOS R band data, with $R_{lim} \approx 29$ (AB, 1$\sigma$, 1\arcsec radius aperture), and image quality $\approx 0.75 \arcsec$. We discuss the strategies for the observations and data reduction, and present the first results from the analysis of the coadded images.

It has recently been shown that galaxy formation models within the LambdaCDM cosmology predict that, compared to the observed population, small galaxies (with stellar masses < 10^{11} M_sun) form too early, are too passive since z ~ 3 and host too old stellar populations at z=0. We then expect an overproduction of small galaxies at z > 4 that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the MORGANA galaxy formation model and GRASIL spectro-photometric + radiative transfer code to generate mock catalogues of deep fields observed with HST-ACS. We add observational noise and the effect of Lyman-alpha emission, and perform color-color selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i-dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of z_{850} ~ 27 V-dropouts at z ~ 5. We quantify the properties of these "excess" galaxies and discuss the implications: these galaxies are hosted in dark matter halos with circular velocities in excess of 100 km s^{-1}, and their suppression may require a deep re-thinking of stellar feedback processes taking place in galaxy formation.

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 quantify the evolution of the spiral, S0 and elliptical fractions in galaxy clusters as a function of cluster velocity dispersion ($\sigma$) and X-ray luminosity ($L_X$) using a new database of 72 nearby clusters from the WIde-Field Nearby Galaxy-cluster Survey (WINGS) combined with literature data at $z=0.5-1.2$. Most WINGS clusters have $\sigma$ between 500 and 1100 $\rm km s^{-1}$, and $L_X$ between 0.2 and $5 \times 10^{44} \rm erg/s$. The S0 fraction in clusters is known to increase with time at the expense of the spiral population. We find that the spiral and S0 fractions have evolved more strongly in lower $\sigma$, less massive clusters, while we confirm that the proportion of ellipticals has remained unchanged. Our results demonstrate that morphological evolution since $z=1$ is not confined to massive clusters, but is actually more pronounced in low mass clusters, and therefore must originate either from secular (intrinsic) evolution and/or from environmental mechanisms that act preferentially in low-mass environments, or both in low- and high-mass systems. We also find that the evolution of the spiral fraction perfectly mirrors the evolution of the fraction of star-forming galaxies. Interestingly, at low-z the spiral fraction anticorrelates with $L_X$. Conversely, no correlation is observed with $\sigma$. Given that both $\sigma$ and $L_X$ are tracers of the cluster mass, these results pose a challenge for current scenarios of morphological evolution in clusters.

We report on observations of Lyman-break galaxies (LBGs) selected from the Great Observatories Origins Deep Survey at mean redshift z~4, 5 and 6, obtained with FORS2 spectrograph at the ESO-VLT.This program has yielded spectroscopic identifications for 114 galaxies of which 51 are at z~4, 31 at z~5, and 32 at z~6. The adopted selection criteria are effective, identifying galaxies at the expected redshift with minimal foreground contamination. Once selection effects are properly accounted for, the rest-frame UV spectra of the higher-redshift LBGs appear to be similar to their counterparts at z~3. As at z~3, LBGs at z~4 and z~5 show Lya both in emission and in absorption; when in absorption, strong interstellar lines and red UV spectra are observed, a fact also measured at z~2 and 3. At z~6, sensitivity issues bias our sample towards galaxies with Lya in emission; nevertheless, these spectra appear to be similar to their lower-redshift counterparts. As in other studies at similar redshifts, we find clear evidence that brighter LBGs tend to have weaker Lya emission lines. At fixed rest-frame UV luminosity, the equivalent width of the Lya emission line is larger at higher redshifts. At all redshifts where the measurements can be reliably made, the redshift of the Lya emission line turns out to be larger than that of the interstellar absorption lines,with a median velocity difference DeltaV~400 km/s at z~4 and 5, consistent with results at lower redshifts. This shows that powerful, large-scale winds are common at high redshift.In general,there is no strong correlation between the morphology of the UV light and the spectroscopic properties. However, galaxies with deep interstellar absorption lines and strong Lya absorption appear to be more diffuse than galaxies with Lya in emission.(abridged)

We investigate the properties of a clump-cluster galaxy at redshift 1.57. The morphology of this galaxy is dominated by eight star-forming clumps in optical observations, and has photometric properties typical of most clump-cluster and chain galaxies. Its complex asymmetrical morphology has led to the suggestion that this system is a group merger of several initially separate proto-galaxies. We performed H_alpha integral field spectroscopy of this system using SINFONI on VLT UT4. These observations reveal a large-scale velocity gradient throughout the system, but with large local kinematic disturbances. Using a numerical model of gas-rich disk fragmentation, we find that clump interactions and migration can account for the observed disturbed rotation. On the other hand, the global rotation would not be expected for a multiply merging system. We further find that this system follows the stellar mass vs. metallicity, star formation rate and size relations expected for a disk at this redshift, and exhibits a disk-like radial metallicity gradient, so that the scenario of internal disk fragmentation is the most likely one. A red and metallic central concentration appears to be a bulge in this proto-spiral clumpy galaxy. A chain galaxy at redshift 2.07 in the same field also shows disk-like rotation. Such systems are likely progenitors of the present-day bright spiral galaxies, forming their exponential disks through clump migration and disruption and fueling their bulges. Our present results show that disturbed morphologies and kinematics are not necessarily signs of galaxy mergers and interactions, and can instead result from the internal evolution of primordial disks.

We present the first results of the VIsible Multiobject Spectrograph (VIMOS) ESO/GOODS program of spectroscopy of faint galaxies in the Chandra Deep Field South (CDF-S). The program complements the FORS2 ESO/GOODS campaign. 3312 spectra have been obtained in service mode with VIMOS at the ESO/VLT UT3. The VIMOS LR-Blue and MR grisms have been used to cover different redshift ranges. Galaxies at 1.8 < z < 3.5 have been observed in the GOODS VIMOS-LR-Blue campaign. Galaxies at z < 1 and Lyman Break Galaxies at z > 3.5 have been observed in the VIMOS MR survey. Here we report results for the first 6 masks (out of 10 total) that have been analyzed from each of the LR-Blue and MR grisms. Spectra of 2344 spectra have been extracted from these 6 LR-Blue masks and 968 from 6 MR masks. 33% of the LR-Blue and 18% of the MR spectra are serendipitous observations. We obtained 1481 redshifts in the LR-Blue campaign and 656 in the MR campaign for a total success rate of 63% and 68%, respectively, which increase to 70% and 75% when only the primary targets are considered. By complementing our VIMOS spectroscopic catalog with all existing spectroscopic redshifts publicly available in the CDF-S, we created a redshift master catalog. By comparing this redshift compilation with different photometric redshift catalogs we estimate the completeness level of the CDF-S spectroscopic coverage in several redshift bins. The completeness level is very high, > 60%, at z < 3.5, and it is very uncertain at higher redshift. The master catalog has been used also to estimate completeness and contamination levels of different galaxy photometric selection techniques, such as the BzK, the so called 'sub'-U-dropout and the drop-out methods and to identify large scale structures in the field.

CODEX and ESPRESSO are concepts for ultra-stable, high-resolution spectrographs at the E-ELT and VLT, respectively. Both instruments are well motivated by distinct sets of science drivers. However, ESPRESSO will also be a stepping stone towards CODEX both in a scientific as well as in a technical sense. Here we discuss this role of ESPRESSO with respect to one of the most exciting CODEX science cases, i.e. the dynamical determination of the cosmic expansion history.

The redshifts of all cosmologically distant sources are expected to experience a small, systematic drift as a function of time due to the evolution of the Universe's expansion rate. A measurement of this effect would represent a direct and entirely model-independent determination of the expansion history of the Universe over a redshift range that is inaccessible to other methods. Here we investigate the impact of the next generation of Extremely Large Telescopes on the feasibility of detecting and characterising the cosmological redshift drift. We consider the Lyman alpha forest in the redshift range 2 < z < 5 and other absorption lines in the spectra of high redshift QSOs as the most suitable targets for a redshift drift experiment. Assuming photon-noise limited observations and using extensive Monte Carlo simulations we determine the accuracy to which the redshift drift can be measured from the Ly alpha forest as a function of signal-to-noise and redshift. Based on this relation and using the brightness and redshift distributions of known QSOs we find that a 42-m telescope is capable of unambiguously detecting the redshift drift over a period of ~20 yr using 4000 h of observing time. Such an experiment would provide independent evidence for the existence of dark energy without assuming spatial flatness, using any other cosmological constraints or making any other astrophysical assumption.

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.

CODEX, a high resolution, super-stable spectrograph to be fed by the E-ELT, the most powerful telescope ever conceived, will for the first time provide the possibility of directly measuring the change of the expansion rate of the Universe with time and much more, from the variability of fundamental constants to the search for other earths. A study for the implementation at the VLT of a precursor of CODEX, dubbed ESPRESSO, is presently carried out by a collaboration including ESO, IAC, INAF, IoA Cambridge and Observatoire de Geneve. The present talk is focused on the cosmological aspects of the experiment.

Aims. We present the full data set of the spectroscopic campaign of the ESO/GOODS program in the GOODS-South field, obtained with the FORS2 spectrograph at the ESO/VLT. Method. Objects were selected as candidates for VLT/FORS2 observations primarily based on the expectation that the detection and measurement of their spectral features would benefit from the high throughput and spectral resolution of FORS2. The reliability of the redshift estimates is assessed using the redshift-magnitude and color-redshift diagrams, and comparing the results with public data. Results. Including the third part of the spectroscopic campaign (12 masks) to the previous work (26 masks, Vanzella et al. 2005, 2006), 1715 spectra of 1225 individual targets have been analyzed. The actual spectroscopic catalog provides 887 redshift determinations. The typical redshift uncertainty is estimated to be sigma(z) ~ 0.001. Galaxies have been selected adopting different color criteria and using photometric redshifts. The resulting redshift distribution typically spans two domains: from z=0.5 to 2 and z=3.5 to 6.3. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/

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 present the first spatial clustering measurements of z~1, 24um-selected, star forming galaxies in the Great Observatories Origins Deep Survey (GOODS). The sample under investigation includes 495 objects in GOODS-South and 811 objects in GOODS-North selected down to flux densities of f_24>20 uJy and z_AB<23.5 mag, for which spectroscopic redshifts are available. The median redshift, IR luminosity and star formation rate (SFR) of the samples are z~0.8, L_IR~4.4 x 10^10 L_sun, and SFR~7.6 M_sun/yr, respectively. We measure the projected correlation function w(r_p) on scales of r_p=0.06-10 h^-1 Mpc, from which we derive a best fit comoving correlation length of r_0 = 4.0 +- 0.4 h^-1 Mpc and slope of gamma=1.5 +- 0.1 for the whole f_24>20uJy sample after combining the two fields. We find indications of a larger correlation length for objects of higher luminosity, with Luminous Infrared Galaxies (LIRGs, L_IR>10^11 L_sun) reaching r_0~5.1 h^-1 Mpc. This would imply that galaxies with larger SFRs are hosted in progressively more massive halos, reaching minimum halo masses of ~3 x 10^12 M_sun for LIRGs. We compare our measurements with the predictions from semi-analytic models based on the Millennium simulation. The variance in the models is used to estimate the errors in our GOODS clustering measurements, which are dominated by cosmic variance. The measurements from the two GOODS fields are found to be consistent within the errors. On scales of the GOODS fields, the real sources appear more strongly clustered than objects in the Millennium-simulation based catalogs, if the selection function is applied consistently. This suggests that star formation at z~0.5-1 is being hosted in more massive halos and denser environments than currently predicted by galaxy formation models.[truncated]

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 study the relationship between the local environment of galaxies and their star formation rate (SFR) in the Great Observatories Origins Deep Survey, GOODS, at z~1, from ultradeep imaging at 24 microns with the MIPS camera onboard Spitzer. We show that the star formation-density relation observed locally was reversed at z~1: the average SFR of an individual galaxy increased with local galaxy density when the universe was less than half its present age. Hierarchical galaxy formation models (simulated lightcones from the Millennium model) predicted such a reversal to occur only at only at earlier epochs (z>2) and at a lower level. We present a remarkable structure at z~1.016, containing X-ray traced galaxy concentrations, which will eventually merge into a Virgo-like cluster. This structure illustrates how the individual SFR of galaxies increases with density at the ~1-2 Mpc scale. The SFR of z~1 galaxies is found to correlate with stellar mass suggesting that mass plays a role in the observed star formation-density trend. However the specific SFR (=SFR/M*) decreases with stellar mass while it increases with galaxy density, which implies that the environment does directly affect the star formation activity of galaxies. Major mergers do not appear to be the unique or even major cause for this effect since nearly half (46%) of the luminous infrared galaxies (LIRGs) at z~1 present the HST-ACS morphology of spirals, while only a third present a clear signature of major mergers. The remaining galaxies are divided into compact (9%) and irregular (14%) galaxies. Moreover, the specific SFR of major mergers is only marginally stronger than that of spirals. Reproducing the SFR-density relation at z ~ 1 is a new challenge for models, requiring a correct balance between mergers and in-situ star formation at early epochs.

We present initial results from our ongoing campaign of spectroscopic identifications of Lyman--break galaxies (LBGs) at z ~ 4, 5 and 6 with FORS2 at the ESO VLT. The sample of LBGs is selected from the HST/ACS images in the GOODS southern field (CDF--S), and here we discuss the spectroscopic properties of 103 LBGs that we have cofirmed in the redshift interval 3.5 < z < 6.3, of which 34 at z>5 and 27 between 5.5 and 6.3. As for LBGs at lower redshifts, we find that LBGs with LyA in emission are characterized by weaker intestellar absorption lines and bluer UV continua than galaxies whose LyA is observed in absorption. No correlation between the LyA properties and the apparent elongation of the UV morphology of the galaxies is observed, suggesting that geometrical effects are not the cause of the spectral differences. A trend of increasing LyA equivalent width with decreasing luminosity is also observed. We also present new measures of the UV luinosity function of LBGs.

We measure the luminosity function of QSOs in the redshift range 3.5 < z < 5.2 for the absolute magnitude interval -21 < M_{145} < -28. Suitable criteria are defined to select faint QSOs in the GOODS fields, checking their effectiveness and completeness in detail. The confirmed sample of faint QSOs is compared with a brighter one derived from the SDSS. Using a Monte-Carlo technique we estimate the properties of the luminosity function. Our results show that models based on pure density evolution show better agreement with observation than models based on pure luminosity evolution, even if a different break magnitude with respect to z ~ 2.1 is required at 3.5 < z < 5.2. According to our modeling a faint-end slope steeper than low-redshift observations is required to reproduce the data, moreover models with a steep bright-end slope score a higher probability than models with a bright-end flattening. Determining the faint-end of the luminosity function at these redshifts provides important constraints on models of the joint evolution of galaxies and AGNs.

Aims: In this work the luminosity function of QSOs is measured in the redshift range 3.5<z<5.2 for the absolute magnitude interval -21<M_{145}<-28. The determination of the faint end of the luminosity function at these redshifts provides important constraints on models of joint evolution of galaxies and AGNs. Methods: We have defined suitable criteria to select faint QSOs in the GOODS fields, checking in detail their effectiveness and completeness. Spectroscopic follow-up of the resulting QSO candidates has been carried out. The confirmed sample of faint QSOs is compared with a brighter one derived from the SDSS. We have used a Monte-Carlo technique to estimate the properties of the luminosity function, checking various parameterizations for its shape and evolution. Results: Models based on Pure Density Evolution show better agreement with observation than models based on Pure Luminosity Evolution. However a different break magnitude with respect to z~2.1 is required at 3.5<z<5.2. Models with a steeper faint end score a higher probability. We do not find any evidence for a flattening of the bright end at redshift z>3.5. Conclusions: The estimated space density evolution of QSOs indicates a suppression of the formation and/or feeding of Supermassive Black Holes at these redshifts. The QSO contribution to the UV background is insufficient to ionize the IGM at 3.5<z<5.2.

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 the second campaign of the ESO/GOODS program of spectroscopy of faint galaxies in the GOODS-South field. Objects were selected as candidates for VLT/FORS2 observations primarily based on the expectation that the detection and measurement of their spectral features would benefit from the high throughput and spectral resolution of FORS2. The reliability of the redshift estimates is assessed using the redshift-magnitude and color-redshift diagrams and comparing the results with public data. 807 spectra of 652 individua targets have been obtained in service mode with the FORS2 spectrograph at the ESO/VLT, providing 501 redshift determinations. The typical redshift uncertainty is estimated to be sigma_z ~ 0.001. Galaxies have been selected adopting three different color criteria and using the photometric redshifts.The resulting redshift distribution typically spans two redshift domains: from z=0.5 to 2 and z=3.5 to 6.2. In particular, 94 B435-,V606-,i775 "dropout" Lyman break galaxies have been observed, yielding redshifts for 65 objects in the interval 3.4<z<6.2. Three sources have been serendipitously discovered in the redshift interval 4.8<z<5.5. Together with the previous release, 930 sources have now been observed and 724 redshift determinations have been carried out. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/. Large scale structures are clearly detected at z=0.666, 0.734, 1.096, 1.221, 1.300, and 1.614. A sample of 34 sources with tilted [OII]3727 emission has been identified, 32 of them in the redshift range 0.9<z<1.5.

We use the GOODS-MUSIC sample, a catalog of ~3000 Ks-selected galaxies based on VLT and HST observation of the GOODS-South field with extended multi-wavelength coverage (from 0.3 to 8 micron) and accurate estimates of the photometric redshifts to select 179 DRGs with J-Ks>1.3 in an area of 135 sq. arcmin. We first show that the J-Ks>1.3 criterion selects a rather heterogeneous sample of galaxies, going from the targeted high-redshift luminous evolved systems, to a significant fraction of lower redshift (1<z<2) and less luminous dusty starbursts. These low-redshift DRGs are significantly less clustered than higher-z DRGs. With the aid of extreme and simplified theoretical models of clustering evolution we show that it is unlikely that the two samples are drawn from the same population observed at two different stages of evolution. High-z DRGs likely represent the progenitors of the more massive and more luminous galaxies in the local Universe and might mark the regions that will later evolve into structures of intermediate mass, like groups or small galaxy clusters. Low-z DRGs, on the other hand, will likely evolve into slightly less massive field galaxies.

We present a high quality multiwavelength (from 0.3 to 8.0 micron) catalog of the large and deep area in the GOODS Southern Field covered by the deep near-IR observations obtained with the ESO VLT. The catalog is entirely based on public data: in our analysis, we have included the F435W, F606W, F775W and F850LP ACS images, the JHKs VLT data, the Spitzer data provided by IRAC instrument (3.6, 4.5, 5.8 and 8.0 micron), and publicly available U-band data from the 2.2ESO and VLT-VIMOS. We describe in detail the procedures adopted to obtain this multiwavelength catalog. In particular, we developed a specific software for the accurate "PSF-matching" of space and ground-based images of different resolution and depth (ConvPhot), of which we analyse performances and limitations. We have included both z-selected, as well as Ks-selected objects, yielding a unique, self-consistent catalog. The largest fraction of the sample is 90% complete at z~26 or Ks~23.8 (AB scale). Finally, we cross-correlated our data with all the spectroscopic catalogs available to date, assigning a spectroscopic redshift to more than 1000 sources. The final catalog is made up of 14847 objects, at least 72 of which are known stars, 68 are AGNs, and 928 galaxies with spectroscopic redshift (668 galaxies with reliable redshift determination). We applied our photometric redshift code to this data set, and the comparison with the spectroscopic sample shows that the quality of the resulting photometric redshifts is excellent, with an average scatter of only 0.06. The full catalog, which we named GOODS-MUSIC (MUltiwavelength Southern Infrared Catalog), including the spectroscopic information, is made publicly available, together with the software specifically designed to this end.

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

[Abridged] We present the rest-frame Js-band and Ks-band luminosity function of a sample of about 300 galaxies selected in the HDF-S at Ks<23 (Vega). We use calibrated photometric redshift together with spectroscopic redshift for 25% of the sample. The sample has allowed to probe the evolution of the LF in the three redshift bins [0;0.8), [0.8;1.9) and [1.9;4) centered at the median redshift z_m ~ [0.6,1.2,3]. The values of alpha we estimate are consistent with the local value and do not show any trend with redshift. We do not see evidence of evolution from z=0 to z_m ~ 0.6 suggesting that the population of local bright galaxies was already formed at z<0.8. On the contrary, we clearly detect an evolution of the LF to z_m ~ 1.2 characterized by a brightening of M* and by a decline of phi*. To z_m ~ 1.2 M* brightens by about 0.4-0.6 mag and phi* decreases by a factor 2-3. This trend persists, even if at a less extent, down to z_m ~ 3 both in the Js-band and in the Ks-band LF. The decline of the number density of bright galaxies seen at z>0.8 suggests that a significant fraction of them increases their stellar mass at 1<z<2-3 and that they underwent a strong evolution in this redshift range. On the other hand, this implies also that a significant fraction of local bright/massive galaxies was already in place at z>3. Thus, our results suggest that the assembly of high-mass galaxies is spread over a large redshift range and that the increase of their stellar mass has been very efficient also at very high redshift at least for a fraction of them.

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

We present the first results of the ESO/GOODS program of spectroscopy of faint galaxies in the Chandra Deep Field South (CDF-S). 399 spectra of 303 unique targets have been obtained in service mode with the FORS2 spectrograph at the ESO/VLT, providing 234 redshift determinations (the median of the redshift distribution is at 1.04). The typical redshift uncertainty is estimated to be sig(z) ~ 0.001. Galaxies have been color selected in a way that the resulting redshift distribution typically spans from z=0.5 to 2. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/ Large scale structure is clearly detected at z ~ 0.67, 0.73, 1.10 and 1.61. Three Lyman-break galaxies have also been included as targets and are confirmed to have redshifts z=4.800, 4.882 and 5.828. In a few cases, we observe clear [OII]3727 rotation curves, even at the relatively low resolution (R = 860) of the present observations. Assuming that the observed velocity structure is due to dynamically-relaxed rotation, this is an indication of large galactic masses (few times 10^(11) solar masses) at z ~ 1.

The Great Observatories Origins Deep Survey provides significant constraints on the space density of less luminous QSOs at high redshift, which is particularly important to understand the interplay between the formation of galaxies and super-massive black holes and to measure the QSO contribution to the UV ionizing background. We present the results of a search for high-z QSOs, identified in the two GOODS fields on the basis of deep imaging in the optical (with HST) and X-ray (Chandra), and discuss the allowed space density of QSOs in the early universe.

[abr] We report the analysis of three galaxies selected in the Hubble Deep Field South at Ks<22 on the basis of their unusually red near-IR color J-K>3. We have used population synthesis models to constrain their redshifts and their stellar masses. One galaxy is at redshift z_p~2.4 while the other two are at z_p~2.9-3.0. All three galaxies have already assembled a stellar mass of about 10^{11} M_sun at the observed redshift placing the possible merging event of their formation at z>3.5. The inferred mass weighted age of their stellar populations implies that the bulk of the stars formed at z_f>3.5. The resulting co-moving density of M_{stars}>10^{11} M_sun galaxies at <z>~2.7 is rho=1.2(+-0.7)x10^{-4} Mpc^{-3}, about a factor two higher than the predictions of hierarchical models. The comparison with the local density of galaxies implies that the three galaxies must have already formed most of their stellar mass and that they cannot follow an evolution significantly different from a passive aging.The comparison with the density of local L>L* early types (passively evolved galaxies) suggests that their co-moving density cannot decrease by more than a factor 2.5-3 from z=0 to z~3 and that up to 40% of the stellar mass content of bright (L>L*) local early type galaxies was already in place at z>2.5.

A census of massive galaxies at redshift increasingly higher than $z\sim$1 may provide strong constraints on the history of mass assembly and of star formation. Here we report on the analysis of three galaxies selected in the Hubble Deep Field South on the basis of their unusually red near-IR color J-K$\ge$3. These objects result to be massive (M$_{star}\sim 10^{11}$M$_{\odot}$) galaxies at redshift 2.5$<$z$<$3 which cannot follow an evolution significantly different from passive aging. Thus, their counterpart can be found among the EROs (i.e., R-K>5) at 1$<$z$<$1.5, and among the local massive old early-type galaxies at $z=0$. In order to find their counterpart also at the intermediate redshift z$\sim$2, we present a new color selection criterion (J-K$<$3 and I-H$>$4) on the basis of which we identified 4 massive evolved galaxies at 1.5$<z<$2.5. We show an example of a passive evolution scenario in which three massive galaxies at different redshift are apparently those we expect from a massive galaxy fully assembled at z$\sim$3 which evolves passively in time down to $z=0$.

We present a technique for the estimation of photometric redshifts based on feed-forward neural networks. The Multilayer Perceptron (MLP) Artificial Neural Network is used to predict photometric redshifts in the HDF-S from an ultra deep multicolor catalog. Various possible approaches for the training of the neural network are explored, including the deepest and most complete spectroscopic redshift catalog currently available (the Hubble Deep Field North dataset) and models of the spectral energy distribution of galaxies available in the literature. The MLP can be trained on observed data, theoretical data and mixed samples. The prediction of the method is tested on the spectroscopic sample in the HDF-S (44 galaxies). Over the entire redshift range, $0.1<z<3.5$, the agreement between the photometric and spectroscopic redshifts in the HDF-S is good: the training on mixed data produces sigma_z(test) ~ 0.11, showing that model libraries together with observed data provide a sufficiently complete description of the galaxy population. The neural system capability is also tested in a low redshift regime, 0<z<0.4, using the Sloan Digital Sky Survey Data Release One (DR1) spectroscopic sample. The resulting accuracy on 88108 galaxies is sigma_z(test) ~ 0.022. Inputs other than galaxy colors - such as morphology, angular size and surface brightness - may be easily incorporated in the neural network technique. An important feature, in view of the application of the technique to large databases, is the computational speed: in the evaluation phase, redshifts of 10^5 galaxies are estimated in few seconds.

This is the third paper of a series describing the Asiago-ESO/RASS QSO survey (AERQS), a project aimed at the construction of an all-sky statistically well-defined sample of relatively bright QSOs (B<15) at z<0.3. We present here the clustering analysis of the full spectroscopically identified database (392 AGN). The clustering signal at 0.02<z<0.22 is detected at a 3-4 sigma level and its amplitude is measured to be r_0=8.6\pm 2.0 h^{-1} Mpc (in a LambdaCDM model). The comparison with other classes of objects shows that low-redshift QSOs are clustered in a similar way to Radio Galaxies, EROs and early-type galaxies in general, although with a marginally smaller amplitude. The comparison with recent results from the 2QZ shows that the correlation function of QSOs is constant in redshift or marginally increasing toward low redshift. We discuss this behavior with physically motivated models, deriving interesting constraints on the typical mass of the dark matter halos hosting QSOs, M_DMH= 10^{12.7} h^{-1} M_sun (10^{12.0}-10^{13.5}h^{-1} M_sun at 1 sigma confidence level). Finally, we use the clustering data to infer the physical properties of local AGN, obtaining M_BH=2 10^8 h^{-1} M_sun (10^7-3 10^9 h^{-1} M_sun) for the mass of the active black holes, tau_{AGN}= 8 10^6 yr (2 10^{6}-5 10^{7} yr) for their life-time and eta = 0.14 for their efficiency (always for a LambdaCDM model).

We use extensive multi-wavelength photometric data from the Great Observatories Origins Deep Survey (GOODS) to estimate photometric redshifts for a sample of 434 galaxies with spectroscopic redshifts in the Chandra Deep Field South. Using the Bayesian method, which incorporates redshift/magnitude priors, we estimate photometric redshifts for galaxies in the range 18 < R (AB) < 25.5, giving an rms scatter of 0.11. The outlier fraction is < 10%, with the outlier-clipped rms being 0.047. We examine the accuracy of photometric redshifts for several, special sub--classes of objects. The results for extremely red objects are more accurate than those for the sample as a whole, with rms of 0.051 and very few outliers (3%). Photometric redshifts for active galaxies, identified from their X-ray emission, have a dispersion of 0.104, with 10% outlier fraction, similar to that for normal galaxies. Employing a redshift/magnitude prior in this process seems to be crucial in improving the agreement between photometric and spectroscopic redshifts.

We report early results on galaxies at z~6, selected from Hubble Space Telescope imaging for the Great Observatories Origins Deep Survey. Spectroscopy of one object with the Advanced Camera for Surveys grism and from the Keck and VLT observatories a shows a strong continuum break and asymmetric line emission, identified as Lya at z=5.83. We detect only five spatially extended, z~6 candidates with signal-to-noise ratios > 10, two of which have spectroscopic confirmation. This is many fewer than would be expected if galaxies at z=6 had the same luminosity function as those at z=3. There are many fainter candidates, but we expect substantial contamination from foreground interlopers and spurious detections. Our best estimates favor a z=6 galaxy population with fainter luminosities, higher space density, and similar co-moving ultraviolet emissivity to that at z=3, but this depends critically on counts at fluxes fainter than those reliably probed by the current data.

We present a sample of 17 high-redshift (3.5<z<5.2) QSO candidates in the 320 sq.arcmin area of the Great Observatories Origins Deep Survey, selected in the magnitude range 22.45<z_{850}<25.25 using deep imaging with the Advanced Camera for Surveys onboard the Hubble Space Telescope and the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. On the basis of seven spectroscopic and ten photometric redshifts we estimate that the final sample will contain between two and four QSOs with 4<z<5.2. A dearth of high-redshift, moderate-luminosity (M_{145}=~-23) QSOs is observed with respect to predictions based on a) the extrapolation of the z~2.7 luminosity function (LF), according to a pure luminosity evolution calibrated by the results of the Sloan Digital Sky Survey; and b) a constant universal efficiency in the formation of super-massive black holes (SMBHs) in dark-matter halos. Evidence is gathered in favor of a density evolution of the LF at high redshift and of a suppression of the formation or feeding of SMBHs in low-mass halos.

We use a deep K(AB)<25 galaxy sample in the Hubble Deep Field South to trace the evolution of the cosmological stellar mass density from z~ 0.5 to z~3. We find clear evidence for a decrease of the average stellar mass density at high redshift, 2<z<3.2, that is 15^{+25}_{-5}% of the local value, two times higher than what observed in the Hubble Deep Field North. To take into account for the selection effects, we define a homogeneous subsample of galaxies with 10^{10}M_\odot \leq M_* \leq 10^{11}M_\odot: in this sample, the mass density at z>2 is 20^{+20}_{-5} % of the local value. In the mass--limited subsample at z>2, the fraction of passively fading galaxies is at most 25%, although they can contribute up to about 40% of the stellar mass density. On the other hand, star--forming galaxies at z>2 form stars with an average specific rate at least <\dot M/M_*> ~4 x10^{-10} yr$^{-1}$, 3 times higher than the z<~1 value. This implies that UV bright star--forming galaxies are substancial contributors to the rise of the stellar mass density with cosmic time. Although these results are globally consistent with $\Lambda$--CDM scenarios, the present rendition of semi analytic models fails to match the stellar mass density produced by more massive galaxies present at z>2.

We present an estimate of the cosmological evolution of the field galaxy luminosity function (LF) in the rest frame 4400 Angstrom B -band up to redshift z=3.5. To this purpose, we use a composite sample of 1541 I--selected galaxies selected down to I_(AB)=27.2 and 138 galaxies selected down to K_(AB)=25 from ground-based and HST multicolor surveys, most notably the new deep JHK images in the Hubble Deep Field South (HDF-S) taken with the ISAAC instrument at the ESO-VLT telescope. About 21% of the sample has spectroscopic redshifts, and the remaining fraction well calibrated photometric redshifts. The resulting blue LF shows little density evolution at the faint end with respect to the local values, while at the bright end (M_B(AB)<-20) a brightening increasing with redshift is apparent with respect to the local LF. Hierarchical CDM models overpredict the number of faint galaxies by about a factor 3 at z=1. At the bright end the predicted LFs are in reasonable agreement only at low and intermediate redshifts (z=1), but fail to reproduce the pronounced brightening observed in the high redshift (z=2-3) LF. This brightening could mark the epoch where a major star formation activity is present in the galaxy evolution.

We report on low-resolution multi-object spectroscopy of 65 objects from I(AB) ~= 20 to I(AB) ~= 25 in the HDF-S obtained with the VLT Focal Reducer/low dispersion Spectrograph (FORS2). 18 objects belong to the HDF-S proper, i.e. the WFPC2 deep area. 15 high-redshift galaxies with 2.0 < z < 3.5 (10 in the HDF-S proper) have been identified. The spectroscopic redshifts are in good agreement with the photometric ones derived from a chi^{2} minimization technique comparing the observed spectral energy distribution with synthetic libraries and with a new neural network (NN) approach. The dispersion with the former method is sigma_z=0.16 whereas the latter provides sigma_z=0.13. No "catastrophic" difference is encountered. The inferred star formation rates of the individual objects range from tens to a few hundreds of M_{\odot} yr^{-1} and the global star formation rate of the Universe at <z>=~2.4 is estimated to be 0.15 M_{\odot} yr^{-1} Mpc^{-3} with a statistical error of 0.04. Evidence for large scale structure is found with two groups' redshifts observed at z ~= 2.1 and z ~= 2.7 and a pronounced low redshift peak around z ~= 0.58. An elliptical galaxy lensing a background object turns out to be at a redshift z=0.577.

This is the second paper of a series describing the Asiago-ESO/RASS QSO survey, a project aimed at the construction of an all-sky statistically well-defined sample of very bright QSOs (B_J < 15). Such a survey is required to remove the present uncertainties about the properties of the local QSO population and constitutes an homogeneous database for detailed evolutionary studies of AGN. We present here the complete Southern Sample, which comprises 243 bright (12.60 < B_J < 15.13) QSO candidates at high galactic latitudes (|b_{gal}| > 30^{\circ}). The area covered by the survey is 5660 sq. deg. Spectroscopy for the 137 still unidentified objects has been obtained. The total number of AGN turns out to be 111, 63 of which are new identifications. The properties of the selection are discussed. The completeness and the success rate for this survey at the final stage are 63% and 46%, respectively.

The addition of deep near infrared images to the database provided by the HDF-S WFPC2 is essential to monitor the SEDs of the objects on a wide baseline and address a number of key issues including the total stellar content of baryonic mass, the effects of dust extinction, the dependence of morphology on the rest frame wavelength, the photometric redshifts, the detection and nature of extremely red objects (EROs). For these reasons deep near infrared images were obtained with the ISAAC instrument at the ESO VLT in the Js, H and Ks bands reaching, respectively, 23.5, 22.0, 22.0 limiting Vega-magnitude. A multi-color (F300, F450, F606, F814, Js, H, Ks) photometric catalog of the HDF-S has been produced. Photometric redshifts have been generated both fitting templates to the observed SEDs and with neural network techniques. Spectroscopic observations of the 9 candidates with I_AB <24.25 have confirmed all of them to be galaxies with 2<z<3.5. The photometric redshifts for all the galaxies brighter than I_AB< 27.5 have been used to study the evolution of galaxy clustering in the interval 0<z<4.5.

We present an analysis of the evolution of galaxy clustering in the redshift interval 0<z<4.5 in the HDF-S. The HST optical data are combined with infrared ISAAC/VLT observations, and photometric redshifts are used for all the galaxies brighter than I_AB<27.5. The clustering signal is obtained in different redshift bins using two different approaches: a standard one, which uses the best redshift estimate of each object, and a second one, which takes into account the redshift probability function of each object. This second method makes it possible to improve the information in the redshift intervals where contamination from objects with insecure redshifts is important. With both methods, we find that the clustering strength up to z~3.5 in the HDF-S is consistent with the previous results in the HDF-N. While at redshift lower than z~1 the HDF galaxy population is un/anti-biased (b<1) with respect to the underlying dark matter, at high redshift the bias increases up to b~2-3, depending on the cosmological model. These results support previous claims that, at high redshift, galaxies are preferentially located in massive haloes, as predicted by the biased galaxy formation scenario. The impact of cosmic errors on our analyses has been quantified, showing that errors in the clustering measurements in the HDF surveys are indeed dominated by shot-noise in most regimes. Future observations with instruments like the ACS on HST will improve the S/N by at least a factor of two and more detailed analyses of the errors will be required. In fact, pure shot-noise will give a smaller contribution with respect to other sources of errors, such as finite volume effects or non-Poissonian discreteness effects.

We present a deep multicolor (UBVIJsHKs) catalog of galaxies in the HDF-S, based on observations obtained with the HST WFPC2 in 1998 and VLT-ISAAC in 1999. The photometric procedures were tuned to derive a catalog optimized for the estimation of photometric redshifts. In particular we adopted a ``conservative'' detection threshold which resulted in a list of 1611 objects. The behavior of the observed source counts is in general agreement with the result of Casertano et al. (2000) in the HDF-S and Williams et al. (1996) in the HDF-N, while the corresponding counts in the HDF-N provided by Fernandez-Soto et al. (1999) are systematically lower by a factor 1.5 beyond I_AB=26. After correcting for the incompleteness of the source counts, the object surface density at I_AB<27.5 is estimated to be 220 per square arcmin, providing an estimate of the Extragalactic Background Light in the I band consistent with the work of Madau & Pozzetti(2000). The comparison between the median V-I color in the HDF-North and South shows a significant difference around I_AB~26, possibly due to the presence of large scale structure at z~1 in the HDF-N. High-z galaxy candidates (90 U dropout and 17 B dropout) were selected by means of color diagrams, down to a magnitude I_AB=27, with a surface density of (21+-1) and (3.9+-0.9) per square arcmin, respectively. 11 EROs (with (I-K)_AB>2.7) were selected down to K_AB=24, plus 3 objects whose upper limit to the Ks flux is still compatible with the selection criterion. The corresponding surface density of EROs is (2.5+-0.8) per sq.arcmin ((3.2+-0.9) per sq.arcmin if we include the three Ks upper limits). They show a remarkably non-uniform spatial distribution and are classified with roughly equal fractions in the categories of elliptical and starburst galaxies.

We present near-IR (J and Ks) number counts and colors of galaxies detected in deep VLT-ISAAC images centered on the Chandra Deep Field and Hubble Deep Field-South for a total area of 13.6 arcmin$^2$. The limiting surface brightness obtained is Ks$\simeq$22.8 mag/arcsec$^2$ and J$\simeq$24.5 (1$\sigma$) on both fields. A d$log$N/dm relation with a slope of $\sim0.34$ in J and $\sim0.28$ in Ks is found in both fields with no evidence of decline near the magnitude limit. The median J-Ks color of galaxies becomes bluer at magnitudes fainter than Ks$\sim18$, in agreement with the different number counts slope observed in the two bands. We find a fraction ($\le5%$ of the total sample) of sources with color redder than J-Ks=2.3 at magnitudes Ks$>20$. Most of them appear as isolated sources, possibly elliptical or dusty starburst galaxies at redshift $z>2$. The comparison of the observed number counts with models shows that our J-band and Ks-band counts are consistent with the prediction of a model based on a small amount of merging in a $\Omega=1$ cosmology. On the other hand, we fail to reproduce the observed counts if we do not consider merging independently of the parameters defining the universe.