Aims. 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. Methods. 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 µm 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. Results. 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 Gyr up to z 1.5, and proceeds much faster thereafter, with an exponential timescale of 0.6 Gyr. 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. Conclusions. The observed behaviour of the Galaxy Stellar Mass Function is consistent with a scenario where about 50% of presentday 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.

Aims. We present a high quality multiwavelength (from 0.3 to 8.0 µm) catalog of the large and deep area in the GOODS Southern Field covered by the deep near-IR observations obtained with the ESO VLT. Methods. 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 µm), 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. Results. The final catalog is made up of 14 847 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.

Aims. The goal of this work is to infer the star formation properties and the mass assembly process of high redshift (0.3 ≤ z < 2.5) galaxies from their IR emission using the 24 μm band of MIPS-Spitzer. Methods. We used an updated version of the GOODS-MUSIC catalog, which has multiwavelength coverage from 0.3 to 24 μm and either spectroscopic or accurate photometric redshifts. We describe how the catalog has been extended by the addition of mid-IR fluxes derived from the MIPS 24 μm image. We compared two different estimators of the star formation rate (SFR hereafter). One is the total infrared emission derived from 24 μm, estimated using both synthetic and empirical IR templates. The other one is a multiwavelength fit to the full galaxy SED, which automatically accounts for dust reddening and age-star formation activity degeneracies. For both estimates, we computed the SFR density and the specific SFR. Results. We show that the two SFR indicators are roughly consistent, once the uncertainties involved are taken into account. However, they show a systematic trend, IR-based estimates exceeding the fit-based ones as the star formation rate increases. With this new catalog, we show that: a) at z > 0.3, the star formation rate is correlated well with stellar mass, and this relationship seems to steepen with redshift if one relies on IR-based estimates of the SFR; 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 M yr −1 . During this epoch, our targeted galaxies assemble a substantial part of their final stellar mass; e) the specific SFR (SSFR) shows a clear bimodal distribution. Conclusions. The analysis of the SFR density and the SSFR seems to support the downsizing scenario, according to which high mass galaxies have formed their stars earlier and more rapidly than their low mass counterparts. A comparison with renditions of theoretical simulations of galaxy formation and evolution indicates that these models follow the global increase in the SSFR with redshift and predict the existence of quiescent galaxies even at z > 1.5. However, the average SSFR is systematically underpredicted by all models considered.

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 9735Å and 9858Å respectively: the lines have fluxes of ∼ 1.6 − 1.2 × 10 −17 ergs −1 cm −2 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α 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 [O ii], [O iii] or Hα emission: however from the spectroscopic and the photometric data we conclude that there are no other plausible identifications, except for Lyα 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.

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