Using the ligthcone from the cosmological hydrodynamical simulation Horizon-AGN, we produced a photometric catalogue over $0<z<4$ with apparent magnitudes in COSMOS, DES, LSST-like, and Euclid-like filters at depths comparable to these surveys. The virtual photometry accounts for the complex star formation history and metal enrichment of Horizon-AGN galaxies, and consistently includes magnitude errors, dust attenuation and absorption by inter-galactic medium. The COSMOS-like photometry is fitted in the same configuration as the COSMOS2015 catalogue. We then quantify random and systematic errors of photometric redshifts, stellar masses, and star-formation rates (SFR). Photometric redshifts and redshift errors capture the same dependencies on magnitude and redshift as found in COSMOS2015, excluding the impact of source extraction. COSMOS-like stellar masses are well recovered with a dispersion typically lower than 0.1 dex. The simple star formation histories and metallicities of the templates induce a systematic underestimation of stellar masses at $z<1.5$ by at most 0.12 dex. SFR estimates exhibit a dust-induced bimodality combined with a larger scatter (typically between 0.2 and 0.6 dex). We also use our mock catalogue to predict photometric redshifts and stellar masses in future imaging surveys. We stress that adding Euclid near-infrared photometry to the LSST-like baseline improves redshift accuracy especially at the faint end and decreases the outlier fraction by a factor $\sim$2. It also considerably improves stellar masses, reducing the scatter up to a factor 3. It would therefore be mutually beneficial for LSST and Euclid to work in synergy.

The role of the cosmic web in shaping galaxy properties is investigated in the GAMA spectroscopic survey in the redshift range $0.03 \leq z \leq 0.25$. The stellar mass, $u - r$ dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. Comparable trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar mass and large-scale density, the traceless component of the tides from the anisotropic large-scale environment also shapes galactic properties. An extension of excursion theory accounting for filamentary tides provides a qualitative explanation in terms of anisotropic assembly bias: at a given mass, the accretion rate varies with the orientation and distance to filaments. It also explains the absence of type/colour gradients in the data on smaller, non-linear scales.

Current constraints on models of galaxy evolution rely on morphometric catalogs extracted from multi-band photometric surveys. However, these catalogs are altered by selection effects that are difficult to model, that correlate in non trivial ways, and that can lead to contradictory predictions if not taken into account carefully. To address this issue, we have developed a new approach combining parametric Bayesian indirect likelihood (pBIL) techniques and empirical modeling with realistic image simulations that reproduce a large fraction of these selection effects. This allows us to perform a direct comparison between observed and simulated images and to infer robust constraints on model parameters. We use a semi-empirical forward model to generate a distribution of mock galaxies from a set of physical parameters. These galaxies are passed through an image simulator reproducing the instrumental characteristics of any survey and are then extracted in the same way as the observed data. The discrepancy between the simulated and observed data is quantified, and minimized with a custom sampling process based on adaptive Monte Carlo Markov Chain methods. Using synthetic data matching most of the properties of a CFHTLS Deep field, we demonstrate the robustness and internal consistency of our approach by inferring the parameters governing the size and luminosity functions and their evolutions for different realistic populations of galaxies. We also compare the results of our approach with those obtained from the classical spectral energy distribution fitting and photometric redshift approach.Our pipeline infers efficiently the luminosity and size distribution and evolution parameters with a very limited number of observables (3 photometric bands). When compared to SED fitting based on the same set of observables, our method yields results that are more accurate and free from systematic biases.

