Julien Devriendt

Abstract:

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. © 2009 ESO.

Abstract:

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

Abstract:

In a recent paper we presented the first semi-analytic model of galaxy formation in which the Thermally-Pulsing Asymptotic Giant Branch phase of stellar evolution has been fully implemented. Here we address the comparison with observations, and show how the TP-AGB recipe affects the performance of the model in reproducing the colours and near-IR luminosities of high-redshift galaxies. We find that the semi-analytic model with the TP-AGB better matches the colour-magnitude and colour-colour relations at z ~ 2, both for nearly-passive and for star-forming galaxies. The model with TP-AGB produces star-forming galaxies with red V-K colours, thus revising the unique interpretation of high-redshift red objects as 'red & dead'. We also show that without the TP-AGB the semi-analytic model fails at reproducing the observed colours, a situation that cannot be corrected by dust reddening. We also explore the effect of nebular emission on the predicted colour-magnitude relation of star-forming galaxies, to conclude that it does not play a significant role in reddening their colours, at least in the range of star-formation rates covered by the model. Finally, the rest-frame K-band luminosity function at z ~ 2.5 is more luminous by almost 1 magnitude. This indicates that the AGN feedback recipe that is adopted to regulate the high-mass end of the luminosity function should be sophisticated to take the effect of the stellar populations into account at high redshifts.