A code computing consistently the evolution of stars, gas and dust, as well as the energy they radiate, is required to derive reliably the history of galaxies by fitting synthetic SEDs to multiwavelength observations. The new code P\'egase.3 described in this paper extends to the far-IR/submm the UV-to-near-IR modeling provided by previous versions of P\'egase. It first computes the properties of single stellar populations at various metallicities. It then follows the evolution of the stellar light of a galaxy and the abundances of the main metals in the ISM, assuming some scenario of mass assembly and star formation. It simultaneously calculates the masses of the various grain families, the optical depth of the galaxy and the attenuation of the SED through the diffuse ISM in spiral and spheroidal galaxies, using grids of radiative transfer precomputed with Monte Carlo simulations taking scattering into account. The code determines the mean radiation field and the temperature probability distribution of stochastically heated individual grains. It then sums up their spectra to yield the overall emission by dust in the diffuse ISM. The nebular emission of the galaxy is also computed, and a simple modeling of the effects of dust on the SED of star-forming regions is implemented. The main outputs are UV-to-submm SEDs of galaxies from their birth up to 20 Gyr, colors, masses of galactic components, ISM abundances of metallic elements and dust species, supernova rates. The temperatures and spectra of individual grains are also available. The paper discusses several of these outputs for a scenario representative of Milky Way-like spirals. P\'egase.3 is fully documented and its Fortran 95 source files are public. The code should be especially useful for cosmological simulations and to interpret future mid- and far-IR data, whether obtained by JWST, LSST, Euclid or e-ELT.

P\'egase.3 is a Fortran 95 code modeling the spectral evolution of galaxies from the far-ultraviolet to submillimeter wavelengths. It also follows the chemical evolution of their stars, gas and dust. For a given scenario (a set of parameters defining the history of mass assembly, the star formation law, the initial mass function...), P\'egase.3 consistently computes the following: * the star formation, infall, outflow and supernova rates from 0 to 20 Gyr; * the stellar metallicity, the abundances of main elements in the gas and the composition of dust; * the unattenuated stellar spectral energy distribution (SED); * the nebular SED, using nebular continua and emission lines precomputed with code Cloudy (Ferland et al. 2017); * the attenuation in star-forming clouds and the diffuse interstellar medium, by absorption and scattering on dust grains, of the stellar and nebular SEDs. For this, the code uses grids of the transmittance for spiral and spheroidal galaxies. We precomputed these grids through Monte Carlo simulations of radiative transfer based on the method of virtual interactions; * the re-emission by grains of the light they absorbed, taking into account stochastic heating. The main innovation compared to P\'egase.2 is the modeling of dust emission and its evolution. The computation of nebular emission has also been entirely upgraded to take into account metallicity effects and infrared lines. Other major differences are that complex scenarios of evolution (derived for instance from cosmological simulations), with several episodes of star formation, infall or outflow, may now be implemented, and that the detailed evolution of the most important elements -- not only the overall metallicity -- is followed.

The z=8 galaxy sample (EGSY-2008532660, EGS-zs8-1, EGS-zs8-2), observed with HST and Spitzer/IRAC, was selected for its atypical [3.6]-[4.5]mum IR excess. The best-fitting of observations is processed with multiple libraries of continuous UV-to-farIR spectral energy distributions (SEDs) with the help of the code Pegase.3 with dust (Fioc and Rocca-Volmerange, 2018, submitted to AA). The main results are: i) the IR excess is due to the Balmer break of red giant stars with galaxy ages of 450-500 Myrs ii) masses are high: 7 10^10 Msun for stellar mass, 10^11 Msun for galaxy mass, 3 10^8 Msun for dust mass iii) they are probing mature (no starbursting) galaxies only with the simulations of spheroidal (and no spiral) geometries and the star formation history (SFH ) of massive ellipticals. Less constrained by observations not covering the Balmer zone, the galaxy GN-z11 is consistent. These results constrain the rapid and efficient formation of massive elliptical galaxies, significantly distinct of early-type disk-spiral galaxies, and the reionization of the Universe

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We present a study of the complete ultraviolet to submillimetre spectral energy distributions (SEDs) of twelve 3CR radio galaxy hosts in the redshift range $1.0 < z < 2.5$, which were all detected in the far-infrared by the Herschel Space Observatory. The study employs the new spectro-chemical evolutionary code P\'EGASE.3, in combination with recently published clumpy AGN torus models. We uncover the properties of the massive host galaxy stellar populations, the AGN torus luminosities, and the properties of the recent starbursts, which had earlier been inferred in these objects from their infrared SEDs. The P\'EGASE.3 fitting yields very luminous (up to 10$^{13}$L$_{\odot}$) young stellar populations with ages of several hundred million years in hosts with masses exceeding 10$^{11}$M$_{\odot}$. Dust masses are seen to increase with redshift, and a surprising correlation -- or better upper envelope behaviour -- is found between the AGN torus luminosity and the starburst luminosity, as revealed by their associated dust components. The latter consistently exceeds the former by a constant factor, over a range of one order of magnitude in both quantities.

High-redshift radio galaxies present signs of both star formation and AGN activity, making them ideal candidates to investigate the connection and coevolution of AGN and star formation in the progenitors of present-day massive galaxies. We make use of a sample of 11 powerful radio galaxies spanning 1<z<4 which have complete coverage of their spectral energy distribution (SED) from UV to FIR wavelengths. Using Herschel data, we disentangle the relative contribution of the AGN and star formation by combining the galaxyevolutioncodePEGASE.3 with an AGN torus model. We find that three components are necessary to reproduce the observed SEDs: an evolved and massive stellar component, a submm bright young starburst, and an AGN torus. We find that powerful radio galaxies form at very high-redshift, but experience episodic and important growth at 1<z<4 as the mass of the associated starburst varies from 5 to 50% of the total mass of the system. The properties of star formation differ from source to source, indicating no general trend of the star formation properties in the most infrared luminous high-redshift radio galaxies and no correlation with the AGN bolometric luminosity. Moreover, we find that AGN scattered light have a very limited impact on broad-band SED fitting on our sample. Finally, our analysis also suggests a wide range in origins for the observed star formation,which we partially constrain for some sources.