The variations of galaxy stellar masses and colour-types with the distance to projected cosmic filaments are quantified using the precise photometric redshifts of the COSMOS2015 catalogue extracted from COSMOS field (2 deg$^{2}$). Realistic mock catalogues are also extracted from the lightcone of the cosmological hydrodynamical simulation Horizon-AGN. They show that the photometric redshift accuracy of the observed catalogue ($\sigma_z<0.015$ at $M_*>10^{10}{\rm M}_{\odot}$ and $z<0.9$) is sufficient to provide 2D filaments that closely match their projected 3D counterparts. Transverse stellar mass gradients are measured in projected slices of thickness 75 Mpc between $0.5< z <0.9$, showing that the most massive galaxies are statistically closer to their neighbouring filament. At fixed stellar mass, passive galaxies are also found closer to their filament while active star-forming galaxies statistically lie further away. The contributions of nodes and local density are removed from these gradients to highlight the specific role played by the geometry of the filaments. We find that the measured signal does persist after this removal, clearly demonstrating that proximity to a filament is not equivalent to proximity to an over-density. These findings are in agreement with gradients measured both in 2D or 3D in the Horizon-AGN simulation and those observed in the spectroscopic VIPERS survey (which rely on the identification of 3D filaments). They are consistent with a picture in which the influence of the geometry of the large-scale environment drives anisotropic tides which impact the assembly history of galaxies, and hence their observed properties.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3<z<6$, using the unique database of spectroscopic redshifts in COSMOS, we find $\sigma_{\Delta z/(1+z_s)}$ = 0.021 and $\eta=13.2\% $. The deepest regions reach a 90\% completeness limit of 10$^{10}M_\odot$ to $z=4$. Detailed comparisons of the color distributions, number counts, and clustering show excellent agreement with the literature in the same mass ranges. COSMOS2015 represents a unique, publicly available, valuable resource with which to investigate the evolution of galaxies within their environment back to the earliest stages of the history of the Universe. The COSMOS2015 catalog is distributed via anonymous ftp (ftp://ftp.iap.fr/pub/from_users/hjmcc/COSMOS2015/) and through the usual astronomical archive systems (CDS, ESO Phase 3, IRSA).

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.

We propose a new method of pushing $Herschel$ to its faintest detection limits using universal trends in the redshift evolution of the far infrared over 24$\mu$m colours in the well-sampled GOODS-North field. An extension to other fields with less multi-wavelength information is presented. This method is applied here to raise the contribution of individually detected $Herschel$ sources to the cosmic infrared background (CIRB) by a factor 5 close to its peak at 250$\mu$m and more than 3 in the 350$\mu$m and 500$\mu$m bands. We produce realistic mock $Herschel$ images of the deep PACS and SPIRE images of the GOODS-North field from the GOODS-$Herschel$ Key Program and use them to quantify the confusion noise at the position of individual sources, i.e., estimate a "local confusion noise". Two methods are used to identify sources with reliable photometric accuracy extracted using 24$\mu$m prior positions. The clean index (CI), previously defined but validated here with simulations, which measures the presence of bright 24$\mu$m neighbours and the photometric accuracy index (PAI) directly extracted from the mock $Herschel$ images. After correction for completeness, thanks to our mock $Herschel$ images, individually detected sources make up as much as 54% and 60% of the CIRB in the PACS bands down to 1.1 mJy at 100$\mu$m and 2.2 mJy at 160$\mu$m and 55, 33, and 13% of the CIRB in the SPIRE bands down to 2.5, 5, and 9 mJy at 250$\mu$m, 350$\mu$m, and 500$\mu$m, respectively. The latter depths improve the detection limits of $Herschel$ by factors of 5 at 250$\mu$m, and 3 at 350$\mu$m and 500$\mu$m as compared to the standard confusion limit. Interestingly, the dominant contributors to the CIRB in all $Herschel$ bands appear to be distant siblings of the Milky Way ($z$$\sim$0.96 for $\lambda$$<$300$\mu$m) with a stellar mass of $M_{\star}$$\sim$9$\times$10$^{10}$M$_{\odot}$.

The kinematic analysis of dark matter and hydrodynamical simulations suggests that the vorticity in large-scale structure is mostly confined to, and predominantly aligned with their filaments, with an excess of probability of 20 per cent to have the angle between vorticity and filaments direction lower than 60 degrees relative to random orientations. The cross sections of these filaments are typically partitioned into four quadrants with opposite vorticity sign, arising from multiple flows, originating from neighbouring walls. The spins of halos embedded within these filaments are consistently aligned with this vorticity for any halo mass, with a stronger alignment for the most massive structures up to an excess of probability of 165 per cent. On large scales, adiabatic/cooling hydrodynamical simulations display the same vorticity in the gas as in the dark matter. The global geometry of the flow within the cosmic web is therefore qualitatively consistent with a spin acquisition for smaller halos induced by this large-scale coherence, as argued in Codis et al. (2012). In effect, secondary anisotropic infall (originating from the vortex-rich filament within which these lower-mass halos form) dominates the angular momentum budget of these halos. The transition mass from alignment to orthogonality is related to the size of a given multi-flow region with a given polarity. This transition may be reconciled with the standard tidal torque theory if the latter is augmented so as to account for the larger scale anisotropic environment of walls and filaments.

We use Horizon-AGN, a hydrodynamical cosmological simulation, to explore the role of mergers in the evolution of massive (M > 10^10 MSun) galaxies around the epoch of peak cosmic star formation (1<z<4). The fraction of massive galaxies in major mergers (mass ratio R<4:1) is around 3%, a factor of ~2.5 lower than minor mergers (4:1<R <10:1) at these epochs, with no trend with redshift. At z~1, around a third of massive galaxies have undergone a major merger, while all such systems have undergone either a major or minor merger. While almost all major mergers at z>3 are 'blue' (i.e. have significant associated star formation), the proportion of 'red' mergers increases rapidly at z<2, with most merging systems at z~1.5 producing remnants that are red in rest-frame UV-optical colours. The star formation enhancement during major mergers is mild (~20-40%) which, together with the low incidence of such events, implies that this process is not a significant driver of early stellar mass growth. Mergers (R < 10:1) host around a quarter of the total star formation budget in this redshift range, with major mergers hosting around two-thirds of this contribution. Notwithstanding their central importance to the standard LCDM paradigm, mergers are minority players in driving star formation at the epochs where the bulk of today's stellar mass was formed.

Long-duration gamma-ray bursts (LGRBs) are produced by the collapse of very massive stars. Due to the short lifetime of their progenitors, LGRBs pinpoint star-forming galaxies. We present here a multi-band search for the host galaxy of the long dark GRB 050219A within the enhanced Swift/XRT error circle. We used spectroscopic observations acquired with VLT/X-shooter to determine the redshift and star-formation rate of the putative host galaxy. We compared the results with the optical/IR spectral energy distribution obtained with different facilities. Surprisingly, the host galaxy is a old and quiescent early-type galaxy at z = 0.211 characterised by an unprecedentedly low specific star-formation rate. It is the first LGRB host to be also an early-type post-starburst galaxy. This is further evidence that GRBs can explode in all kind of galaxies, with the only requirement being an episode of star-formation.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

We study the mid-infrared (IR) excess emission of early-type galaxies (ETGs) on the red-sequence at $z <$ 1 using a spectroscopic sample of galaxies in the fields of Great Observatories Origins Deep Survey (GOODS). In the mass-limited sample of 1025 galaxies with $M_{star}$ $>$ 10$^{10.5}$ $M_{\odot}$ and $0.4<z<1.05$, we identify 696 $Spitzer$ 24 $\mu$m detected (above the 5$\sigma$) galaxies and find them to have a wide range of NUV-$r$ and $r$-[12 $\mu$m] colors despite their red optical $u-r$ colors. Even in the sample of very massive ETGs on the red sequence with $M_{star}$ $>$ 10$^{11.2}$ $M_{\odot}$, more than 18% show excess emission over the photospheric emission in the mid-IR. The combination with the results of red ETGs in the local universe suggests that the recent star formation is not rare among quiescent, red ETGs at least out to $z \sim 1$ if the mid-IR excess emission results from intermediate-age stars or/and from low-level ongoing star formation. Our color$-$color diagram including near-UV and mid-IR emissions are efficient not only for identifying ETGs with recent star formation, but also for distinguishing quiescent galaxies from dusty star-forming galaxies.

Gamma-ray bursts (GRBs) have been proposed as a tool for studying star formation in the Universe, so it is crucial to investigate whether their host galaxies and immediate environments are in any way special compared with other star-forming galaxies. Here we present spatially resolved maps of dust emission of the host galaxy of the closest known GRB 980425 at z=0.0085 using our new high-resolution observations from Herschel, Atacama Pathfinder Experiment (APEX), Atacama Large Millimeter Array (ALMA) and Australia Telescope Compact Array (ATCA). We modelled the spectral energy distributions of the host and of the star-forming region displaying the Wolf-Rayet signatures in the spectrum (WR region), located 800 pc from the GRB position. The host is characterised by low dust content and a high fraction of UV-visible star formation, similar to other dwarf galaxies. These galaxies are abundant in the local universe, so it is not surprising to find a GRB in one of them, assuming the correspondence between the GRB rate and star formation. The WR region contributes substantially to the host emission at the far-infrared, millimetre, and radio wavelengths and we propose that this is a consequence of its high gas density. If dense environments are also found close to the positions of other GRBs, then the ISM density should also be considered, along with metallicity, an important factor influencing whether a given stellar population can produce a GRB.

Until recently, dust emission has been detected in very few host galaxies of gamma-ray bursts (GRBHs). With Herschel, we have now observed 17 GRBHs up to redshift z~3 and detected seven of them at infrared (IR) wavelengths. This relatively high detection rate (41%) may be due to the composition of our sample which at a median redshift of 1.1 is dominated by the hosts of dark GRBs. Although the numbers are small, statistics suggest that dark GRBs are more likely to be detected in the IR than their optically-bright counterparts. Combining our IR data with optical, near-infrared, and radio data from our own datasets and from the literature, we have constructed spectral energy distributions (SEDs) which span up to 6 orders of magnitude in wavelength. By fitting the SEDs, we have obtained stellar masses, dust masses, star-formation rate (SFR), and extinctions for our sample galaxies. We find that GRBHs are galaxies that tend to have a high specfic SFR (sSFR), and like other star-forming galaxies, their ratios of dust-to-stellar mass are well correlated with sSFR. We incorporate our Herschel sample into a larger compilation of GRBHs, and compare this combined sample to SFR-weighted median stellar masses of the widest, deepest galaxy survey to date. This is done in order to establish whether or not GRBs can be used as an unbiased tracer of cosmic comoving SFR density (SFRD) in the universe. In contrast with previous results, this comparison shows that GRBHs are medium-sized galaxies with relatively high sSFRs; stellar masses and sSFRs of GRBHs as a function of redshift are similar to what is expected for star-forming galaxy populations at similar redshifts. We conclude that there is no strong evidence that GRBs are biased tracers of SFRD; thus they should be able to reliably probe the SFRD to early epochs.

We seek to derive star formation rates (SFR) and stellar masses (M_star) in distant galaxies and to quantify the main uncertainties affecting their measurement. We explore the impact of the assumptions made in their derivation with standard calibrations or through a fitting process, as well as the impact of the available data, focusing on the role of IR emission originating from dust. We build a sample of galaxies with z>1, all observed from the UV to the IR (rest frame). The data are fitted with the code CIGALE, which is also used to build and analyse a catalogue of mock galaxies. Models with different SFHs are introduced. We define different set of data, with or without a good sampling of the UV range, NIR, and thermal IR data. The impact of these different cases on the determination of M_star and SFR are analysed. Exponentially decreasing models with a redshift formation of the stellar population z ~8 cannot fit the data correctly. The other models fit the data correctly at the price of unrealistically young ages when the age of the single stellar population is taken to be a free parameter. The best fits are obtained with two stellar populations. As long as one measurement of the dust emission continuum is available, SFR are robustly estimated whatever the chosen model is, including standard recipes. M_star measurement is more subject to uncertainty, depending on the chosen model and the presence of NIR data, with an impact on the SFR-M_star scatter plot. Conversely, when thermal IR data from dust emission are missing, the uncertainty on SFR measurements largely exceeds that of stellar mass. Among all physical properties investigated here, the stellar ages are found to be the most difficult to constrain and this uncertainty acts as a second parameter in SFR measurements and as the most important parameter for M_star measurements.

We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 {\mu}Jy. Material of the model (predictions, SED library, mock catalogs...) is available online at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

We study the evolution of the galactic spin using data of high redshift galaxies in the fields of the Great Observatories Origins Deep Survey (GOODS). Through simple dynamical considerations we estimate the spin for the disc galaxies in our sample and find that its distribution is consistent with that found for nearby galaxies. Defining a dimensionless angular momentum parameter for the disc component of the galaxies ($\lambda_{d}$), we do not find signs of evolution in the redshift range $0.4 \leq z \leq 1.2$. We find that the mass and environmental dependence of the spin of our high redshift galaxies are similar to that of low-$z$ galaxies; showing a strong dependence on mass, in the sense that low-mass systems present higher $\lambda_{d}$ values than high-mass galaxies, with no significant dependence on the environmental density. These results lead us to conclude that, although individual disc galaxies might occasionally suffer from strong evolution, they evolve in such a way that the overall spin distribution of the galactic population remains constant from $z\sim1$ to the present epoch.

Using a mass-limited sample of 24um-detected, star-forming galaxies at 0.5<z<1.3, we study the mass-star formation rate (SFR) correlation and its tightness. The correlation is well defined (sigma=0.28dex) for disk galaxies (n_sersic<1.5), while more bulge-dominated objects often have lower specific SFRs. For disk galaxies, a much tighter correlation (sigma=0.19dex) is obtained if the rest-frame H-band luminosity is used instead of stellar mass derived from multicolor photometry. The specific SFR (sSFR) correlates strongly with rest-frame optical colors (hence luminosity-weighted stellar age) and also with clumpiness (which likely reflects the molecular gas fraction). This implies that most of the observed scatter is real, despite its low level, and not dominated by random measurement errors. After correcting for these differential effects a remarkably small dispersion remains (sigma=0.14dex), suggesting that measurement errors in mass or SFR are ~0.10dex, excluding systematic uncertainties. Measurement errors in stellar masses, the thickening of the correlation due to real sSFR variations, and varying completeness with stellar mass, can spuriously bias the derived slope to lower values due to the finite range over which observables (mass and SFR) are available. When accounting for these effects, the intrinsic slope for the main sequence for disk galaxies gets closer to unity.

Long duration gamma-ray bursts (GRBs) allow us to explore the distant Universe, and are potentially the most effective tracer of the most distant objects. Our current knowledge of the properties of GRB host galaxies at redshifts >5 is very scarce. We propose to improve this situation by obtaining more observations of high-redshift hosts to better understand their properties and help enable us to use GRBs as probes of the high-redshift universe. We performed very deep photometric observations of three high-redshift GRB host galaxies, GRB 080913 at z =6.7, GRB 060927 at z =5.5 and GRB 060522 at z =5.1. In addition, we completed deep spectroscopic observations of the GRB080913 host galaxy with X-Shooter at the VLT to search for Ly-alpha emission. For the sake of the discussion, we use published results on another high-redshift GRB host, GRB 050904 at z = 6.3. The sample of GRB host galaxies studied in this paper consists of four out of the five spectroscopically confirmed GRBs at z>5. Despite our presented observations being the deepest ever reported of high-redshift GRB host galaxies, we do not detect any of the hosts, neither in photometry nor in spectroscopy in the case of GRB 080913. These observations indicate that the GRB host galaxies seem to evolve with time and to have lower SFRs at z >5 than they have at z<1. In addition, the host galaxy of GRB 080913 at z =6.7 does not show Ly-alpha emission. While the measured properties of the galaxies in our sample agree with the properties of the general galaxy population at z>5, our observations are not sufficiently sensitive to allow us to infer further conclusions on whether this specific population is representative of the general one. The characterization of high-redshift GRB host galaxies is a very challenging endeavor requiring a lot of telescope time, but is necessary to improve our understanding of the high-redshift universe.

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

We study the impact of galaxy-galaxy interactions on the FIR properties of galaxies and its evolution at 0<z<1.2. Using the high-z galaxies in the fields of GOODS observed by Herschel in the framework of the GOODS-Herschel key program and the local IRAS or AKARI-selected galaxies in the field of SDSS DR7, we investigate the dependence of galaxy properties on the morphology of and the distance to the nearest neighbor galaxy. We find that the SFRs and the SSFRs of galaxies, on average, depend on the morphology of and the distance to the nearest neighbor galaxy in this redshift range. When a late-type galaxy has a close neighbor galaxy, the SFR and the SSFR increase as it approaches a late-type neighbor, which is supported by K-S and MC tests with a significance level of >99%. However, the SFR and the SSFR decrease or do not change much as it approaches an early-type neighbor. The bifurcations of SFRs and SSFRs depending on the neighbor's morphology seem to occur at R_n~0.5r_vir,nei, which is supported by K-S and MC tests with a significance level of >98%. For all redshift bins, the SSFRs of late-type galaxies interacting with late-type neighbors are increased by factors of about 1.8+/-0.7 and 4.0+/-1.2 compared to those of non-interacting galaxies when the pair separation is smaller than 0.5r_vir,nei and 0.1r_vir,nei, respectively. The Tdust of both local and high-z late-type galaxies strongly interacting with late-type neighbors (i.e. R_n < 0.1r_vir,nei) appears to be higher than that of non-interacting galaxies with a significance level of 96-99%. However, Tdust of local late-type galaxies strongly interacting with early-type neighbors seems to be lower than or similar to that of non-interacting galaxies. Our results suggest that galaxy-galaxy interactions and mergers have been strongly affecting the SFR and the dust properties of star-forming galaxies over at least 8 billion years.

We present an analysis of the size growth seen in early-type galaxies over 10 Gyr of cosmic time. Our analysis is based on a homogeneous synthesis of published data from 17 spectroscopic surveys observed at similar spatial resolution, augmented by new measurements for galaxies in the Gemini Deep Deep Survey. In total, our sample contains structural data for 465 galaxies (mainly early-type) in the redshift range 0.2<z<2.7. The size evolution of passively-evolving galaxies over this redshift range is gradual and continuous, with no evidence for an end or change to the process around z~1, as has been hinted at by some surveys which analyze subsets of the data in isolation. The size growth appears to be independent of stellar mass, with the mass-normalized half-light radius scaling with redshift as R_e (1+z)^(-1.62 +/- 0.34). Surprisingly, this power law seems to be in good agreement with the recently reported continuous size evolution of UV-bright galaxies in the redshift range z~0.5-3.5. It is also in accordance with the predictions from recent theoretical models.

We aim at modeling the infrared galaxy evolution in an as simple as possible way and reproduce statistical properties among which the number counts between 15 microns and 1.1 mm, the luminosity functions, and the redshift distributions. We then aim at using this model to interpret the recent observations (Spitzer, Akari, BLAST, LABOCA, AzTEC, SPT and Herschel), and make predictions for future experiments like CCAT or SPICA. This model uses an evolution in density and luminosity of the luminosity function with two breaks at redshift ~0.9 and 2 and contains the two populations of the Lagache et al. (2004) model: normal and starburst galaxies. We also take into account the effect of the strong lensing of high-redshift sub-millimeter galaxies. It has 13 free parameters and 8 additional calibration parameters. We fit the parameters to the IRAS, Spitzer, Herschel and AzTEC measurements with a Monte-Carlo Markov chain. The model ajusted on deep counts at key wavelengths reproduces the counts from the mid-infrared to the millimeter wavelengths, as well as the mid-infrared luminosity functions. We discuss the contribution to the cosmic infrared background (CIB) and to the infrared luminosity density of the different populations. We also estimate the effect of the lensing on the number counts, and discuss the recent discovery by the South Pole Telescope (SPT) of a very bright population lying at high-redshift. We predict confusion level for future missions using a P(D) formalism, and the Universe opacity to TeV photons due to the CIB.

We derive the evolution of the infrared (IR) luminosity function (LF) over the last 4/5ths of cosmic time, using deep 24um and 70um imaging of the GOODS North and South fields. We use an extraction technique based on prior source positions at shorter wavelengths to build the 24 and 70um source catalogs. The majority (93%) of the sources have a spectroscopic (39%) or a photometric redshift (54%) and, in our redshift range of interest (i.e., 1.3<z<2.3) ~20% of the sources have a spectroscopic redshifts. To extend our study to lower 70um luminosities we perform a stacking analysis and we characterize the observed L_24/(1+z) vs L_70/(1+z) correlation. Using spectral energy distribution templates which best fit this correlation, we derive the IR luminosity of sources from their 24 and 70 um fluxes. We then compute the IR LF at z=1.55+/-0.25 and z=2.05+/-0.25. The redshift evolution of the IR LF from z=1.3 to z=2.3 is consistent with a luminosity evolution proportional to (1+z)^1.0+/-0.9 combined with a density evolution proportional to (1+z)^-1.1+/-1.5. At z~2, luminous IR galaxies (LIRGs: 10^11Lsun< LIR <10^12Lsun) are still the main contributors to the total comoving IR luminosity density (IR LD) of the Universe. At z~2, LIRGs and ultra-luminous IR galaxies (ULIRGs: 10^12Lsun< LIR) account for ~49% and ~17% respectively of the total IR LD of the Universe. Combined with previous results for galaxies at z<1.3 and assuming a constant conversion between the IR luminosity and star-formation rate (SFR) of a galaxy, we study the evolution of the SFR density of the Universe from z=0 to z=2.3. We find that the SFR density of the Universe strongly increased with redshift from z=0 to z=1.3, but is nearly constant at higher redshift out to z=2.3. As part of the online material accompanying this article, we present source catalogs at 24um and 70um for both the GOODS-North and -South fields.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

The correlation between the large-scale distribution of galaxies and their spectroscopic properties at z=1.5 is investigated using the Horizon MareNostrum cosmological run. We have extracted a large sample of 10^5 galaxies from this large hydrodynamical simulation featuring standard galaxy formation physics. Spectral synthesis is applied to these single stellar populations to generate spectra and colours for all galaxies. We use the skeleton as a tracer of the cosmic web and study how our galaxy catalogue depends on the distance to the skeleton. We show that galaxies closer to the skeleton tend to be redder, but that the effect is mostly due to the proximity of large haloes at the nodes of the skeleton, rather than the filaments themselves. This effects translate into a bimodality in the colour distribution of our sample. The origin of this bimodality is investigated and seems to follow from the ram pressure stripping of satellite galaxies within the more massive clusters of the simulation. The virtual catalogues (spectroscopical properties of the MareNostrum galaxies at various redshifts) are available online at http://www.iap.fr/users/pichon/MareNostrum/catalogues

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

A broad continuum excess in the near-infrared, peaking in the rest-frame at 2-5 micron, is detected in a spectroscopic sample of 88 galaxies at 0.5<z<2.0 taken from the Gemini Deep Deep Survey. Line emission from polycyclic aromatic hydrocarbons (PAHs) at 3.3 micron alone cannot explain the excess, which can be fit by a spectral component consisting of a template of PAH emission lines superposed on a modified blackbody of temperature T~850 K. The luminosity of this near-infrared excess emission at 3 micron is found to be correlated with the star formation rate of the galaxy. The origin of the near-infrared excess is explored by examining similar excesses observed locally in massive star forming regions, reflection and planetary nebulae, post-asymptotic giant branch stars and in the galactic cirrus. We also consider the potential contribution from dust heated around low-luminosity active galactic nuclei. We conclude that the most likely explanation for the 2-5 micron excess is the contribution from circumstellar disks around massive young stellar objects seen in the integrated light of high-redshift galaxies. Assuming circumstellar disks extend down to lower masses, as they do in our own Galaxy, the excess emission presents us with an exciting opportunity to measure the formation rate of planetary systems at cosmic epochs before our own solar system formed.

We discuss observational evidence that quasars play a key role in the formation of galaxies starting from the detailed study of the quasar HE0450-2958 and extending the discussion to a series of converging evidence that radio jets may trigger galaxy formation. The direct detection with VISIR at the ESO-VLT of the 7 kpc distant companion galaxy of HE0450-2958 allows us to spatially separate the sites of quasar and star formation activity in this composite system made of two ultra-luminous infrared galaxies (ULIRGs). No host galaxy has yet been detected for this quasar, but the companion galaxy stellar mass would bring HE0450-2958 in the local M(BH)-M(stellar bulge) relation if it were to merge with the QSO. This is bound to happen because of their close distance (7 kpc) and small relative velocity (~60-200 km/s). We conclude that we may be witnessing the building of the M(BH)-M(stellar bulge) relation, or at least of a major event in that process. The star formation rate (~340 Msun/yr), age (40-200 Myr) and stellar mass ([5-6]x10^10 Msun) are consistent with jet-induced formation of the companion galaxy. We suggest that HE0450-2958 may be fueled in fresh material by cold gas accretion from intergalactic filaments. We map the projected galaxy density surrounding the QSO as a potential tracer of intergalactic filaments and discuss a putative detection. Comparison to other systems suggests that inside-out formation of quasar host galaxies and jet-induced galaxy formation may be a common process. Two tests are proposed for this new paradigm: (1) the detection of offset molecular gas or dust emission with respect to the position of distant QSOs, (2) the delayed formation of host galaxies as a result of QSO activity, hence the two step building of the M(BH)/M(stellar bulge) ratio.

[Abridged] We use the deepest existing mid- and far-infrared observations (reaching ~3 mJy at 70 um) obtained with Spitzer in the GOODS and FIDEL fields to derive the evolution of the rest-frame 15 um, 35 um, and total infrared luminosity functions of galaxies spanning z < 1.3. In comparison with previous studies, the present one takes advantage of deep 70 um observations that provide a more robust infrared luminosity indicator than 24 um affected by the emission of PAHs at high redshift (z~1), and we use several independent fields to control cosmic variance. We use a new extraction technique based on the well-determined positions of galaxies at shorter wavelengths to extract the 24 and 70 um flux densities of galaxies. Using a combination of photometric and spectroscopic redshifts that exist for ~80% of the sources in our sample, we are able to estimate the rest-frame luminosities of galaxies at 15um and 35um. Then using a careful stacking analysis to validate the infrared bolometric corrections we compute the infrared luminosity function in three redshift bins z~0.55, z~0.85, and z~1.15. We find that the average infrared spectral energy distribution of galaxies over the last 2/3 of the cosmic time is consistent with that of local galaxies, although individual sources do present significant scatter. We also measured both the bright and faint ends of the infrared luminosity functions and find no evidence for a change in the slope of the double power law used to characterize the luminosity function. The redshift evolution of infrared luminous galaxies is consistent with pure luminosity evolution proportional to (1+z)^3.6 up to z~1.3. The comoving number density of infrared luminous galaxies has increased by a factor of ~100 between 0< z <1. By z~1.0, LIRGs produce half of the total comoving infrared luminosity density.

[Abridged] This paper aims at providing new conservative constraints to the cosmic star-formation history from the empirical modeling of mid- and far-infrared data. We perform a non-parametric inversion of galaxy counts at 15, 24, 70, 160, and 850 microns simultaneously. It is a "blind" search (no redshift information is required) of all possible evolutions of the infrared luminosity function of galaxies, from which the evolution of the star-formation rate density and its uncertainties are derived. The cosmic infrared background (CIRB) measurements are used a posteriori to tighten the range of solutions. The inversion relies only on two hypotheses: (1) the luminosity function remains smooth both in redshift and luminosity, (2) a set of infrared spectral energy distributions (SEDs) of galaxies must be assumed. The range of star-formation histories that we derive is well constrained and consistent with redshift-based measurements from deep surveys. The redshift decompositions of the counts are also recovered successfully. Therefore, multi-wavelength counts and CIRB (both projected observations) alone seem to contain enough information to recover the cosmic star-formation history with quantifiable errors. A peak of the SFRD at z~2 is preferred, although higher redshifts are not excluded. We also find a good consistency between the observed evolution of the stellar mass density and the prediction from our model. Finally, the inability of the inversion to model perfectly and simultaneously all the multi-wavelength infrared counts (especially at 160 microns where an excess is seen around 20 mJ) implies either (i) the existence of a sub-population of colder galaxies, (ii) a larger dispersion of dust temperatures among local galaxies than expected, (iii) or a redshift evolution of the infrared SEDs of galaxies.

We present the results of NICMOS imaging of a sample of 16 high mass passively evolving galaxies with 1.3<z<2, taken primarily from the Gemini Deep Deep Survey. Around 80% of galaxies in our sample have spectra dominated by stars with ages >1 Gyr. Our rest-frame R-band images show that most of these objects have compact regular morphologies which follow the classical R^1/4 law. These galaxies scatter along a tight sequence in the Kormendy relation. Around one-third of the massive red objects are extraordinarily compact, with effective radii under one kiloparsec. Our NICMOS observations allow the detection of such systems more robustly than is possible with optical (rest-frame UV) data, and while similar systems have been seen at z>2, this is the first time such systems have been detected in a rest-frame optical survey at 1.3<z<2. We refer to these compact galaxies as "red nuggets". Similarly compact massive galaxies are completely absent in the nearby Universe. We introduce a new "stellar mass Kormendy relation" (stellar mass density vs size) which isolates the effects of size evolution from those of luminosity and color evolution. The 1.1 < z < 2 passive galaxies have mass densities that are an order of magnitude larger then early type galaxies today and are comparable to the compact distant red galaxies at 2 < z < 3. We briefly consider mechanisms for size evolution in contemporary models focusing on equal-mass mergers and adiabatic expansion driven by stellar mass loss. Neither of these mechanisms appears able to transform the high-redshift Kormendy relation into its local counterpart. <ABRIDGED>

We present the most extensive and complete study of the properties for the largest sample (46 objects) of gamma-ray burst (GRB) host galaxies. The redshift interval and the mean redshift of the sample are 0<z<6.3 and z=0.96 (look-back time: 7.2 Gyr), respectively; 89% of the hosts are at z <~ 1.6. Optical-near-infrared (NIR) photometry and spectroscopy are used to derive stellar masses, star formation rates (SFRs), dust extinctions and metallicities. The average stellar mass is 10^9.3 M_sun, with a 1 sigma dispersion of 0.8 dex. The average metallicity for a subsample of 17 hosts is about 1/6 solar and the dust extinction in the visual band (for a subsample of 10 hosts) is A_V=0.5. We obtain new relations to derive SFR from [OII] or UV fluxes, when Balmer emission lines are not available. SFRs, corrected for dust extinction, aperture-slit loss and stellar Balmer absorption are in the range 0.01-36 M_sun yr^-1. The median SFR per unit stellar mass (specific SFR) is 0.8 Gyr^-1. Equivalently the inverse quantity, the median formation timescale is 1.3 Gyr. Most GRBs are associated with the death of young massive stars, more common in star-forming galaxies. Therefore GRBs are an effective tool to detect star-forming galaxies in the universe. Star-forming galaxies at z<1.6 are a faint and low-mass population, hard to detect by conventional optical-NIR surveys, unless a GRB event occurs. There is no compelling evidence that GRB hosts are peculiar galaxies. More data on the subclass of short GRB are necessary to establish the nature of their hosts.

At redshifts around 0.1 the CFHT Legacy Survey Deep fields contain some 6x10^4 galaxies spanning the mass range from 10^5 to 10^12 Msun. We measure the stellar mass dependence of the two point correlation using angular measurements to largely bypass the errors, approximately 0.02 in the median, of the photometric redshifts. Inverting the power-law fits with Limber's equation we find that the auto-correlation length increases from a very low 0.4hMpc at 10^5.5 Msun to the conventional 4.5hMpc at 10^10.5 Msun. The power law fit to the correlation function has a slope which increases from gamma approximately 1.6 at high mass to gamma approximately 2.3 at low mass. The spatial cross-correlation of dwarf galaxies with more massive galaxies shows fairly similar trends, with a steeper radial dependence at low mass than predicted in numerical simulations of sub-halos within galaxy halos. To examine the issue of missing satellites we combine the cross-correlation measurements with our estimates of the low mass galaxy number density. We find on the average there are 60+/-20 dwarfs in sub-halos with M(total) > 10^7 Msun for a typical Local Group M(total)/M(stars)=30, corresponding to M/L_V approximately 100 for a galaxy with no recent star formation. The number of dwarfs per galaxy is about a factor of two larger than currently found for the Milky Way. Nevertheless, the average dwarf counts are about a factor of 30 below LCDM simulation results. The divergence from LCDM predictions is one of slope of the relation, approximately dN/dlnM approximately -0.5 rather than the predicted -0.9, not sudden onset at some characteristic scale. The dwarf galaxy star formation rates span the range from passive to bursting, which suggests that there are few completely dark halos.

Timmes, Brown & Truran found that metallicity variations could theoretically account for a 25% variation in the mass of 56Ni synthesized in Type Ia supernovae (SNe Ia), and thus account for a large fraction of the scatter in observed SN Ia luminosities. Higher-metallicity progenitors are more neutron-rich, producing more stable burning products relative to radioactive 56Ni. We develop a new method for estimating bolometric luminosity and 56Ni yield in SNe Ia and use it to test the theory with data from the Supernova Legacy Survey. We find that the average 56Ni yield does drop in SNe Ia from high metallicity environments, but the theory can only account for 7%--10% of the dispersion in SN Ia 56Ni mass, and thus luminosity. This is because the effect is dominant at metallicities significantly above solar, whereas we find that SN hosts have predominantly subsolar or only moderately above-solar metallicities. We also show that allowing for changes in O/Fe with the metallicity [Fe/H] does not have a major effect on the theoretical prediction of Timmes, Brown & Truran, so long as one is using the O/H as the independent variable. Age may have a greater effect than metallicity -- we find that the luminosity weighted age of the host galaxy is correlated with 56Ni yield, and thus more massive progenitors give rise to more luminous explosions. This is hard to understand if most SNe Ia explode when the primaries reach the Chandrasekhar mass. Finally, we test the findings of Gallagher et al., that the residuals of SNe Ia from the Hubble diagram are correlated with host galaxy metallicity, and we find no such correlation.

For the first time the cross-correlation between type Ia supernova host galaxies and surrounding field galaxies is measured using the Supernova Legacy Survey sample. Over the z=0.2 to 0.9 redshift range we find that supernova hosts are correlated an average of 60% more strongly than similarly selected field galaxies over the 3-100 arcsec range and about a factor of 3 more strongly below 10 arcsec. The correlation errors are empirically established with a jackknife analysis of the four SNLS fields. The hosts are more correlated than the field at a significance of 99% in the fitted amplitude and slope, with the point-by-point difference of the two correlation functions having a reduced $\chi^2$ for 8 degrees of freedom of 4.3, which has a probability of random occurrence of less than 3x10^{-5}. The correlation angle is 1.5+/-0.5 arcsec, which deprojects to a fixed co-moving correlation length of approximately 6.5+/- 2/h mpc. Weighting the field galaxies with the mass and star formation rate supernova frequencies of the simple A+B model produces good agreement with the observed clustering. We conclude that these supernova clustering differences are primarily the expected outcome of the dependence of supernova rates on galaxy masses and stellar populations with their clustering environment.

Fitting whole spectra at intermediate spectral resolution (R = 1000 -- 3000), to derive physical properties of stellar populations, appears as an optimized alternative to methods based on spectrophotometric indices: it uses all the redundant information contained in the signal. This paper addresses the validation of the method and it investigates the quality of the population models together with the reliability of the fitting procedures. We are using two algorithms: STECKMAP, a non-parametric regularized program and NBURSTS a parametric non-linear minimization. We compare three spectral synthesis models for single stellar populations: Pegase-HR, Galaxev (BC03) and Vazdekis/Miles, and we analyse spectra of Galactic clusters whose populations are known from studies of color-magnitude diagrams (CMD) and spectroscopy of individual stars. We find that: (1) The quality of the models critically depends on the stellar library they use. Pegase-HR and Vazdekis/Miles are consistent, while the comparison between Pegase-HR and BC03 shows some systematics reflecting the limitations of the stellar library (STELIB) used to generate the latter models; (2) The two fitting programs are consistent; (3) For globular clusters and M67 spectra, the method restitutes metallicities in agreement with spectroscopy of stars within 0.14 dex; (4) The spectroscopic ages are very sensitive to the presence of a blue horizontal branch (BHB) or of blue stragglers. A BHB morphology results in a young SSP-equivalent age. Fitting a free amount of blue stars in addition to the SSP model to mimic the BHB improves and stabilizes the fit and restores ages in agreement with CMDs studies. This method is potentially able to disentangle age or BHB effects in extragalactic clusters.

We have used the Hubble Space Telescope's Advanced Camera for Surveys (Ford et al. 2003) to measure the cumulative mass density in morphologically-selected early-type galaxies over the redshift range 0.8 < z < 1.7. Our imaging data set covers four well-separated sight-lines, and is roughly intermediate (in terms of both depth and area) between the GOODS/GEMS imaging data, and the images obtained in the Hubble Deep Field campaigns. Our images contain 144 galaxies with ultra-deep spectroscopy obtained as part of the Gemini Deep Deep Survey. These images have been analyzed using a new purpose-written morphological analysis code which improves the reliability of morphological classifications by adopting a 'quasi-Petrosian' image thresholding technique. We find that at z \~ 1 about 80% of the stars living in the most massive galaxies reside in early-type systems. This fraction is similar to that seen in the local Universe. However, we detect very rapid evolution in this fraction over the range 0.8 < z < 1.7, suggesting that over this redshift range the strong morphology-mass relationship seen in the nearby Universe is beginning to fall into place. By comparing our images to published spectroscopic classifications, we show that little ambiguity exists in connecting spectral classes to morphological classes for spectroscopically quiescent systems. However, the mass density function of early-type galaxies is evolving more rapidly than that of spectroscopically quiescent systems, which we take as further evidence that we are witnessing the formation of massive early-type galaxies over the 0.8 < z < 1.7 redshift range.

We present ELODIE.3.1, an updated release of the library published in Prugniel & Soubiran (2001, 2004). The library includes 1962 spectra of 1388 stars obtained with the ELODIE spectrograph at the Observatoire de Haute-Provence 193cm telescope in the wavelength range 390 to 680 nm. It provides a wide coverage of atmospheric parameters : T_eff from 3100 K to 50000 K, log g from -0.25 to 4.9 and [Fe/H] from -3 to +1. The library is given at two resolutions: R~42000, with the flux normalized to the pseudo-continuum, FWHM=0.55\AA (R~10000) calibrated in physical flux (reduced above earth atmosphere) with a broad-band photometric precision of 2.5% and narrow-band precision of 0.5%. In this new release the data-reduction (flux calibration, reconnection of the echelle orders) has been improved, and in particular the blue region, between 390 and 400 nm has been added. The FITS files for each spectra, and the measured atmospheric parameters are publicly available. See the ELODIE.3.1 page for more details: http://www.obs.u-bordeaux1.fr/m2a/soubiran/elodie_library.html

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

The stellar population models dramatically progressed with the arrival of large and complete libraries, ELODIE, CFLIB (=Indo-US) and MILES at a relatively high resolution. We show that the quality of the fits is not anymore limited by the size of the stellar libraries in a large range of ages (0.1 to 10 Gyrs) and metallicities (-2 to +0.4 dex). The main limitations of the empirical stellar libraries are (i) the coverage of the parameters space (lack of hot stars of low metallicity), (ii) the precision and homogeneity of the atmospheric parameters and (iii) the non-resolution of individual element abundances (in particular [$\alpha$/Fe]). Detailed abundance measurements in the large libraries, and usage of theoretical libraries are probably the next steps, and we show that a combination between an empirical (ELODIE) and a theoretical library (Coelho et al. 2005) immediately improves the modeling of ($\alpha$-enhanced) globular clusters.

We have compared simple stellar populations (SSPs) generated with different population synthesis tools: BC03, Vazdekis and Pegase.HR and different stellar libraries: ELODIE3.1, SteLib and MILES. We find that BC03/SteLib SSPs are biased toward solar metallicity, however Pegase.HR/ELODIE3.1 and Vazdekis/MILES are extremely consistent. The extensive coverage of the space of atmospheric parameters in the large stellar libraries allows precise synthesis for a large range of ages (0.1 .. 10 Gyr) and metallicities (-2 .. +0.4 dex) limited by the quality of the determination of stellar parameters (like temperature scale of the giants).

We are using full spectrum fitting to determine the ages and metallicities of Galactic clusters (M67 and globular clusters). We find that the method is very accurate to measure the metallicity. Blue horizontal branches are well identified as a 'young' sub-component, and the age of the 'old' component is in fair agreement with CMD determinations.

We describe a compact cluster of massive red galaxies at z=1.51 discovered in one of the Gemini Deep Deep Survey (GDDS) fields. Deep imaging with the Near Infrared Camera and Multi Object Spectrometer (NICMOS) on the Hubble Space Telescope reveals a high density of galaxies with red optical to near-IR colors surrounding a galaxy with a spectroscopic redshift of 1.51. Mid-IR imaging with Infrared Array Camera (IRAC) on the Spitzer Space telescope shows that these galaxies have spectral energy distributions that peak between 3.6 and 4.5 microns. Fits to 12-band photometry reveal 12 or more galaxies with spectral shapes consistent with z = 1.51. Most are within ~170 co-moving kpc of the GDDS galaxy. Deep F814W images with the Advanced Camera for Surveys (ACS) on HST reveal that these galaxies are a mix of early-type galaxies, disk galaxies and close pairs. The total stellar mass enclosed within a sphere of 170 kpc in radius is > 8E+11 solar masses. The colors of the most massive galaxies are close to those expected from passive evolution of simple stellar populations (SSP) formed at much higher redshifts. We suggest that several of these galaxies will merge to form a single, very massive galaxy by the present day. This system may represent an example of a short-lived dense group or cluster core typical of the progenitors of massive clusters in the present day and suggests the red sequence was in place in over-dense regions at early times.

We have used the Hubble Space Telescope's Advanced Camera for Surveys to measure the mass density function of morphologically selected early-type galaxies in the Gemini Deep Deep Survey fields, over the redshift range 0.9 < z < 1.6. Our imaging data set covers four well-separated sight-lines, and is roughly intermediate (in terms of both depth and area) between the GOODS/GEMS imaging data, and the images obtained in the Hubble Deep Field campaigns. Our images contain 144 galaxies with ultra-deep spectroscopy, and they have been analyzed using a new purpose-written morphological analysis code which improves the reliability of morphological classifications by adopting a 'quasi-petrosian' image thresholding technique. We find that at z = 1 approximately 70% of the stars in massive galaxies reside in early-type systems. This fraction is remarkably similar to that seen in the local Universe. However, we detect very rapid evolution in this fraction over the range 1.0 < z < 1.6, suggesting that in this epoch the strong color-morphology relationship seen in the nearby Universe is beginning to fall into place.

The acceleration of the expansion of the universe, and the need for Dark Energy, were inferred from the observations of Type Ia supernovae (SNe Ia). There is consensus that SNe Ia are thermonuclear explosions that destroy carbon-oxygen white dwarf stars that accrete matter from a companion star, although the nature of this companion remains uncertain. SNe Ia are thought to be reliable distance indicators because they have a standard amount of fuel and a uniform trigger -- they are predicted to explode when the mass of the white dwarf nears the Chandrasekhar mass -- 1.4 solar masses. Here we show that the high redshift supernova SNLS-03D3bb has an exceptionally high luminosity and low kinetic energy that both imply a super-Chandrasekhar mass progenitor. Super-Chandrasekhar mass SNe Ia should preferentially occur in a young stellar population, so this may provide an explanation for the observed trend that overluminous SNe Ia only occur in young environments. Since this supernova does not obey the relations that allow them to be calibrated as standard candles, and since no counterparts have been found at low redshift, future cosmology studies will have to consider contamination from such events.

We compare the rise times of nearby and distant Type Ia supernovae (SNe Ia) as a test for evolution using 73 high-redshift spectroscopically-confirmed SNe Ia from the first two years of the five year Supernova Legacy Survey (SNLS) and published observations of nearby SN. Because of the ``rolling'' search nature of the SNLS, our measurement is approximately 6 times more precise than previous studies, allowing for a more sensitive test of evolution between nearby and distant supernovae. Adopting a simple $t^2$ early-time model (as in previous studies), we find that the rest-frame $B$ rise times for a fiducial SN Ia at high and low redshift are consistent, with values $19.10^{+0.18}_{-0.17}({stat}) \pm 0.2 ({syst})$ and $19.58^{+0.22}_{-0.19}$ days, respectively; the statistical significance of this difference is only 1.4 \sg . The errors represent the uncertainty in the mean rather than any variation between individual SN. We also compare subsets of our high-redshift data set based on decline rate, host galaxy star formation rate, and redshift, finding no substantive evidence for any subsample dependence.

(ABRIDGED) We show that Type Ia supernovae (SNe Ia) are formed within both very young and old stellar populations, with observed rates that depend on the stellar mass and mean star-formation rates (SFRs) of their host galaxies. Models where the SN Ia rate depends solely on host galaxy stellar mass are ruled out with 99% confidence. Our analysis is based on 100 spectroscopically-confirmed SNe Ia, plus 24 photometrically-classified events, all from the Supernova Legacy Survey (SNLS) and distributed over 0.2<z<0.75. Using multi-band photometry, we estimate stellar masses and SFRs for the SN Ia host galaxies by fitting their broad-band spectral energy distributions with the galaxy spectral synthesis code, PEGASE.2. We show that the SN Ia rate per unit mass is proportional to the specific SFR of the parent galaxies -- more vigorously star-forming galaxies host more SNe Ia per unit stellar mass, broadly equivalent to the trend of increasing SN Ia rate in later-type galaxies seen in the local universe. Following earlier suggestions for a simple "two-component" model approximating the SN Ia rate, we find bivariate linear dependencies of the SN Ia rate on both the stellar masses and the mean SFRs of the host systems. We also demonstrate a dependence of distant SN Ia light-curve shapes on star-formation in the host galaxy, similar to trends observed locally. Passive galaxies, with no star-formation, preferentially host faster-declining/dimmer SNe Ia, while slower-declining/brighter events are only found in systems with ongoing star-formation. We model the light-curve width distribution in star-forming galaxies as the sum of a young component, and an old component taken from the distribution in non-star-forming galaxies.

The GRB Host Studies (GHostS) is a public archive collecting observed quantities of GRB host galaxies. At this time (January 2006) it contains information on 32 GRB hosts, i.e. about half of the total number of GRBs with known redshift. Here we present some preliminary statistical analysis of the sample, e.g. the total stellar mass, metallicity and star formation rate for the hosts. We found that these are generally low-mass objects, with 79% having M_*<10^10 M_sun. The total stellar mass and the metallicity for a subsample of 7 hosts at 0.4<z<1 are consistent with the mass-metallicity relation recently found for normal star-forming galaxies in the same redshift interval. At least 56% of the total sample are bursty galaxies: their growth time-scale (the time required to form the observed stellar mass assuming that the observed SFR is constant over the entire life of the galaxy) is shorter than 400 Myr.

We show that there has been a dramatic decline in the abundance of massive galaxies with strong Hdelta stellar absorption lines from z=1.2 to the present. These ``Hdelta-strong'', or HDS, galaxies have undergone a recent and rapid break in their star-formation activity. Combining data from the Gemini Deep Deep and the Sloan Digital Sky Surveys to make mass-matched samples (M*>=10^10.2 Msun), with 25 and 50,255 galaxies, respectively), we find that the fraction of galaxies in an HDS phase has decreased from about 50% at z=1.2 to a few percent today. This decrease in fraction is due to an actual decrease in the number density of massive HDS systems by a factor of 2-4, coupled with an increase in the number density of massive galaxies by about 30 percent. We show that this result depends only weakly on the threshold chosen for the Hdelta equivalent width to define HDS systems (if greater than 4 A) and corresponds to a (1+z)^{2.5\pm 0.7} evolution. Spectral synthesis studies of the high-redshift population using the PEGASE code, treating Hdelta_A, EW[OII], Dn4000, and rest-frame colors, favor models in which the Balmer absorption features in massive Hdelta-strong systems are the echoes of intense episodes of star-formation that faded about 1 Gyr prior to the epoch of observation. The z=1.4-2 epoch appears to correspond to a time at which massive galaxies are in transition from a mode of sustained star formation to a relatively quiescent mode with weak and rare star-formation episodes. We argue that the most likely local descendants of the distant massive HDS galaxies are passively evolving massive galaxies in the field and small groups.

