The bright star $\pi$ Men was chosen as the first target for a radial velocity follow-up to test the performance of ESPRESSO, the new high-resolution spectrograph at the ESO's Very-Large Telescope (VLT). The star hosts a multi-planet system (a transiting 4 M$_\oplus$ planet at $\sim$0.07 au, and a sub-stellar companion on a $\sim$2100-day eccentric orbit) which is particularly appealing for a precise multi-technique characterization. With the new ESPRESSO observations, that cover a time span of 200 days, we aim to improve the precision and accuracy of the planet parameters and search for additional low-mass companions. We also take advantage of new photometric transits of $\pi$ Men c observed by TESS over a time span that overlaps with that of the ESPRESSO follow-up campaign. We analyse the enlarged spectroscopic and photometric datasets and compare the results to those in the literature. We further characterize the system by means of absolute astrometry with Hipparcos and Gaia. We used the spectra of ESPRESSO for an independent determination of the stellar fundamental parameters. We present a precise characterization of the planetary system around $\pi$ Men. The ESPRESSO radial velocities alone (with typical uncertainty of 10 cm/s) allow for a precise retrieval of the Doppler signal induced by $\pi$ Men c. The residuals show an RMS of 1.2 m/s, and we can exclude companions with a minimum mass less than $\sim$2 M$_\oplus$ within the orbit of $\pi$ Men c). We improve the ephemeris of $\pi$ Men c using 18 additional TESS transits, and in combination with the astrometric measurements, we determine the inclination of the orbital plane of $\pi$ Men b with high precision ($i_{b}=45.8^{+1.4}_{-1.1}$ deg). This leads to the precise measurement of its absolute mass $m_{b}=14.1^{+0.5}_{-0.4}$ M$_{Jup}$, and shows that the planetary orbital planes are highly misaligned.

Using the GAEA semi-analytic model, we analyse the connection between Damped Lyman-$\alpha$ systems (DLAs) and HI in galaxies. Our state-of-the-art semi-analytic model is tuned to reproduce the local galaxy HI mass function, and that also reproduces other important galaxy properties, including the galaxy mass - gas metallicity relation. To produce catalogs of simulated DLAs we throw $10^5$ random lines of sight in a composite simulated volume: dark matter haloes with log$(\frac{M_{200}}{ M_{\odot}}) \geq 11.5$ are extracted from the Millennium Simulation, while for $9.2 \leq \log(\frac{M_{200}}{ M_{\odot}})<11.5$ we use the Millennium II, and for $8 \leq \log(\frac{M_{200}}{M_{\odot}}) < 9.2$ a halo occupation distribution model. At $2 < z < 3$, where observational data are more accurate, our fiducial model predicts the correct shape of the column density distribution function, but its normalization falls short of the observations, with the discrepancy increasing at higher redshift. The agreement with observations is significantly improved increasing both the HI masses and the disk radii of model galaxies by a factor 2, as implemented 'a posteriori' in our $2M-2R$ model. In the redshift range of interest, haloes with $M_{200} \geq {10}^{11} M_{\odot}$ give the major contribution to $\Omega_{\rm DLA}$, and the typical DLA host halo mass is $\sim {10}^{11} M _{\odot}$. The simulated DLA metallicity distribution is in relatively good agreement with observations, but our model predicts an excess of DLAs at low metallicities. Our results suggest possible improvements for the adopted modelling of the filtering mass and metal ejection in low-mass haloes.

We characterized the transiting planetary system orbiting the G2V star K2-38 using the new-generation echelle spectrograph ESPRESSO. We carried out a photometric analysis of the available K2 photometric light curve of this star to measure the radius of its two known planets. Using 43 ESPRESSO high-precision radial velocity measurements taken over the course of 8 months along with the 14 previously published HIRES RV measurements, we modeled the orbits of the two planets through a MCMC analysis, significantly improving their mass measurements. Using ESPRESSO spectra, we derived the stellar parameters, $T_{\rm eff}$=5731$\pm$66, $\log g$=4.38$\pm$0.11~dex, and $[Fe/H]$=0.26$\pm$0.05~dex, and thus the mass and radius of K2-38, $M_{\star}$=1.03 $^{+0.04}_{-0.02}$~M$_{\oplus}$ and $R_{\star}$=1.06 $^{+0.09}_{-0.06}$~R$_{\oplus}$. We determine new values for the planetary properties of both planets. We characterize K2-38b as a super-Earth with $R_{\rm P}$=1.54$\pm$0.14~R$_{\rm \oplus}$ and $M_{\rm p}$=7.3$^{+1.1}_{-1.0}$~M$_{\oplus}$, and K2-38c as a sub-Neptune with $R_{\rm P}$=2.29$\pm$0.26~R$_{\rm \oplus}$ and $M_{\rm p}$=8.3$^{+1.3}_{-1.3}$~M$_{\oplus}$. We derived a mean density of $\rho_{\rm p}$=11.0$^{+4.1}_{-2.8}$~g cm$^{-3}$ for K2-38b and $\rho_{\rm p}$=3.8$^{+1.8}_{-1.1}$~g~cm$^{-3}$ for K2-38c, confirming K2-38b as one of the densest planets known to date. The best description for the composition of K2-38b comes from an iron-rich Mercury-like model, while K2-38c is better described by a rocky model with a H2 envelope. The maximum collision stripping boundary shows how giant impacts could be the cause for the high density of K2-38b. The irradiation received by each planet places them on opposite sides of the radius valley. We find evidence of a long-period signal in the radial velocity time-series whose origin could be linked to a 0.25-3~M$_{\rm J}$ planet or stellar activity.

We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} M_\odot$) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at $< 10^{10} M_\odot$. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 M_\odot$ are already in a self-regulated state as inferred from observations.

It is a widespread opinion that hydrogen reionization is mainly driven by primeval star-forming galaxies, with a minor role of high-z active galactic nuclei. Recent observations, however, challenge this notion, indicating a number of issues related to a galaxy-driven reionization scenario. We provide here an updated assessment of the space density of relatively faint (M1450~-22.5) AGNs at zspec~5.5 in order to improve the estimate of the photo-ionization rate contribution from accreting super massive black holes. Exploiting deep UV rest-frame ground-based spectra collected at the Very Large Telescope on the CANDELS/GOODS-South field and deep Chandra X-ray images in the CANDELS/GOODS-North and EGS areas, we find two relatively bright (M1450~-22.5) AGNs at zspec~5.5. We derive an AGN space density of Phi=1.29x10^-6 cMpc^-3 at z~5.5 and M1450~-22.5 by simply dividing their observed number by the cosmological volume in the range 5.0<z<6.1. Our estimate does not consider corrections for incompleteness, therefore it represents a lower limit, although uncertainties due to cosmic variance can still be significant. This value supports a high space density of AGNs at z>5, in contrast with previous claims mostly based on standard color selection, possibly affected by significant incompleteness. Our estimate for the AGN photo-ionization rate at z~5.5 is in agreement with the observed values at similar redshifts, which are needed to keep the intergalactic medium highly ionized. Upcoming JWST and giant ground based telescopes observations will improve the study of high-z AGNs and their contribution to the reionization of the Universe.

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.

Ultra-hot giant exoplanets receive thousands of times Earth's insolation. Their high-temperature atmospheres (>2,000 K) are ideal laboratories for studying extreme planetary climates and chemistry. Daysides are predicted to be cloud-free, dominated by atomic species and substantially hotter than nightsides. Atoms are expected to recombine into molecules over the nightside, resulting in different day-night chemistry. While metallic elements and a large temperature contrast have been observed, no chemical gradient has been measured across the surface of such an exoplanet. Different atmospheric chemistry between the day-to-night ("evening") and night-to-day ("morning") terminators could, however, be revealed as an asymmetric absorption signature during transit. Here, we report the detection of an asymmetric atmospheric signature in the ultra-hot exoplanet WASP-76b. We spectrally and temporally resolve this signature thanks to the combination of high-dispersion spectroscopy with a large photon-collecting area. The absorption signal, attributed to neutral iron, is blueshifted by -11+/-0.7 km s-1 on the trailing limb, which can be explained by a combination of planetary rotation and wind blowing from the hot dayside. In contrast, no signal arises from the nightside close to the morning terminator, showing that atomic iron is not absorbing starlight there. Iron must thus condense during its journey across the nightside.

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.

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

One of the major challenges in observational cosmology is related to the redshift evolution of the average hydrogen ionization in the Universe. In order to probe the ionization level of the IGM and the ionization capabilities of bright QSOs at z=4, we have selected a sample of 2508 QSOs from the SDSS survey (DR14) at 3.6<z<4.6 and -29.0<M1450<-26.0. Starting from UV/optical rest-frame spectra of the whole QSO sample, we estimate the escape fraction and free path individually for each of them. We calculate the Lyman Continuum (LyC) escape fraction as the flux ratio blue-ward (~900 A rest-frame) and red-ward (~930 A rest-frame) of the Lyman limit. We then obtain the probability distribution function (PDF) of the individual free paths of the QSOs in the sample and study its evolution in luminosity and redshift, comparing our results with the ones found in literature so far. We find a lower limit to the mean LyC escape fraction of 0.49, in agreement with the values obtained for both brighter and fainter sources at the same redshift. We show that the free paths of ionizing photons are characterized by a skewed distribution function peaked at low values, with an average of ~49-59 proper Mpc at z~4, after excluding possible associated absorbers. This value is larger than the one obtained at the same redshift by many authors in the literature using different techniques. Finally, we also find that the redshift evolution of this parameter results to be possibly milder than previously thought. Our new determination of the mean free path at z~4 implies that previous estimates of the HI photo-ionization rate Gamma_HI available in the literature should be corrected by a factor of 1.2-1.7. These results have important implications if extrapolated at the epoch of reionization.

We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight, including two bright z ~ 3.2 quasars with close separation (~ 500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, log(Mh/Msun) ~ 11 - 13.5. For four of the groups, we find associated MgII absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ~ 480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionised gas seen in emission in lower-redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross section of cool gas and leading to the high fraction of strong MgII absorbers that we detect.

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

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

One of the main scientific drivers for ESPRESSO,\'Echelle SPectrograph, is the detection and characterization of Earth-class exoplanets. With this goal in mind, the ESPRESSO Guaranteed Time Observations (GTO) Catalog identifies the best target stars for a blind search for the radial velocity (RV) signals caused by Earth-class exoplanets. Using the most complete stellar catalogs available, we screened for the most suitable G, K, and M dwarf stars for the detection of Earth-class exoplanets with ESPRESSO. For most of the stars, we then gathered high-resolution spectra from new observations or from archival data. We used these spectra to spectroscopically investigate the existence of any stellar binaries, both bound or background stars. We derived the activity level using chromospheric activity indexes using $log(R'_{HK})$, as well as the projected rotational velocity $\textit{v sin i}$. For the cases where planet companions are already known, we also looked at the possibility that additional planets may exist in the host's habitable zone using dynamical arguments. We estimated the spectroscopic contamination level, $\textit{v sin i}$, activity, stellar parameters and chemical abundances for 249 of the most promising targets. Using these data, we selected 45 stars that match our criteria for detectability of a planet like Earth. The stars presented and discussed in this paper constitute the ESPRESSO GTO catalog for the RV blind search for Earth-class planets. They can also be used for any other work requiring a detailed spectroscopic characterization of stars in the solar neighborhood.

We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations.

We report the discovery of Serenity-18, a galaxy at z=5.939 for which we could measure the content of molecular gas, M(H_2)~ 5 x10^9 M_sun, traced by the CO(6-5) emission, together with the metal-poor ([Fe/H]=-3.08 +- 0.12, [Si/H]=-2.86 +- 0.14) gas clump/filament which is possibly feeding its growth. The galaxy has an estimated star formation rate of ~100 M_sun yr^{-1}, implying that it is a typical main sequence galaxy at these redshifts. The metal-poor gas is detected through a damped Lyman-alpha absorber (DLA) observed at a spatial separation of 40 kpc and at the same redshift of Serenity-18, along the line of sight to the quasar SDSS J2310+1855 (z_em = 6.0025). The chemical abundances measured for the damped Lyman-alpha system are in very good agreement with those measured for other DLAs discovered at similar redshifts, indicating an enrichment due to massive PopII stars. The galaxy/Damped system we discovered is a direct observational evidence of the assembly of a galaxy at the edge of the reionization epoch.

We investigate the NV absorption signal along the line of sight of background quasars, in order to test the robustness of the use of this ion as criterion to select intrinsic (i.e. physically related to the quasar host galaxy) narrow absorption lines (NALs). We build composite spectra from a sample of $\sim$ 1000 CIV absorbers, covering the redshift range 2.55 < z < 4.73, identified in 100 individual sight lines from the XQ-100 Legacy Survey. We detect a statistical significant NV absorption signal only within 5000 km s$^{-1}$ of the systemic redshift, z$\rm_{em}$. This absorption trough is $\sim$ 15$\sigma$ when only CIV systems with N(CIV) > 10$^{14}$ cm$^{-2}$ are included in the composite spectrum. This result confirms that NV offers an excellent statistical tool to identify intrinsic systems. We exploit the stacks of 11 different ions to show that the gas in proximity to a quasar exhibits a considerably different ionization state with respect to gas in the transverse direction and intervening gas at large velocity separations from the continuum source. Indeed, we find a dearth of cool gas, as traced by low-ionization species and in particular by MgII, in the proximity of the quasar. We compare our findings with the predictions given by a range of Cloudy ionization models and find that they can be naturally explained by ionization effects of the quasar.

We present ALMA observations of the CO(6-5) and [CII] emission lines and the sub-millimeter continuum of the $z\sim6$ quasi-stellar object (QSO) SDSS J231038.88+185519.7. Compared to previous studies, we have analyzed a synthetic beam that is ten times smaller in angular size, we have achieved ten times better sensitivity in the CO(6-5) line, and two and half times better sensitivity in the [CII] line, enabling us to resolve the molecular gas emission. We obtain a size of the dense molecular gas of $2.9\pm0.5$ kpc, and of $1.4\pm0.2$ kpc for the 91.5 GHz dust continuum. By assuming that CO(6-5) is thermalized, and by adopting a CO--to--$H_2$ conversion factor $\rm \alpha_{CO} = 0.8~ M_{\odot}~K^{-1}~ (km/s)^{-1} ~pc^{2}$, we infer a molecular gas mass of $\rm M(H_2)=(3.2 \pm0.2) \times 10^{10}\rm M_{\odot}$. Assuming that the observed CO velocity gradient is due to an inclined rotating disk, we derive a dynamical mass of $\rm M_{dyn}~sin^2(i) = (2.4\pm0.5) \times 10^{10}~ M_{\odot}$, which is a factor of approximately two smaller than the previously reported estimate based on [CII]. Regarding the central black hole, we provide a new estimate of the black hole mass based on the C~IV emission line detected in the X-SHOOTER/VLT spectrum: $\rm M_{BH}=(1.8\pm 0.5) \times 10^{9}~ M_{\odot}$. We find a molecular gas fraction of $\rm \mu=M(H_2)/M^*\sim4.4$, where $\rm M^*\approx M_{dyn} - M(H_2)-M(BH)$. We derive a ratio $v_{rot}/\sigma \approx 1-2$ suggesting high gas turbulence, outflows/inflows and/or complex kinematics due to a merger event. We estimate a global Toomre parameter $Q\sim 0.2-0.5$, indicating likely cloud fragmentation. We compare, at the same angular resolution, the CO(6-5) and [CII] distributions, finding that dense molecular gas is more centrally concentrated with respect to [CII]. We find that the current BH growth rate is similar to that of its host galaxy.