In this paper, we show that the supernova remnant (SNR) masses cumulated from core-collapse supernovae along the star formation history of two powerful z=3.8 radio galaxies 4C41.17 and TN J2007-1316 reach up to > 10^9 Msun, comparable with supermassive black hole (SMBH) masses measured from the SDSS sample at similar redshifts. The SNR mass is measured from the already exploded supernova mass after subtraction of ejecta at the galaxy age where the mass of still luminous stars fits at best the observed spectral energy distribution (SED), continuously extended to the optical-Spitzer-Herschel-submm domains, with the help of the galaxy evolution model P\'egase.3. For the recent and old stellar populations, SNR masses vary on 10^(9 to 10) Msun and the SNR-to-star mass ratio between 1 and 0.1 percent is comparable to the observed low-z SMBH-to-star mass ratio. For the template radio galaxy 4C41.17, SNR and stellar population masses estimated from large aperture (>4arcsec=30kpc) observations are compatible, within one mass order, with the total mass of multiple optical HST (~700pc) structures, associated with VLA radio emissions, both at 0.1 arcsec. Probing the SNR accretion by central black holes is a simple explanation for SMBH growth, requiring physics on star formation, stellar and galaxy dynamics with consequences on various processes (quenching, mergers, negative feedback) and a key to the relation bulge-SMBH.

We present results from a survey of 70 radio galaxies (RGs) at redshifts 1<z<5.2 using the PACS and SPIRE on-board Herschel. Combined with existing mid-IR photometry from Spitzer and observations obtained with LABOCA, the SEDs of galaxies in our sample are continuously covered across 3.6-870um. The total infrared luminosities of these RGs are such that they almost all are either ultra-or hyper-luminous infrared galaxies. We fit the infrared SEDs with a set of empirical templates which represent dust heated (1) by a variety of SB and (2) by a AGN. We find that the SEDs of RGs require the dust to be heated by both AGN and SB, but the luminosities of these two components are not strongly correlated. Assuming empirical relations and simple physical assumptions, we calculate the SFR, the black hole mass accretion rate (MdotBH), and the black hole mass (MBH) for each RG. We find that the host galaxies and their BHs are growing extremely rapidly, having SFR~100-5000 Msun/yr and MdotBH~1-100 Msun/yr. The mean sSFR of RGs at z>2.5 are higher than the sSFR of typical star-forming galaxies over the same redshift range but are similar or perhaps lower than the galaxy population for RGs at z<2.5. By comparing the sSFR and the specific black hole mass accretion rate, we conclude that BHs in radio loud AGN are already, or soon will be, overly massive compared to their host galaxies in terms of expectations from the local MBH-MGal relation. In order to ``catch up'' with the BH, the galaxies require about an order-of magnitude more time to grow in mass, at the observed SFRs, compared to the time the BH is actively accreting. However, during the current cycle of activity, we argue that this catching-up is likely to be difficult due to the short gas depletion times. Finally, we speculate on how the host galaxies might grow sufficiently in stellar mass to ultimately fall onto the local MBH-MGal relation.

The Gaia satellite will survey the entire celestial sphere down to 20th magnitude, obtaining astrometry, photometry, and low resolution spectrophotometry on one billion astronomical sources, plus radial velocities for over one hundred million stars. Its main objective is to take a census of the stellar content of our Galaxy, with the goal of revealing its formation and evolution. Gaia's unique feature is the measurement of parallaxes and proper motions with hitherto unparalleled accuracy for many objects. As a survey, the physical properties of most of these objects are unknown. Here we describe the data analysis system put together by the Gaia consortium to classify these objects and to infer their astrophysical properties using the satellite's data. This system covers single stars, (unresolved) binary stars, quasars, and galaxies, all covering a wide parameter space. Multiple methods are used for many types of stars, producing multiple results for the end user according to different models and assumptions. Prior to its application to real Gaia data the accuracy of these methods cannot be assessed definitively. But as an example of the current performance, we can attain internal accuracies (RMS residuals) on F,G,K,M dwarfs and giants at G=15 (V=15-17) for a wide range of metallicites and interstellar extinctions of around 100K in effective temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]), and 0.25dex in surface gravity (logg). The accuracy is a strong function of the parameters themselves, varying by a factor of more than two up or down over this parameter range. After its launch in November 2013, Gaia will nominally observe for five years, during which the system we describe will continue to evolve in light of experience with the real data.

The new evolutionary code Pegase.3 is used for the spectral synthesis of continuous optical-Spitzer-Herschel-submm spectral energy distributions (SEDs) of two z=3.8 radio galaxies 4C 41.17 and TN J2007-1316. Both galaxies are selected from the Herschel Radio Galaxy Evolution Project for their faint AGN contribution and because they show evidence of a large stellar contribution to their bolometric luminosity. Pegase.3 coherently models the reprocessing of the stellar luminosity to dust emission, allowing to build UV to IR-submm libraries by types. These templates allow to predict SEDs at redshifts z in the observer's frame by using z=0 templates corrected for cosmology and evolution. Our best fits are a sum of two distinct evolving stellar populations: i) a very massive starburst observed ~30 Myrs after formation from optical and far-IR peaks ii) a significantly older (~1Gyr old) massive early-type population revealed from mid-IR Spitzer data. The AGN emission fits a roughly thermal model. These results confirm that many of the stellar populations in high-z radio galaxies were already formed by massive (a few thousands of billion solar masses) starbursts at z=4, confirming the previous K-z diagram interpretation. Gas-rich mergers and/or jet-cloud interactions are favored for triggering the intense star formation, possibly in relation with massive super black holes of these radio galaxies. These multiple stellar populations may be a generic feature of the luminous infrared radio galaxy population as a whole and strongly constrain mass accumulation and galaxy formation models.

We examine the relative orientation of radio jets and dusty tori surrounding the AGN in powerful radio galaxies at z > 1. The radio core dominance R = P(20 GHz) /P(500 MHz) serves as an orientation indicator, measuring the ratio between the anisotropic Doppler-beamed core extended core emission and the isotropic lobe emission. Assuming a fixed cylindrical geometry for the hot, dusty torus, we derive its inclination i by fitting optically-thick radiative transfer models to spectral energy distributions obtained with the Spitzer Space Telescope. We find a highly significant anti-correlation (p < 0.0001) between R and i in our sample of 35 type 2 AGN combined with a sample of 18 z = 1 3CR sources containing both type 1 and 2 AGN. This analysis provides observational evidence both for the Unified scheme of AGN and for the common assumption that radio jets are in general perpendicular to the plane of the torus. The use of inclinations derived from mid-infrared photometry breaks several degeneracies which have been problematic in earlier analyses. We illustrate this by deriving the core Lorentz factor Gamma from the R-i anti-correlation, finding Gamma > 1.3.

