We present the first results of a project called SAGAN, dedicated solely to the studies of relatively rare megaparsec-scale radio galaxies in the Universe, called the Giant Radio Galaxies (GRGs). We have identified 162 new GRGs primarily from the NVSS with sizes ranging from ~0.71 Mpc to ~2.82 Mpc in the redshift range of ~0.03 - 0.95. As part of the project SAGAN, we have created a database of all known GRGs-`GRG-catalogue' from literature (including our new sample) of 820 sources. For the first time, we present the multi-wavelength properties of the largest sample of GRGs, providing new insights about their nature. Our results firmly establish that the distributions of radio spectral index and the black hole mass of GRGs do not differ from the corresponding distributions of normal sized radio galaxies (RGs). However, GRGs have lower Eddington ratio than RGs. Using the mid-infrared data, we have classified GRGs in terms of their accretion mode: either high-power radiatively-efficient, high-excitation state or a radiatively-inefficient low-excitation state. This enables us to compare key physical properties of their AGN like the black hole mass, spin, Eddington ratio, jet kinetic power, total radio power, magnetic field and size. We find that GRGs in high excitation state statistically have larger sizes, radio power, jet kinetic power and Eddington ratio than those in low excitation state. Our analysis reveals a strong correlation between black hole's accretion efficiency and jet kinetic power, thus suggesting a disk-jet coupling. Our environmental study reveals that ~10% of all GRGs may reside at the centres of galaxy clusters, in a denser galactic environment while majority seem to reside in sparse environment. The probability of finding the brightest cluster galaxy (BCG) as GRG is quite low and even lower for high mass clusters. Abridged.

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged)

(Abridged) We investigate the Lyman $\alpha$ emitter luminosity function (LAE LF) within the redshift range $2.9 \leq z \leq 6$ from the first instalment of the blind integral field spectroscopic survey MUSE-Wide. This initial part of the survey probes a region of 22.2 arcmin$^2$ in the CANDELS/GOODS-S field. The dataset provided us with 237 LAEs from which we construct the LAE LF in the luminosity range $42.2 \leq \log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}]\leq 43.5$ within a volume of $2.3\times10^5$ Mpc$^3$. For the LF construction we utilise three different non-parametric estimators: The classical $1/V_\mathrm{max}$ method, the $C^{-}$ method, and an improved binned estimator for the differential LF. All three methods deliver consistent results, with the cumulative LAE LF being $\Phi(\log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}] = 43.5) \simeq 3\times 10^{-6}$ Mpc$^{-3}$ and $\Phi(\log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}] = 42.2) \simeq 2 \times 10^{-3}$ Mpc$^{-3}$ towards the bright- and faint-end of our survey, respectively. By employing a non-parametric statistical test, as well as by comparing the full sample to sub-samples in redshift bins, we find no supporting evidence for an evolving LAE LF over the probed redshift and luminosity range. We determine the best-fitting Schechter function parameters $\alpha = -1.84^{+0.42}_{-0.41}$ and $\log L^* [\mathrm{erg\,s}^{-1}] = 42.2^{+0.22}_{-0.16}$ with the corresponding normalisation $\log \phi^* [\mathrm{Mpc}^{-3}] = -2.71$. When correcting for completeness in the LAE LF determinations, we take into account that LAEs exhibit diffuse extended low surface-brightness haloes. We compare the resulting LF to one obtained where we apply a correction assuming compact point-like emission. We find that the standard correction underestimates the LAE LF at the faint end of our survey by a factor of 2.5.

Galaxies are surrounded by large reservoirs of gas, mostly hydrogen, fed by inflows from the intergalactic medium and by outflows due to galactic winds. Absorption-line measurements along the sightlines to bright and rare background quasars indicate that this circumgalactic medium pervades far beyond the extent of starlight in galaxies, but very little is known about the spatial distribution of this gas. A new window into circumgalactic environments was recently opened with the discovery of ubiquitous extended Lyman-alpha emission from hydrogen around high-redshift galaxies, facilitated by the extraordinary sensitivity of the MUSE instrument at the ESO Very Large Telescope. Due to the faintness of this emission, such measurements were previously limited to especially favourable systems or to massive statistical averaging. Here we demonstrate that low surface brightness Lyman-alpha emission surrounding faint galaxies at redshifts between 3 and 6 adds up to a projected sky coverage of nearly 100%. The corresponding rate of incidence (the mean number of Lyman-alpha emitters penetrated by any arbitrary line of sight) is well above unity and similar to the incidence rate of high column density absorbers frequently detected in the spectra of distant quasars. This similarity suggests that most circumgalactic atomic hydrogen at these redshifts has now been detected also in emission.

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6.

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

We present the MUSE Hubble Ultra Deep Survey, a mosaic of nine MUSE fields covering 90\% of the entire HUDF region with a 10-hour deep exposure time, plus a deeper 31-hour exposure in a single 1.15 arcmin2 field. The improved observing strategy and advanced data reduction results in datacubes with sub-arcsecond spatial resolution (0.65 arcsec at 7000 A) and accurate astrometry (0.07 arcsec rms). We compare the broadband photometric properties of the datacubes to HST photometry, finding a good agreement in zeropoint up to mAB=28 but with an increasing scatter for faint objects. We have investigated the noise properties and developed an empirical way to account for the impact of the correlation introduced by the 3D drizzle interpolation. The achieved 3 sigma emission line detection limit for a point source is 1.5 and 3.1 10-19 erg.s-1.cm-2 for the single ultra-deep datacube and the mosaic, respectively. We extracted 6288 sources using an optimal extraction scheme that takes the published HST source locations as prior. In parallel, we performed a blind search of emission line galaxies using an original method based on advanced test statistics and filter matching. The blind search results in 1251 emission line galaxy candidates in the mosaic and 306 in the ultradeep datacube, including 72 sources without HST counterparts (mAB>31). In addition 88 sources missed in the HST catalog but with clear HST counterparts were identified. This data set is the deepest spectroscopic survey ever performed. In just over 100 hours of integration time, it provides nearly an order of magnitude more spectroscopic redshifts compared to the data that has been accumulated on the UDF over the past decade. The depth and high quality of these datacubes enables new and detailed studies of the physical properties of the galaxy population and their environments over a large redshift range.

We present a clustering analysis of a sample of 238 Ly{$\alpha$}-emitters at redshift 3<z<6 from the MUSE-Wide survey. This survey mosaics extragalactic legacy fields with 1h MUSE pointings to detect statistically relevant samples of emission line galaxies. We analysed the first year observations from MUSE-Wide making use of the clustering signal in the line-of-sight direction. This method relies on comparing pair-counts at close redshifts for a fixed transverse distance and thus exploits the full potential of the redshift range covered by our sample. A clear clustering signal with a correlation length of r0 = 2.9(+1.0/-1.1) Mpc (comoving) is detected. Whilst this result is based on only about a quarter of the full survey size, it already shows the immense potential of MUSE for efficiently observing and studying the clustering of Ly{$\alpha$}-emitters.

GalICS 2.0 is a new semianalytic code to model the formation and evolution of galaxies in a cosmological context. N-body simulations based on a Planck cosmology are used to construct halo merger trees, track subhaloes, compute spins and measure concentrations. The accretion of gas onto galaxies and the morphological evolution of galaxies are modelled with prescriptions derived from hydrodynamic simulations. Star formation and stellar feedback are described with phenomenological models (as in other semianalytic codes). GalICS 2.0 computes rotation speeds from the gravitational potential of the dark matter, the disc and the central bulge. As the rotation speed depends not only on the virial velocity but also on the ratio of baryons to dark matter within a galaxy, our calculation predicts a different Tully-Fisher relation from models in which the rotation speed is proportional to the virial velocity. This is why GalICS 2.0 is able to reproduce the galaxy stellar mass function and the Tully-Fisher relation simultaneously. Our results are also in agreement with halo masses from weak lensing and satellite kinematics, gas fractions, the relation between star formation rate (SFR) and stellar mass, the evolution of the cosmic SFR density, bulge-to-disc ratios, disc sizes and the Faber-Jackson relation.

Massive galaxy clusters are the most violent large scale structures undergoing merger events in the Universe. Based upon their morphological properties in X-rays, they are classified as un-relaxed and relaxed clusters and often host (a fraction of them) different types of non-thermal radio emitting components, viz., haloes, mini-haloes, relics and phoenix within their Intra Cluster Medium (ICM). The radio haloes show steep (alpha = -1.2) and ultra steep (alpha < -1.5) spectral properties at low radio frequencies, giving important insights on the merger (pre or post) state of the cluster. Ultra steep spectrum radio halo emissions are rare and expected to be the dominating population to be discovered via LOFAR and SKA in the future. Further, the distribution of matter (morphological information), alignment of hot X-ray emitting gas from the ICM with the total mass (dark + baryonic matter) and the bright cluster galaxy (BCG) is generally used to study the dynamical state of the cluster. We present here a multi wavelength study on 14 massive clusters from the CLASH survey and show the correlation between the state of their merger in X-ray and spectral properties (1.4 GHz - 150 MHz) at radio wavelengths. Using the optical data we also discuss about the gas-mass alignment, in order to understand the interplay between dark and baryonic matter in massive galaxy clusters.

We report radio imaging and monitoring observations in the frequency range 0.235 - 2.7 GHz during the flaring mode of PKS 2155-304, one of the brightest BL Lac objects. The high sensitivity GMRT observations not only reveal extended kpc-scale jet and FRI type lobe morphology in this erstwhile `extended-core' blazar but also delineate the morphological details, thanks to its arcsec scale resolution. The radio light curve during the end phase of the outburst measured in 2008 shows high variability (8.5%) in the jet emission in the GHz range, compared to the lower core variability (3.2%) seen at the lowest frequencies. The excess of flux density with a very steep spectral index in the MHz range supports the presence of extra diffuse emission at low frequencies. The analysis of multi wavelength (radio/ optical/ gamma-ray) light curves at different radio frequencies confirms the variability of the core region and agrees with the scenario of high energy emission in gamma-rays due to inverse Compton emission from a collimated relativistic plasma jet followed by synchrotron emission in radio. Clearly, these results give an interesting insight of the jet emission mechanisms in blazars and highlight the importance of studying such objects with low frequency radio interferometers like LOFAR and the SKA and its precursor instruments.

