Search

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

Number Counts of galaxies in two GALEX bands (FUV: 1530A and NUV: 2310A, both in AB magnitudes) are reported. They provide for the first time in the literature homogeneously calibrated number counts of UV galaxies covering continuously a very wide range of UV magnitude (14 -- 23.8). Both the FUV and NUV counts are inconsistent with a non-evolution model, while they are in good agreement with evolution models (essentially luminosity evolution) derived from the high-z UV luminosity functions of Arnouts et al. (2004). It is found that the contribution from galaxies detected by GALEX to the UV background is 0.68+-0.10 nW m-2 sr-1 at 1530A and 0.99+-0.15 nW m-2 sr-1 at 2310A. These are 66+-9% and 44+-6% of the total contributions of galaxies to the the UV background at 1530A, respectively, as estimated using the evolution models. ...

Non-resonant FeII* 2365, 2396, 2612, 2626 emission can potentially trace galactic winds in emission and provide useful constraints to wind models. From the 3'x3' mosaic of the Hubble Ultra Deep Field (UDF) obtained with the VLT/MUSE integral field spectrograph, we identify a statistical sample of 40 FeII* emitters and 50 MgII 2796, 2803 emitters from a sample of 271 [OII] 3726, 3729 emitters with reliable redshifts from z = 0.85 - 1.5 down to 2E-18 (3 sigma) ergs/s/cm^2 (for [OII]), covering the stellar mass range 10^8 - 10^11 Msun. The FeII* and MgII emitters follow the galaxy main sequence, but with a clear dichotomy. Galaxies with masses below 10^9 Msun and star formation rates (SFRs) of 10 Msun/year have FeII* emission without accompanying MgII emission. Between these two regimes, galaxies have both MgII and FeII* emission, typically with MgII P-Cygni profiles. Indeed, the MgII profile shows a progression along the main sequence from pure emission to P-Cygni profiles to strong absorption, due to resonant trapping. Combining the deep MUSE data with HST ancillary information, we find that galaxies with pure MgII emission profiles have lower star formation rate surface densities than those with either MgII P-Cygni profiles or FeII* emission. These spectral signatures produced through continuum scattering and fluorescence, MgII P-Cygni profiles and FeII* emission, are better candidates for tracing galactic outflows than pure MgII emission, which may originate from HII regions. We compare the absorption and emission rest-frame equivalent widths for pairs of FeII transitions to predictions from outflow models and find that the observations consistently have less total re-emission than absorption, suggesting either dust extinction or non-isotropic outflow geometries., 16 pages, 12 figures, 3 tables, accepted for publication in A&A (14 October 2017)

We present rest-frame Lyα equivalent widths (E W 0 ) of 417 Lyα emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at 2.9 Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured E W 0 values taking into account extended Lyα emission and UV continuum slopes (β ). Our LAEs reach unprecedented depths, both in Lyα luminosities and UV absolute magnitudes, from log (L Lyα /erg s-1 ) ~ 41.0 to 43.0 and from M UV ~ −16 to −21 (0.01−1.0 L * z =3 ). The E W 0 values span the range of ~ 5 to 240 A or larger, and their distribution can be well fitted by an exponential law N = N 0 exp(−E W 0 /w 0 ). Owing to the high dynamic range in M UV , we find that the scale factor, w 0 , depends on M UV in the sense that including fainter M UV objects increases w 0 , i.e., the Ando effect. The results indicate that selection functions affect the E W 0 scale factor. Taking these effects into account, we find that our w 0 values are consistent with those in the literature within 1σ uncertainties at 2.9 UV and L Lyα . Interestingly, we find 12 objects with E W 0 > 200 A above 1σ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak Civλ 1549 emission line. For the remaining 10 very large E W 0 LAEs, we find that the E W 0 values can be reproduced by young stellar ages ( ⊙ ). Otherwise, at least part of the Lyα emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

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