We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.

Context. This study has been developed in the framework of the computational simulations executed for the preparation of the ESA Gaia astrometric mission. Aims. We focus on describing the objects and characteristics that Gaia will potentially observe without taking into consideration instrumental effects (detection efficiency, observing errors). Methods. The theoretical Universe Model prepared for the Gaia simulation has been statistically analyzed at a given time. Ingredients of the model are described, giving most attention to the stellar content, the double and multiple stars, and variability. Results. In this simulation the errors have not been included yet. Hence we estimate the number of objects and their theoretical photometric, astrometric and spectroscopic characteristics in the case that they are perfectly detected. We show that Gaia will be able to potentially observe 1.1 billion of stars (single or part of multiple star systems) of which about 2% are variable stars, 3% have one or two exoplanets. At the extragalactic level, observations will be potentially composed by several millions of galaxies, half million to 1 million of quasars and about 50,000 supernovas that will occur during the 5 years of mission. The simulated catalogue will be made publicly available by the DPAC on the Gaia portal of the ESA web site http://www.rssd.esa.int/gaia/.

Aims:This paper is the third in a series implementing a classification system for Gaia observations of unresolved galaxies. The system makes use of template galaxy spectra in order to determine spectral classes and estimate intrinsic astrophysical parameters. In previous work we used synthetic galaxy spectra produced by PEGASE.2 code to simulate Gaia observations and to test the performance of Support Vector Machine (SVM) classifiers and parametrizers. Here we produce a semi-empirical library of galaxy spectra by fitting SDSS spectra with the previously produced synthetic libraries. We present (1) the semi-empirical library of galaxy spectra, (2) a comparison between the observed and synthetic spectra, and (3) first results of claassification and parametrization experiments with simulated Gaia spectrophotometry of this library. Methods: We use chi2-fitting to fit SDSS galaxy spectra with the synthetic library in order to construct a semi-empirical library of galaxy spectra in which (1) the real spectra are extended by the synthetic ones in order to cover the full wavelength range of Gaia, and (2) astrophysical parameters are assigned to the SDSS spectra by the best fitting synthetic spectrum. The SVM models were trained with and applied to semi-empirical spectra. Tests were performed for the classification of spectral types and the estimation of the most significant galaxy parameters (in particular redshift, mass to light ratio and star formation history). Results: We produce a semi-empirical library of 33670 galaxy spectra covering the wavelength range 250 to 1050 nm at a sampling of 1 nm or less. Using the results of the fitting of the SDSS spectra with our synthetic library, we investigate the range of the input model parameters that produces spectra which are in good agreement with observations. (abridged)