We present the first estimate of the Ly{\alpha} luminosity function using blind spectroscopy from the Multi Unit Spectroscopic Explorer, MUSE, in the Hubble Deep Field South. Using automatic source-detection software, we assemble a homogeneously-detected sample of 59 Ly{\alpha} emitters covering a flux range of -18.0 < log10 (F) < -16.3 (erg s^-1 cm^-2), corresponding to luminosities of 41.4 < log10 (L) < 42.8 (erg s^-1). As recent studies have shown, Ly{\alpha} fluxes can be underestimated by a factor of two or more via traditional methods, and so we undertake a careful assessment of each object's Ly{\alpha} flux using a curve-of-growth analysis to account for extended emission. We describe our self-consistent method for determining the completeness of the sample, and present an estimate of the global Ly{\alpha} luminosity function between redshifts 2.91 < z < 6.64 using the 1/Vmax estimator. We find the luminosity function is higher than many number densities reported in the literature by a factor of 2 - 3, although our result is consistent at the 1{\sigma} level with most of these studies. Our observed luminosity function is also in good agreement with predictions from semi-analytic models, and shows no evidence for strong evolution between the high- and low-redshift halves of the data. We demonstrate that one's approach to Ly{\alpha} flux estimation does alter the observed luminosity function, and caution that accurate flux assessments will be crucial in measurements of the faint end slope. This is a pilot study for the Ly{\alpha} luminosity function in the MUSE deep-fields, to be built on with data from the Hubble Ultra Deep Field which will increase the size of our sample by almost a factor of 10.

The VLT Multi Unit Spectroscopic Explorer (MUSE) integral-field spectrograph can detect Ly\alpha{} emitters (LAE) in the redshift range $2.8 \lesssim z \lesssim 6.7$ in a homogeneous way. Ongoing MUSE surveys will notably probe faint Ly\alpha{} sources that are usually missed by current narrow-band surveys. We provide quantitative predictions for a typical wedding-cake observing strategy with MUSE based on mock catalogs generated with a semi-analytic model of galaxy formation coupled to numerical Ly\alpha{} radiation transfer models in gas outflows. We expect $\approx$ 1500 bright LAEs ($F_{Ly\alpha}$ $\gtrsim$ $10^{-17}$ erg s$^{-1}$ cm$^{-2}$) in a typical Shallow Field (SF) survey carried over $\approx$ 100 arcmin$^2$, and $\approx$ 2,000 sources as faint as $10^{-18}$ erg s$^{-1}$ cm$^{-2}$ in a Medium-Deep Field (MDF) survey over 10 arcmin$^2$. In a typical Deep Field (DF) survey of 1 arcmin$^2$, we predict that $\approx$ 500 extremely faint LAEs ($F_{Ly\alpha}$ $\gtrsim$ $4 \times 10^{-19}$ erg s$^{-1}$ cm$^{-2}$) will be found. Our results suggest that faint Ly\alpha{} sources contribute significantly to the cosmic Ly\alpha{} luminosity and SFR budget. While the host halos of bright LAEs at z $\approx$ 3 and 6 have descendants with median masses of $2 \times 10^{12}$ and $5 \times 10^{13}$ $M_{\odot}$ respectively, the faintest sources detectable by MUSE at these redshifts are predicted to reside in halos which evolve into typical sub-$L^{*}$ and $L^{*}$ galaxy halos at z = 0. We expect typical DF and MDF surveys to uncover the building blocks of Milky Way-like objects, even probing the bulk of the stellar mass content of LAEs located in their progenitor halos at z $\approx$ 3.

We explore the properties of high-redshift Lyman-alpha emitters (LAE), and their link with the Lyman-Break galaxy population (LBG), using a semi-analytic model of galaxy formation that takes into account resonant scattering of Lya photons in gas outflows. We can reasonably reproduce the abundances of LAEs and LBGs from redshift 3 to 7, as well as most UV LFs of LAEs. The stronger dust attenuation for (resonant) Lya photons compared to UV continuum photons in bright LBGs provides a natural interpretation to the increase of the LAE fraction in LBG samples, X_LAE, towards fainter magnitudes. The redshift evolution of X_LAE seems however very sensitive to UV magnitudes limits and EW cuts. In spite of the apparent good match between the statistical properties predicted by the model and the observations, we find that the tail of the Lya equivalent width distribution (EW > 100 A) cannot be explained by our model, and we need to invoke additional mechanisms. We find that LAEs and LBGs span a very similar dynamical range, but bright LAEs are about 4 times rarer than LBGs in massive halos. Moreover, massive halos mainly contain weak LAEs in our model, which might introduce a bias towards low-mass halos in surveys which select sources with high EW cuts. Overall, our results are consistent with the idea that LAEs and LBGs make a very similar galaxy population. Their apparent differences seem mainly due to EW selections, UV detection limits, and a decreasing Lya-to-UV escape fraction ratio in high SFR galaxies.

Semi-analytical models (SAMs) are currently the best way to understand the formation of galaxies within the cosmic dark-matter structures. While they fairly well reproduce the local stellar mass functions, correlation functions and luminosity functions, they fail to match observations at high redshift (z > 3) in most cases, particularly in the low-mass range. The inconsistency between models and observations indicates that the history of gas accretion in galaxies, within their host dark-matter halo, and the transformation of gas into stars, are not well followed. Hereafter, we briefly present a new version of the GalICS semi-analytical model. We explore the impacts of classical mechanisms, such as supernova feedback or photoionization, on the evolution of the stellar mass assembly. Even with a strong efficiency, these two processes cannot explain the observed stellar mass function and star formation rate distribution and some other relations. We thus introduce an ad-hoc modification of the standard paradigm, based on the presence of a \textit{no-star-forming} gas component, and a concentration of the star-forming gas in galaxy discs. The main idea behind the existence of the no-star-forming gas reservoir is that only a fraction of the total gas mass in a galaxy is available to form stars. The reservoir generates a delay between the accretion of the gas and the star formation process. This new model is in much better agreement with the observations of the stellar mass function in the low-mass range than the previous models, and agrees quite well with a large set of observations, including the redshift evolution of the specific star formation rate. However, it predicts a large fraction of no-star-forming baryonic gas, potentially larger than observed, even if its nature has still to be examined in the context of the missing baryon problem.