The ESPRESSO instrument, to be commissioned in the next months at the ESO VLT, is bound to became a landmark in the field of high-resolution optical spectroscopy, both for its ground-breaking science objectives (search for Earth-like exoplanets; measure of a possible variation of fundamental constants) and for its novel approach to data treatment. For the first time for an ESO instrument, scientific information will be extracted in real time by a dedicated Data Analysis Software (DAS), which includes several interactive workflows to handle the typical analysis cases in stellar and QSO spectroscopy. Data analysis tools in the oncoming ELT era will face very demanding requirements from compelling science case, such as the Sandage Test: the need of handling larger data sizes with a higher degree of accuracy, and the possibility to compare observations and simulated data on the fly. To this purpose, we are currently porting the solutions developed for ESPRESSO to a wider framework, integrating the algorithms within a full-fledged set of Python modules. The project, named "Astrocook", is aimed to provide a set of high-level, instrument-agnostic procedures to automatically extract physical information from the data.

The data analysis software (DAS) for VLT ESPRESSO is aimed to set a new benchmark in the treatment of spectroscopic data towards the extremely-large-telescope era, providing carefully designed, fully interactive recipes to take care of complex analysis operations (e.g. radial velocity estimation in stellar spectra, interpretation of the absorption features in quasar spectra). A few months away from the instrument's first light, the DAS is now mature for science validation, with most algorithms already implemented and operational. In this paper, I will showcase the DAS features which are currently employed on high-resolution HARPS and UVES spectra to assess the scientific reliability of the recipes and their range of application. I will give a glimpse on the science that will be possible when ESPRESSO data become available, with a particular focus on the novel approach that has been adopted to simultaneously fit the emission continuum and the absorption lines in the Lyman-alpha forest of quasar spectra.

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.

In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at far-blue visible wavelengths for the foreseeable future. High-efficiency (>20%) observations of the ground UV (300-400 nm) at medium resolving power (R~20,000) are required to address a number of exciting topics in stellar astrophysics, while also providing new insights in extragalactic science. Anticipating strong demand to better exploit this diagnostic-rich wavelength region, we revisit the science case and instrument requirements previously assembled for the CUBES concept for the Very Large Telescope.

We present new measurements of the free-streaming of warm dark matter (WDM) from Lyman-$\alpha$ flux-power spectra. We use data from the medium resolution, intermediate redshift XQ-100 sample observed with the X-shooter spectrograph ($z=3 - 4.2$) and the high-resolution, high-redshift sample used in Viel et al. (2013) obtained with the HIRES/MIKE spectrographs ($z=4.2 - 5.4$). Based on further improved modelling of the dependence of the Lyman-$\alpha$ flux-power spectrum on the free-streaming of dark matter, cosmological parameters, as well as the thermal history of the intergalactic medium (IGM) with hydrodynamical simulations, we obtain the following limits, expressed as the equivalent mass of thermal relic WDM particles. The XQ-100 flux power spectrum alone gives a lower limit of 1.4 keV, the re-analysis of the HIRES/MIKE sample gives 4.1 keV while the combined analysis gives our best and significantly strengthened lower limit of 5.3 keV (all 2$\sigma$ C.L.). The further improvement in the joint analysis is partly due to the fact that the two data sets have different degeneracies between astrophysical and cosmological parameters that are broken when the data sets are combined, and more importantly on chosen priors on the thermal evolution. These results all assume that the temperature evolution of the IGM can be modelled as a power law in redshift. Allowing for a non-smooth evolution of the temperature of the IGM with sudden temperature changes of up to 5000K reduces the lower limit for the combined analysis to 3.5 keV. A WDM with smaller thermal relic masses would require, however, a sudden temperature jump of $5000\,K$ or more in the narrow redshift interval $z=4.6-4.8$, in disagreement with observations of the thermal history based on high-resolution resolution Lyman-$\alpha$ forest data and expectations for photo-heating and cooling in the low density IGM at these redshifts.

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

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.

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)

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.

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

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

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 present the Lyman-$\alpha$ flux power spectrum measurements of the XQ-100 sample of quasar spectra obtained in the context of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high redshift universe with VLT/XSHOOTER". Using $100$ quasar spectra with medium resolution and signal-to-noise ratio we measure the power spectrum over a range of redshifts $z = 3 - 4.2$ and over a range of scales $k = 0.003 - 0.06\,\mathrm{s\,km^{-1}}$. The results agree well with the measurements of the one-dimensional power spectrum found in the literature. The data analysis used in this paper is based on the Fourier transform and has been tested on synthetic data. Systematic and statistical uncertainties of our measurements are estimated, with a total error (statistical and systematic) comparable to the one of the BOSS data in the overlapping range of scales, and smaller by more than $50\%$ for higher redshift bins ($z>3.6$) and small scales ($k > 0.01\,\mathrm{s\,km^{-1}}$). The XQ-100 data set has the unique feature of having signal-to-noise ratios and resolution intermediate between the two data sets that are typically used to perform cosmological studies, i.e. BOSS and high-resolution spectra (e.g. UVES/VLT or HIRES). More importantly, the measured flux power spectra span the high redshift regime which is usually more constraining for structure formation models.

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.

ESPRESSO is a high-resolution-ultra-stable spectrograph for the VLT, whose commissioning will start in 2017. One of its key science goals is to test the stability of nature's fundamental couplings with unprecedented accuracy and control of possible systematics. A total of 27 nights of the ESPRESSO Consortium's guaranteed time observations (GTO) will be spent in testing the stability of the fine-structure constant and other fundamental couplings. A set of 14 priority optimal targets have been selected for the GTO period. Here we briefly discuss the criteria underlying this selection and describe the selected targets, and then present detailed forecasts of the impact of these measurements on fundamental physics and cosmology, focusing on dark energy constraints and using future supernova type Ia surveys as a comparison point. We show how canonical reconstructions of the dark energy equation of state are improved by the extended redshift range enabled by these spectroscopic measurements, and also quantify additional improvements foreseen for a future ELT-HIRES instrument.

At low densities the standard ionisation history of the intergalactic medium (IGM) predicts a decreasing temperature of the IGM with decreasing density once hydrogen (and helium) reionisation is complete. Heating the high-redshift, low-density IGM above the temperature expected from photo-heating is difficult, and previous claims of high/rising temperatures in low density regions of the Universe based on the probability density function (PDF) of the opacity in Lyman-$\alpha$ forest data at $2<z<4$ have been met with considerable scepticism, particularly since they appear to be in tension with other constraints on the temperature-density relation (TDR). We utilize here an ultra-high signal-to-noise spectrum of the QSO HE0940-1050 and a novel technique to study the low opacity part of the PDF. We show that there is indeed evidence (at 90% confidence level) that a significant volume fraction of the under-dense regions at $z \sim 3$ has temperatures as high or higher than those at densities comparable to the mean and above. We further demonstrate that this conclusion is nevertheless consistent with measurements of a slope of the TDR in over-dense regions that imply a decreasing temperature with decreasing density, as expected if photo-heating of ionised hydrogen is the dominant heating process. We briefly discuss implications of our findings for the need to invoke either spatial temperature fluctuations, as expected during helium reionization, or additional processes that heat a significant volume fraction of the low-density IGM.

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.

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

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 .

We describe the execution and data reduction of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high-redshift universe with VLT/XSHOOTER" (hereafter `XQ-100'). XQ-100 has produced and made publicly available an homogeneous and high-quality sample of echelle spectra of 100 QSOs at redshifts z~3.5-4.5 observed with full spectral coverage from 315 to 2500 nm at a resolving power ranging from R~4000 to 7000, depending on wavelength. The median signal-to-noise ratios are 33, 25 and 43, as measured at rest-frame wavelengths 1700, 3000 and 3600 Angstrom, respectively. This paper provides future users of XQ-100 data with the basic statistics of the survey, along with details of target selection, data acquisition and data reduction. The paper accompanies the public release of all data products, including 100 reduced spectra. XQ-100 is the largest spectroscopic survey to date of high-redshift QSOs with simultaneous rest-frame UV/optical coverage, and as such enables a wide range of extragalactic research, from cosmology and galaxy evolution to AGN astrophysics.

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 statistically study the physical properties of a sample of narrow absorption line (NAL) systems looking for empirical evidences to distinguish between intrinsic and intervening NALs without taking into account any a priori definition or velocity cut-off. We analyze the spectra of 100 quasars with 3.5 < z$\rm_{em}$ < 4.5, observed with X-shooter/VLT in the context of the XQ-100 Legacy Survey. We detect a $\sim$ 8 $\sigma$ excess in the number density of absorbers within 10,000 km/s of the quasar emission redshift with respect to the random occurrence of NALs. This excess does not show a dependence on the quasar bolometric luminosity and it is not due to the redshift evolution of NALs. It extends far beyond the standard 5000 km/s cut-off traditionally defined for associated absorption lines. We propose to modify this definition, extending the threshold to 10,000 km/s when also weak absorbers (equivalent width < 0.2 \AA) are considered. We infer NV is the ion that better traces the effects of the quasar ionization field, offering the best statistical tool to identify intrinsic systems. Following this criterion we estimate that the fraction of quasars in our sample hosting an intrinsic NAL system is 33 percent. Lastly, we compare the properties of the material along the quasar line of sight, derived from our sample, with results based on close quasar pairs investigating the transverse direction. We find a deficiency of cool gas (traced by CII) along the line of sight associated with the quasar host galaxy, in contrast with what is observed in the transverse direction.

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.

Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and could thus be a potentially powerful tool to trace cosmic star formation. However, especially at low redshifts (z < 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the impact of the environment on GRB formation. We study host galaxy spectra of the Swift/BAT6 complete sample of 14 z < 1 bright LGRBs. We use the detected nebular emission lines to measure the dust extinction, star formation rate (SFR) and nebular metallicity (Z) of the hosts and supplement the data set with previously measured stellar masses M$_{\star}$. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M$_{\star}$ relations) are compared to samples of field star-forming galaxies.We find that LGRB hosts at z < 1 have on average lower SFRs than if they were direct star-formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to Z $\sim 8.4-8.5$, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are not inconsistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of increased incidence of starburst galaxies in the Swift/BAT6 z < 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevent us from reaching more quantitative conclusions.

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.

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.

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

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

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.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

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.

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

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

(Abridged) We explore new observationally-constrained sub-resolution models of galactic outflows and investigate their impact on the circumgalactic medium (CGM) over redshifts z = 2 - 4. We perform cosmological hydrodynamic simulations, including star formation, chemical enrichment, and four cases of SNe-driven outflows: no wind (NW), an energy-driven constant velocity wind (CW), a radially varying wind (RVWa) where the outflow velocity has a positive correlation with galactocentric distance (r), and a RVW with additional dependence on halo mass (RVWb). Overall, we find that the outflows expel metal-enriched gas away from galaxies, significantly quench star formation, and enrich the CGM. At z = 2, the radial profiles of gas properties around galaxy centers are most sensitive to the choice of the wind model for halo masses (10^9 - 10^11) M_sun. We infer that the RVWb model is similar to the NW case, except that it substantially enriches the CGM: the carbon metallicity (Z_C) is 10 times higher in RVWb than in NW at r > R_200; and the warm gas of 10^4 - 10^5 K and delta < 100 is enriched to 50 times higher than in NW. We also find that the impact of models CW and RVWa are similar, with the following differences. RVWa causes a greater suppression of star formation rate at z < 5, and has a higher fraction of low-density (delta < 10), warm-hot (10^4 - 10^6 K) gas than in CW. Outflows in CW produce a higher and earlier enrichment of some IGM phases than RVWa. By visual inspection, we note that the RVWa model shows more pronounced bipolar outflows and galactic disks. We present fitting formulae for [Z_C-delta] and [Z_C-r], also for the abundance of CIV as a function of r. We predict observational diagnostics to distinguish between different outflow scenarios: Z_C of the CGM gas at r = (30 - 300) kpc/h comoving, and CIV fraction of the inner gas at r < (4 - 5) kpc/h comoving.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

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 present a measurement of the Lyman alpha flux probability distribution function (PDF) measured from a set of eight high resolution quasar spectra with emission redshifts at 3.3 < z < 3.8. We carefully study the effect of metal absorption lines on the shape of the PDF. Metals have a larger impact on the PDF measurements at lower redshift, where there are fewer Lyman alpha absorption lines. This may be explained by an increase in the number of metal lines which are blended with Lyman alpha absorption lines toward higher redshift, but may also be due to the presence of fewer metals in the intergalactic medium with increasing lookback time. We also provide a new measurement of the redshift evolution of the effective optical depth, tau_eff, at 2.8 < z < 3.6, and find no evidence for a deviation from a power law evolution in the log(tau_eff)-log(1+z) plane. The flux PDF measurements are furthermore of interest for studies of the thermal state of the intergalactic medium (IGM) at z ~ 3 . By comparing the PDF to state-of-the-art cosmological hydrodynamical simulations, we place constraints on the temperature of the IGM and compare our results with previous measurements of the PDF at lower redshift. At redshift z=3, our new PDF measurements are consistent with an isothermal temperature-density relation, T=T_0 Delta^{gamma-1}, with a temperature at the mean density of T_0 = 19250 +/- 4800 K and a slope gamma=0.90+/-0.21 (1 sigma uncertainties). In comparison, joint constraints with previous PDF measurements at z<3 favour an inverted (gamma<1) temperature-density relation with T_0=17900 +/- 3500 K and gamma=0.70 +/- 0.12, in broad agreement with previous analyses.

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