We present results from a comprehensive imaging survey of 70 radio galaxies at redshifts 1<z<5.2 using all three cameras onboard the Spitzer Space Telescope. The resulting spectral energy distributions unambiguously show a stellar population in 46 sources and hot dust emission associated with the active nucleus in 59. Using a new restframe S_3um/S_1.6um versus S_um/S_3um criterion, we identify 42 sources where the restframe 1.6um emission from the stellar population can be measured. For these radio galaxies, the median stellar mass is high, 2x10^11 M_sun, and remarkably constant within the range 1<z<3. At z>3, there is tentative evidence for a factor of two decrease in stellar mass. This suggests that radio galaxies have assembled the bulk of their stellar mass by z~3, but confirmation by more detailed decomposition of stellar and AGN emission is needed. The restframe 500 MHz radio luminosities are only marginally correlated with stellar mass but are strongly correlated with the restframe 5um hot dust luminosity. This suggests that the radio galaxies have a large range of Eddington ratios. We also present new Very Large Array 4.86 and 8.46 GHz imaging of 14 radio galaxies and find that radio core dominance --- an indicator of jet orientation --- is strongly correlated with hot dust luminosity. While all of our targets were selected as narrow-lined, type 2 AGNs, this result can be understood in the context of orientation-dependent models if there is a continuous distribution of orientations from obscured type 2 to unobscured type 1 AGNs rather than a clear dichotomy. Finally, four radio galaxies have nearby (<6") companions whose mid-IR colors are suggestive of their being AGNs. This may indicate an association between radio galaxy activity and major mergers.

This paper is the second in a series, implementing a classification system for Gaia observations of unresolved galaxies. Our goals are to determine spectral classes and estimate intrinsic astrophysical parameters via synthetic templates. Here we describe (1) a new extended library of synthetic galaxy spectra, (2) its comparison with various observations, and (3) first results of classification and parametrization experiments using simulated Gaia spectrophotometry of this library. Using the PEGASE.2 code, based on galaxy evolution models that take account of metallicity evolution, extinction correction, and emission lines (with stellar spectra based on the BaSeL library), we improved our first library and extended it to cover the domain of most of the SDSS catalogue. We produce an extended library of 28885 synthetic galaxy spectra at zero redshift covering four general Hubble types of galaxies, over the wavelength range between 250 and 1050 nm at a sampling of 1 nm or less. The library is also produced for 4 random values of redshift in the range of 0-0.2. It is computed on a random grid of four key astrophysical parameters (infall timescale and 3 parameters defining the SFR) and, depending on the galaxy type, on two values of the age of the galaxy. The synthetic library was compared and found to be in good agreement with various observations. The first results from the SVM classifiers and parametrizers are promising, indicating that Hubble types can be reliably predicted and several parameters estimated with low bias and variance.

We present a detailed reduction of a mid-infrared 12um (LW10 filter) ISOCAM open time observation performed on the ESO-Sculptor Survey field (Arnouts et al. 1997). A complete catalogue of 142 sources (120 galaxies and 22 stars), detected with high significance (equivalent to 5sigma), is presented above an integrated flux density of 0.24mJy. Star/galaxy separation is performed by a detailed study of colour-colour diagrams. The catalogue is complete to 1mJy and below this flux density the incompleteness is corrected using two independent methods. The first method uses stars and the second uses optical counterparts of the ISOCAM galaxies; these methods yield consistent results. We also apply an empirical flux density calibration using stars in the field. For each star, the 12um flux density is derived by fitting optical colours from a multi-band chi^2 to stellar templates (BaSel-2.0) and using empirical optical-IR colour-colour relations. This article is a companion analysis to Rocca-Volmerange 2007 et al. where the 12um faint galaxy counts are presented and analysed by galaxy type with the evolutionary code PEGASE.3.

Aims:The Gaia astrometric survey mission will, as a consequence of its scanning law, obtain low resolution optical (330-1000 nm) spectrophotometry of several million unresolved galaxies brighter than V=22. We present the first steps in a project to design and implement a classification system for these data. The goal is both to determine morphological classes and to estimate intrinsic astrophysical parameters via synthetic templates. Here we describe (1) a new library of synthetic galaxy spectra, and (2) first results of classification and parametrization experiments using simulated Gaia spectrophotometry of this library. Methods:We have created a large grid of synthetic galaxy spectra using the PEGASE.2 code, which is based on galaxy evolution models that take into account metallicity evolution, extinction correction, emission lines (with stellar spectra based on the BaSeL library). Our classification and regression models are Support Vector Machines (SVMs), which are kernel-based nonlinear estimators. Results:We produce a basic library of about 4000 zero redshift galaxy spectra covering the main Hubble types over wavelength range 250 to 1050 nm at a sampling of 1 nm or less. It is computed on a regular grid of four key astrophysical parameters for each type and for intermediate random values of the same parameters. An extended library reproduces this at a series of redshifts. Initial results from the SVM classifiers and parametrizers are promising, indicating that Hubble types can be reliably predicted and several parameters estimated with low bias and variance. Comparing the colours of our synthetic library with Sloan Digital Sky Survey (SDSS) spectra we find good agreement over the full range of Hubble types and parameters.

Multi-wavelength galaxy number counts provide clues on the nature of galaxy evolution. The interpretation per galaxy type with the new code PEGASE.3 of the 12mu (ISO) and 24mu (Spitzer) faint galaxy counts provide new constraints on the dust and stellar emission. It also reveals the nature of ULIRGs (L/Lsun > 10^12) and informs on the cosmic star formation history and the time-scales for mass build-up. We firstly present the faint galaxy counts at 12um derived from the catalogue of the ISO-ESO-Sculptor Survey (ISO-ESS) in a companion paper (Seymour et al. 2007) which go down to 0.24 mJy after corrections for incompleteness. We check that they are consistent with ISO number counts at 15um. Secondly we show that the ``normal'' scenarios which fit the deep UV-optical-near-IR counts with PEGASE.2 (Fioc et al. 1999), are unsuccessful in modelling the strong excess simultaneously observed in the cumulative and differential counts at 12um, 15um and 24um. Based on observed 12mu and 25mu IRAS luminosity functions and optical/mid-IR colors, we finally succeed in modelling cumulative and differential counts by only changing 9% of normal galaxies (1/3 of the ellipticals) into ultra-bright dusty elliptical galaxies, interpreted as ULIRGs. This ULIRG population has similarities with high-z radio-galaxy hosts. No number density evolution is included. The Herschel observatory will hopefully confirm these results.

We present the results of a comprehensive Spitzer survey of 69 radio galaxies across 1<z<5.2. Using IRAC (3.6-8.0um), IRS (16um) and MIPS (24-160um) imaging, we decompose the rest-frame optical to infrared spectral energy distributions into stellar, AGN, and dust components and determine the contribution of host galaxy stellar emission at rest-frame H-band. Stellar masses derived from rest-frame near-IR data, where AGN and young star contributions are minimized, are significantly more reliable than those derived from rest-frame optical and UV data. We find that the fraction of emitted light at rest-frame H-band from stars is >60% for ~75% the high redshift radio galaxies. As expected from unified models of AGN, the stellar fraction of the rest-frame H-band luminosity has no correlation with redshift, radio luminosity, or rest-frame mid-IR (5um) luminosity. Additionally, while the stellar H-band luminosity does not vary with stellar fraction, the total H-band luminosity anti-correlates with the stellar fraction as would be expected if the underlying hosts of these radio galaxies comprise a homogeneous population. The resultant stellar luminosities imply stellar masses of 10^{11-11.5}Msun even at the highest redshifts. Powerful radio galaxies tend to lie in a similar region of mid-IR color-color space as unobscured AGN, despite the stellar contribution to their mid-IR SEDs at shorter-wavelengths. The mid-IR luminosities alone classify most HzRGs as LIRGs or ULIRGs with even higher total-IR luminosities. As expected, these exceptionally high mid-IR luminosities are consistent with an obscured, highly-accreting AGN. We find a weak correlation of stellar mass with radio luminosity.

The purpose of this study is to explore the relationship between galaxy stellar masses, based on multiwavelength photometry spectral template fitting and dynamical masses based on published velocity dispersion measurements, for a sample of 48 early-type galaxies at z ~ 1 with HST/ACS morphological information. We determine photometric-stellar masses and perform a quantitative morphological analysis of cluster and field galaxies at redshift 0.6 < z < 1.2, using ground- and space-based multiwavelegth data available on the GOODS-S field and on the field around the X-ray luminous cluster RDCS1252.9-2927 at z = 1.24. We use multi-band photometry over 0.4-8um from HST/ACS, VLT/ISAAC and Spitzer/IRAC to estimate photometric-stellar masses using Composite Stellar Population (CSP) templates computed with PEGASE.2 models. We compare stellar masses with those obtained using CSPs built with Bruzual & Charlot and Maraston models. We then compare photometric-stellar mass and dynamical mass estimates as a function of morphological parameters obtained from HST/ACS imaging. Based on our sample, which spans the mass range log(Mphot)=[10, 11.5], we find that 1) PEGASE.2, BC03, M05 yield consistent photometric-stellar masses for early-type galaxies at z ~ 1 with a small scatter (0.15 dex rms); 2) adopting a Kroupa IMF, photometric-stellar masses match dynamical mass estimates for early-type galaxies with an average offset of 0.27 dex; 3) assuming a costant IMF, increasing dark matter fraction with the increasing galaxy mass can explain the observed trend.