We have investigated the mass-metallicity (M-Z) relation using galaxies at 0.4<z<1.0 from the Gemini Deep Deep Survey and Canada-France Redshift Survey. Deep K and z' band photometry allowed us to measure stellar masses for 69 galaxies. From a subsample of 56 galaxies, for which metallicity of the interstellar medium is also measured, we identified a strong correlation between mass and metallicity, for the first time in the distant Universe. This was possible because of the larger base line spanned by the sample in terms of metallicity (a factor of 7) and mass (a factor of 400) than in previous works. This correlation is much stronger and tighter than the luminosity-metallicity, confirming that stellar mass is a more meaningful physical parameter than luminosity. We find clear evidence for temporal evolution in the M-Z relation in the sense that at a given mass, a galaxy at z=0.7 tends to have lower metallicity than a local galaxy of similar mass. We use the z=0.1 SDSS M-Z relation, and a small sample of z=2.3 Lyman break galaxies with known mass and metallicity, to propose an empirical redshift-dependent M-Z relation, according to which the stellar mass and metallicity in small galaxies evolve for a longer time than in massive galaxies. This relation predicts that the generally metal poor damped Lyman-alpha galaxies have stellar masses of the order of 10^8.8 M_sun (with a dispersion of 0.7 dex) all the way from z=0.2 to z=4. The observed redshift evolution of the M-Z relation can be reproduced remarkably well by a simple closed-box model where the key assumption is an e-folding time for star formation which is higher or, in other words, a period of star formation that lasts longer in less massive galaxies than in more massive galaxies. Such a picture supports the downsizing scenario for galaxy formation.