[abridged] We study the evolution of the Star Formation Rate Function (SFRF) of massive galaxies over the 0.4<z<1.8 redshift range and its implications for our understanding of the physical processes responsible for galaxy evolution. We use multiwavelength observations included in the GOODS-MUSIC catalogue, which provides a suitable coverage of the spectral region from 0.3 to 24 micron and either spectroscopic or photometric redshifts for each object. Individual SFRs have been obtained by combining UV and 24 micron observations, when the latter were available. For all other sources an "SED fitting" SFR estimate has been considered. We then define a stellar mass limited sample, complete in the Mstar>1.e10 Msun range and determine the SFRF using the 1/Vmax algorithm. We define simulated galaxy catalogues based on three different semi-analytical models of galaxy formation and evolution. We show that the theoretical SFRFs are well described by a double power law functional form and its redshift evolution is approximated with high accuracy by a pure evolution of the typical SFR. We find good agreement between model predictions and the high-SFR end of the SFRF, when the observational errors on the SFR are taken into account. However, the observational SFRF is characterised by a double peaked structure, which is absent in its theoretical counterparts. At z>1.0 the observed SFRF shows a relevant density evolution, which is not reproduced by SAMs, due to the well known overprediction of intermediate mass galaxies at z~2. The agreement at the low-SFR end is poor: all models overpredict the space density of SFR~1 Msun/yr and no model reproduces the double peaked shape of the observational SFRF. If confirmed by deeper IR observations, this discrepancy will provide a key constraint on theoretical modelling of star formation and stellar feedback.

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 report on our serendipitous discovery that the objects Q 01323-4037 and Q 0132-4037, listed in the V\'eron-Cetty & V\'eron catalog (2006) as two different quasars, are actually a quasar and a star. We briefly discuss the origin of the misidentification, and provide a refined measurement of the quasar redshift.

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

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

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 investigate the properties of triply ionized Carbon (CIV) in the Intergalactic Medium using a set of high-resolution and large box-size cosmological hydrodynamic simulations of a $\Lambda$CDM model. We rely on a modification of the GADGET-2 code that self-consistently follows the metal enrichment mechanism by means of a detailed chemical evolution model. We focus on several numerical implementations of galactic feedback: galactic winds in the energy driven and momentum driven prescriptions and Active Galactic Nuclei (AGN) powered by gas accretion onto massive black holes. We extract mock IGM transmission spectra in neutral hydrogen (HI) and CIV and perform Voigt profile fitting. The results are then compared with high-resolution quasar (QSO) spectra obtained with the UVES spectrograph at the VLT and the HIRES spectrograph at Keck. We find that feedback has little impact on statistics related to the neutral hydrogen, while CIV is more affected by galactic winds and/or AGN feedback. When the same analysis is performed over observed and simulated CIV lines, we find reasonables good agreement between data and simulations over the column density range $N_{\rm CIV}=10^{12.5-15}$ cm$^{-2}$. Also the CIV line-widths distribution appears to be in agreement with the observed values, while the HI Doppler parameters, $b_{\rm HI}$, are in general too large showing that the diffuse cosmic web is heated more than what is inferred by observations. The simulation without feedback fails in reproducing the CIV systems at high column densities at all redshift, while the AGN feedback case agrees with observations only at $z<3$, when this form of feedback is particularly effective. We also present scatter plots in the $b-N$ and in the $N_{\rm CIV}-N_{\rm HI}$ planes, showing that there is rough agreement between observations and simulations only when feedback is taken into account.

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 study the three-dimensional distribution of matter at z~2 using high resolution spectra of QSO pairs and simulated spectra drawn from cosmological hydro-dynamical simulations. We present a sample of 15 QSOs, corresponding to 21 baselines of angular separations evenly distributed between ~1 and 14 arcmin, observed with the Ultraviolet and Visual Echelle Spectrograph (UVES) at the European Southern Observatory-Very Large Telescope (ESO-VLT). The observed correlation functions of the transmitted flux in the HI Lya forest transverse to and along the line of sight are in agreement, implying that the distortions in redshift space due to peculiar velocities are relatively small and - within the relatively large error bars - not significant. The clustering signal is significant up to velocity separations of ~300 km/s, corresponding to about 5 h^{-1} comoving Mpc. Compatibility at the 2 sigma level has been found both for the Auto- and Cross-correlation functions and for the set of the Cross correlation coefficients. The analysis focuses in particular on two QSO groups of the sample. Searching for alignments in the redshift space between Lya absorption lines belonging to different lines of sight, it has been possible to discover the presence of a wide HI structures extending over about ten Mpc in comoving space, and give constraints on the sizes of two cosmic under-dense regions in the intergalactic medium.

The cosmic evolution of the metal content of the intergalactic medium puts stringent constraints on the properties of galactic outflows and on the nature of UV background. In this paper, we present a new measure of the redshift evolution of the mass density of C IV, Omega_CIV, in the interval 1.5 < z < 4 based on a sample of more than 1500 C IV lines with column densities 10^{12} < N(CIV) < 10^{15} cm^{-2}. This sample more than doubles the absorption redshift path covered in the range z<2.5 by previous samples. The result shows a significant increase of Omega_CIV towards the lower redshifts at variance with the previously pictured constant behaviour.

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]

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 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 revisit the proximity effect produced by QSOs at redshifts 2.1-3.3 applying the FLO approach (Saitta et al. 2008) to a sample of ~6300 Ly-alpha lines fitted in 21 high resolution, high signal-to-noise spectra. This new technique allows to recover the hydrogen density field from the HI column densities of the lines in the Ly-alpha forest, on the basis of simple assumptions on the physical state of the gas. To minimize the systematic uncertainties that could affect the density recovering in the QSO vicinity, we carefully determined the redshifts of the QSOs in our sample and modelled in detail their spectra to compute the corresponding ionising fluxes. The mean density field obtained from the observed spectra shows a significant over-density in the region within 4 proper Mpc from the QSO position, confirming that QSOs are hosted in high density peaks. The absolute value of rho/<rho> for the peak is uncertain by a factor of ~3, depending on the assumed QSO spectral slope and the minimum HI column density detectable in the spectra. We do not confirm the presence of a significant over-density extending to separations of ~15 proper Mpc from the QSO, claimed in previous works at redshifts <z>=2.5 and 3.8. Our best guess for the UV background ionisation rate based on the IGM mean density recovered by FLO is Gamma_UVB ~ 10^{-12} s^{-1}. However, values of Gamma_UVB ~ 3x10^{-12} s^{-1} could be viable if an inverted temperature-density relation with index alpha=-0.5 is adopted.

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.

[Abridged] We present FLO (From Lines to Over-densities), a new technique to reconstruct the hydrogen density field for the Lya forest lines observed in high resolution QSO spectra. The method is based on the hypothesis that the Lya lines arise in the low to intermediate density intergalactic gas and that the Jeans length is the typical size of the Lya absorbers. The reliability of FLO is tested against mock spectra obtained from cosmological simulations. The recovering algorithm gives satisfactory results in the range from the mean density to over-densities of ~30 and reproduces correctly the correlation function of the density field and the 1D power spectrum on scales between ~20 and 60 comoving Mpc. A sample of Lya forests from 22 high resolution QSO spectra is analysed, covering the redshift range 1.7<z<3.5. For each line of sight, we fit Voigt profiles to the lines of the Lya forest, providing the largest, homogeneous sample of fitted Lya lines ever studied. The line number density evolution with redshift follows a power-law relation: dn/dz=(166 +/- 4) [(1+z)/3.5]^{(2.8 +/- 0.2)} (1 sigma errors). The two-point correlation function of lines shows a signal up to separations of ~2 comoving Mpc; weak lines (log N(HI)<13.8) also show a significant clustering but on smaller scales (r<1.5 comoving Mpc). We estimate with FLO the hydrogen density field toward the 22 observed lines of sight. The redshift distribution of the average densities computed for each QSO is consistent with the cosmic mean hydrogen density in the analysed redshift range. The two-point correlation function and the 1D power spectrum of the delta field are estimated. The correlation function shows clustering signal up to ~4 comoving Mpc.

We have conducted a long slit search for low surface brightness Lyman-alpha emitters at redshift 2.67 < z < 3.75. A 92 hour long exposure with VLT/FORS2 down to a 1-sigma surface brightness detection limit of 8x10^-20 erg/cm2/s/sqarcsec yielded a sample of 27 single line emitters with fluxes of a few times 10^-18 erg/s/cm2. We present arguments that most objects are indeed Lyman-alpha. The large comoving number density, the large covering factor, dN/dz ~ 0.2-1, and the often extended Lyman-alpha emission suggest that the emitters be identified with the elusive host population of damped Lyman-alpha systems (DLAS) and high column density Lyman limit systems. A small inferred star formation rate, perhaps supplanted by cooling radiation, appears to energetically dominate the Lyman-alpha emission, and is consistent with the low metallicity, low dust content, and theoretically inferred low masses of DLAS, and with the relative lack of success of earlier searches for their optical counterparts. (abridged)

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

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

In order to constrain the physical processes that regulate and downsize the AGN population, the predictions of the MOdel for the Rise of GAlaxies aNd Active nuclei (MORGANA) are compared to luminosity functions (LFs) of AGNs in the optical, soft X-ray and hard X-ray bands, to the local BH-bulge mass relation, and to the observed X-ray number counts and background. We also give predictions on the accretion rate of AGNs in units of the Eddington rate and on the BH--bulge relation expected at high redshift. We find that it is possible to reproduce the downsizing of AGNs within the hierarchical LambdaCDM cosmogony, and that the most likely responsible for this downsizing is the stellar kinetic feedback that arises in star-forming bulges as a consequence of the high level of turbulence and leads to a massive removal of cold gas in small elliptical galaxies. At the same time, to obtain good fits to the number of bright quasars we need to require that quasar-triggered galactic winds self-limit the accretion onto BHs. In all cases, the predicted BH--bulge relation steepens considerably with respect to the observed one at bulge masses <10^{11} Msun; this problem is related to a known excess in the predicted number of small bulges, common to most similar models, so that the reproduction of the correct number of faint AGNs is done at the cost of underestimating their BH masses. This highlights an insufficient downsizing of elliptical galaxies, and hints for another feedback mechanism able to act on the compact discs that form and soon merge at high redshift. The results of this paper reinforce the need for direct investigations of the feedback mechanisms in active galaxies, that will be possible with the next generation of astronomical telescopes from sub-mm to X-rays.

We study the three-dimensional distribution of non virialised matter at z~2 using high resolution spectra of QSO pairs and simulated spectra drawn from cosmological hydrodynamical simulations. We have collected the largest sample of QSO pairs ever observed with UVES at the ESO-VLT, with angular separations between ~1 and 14 arcmin. The observed correlation functions of the transmitted flux in the HI Lyman alpha forest along and transverse to the lines of sight are in good agreement implying that the distortions in redshift space due to peculiar velocities are small. The clustering signal is significant up to velocity separations of ~200 km/s, or about 3 h^{-1} comoving Mpc. The regions at lower overdensity (rho/<rho> < 6.5) are still clustered but on smaller scales (Delta v < 100 km/s). The observed and simulated correlation functions are compatible at the 3 sigma level. A better concordance is obtained when only the low overdensity regions are selected for the analysis or when the effective optical depth of the simulated spectra is increased artificially, suggesting a deficiency of strong lines in the simulated spectra. We found that also a lower value of the power-law index of the temperature-density relation for the Lyman alpha forest gas improves the agreement between observed and simulated results. If confirmed, this would be consistent with other observations favouring a late HeII reionization epoch (at z~3). We remark the detection of a significant clustering signal in the cross correlation coefficient at a transverse velocity separation Delta v_{\perp} ~500 km/s whose origin needs further investigation.

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 deep radio observations of the most distant complete quasar sample drawn from the Sloan Digital Sky Survey. Combining our new data with those from literature we obtain a sample which is ~100 per cent complete down to S_1.4GHz = 60 \mu Jy over the redshift range 3.8 < z < 5. The fraction of radio detections is relatively high (~43 per cent), similar to what observed locally in bright optical surveys. Even though the combined radio and optical properties of quasars remain overall unchanged from z ~ 5 to the local Universe, there is some evidence for a slight over-abundance of radio-loud sources at the highest redshifts when compared with the lower-z regime. Exploiting the deep radio VLA observations we present the first attempt to directly derive the radio luminosity function of bright quasars at z ~ 4. The unique depth -- both in radio and optical -- allows us to thoroughly explore the population of optically bright FR~II quasars up to z ~ 5 and opens a window on the behaviour of the brightest FR~I sources. A close investigation of the space density of radio loud quasars also suggests a differential evolution, with the more luminous sources showing a less pronounced cut-off at high z when compared with the less luminous ones.

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

We investigate the 3-D matter distribution at z~2 with high resolution (R ~ 40000) spectra of QSO pairs and groups obtained with the UVES spectrograph at ESO VLT. Our sample is unique for the number density of objects and the variety of separations, between 0.5 and 7 proper Mpc. We compute the real space cross-correlation function of the Lyman-alpha forest transmitted fluxes. There is a significant clustering signal up to ~2 proper Mpc, which is still present when absorption lines with high column density (log N > 13.8) are excluded.

The K20 survey is a near infrared-selected, deep (Ks < 20) redshift survey targeting galaxies in two independent regions of the sky, the CDFS and the q0055-2659 field. The total Ks-selected sample includes 545 objects. Optical spectra for 525 of them have been obtained, providing 501 spectroscopic identifications (including 12 type-1 AGN and 45 stars). This paper describes the final K20 spectroscopic catalogue, along with the technique used to determine redshifts, measure the spectral features and characterize the spectra. The classification of the galaxy spectra has been performed according to a simple parametric recipe that uses the equivalent widths of the two main emission lines (O[II] and Halpha+[N II]) and two continuum indices (the 4000A break index, and a near-UV color index, C(28-39)). We defined three main spectroscopic classes: red early-type galaxies, blue emission-line galaxies and the intermediate galaxies, which show emission lines but a red continuum. More than 95% of the examined galaxies is included in one of these spectral types and a composite spectrum is built for each of the three galaxy classes. The blue emission-line and the early-type galaxies have been divided in redshift bins, and the corresponding composite spectra have been built, in order to investigate the evolution of the spectral properties of the K20 galaxies with redshift.

We have explored the redshift evolution of the luminosity function of red and blue galaxies up to $z=3.5$. This was possible joining a deep I band composite galaxy sample, which includes the spectroscopic K20 sample and the HDFs samples, with the deep $H_{AB}=26$ and $K_{AB}=25$ samples derived from the deep NIR images of the Hubble Deep Fields North and South, respectively. About 30% of the sample has spectroscopic redshifts and the remaining fraction well-calibrated photometric redshifts. This allowed to select and measure galaxies in the rest-frame blue magnitude up to $z\sim 3$ and to derive the redshift evolution of the B-band luminosity function of galaxies separated by their rest-frame $U-V$ color or specific (i.e. per unit mass) star-formation rate. The class separation was derived from passive evolutionary tracks or from their observed bimodal distributions. Both distributions appear bimodal at least up to $z\sim 2$ and the locus of red/early galaxies is clearly identified up to these high redshifts. Both luminosity and density evolutions are needed to describe the cosmological behaviour of the red/early and blue/late populations. The density evolution is greater for the early population with a decrease by one order of magnitude at $z\sim 2-3$ with respect to the value at $z\sim 0.4$. The luminosity densities of the early and late type galaxies with $M_B<-20.6$ appear to have a bifurcation at $z>1$. Indeed while star-forming galaxies slightly increase or keep constant their luminosity density, "early" galaxies decrease in their luminosity density by a factor $\sim 5-6$ from $z\sim 0.4$ to $z\sim 2.5-3$. A comparison with one of the latest versions of the hierarchical CDM models shows a broad agreement with the observed number and luminosity density evolutions of both populations.