This is the second in a series of papers presenting observations and results for a sample of 76 ultra-steep-spectrum (USS) radio sources in the southern hemisphere designed to find galaxies at high redshift. Here we focus on the optical spectroscopy program for 53 galaxies in the sample. We report 35 spectroscopic redshifts, based on observations with the Very Large Telescope (VLT), the New Technology Telescope (NTT) and the Australian National University's 2.3m telescope; they include five radio galaxies with z>3. Spectroscopic redshifts for the remaining 18 galaxies could not be confirmed: three are occulted by Galactic stars, eight show continuum emission but no discernible spectral lines, whilst the remaining seven galaxies are undetected in medium-deep VLT integrations. The latter are either at very high redshift (z >~7) or heavily obscured by dust. A discussion of the efficiency of the USS technique is presented. Based on the similar space density of z>3 radio galaxies in our sample compared with other USS-selected samples, we argue that USS selection at 843-1400 MHz is an efficient and reliable technique for finding distant radio galaxies.

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 spectra of the powerful 3CR radio galaxies present a typical distribution in the far-infrared. From the observed radio to X-ray spectral energy distribution (SED) templates, we propose to subtract the typical energy distributions of, respectively, the elliptical galaxy host and the synchrotron radiation. The resulting SED reveals that the main dust emission is well fitted by the sum of two blackbody components at the respective temperatures 340K+-50K and 40K+-16K. When the AGN is active, the energy rate released by hot dust is much more dissipative than cold dust and stellar emission, even when the elliptical galaxy emission is maximum at age of ~ 90 Myr. Hot dust appears as a huge cooling source which implies an extremely short time-scale t(cool). In balance, with the short gravitational time-scale t(grav) of massive galaxies, the dissipative self-gravitational models (Rees & Ostriker, 1977) are favoured for radio sources. They justify the existence of massive radio galaxies at z=4 (Rocca-Volmerange et al., 2004). The synchrotron emission is emitting up to the X-ray wavelength range, so that strong "EXOs" sources could be assimilated to 3CR radio sources. This analysis applied to ISO and SPITZER data on a larger sample will statistically confirm these results.

We present a re--analysis of the European Large Area ISO Survey (ELAIS) 90mum observations carried out with ISOPHOT, an instrument on board the ESA's Infrared Space Observatory (ISO). With more than 12 sq. deg., the ELAIS survey is the largest area covered by ISO in a single program and is about one order of magnitude deeper than the IRAS 100mum survey. The data analysis is presented and was mainly performed with the Phot Interactive Analysis software but using the pairwise method of Stickel et al. (2003) for signal processing from ERD (Edited Raw Data) to SCP (Signal per Chopper Plateau). The ELAIS 90mum catalogue contains 229 reliable sources with fluxes larger than 70 mJy and is available at www.blackwell-synergy.com. Number counts are presented and show an excess above the no-evolution model prediction. This confirms the strong evolution detected at shorter(15mum) and longer (170mum) wavelengths in other ISO surveys. The ELAIS counts are in agreement with previous works at 90mum and in particular with the deeper counts extracted from the Lockman hole observations. Comparison with recent evolutionary models show that the models of Franceschini et al. and Guiderdoni et al. which includes a heavily-extinguished population of galaxies give the best fit to the data. Deeper observations are nevertheless required to better discriminate between the model predictions in the far-infrared and are scheduled with the Spitzer Space Telescope which already started operating and will also be performed by ASTRO-F.

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.

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.

We present the first results from a pilot study to search for distant radio galaxies in the southern hemisphere (delta < -32). Within a 360 deg^2 region of sky, we define a sample of 76 ultra-steep spectrum (USS) radio sources from the 843 MHz Sydney University Molonglo Sky Survey (SUMSS) and 1.4 GHz NRAO VLA Sky Survey (NVSS) radio surveys with alpha_843^1400 < -1.3 and S_1400 > 15 mJy. We observed 71 sources without bright optical or near-infrared counterparts at 1.385 GHz with the ATCA, providing ~5" resolution images and sub-arcsec positional accuracy. To identify their host galaxies, we obtained near-IR K-band images with IRIS2 at the AAT and SofI at the NTT. We identify 92% of the USS sources down to K~20.5. The SUMSS-NVSS USS sources have a surface density more than 4 times higher than USS sources selected at lower frequencies. This is due to the higher effective selection frequency, and the well-matched resolutions of both surveys constructed using the same source fitting algorithm. The scattering of alpha >-1.3 sources into the USS sample due to spectral index uncertainties can account for only 35% of the observed USS sources. Since our sample appears to contain a similar fraction of very distant (z>3) galaxies, selecting USS sources from SUMSS-NVSS should allow us to identify large numbers of massive galaxies at high redshift.

We present STELIB, a new spectroscopic stellar library, available at http://webast.ast.obs-mip.fr/stelib . STELIB consists of an homogeneous library of 249 stellar spectra in the visible range (3200 to 9500A), with an intermediate spectral resolution (~3A) and sampling (1A). This library includes stars of various spectral types and luminosity classes, spanning a relatively wide range in metallicity. The spectral resolution, wavelength and spectral type coverage of this library represents a substantial improvement over previous libraries used in population synthesis models. The overall absolute photometric uncertainty is 3%.

We present results of a comprehensive multi-frequency study of the radio galaxy B3 J2330+3927. The 1.9" wide radio source, consisting of 3 components, is bracketed by 2 objects in our Keck K-band image. Optical and near-IR Keck spectroscopy of these two objects yield z=3.087+-0.004. The brightest (K=18.8) object has a standard type II AGN spectrum, and is the most likely location of the AGN, which implies a one-sided jet radio morphology. Deep 113 GHz observations with the IRAM Plateau de Bure Interferometer reveal CO J=4-3 emission, which peaks at the position of the AGN. The CO line is offset by 500 km/s from the systemic redshift of the AGN, but corresponds very closely to the velocity shift of an associated HI absorber seen in Lya. This strongly suggests that both originate from the same gas reservoir surrounding the AGN host galaxy. Simultaneous 230 GHz interferometer observations find a ~3x lower integrated flux density when compared to single dish 250 GHz observations with MAMBO at the IRAM 30m telescope. This can be interpreted as spatially resolved thermal dust emission at scales of 0.5" to 6". Finally, we present a tau <1.3% limit to the HI 21 cm absorption against the radio source, which represents the seventh non-detection out of 8 z>2 radio galaxies observed to date with the WSRT. We present mass estimates for the atomic, neutral, and ionized hydrogen, and for the dust, ranging from M(HI)=2x10^7 M_Sun derived from the associated HI absorber in Lya up to M(H_2)=7x10^{10} M_Sun derived from the CO emission. This indicates that the host galaxy is surrounded by a massive reservoir of gas and dust. The K-band companion objects may be concentrations within this reservoir, which will eventually merge with the central galaxy hosting the AGN.

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