One way to extract more information from the integrated light of galaxies is to improve the spectral resolution at which observations and analysis are carried out. The population synthesis code currently providing the highest spectral resolution is Pegase-HR, which was made available by D. Le Borgne et al. in 2004. Based on an empirical stellar library, it provides synthetic spectra between 4000 and 6800 A at lambda/d(lambda)=10000 for any star formation history, with or without chemical evolution. Such a resolution is particularly useful for the study of low mass galaxies, massive star clusters, or other galaxy regions with low internal velocity dispersions. After a summary of the main features of Pegase-HR and comparisons with other population synthesis codes, this paper focuses on the inversion of optical galaxy spectra. We explore the limits of what information can or can not be recovered, based on theoretical principles and extensive simulations. First applications to extragalactic objects are shown.

The ISM metallicity and the stellar mass are examined in a sample of 66 galaxies at 0.4<z<1, selected from the Gemini Deep Deep Survey (GDDS) and the Canada-France Redshift Survey (CFRS). We observe a mass-metallicity relation similar to that seen in z~0.1 SDSS galaxies, but displaced towards higher masses and/or lower metallicities. Using this sample, and a small sample of z~2.3 LBGs, a redshift dependent mass-metallicity relation is proposed which describes the observed results.

We examine the cosmic star formation rate (SFR) and its dependence on galaxy stellar mass over the redshift range 0.8 < z < 2 using data from the Gemini Deep Deep Survey (GDDS). The SFR in the most massive galaxies (M > 10^{10.8} M_sun) was six times higher at z = 2 than it is today. It drops steeply from z = 2, reaching the present day value at z ~ 1. In contrast, the SFR density of intermediate mass galaxies (10^{10.2} < M < 10^{10.8} M_sun) declines more slowly and may peak or plateau at z ~ 1.5. We use the characteristic growth time t_SFR = rho_M / rho_SFR to provide evidence of an associated transition in massive galaxies from a burst to a quiescent star formation mode at z ~ 2. Intermediate mass systems transit from burst to quiescent mode at z ~ 1, while the lowest mass objects undergo bursts throughout our redshift range. Our results show unambiguously that the formation era for galaxies was extended and proceeded from high to low mass systems. The most massive galaxies formed most of their stars in the first ~3 Gyr of cosmic history. Intermediate mass objects continued to form their dominant stellar mass for an additional ~2 Gyr, while the lowest mass systems have been forming over the whole cosmic epoch spanned by the GDDS. This view of galaxy formation clearly supports `downsizing' in the SFR where the most massive galaxies form first and galaxy formation proceeds from larger to smaller mass scales.

We present PEGASE-HR, a new stellar population synthesis program generating high resolution spectra (R=10 000) over the optical range lambda=400--680 nm. It links the spectro-photometric model of galaxy evolution PEGASE.2 (Fioc & Rocca-Volmerange 1997) to an updated version of the ELODIE library of stellar spectra observed with the 193 cm telescope at the Observatoire de Haute-Provence (Prugniel & Soubiran 2001a). The ELODIE star set gives a fairly complete coverage of the Hertzprung-Russell (HR) diagram and makes it possible to synthesize populations in the range [Fe/H]=-2 to +0.4. This code is an exceptional tool for exploring signatures of metallicity, age, and kinematics. We focus on a detailed study of the sensitivity to age and metallicity of the high-resolution stellar absorption lines and of the classical metallic indices proposed until now to solve the age-metallicity degeneracy. Validity tests on several stellar lines are performed by comparing our predictions for Lick indices to the models of other groups. The comparison with the lower resolution library BaSeL (Lejeune et al. 1997) confirms the quality of the ELODIE library when used for simple stellar populations (SSPs) from 10 Myr to 20 Gyr. Predictions for the evolved populations of globular clusters and elliptical galaxies are given and compared to observational data. Two new high-resolution indices are proposed around the Hgamma line. They should prove useful in the analysis of spectra from the new generation of telescopes and spectrographs.

We present spectroscopic evidence from the Gemini Deep Deep Survey (GDDS) for a significant population of color-selected red galaxies at 1.3 < z < 2.2 whose integrated light is dominated by evolved stars. Unlike radio-selected objects, the z > 1.5 old galaxies have a sky density > 0.1 per sq. arcmin. Conservative age estimates for 20 galaxies with z > 1.3; <z> = 1.49, give a median age of 1.2 Gyr and <z_f> = 2.4. One quarter of the galaxies have inferred z_f > 4. Models restricted to abundances less than or equal to solar give median ages and z_f of 2.3 Gyr and 3.3, respectively. These galaxies are among the most massive and contribute approximately 50% of the stellar mass density at 1 < z < 2. The derived ages and most probable star formation histories suggest a high star-formation-rate (300-500 solar masses per year) phase in the progenitor population. We argue that most of the red galaxies are not descendants of the typical z=3 Lyman break galaxies. Galaxies associated with luminous sub-mm sources have the requisite star formation rates to be the progenitor population. Our results point toward early and rapid formation for a significant fraction of present day massive galaxies.

The narrow K-z relation of powerful radio galaxies in the Hubble K diagram is often attributed to the stellar populations of massive elliptical galaxies. Extended over a large range of redshifts (0<z<4), it is difficult to estimate masses at high z from galaxy evolution. In the present paper, we propose to estimate the stellar masses of galaxies using the galaxy evolution model PEGASE. We use star formation scenarios that successfully fit faint galaxy counts and z=0 galaxy templates. They also allow to estimate valid photometric redshifts. The baryonic mass of the initial gas cloud M_bar is then derived. The K-z relation is remarkably reproduced by our scenario for elliptical galaxies of baryonic mass M_(bar,max)=10^12 Msol, at all z up to z=4. M_(bar,max) is also the maximum mass limit of all types of galaxies. Using another IMF, even a top-heavy one, does not alter our conclusions. The high value of M_(bar,max) observed at z > 4 implies that massive clouds were already formed at early epochs. We also find that the M_(bar,max) limit is similar to the critical mass M_crit of a self-gravitating cloud regulated by cooling (Rees & Ostriker, 1977; Silk, 1977). Moreover, the critical size r_crit = 75 Kpc is remarkably close to the typical diameter of Ly_alpha haloes surrounding distant radio galaxies. This confirms the validity of the method of baryonic mass determination based on the K-band luminosity. A puzzling question that remains to be answered is the short time-scale of mass-accumulation required to form such massive galaxies at z=4. We discuss the dispersion of the K-z relation and the link between the active nucleus and a large stellar mass.

Photometric redshifts are estimated on the basis of template scenarios with the help of the code ZPEG, an extension of the galaxy evolution model PEGASE.2 and available on the PEGASE web site. The spectral energy distribution (SED) templates are computed for nine spectral types including starburst, irregular, spiral and elliptical. Dust, extinction and metal effects are coherently taken into account, depending on evolution scenarios. The sensitivity of results to adding near-infrared colors and IGM absorption is analyzed. A comparison with results of other models without evolution measures the evolution factor which systematically increases the estimated photometric redshift values by $\Delta z$ > 0.2 for z > 1.5. Moreover we systematically check that the evolution scenarios match observational standard templates of nearby galaxies, implying an age constraint of the stellar population at z=0 for each type. The respect of this constraint makes it possible to significantly improve the accuracy of photometric redshifts by decreasing the well-known degeneracy problem. The method is applied to the HDF-N sample. From fits on SED templates by a $\chi^2$-minimization procedure, not only is the photometric redshift derived but also the corresponding spectral type and the formation redshift $z_for$ when stars first formed. Early epochs of galaxy formation z > 5 are found from this new method and results are compared to faint galaxy count interpretations. The new tool is available at: http://www.iap.fr/pegase