The deep X-ray, optical, and far-infrared fields that constitute GOODS are sensitive to obscured AGN (N_H>10^{22} cm^{-2}) at the quasar epoch (z~2-3), as well as to unobscured AGN as distant as z~7. Luminous X-ray emission is a sign of accretion onto a supermassive black hole and thus reveals all but the most heavily obscured AGN. We combine X-ray luminosity functions with appropriate spectral energy distributions for AGN to model the X-ray, optical and far-infrared flux distributions of the X-ray sources in the GOODS fields. A simple model based on the unified paradigm for AGN, with ~3 times as many obscured AGN as unobscured, successfully reproduces the z-band flux distributions measured in the deep HST ACS observations on the GOODS North and South fields. This model is also consistent with the observed spectroscopic and photometric redshift distributions once selection effects are considered. The previously reported discrepancy between observed spectroscopic redshift distributions and the predictions of population synthesis models for the X-ray background can be explained by bias against the most heavily obscured AGN generated both by X-ray observations and the identification of sources via optical spectroscopy. We predict the AGN number counts for Spitzer MIPS 24 um and IRAC 3.6-8 um observations in the GOODS fields, which will verify whether most AGN in the early Universe are obscured in the optical. Such AGN should be very bright far-infrared sources and include some obscured AGN missed even by X-ray observations.

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.

We use the LUQAS sample of 27 high resolution, high signal-to-noise QSO absorption spectra (Kim et al. 2004) and the results from Croft et al. (2002) together with a suite of high-resolution hydro-dynamical simulations run with the GADGET-II code, to infer the linear dark matter power spectrum on scales of \~ 0.3-30 Mpc at z=2.125 and at z=2.72.

The Damped Lyman-alpha systems (DLAs), seen in absorption in the spectrum of quasars, are believed to contain a large fraction of the neutral gas in the Universe. Paradoxically, these systems are more difficult to observe at z_abs<1.7, since they are rare and their HI feature then falls in UV spectra. Rao & Turnshek (2000) pioneered a method based on MgII-selected DLAs, that is absorbers discovered thanks to our knowledge of their MgII feature in optical spectra. We use new observations undertaken at the TNG as well as a careful literature & archival search to build samples of low redshift absorbers classified according to the technique used for their discovery. We successfully recover N(HI) and equivalent widths of FeII 2600, MgII 2796, MgII 2803 and MgII 2852 for a sample of 36 absorbers, 21 of which are MgII-selected. We find that the MgII-selected sample contains a marginally larger fraction of absorbers with log N(HI)>21.0 than seen otherwise at low redshift. If confirmed, this property will in turn affect estimates of Omega_HI which is dominated by the highest HI column densities. We find that log N(HI) does not correlate significantly with metal equivalent widths. Similarly, we find no evidence that gravitational lensing, the fraction of associated systems or redshift evolution affect the absorber samples in a different way. We conclude that the hint of discrepancies in N(HI) distributions most likely arises from small number statistics. Therefore, further observations are required to better clarify the impact of this discrepancy on estimates of Omega_HI at low redshift.

The flux power spectra of the Lyman-alpha forest from a sample of 27 QSOs taken with the high resolution echelle spectrograph UVES on VLT are presented. We find a similar fluctuation amplitude at the peak of the ``3D'' flux power spectrum at k ~ 0.03 (km/sec)^(-1) as the study by Croft et al. (2002), in the same redshift range. The amplitude of the flux power spectrum increases with decreasing redshift if corrected for the increase in the mean flux level as expected if the evolution of the flux power spectrum is sensitive to the gravitational growth of matter density fluctuations. This is in agreement with the findings of McDonald et al. (2000) at larger redshift. The logarithmic slope of the "3D" flux power spectrum, P_F(k), at large scales k < 0.03 (km/sec)^(-1), is 1.4 +- 0.3, i.e. 0.3 shallower than that found by Croft et al. (2002) but consistent within the errors.

We present a detailed analysis of the stellar mass content of galaxies up to z=2.5 in the K20 galaxy sample, that has a 92% spectroscopic completeness and a complete $UBVRIzJK_s$ multicolor coverage. We find that the M/L ratio decreases with redshift: in particular, the average M/L ratio of early type galaxies decreases with $z$, with a scatter that is indicative of a range of star--formation time-scales and redshift of formation. More important, the typical M/L of massive early type galaxies is larger than that of less massive ones, suggesting that their stellar population formed at higher z. The final K20 galaxy sample spans a range of stellar masses from M*=10^9Msun to M*=10^12Msun, with massive galaxies ($M*>10^11Msun) detected up to z~2. We compute the Galaxy Stellar Mass Function at various z, of which we observe only a mild evolution (i.e. by 20-30%) up to z~1. At z>1, the evolution of the GSMF appears to be much faster: at z~2, about 35% of the present day stellar mass in objects with M*~10^11Msun appear to have assembled. We also detect a change in the physical nature of the most massive galaxies, since at z>1 a population of massive star--forming galaxies progressively appears. We finally analyze our results in the framework of Lambda-CDM hierarchical models. First, we show that the large number of massive galaxies detected at high z does not violate any fundamental Lambda-CDM constraint based on the number of massive DM halos. Then, we compare our results with the predictions of renditions of both semianalytic and hydro-dynamical models, that range from severe underestimates to slight overestimates of the observed mass density at z<~2. We discuss how the differences among these models are due to the different implementation of the main physical processes. (Abridged)

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.

We determined C, N and alpha-element relative abundances in the gas surrounding six QSOs at an average redshift of < z > ~ 2.4, by studying six narrow associated absorption systems in UVES high-resolution spectra. We found five systems with a metallicity (measured by C/H) consistent or above the solar value. The ionization structure observed in the associated systems is clearly different from that of the intervening ones, indicating that the associated systems are influenced by the strong UV flux from the QSO. There is a possible correlation (anticorrelation) between [N/C] ([Si/C]) and [C/H] of the studied associated systems, and [N/C] >= 0 when [C/H] >= 0. We have compared these observational results with the predictions of a model simulating the joint evolution of QSOs and their spheroidal hosts. The agreement turns out to be very good, in particular, the case envisaging massive haloes and high star-formation rates recovers both the correlation between [N/C] and [C/H] and the anticorrelation for [Si/C] vs. [C/H]. Narrow associated absorption systems prove to be powerful tracers of the chemical abundances in gas belonging to high redshift spheroidal galaxies. The outflow of this same gas, triggered by the QSO feedback, is probably going to contribute to the early enrichment of the surrounding intergalactic medium. A larger statistics, possibly increasing the number of ionisation stages, chemical elements and the redshift range, would allow us to put firm constraints on detailed chemical evolution models of galaxies at high redshifts.

