Building merger trees from cosmological N-body simulations
Astronomy and Astrophysics 506:2 (2009) 647-660
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.Influence of AGN jets on the magnetized ICM
ArXiv 0905.3345 (2009)
Abstract:
Galaxy clusters are the largest structures for which there is observational evidence of a magnetised medium. Central cores seem to host strong magnetic fields ranging from a few 0.1 microG up to several 10 microG in cooling flow clusters. Numerous clusters harbor central powerful AGN which are thought to prevent cooling flows in some clusters. The influence of such feedback on the magnetic field remains unclear: does the AGN-induced turbulence compensate the loss of magnetic amplification within a cool core? And how is this turbulence sustained over several Gyr? Using high resolution magneto-hydrodynamical simulations of the self-regulation of a radiative cooling cluster, we study for the first time the evolution of the magnetic field within the central core in the presence of a powerful AGN jet. It appears that the jet-induced turbulence strongly amplifies the magnetic amplitude in the core beyond the degree to which it would be amplified by pure compression in the gravitational field of the cluster. The AGN produces a non-cooling core and increases the magnetic field amplitude in good agreement with microG field observations.Building Merger Trees from Cosmological N-body Simulations
ArXiv 0902.0679 (2009)
Abstract:
Although a fair amount of work has been devoted to growing Monte-Carlo merger trees which 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. The first step to 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 extend to most if not all (sub)structure finders. We then highlight different ways to build merger histories from AdaptaHOP haloes and subhaloes, contrasting their various advantages and drawbacks. 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 which conducts a straightforward sequential treatment of these two steps. This is rooted in the complexity of the merging trees which 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.Cooling, gravity, and geometry: Flow-driven massive core formation
Astrophysical Journal 674:1 (2008) 316-328
Abstract:
We study numerically the formation of molecular clouds in large-scale colliding flows including self-gravity. The models emphasize the competition between the effects of gravity on global and local scales in an isolated cloud. Global gravity builds up large-scale filaments, while local gravity, triggered by a combination of strong thermal and dynamical instabilities, causes cores to form. The dynamical instabilities give rise to a local focusing of the colliding flows, facilitating the rapid formation of massive protostellar cores of a few hundred M⊙. The forming clouds do not reach an equilibrium state, although the motions within the clouds appear to be comparable to virial. The self-similar core mass distributions derived from models with and without self-gravity indicate that the core mass distribution is set very early on during the cloud formation process, predominantly by a combination of thermal and dynamical instabilities rather than by self-gravity. © 2008. The American Astronomical Society. All rights reserved.Cooling, Gravity and Geometry: Flow-driven Massive Core Formation
ArXiv 0709.2451 (2007)