We have analyzed the ultraviolet and optical emission line ratios of a large sample of extragalactic radio sources (QSOs and radio galaxies), with the help of models combining AGN photoionization and shocks. The results strongly suggest that the two ionizing mechanisms frequently coexist. The model sequences obtained by varying the balance between shocks and AGN photoionization account for most emission line data in the 12 line ratio diagrams we have considered. In the frequently used diagrams involving [OIII]5007/Hbeta, [OI]6300/Halpha, [NII]6584/Halpha and [SII]6716,6731, the effect of varying the shock-photoionization balance mimics a variation of the ionization parameter (U) in traditional photoionization sequences. In most of the remaining diagrams, such as [OI]6300/[OIII]5007 vs.[OIII]4363/[OIII]5007, [OIII]5007/H vs. [OIII]4363/[OIII]5007 and CIII]1909/[CII]2326 vs. CIV1549/[CII]2326, the data can only be accounted for if both photoionization and shocks contribute to the line fluxes. The coexistence of shocks and AGN photoionization also provides an explanation for the most extreme objects in the NV1240/HeII1640 vs. NV1240/CIV1549 diagram without requiring largely super-solar metallicities. In addition, we show that there is a relationship between the [OII]3727/[OIII]5007 ratio and the radio size in radio galaxies. This strongly supports the hypothesis that the most compact (< 2kpc) and the largest (> 150 kpc) sources are dominated by photoionization, while intermediate-sized radio galaxies are dominated by shocks. We briefly discuss the possible origin of the relation between the shock-ionization balance and the radio size.

We analyze the stellar and nebular energy distributions of starbursts using our evolutionary synthesis model PEGASE coupled to the photoionization code CLOUDY. The originality of this study is to relate the evolution and the metallicity of the starburst to the past star formation history of the host galaxy.We compare our model predictions to an observed sample of ~750 starbursts. the most striking feature from line ratio fitting is the decreasing spread in ionization parameter U at high metallicity Z. The best fits of emission line ratios are obtained with a combination of a high- and a low-ionization components. No additional source of ionizing photons --shocks or hidden AGN -- is needed. Colors and equivalent widths are fitted in coherency with emission line ratios. An underlying population is needed, even for small-aperture observations. This evolved population not only reddens the continuum and dilutes the equivalent width of the emission lines, but also participates in the ionization process. Its main effect on line ratios is to maintain a high level of excitation when the burst stops. Models combining underlying populations typical of Hubble sequence galaxies and instantaneous starbursts with ages between 0 and 8 Myr agree satisfactorily with all the data.

We present results and source counts at 90micron extracted from the Preliminary Analysis of the European Large Area ISO Survey (ELAIS). The survey covered about 11.6 square degrees of the sky in four main areas and was carried out with the PHOT instrument onboard the Infrared Space Observatory (ISO). The survey is at least an order of magnitude deeper than the IRAS 100micron survey and is expected to provide constraints on the formation and evolution of galaxies. The majority of the detected sources are associated with galaxies on optical images. In some cases the optical associations are interacting pairs or small groups of galaxies suggesting the sample may include a significant fraction of luminous infrared galaxies. The source counts extracted from a reliable subset of the detected sources are in agreement with strongly evolving models of the starburst galaxy population.

We report the results of near-infrared VLT-ISAAC spectroscopy of a sample of 12 galaxies at z = 0.4-1.4, drawn from the ISOCAM survey of the Hubble Deep Field South. We find that the rest frame R-band spectra of the ISOCAM galaxies resemble those of powerful dust-enshrouded starbursts. Halpha emission is detected in 11 out of 12 objects down to a flux limit of 7x10^(-17) erg/cm^2/s, corresponding to a luminosity limit of 10^41 erg/s at z = 0.6, (for an Ho = 50 and Omega = 0.3 cosmology). From the Halpha luminosities in these galaxies we derive estimates of the star formation rate in the range 2--50 Mo/yr for stellar masses 1--100 Mo. The raw Halpha-based star formation rates are an order of magnitude or more lower than SFR(FIR) estimates based on ISOCAM LW3 fluxes. If the Halpha emission is corrected for extinction the median offset is reduced to a factor of 3. The sample galaxies are part of a new population of optically faint but infrared--luminous active starburst galaxies, which are characterized by an extremely high rate of evolution with redshift up to z~1.5 and expected to contribute significantly to the cosmic far-IR extragalactic background.

We present preliminary source counts at 6.7um and 15um from the Preliminary Analysis of the European Large Area ISO survey, with limiting flux densities of \~2mJy at 15um & ~1mJy at 6.7um. We separate the stellar contribution from the extragalactic using identifications with APM sources made with the likelihood ratio technique. We quantify the completeness & reliability of our source extraction using (a) repeated observations over small areas, (b) cross-IDs with stars of known spectral type, (c) detections of the PSF wings around bright sources, (d) comparison with independent algorithms. Flux calibration at 15um was performed using stellar IDs; the calibration does not agree with the pre-flight estimates, probably due to effects of detector hysteresis and photometric aperture correction. The 6.7um extragalactic counts are broadly reproduced in the Pearson & Rowan-Robinson model, but the Franceschini et al. (1997) model underpredicts the observed source density by ~0.5-1 dex, though the photometry at 6.7um is still preliminary. At 15um the extragalactic counts are in excellent agreement with the predictions of the Pearson & Rowan-Robinson (1996), Franceschini et al. (1994), Guiderdoni et al. (1997) and the evolving models of Xu et al. (1998), over 7 orders of magnitude in 15um flux density. The counts agree with other estimates from the ISOCAM instrument at overlapping flux densities (Elbaz et al. 1999), provided a consistent flux calibration is used. Luminosity evolution at a rate of (1+z)^3, incorporating mid-IR spectral features, provides a better fit to the 15um differential counts than (1+z)^4 density evolution. No-evolution models are excluded, and implying that below around 10mJy at 15um the source counts become dominated by an evolving cosmological population of dust-shrouded starbursts and/or active galaxies.