[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 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 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 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 describe a possible new class of X-ray sources that have robust detections in ultra-deep Chandra data, yet have no detections at all in our deep multi-band GOODS Hubble Space Telescope (HST) ACS images, which represent the highest quality optical imaging obtained to date on these fields. These extreme X-ray / Optical ratio sources ("EXO"s) have values of Fx/Fopt at least an order of magnitude above those generally found for other AGN, even those that are harbored by reddened hosts. We thus infer two possible scenarios: (1) if these sources lie at redshifts z < 6, then their hosts need to be exceedingly underluminous, or more reddened, compared with other known sources; (2) if these sources lie above z ~ 6-7, such that even their Lyman-alpha emission is redshifted out of the bandpass of our ACS z(850) filter, then their optical and X-ray fluxes can be accounted for in terms of relatively normal L* hosts and moderate-luminosity AGN.

We present the galaxy rest-frame near-IR Luminosity Function (LF) and its cosmic evolution to z=1.5 based on a spectroscopic survey of a magnitude limited sample of galaxies with Ks<20 (the K20 survey, Cimatti et al. 2002b). The LFs have been derived in the rest-frame J and Ks bands using 3 z bins (z_mean= 0.5, 1, 1.5) and compared to the local near-IR LF. The faint-end of the LFs is consistent with the local estimates, with no evidence for a change either in the slope or normalization up to z<1.3. Viceversa, the density of luminous galaxies (M_Ks-5logh_70<-25.5) is higher than locally at all z and relatively constant or mildly increasing with z within our sample. The data are consistent with a mild luminosity evolution both in the J and Ks-band up to z=1.5 (DeltaM_J=-0.69+-0.12 and DeltaM_K=-0.54+-0.12 at z=1). Moreover, we find that red and early-type galaxies dominate the bright-end of the LF, and that their number density shows at most a small decrease (<30%) up to z=1, thus suggesting that massive elliptical galaxies were already in place at z=1 and they should have formed their stars and assembled their mass at higher z. There appears to be a correlation of the optical/near-IR colors with near-IR luminosities, the most luminous/massive galaxies being red/old. We find a slow evolution with z of the near-IR comoving luminosity density to z=1.5. Finally, we show that hierarchical models overpredict significantly the density of low luminosity galaxies at z<=1 and underpredict the density of luminous galaxies at z>=1, whereas PLE models are more consistent with the data up to z=1.5. The GIF model (Kaufmann et al. 1999) shows a clear deficiency of red luminous galaxies at z=1 compared to our observations and predicts a decrease of luminous galaxies with z not observed in our sample.

We present the results of a VLT and HST imaging survey aimed at the identification of 4.5<z<6 galaxies. In the VLT data, a set of broad and intermediate band filters has been used to select 13 high--$z$ candidates in a Z(AB)<25 catalog, over an area of about 30 arcmin2. Discrimination against lower redshift interlopers (mainly early--type galaxies at high redshift and cool Galactic stars) has been done combining morphological and spectral classification. This sample has been combined with a deeper I<27.2 sample obtained from the Hubble Deep Field campaigns. The VLT final sample consists of 13 high z candidates, 4 of which are identified with high confidence as z>4.5 galaxies. The resulting integral surface density of the Z<25 candidates at z>4.5z is in the range 0.13-0.44/arcmin^2 and that in the highest redshift bin 5<z<6 is between 0.07-0.13/arcmin^2. In the two HDFs, we identify at I<27.2 25 galaxies in the range 4.5<z<5 and 16 at 5<z<6$, corresponding to surface densities of 3.1 arcmin^-2 and 2 arcmin^-2, respectively. We show that the observed Z<25 UV luminosity density appear to drop by about one order of magnitude from z~3 to z~6. However, if we apply a threshold to obtain an absolute--magnitude limited sample, the UV luminosity density results to be roughly constant up to z~6. We finally show that recent semi-analytic hierarchical models for galaxy formation, while predicting a nearly constant total UV luminosity density up to z~6, under-predict the observed UV luminosity density at Z<25 and over-predicts the I<27.2 one. This behaviour can be understood in term of a poor match to the slope of the UV luminosity function.

We present the ISO far-IR photometry of a complete sub-sample of optically selected bright quasars belonging to two complete surveys selected through multicolour (U,B,V,R,I) techniques. The ISOPHOT camera on board of the ISO Satellite was used to target these quasars at wavelengths of 7.3, 11.5, 60, 100 and 160 micron. Almost two thirds of the objects were detected at least in one ISOPHOT band. The detection rate is independent of the source redshift, very likely due to the negative K-correction of the far-IR thermal emission. More than a half of the optically selected QSOs show significant emission between 4 and 100 micron in the quasar rest-frame. These fluxes have a very likely thermal origin, although in a few objects an additional contribution from a non-thermal component is plausible in the long wavelength bands. In a colour-colour diagram these objects span a wide range of properties from AGN-dominated to ULIRG-like. The far-IR composite spectrum of the quasar population presents a broad far-IR bump between 10 and 30 micron and a sharp drop at wavelengths greater than 100 micron in the quasar restframe. The amount of energy emitted in the far-IR, is on average a few times larger than that emitted in the blue and the ratio L(FIR)/L(B) increases with the bolometric luminosity. Objects with fainter blue magnitudes have larger ratios between the far-IR (wavelengths > 60 micron) fluxes and the blue band flux, which is attributed to extinction by dust around the central source. No relation between the blue absolute magnitude and the dust colour temperature is seen, suggesting that the dominant source of FIR energy could be linked to a concurrent starburst rather than to gravitational energy produced by the central engine.

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.

The AERQS is a project aimed at the construction of an all-sky statistically well-defined sample of local (z < 0.3) and bright (R < 15.5) QSOs. Present uncertainties on the theoretical modeling of the QSO evolution are to be ascribed to the relatively small number of objects observed in this magnitude and redshift domain. We are filling this gap with a sample of 392 AGN covering all the sky at high Galactic latitude: we present here the Clustering analysis and the Luminosity Function of AGN at low redshift. The AERQS sets an important zero point about the properties of the local QSO population and the evolutionary pattern of QSOs between the present epoch and the highest redshift.

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.

We present the redshift distribution of a complete sample of 480 galaxies with Ks<20 distributed over two independent fields covering a total area of 52 arcmin^2. The redshift completeness is 87% and 98% respectively with spectroscopic and high-quality and tested photometric redshifts. The redshift distribution of field galaxies has a median redshift z_{med}\sim 0.80, with \sim 32% and \sim 9% of galaxies at z>1 and z>1.5 respectively. A ``blind'' comparison is made with the predictions of a set of the most recent LambdaCDM hierarchical merging and pure luminosity evolution (PLE) models. The hierarchical merging models overpredict and underpredict the number of galaxies at low-z and high-z respectively, whereas the PLE models match the median redshift and the low-z distribution, still being able to follow the high-z tail of N(z). We briefly discuss the implications of this comparison and the possible origins of the observed discrepancies. We make the redshift distribution publicly available.

The K20 survey is an ESO VLT optical and near-infrared spectroscopic survey aimed at obtaining spectral information and redshifts of a complete sample of about 550 objects to K_s\leq20.0 over two independent fields with a total area of 52 arcmin^2. In this paper we discuss the scientific motivation of such a survey, we describe the photometric and spectroscopic properties of the sample, and we release the $K_s$-band photometric catalog. Extensive simulations showed that the sample is photometrically highly complete to K_s=20. The observed galaxy counts and the R-K_s color distribution are consistent with literature results. We observed spectroscopically 94% of the sample, reaching a spectroscopic redshift identification completeness of 92% to K_s\leq20.0 for the observed targets, and of 87% for the whole sample (i.e. counting also the unobserved targets). Deep spectroscopy was complemented with multi-band deep imaging in order to derive tested and reliable photometric redshifts for the galaxies lacking spectroscopic redshifts. The results show a very good agreement between the spectroscopic and the photometric redshifts with <z_{spe}-z_{phot}>=0.01 and with a dispersion of \sigma_{\Delta z}=0.09. Using both the spectroscopic and the photometric redshifts, we reached an overall redshift completeness of about 98%. The size of the sample, the redshift completeness, the availability of high quality photometric redshifts and multicolor spectral energy distributions make the K20 survey database one of the most complete samples available to date for constraining the currently competing scenarios of galaxy formation and for a variety of other galaxy evolution studies.

We present HST STIS observations of four quasar pairs with redshifts 0.84<z_em<1.56 and angular separation 2-3 arcmin corresponding to ~ 1-1.5 h^-1_50 Mpc transverse proper distance at z ~ 0.9. We study the distribution of velocity differences between nearest neighbor HI Lyman_alpha absorption lines detected in the spectra of adjacent QSOs in order to search for the possible correlation caused by the extent or the clustering properties of the structures traced by the absorption lines over such a scale. The significance of the correlation signal is determined by comparison with Monte-Carlo simulations of spectra with randomly distributed absorption lines. We find an excess of lines with a velocity separation smaller than Delta(V) = 500 km/s significant at the 99.97% level. This clearly shows that the Lyman_alpha forest is correlated on scales larger than 1 h^-1_50 Mpc at z ~ 1. However, out of the 20 detected coincidences within this velocity bin, 12 have Delta(V) > 200 km/s. This probably reflects the fact that the scale probed by our observations is not related to the real size of individual absorbers but rather to large scale correlation. Statistics are too small to conclude about any difference between pairs separated by either 2 or 3 arcmin. A damped Lyman_alpha system is detected at z_abs = 1.2412 toward LBQS0019-0145A with log[N(HI)] ~ 20.5. From the absence of ZnII absorption, we derive a metallicity relative to solar [Zn/H]<-1.75.

We study the transverse clustering properties of high matter density peaks as traced by high column density absorption systems (either Lyman limit systems characterized by N(HI)> 2 x 10^{17} cm^{-2} or CIV systems with W_{r}> 0.5 A) at redshifts between 2 and 3 with UVES spectra of two QSO pairs (UM680/UM681 at 56 arcsec angular separation and Q2344+1228/Q2343+1232 at 5 arcmin angular separation) and a QSO triplet (Q2139-4433/Q2139-4434/Q2138-4427 at 1, 7 and 8 arcmin angular separation). We find 3 damped Ly-alpha systems (N(HI)> 2 x 10^{20} cm^{-2}): 2 coinciding with strong metal systems in the nearby line of sight and 1 matching the emission redshift of the paired QSO; plus 7 Lyman limit systems: 4 forming two matching couples and 3 without a corresponding metal system within ~3000 km/s in the coupled line of sight. In summary, we detect five out of ten matching systems within 1000 km/s, indicating a highly significant overdensity of strong absorption systems over separation lengths from ~1 to 8 h^{-1} Mpc. The observed coincidences could arise in gas due to starburst-driven superwinds associated with a quasar or a galaxy, or gas belonging to large scale structures like filaments or sheets. We also determine chemical abundance ratios for three damped Ly-alpha systems. In particular, for the damped system at z ~ 2.53788 in the spectrum of Q2344+1228, new estimates of the ratios O/Fe, C/Fe are obtained: [C/Fe]<0.06, [O/Fe]<0.2. They indicate that O and C are not over-solar in this system.

Combining a new, increased dataset of 8 QSOs covering the Ly-alpha forest at redshifts 1.5 < z < 3.6 from VLT/UVES observations with previously published results, we have investigated the properties of the Ly-alpha forest at 1.5 < z < 4. We have fitted Voigt profiles to the Ly-alpha absorption lines as in previous studies, and have, for two QSOs, fitted Ly-alpha at z ~ 2.2 and higher order of Lyman lines down to 3050 Angstrom simultaneously. For a given HI column density range, the Ly-alpha forest at 1.5 < z < 4 shows the monotonic evolution, in terms of the line number density and the HI opacity. For higher column density clouds at 10^{14} cm^(-2),there is a variation in the line number density from sightline to sightline at z < 2.5. This variation is stronger for higher column density systems, probably due to more gravitationally evolved structures at lower z. The correlation strength of the step optical depth correlation function shows the strong evolution from z ~ 3.3 to z ~ 2.1, although there is a large scatter along different sightlines. The analyses of the Ly-beta forest at z ~ 2.2 are, in general, in good agreement with those of the Ly-alpha forest.

Currently public ESO data sets pertinent to the CDFS/GOODS field are briefly illustrated along with an indication on how to get access to them. Future ESO plans for complementing the GOODS database with optical/IR imaging and optical spectroscopy are also described.

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.

Long exposure, long-slit spectra have been obtained in the UV/optical bands for two galaxies at z=2.96 and z=3.32 to investigate the fraction of ionizing UV photons escaping from high redshifts galaxies. The two targets are among the brightest galaxies discovered by Steidel and collaborators and they have different properties in terms of Lyman-alpha emission and dust reddening. No significant Lyman continuum emission has been detected. The noise level in the spectra implies an upper limit of f_{rel,esc}\equiv 3 f(900)/f(1500)< 16% for the relative escape fraction of ionizing photons, after correction for absorption by the intervening intergalactic medium. This upper limit is 4 times lower than the previous detection derived from a composite spectrum of 29 Lyman break galaxies at z 3.4. If this value is typical of the escape fraction of the z 3 galaxies, and is added to the expected contribution of the QSO population, the derived UV background is in good agreement with the one derived by the proximity effect.

We compare the 3D clustering of old passively-evolving and dusty star-forming z~1 EROs from the K20 survey. With detailed simulations of clustering, the comoving correlation length of dusty star-forming EROs is constrained to be less than r_0~2.5 Mpc/h. In contrast, the old EROs are much more positively correlated, with 5.5<~r_0/(Mpc/h)<~16, consistent with previous claims for z~1 field early-type galaxies based on analyses of ERO angular clustering. The low level of clustering of dusty star-forming EROs does not support these to be major mergers building up an elliptical galaxy, or typical counterparts of SCUBA sources, but it is instead consistent with the weak clustering of high redshift blue galaxies and of luminous local IRAS galaxies. Current hierarchical merging models can explain the large r_0 for z~1 field early-type galaxies, but fail in matching their high number density and overall old ages.

Using a new, increased dataset of 7 QSOs from VLT/UVES observations combined with one QSO from the literature, the minimum Doppler parameters as a function of neutral hydrogen column density N(HI), b_c(N), of the Lyman-alpha forest has been derived at three redshifts z = 2.1, 3.3 and 3.8. When incompleteness of the number of the observed lines towards lower N(HI) is accounted for, the derived slopes of b_c(N), (Gamma-1), are consistent with no-redshift evolution with an indication of lower value at z = 3.3, while b_c(N) at a fixed column density N(HI) = 10^(13.6) cm^(-2), b_c(13.6), increases as z decreases. Assuming a QSO-dominated UV background, the slope of the equation of state (gamma-1) shows no z-evolution within large uncertainties and the temperature at the mean density, T_0, decreases as z decreases at three redshift ranges. There is a large fluctuation of (Gamma-1) and b_c(13.6) even at the similar redshifts, in particular at z = 3.3 and 3.8. The lower (Gamma-1) and higher b_c(13.6) values at z ~ 3.1 and 3.6 compared to ones at z ~ 3.4 and 3.9 are caused by a lack of lower-N(HI) and lower-b lines, probably due to the HeII reionization. This result suggests that an impact from the HeII reionization on the forest might be mainly on the lower-N(HI) forest. We find some forest clouds with a high ratio of SiIV column density to CIV column density at z < 2.5, although the bulk of the forest clouds shows lower N(SiIV)/N(CIV). This high N(SiIV)/N(CIV) at z < 2.5 suggests that some forest clouds are exposed to a soft UV background.

The remarkable efficiency of the UVES spectrograph at the VLT has made it possible to push high-resolution, high-S/N ground observations of the Ly-a forest down to z~1.5, gaining new insight into the physical conditions of the intergalactic medium and its evolution over more than 90% of the cosmic time. The universal expansion, the UV ionizing background and the gravitational condensation of structures are the driving factors shaping the number density and the column density distribution of the absorbers. A (limited) contribution of UV photons produced by galaxies is found to be important to reproduce the observed evolutionary pattern at very high and low redshift. The Lyman forest contains most of the baryons, at least at z>1.5, and acts as a reservoir for galaxy formation. The typical Doppler parameter at a fixed column density is measured to slightly increase with decreasing redshift, but the inferred temperature at the mean density is increasing with redshift. The signatures of HeII reionization and feedback from the formation of galactic structures have possibly been detected in the Lyman forest.

We present the results of VLT optical spectroscopy of a complete sample of 78 EROs with R-Ks\geq5 over a field of 52 arcmin^2. About 70% of the 45 EROs with Ks\leq19.2 have been spectroscopically identified with old passively evolving and dusty star-forming galaxies at 0.7<z<1.5. The two classes are about equally populated and for each of them we present and discuss the average spectrum. From the old ERO average spectrum and for Z=Z_{\odot} we derive a minimum age of \sim 3 Gyr, corresponding to a formation redshift of z_f \gtsima 2.4. PLE models with such formation redshifts well reproduce the density of old EROs (consistent with being passively evolving ellipticals), whereas the predictions of the current hierarchical merging models are lower than the observed densities by large factors (up to an order of magnitude). From the average spectrum of the star-forming EROs we estimate a substantial dust extinction with E(B-V) \gtsima 0.5. The star formation rates, corrected for the average reddening, suggest a significant contribution from EROs to the cosmic star-formation density at z \sim 1.

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 preliminary results of an investigation of the clustering properties of high matter density peaks between redshift ~2 and ~3, as traced by Lyman limit and Damped Ly-alpha systems in spectra of close QSO pairs and groups.

A sample of 8 QSOs with 1.7 < z_em < 3.7 has been observed with VLT/UVES at a typical resolution 45000 and S/N~40-50. Thanks to the two-arm design of the spectrograph, a remarkable efficiency has been achieved below 400nm and above 800nm, which translates immediately in the possibility of obtaining new results, especially at z<2.5. We report here new insight gained about the evolution of the number density of Ly-a lines, their column density distribution, the ionizing UV background and the cosmic baryon density.

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 have modelled the evolution of the number of Ly_alpha absorbers with redshift, resulting from the evolution of the ionising background and the Hubble expansion. The contribution of quasars (QSOs) and galaxies to the HI-ionising UV background has been estimated. The QSOs emissivity is derived from recent fits of their luminosity function. The galaxy emissivity is computed using a stellar population synthesis model, with a star-formation history scaled on observations of faint galaxies at lambda <= 1500 A. We allow for three values of the fraction of ionising photons that can escape the interstellar medium, f_esc=0.05, 0.1 and 0.4. The Intergalactic Medium is modelled as made of purely-absorbing clouds with the distribution in redshift and column density obtained from QSOs absorption lines. For the adopted values of f_esc, the contribution of galaxies to the ionising UV background is comparable or greater than that of QSOs. Accounting for the contribution of clouds to the UV emission, all models with f_esc <~ 0.1 provide an ionising flux compatible with local and high-z determination, including those with a pure QSOs background. The observed z ~ 1 break in the evolution can be better explained by a dominant contribution from galaxies. We find that models in Lambda-cosmology with Omega_matter=0.3, Omega_Lambda=0.7 describe the flat absorbers evolution for z <~ 1.0 better than models for Omega_matter=1.0.

A detailed analysis of NGST multiband photometry applied to the reference case of the study of high-redshift galaxies has been carried out with simulations based on galaxy SEDs derived from the currently available empirical and model templates and on plausible standard filter-sets. In order to correctly identify star forming galaxies in the redshift range 5<z<9 and early-type galaxies above z>2 and avoid confusion with other SEDs it is mandatory to have photometric information in optical bands, besides a standard IR filter-set like F110 F160 K L M. In particular by adding the V606, I814 and z-Gunn filters a good discrimination is obtained above z>5 for star forming galaxies and z>1 for early-types. The case for an extension of the NGST wavelength domain to the optical range is therefore strongly supported by this analysis. The effects of leaks in the filter blocking have also been investigated. In spite of rather pessimistic assumptions (a constant leak at a level of 10^{-4} of the peak transmission over the whole spectral range or a leak of Gaussian shape placed at 1.5 times the effective wavelength of the filter with an amplitude of 10^{-3} of the peak transmission and a width of 5% of the central wavelength) the effects are not dramatic: the accuracy of the determination of photometric redshifts with a standard V606 I814 z-Gunn F110 F160 K L M filter set does not significantly deteriorate for a sample limited to M_AB~30. In the range 5.5<z<12.5 a rather high accuracy with a typical error |z_sim-z_phot|<=0.2 is achieved. With the filter blocking achieved by the present standard technology, extended and continuous spectral coverage appears to be the driving factor to maximize the scientific output.

Using high resolution (R ~ 45000), high S/N (20-50) VLT/UVES observations covering the Lyman-alpha forest at 1.5 < z < 2.4 combined with similar observations at z > 2.4 from the literature, we have presented the redshift evolution of the Lyman-alpha forest at 1.5 < z < 4. We applied two parallel analyses: the tradiational Voigt profile fitting and a statistical measure of the flux. The general conclusions from both analyses are in good agreement. In general, the properties of the Lyman-alpha forest at 1.5 < z < 2.4, such as the different density distribution function, the correlation strengths and the mean neutral hydrogen opacity, evolve as expected extrapolating from at z > 2.4. These results suggest the main drive of the redshift evolution of the forest at 1.5 < z < 4 is the Hubble expansion. When combined with the HST observations from the QSO key project at z < 1.7, a flattening of the line number density evolution does seem to occur at z ~ 1.2. The result suggests that the ionizing UV background might not evolve as strongly as a QSO-dominated background, indicating a contribution to the UV background from galaxies at z < 2. The lower cut-off Doppler paramter at the fixed column denisty log N = 13.5 shows a slight increase as z decreases, with a possible local maximum at z ~ 2.9. The slope of the lower cut-off Doppler parameter is not well-defined with z. However, it is flatter at z < 3.1 than at z ~ 3.75.

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.

The galaxy Luminosity Function (LF) has been estimated in the rest frame B luminosity at 0<z<1.25 and at 1700 {\AA} for 2.5<z<4.5 from deep multicolor surveys in the HDF-N, HDF-S, NTT-DF. The results have been compared with a recent version of galaxy formation models in the framework of hierarchical clustering in a flat Cold Dark Matter Universe with cosmological constant. The results show a general agreement for z<= 1, although the model LF has a steeper average slope at the faint end; at z~3 such feature results in an overprediction of the number of faint (I_{AB}~ 27) galaxies, while the agreement at the bright end becomes critically sensitive to the details of dust absorption at such redshifts. The discrepancies at the faint end show that a refined treatement of the physical processes involving smaller galaxies is to be pursued in the models, in terms of aggregation processes and/or stellar feedback heavily affecting the luminosity of the low luminosity objects. The implications of our results on the evolution of the cosmological star formation rate are discussed.

We present progress results of a new survey for bright QSOs (V<14.5, R<15.4, B_J<15.2) covering the whole sky at high galactic latitudes, |b|>30. The surface density of QSOs brighter than B_J=14.8 turns out to be 2.9 +- 0.8 10^{-3} deg^{-2}. The optical Luminosity Function at 0.04 < z < 0.3 shows significant departures from the standard pure luminosity evolution, providing new insights in the modelling of the QSO phenomenon.

We present and compare in this paper new photometric redshift catalogs of the galaxies in three public fields: the NTT Deep Field, the HDF-N and the HDF-S. Photometric redshifts have been obtained for thewhole sample, by adopting a $\chi^2$ minimization technique on a spectral library drawn from the Bruzual and Charlot synthesis models, with the addition of dust and intergalactic absorption. The accuracy, determined from 125 galaxies with known spectroscopic redshifts, is $\sigma_z\sim 0.08 (0.3)$ in the redshift intervals $z=0-1.5 (1.5-3.5)$. The global redshift distribution of I-selected galaxies shows a distinct peak at intermediate redshifts, z~0.6 at I_{AB}<26 and z~0.8 at I_{AB}<27.5 followed by a tail extending to z~6. We also present for the first time the redshift distribution of the total IR-selected sample to faint limits ($Ks \leq 21$ and $J\leq22$). It is found that the number density of galaxies at 1.25<z<1.5 is ~ 0.1 /arcmin^22 at J<=21 and ~1./arcmin^2} at J<22, and drops to 0.3/arcmin^2 (at J<22) at 1.5<z<2. The HDFs data sets are used to compare the different results from color selection criteria and photometric redshifts in detecting galaxies in the redshift range 3.5<z<4.5 Photometric redshifts predict a number of high z candidates in both the HDF-N and HDF-S that is nearly 2 times larger than color selection criteria, and it is shown that this is primarily due to the inclusion of dusty models that were discarded in the original color selection criteria by Madau et al 1998. In several cases, the selection of these objects is made possible by the constraints from the IR bands. Finally, it is shown that galactic M stars may mimic z>5 candidates in the HDF filter set and that the 4 brightest candidates at $z>5$ in the HDF-S are indeed most likely M stars. (ABRIDGED)