[abridged] We describe LyMAS (Ly-alpha Mass Association Scheme), a method of predicting clustering statistics in the Ly-alpha forest on large scales from moderate resolution simulations of the dark matter distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the "Horizon MareNostrum" simulation, a 50 Mpc/h comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F_s|delta_s) of the transmitted flux F_s, smoothed (1-dimensionally) over the spectral resolution scale, on the dark matter density contrast delta_s, smoothed (3-dimensionally) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III BOSS at z=2.5, and we find optimal results for a dark matter smoothing length sigma=0.3 Mpc/h (comoving). In extended form, LyMAS exactly reproduces both the 1-dimensional power spectrum and 1-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum dark matter field, LyMAS accurately predicts the 2-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 Mpc/h, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ("Fluctuating Gunn-Peterson Approximation"), often used for large volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024^3 N-body simulations of a 300 Mpc/h and 1.0 Gpc/h cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3-d Ly-alpha forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, dark matter, space geometry, and IGM physics., Accepted for publication in ApJ (minor corrections from the previous version). Catalogs of mock BOSS spectra and relevant data can be found at: http://www2.iap.fr/users/peirani/lymas/lymas.htm

We present a new implementation of radiation hydrodynamics (RHD) in the adaptive mesh refinement (AMR) code RAMSES. The multi-group radiative transfer (RT) is performed on the AMR grid with a first-order Godunov method using the M1 closure for the Eddington tensor, and is coupled to the hydrodynamics via non-equilibrium thermochemistry of hydrogen and helium. This moment-based approach has the large advantage that the computational cost is independent of the number of radiative sources - it can even deal with continuous regions of emission such as bound-free emission from gas. As it is built directly into RAMSES, the RT takes natural advantage of the refinement and parallelization strategies already in place. Since we use an explicit advection solver for the radiative transport, the time step is restricted by the speed of light - a severe limitation that can be alleviated using the so--called "reduced speed of light" approximation. We propose a rigorous framework to assess the validity of this approximation in various conditions encountered in cosmology and galaxy formation. We finally perform with our newly developed code a complete suite of RHD tests, comparing our results to other RHD codes. The tests demonstrate that our code performs very well and is ideally suited for exploring the effect of radiation on current scenarios of structure and galaxy formation., Comment: 48 pages, 37 figures, accepted for publication in MNRAS

