To address the problem concerning the early formation of stars in massive galaxies, we present the results of a semi-analytic model of galaxy formation which includes a physical description of starbursts triggered by galaxy interactions. These originate from the destabilization of cold galactic gas occurring in galaxy encounters, which in part feeds the accretion onto black holes powering quasars, and in part drives circumnuclear starsbursts at redshifts $z\approx 2-4$, preferentially in massive objects. This speeds up the formation of stars in massive galaxies at high redshifts without altering it in low mass galactic halos. Thus, at intermediate $z\approx 1.5-2$ we find that a considerable fraction of the stellar content of massive galaxies is already in place, at variance with the predictions of previous hierarchical models. The resulting high-$z$ star formation rate and B-band luminosity functions, and the luminosity and redshift distribution of galaxies in K-band at $z\lesssim 2$ are all in good agreement with the existing observations concerning the bright galaxy population.

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 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 link the evolution of the galaxies in the hierarchical clustering scenario with the changing accretion rates of cold gas onto the central massive black holes that power the quasars. We base on galaxy interactions as main triggers of accretion; the related scaling laws are taken up from Cavaliere & Vittorini (2000), and grafted to a semi-analytic code for galaxy formation. As a result, at high $z$ the protogalaxies grow rapidly by hierarchical merging; meanwhile, much fresh gas is imported and also destabilized, so the holes are fueled at their full Eddington rates. At lower $z$ the galactic dynamical events are mostly encounters in hierarchically growing groups; now the refueling peters out, as the residual gas is exhausted while the destabilizing encounters dwindle. So, with no parameter tuning other than needed for stellar observables, our model uniquely produces at $z>3$ a rise, and at $z\lesssim 2.5 $ a decline of the bright quasar population as steep as observed. In addition, our results closely fit the observed luminosity functions of quasars, their space density at different magnitudes from $z\approx 5$ to $z\approx 0$, and the local $m_{BH}-\sigma$ relation.

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 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 develop a semi-analytic model of hierarchical galaxy formation with an improved treatment of the evolution of galaxies inside dark matter haloes. We take into account not only dynamical friction processes building up the central dominant galaxy, but also binary aggregations of satellite galaxies inside a common halo described using the kinetic Smoluchowski equation. The description of gas cooling, star formation and evolution, and Supernova feedback follows the standard prescriptions widely used in semi-analytic modelling. We find that binary aggregations are effective in depleting the number of small/intermediate mass galaxies over the redshift range 1<z<3, thus flattening the slope of the luminosity function at the faint end. At z~0 the flattening occurs for -20<M_B<-18, but an upturn is obtained at the very faint end for M_B>-16. We compare our predicted luminosity functions with those obtained from deep multicolor surveys in the rest-frame B and UV bands for the redshift ranges 0<z<1 and 2.5<z<3.5, respectively. The comparison shows that the discrepancy of the predictions of other semi-analytic models with the observations is considerably reduced at z>1 and even more at z ~ 3 by the effect of binary aggregations. The predictions from our dynamical model are discussed and compared with the effects of complementary processes which may conspire in affecting the shape of the luminosity function.

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.

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

The intrinsic sizes of the field galaxies with I<26 in the Hubble and ESO-NTT Deep Fields are shown as a function of their redshifts and absolute magnitudes using photometric redshifts derived from the multicolor catalogs and are compared with the CDM predictions. Extending to lower luminosities and to higher z our previous analysis performed on the NTT field alone, we find that the distribution of the galaxy disk sizes at different cosmic epochs is within the range predicted by typical CDM models. However, the observed size distribution of faint (M_B>-19) galaxies is skewed with respect to the CDM predictions and an excess of small-size disks (R_d<2 kpc) is already present at z~ 0.5. The excess persists up to z~3 and involves brighter galaxies . Such an excess may be reduced if luminosity-dependent effects, like starburst activity in interacting galaxies, are included in the physical mechanisms governing the star formation history in CDM models.

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.

The sizes of the field galaxies with I<25 have been measured in the NTT Deep Field. Intrinsic sizes have been obtained after deconvolution of the PSF with a multigaussian method. The reliability of the method has been tested using both simulated data and HST observations of the same field. The distribution of the half light radii is peaked at r_{hl} 0.3 arcsec, in good agreement with that derived from HST images at the same magnitude. An approximate morphological classification has been obtained using the asymmetry and concentration parameters. The intrinsic sizes of the galaxies are shown as a function of their redshifts and absolute magnitudes using photometric redshifts derived from the multicolor catalog. While the brighter galaxies with morphological parameters typical of the normal spirals show a flat distribution in the range r_{d}=1-6 kpc, the fainter population at 0.4<z<0.8 dominates at small sizes. To explore the significance of this behaviour, an analytical rendition of the standard CDM model for the disc size evolution has been computed. The model showing the best fit to the local luminosity function and the Tully-Fisher relation is able to reproduce at intermediate redshifts a size distribution in general agreement with the observations, although it tends to underestimate the number of galaxies fainter than M_B~ -19 with disk sizes r_d~ 1-2 kpc.

We present the deep UBVRI observations of the HDF-S NICMOS field obtained as part of the Science Verification of the VLT Unit 1 telescope. The images have been used to construct object catalogs and to obtain photometric redshifts. The effective field of view is $\simeq 70\times70$ arcsec$^2$, and the formal $5\sigma$ limiting magnitudes (in a 2 FWHM aperture) are 26.3, 27.8, 27.5, 26.9, 25.2 in the $U$, $B$, $V$, $R$ and $I$ bands, respectively. Thanks to the sub-arcsecond image quality, relatively long exposure time, and large collecting area of the VLT, this is the deepest set of multicolor images ever obtained from a ground-based telescope. Galaxy counts have been derived independently in each band, and show no significant departures from previous data from wider areas. A multicolor photometric catalog of all the galaxies selected in the $R$ band has also been obtained and used to derive photometric redshifts for all galaxies with $R \leq 26.5 $, using also the J,H and K magnitudes from the NICMOS deep observations. A significant fraction (\simeq 28%)of the galaxies is placed at $z\geq 2$. Among them, six robust candidates are found at $z \geq 2.5$ .