We present progress results of a new survey for bright QSOs (V<14.5, R<15.4, B_J<15.2) covering the whole sky at |b|>30. The surface density of QSOs brighter than B_J=14.8 turns out to be $2.9 \pm 0.8 \cdot 10^{-3} deg^{-2}$. The optical Luminosity Function at 0.04 < z < 0.3 shows significant departures from the standard pure luminosity evolution, providing new insights in the modelling of the QSO phenomenon.

We report on high-resolution observations ($\Re \simeq 45000$) of the Hubble Deep Field South QSO J2233-606 obtained with the VLT UV-Visual Echelle Spectrograph (UVES). We present spectral data for the wavelength region $3050 < \lambda < 10000$ \AA. The $S/N$ ratio of the final spectrum is about 50 per resolution element at 4000 \AA, 90 at 5000 \AA, 80 at 6000 \AA, 40 at 8000 \AA. Redshifts, column densities and Doppler widths of the absorption features have been determined with Voigt-profile fitting. A total of 621 lines have been measured. In particular 270 Ly-alpha lines, 41 Ly-beta and 24 systems containing metal lines have been identified. Together with other data in the literature, the present spectrum confirms that the evolution of the number density of Ly-alpha lines with $\log N($\huno$) > 14$ has an upturn at $z \sim 1.4-1.6$.

We report on low-resolution multi-object spectroscopy of 30 faint targets (R \~ 24-25) in the HDF-S and AXAF deep field obtained with the VLT Focal Reducer/low dispersion Spectrograph (FORS1). Eight high-redshift galaxies with 2.75< z < 4 have been identified. The spectroscopic redshifts are in good agreement with the photometric ones with a dispersion $\sigma_z = 0.07$ at z<2 and $\sigma_z = 0.16$ at z>2. The inferred star formation rates of the individual objects are moderate, ranging from a few to a few tens solar masses per year. Five out of the eight high-z objects do not show prominent emission lines. One object has a spectrum typical of an AGN. In the AXAF field two relatively close pairs of galaxies have been identified, with separations of 8.7 and 3.1 proper Mpc and mean redshifts of 3.11 and 3.93, respectively.

This paper presents the first results of a survey for bright quasars (V < 14.5 and R<15.4) covering the North Hemisphere at galactic latitudes |b|>30. The photometric database is derived from the GSC and USNO catalogs. Quasars are identified on the basis of their X-ray emission measured in the ROSAT All Sky Survey. The surface density of quasars brighter than 15.5 mag turns out to be $9 \pm 1 \cdot 10^{-3} deg^{-2}$, about 3 times higher than that estimated by the PG survey. The quasar optical Luminosity Function (LF) at $0.04 < z \le 0.3$ is computed and shown to be consistent with a Luminosity Dependent Luminosity Evolution of the type derived by La Franca and Cristiani (1997) in the range $0.3 < z \le 2.2$. The predictions of semi-analytical models of hierarchical structure formation agree remarkably well with the present observations.

A combination of ground-based (NTT and VLT) and HST (HDF-N and HDF-S) public imaging surveys have been used to collect a sample of 1712 I-selected and 319 $K\leq 21$ galaxies. Photometric redshifts have been obtained for all these galaxies. The results have been compared with the prediction of an analytic rendition of the current CDM hierarchical models for galaxy formation. We focus in particular on two observed quantities: the galaxy redshift distribution at K<21 and the evolution of the UV luminosity density. The derived photometric redshift distribution is in agreement with the hierarchical CDM prediction, with a fraction of only 5% of galaxies detected at z>2. This result strongly supports hierarchical scenarios where present-day massive galaxies are the result of merging processes. The observed UV luminosity density in the I-selected sample is confined within a factor of 4 over the whole range 0<z<4.5. CDM models in a critical Universe are not able to produce the density of UV photons that is observed at z>3. CDM models in $\Lambda$-dominated universe are in better agreement at 3<z<4.5, but predict a pronounced peak at z~1.5 and a drop by a factor of 8 from z=1.5 to z=4 that is not observed in the data. We conclude that improvements are required in the treatment of the physical processes directly related to the SFR, e.g. the starbust activity in merger processes and/or different feedback to the star formation activity.

Spectra, position, magnitudes and colors are presented for 485 faint (B<20.5) emission line objects selected with the ultraviolet-excess (UVX) criterion on a area of 24.6 sq. deg in the South Galactic Pole. The objects were selected from the analysis of pixel-to-pixel stacking of COSMOS scans of UKST U, J and R plates. The candidates were observed with the Meudon-ESO Fiber Optics System (MEFOS) at the ESO 3.6m telescope. 429 type 1 AGNs have been identified (373 in the redshift range 0.3<z<2.2). This sample has allowed the measure of a difference on the QSO clustering evolution in comparison with that found for galaxies (La Franca et al 1998). The region is part of the ESO Imaging Survey (EIS) and of the 2dF QSO redshift survey.

The results of the FORS and ISAAC Science Verification of the FORS and ISAAC instruments at the VLT ANTU/UT1 are described. The following observations have been carried out: 1) the Cluster Deep Field MS1008.1-1224 2) the Antlia dwarf spheroidal Galaxy 3) multiple Object Spectroscopy of Lyman break galaxies in the AXAF and Hubble Deep Field South 4) ISAAC IR Spectroscopy of a gravitationally magnified galaxy at z=2.72. The data have been made public for the ESO community and, in the case of HDF-S, worldwide.

(abridged) The evolution of galaxy clustering from z=0 to z=4.5 is analyzed using the angular correlation function and the photometric redshift distribution of galaxies brighter than I_{AB}\le 28.5 in the HDF North. The reliability of the photometric redshift estimates is discussed on the basis of the available spectroscopic redshifts, comparing different codes and investigating the effects of photometric errors. The redshift bins in which the clustering properties are measured are then optimized to take into account the uncertainties of the photometric redshifts. The results show that the comoving correlation length has a small decrease in the range 0<z<1 followed by an increase at higher z. We compare these results with the theoretical predictions of a variety of cosmological models belonging to the general class of CDM. The comparison with the expected mass clustering evolution indicates that the observed high-redshift galaxies are biased tracers of the dark matter with an effective bias b strongly increasing with redshift. Assuming an Einstein-de Sitter universe, we obtain b\simeq 2 at z=2 and b\simeq 5 at z=4. A comparison of the clustering amplitudes that we measured at z=3 with those reported for LBG suggests that the clustering depends on the abundance of the objects: more abundant objects are less clustered, as expected in the paradigm of hierarchical galaxy formation. The strong clustering and high bias measured at z=3 are consistent with the expected density of massive haloes predicted for the various cosmologies here considered. At z=4, the strong clustering observed in the HDF requires a significant fraction of massive haloes to be already formed by that epoch. This feature could be a discriminant test for the cosmological parameters if confirmed by future observations.

Together with the CMB, the three sources of information that astronomers have at high redshift as probes of the formation and evolution of the LSS are QSOs, galaxies and absorbers observed in the spectrum of distant background objects. In this contribution I try to give a hint of historical perspective, following how the technological advances have driven the emphasis from one class to another, in order to show what are the likely forthcoming milestones.

We present spectroscopic observations of galaxies lying within 1' from the QSO J2233-606 in the HDF-S. Several are found to be coincident in redshift with absorption line systems seen in the HST spectrum of the QSO. We detect a new quasar with z=1.3360 at only 45" angular separation from this QSO. This quasar pair is thus an ideal target for QSO environment studies: indeed, strong HI absorptions are seen at z(abs) = 1.3367 in the QSO spectrum, consistent with neutral hydrogen column densities log N(HI) > 16. A bright galaxy at 4"7 from the QSO is at z = 0.570, corresponding also to HI absorption in the QSO's spectrum. Other galaxies are also discussed.

QSOs and radio-galaxies, together with the CMB, ``normal'' galaxies and clusters, represent the main source of information about the origin and evolution of the Large Scale Structure. They can be used either directly, as tracers of the density peaks, or as cosmic lighthouses, probing the Universe along the line of sight with their conspicuous flow of photons. The amplitude of the QSO clustering has been measured to be $r_o \simeq 6 \hmpc$ at $z\sim 2$ and appears to increase with redshift. These properties are consistent with a scenario in which the quasars, with a lifetime of $t \sim 10^8 yr$, sparsely sample halos of mass greater than $M_{min} \sim 10^{12}-10^{13} h^{-1}$ M$_{\odot}$. Radio-galaxies are found to cluster on similar scales. The study of QSO and radio-galaxy environments is in qualitative agreement with their clustering properties. Absorption-line systems observed in high-resolution spectra of QSOs show structure, including voids/underdensities and overdensities, from the smallest scales to tens of Mpc. QSO pairs, groups and lenses allow us to get a 3-D glimpse of the absorbers, providing invaluable information about their size, ionization, mass. Great progress in this field is expected from the new instrumentation coming into operation at the 8-m class telescopes.

We present a deep BVrI multicolor catalog of galaxies in a 5.62 sq.arcmin field 80 arcsec south of the high redshift (z=4.7) quasar BR 1202-0725, derived from observations with the direct CCD camera SUSI at the ESO NTT. The formal 5$\sigma$ magnitude limits (in 2 x FWHM apertures) are 26.9, 26.5, 25.9 and 25.3 in B, V, r and I respectively. Counts, colors for the star and galaxy samples are discussed and a comparison with a deep HST image in the I band is presented. The percentage of merged or blended galaxies in the SUSI data to this magnitude limit is estimated to be not higher than 1%. At the same galactic latitude of the HDF but pointing toward the galactic center, the star density in this field is found to be ~3 times higher, with ~20% of the objects with V-I > 3.0. Reliable colors have been measured for galaxies selected down to r = 26. The choice of the optical filters has been optimized to define a robust multicolor selection of galaxies at 3.8 <= z <= 4.4. Within this interval the surface density of galaxy candidates with r < 26 in this field is 2.7 +/- 0.4 arcmin^-2 corresponding to a comoving density of Star Formation Rate at 3.8 <= z <= 4.4 of 10^{-2.00} - 10^{-1.82} h^3 Mo /yr/ Mpc^3.

The quasar Q1208+1011 (z_{em}=3.8) is the second highest redshift double quasar ever detected. Several indications point toward it being a gravitational lensed system, although a definitive proof is still lacking. We present new evidence of its lensed nature based on the weakness of the ``proximity effect'' measured in the high resolution Lyman absorption spectrum of the QSO. A luminosity amplification as large as 22 has been derived from this analysis. Indications on the redshift of the lensing galaxy can be obtained from the analysis of the intervening heavy element absorption systems discovered in the QSO high resolution spectrum. On statistical and dynamical grounds a MgII system present at z=1.13 appears as the most likely candidate for the lensing galaxy. We compare the observed parameters with a simple isothermal model for the lens to derive the properties of the lensing galaxy. The resulting magnification factor is smaller, although marginally consistent with that derived by the analysis of the proximity effect.

Spectra of the QSO pair Q0307-195A,B have been obtained in the Ly-alpha Forest (3660 - 3930 A) and C IV (4720 - 4850 A) regions with a FWHM resolution between 0.7 and 0.5 A. 46 lines have been detected in the spectrum of object A while 36 in the spectrum of object B, of them 29 and 20 were identified as Ly-alpha absorptions respectively. The present observations have been supplemented with data of comparable quality on other 7 QSO pairs available in the literature to give an enlarged sample of 217 Ly-alpha lines with rest equivalent width W_o > 0.3 A. The analysis of the hits (i.e. when an absorption line appears in both QSO spectra) and misses (i.e. when a line is seen in any of the QSO spectra, but no line is seen in the other), carried out with an improved statistical approach, indicates that the absorbers have typically a large size (h = H_0/100 km/s/Mpc, q_0 = 0.5): R = 362 h^{-1} kpc, with 95 % confidence limits 298< R <426 h^{-1} kpc and R = 412 h^{-1} kpc, with 95 % confidence limits 333< R <514 h^{-1} kpc for the radius in idealized spherical and disc geometries, respectively. The present data do not allow to establish any correlation of the typical inferred size with the proper separation or with the redshift of the pairs. The correlation between the observed equivalent widths of the absorbers in the adjacent lines of sight becomes poorer and poorer with increasing proper separation. A disc geometry with a column density profile N(r) \propto (r/R_0)^{-gamma}, gamma=4, is found to reasonably reproduce the data with R_0 ~ 100-200 h^{-1} kpc, but also spherical clouds with the same column density profile and a power-law distribution of radii may give a satisfactory representation of the observations.

We present a deep BVrIK multicolor catalog of galaxies in the field of the high redshift (z=4.7) quasar BR 1202-0725. Reliable colors have been measured for galaxies selected down to R=25. The choice of the optical filters has been optimized to define a robust multicolor selection of galaxies at 3.8<z< 4.5. Within this interval the surface density of galaxy candidates with z~4 in this field is 1 arcmin^{-2}. Photometric redshifts have been derived for the galaxies in the field with the maximum likelihood analysis using the GISSEL library of ~10^6 synthetic spectra. The accuracy of the method used has been discussed and tested using galaxies in the Hubble Deep Field (HDF) with known spectroscopic redshifts and accurate photometry. A peak in the redshift distribution is present at z~0.6 with relatively few galaxies at z>1.5. At variance with brighter surveys (I<22.5) there is a tail in the distribution towards high redshifts up to z~4. The luminosity function at z~0.6 shows a steepening for M_B>-19. The observed cosmological ultraviolet luminosity density is computed in the overall redshift interval z=0.3-4.5 reaching a value 2x10^{19} W/Hz/Mpc^3 at z~0.8. We have derived in a homogeneous way, using the GISSEL libraries, the physical parameters connected with the fitted spectral energy distributions. The bulk of the blue intermediate redshift population with z=0.4-1 mostly consists of very young star-forming galaxies with a median starburst age of the order of a few 10^8 yr and typical mass in luminous stars ~2x10^8 Mo.