We present an extended version of the semi-analytical model, GalICS. Like its predecessor, eGalICS applies a post-treatment of the baryonic physics on pre-computed dark-matter merger trees extracted from an N-body simulation. We review all the mechanisms that affect, at any given time, the formation and evolution of a galaxy in its host dark-matter halo. We mainly focus on the gas cycle from the smooth cosmological accretion to feedback processes. To follow this cycle with high accuracy we introduce some novel prescriptions: i) a smooth baryonic accretion with two phases: a cold mode and a hot mode built on the continuous dark-matter accretion. In parallel to this smooth accretion, we implement the standard photoionisation modelling to reduce the input gas flow on the smallest structures. ii) a complete monitoring of the hot gas phase. We compute the evolution of the core density, the mean temperature and the instantaneous escape fraction of the hot atmosphere by considering that the hot gas is in hydrostatic equilibrium in the dark-matter potential well, and by applying a principle of conservation of energy on the treatment of gas accretion, supernovae and super massive black hole feedback iii) a new treatment for disc instabilities based on the formation, the migration and the disruption of giant clumps. The migration of such clumps in gas-rich galaxies allows to form pseudo-bulges. The different processes in the gas cycle act on different time scales, and we thus build an adaptive time-step scheme to solve the evolution equations. The model presented here is compared in detail to the observations of stellar-mass functions, star formation rates, and luminosity functions, in a companion paper.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

[abridged] Aims. The aim of this paper is to assess the impact of the interstellar medium (ISM) physics on Lyman-alpha (Lya) radiation transfer and to quantify how galaxy orientation with respect to the line of sight alters observational signatures. Methods. We compare the results of Lya radiation transfer calculations through the ISM of a couple of idealized galaxy simulations with different ISM models. Results. First, the small-scale structuration of the ISM plays a determinant role in shaping a galaxys Lya properties.The artificially warm, and hence smooth, ISM of G1 yields an escape fraction of 50 percent at the Lya line center, and produces symmetrical double-peak profiles. On the contrary, in G2, most young stars are embedded in thick star-forming clouds, and the result is a 10 times lower escape fraction. G2 also displays a stronger outflowing velocity field, which favors the escape of red-shifted photons, resulting in an asymmetric Lya line. Second, the Lya properties of G2 strongly depend on the inclination at which it is observed: From edge-on to face-on, the line goes from a double-peak profile with an equivalent width of -5 Angstrom to a 15 times more luminous red-shifted asymmetric line with EW 90 Angstrom. Conclusions. Lya radiation transfer calculations can only lead to realistic properties in simulations where galaxies are resolved into giant molecular clouds, putting these calculations out of reach of current large scale cosmological simulations. Finally, we find inclination effects to be much stronger for Lya photons than for continuum radiation. This could potentially introduce severe biases in the selection function of narrow-band Lya emitter surveys, which could indeed miss a significant fraction of the high-z galaxy population.

We present a new model for high redshift Lyman-Alpha Emitters (LAEs) in the cosmological context which takes into account the resonant scattering of Ly-a photons through expanding gas. The GALICS semi-analytic model provides us with the physical properties of a large sample of high redshift galaxies. We implement a gas outflow model for each galaxy based on simple scaling arguments. The coupling with a library of numerical experiments of Ly-a transfer through expanding or static dusty shells of gas allows us to derive the Ly-a escape fractions and profiles. The predicted distribution of Ly-a photons escape fraction shows that galaxies with a low star formation rate have a f_esc of the order of unity, suggesting that, for those objects, Ly-a may be used to trace the star formation rate assuming a given conversion law. In galaxies forming stars intensely, the escape fraction spans the whole range from 0 to 1. The model is able to get a good match to the UV and Ly-a luminosity function (LF) data at 3 < z < 5. We find that we are in good agreement with both the bright Ly-a data and the faint population observed by Rauch et al. (2008) at z=3. Most of the Ly-a profiles of our LAEs are redshifted by the diffusion in the outflow which suppresses IGM absorption. The bulk of the observed Ly-a equivalent width (EW) distribution is recovered by our model, but we fail to obtain the very large values sometimes detected. Predictions for stellar masses and UV LFs of LAEs show a satisfactory agreement with observational estimates. The UV-brightest galaxies are found to show only low Ly-a EWs in our model, as it is reported by many observations of high redshift LAEs. We interpret this effect as the joint consequence of old stellar populations hosted by UV-bright galaxies, and high HI column densities that we predict for these objects, which quench preferentially resonant Ly-a photons via dust extinction.

Intermediate mass galaxies (logM(Msun)>10) at z~0.6 are the likeliest progenitors of the present-day numerous population of spirals. There is growing evidence that they have evolved rapidly since the last 6 to 8 Gyr ago, and likely have formed a significant fraction of their stellar mass, often showing perturbed morphologies and kinematics. We have gathered a representative sample of 88 such galaxies and have provided robust estimates of their gas phase metallicity. For doing so, we have used moderate spectral resolution spectroscopy at VLT/FORS2 with unprecedented high S/N allowing to remove biases coming from interstellar absorption lines and extinction to establish robust values of R23=([OII]3727 + [OIII]4959,5007)/Hbeta. We definitively confirm that the predominant population of z~0.6 starbursts and luminous IR galaxies (LIRGs) are on average, two times less metal rich than the local galaxies at a given stellar mass. We do find that the metal abundance of the gaseous phase of galaxies is evolving linearly with time, from z=1 to z=0 and after comparing with other studies, from z=3 to z=0. Combining our results with the reported evolution of the Tully Fisher relation, we do find that such an evolution requires that ~30% of the stellar mass of local galaxies have been formed through an external supply of gas, thus excluding the close box model. Distant starbursts & LIRGs have properties (metal abundance, star formation efficiency & morphologies) similar to those of local LIRGs. Their underlying physics is likely dominated by gas infall probably through merging or interactions. Our study further supports the rapid evolution of z~0.4-1 galaxies. Gas exchanges between galaxies is likely the main cause of this evolution.

Using the multi-integral field spectrograph GIRAFFE at VLT, we have derived the K-band Tully-Fisher relation (TFR) at z~0.6 for a representative sample of 65 galaxies with emission lines. We confirm that the scatter in the z~0.6 TFR is caused by galaxies with anomalous kinematics, and find a positive and strong correlation between the complexity of the kinematics and the scatter that they contribute to the TFR. Considering only relaxed-rotating disks, the scatter, and possibly also the slope of the TFR, do not appear to evolve with z. We detect an evolution of the K-band TFR zero point between z~0.6 and z=0, which, if interpreted as an evolution of the K-band luminosity of rotating disks, would imply that a brightening of 0.66+/-0.14 mag occurs between z~0.6 and z=0. Any disagreement with the results of Flores et al. (2006) are attributed to both an improvement of the local TFR and the more detailed accurate measurement of the rotation velocities in the distant sample. Most of the uncertainty can be explained by the relatively coarse spatial-resolution of the kinematical data. Because most rotating disks at z~0.6 are unlikely to experience further merging events, one may assume that their rotational velocity does not evolve dramatically. If true, our result implies that rotating disks observed at z~0.6 are rapidly transforming their gas into stars, to be able to double their stellar masses and be observed on the TFR at z=0. The rotating disks observed are indeed emission-line galaxies that are either starbursts or LIRGs, which implies that they are forming stars at a high rate. Thus, a significant fraction of the rotating disks are forming the bulk of their stars within 6 to 8 Gyr, in good agreement with former studies of the evolution of the M-Z relation.

We present a first combined analysis of the morphological and dynamical properties for the Intermediate MAss Galaxy Evolution Sequence (IMAGES) sample. It is a representative sample of 52 z~0.6 galaxies with Mstell from 1.5 to 15 10^10Msun and possessing 3D resolved kinematics and HST deep imaging in at least two broad band filters. We aim at evaluate robustly the evolution of rotating spirals since z~0.6, as well as to test the different schemes for classifying galaxies morphologically. We used all the information provided by multi-band images, color maps and 2 dimensional light fitting to assign to each object a morphological class. We divided our sample between spiral disks, peculiar objects, compact objects and mergers. Using our morphological classification scheme, 4/5 of identified spirals are rotating disks and more than 4/5 of identified peculiar galaxies show complex kinematics, while automatic classification methods such as Concentration-Asymmetry and GINI-M20 severely overestimate the fraction of relaxed disk galaxies. Using this methodology, we find that the fraction of rotating spirals has increased by a factor ~ 2 during the last 6 Gyrs, a much higher fraction that found previously based on morphologies alone. These rotating spiral disks are forming stars very rapidly, doubling their stellar masses over the last 6 Gyrs, while most of their stars have been formed few Gyrs earlier, which reveals the presence of a large gas supply. Because they are likely the progenitors of local spirals, we can conjecture how their properties are evolving. Their disks show some evidence for an inside-out growth and the gas supply/accretion is not made randomly as the disk need to be stable in order to match the local disk properties.