Deep galaxy counts are among the best constraints on the cosmic star formation history (SFH) of galaxies. The evolution of the star formation activity is followed on a wide range of redshifts (0< z < 4) covering most of the history of the Universe. Two incompatible interpretations of the observations are currently competing. After applying star formation rate (SFR) conversion factors to the CFRS, Halpha or ISO samples, many authors conclude to a strong increase (~ a factor 10) of the SFR from z=0 to z=1. They also find some evidence for a peak at z~1 and for a rapid decrease at higher redshifts. On the other side, the Hawaii deep surveys favor only a mild increase between z=0 and 1 (Cowie et al., 1996, 1999). To understand the reason for these discrepant interpretations, we consider three classes of galaxies: E/S0, Sa-Sb-Sbc, Sc-Sd-Im and bursting dwarfs. We use the new version of our evolutionary synthesis code, PEGASE which takes into account metallicity and dust effects. The main results are: i) Late-type galaxies contribute significantly to the local SFR, especially bursting dwarfs (Fioc and Rocca-Volmerange, 1999). Because of that, the cosmic SFR can not decrease by a factor 10 from z=0 to 1. This is in agreement with Cowie's result. ii) The SFR of intermediate-type galaxies has strongly decreased since z=1. Though the decrease is less than what find Lilly et al., 1996, this suggests that the CFRS and Halpha surveys are dominated by such bright early spirals. iii) The contribution of early-type galaxies increases rapidly from z=1 to their redshift of formation (> 2-3 for cosmological reasons). Their intense star formation rates at high-z give strong constraints on early ionization phases, primeval populations or metal enrichments.(abridged)

We provide here the documentation of the new version of the spectral evolution model PEGASE. PEGASE computes synthetic spectra of galaxies in the UV to near-IR range from 0 to 20 Gyr, for a given stellar IMF and evolutionary scenario (star formation law, infall, galactic winds). The radiation emitted by stars from the main sequence to the pre-supernova or white dwarf stage is calculated, as well as the extinction by dust. A simple modeling of the nebular emission (continuum and lines) is also proposed. PEGASE may be used to model starbursts as well as old galaxies. The main improvements of PEGASE.2 relative to PEGASE.1 (Fioc & Rocca-Volmerange 1997) are the following: (1)The stellar evolutionary tracks of the Padova group for metallicities between 0.0001 and 0.1 have been included; (2)The evolution of the metallicity of the interstellar medium (ISM) due to SNII, SNIa and AGB stars is followed. Stars are formed with the same metallicity as the ISM (instead of a solar metallicity in PEGASE.1), providing thus a metallicity-consistent model; (3)Lejeune et al.'s library of stellar spectra is used; (4)The extinction by dust is computed for geometries corresponding to disk and spheroidal galaxies using a radiative transfer code taking into account the scattering. The main outputs (as a function of time) are spectra, colors and magnitudes in various photometric systems, luminosities, type II and Ia supernovae rates, line intensities and equivalent widths, amount and metallicity of stars and gas, mass locked in stellar remnants, optical depth and total dust emission. The corresponding article (Fioc & Rocca-Volmerange 2000) will be submitted soon. A detailed modeling of the spectrum of the dust emission and of HII regions (Moy, Rocca-Volmerange & Fioc 2000) will be included in future versions.