Recent results have questioned the description of the QSO luminosity function in terms of a pure luminosity evolution and call for a luminosity dependent luminosity evolution. Measurements of the QSO clustering amplitude and evolution allow further distinguishing among the various physical scenarios proposed to interpret the QSO phenomenon. The general properties of the QSO population would arise naturally if quasars are short-lived events connected to a characteristic halo mass $\sim 5 \cdot 10^{12}$ M$_\odot$. This is the typical mass of groups of galaxies in which the interactions triggering the QSO activity preferentially take place.

An ESO Key programme dedicated to an Homogeneous Bright QSO Survey (HBQS) has been completed. 327 QSOs (Mb<-23, 0.3<z<2.2) have been selected over 555 deg^2 with 15<B<18.75. For B<16.4 the QSO surface density turns out to be a factor 2.2 higher than what measured by the PG survey, corresponding to a surface density of 0.013+/-.006 deg^{-2}. If the Edinburgh QSO Survey is included, an overdensity of a factor 2.7 is observed, corresponding to a density of 0.016+/-0.005 deg^{-2}. In order to derive the QSO optical luminosity function (LF) we used Monte Carlo simulations that take into account of the selection criteria, photometric errors and QSO spectral slope distribution. The LF can be represented with a Pure Luminosity Evolution (L(z)\propto(1+z)^k) of a two power law both for q_0=0.5 and q_0=0.1. For q_0=0.5 k=3.26, slower than the previous Boyle's (1992) estimations of k=3.45. A flatter slope beta=-3.72 of the bright part of the LF is also required. The observed overdensity of bright QSOs is concentrated at z<0.6. It results that in the range 0.3<z<0.6 the luminosity function is flatter than observed at higher redshifts. In this redshift range, for Mb<-25, 32 QSOs are observed instead of 19 expected from our best-fit PLE model. This feature requires a luminosity dependent luminosity evolution in order to satisfactorily represent the data in the whole 0.3<z<2.2 interval.

The evolution of QSO clustering is investigated with a new sample of 388 QSOs with 0.3<z<2.2, B<20.5 and M_B<-23. Evidence is found for an increase of the clustering amplitude with increasing redshift. These measurements allow to further distinguish among the various physical scenarios proposed to interpret the QSO phenomenon. A single population model is inconsistent with the observations. The general properties of the QSO population would arise naturally if quasars are short-lived events connected to a characteristic halo mass ~5x10^{12} M_sun.

We present here the spectrum of the galaxy companion to the z=4.7 quasar BR1202-0725, in the optical range 6000-9000 A, corresponding to 1050-1580 A rest-frame. We detect a strong Lya emission line at z=4.702, with an integrated flux of 2x10^{-16} ergs/cm^2/s, and a UV continuum longward of the Lya emission at a flux level of 3x10^{-19}ergs/cm^2/s/A. We fail to detect any CIV[1550] emission with a 3 sigma upper limit of 3x10^{-17} ergs/cm^2/s $. We show that the ratio between Lya and continuum intensity and the absence of a strong CIV emission imply that the UV continuum radiation is the result of an intense star-formation activity rather than of a reprocessing of the QSO flux. The total estimated SFR of this z=4.7 star-forming region is about 15-54 Msun/yr, depending on the IMF and the metallicity. The present data suggest that the Lya emission has a velocity and spatial structure, with possible velocity differences of about 500 km/s on scales of few kpc. These velocity patterns may be a signature of collapsing or merging phenomena in the QSO and its environment, as expected from current models of galaxy formation at high z.

We present the results of the detailed surface photometry of a sample of early-type galaxies in the Hubble Deep Field. Effective radii, surface brightnesses and total V_606 magnitudes have been obtained, as well as U_300, B_450, I_814, J, H and K colors, which are compared with the predictions of chemical-spectrophotometric models of population synthesis. Spectroscopic redshifts are available for 23 objects. For other 25 photometric redshifts are given. In the <mu_e>-r_e plane the early-type galaxies of the HDF, once the appropriate K+E corrections are applied, turn out to follow the `rest frame' Kormendy relation. This evidence, linked to the dynamical information gathered by Steidel et al.(1996), indicates that these galaxies, even at z~2-3, lie in the Fundamental Plane, in a virial equilibrium condition. At the same redshifts a statistically significant lack of large galaxies [i.e. with Log r_e(kpc) > 0.2] is observed.

The observed clustering of Lyman-$\alpha$ lines is reviewed and compared with the clustering of CIV systems. We argue that a continuity of properties exists between Lyman-$\alpha$ and metal systems and show that the small-scale clustering of the absorbers is consistent with a scenario of gravitationally induced correlations. At large scales statistically significant over and under-densities (including voids) are found on scales of tens of Mpc.

The average optical-UV continuum shape of QSOs has been investigated using spectra of 62 QSOs having good relative photometric calibrations. The QSO spectra were extracted from two complete color selected samples in the magnitude intervals B ~ 16-20. The analysis was performed fitting power-law continua (f proportional to nu^(alpha)) in well defined rest-frame wavelength intervals after removing regions of the spectrum affected by strong emission lines or weak emission bumps. The average slope in the rest-frame optical-UV region 1200 - 5500 A shows a rapid change around the 3000 A emission bump with alpha=0.15 longward of it and alpha=-0.65 at shorter wavelengths. Although these average slopes have been obtained using spectra of QSOs with different luminosities and redshifts, there are no significant correlations of the average spectral index with these quantities. For a few QSOs in the sample we were able to measure the same softening of the spectral shape within the individual spectrum. These results have significant consequences on the estimate of the cosmological evolution of the optically selected QSOs as they affect, for instance, the k-corrections. New k-corrections in the B, V, R and Gr bands were computed. The derived average spectral shape in the optical-UV band puts interesting constraints on the expected emission mechanisms.

The evolution of quasar clustering is investigated with a new sample of 388 quasars with 0.3<z<=2.2, B<=20.5 and Mb<-23, selected over an area of 24.6 sq. deg. in the South Galactic Pole. Assuming a two-point correlation function of the form xi(r) = (r/r_o)^-1.8, we detect clustering with r_0=6.2 +/- 1.6 h^-1 comoving Mpc at an average redshift of z=1.3. We find a 2 sigma significant increase of the quasar clustering between z=0.95 and z=1.8, independent of the quasar absolute magnitude and inconsistent with recent evidence on the evolution of galaxy clustering. If other quasar samples are added (resulting in a total data-set of 737 quasars) the increase of the quasar clustering is still favoured although it becomes less significant. We find epsilon=-2.5. Evolutionary parameters epsilon>0.0 are excluded at a 0.3% probability level, to be compared with epsilon=0.8 found for galaxies. The observed clustering properties appear qualitatively consistent with a scenario of Omega=1 CDM in which a) the difference between the quasar and the galaxy clustering can be explained as a difference in the effective bias and redshift distributions, and b) the quasars, with a lifetime of t~10^8 yr, sparsely sample halos of mass greater than M_min~10^12-10^13 h^-1 M_sun. We discuss also the possibility that the observed change in the quasar clustering is due to an increase in the fraction of early-type galaxies as quasar hosts at high z.

The estimate of the ionizing UV background through the proximity effect analysis is discussed. Taking into account the correct bending of the column density distribution which appears in high resolution high s/n data, a value of the UV background ~5x10^{-22} cgs is obtained in the redshift interval z=1.7-4 without indication of any appreciable redshift evolution. Inferences on the spectral shape of the UV background between the HI and HeII edges are provided by the study of metal line ratios like SiIV/CIV in optically thin systems. The presence of a jump near the HeII edge seems favoured at very high redshifts z>3. An increasing HeII jump with increasing redshift can be responsible for the decreasing in the minimum temperatures of the Lyman alpha clouds with increasing z. New studies of the cosmic evolution of star-forming galaxies allow the estimate of a galaxy-dominated UV background intensity ~1.3 x10^{-21} f_{esc} cgs at z=0.5 where f_{esc}<20% is the escape fraction of ionizing UV photons from the galaxies. Finally it is also shown that in a UV background dominated by QSOs and/or star-forming galaxies the cosmological baryon density of the Lyman alpha clouds decreases rapidly with cosmic time possibly due to increasing bulk heating in the intergalactic medium with cosmic time.

We present the results of the detailed surface photometry of a sample of elliptical galaxies in the Hubble Deep Field. In the (surface_brightness-effective_radius) plane, the elliptical galaxies of the HDF turn out to follow a `rest frame' Kormendy relation, once the appropriate K+E corrections are applied. This evidence, linked to the dynamical information gathered by Steidel et al. (1996), indicates that these galaxies, even at z ~ 2-3, lie in the Fundamental Plane, in a virial equilibrium condition. At the same redshifts a statistically signifcant lack of large galaxies [i.e. with log r_e(kpc) > 0.2] is observed.

Spectra of the z_em = 3.268 quasar PKS 2126-158 have been obtained in the range lambda lambda 4300-6620 A with a resolution R = 27000 and an average signal-to-noise ratio s/n = 25 per resolution element. The list of the identified absorption lines is given together with their fitted column densities and Doppler widths. The modal value of the Doppler parameter distribution for the Ly-alpha lines is = 25 km/s. The column density distribution can be described by a power-law dn/dN \propto N^{-beta} with beta = 1.5. 12 metal systems have been identified, two of which were previously unknown. In order to make the column densities of the intervening systems compatible with realistic assumptions about the cloud sizes and the silicon to carbon overabundance, it is necessary to assume a jump beyond the He II edge in the spectrum of the UV ionizing background at z = 3 a factor 10 larger than the standard predictions for the integrated quasar contribution. An enlarged sample of C IV absorptions (71 doublets) has been used to analyze the statistical properties of this class of absorbers strictly related to galaxies. The column density distribution is well described by a single power-law, with beta=1.64 and the Doppler parameter distribution shows a modal value b = 14 km/s. The two point correlation function has been computed in the velocity space for the individual components of C IV features. A significant signal is obtained for scales smaller than 200-300 km/s, xi(30< Delta v < 90 km/s) = 32.71 +- 2.89. A trend of decreasing clustering amplitude with decreasing column density is apparent, analogously to what has been observed for Ly-alpha lines.

An ESO Key programme dedicated to an Homogeneous Bright QSO Survey (HBQS) has been completed. 327 QSOs (Mb<-23, 0.3<z<2.2) have been selected over 555 deg^2 with 15<B<18.75. For B<16.4 the QSO surface density turns out to be a factor 2.2 higher than what measured by the PG survey, corresponding to a surface density of 0.013+/-.006 deg^{-2}. If the Edinburgh QSO Survey is included, an overdensity of a factor 2.7 is observed, corresponding to a density of 0.016+/-0.005 deg^{-2}. In order to derive the QSO optical luminosity function (LF) we used Monte Carlo simulations that take into account of the selection criteria, photometric errors and QSO spectral slope distribution. The LF can be represented with a Pure Luminosity Evolution (L(z)\propto(1+z)^k) of a two power law both for q_0=0.5 and q_0=0.1. For q_0=0.5 k=3.26, slower than the previous Boyle's (1992) estimations of k=3.45. A flatter slope beta=-3.72 of the bright part of the LF is also required. The observed overdensity of bright QSOs is concentrated at z<0.6. It results that in the range 0.3<z<0.6 the luminosity function is flatter than observed at higher redshifts. In this redshift range, for Mb<-25, 32 QSOs are observed instead of 19 expected from our best-fit PLE model. This feature requires a luminosity dependent luminosity evolution in order to satisfactorily represent the data in the whole 0.3<z<2.2 interval.

The long-term variability of a sample of 149 optically selected QSOs in the field of the Selected Area 94 has been studied in the R-band. The relations between variability and luminosity and between variability and redshift have been investigated by means of ``robust'' statistical estimators, allowing to disentangle the effects of the measurement errors. The results are compared with the corresponding properties of the variability in the B-band for the same sample. An anti-correlation between the R-band variability and the intrinsic luminosity is found, analogously to what is observed in the B-band. The amplitude of the R-band variability turns out to be smaller (of a factor $1.13 \pm 0.05$) than the B-band variability. The implications in terms of the black-hole, starburst and microlensing models are discussed.

An ESO Key programme dedicated to an Homogeneous Bright QSO Survey (HBQS) has been completed. 327 QSOs (Mb<-23, 0.3<z<2.2) have been selected over 555 deg^2 with 15<B<18.75. For B<16.4 the QSO surface density turns out to be a factor 2.2 higher than what measured by the PG survey, corresponding to a surface density of 0.013+/-.006 deg^{-2}. If the Edinburgh QSO Survey is included, an overdensity of a factor 2.5 is observed, corresponding to a density of 0.016+/-0.005 deg^{-2}. In order to derive the QSO optical luminosity function (LF) we used Monte Carlo simulations that take into account of the selection criteria, photometric errors and QSO spectral slope distribution. The LF can be represented with a Pure Luminosity Evolution (L(z)\propto(1+z)^k) of a two power law both for q_0=0.5 and q_0=0.1. For q_0=0.5 k=3.26, slower than the previous Boyle's (1992) estimations of k=3.45. A flatter slope beta=-3.72 of the bright part of the LF is also required. The observed overdensity of bright QSOs is concentrated at z<0.6. It results that in the range 0.3<z<0.6 the luminosity function is flatter than observed at higher redshifts. In this redshift range, for Mb<-25, 32 QSOs are observed instead of 19 expected from our best-fit PLE model. This feature requires a luminosity dependent luminosity evolution in order to satisfactorily represent the data in the whole 0.3<z<2.2 interval.

We analyze the clustering properties of a high-resolution (~10 km/s) sample of about 1600 Lyman-alpha lines in the spectra of 15 quasars, obtained in the framework of an ESO key-programme with the addition of literature data. The two-point correlation function in the velocity space shows a significant signal on small velocity scales (Delta v < 300 km/s) with amplitude and significance increasing with increasing column density. The correlation scale at z~3 turns out to be 200 - 300 h_{50}^{-1} kpc. A trend of increasing correlation with decreasing redshift is apparent. The existence of over- and under-densities on scales of a few tens of Megaparsec is confirmed with a high confidence level and a number of possible structures are identified. The present observations are found to be consistent with models of gravitationally induced correlations. A continuity scenario between Lyman-alpha and metal systems emerges, with a suggested physical association between the Lyman-$\alpha$ clouds with log N_{HI} > 14 and the halos of protogalactic systems.