We address the origin of the `downsizing' of elliptical galaxies, according to which the stars in more massive galaxies formed earlier and over a shorter period than those in less massive galaxies. We show that this could be the natural result of a shutdown of star formation in dark matter haloes above a critical mass of 10^12MSun. This is demonstrated using a semianalytic simulation of galaxy formation within the standard hierarchical scenario of structure formation. The assumed threshold mass is motivated by the prediction of stable shock heating above this mass and the finding that such a shutdown reproduces the observed distribution of galaxies in luminosity and colour.The shutdown at a critical halo mass introduces a characteristic stellar mass for the transition of galaxies into the `red sequence' of the galaxy colour-magnitude diagram. Central galaxies of haloes that are more massive today have reached this mass earlier and can therefore grow further along the red sequence by dry mergers, ending up more massive and containing older stars. Small galaxies formed in haloes below the critical mass can shutdown late, when they fall into haloes above the critical mass and become satellites. While our semianalytic simulation that incorporates an explicit shutdown reproduces downsizing as inferred from the stellar ages of ellipticals, we explain why it is much harder to detect downsizing using the mass functions of different galaxy types.

(abbreviated) We present the first results of the ESO large program, ``IMAGES'' which aims at obtaining robust measurements of the kinematics of distant galaxies using the multi-IFU mode of GIRAFFE on the VLT. 3D spectroscopy is essential to robustly measure the often distorted kinematics of distant galaxies (e.g., Flores et al. 2006). We derive the velocity fields and $\sigma$-maps of 36 galaxies at 0.4<z<0.75 from the kinematics of the [OII] emission line doublet, and generate a robust technique to identify the nature of the velocity fields based on the pixels of the highest signal-to-noise ratios (S/N). We have gathered a unique sample of 63 velocity fields of emission line galaxies (W0([OII]) > 15 A) at z=0.4-0.75, which are a representative subsample of the population of M_stellar>1.5x10^{10} M_sun emission line galaxies in this redshift range, and are largely unaffected by cosmic variance. Taking into account all galaxies -with or without emission lines- in that redshift range, we find that at least 41+/-7% of them have anomalous kinematics, i.e., they are not dynamically relaxed. This includes 26+/-7% of distant galaxies with complex kinematics, i.e., they are not simply pressure or rotationally supported. Our result implies that galaxy kinematics are among the most rapidly evolving properties, because locally, only a few percent of the galaxies in this mass range have complex kinematics.

The dust produced by the first generation of stars will be a foreground to cosmic microwave background. In order to evaluate the effect of this early dust, we calculate the power spectrum of the dust emission anisotropies and compare it with the sensitivity limit of the Planck satellite. The spatial distribution of the dust is estimated through the distribution of dark matter. At small angular scales ($\ell \gtrsim 1000$) the dust signal is found to be noticeable with the Planck detector for certain values of dust lifetime and production rates. The dust signal is also compared to sensitivities of other instruments. The early dust emission anisotropies are finally compared to those of local dust and they are found to be similar in magnitude at mm wavelengths.

(Abbrev.) This paper prepares a series of papers analysing the Intermediate MAss Galaxy Evolution Sequence (IMAGES) up to z=1. Intermediate mass galaxies (MJ <=-20.3) are selected from the Chandra Deep Field South (CDFS) for which we identify a serious lack of spectroscopically determined redshifts..... We have spectroscopically identified 691 objects including 580 gal., 7 QSOs, and 104 stars. This study provides 531 new redshifts in the CDFS. It confirms the presence of several large scale structures in the CDFS. To test the impact of these structures in the GOODS-South field, we ... compare the evolution of rest-frame U, B, V and K galaxy luminosity densities to that derived from the CFRS. The CDFS field shows a significant excess of luminosity densities in the z=0.5-0.75 range, which increases with the wavelength, reaching up to 0.5 dex at 2.1 um. Stellar mass and specific star formation evolutions might be significantly affected by the presence of the peculiar large scale structures at z= 0.668 and at z= 0.735, that contain a significant excess of evolved, massive galaxies when compared to other fields. This leads to a clear warning to results based on the CDFS/GOODS South fields, especially those related to the evolution of red luminosity densities, i.e. stellar mass density and specific star formation rate. Photometric redshift techniques, when applied to that field, are producing quantities which are apparently less affected by cosmic variance (0.25 dex at 2.1 um), however at the cost of the difficulty in disentangling between evolutionary and cosmic variance effects.

We present the Lensed Mock Map Facility (LeMoMaF), a tool designed to perform mock weak lensing measurements on numerically simulated chunks of the universe. Coupling N-body simulations to a semi-analytical model of galaxy formation, LeMoMaF can create realistic lensed images and mock catalogues of galaxies, at wavelengths ranging from the UV to the submm. To demonstrate the power of such a tool we compute predictions of the source-lens clustering effect on the convergence statistics, and quantify the impact of weak lensing on galaxy counts in two different filters. We find that the source-lens clustering effect skews the probability density function of the convergence towards low values, with an intensity which strongly depends on the redshift distribution of galaxies. On the other hand, the degree of enhancement or depletion in galaxy counts due to weak lensing is independent of the source-lens clustering effect. We discuss the impact on the two-points shear statistics to be measured by future missions like SNAP and LSST. The source-lens clustering effect would bias the estimation of sigma_8 from two point statistics by 2% -5%. We conclude that accurate photometric redshifts for individual galaxies are necessary in order to quantify and isolate the source-lens clustering effect.

IRAS observations show the existence of a correlation between the infrared luminosity Lir and dust temperature Td in star-forming galaxies, in which larger Lir leads to higher dust temperature. The Lir-Td relation is commonly seen as reflecting the increase in dust temperature in galaxies with higher star formation rate. Even though the correlation shows a significant amount of dispersion, a unique relation has been commonly used to construct spectral energy distributions of galaxies in distant universe studies, such as source number counting or photometric redshift determination. In this work, we introduce a new parameter, namely the size of the star-forming region Rir and lay out the empirical and modelled relation between the global parameters Lir, Td and Rir of IR-bright non-AGN galaxies. IRAS 60-to-100um color is used as a proxy for the dust temperature and the 1.4GHz radio contiuum emission for the infrared spatial distribution. The analysis has been carried out on two samples. The first one is made of the galaxies from the 60um flux-limited IRAS Revised Bright Galaxy Samples which have a reliable RC size estimate from the VLA follow-ups of the IRAS Bright Galaxy Samples. The second is made of the sources from the 170um ISOPHOT Serendipity Sky Survey which are resolved by the NVSS or FIRST surveys. We show that the dispersion in the Lir-Td diagram can be reduced to a relation between the infrared surface brightness and the dust temperature, a relation that spans 5 orders of magnitude in surface brightness. We explored the physical processes giving rise to the Sir-Td relation, and show that it can be derived from the Schmidt law, which relates the star formation rate to the gas surface density.

The Multi Unit Spectroscopic Explorer (MUSE) is a second-generation VLT panoramic integral-field spectrograph under preliminary design study. MUSE has a field of 1x1 arcmin**2 sampled at 0.2x0.2 arcsec**2 and is assisted by the VLT ground layer adaptive optics ESO facility using four laser guide stars. The simultaneous spectral range is 465-930 nm, at a resolution of R~3000. MUSE couples the discovery potential of a large imaging device to the measuring capabilities of a high-quality spectrograph, while taking advantage of the increased spatial resolution provided by adaptive optics. This makes MUSE a unique and tremendously powerful instrument for discovering and characterizing objects that lie beyond the reach of even the deepest imaging surveys. MUSE has also a high spatial resolution mode with 7.5x7.5 arcsec**2 field of view sampled at 25 milli-arcsec. In this mode MUSE should be able to obtain diffraction limited data-cubes in the 600-930 nm wavelength range. Although the MUSE design has been optimized for the study of galaxy formation and evolution, it has a wide range of possible applications; e.g. monitoring of outer planets atmosphere, environment of young stellar objects, super massive black holes and active nuclei in nearby galaxies or massive spectroscopic surveys of stellar fields in the Milky Way and nearby galaxies.