Simulations of galaxy evolution aim to capture our current understanding as well as to make predictions for testing by future experiments. Simulations and observations are often compared in an indirect fashion: physical quantities are estimated from the data and compared to models. However, many applications can benefit from a more direct approach, where the observing process is also simulated and the models are seen fully from the observer's perspective. To facilitate this, we have developed the Millennium Run Observatory (MRObs), a theoretical virtual observatory which uses virtual telescopes to `observe' semi-analytic galaxy formation models based on the suite of Millennium Run dark matter simulations. The MRObs produces data that can be processed and analyzed using the standard software packages developed for real observations. At present, we produce images in forty filters from the rest-frame UV to IR for two stellar population synthesis models, three different models of IGM absorption, and two cosmologies (WMAP1/7). Galaxy distributions for a large number of mock lightcones can be `observed' using models of major ground- and space-based telescopes. The data include lightcone catalogues linked to structural properties of galaxies, pre-observation model images, mock telescope images, and Source Extractor products that can all be traced back to the higher level dark matter, semi-analytic galaxy, and lightcone catalogues available in the Millennium database. Here, we describe our methods and announce a first public release of simulated surveys (e.g., SDSS, CFHT-LS, GOODS, GOODS/ERS, CANDELS, and HUDF). The MRObs browser, an online tool, further facilitates exploration of the simulated data. We demonstrate the benefits of a direct approach through a number of example applications (galaxy number counts in CANDELS, clusters, morphologies, and dropout selections)., MNRAS, in press. Millennium Run Observatory data products, online tools, and more available through http://galformod.mpa-garching.mpg.de/mrobs/

We develop a Principal Component Analysis aimed at classifying a sub-set of 27,350 spectra of galaxies in the range 0.4 < z < 1.0 collected by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We apply an iterative algorithm to simultaneously repair parts of spectra affected by noise and/or sky residuals, and reconstruct gaps due to rest-frame transformation, and obtain a set of orthogonal spectral templates that span the diversity of galaxy types. By taking the three most significant components, we find that we can describe the whole sample without contamination from noise. We produce a catalogue of eigen-coefficients and template spectra that will be part of future VIPERS data releases. Our templates effectively condense the spectral information into two coefficients that can be related to the age and star formation rate of the galaxies. We examine the spectrophotometric types in this space and identify early, intermediate, late and starburst galaxies., 15 pages, 20 images, accepted for publication in MNRAS: MN-12-1739-MJ.R1

[Abridged] We perform on galaxy mock catalogues the same colour-density analysis made by Cucciati et al. (2006) on a 5 Mpc/h scale using the VVDS-Deep survey, and compare the results from mocks with observed data. We use mocks with the same flux limits (I=24) as the VVDS (CMOCKS), built using the semi- analytic model by De Lucia & Blaizot (2007) applied to the Millennium Simulation. From CMOCKS, we extracted samples of galaxies mimicking the VVDS observational strategy (OMOCKS). We computed the B-band Luminosity Function LF and the colour-density relation (CDR) in the mocks. We find that the LF in mocks roughly agrees with the observed LF, but at z, Comment: 15 pages, 12 figures, accepted for publication in A&A

Lyman-$\alpha$ nebulae are usually found in massive environments at high redshift ($z > 2$). The origin of their Lyman-$\alpha$ (Lya) emission remains debated. Recent polarimetric observations showed that at least some Lya sources are polarized. This is often interpreted as a proof that the photons are centrally produced, and opposed to the scenario in which the Lya emission is the cooling radiation emitted by gas heated during the accretion onto the halo. We suggest that this scenario is not incompatible with the polarimetric observations. In order to test this idea, we post-process a radiative hydrodynamics simulation of a blob with the MCLya Monte Carlo transfer code. We compute radial profiles for the surface brightness and the degree of polarization and compare them to existing observations. We find that both are consistent with a significant contribution of the extragalactic gas to the Lya emission. Most of the photons are centrally emitted and scattered inside the filament afterwards, producing the observed high level of polarization. We argue that the contribution of the extragalactic gas to the Lya emission does not prevent polarization to arise. On the contrary, we find that pure galactic emission causes the polarization profile to be too steep to be consistent with observations., Comment: 14 pages, 9 figures (+4 in appendices). Accepted for publication in Astronomy & Astrophysics

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

[Abridged] We present a homogeneous and complete catalogue of optical groups identified in the purely flux limited (17.5, Submitted to A&A, revised version after referee comments, Table 5 added

The Multi Unit Spectroscopic Explorer (MUSE) is a second-generation VLT panoramic integral-field spectrograph currently in manufacturing, assembly and integration phase. MUSE has a field of 1x1 arcmin2 sampled at 0.2x0.2 arcsec2 and is assisted by the VLT ground layer adaptive optics ESO facility using four laser guide stars. The instrument is a large assembly of 24 identical high performance integral field units, each one composed of an advanced image slicer, a spectrograph and a 4kx4k detector. In this paper we review the progress of the manufacturing and report the performance achieved with the first integral field unit., Comment: 10 pages, 5 figures

(Abridged) We assemble a sample of 3258 low-redshift galaxies from the SDSS DR6 with complementary photometric observations by GALEX, 2MASS and IRAS at far-ultraviolet and infrared wavelengths. We use a recent, simple but physically motivated model to interpret the observed spectral energy distributions of the galaxies in this sample in terms of statistical constraints on physical parameters describing the star formation history and dust content. The focus on a subsample of 1658 galaxies with highest S/N observations enables us to investigate most clearly several strong correlations between various derived physical properties of galaxies. We find that the typical dust mass of a star-forming correlates remarkably well with the star formation rate (SFR). We also find that the dust-to-stellar mass ratio, the ratio of dust mass to star formation rate and the fraction of dust luminosity contributed by the diffuse interstellar medium all correlate strongly with specific SFR. A comparison with recent models of chemical and dust evolution of galaxies suggests that these correlations could arise, at least in part, from an evolutionary sequence. As galaxies form stars, their ISM becomes enriched in dust, while the drop in gas supply makes the specific SFR decrease. Interestingly, as a result, a young, actively star-forming galaxy with low dust-to-gas ratio may still be highly dusty because it contains large amounts of interstellar gas. This may be important for the interpretation of the infrared emission from young, gas-rich star-forming galaxies at high redshift. Our study provides a useful local reference for future statistical studies of the star formation and dust properties of galaxies at high redshifts., Comment: 15 pages, 10 figures, accepted for publication in MNRAS. Full-resolution figures available from http://users.physics.uoc.gr/~dacunha/paper_dacunha.pdf

Several observational works have attempted to isolate the effects of galaxy interactions by comparing galaxies in pairs with isolated galaxies. However, different authors have proposed different ways to build these so-called control samples (CS). By using mock galaxy catalogues of the Sloan Digital Sky Survey Data Release 4 buildup from the Millennium Simulation, we explore how the way of building a CS might introduce biases which could affect the interpretation of results. We make use of the fact that the physics of interactions are not included in the semi-analytic model, to infer that any difference between the mock control and pair samples can be ascribed to selection biases. Thus, we find that galaxies in pairs artificially tend to be older and more bulge dominated, and to have less cold gas and different metallicities than their isolated counterparts. Also because of a biased selection, galaxies in pairs tend to live in higher density environments and in haloes of larger masses. We find that imposing constraints on redshift, stellar masses and local densities diminishes the selection biases by ≈70 per cent. Based on these findings, we suggest observers how to build a unique and unbiased CS in order to reveal the effect of galaxy interactions., Facultad de Ciencias Astronómicas y Geofísicas

International audience; Context: In the past decade or so, using numerical N-body simulations to describe the gravitational clustering of dark matter (DM) in an expanding universe has become the tool of choice for tackling the issue of hierarchical galaxy formation. As mass resolution increases with the power of supercomputers, one is able to grasp finer and finer details of this process, resolving more and more of the inner structure of collapsed objects. This begs one to revisit time and again the post-processing tools with which one transforms particles into “invisible” dark matter haloes and from thereon into luminous galaxies. Aims: Although a fair amount of work has been devoted to growing Monte-Carlo merger trees that resemble those built from an N-body simulation, comparatively little effort has been invested in quantifying the caveats one necessarily encounters when one extracts trees directly from such a simulation. To somewhat revert the tide, this paper seeks to provide its reader with a comprehensive study of the problems one faces when following this route. Methods: The first step in building merger histories of dark matter haloes and their subhaloes is to identify these structures in each of the time outputs (snapshots) produced by the simulation. Even though we discuss a particular implementation of such an algorithm (called AdaptaHOP) in this paper, we believe that our results do not depend on the exact details of the implementation but instead extend to most if not all (sub)structure finders. To illustrate this point in the appendix we compare AdaptaHOP's results to the standard friend-of-friend (FOF) algorithm, widely utilised in the astrophysical community. We then highlight different ways of building merger histories from AdaptaHOP haloes and subhaloes, contrasting their various advantages and drawbacks. Results: We find that the best approach to (sub)halo merging histories is through an analysis that goes back and forth between identification and tree building rather than one that conducts a straightforward sequential treatment of these two steps. This is rooted in the complexity of the merging trees that have to depict an inherently dynamical process from the partial temporal information contained in the collection of instantaneous snapshots available from the N-body simulation. However, we also propose a simpler sequential “Most massive Substructure Method” (MSM) whose trees approximate those obtained via the more complicated non sequential method. Appendices are only available in electronic form at: http://www.aanda.org

We present a new method to estimate the average star formation rate per unit stellar mass (SSFR) of a stacked population of galaxies. We combine the spectra of 600-1000 galaxies with similar stellar masses and parameterise the star formation history of this stacked population using a set of exponentially declining functions. The strength of the Hydrogen Balmer absorption line series in the rest-frame wavelength range 3750-4150\AA is used to constrain the SSFR by comparing with a library of models generated using the BC03 stellar population code. Our method, based on a principal component analysis (PCA), can be applied in a consistent way to spectra drawn from local galaxy surveys and from surveys at $z \sim 1$, and is only weakly influenced by attenuation due to dust. We apply our method to galaxy samples drawn from SDSS and DEEP2 to study mass-dependent growth of galaxies from $z \sim 1$ to $z \sim 0$. We find that, (1) high mass galaxies have lower SSFRs than low mass galaxies; (2) the average SSFR has decreased from $z=1$ to $z=0$ by a factor of $\sim 3-4$, independent of galaxy mass. Additionally, at $z \sim 1$ our average SSFRs are a factor of $2-2.5$ lower than those derived from multi-wavelength photometry using similar datasets. We then compute the average time (in units of the Hubble time, $t_{\rm H}(z)$) needed by galaxies of a given mass to form their stars at their current rate. At both $z=0$ and at $z=1$, this timescale decreases strongly with stellar mass from values close to unity for galaxies with masses $\sim 10^{10} M_{\odot}$, to more than ten for galaxies more massive than $ 10^{11} M_{\odot}$. Our results are in good agreement with models in which AGN feedback is more efficient at preventing gas from cooling and forming stars in high mass galaxies., Comment: 14 pages, 11 figures, accepted by MNRAS

We analyze the clustering properties of ultraviolet selected galaxies by using GALEX-SDSS data at z 2) to low redshift (b ~ 0.79^{+0.1}_{-0.08}). When accounting for the fraction of the star formation activity enclosed in the different samples, our results suggest that the bulk of star formation migrated from high mass dark matter halos at z>2 (10^12 < M_min < 10^13 M_sun, located in high density regions), to less massive halos at low redshift (M_min < 10^12 M_sun, located in low density regions). This result extends the ``downsizing'' picture (shift of the star formation activity from high stellar mass systems at high z to low stellar mass at low z) to the dark matter distribution., Accepted for Publication in the Special GALEX Ap. J. Supplement, December 2007 Version with full resolution fig1 available at http://taltos.pha.jhu.edu/~sebastien/papers/Galex_p2.ps.gz

We present the first measurements of the angular correlation function of galaxies selected in the far (1530 A) and near (2310 A) Ultraviolet from the GALEX survey fields overlapping SDSS DR5 in low galactic extinction regions. The area used covers 120 sqdeg (GALEX - MIS) down to magnitude AB = 22, yielding a total of 100,000 galaxies. The mean correlation length is ~ 3.7 \pm 0.6 Mpc and no significant trend is seen for this value as a function of the limiting apparent magnitude or between the GALEX bands. This estimate is close to that found from samples of blue galaxies in the local universe selected in the visible, and similar to that derived at z ~ 3 for LBGs with similar rest frame selection criteria. This result supports models that predict anti-biasing of star forming galaxies at low redshift, and brings an additional clue to the downsizing of star formation at z, Accepted for publication in GALEX Special ApJs, December 2007

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

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

We have developed a method for detecting clusters in large imaging surveys, based on the detection of structures in galaxy density maps made in slices of photometric redshifts. This method was first applied to the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Deep 1 field by Mazure et al. (2007), then to all the Deep and Wide CFHTLS fields available in the T0004 data release by Adami et al. (2010). The validity of the cluster detection rate was estimated by applying the same procedure to galaxies from the Millennium simulation. Here we analyse with the same method the full CFHTLS Wide survey, based on the T0006 data release. In a total area of 154 deg2, we have detected 4061 candidate clusters at 3sigma or above (6802 at 2sigma and above), in the redshift range 0.1, Accepted for publication in Astronomy \& Astrophysics

In the last decade we have seen an enormous increase in the size and quality of spectroscopic galaxy surveys, both at low and high redshift. New statistical techniques to analyse large portions of galaxy spectra are now finding favour over traditional index based methods. Here we will review a new robust and iterative Principal Component Analysis (PCA) algorithm, which solves several common issues with classic PCA. Application to the 4000AA break region of galaxies in the VIMOS VLT Deep Survey (VVDS) and Sloan Digital Sky Survey (SDSS) gives new high signal-to-noise ratio spectral indices easily interpretable in terms of recent star formation history. In particular, we identify a sample of post-starburst galaxies at z~0.7 and z~0.07. We quantify for the first time the importance of post-starburst galaxies, consistent with being descendants of gas-rich major mergers, for building the red sequence. Finally, we present a comparison with new low and high redshift "mock spectroscopic surveys" derived from a Millennium Run semi-analytic model., Comment: 7 pages, 3 figures. Conference proceedings in "Classification and Discovery in Large Astronomical Surveys", 2008, C.A.L. Bailer-Jones (ed.)