We present UV spectra of a small sample of high redshift active galaxies: the hyperluminous, gravitationally lensed system SMM J02399-0136 (z=2.8) and three powerful radio galaxies (z~2.5). Extended continuum and emission lines are detected in all objects. The three radio galaxies present typical spectra with dominant Lya and weaker CIV1550, HeII1640 and CIII]1909. The spectrum of SMM J02399-0136 is strikingly different, showing relatively weak Lya and HeII and strong NV1240 relative to the C lines. We find that the weakness of HeII can be explained if the intermediate density narrow line region (n~10^{5-6} cm^-3 dominates the emission line spectrum, rather than the more extended low density gas (n<100 cm^-3). The line ratios of MRC2025-218 suggest that this is also the case in this object. MRC2025-218 and SMM J02399-0136 show the largest NV/HeII and NV/CIV values found in high redshift radio galaxies. The failure of solar abundance models to explain these line ratios and studies of metal abundances in high redshift quasars and radio galaxies, suggest that nitrogen is overabundant in the ionized gas of these objects. Alternative possibilities which we discuss include NV emission from the broad line region and differential amplification by a gravitational lens. We report the discovery of a PCygni profile in the CIV line in the spectrum of MRC2025-218. We detect also CII1334.5 in absorption and PCygni profiles for the lines SiIV1393.8,1402.8. We do not detect any purely photospheric stellar lines. The nature of the absorption features is not clear. It could be due to stars or an associated absorption line system.

The European Large Area ISO Survey (ELAIS) has surveyed 12 square degrees of the sky at 15 and 90 microns, and subsets of this area at 6.75 and 175 microns, using the Infrared Space Observatory (ISO). This project was the largest single open time programme executed by ISO, taking 375 hours of data. A preliminary catalogue of more than 1000 galaxies has been produced. In this talk we describe the goals of the project, describe the follow-up programmes that are in progress, and present some first scientific results including a provisional number count analysis at 15 and 90 microns.

Growth curves of the near-infrared (NIR) magnitude as a function of the aperture have been built and used to derive NIR total magnitudes from aperture data taken from the literature. By cross-correlating with optical and redshift data, absolute magnitudes and optical-to-NIR colors have been computed for some 1000 galaxies of different types. Significant color gradients are observed, underlining that small aperture colors may lead to a biased picture of the stellar populations of galaxies. A statistical analysis using various estimators taking into account the intrinsic scatter has been performed to establish relations between the colors, the morphological type, the inclination or the shape, and the intrinsic luminosity. The combination of the optical and the NIR should obviously improve our understanding of the evolution of galaxies. Despite the intrinsic scatter, especially among star-forming galaxies, optical-to-NIR colors show a very well defined sequence with type, blueing by 1.3 mag from ellipticals to irregulars. The colors of spiral galaxies strongly redden with increasing inclination and put new constraints on the modeling of the extinction. No such effect is observed for lenticular galaxies. We also find that rounder ellipticals tend to be redder. A color-absolute magnitude relation is observed inside each type, with a slope significantly steeper for early and intermediate spirals than for ellipticals or late spirals. This stresses the importance of considering both the mass and the type to study the star formation history of galaxies.

Galaxy counts from bright ultraviolet (UV) and deep optical spectroscopic surveys have revealed an unexpectedly large number of very blue galaxies. The colors and luminosities of these objects indicate that they are dwarf galaxies undergoing bursts of star formation. We use a galaxy evolution model (PEGASE, Fioc & Rocca-Volmerange 1997) to describe this population as galaxies undergoing cyclical bursts of star formation, thereby determining the luminosity function of these galaxies. When these bursting galaxies are added to normally evolving populations, the combination reproduces the UV number counts, color distributions and deep optical redshift distributions fairly well. Optical (including the Hubble Deep Field) and near-infrared number counts are fitted assuming an open or a flat, Lambda-dominated, Universe. The high amplitude of the angular correlation function of very blue galaxies discovered by Landy et al. (1996) is also recovered in this modelling. The number of bursting galaxies is only a small fraction of the total number of galaxies at optical and near-infrared wavelengths, even at faintest magnitudes. In our evolution modelling, normal galaxies explain most of the blue excess in a low-Omega Universe. The problem of the blue excess remains in a flat Universe without a cosmological constant.

A deep survey of the European Large Area ISO Survey (ELAIS) field in the southern celestial hemisphere (hereinafter S1) has been carried out with the Australia Telescope Compact Array (ATCA) at 1.4 GHz. The S1 region, covering about 4 square degrees, has been also surveyed in the mid- and far-infrared (5-200 micron) with the Infrared Space Observatory (ISO). The radio survey provides uniform coverage of the entire S1 region, with a sensitivity (5 sigma) of 0.4 mJy over the whole area and 0.2 mJy in the center. To this sensitivity, virtually all the radio counterparts of the far-infrared extragalactic ISO sources should be detected. This paper presents a radio sample - complete at the 5 sigma level - consisting of 581 sources detected at 1.4 GHz. Of these radio sources, 349 have peak flux density in the range 0.2-1 mJy, forming a new homogeneous sample of sub-mJy radio sources. Due to its size, depth and multi-waveband coverage, the sample will allow us to study in greater detail the sub-mJy radio source population. The full catalogue will be available from http://athena.ph.ac.uk/

We have used the Very Large Array(VLA) in C configuration to carry out a sensitive 20cm radio survey of regions of sky that have been surveyed in the Far Infra-Red over the wavelength range 5-200 microns with ISO as part of the European Large Area ISO Survey(ELAIS). As usual in surveys based on a relatively small number of overlapping VLA pointings the flux limit varies over the area surveyed. The survey has a flux limit that varies from a 5$\sigma$ limit of 0.135mJy over an area of 0.12deg$^2$ to a 5$\sigma$ limit of 1.15mJy or better over the whole region covered of 4.22 deg$^2$. In this paper we present the radio catalogue of 867 sources. These regions of sky have previously been surveyed to shallow flux limits at 20cm with the VLA as part of the VLA D configuration NVSS(FWHM=45 arcsec) and VLA B configuration FIRST(FWHM=5 arcsec) surveys. We have carried out a a detailed comparison of the reliability of our own survey and these two independent surveys in order to assess the reliability and completeness of each survey.

The far-ultraviolet (UV) counts and the deep optical spectroscopic surveys have revealed an unexpected number of very blue galaxies (vBG). Using constraints from the UV and optical, we apply the galaxy evolution model PEGASE (Fioc & Rocca-Volmerange 1997, hereafter FRV) to describe this population with a cycling star formation. When added to normally evolving galaxy populations, vBG are able to reproduce UV number counts and color distributions as well as deep optical redshift distributions fairly well. Good agreement is also obtained with optical counts (including the Hubble Deep Field). The number of vBG is only a small fraction of the number of normal galaxies, even at faintest magnitudes. In our modelling, the latter explain the bulk of the excess of faint blue galaxies in an open Universe. The problem of the blue excess remains in a flat Universe without cosmological constant.

A method of deriving and using merging history trees of dark matter galaxy haloes directly from pure gravity N-body simulations is presented. This combines the full non-linearity of N-body simulations with the flexibility of the semi-analytical approach. Merging history trees derived from power-law initial perturbation spectrum simulations (for indices $n=-2$ and $n=0$) by Warren et al. (1992) are shown. As an example of a galaxy formation model, these are combined with evolutionary stellar population synthesis, via simple scaling laws for star formation rates, showing that if most star formation occurs during merger-induced bursts, then a nearly flat faint-end slope of the galaxy luminosity function may be obtained in certain cases. Interesting properties of hierarchical halo formation are noted: (1) In a given model, merger rates may vary widely between individual haloes, and typically 20%-30% of a halo's mass may be due to infall of uncollapsed material. (2) Small mass haloes continue to form at recent times: as expected, the existence of young, low redshift, low metallicity galaxies (e.g., Izotov et al. 1997) is consistent with hierarchical galaxy formation models. (3) For $n=-2,$ the halo spatial correlation function can have a very high initial bias due to the high power on large scales.

PEGASE is a new spectrophotometric evolution model for starbursts and evolved galaxies of the Hubble sequence. Its main originality is the extension to the near-infrared (NIR) of the atlas of synthetic spectra of Rocca-Volmerange & Guiderdoni (1988) with a revised stellar library including cold star parameters and stellar tracks extended to the TP-AGB and the post-AGB phase. The NIR is coherently linked to the visible and the ultraviolet, so that the model is continuous on an exceptionally large wavelength range from 220 A up to 5 microns. Moreover, a precise algorithm allows to follow very rapid evolutionary phases such as red supergiants or AGB crucial in the NIR. The nebular component is also computed in the NIR. The extinction correction is gas-dependent for spirals and ellipticals. A set of reference synthetic spectra at z=0, to which apply cosmological k- and evolution e- corrections for high-redshift galaxies, is built from fits of observational templates. Because of the lack of visible to NIR spectral templates for each Hubble type, we adopt statistical samples of colors, not fitted by previous models. A first application of this continuous model is to solve the problem of the slope of the bright galaxy counts from B=15 to 19 and of the normalization parameter of the Schechter luminosity function. Code sources, input and output data are available by anonymous ftp or at the WEB address of the authors.