We compare the galaxy population of an SPH simulation to those predicted by the GalICS semi-analytic model and a stripped down version without supernova and AGN feedback. The SPH simulation and the no-feedback GalICS model make similar predictions for the baryonic mass functions of galaxies and for the dependence of these mass functions on environment and redshift. The two methods also make similar predictions for the galaxy content of dark matter haloes as a function of halo mass and for the gas accretion history of galaxies. Both the SPH and no-feedback GalICS models predict a bimodal galaxy population at z=0. The "red'' sequence of gas poor, old galaxies is populated mainly by satellite systems while, contrary to observations, the central galaxies of massive haloes lie on the "blue'' star-forming sequence as a result of continuing hot gas accretion at late times. Furthermore, both models overpredict the observed baryonic mass function, especially at the high mass end. In the full GalICS model, supernova-driven outflows reduce the masses of low and intermediate mass galaxies by about a factor of two. AGN feedback suppresses gas cooling in large haloes, producing a sharp cut-off in the baryonic mass function and moving the central galaxies of these massive haloes to the red sequence. Our results imply that the observational failings of the SPH simulation and the no-feedback GalICS model are a consequence of missing input physics rather than computational inaccuracies, that truncating gas accretion by satellite galaxies automatically produces a bimodal galaxy distribution with a red sequence, but that explaining the red colours of the most massive galaxies requires a mechanism like AGN feedback that suppresses the accretion onto central galaxies in large haloes.

[Abridged] We use mock catalogues based on the GALICS model (Hatton et al. 03) to explore the nature of galaxy clustering observed in the SDSS. We measure low and high order angular clustering statistic from these mock catalogues, after selecting galaxies the same way as for observations, and compare them directly to estimates from SDSS data. Note that we also present measurements of S3-S5 on the SDSS DR1. We find that our model is in general good agreement with observations in the scale/luminosity range where we can trust the predictions. This range is found to be limited (i) by the size of the dark matter simulation used -- which introduces finite volume effects at large scales -- and by the mass resolution of this simulation -- which introduces incompleteness at apparent magnitudes fainter than $r\sim 20$. We then focus on the small scale clustering properties of galaxies and investigate the behaviour of three different prescriptions for positioning galaxies within haloes of dark matter. We show that galaxies are poor tracers both of DM particles or DM sub-structures, within groups and clusters. Instead, SDSS data tells us that the distribution of galaxies lies somewhat in between these two populations. This confirms the general theoretical expectation from numerical simulations and semi-analytic modelling.

We reproduce the blue and red sequences in the observed joint distribution of colour and magnitude for galaxies at low and high redshifts using hybrid N-body/semi-analytic simulations of galaxy formation. The match of model and data is achieved by mimicking the effects of cold flows versus shock heating coupled to feedback from active galactic nuclei (AGNs), as predicted by Dekel & Birnboim (2006). After a critical epoch z=3, only haloes below a critical shock-heating mass 10^12MSun enjoy gas supply by cold flows and form stars, while cooling and star formation are shut down abruptly above this mass. The shock-heated gas is kept hot because being dilute it is vulnerable to feedback from energetic sources such as AGNs in their self-regulated mode. The shutdown explains in detail the bright-end truncation of the blue sequence at ~L*, the appearance of luminous red-and-dead galaxies on the red sequence starting already at z~2, the colour bimodality, its strong dependence on environment density and its correlations with morphology and other galaxy properties. Before z~2-3, even haloes above the shock-heating mass form stars by cold streams penetrating through the hot gas. This explains the bright star-forming galaxies at z~3-4, the early appearance of massive galaxies on the red sequence, the high cosmological star-formation rate at high redshifts and the subsequent low rate at low redshifts.

The Active Galactic Nuclei In Cosmological Simulations (AGNICS) project includes AGN into the GalICS hybrid (N-body + semi-analytic) model to explore the quasar-galaxy link in a cosmological perspective. The key problems are the quasar fuelling mechanism, the origin of the BH to bulge mass relation, the causal and chronological link between BH growth and galaxy formation, the properties of quasar hosts and the role of AGN feedback in galaxy formation. This first paper has two goals: to describe the general structure and the assumptions of the method, and to apply AGNICS to studying the joint formation of BHs and spheroids in galaxy mergers. We investigate under what conditions this scenario can reproduce the local BH distribution and the evolution of the quasar population. AGNICS contains two star formation modes: a quiescent one in discs and a starburst one in protospheroids, the latter triggered by mergers and disc instabilities. We assume that BH growth is linked to the starburst mode. The simplest version of this scenario, in which the BH accretion rate and the star formation rate in the starburst component are simply related by a constant of proportionality, does not to reproduce the cosmic evolution of the quasar population. A model in which this relation is modulated by a scaling with the square root of the starburst gas density can explain the evolution of the quasar luminosity function in B-band and X-rays (accounting for the obscuration inferred from X-ray studies). The scatter and the tilt that this model introduces in the BH-to-bulge mass relation are within the observational constraints. The quasar contribution grows with bolometric luminosity and for a given bulge mass the most massive BHs are in the bulges with the oldest stars.