Spectra of the \ze~$=4.12$ quasar Q$0000-2619$ have been obtained in the range $\lambda\lambda = 4880-8115$ \AA~with a resolution of 13 \kms~and signal--to--noise ratio of $s/n=15-60$ per resolution element. The mode of the distribution of the Doppler parameters for the \lya~lines is $\simeq 25$ \kms. The fraction of lines with $10<b<20$ \kms~is 17\%. The Doppler values derived from uncontaminated \lyb~lines are smaller than those obtained from the corresponding Ly$\alpha$ lines, indicating the contribution of non saturated, non resolved components in the \lya~profiles. The integrated UV background estimated from the proximity effect is found to be $J \sim 7 \times 10^{-22}$ erg s$^{-1}$ cm$^{-2}$~Hz$^{-1}$ sr$^{-1}$. This value is consistent with previous estimates obtained at a lower $z$, implying no appreciable redshift evolution of the UVB up to $z=4$. 13 metal systems are identified, five of which previously unknown. The analysis of the associated metal systems suggests abundances generally below the solar value with an average [C/H] $\sim -0.5$. This value is about one order of magnitude higher than that found in intervening systems at about the same redshift. The analysis of the intervening metal line systems has revealed in particular the presence of three optically thin systems with $\log N_{HI}\sim 15$ showing associated CIV and SiIV absorptions. In order to make the observed column densities consistent with [Si/C] ratios lower than 10 times the solar value, it is necessary to assume a large jump in the spectrum of the ionizing UV background beyond the HeII edge ($J_{912}/J_{228}\gsim 1000$). This result, if confirmed in other spectra at the same redshift is suggestive of a possible dominance of a stellar ionizing emissivity over the declining quasar one at $z > 3$.

We present deep BVrI multicolor photometry in the field of the quasar BR1202-07 (z_{em}=4.694) aimed at selecting field galaxies at z>4. We compare the observed colors of the galaxies in the field with those predicted by spectral synthesis models including UV absorption by the intergalactic medium and we define a robust multicolor selection of galaxies at z>4. We provide spectroscopic confirmation of the high redshift QSO-companion galaxy (z=4.702) selected by our method. The first estimate of the surface density of galaxies in the redshift interval 4<z<4.5 is obtained for the same field, corresponding to a comoving volume density of ~ 10^{-3}$ Mpc$^{-3}. This provides a lower limit to the average star formation rate of the order of 10^{-2} Mo/yr/Mpc^{-3} at z~ 4.25.

We present results from high resolution (R\simeq 28000) spectra of six high- redshift QSOs taken at the ESO NTT telescope that allow the detailed study of the Ly-alpha population in the redshift interval z=2.8-4.1. The typical Doppler parameters found for the Ly-alpha lines lie in the interval b=20\div30 km/s, corresponding to temperatures T>24000K, with a fraction of the order 15% in the range 10<b<20km/s. These values are still consistent with models of low density, highly ionized clouds. The observed redshift and column density distributions obtained from these spectra and from the observations of 4 a dditional QSOs taken in the literature allow an accurate estimate of the proximity effect from a relatively large Ly-alpha sample (more than 1100 lines with logN_{HI}>=13.3) in the redshift interval z=1.7-4.1. A Maximum Likelihood analysis has been applied to estimate SIMULTANEOUSLY the best fit parameters of the Ly-alpha$ statistics AND of the UV background. After correcting for the blanketing of weak lines, we confirm that the column density distribution is best represented by a double power-law with a break at logN_{HI}\simeq14, with a slope beta_s=1.8 for higher column densities and a flatter slope beta_f=1.4 below the break. A value J_{LL}= 5 \pm1 \times10^{-22} erg/cm^2/s/Hz/sr is derived for the UV background in the redshift interval z=1.7-4.1, consistent with the predicted QSO contribution. No evidence is found for redshift evolution of the UVB in the same redshift interval. The comoving volume density distributions of protogalactic Damped and Lyman Lymit systems and Ly-alpha clouds with log N_{HI}>= 14 and radii R\simeq 200 kpc are found to be similar, suggesting a possible common association with galaxies.

We have obtained with the SUSI CCD camera on the ESO 3.5m NTT deep images in the BVRI bands of the field centered on the QSO BRI 1202-0725 ($z_{\it em}=4.694$). In the final combined frames the stellar images have FWHM of 1,1,0.6 and 0.65 arcsec respectively. The R and I images show clearly a galaxy $2.2''$ from the QSO, corresponding to $13h^{-1}_{50}$ kpc at $z\sim 4.5$. Possible identification with the metal absorption systems seen in the line of sight to the QSO, including the highest redshift damped system known to date at $z=4.383$, are discussed. We conclude that its colours can be reconciled only with the spectrum of a primeval galaxy at z >~ 4.4, making it the most distant galaxy detected so far. From its magnitudes and models of young galaxy evolution we deduce that it is forming stars at a rate $\sim 30 M_{\odot}$ yr$^{-1}$ and has an estimated age of the order of $10^8$ yr or less, implying that the bulk of the stellar population formed at $z < 6$.

The properties of the Lyman-$\alpha$ absorptions observed in the spectra of QSOs are reviewed: the distribution of column densities and Doppler widths, the redshift evolution, the ``inverse effect'', the clustering. By interpreting the statistics of the line parameters insight is gained about the nature of the absorbers, their sizes, temperatures, confining agents. On the basis of the ``inverse effect'' it is possible to estimate the ultraviolet background radiation at high redshift and, linking this information with the limits on the Gunn-Peterson optical depth, the density of the diffuse part of the intergalactic medium is derived. Future prospects are briefly discussed.

The long-term variability of a sample of 486 optically selected QSOs in the fields of the SA94, SA57 and SGP has been studied. The relations of variability with luminosity and redshift have been investigated by means of statistical estimators that are ``robust'' and allow at the same time to eliminate the influence of the measurement errors. The analysis in the QSOs rest frame of both the ensemble structure function ($SF$) and individual variability indices show: 1) a negative correlation between variability and luminosity; 2) a significant positive correlation between variability and $z$; 3) such correlations may be equally well parameterized either with a model in which the timescale of the variability is fixed for all the QSOs ($\tau \sim 2.4$ yr), while the amplitude linearly increases with the absolute magnitude and $z$, or with a model in which the timescale of the variability linearly depends on the absolute magnitude and the amplitude is only a function of $z$. The same analysis carried out in the observer's frame shows: 1) a negative correlation between variability and luminosity; 2) the timescale of variability does not depend significantly either on the absolute magnitude or on $z$; 3) the ensemble structure function is well represented by a parameterization in which, with a fixed timescale of about 5.5 yr, the amplitude linearly increases with the absolute magnitude; 4) although the general behaviour of the $SF$ does not show a systematic variation of the timescale and/or amplitude with $z$, if we examine the average variability index for objects with $-25>M_B>-27$, we find that below $z=1$ QSOs are significantly less variable than at higher redshift. Implications in terms of the black-hole, starburst and microlensing models are briefly discussed.

The spectra of the QSOs PKS2126-158 ($z_{em}=3.27$) and Q0055-269 ($z_{em}=3.66$) have been observed at high resolution ($2.2 - 2.7 \times 10^4$). Significant clustering, with $\xi \simeq 1$ at $\Delta v=100$ km s$^{-1}$, is detected in both cases for lines with $\log N_{HI} \geq 13.8$. A similar result is obtained analysing the data available in the literature for Q0014+81. An anticlustering signal is also detected on scales $\Delta v \sim 600$ km s$^{-1}$. Two voids of size $\sim 2000$ km s$^{-1}$ are found in the spectrum of Q0055-269 with a probability of $2\times 10^{-4}$.

We have investigated the spatial distribution of quasars and its relationship with redshift by using the two-point correlation function, the variance of cell counts and the conditional density as a function of redshift. By comparing our results with those found by other authors we conclude that quasars have a correlation function with amplitude similar to or greater than that of galaxies, but the issue of its evolution still needs to be settled. We discuss how this can be achieved by adding new observations and comparing real data with survey simulations.

This is the first paper in a series aimed at defining a statistically significant sample of QSOs in the range $ 15 < B < 18.75$ and $ 0.3 < z < 2.2$. The selection is carried out using direct plates obtained at the ESO and UK Schmidt Telescopes, scanned with the COSMOS facility and searched for objects with an ultraviolet excess. Follow-up spectroscopy, carried out at ESO La Silla, is used to classify each candidate. In this initial paper, we describe the scientific objectives of the survey; the selection and observing techniques used. We present the first sample of 285 QSOs ($M_B < -23$) in a 153 deg$^2$ area, covered by the six ``deep'' fields, intended to obtain significant statistics down $B \simeq 18.75$ with unprecedented photometric accuracy. From this database, QSO counts are determined in the magnitude range $ 17 < B < 18.75$.

Attempts to measure the optical depth due to a diffuse intergalactic medium (IGM) in the spectra of high redshift quasars (GP test) within the ESO key program on the intergalactic medium are reviewed. It is shown that there is no evidence for any Gunn-Peterson effect up to the highest redshifts observable in quasar spectra. The Ly$\alpha$ line statistics consistent with this limit is not able to explain the strong absorption observed at the HeII forest in one QSO, leaving room for a true HeII GP effect. The HI/HeII GP ratio implies a steep ionizing UV background at $z\sim 5$ and ionizing sources of stellar origin.

The spectrum of the Q0055-269 ($z=3.66$) has been observed at the resolution of 14 km s$^{-1}$ in the wavelength interval 4750--6300 \AA. The statistical distribution of the Doppler parameter for the Ly$_\alpha$ lines is peaked at $b\simeq 23$ km s$^{-1}$. The column density distribution is described by a power-law with a break at $\log N_{HI} \simeq 14.5$. Significant clustering, with $\xi \simeq 1$ at $\Delta v=100$ km s$^{-1}$, is detected for lines with $\log N_{HI} \geq 13.8$. Two voids of size $\sim 2000$ km s$^{-1}$ are found in the spectrum with a probability of $2\times 10^{-4}$.

The issue of the X-ray to optical luminosity relationship ($L_o-L_x$) is addressed for both optically and X-ray selected quasar samples. We have applied a generalized regression algorithm for the case of samples involving censored data, with errors on both the dependent and the independent variable. Contrary to some previous results, we find that such relationship is consistent with being a linear one ($L_x\propto L_o$). We argue that previous reports of non-linear relationships (i.e. $L_x\propto L_o^e$ with $e < 1$) are due to the neglect of the influence of the photometric errors, the precise knowledge of which strongly influences the reliability of the results. Further progresses in the determination of the $L_o-L_x$ relationship can be achieved with ROSAT observations of the new generation of large bright quasar surveys.

The spectrum of the quasar Q0055-269 ($z=3.66$) has been observed at the resolution of 14 km s$^{-1}$ in the wavelength interval 4750--6300 \AA. We give a list of Lyman$-\alpha$ lines and metal-line systems for which column densities and Doppler widths have been derived by a fitting procedure. The statistical distribution of the Doppler parameter for the Lyman$-\alpha$ lines is peaked at $b\simeq 23$ km s$^{-1}$, with 13\% of lines with $15\leq b\leq 20$ km s$^{-1}$. The column density distribution of the Lyman$-\alpha$ lines is described by a power-law with a break or cutoff at $\log N_{HI} \simeq 14.5$. A featureless power-law distribution is rejected with a probability of $99.94$\%. Significant clustering, with $\xi \simeq 1$ at $\Delta v=100$ km s$^{-1}$, is detected only for lines with $\log N_{HI}\magcir 13.8$. In addition, two voids of size $\sim 2000$ km s$^{-1}$ are found in the spectrum with a random probability of $2\times 10^{-4}$. {}From the proximity effect for lines in the interval $\log N_{HI}=13.3-14.2$ it is possible to infer an UV ionizing background $J=5\times 10^{-22}$ ergs s$^{-1}$ cm$^{-2}$ Hz$^{-1}$ sr$^{-1}$, within a factor of two from the integrated quasar contribution.

The present status of the ESO key-programme ``A Homogeneous Bright QSO Survey'' is described together with the first results concerning the QSO counts and the optical luminosity function. The related analysis of the evolution of the radio and X-ray selected QSOs is presented.

The sky subtraction performances of multi-fiber spectrographs are discussed, analysing in detail the case of the OPTOPUS system at the 3.6 meter ESO telescope at La Silla. A standard technique, based on flat-fields obtained with a uniformly illuminated screen on the dome, provides poor results. A new method has been developed, using the [OI] emission line at 5577 \AA\ as a calibrator of the fiber transmittance, taking into account the diffuse light and the influence of each fiber on the adjacent ones, and correcting for the effects of the image distortions on the sky sampling. In this way the accuracy of the sky subtraction improves from $2-8 \%$ to $1.3-1.6 \%$.

A measure of the average depression between Lyman absorption lines in the spectrum of the faint quasar BR1202-0725 ($z_{em}=4.695$) is presented. The relatively high resolution of the spectrum ($\sim 40$ km s$^{-1}$) allows the selection of regions free of strong absorption lines in the Lyman alpha forest. A reliable evaluation of the continuum shape is based on the careful flux calibration and on the large wavelength interval covered (5000--9300 ~\AA). A best estimate of $\tau _{_{GP}}\leq 0.02\pm 0.03$ has been found for the Gunn-Peterson optical depth at the highest absorption redshift observed at this resolution, $z \simeq 4.3$. The derived baryon density of the intergalactic medium is $\Omega _{IGM}\mincir 0.01$ if the observed quasars are the major contributor to the ionizing UV background flux. This limit, when compared with the total baryon density deduced from the nucleosynthesis, could imply that most of the baryons are already in bound systems at $z\sim 5$.

{}From three quasar samples with a total of 1038 objects in the redshift range $1.0 \div 2.2$ we measure the variance $\sigma^2$ of counts in cells of volume $V_u$. By a maximum likelihood analysis applied separately on these samples we obtain estimates of $\sigma^2(\ell)$, with $\ell \equiv V_u^{1/3}$. The analysis from a single catalog for $\ell = ~40~h^{-1}$ Mpc and from a suitable average over the three catalogs for $\ell = ~60,~80$ and $100~h^{-1}$ Mpc, gives $\sigma^2(\ell) = 0.46^{+0.27}_{-0.27}$, $0.18^{+0.14}_{-0.15}$, $0.05^{+0.14}_{-0.05}$ and $0.12^{+0.13}_{-0.12}$, respectively, where the $70\%$ confidence ranges account for both sampling errors and statistical fluctuations in the counts. This allows a comparison of QSO clustering on large scales with analogous data recently obtained both for optical and IRAS galaxies: QSOs seem to be more clustered than these galaxies by a biasing factor $b_{QSO}/b_{gal} \sim 1.4 - 2.3$.

The average flux decrement shortward the Ly$_{\alpha}$ emission, due to the well-known ``forest'' of absorptions, has been measured in the spectra of 8 quasars. Quasi-simultaneous optical and IUE observations of the two low redshift quasars PKS 0637--75 (z=0.654) and MC 1104+16 (z=0.632) have been carried out, obtaining relatively high S/N, spectrophotometrically calibrated data on their energy distribution from the rest frame H$_{\beta}$ to the Lyman continuum. Six more quasars in the redshift range 2.5-3.4 have been observed in the optical domain. For all the quasars the ``intrinsic'' continuum slope and normalization have been estimated longward the Ly$_{\alpha}$ emission and extrapolated towards the Lyman continuum to measure the average depressions, which have been compared with the model statistics of the Ly$_{\alpha}$ clouds. When all the known classes of absorbers are taken into account with plausible values for their equivalent width distribution and evolution, a good agreement is obtained with the observations. The results for the observed continuum decrement at $z \sim 0.65$ are identical to those predicted by the evolution with redshift of the number of Ly$_{\alpha}$ forest systems including the HST data and within $2\sigma$ of the predicted value using the ``standard'' Ly$_{\alpha}$ evolution (as determined only at high z).