We present the Mock Map Facility, a powerful tool to generate mock catalogues or images from semi-analytically post-processed snapshots of cosmological N-body simulations. The paper describes in detail an efficient technique to create such mocks from the GALICS semi-analytic model, providing the reader with an accurate quantification of the artifacts it introduces at every step. We show that replication effects introduce a negative bias on the clustering signal -- typically peaking at less than 10 percent around the correlation length. We also thoroughly discuss how the clustering signal is affected by finite volume effects, and show that it vanishes at scales larger than about a tenth of the simulation box size. For the purpose of analysing our method, we show that number counts and redshift distributions obtained with GALICS and MOMAF compare well to K-band observations and to the 2dFGRS. Given finite volume effects, we also show that the model can reproduce the APM angular correlation function. The MOMAF results discussed here are made publicly available to the astronomical community through a public database. Moreover, a user-friendly Web interface (http://galics.iap.fr) allows any user to recover her/his own favourite galaxy samples through simple SQL queries. The flexibility of this tool should permit a variety of uses ranging from extensive comparisons between real observations and those predicted by hierarchical models of galaxy formation, to the preparation of observing strategies for deep surveys and tests of data processing pipelines.

High-resolution N-body re-simulations of 15 massive (10^{14}-10^{15} Msun) dark matter haloes have been combined with the hybrid galaxy formation model GalICS (Hatton et al. 2003), to study the formation and evolution of galaxies in clusters, within the framework of the hierarchical merging scenario. New features in GalICS include a better description of galaxy positioning within dark matter haloes, a more reliable computation of the temperature of the inter-galactic medium as a function of redshift, and a description of the ram pressure stripping process. We focus on the luminosity functions, morphological fractions and colour distributions of galaxies in clusters and in cluster outskirts, at z=0. No systematic dependency on cluster richness is found either for the galaxy luminosity functions, morphological mixes, or colour distributions. Moving from higher density (cluster cores), to lower density environments (cluster outskirts), we detect a progressive flattening of the luminosity functions, an increase of the fraction of spirals and a decrease of that of ellipticals and S0s, and the progressive emergence of a bluer tail in the distributions of galaxy colours, especially for spirals. As compared to cluster spirals, early-type galaxies show a flatter luminosity function, and more homogeneous and redder colours. An overall good agreement is found between our results and the observations, particularly in terms of the cluster luminosity functions and morphological mixes. However, some discrepancies are also apparent, with too faint magnitudes of the brightest cluster members, especially in the B band, and galaxy colours tendentially too red (or not blue enough) in the model, with respect to the observations. Finally, ram pressure stripping appears to affect very little our results.

We present predictions for the abundance of Ly-$\alpha$ emitters in hierarchical structure formation models. We use the {\tt GALFORM} semi-analytical model to explore the impact on the predicted counts of varying assumptions about the escape fraction of Ly-$\alpha$ photons, the redshift at which the universe reionised and the cosmological density parameter. A model with a fixed escape fraction gives a remarkably good match to the observed counts over a wide redshift interval. We present predictions for the expected counts in a typical observation with the Multi Unit Spectroscopic Explorer instrument proposed for the Very Large Telescope.

We present a comparison between two optical cluster finding methods: a matched filter algorithm using galaxy angular coordinates and magnitudes, and a percolation algorithm using also redshift information. We test the algorithms on two mock catalogues. The first mock catalogue is built by adding clusters to a Poissonian background, while the other is derived from N-body simulations. Choosing the physically most sensible parameters for each method, we carry out a detailed comparison and investigate advantages and limits of each algorithm, showing the possible biases on final results. We show that, combining the two methods, we are able to detect a large part of the structures, thus pointing out the need to search for clusters in different ways in order to build complete and unbiased samples of clusters, to be used for statistical and cosmological studies. In addition, our results show the importance of testing cluster finding algorithms on different kinds of mock catalogues to have a complete assessment of their behaviour.

This paper illustrates how mock observational samples of high-redshift galaxies with sophisticated selection criteria can be extracted from the predictions of GALICS, a hybrid model of hierarchical galaxy formation that couples the outputs of large cosmological simulations and semi-analytic recipes to describe dark matter collapse and the physics of baryons respectively. As an example of this method, we focus on the properties of Lyman Break Galaxies at redshift 3. With the MOMAF software package described in a companion paper, we generate a mock observational sample with selection criteria as similar as possible to those implied in the actual observations of z = 3 LBGs by Steidel et al.(1995). Our model predictions are in good agreement with the observed number density and 2D correlation function. We investigate the optical/IR luminosity budget as well as several other physical properties of LBGs and find them to be in general agreement with observed values. Looking into the future of these LBGs we predict that 75% of them end up as massive ellipticals today, even though only 35% of all our local ellipticals are predicted to have a LBG progenitor. In spite of some shortcomings, this new 'mock observation' method clearly represents a necessary first step toward a more accurate comparison between hierarchical models of galaxy formation and real observational surveys.

This is the first paper of a series that describes the methods and basic results of the GalICS model (for Galaxies In Cosmological Simulations). GalICS is a hybrid model for hierarchical galaxy formation studies, combining the outputs of large cosmological N-body simulations with simple, semi-analytic recipes to describe the fate of the baryons within dark matter halos. The simulations produce a detailed merging tree for the dark matter halos including complete knowledge of the statistical properties arising from the gravitational forces. We intend to predict the overall statistical properties of galaxies, with special emphasis on the panchromatic spectral energy distribution emitted by galaxies in the UV/optical and IR/submm wavelength ranges. In this paper, we outline the physically motivated assumptions and key free parameters that go into the model, comparing and contrasting with other parallel efforts. We specifically illustrate the success of the model in comparison to several datasets, showing how it is able to predict the galaxy disc sizes, colours, luminosity functions from the ultraviolet to far infrared, the Tully--Fisher and Faber--Jackson relations, and the fundamental plane in the local universe. We also identify certain areas where the model fails, or where the assumptions needed to succeed are at odds with observations, and pay special attention to understanding the effects of the finite resolution of the simulations on the predictions made. Other papers in this series will take advantage of different data sets available in the literature to extend the study of the limitations and predictive power of GalICS, with particular emphasis put on high-redshift galaxies.

We describe MIFS, a second generation integral-field spectrograph for the VLT, operating in the visible wavelength range. It combines a 1'x1' field of view with the improved spatial resolution provided by multi-conjugate adaptive optics and covers a large simultaneous spectral range (0.6-1.0 microns). A separate mode exploits the highest spatial resolution provided by adaptive optics. With this unique combination of capabilities, MIFS has a wide domain of application and a large discovery potential. The MIFS low-spatial resolution mode (sampled at 0.2") combined with the initial MCAO capabilities planned for the VLT will provide ultra-deep fields with a limiting magnitude for spectroscopy of R~28. MIFS will improve the present day detection limit of Lyman alpha emitters by a factor of 100, and will detect low-mass star-forming galaxies to z~7. The MIFS high-spatial resolution mode (3"x3" field sampled at 0.01") is optimized for the next step in (MC)AO. It will probe, e.g., the relationship between supermassive central black holes and their host galaxy and the physics of winds from accretion disks in young stellar objects at unprecedented spatial resolution. MIFS will extend Europe's lead in integral-field spectroscopy. It capitalizes on new developments in adaptive optics, and is a key step towards instrumentation for OWL.

Bottom-up hierarchical formation of dark matter haloes is not as monotonic as implicitly assumed in the Press-Schechter formalism: matter can be ejected into tidal tails, shells or low density `atmospheres'. The implications that the possible truncation of star formation in these tidal `debris' may have for observational galaxy statistics are examined here using the ArFus N-body plus semi-analytical galaxy modelling software. Upper and lower bounds on stellar losses implied by a given set of N-body simulation output data can be investigated by choice of the merging/identity criterion of haloes between successive N-body simulation output times. A median merging/identity criterion is defined and used to deduce an upper estimate of possible star formation and stellar population losses. A largest successor merging/identity criterion is defined to deduce an estimate which minimises stellar losses. In the N-body simulations studied, the debris losses are short range in length and temporary; maximum loss is around 16%. The induced losses for star formation and luminosity functions are strongest (losses of 10%-30%) for low luminosity galaxies and at intermediate redshifts (1 < z < 3). This upper bound on likely losses is smaller than present observational uncertainties. Hence, Press-Schechter based galaxy formation models are approximately valid despite ignoring loss of debris, provided that dwarf galaxy statistics are not under study.

The FIRBACK (Far Infrared BACKground) survey is one of the deepest imaging surveys carried out at 170 microns with ISOPHOT onboard ISO, and is aimed at the study of the structure of the Cosmic Far Infrared Background. This paper provides the analysis of resolved sources. After a validated process of data reduction and calibration, we perform intensive simulations to optimize the source extraction, measure the confusion noise (sigma_c = 45 mJy), and give the photometric and astrometric accuracies. 196 galaxies with flux S > 3 sigma_c are detected in the area of 3.89 square degrees. Counts of sources with flux S > 4 sigma_c present a steep slope of 3.3 +/- 0.6 on a differential "logN-logS" plot between 180 and 500 mJy. As a consequence, the confusion level is high and will impact dramatically on future IR deep surveys. This strong evolution, compared with a slope of 2.5 from Euclidian geometry, is in line with models implying a strongly evolving Luminous Infrared Galaxy population. The resolved sources account for less than 10% of the Cosmic Infrared Background at 170 microns, which is expected to be resolved into sources in the 1 to 10 mJy range.

The next generation of submillimetre/millimetre instruments will provide us with a deeper insight into the mechanisms that rule galaxy formation. As the brightest starbursts are thought to be heavily obscured at optical wavelengths, the opening of this new window will complement the present observations, and enable a detailed investigation of the hierarchical merging of galaxies at remote epochs. In this context, we will present preliminary results of simulations currently developed for deep galaxy surveys with FIRST. N-body simulations of galaxy mergers are being developed to produce realistic morphologies and star formation histories. A Schmidt law is used to account for the evolution of the stellar and gas content of each particle. These simulations will be gathered in a library, including the temporal information with a timestep of 10 Myr. We then build consistent spectra accounting for this star formation history at each resolution element of the simulations. These morphologies will be included in the framework of the GalICS hybrid model of hierarchical galaxy formation, which reproduces the main observational constraints. We intend to present synthetic maps with the characteristics of the FIRST PACS/SPIRE instruments and discuss the optimal strategy for deep surveys complemented with the ground-based ALMA project.

We study the evolution of stellar content and the chemical enrichment of the universe averaged over the whole population of galaxies by means of a series of chemo-spectrophotometric models that take into account the metallicity and dust obscuration effects. We investigate various classes of cosmic star formation rates (CSFR) histories consistent with current estimations. We are able to reproduce a variety of observational constraints such as the emissivities at rest-frame 0.44, 25, 60 and 100 $\mu$m of the local universe and also the overall shape of the extragalactic background light from UV/NIR galaxy counts and the cosmic infrared background (CIB) from DIRBE/FIRAS measurements. We find that the CIB at 140 $\mu$m is crucial for discriminating between the CSFR histories. The best-fit model to this constraint seems to favour mid-infrared derived CSFR at low-$z$ and a flat CSFR at higher $z$ consistent with the most recent estimations, although the shape of the CSFR at high redshifts has little impact on the FIR/submm part of extragalactic background. We suggest that the bulk of the CIB energy is produced by a population of moderately obscured normal galaxies lying at $0\le z \le 1$. We then derive the global chemical enrichment and stellar content of the universe and find that this model predicts metallicities in good agreement with the metallicity in DLAs given that some outflow of metal-enriched gas from galaxies is assumed, but it overproduces the current present-day stellar mass density and NIR luminosity density.

(Abridged) This paper predicts self-consistent faint galaxy counts from the UV to the submm wavelength range. The STARDUST spectral energy distributions described in Devriendt et al. (1999) are embedded within the explicit cosmological framework of a simple semi-analytic model of galaxy formation and evolution. We build a class of models which capture the luminosity budget of the universe through faint galaxy counts and redshift distributions in the whole wavelength range spanned by our spectra. In contrast with a rather stable behaviour in the optical and even in the far-IR, the submm counts are dramatically sensitive to variations in the cosmological parameters and changes in the star formation history. Faint submm counts are more easily accommodated within an open universe with a low value of $\Omega_0$, or a flat universe with a non-zero cosmological constant. This study illustrates the implementation of multi-wavelength spectra into a semi-analytic model. In spite of its simplicity, it already provides fair fits of the current data of faint counts, and a physically motivated way of interpolating and extrapolating these data to other wavelengths and fainter flux levels.

FIRBACK is a one of the deepest surveys performed at 170 microns with ISOPHOT onboard ISO, and is aimed at the study of cosmic far infrared background sources. About 300 galaxies are detected in an area of four square degrees, and source counts present a strong slope of 2.2 on an integral "logN-logS" plot, which cannot be due to cosmological evolution if no K-correction is present. The resolved sources account for less than 10% of the Cosmic Infrared Background at 170 microns. In order to understand the nature of the sources contributing to the CIB, and to explain deep source counts at other wavelengths, we have developed a phenomenological model, which constrains in a simple way the luminosity function evolution with redshift, and fits all the existing deep source counts from the mid-infrared to the submillimetre range. Images, materials and papers available on the FIRBACK web: http://wwwfirback.ias.u-psud.fr wwwfirback.ias.u-psud.fr

A Cosmic Far-InfraRed Background (CFIRB) has long been predicted that would traces the intial phases of galaxy formation. It has been first detected by Puget et al.(1996) using COBE data and has been later confirmed by several recent studies (Fixsen et al. 1998, Hauser et al. 1998, Lagache et al. 1999). We will present a new determination of the CFIRB that uses for the first time, in addition to COBE data, two independent gas tracers: the HI survey of Leiden/Dwingeloo (hartmann, 1998) and the WHAM H$_{\alpha}$ survey (Reynolds et al 1998). We will see that the CFIRB above 100 micron is now very well constrained. The next step is to see if we can detect its fluctuations. To search for the CFIRB fluctuations, we have used the FIRBACK observations. FIRBACK is a deep cosmological survey conducted at 170 micron with ISOPHOT (Dole et al., 2000). We show that the emission of unresolved extra-galactic sources clearly dominates, at arcminute scales, the background fluctuations in the lowest galactic emission regions. This is the first detection of the CFIRB fluctuations.

The current and forthcoming observations of large samples of high-redshift galaxies selected according to various photometric and spectroscopic criteria can be interpreted in the context of galaxy formation, by means of models of evolving spectral energy distributions (SEDs). We hereafter present STARDUST which gives synthetic SEDs from the far UV to the submm wavelength range. These SEDs are designed to be implemented into semi-analytic models of galaxy formation.

The merging history of dark matter halos is computed with the Merging Cell Model proposed by Rodrigues & Thomas (1996). While originally discussed in the case of scale-free power spectra, it is developed and tested here in the framework of the cold dark matter cosmology. The halo mass function, the mass distribution of progenitors and child halos, as well as the probability distribution of formation times, have been computed and compared to the available analytic predictions. The halo auto-correlation function has also been obtained (a first for a semi-analytic merging tree), and tested against analytic formulae. An overall good agreement is found between results of the model, and the predictions derived from the Press & Schechter theory and its extensions. More severe discrepancies appear when formulae that better describe N-body simulations are used for comparison. In many instances, the model can be a useful tool for following the hierarchical growth of structures. In particular, it is suitable for addressing the issue of the formation and evolution of galaxy clusters, as well as the population of Lyman-break galaxies at high redshift, and their clustering properties.

We have constructed the merging history of dark matter halos in a SCDM cosmology, by means of a the {``Merging Cell Model''} (Rodrigues & Thomas 1996). It is based on the linear theory of growth of density fluctuations, and the Top-Hat model, and it takes into account mergers between halos through simplified criteria. The halo mass function, the progenitors and children mass distribution, the distribution of formation times, and the halo autocorrelation function (in the Lagrangian space of initial conditions) have been computed. To test the reliability of the model, results have been compared to the analytical predictions of the extended Press & Schechter theory, that in turn well describes results of N-body simulations. A general good agreement has been found, thus providing us with a reliable tool for rapidly addressing the issue of galaxy formation and evolution in hierarchical scenarios. The model is particularly suitable for studying galaxy clusters at low redshift, and the population of Lyman-break galaxies at high z.

(abridged) We present STARDUST, a new self-consistent modelling of the spectral energy distributions (SEDs) of galaxies from far-UV to radio wavelengths. In order to derive the SEDs in this broad spectral range, we first couple spectrophotometric and (closed-box) chemical evolutions to account for metallicity effects on the spectra of synthetic stellar populations. We then use a phenomenological fit for the metal-dependent extinction curve and a simple geometric distribution of the dust to compute the optical depth of galaxies and the corresponding obscuration curve. This enables us to calculate the fraction of stellar light reprocessed in the infrared range. In a final step, we define a dust model with various components and we fix the weights of these components in order to reproduce the IRAS correlation of IR colours with total IR luminosities. This allows us to compute far-IR SEDs that phenomenologically mimic observed trends. We are able to predict the spectral evolution of galaxies in a broad wavelength range, and we can reproduce the observed SEDs of local spirals, starbursts, luminous infrared galaxies (LIRGs) and ultra luminous infrared galaxies (ULIRGs). This modelling is so far kept as simple as possible and depends on a small number of free parameters, namely the initial mass function (IMF), star formation rate (SFR) time scale, gas density, and galaxy age, as well as on more refined assumptions on dust properties and the presence (or absence) of gas inflows/outflows.

In spite of their overall success, semi-analytic models of galaxy formation and evolution predict slopes of luminosity functions which are steeper than the observed ones. This discrepancy has generally been explained by subtle surface brightness effects acting on the observational samples. In this paper, we explicitly implement the computation of surface brightness in a simple semi-analytic model (with standard CDM), and we estimate the effect of observational surface brightness thresholds on the predicted luminosity functions. The crucial free parameter in this computation is the efficiency epsilon of supernova feedback which is responsible for the triggering of galactic winds. With the classical formalism for this process, it is difficult to reproduce simultaneously the Tully-Fisher relation and the flat slope of the observational luminosity function with the same value of epsilon. This suggests that the triggering of galactic winds is a complex phenomenon. [...] The conclusion of this heuristic paper is that a more realistic description of SN feedback is needed, and that surface brightness effects should not be neglected in the modelling of galaxy formation.

The analysis of the submm anisotropies that will be mapped by the forthcoming MAP and PLANCK satellites requires careful foreground subtraction before measuring CMB fluctuations. Among these, the foreground due to IR/submm thermal radiation from dusty sources was poorly known until recent observational breakthroughs began unveiling the properties of these objects. We hereafter briefly review the observational evidence for a strong evolution of IR/submm sources with respect to the local universe explored by IRAS. We present the basic principles of a new modeling effort where consistent spectral energy distributions of galaxies are implemented into the paradigm of hierarchical clustering with the fashionable semi-analytic approach. This model provides us with specific predictions in IR/submm wavebands, that seem to reproduce the current status of the observations and help assessing the capabilities of forthcoming instruments to pursue the exploration of the deep universe at IR/submm wavelengths. Finally, the ability of the PLANCK High Frequency Instrument all-sky survey to produce a catalogue of dusty sources at submm wavelengths is briefly described.

FIRBACK is one of the deepest cosmological surveys performed in the far infrared, using ISOPHOT. We describe this survey, its data reduction and analysis. We present the maps of fields at 175 microns. We point out some first results: source identifications with radio and mid infrared, and source counts at 175 microns. These two results suggest that half of the FIRBACK sources are probably at redshifts greater than 1. We also present briefly the large follow-up program.

Our view on the deep universe has been so far biased towards optically bright galaxies. Now, the measurement of the Cosmic Infrared Background in FIRAS and DIRBE residuals, and the observations of FIR/submm sources by the ISOPHOT and SCUBA instruments begin unveiling the ``optically dark side'' of galaxy formation. Though the origin of dust heating is still unsolved, it appears very likely that a large fraction of the FIR/submm emission is due to heavily-extinguished star formation. Consequently, the level of the CIRB implies that about 2/3 of galaxy/star formation in the universe is hidden by dust shrouds. In this review, we introduce a new modeling of galaxy formation and evolution that provides us with specific predictions in FIR/submm wavebands. These predictions are compared with the current status of the observations. Finally, the capabilities of current and forthcoming instruments for all-sky and deep surveys of FIR/submm sources are briefly described.

FIRBACK is a deep survey conducted with the ISOPHOT instrument aboard the Infrared Space Observatory (ISO) at an effective wavelength of 175 $\mu$m. We present here results we have obtained on the first field, the so-called Marano1 which covers around 0.25 square degree. We find that the source density for objects with a flux above 200 mJy exceeds the counts expected for sources found in the IRAS deep surveys with a similar flux by about an order of magnitude. Such an excess was expected on the basis of the high far infrared background detected with the FIRAS and DIRBE instruments aboard the Cosmic Background Explorer (COBE). These sources are likely to be redshifted infrared galaxies. The steep number counts indicate strong cosmological evolution in this population. The detected sources account for only 10 % of the cosmic IR background. An extrapolation of the counts down to about 10 mJy would be needeed to account for the whole background at this wavelength.

The recent detection of the Cosmic Infrared Background in FIRAS and DIRBE residuals, and the observations of IR/submm sources by the ISOPHOT and SCUBA instruments have shed new light on the optically dark side of galaxy formation. It turns out that our view on the deep universe has been so far biassed towards optically bright galaxies. We now know that a significant fraction of galaxy/star formation in the universe is hidden by dust shrouds. In this paper, we introduce a new modelling of galaxy formation and evolution that provides us with specific predictions in the IR/submm wavebands. These predictions are compared with the current status of the observations. Finally, the future all-sky and deep surveys with the PLANCK and FIRST missions are briefly described.

The FIRBACK project is a collaboration of eight institutes aimed at using long wavelength surveys with ISO to uncover the nature of the Cosmic Far-Infrared Background. We have surveyed nearly 4 sq. deg. in three regions in the northern and southern hemispheres with the ISOPHOT instrument on ISO, using the 175 micron filter. As reported by Puget et al. (1998) and Reach et al. (these proceedings), the first 0.25 sq. deg. region studied in this programme contains many more 175 micron sources than expected. Initial reduction of the remaining survey areas appears to confirm this result. The area surveyed by FIRBACK represents the largest and deepest long wavelength far-IR survey available before the launch of SIRTF.

In this paper, we report on a first estimate of the contribution of galaxies to the diffuse extragalactic background from the far-UV to the submm, based on semi--analytic models of galaxy formation and evolution. We conclude that the global multi--wavelength picture seems to be consistent provided a quite important fraction of star--formation be hidden in dust--enshrouded systems at intermediate and high--redshift. We show that, according to such models, galaxies cannot stand as important contributors to the background hydrogen-ionizing flux at high-redshift unless neutral hydrogen absorption sites are clumpy and uncorrelated with star forming regions.We briefly discuss the robustness of such a result.

We estimate the evolution of the contribution of galaxies to the cosmic background flux at $912 \AA$ by means of a semi-analytic model of galaxy formation and evolution. Such a modelling has been quite successful in reproducing the optical properties of galaxies. We assume hereafter the high-redshift damped Lyman-$\alpha$ (DLA) systems to be the progenitors of present day galaxies, and we design a series of models which are consistent with the evolution of cosmic comoving emissivities in the available near infrared (NIR), optical, ultraviolet (UV), and far infrared (FIR) bands along with the evolution of the neutral hydrogen content and average metallicity of damped Lyman-$\alpha$ systems (DLA). We use these models to compute the galactic contribution to the Lyman-limit emissivity and background flux for $0 \simeq z \le 4$. We take into account the absorption of Lyman-limit photons by HI and dust in the interstellar medium (ISM) of the galaxies. We find that the background Lyman-limit flux due to galaxies might dominate (or be comparable to) the contribution from quasars at almost all redshifts if the absorption by HI in the ISM is neglected. The ISM HI absorption results in a severe diminishing of this flux--by almost three orders of magnitude at high redshifts to between one and two orders at $z \simeq 0$. Though the resulting galaxy flux is completely negligible at high redshifts, it is comparable to the quasar flux at $z \simeq 0$.

We study how Type Ia/II supernova rates at various redshifts can be used to constrain the cosmic star formation rate (CSFR) history. We use a spectrophotometric model of galaxy evolution which provides a self-consistent description of global properties of local galaxies. We have chosen two CSFR histories, with and without dust extinction, both normalized in order to reproduce the observed comoving luminosity densities at different wavelengths. We show that as expected, the SNII rates provide a direct measurement of the instantaneous CSFR, whereas the SNIa rate depends on the CSFR history during the past few Gyrs. Therefore the SNIa rate detected up to $z {\sim}$ 1 is able to provide new insight on the CSFR at higher redshift. The comparison with the data are in good agreement with both local and high-$z$ measurement and provides an additional evidence of the validity of the CSFR up to $z$ = 1, but suggests values at $z >$ 1 higher than those of Madau et al. (1996) and consistent with the presence of extinction. We also show that Type Ia supernova observations at $z{\sim}$ 1 can already put constraints on the still puzzling nature of Type Ia progenitors.

Recent optical observations have led to a significant progress in our understanding of galaxy formation and evolution. However, our view on the deep universe is currently limited to the starlight which directly escapes from high-redshift galaxies, since we so far ignore the fraction of luminosity absorbed by dust and released in the IR/submm wavelength range. A new constraint is set by the possible detection of the Cosmic Infrared Background. We briefly review the observations and use a semi-analytic model of galaxy formation and evolution to predict number counts consistent with the level of the background. It turns out that the predictions fairly accomodate preliminary data at 175 and 850 microns. This suggests that a significant fraction of star/galaxy formation at high z is hidden by dust.

This paper proposes a new semi-analytic modelling of galaxy properties in the IR/submm wavelength range, which is explicitly set in a cosmological framework. This type of approach has had some success in reproducing the optical properties of galaxies. We hereafter propose a simple extension to the IR/submm range. We estimate the IR/submm luminosities of ``luminous UV/IR galaxies'', and we explore how much star formation could be hidden in heavily--extinguished, ``ultraluminous IR galaxies'' by designing a family of evolutionary scenarios which are consistent with the current status of the ``cosmic constraints'', as well as with the IRAS luminosity function and faint counts, but with different high-z IR luminosity densities. However, these scenarios generate a Cosmic Infrared Background whose spectrum falls within the range of the isotropic IR component detected by Puget et al. (1996) and revisited by Guiderdoni et al. (1997). We give predictions for the faint galaxy counts and redshift distributions at IR and submm wavelengths. The submm range is very sensitive to the details of the evolutionary scenarios. As a result, the on-going and forthcoming observations with ISO and SCUBA (and later with SIRTF, SOFIA, FIRST and PLANCK) will put strong constraints on the evolution of galaxies at z=1 and beyond.

In this letter, we estimate how much starlight was absorbed by dust and released in the infrared during the formation of galaxies. New constraints are set on this so-far hidden history by the recent detection of the Cosmic Infrared Background (CIB) built up from the accumulated light of faint sources (and revisited here). By means of a semi-analytical model of hierarchical galaxy formation, we estimate the population of dust-enshrouded starburst galaxies fitting the CIB. It is argued that a deep ISO survey at 175 micron could already detect these sources and help close the gap in our understanding of the history of light in the Universe.

A new semi-analytic modelling of galaxy evolution in the IR/submm is hereafter outlined. This type of approach successfully reproduces the optical properties of galaxies. We illustrate a simple extension to the IR/submm wavelength range by taking the case of the SCDM model. We design a family of evolutionary scenarios with star formation histories reproducing the evolution of the overall gas and luminosity densities in the universe. These scenarios are also consistent with IRAS data and preliminary ISO counts, and we use them to disentangle the ``Cosmic Infared Background'' detected by Puget et al. (1996) into discrete sources. We finally give predictions for the faint galaxy counts and redshift distributions in various submm wavelengths, which seem to be very sensitive to the details of the evolutionary scenarios. As a consequence, it is easy to anticipate that the current and forthcoming instruments, such as SCUBA, FIRST and PLANCK, will strongly constrain the evolution of high-redshift galaxies.