Adrianne Slyz

Dancing in the dark: galactic properties trace spin swings along the cosmic web

A large-scale hydrodynamical cosmological simulation, Horizon-AGN , is used to investigate the alignment between the spin of galaxies and the large-scale cosmic filaments above redshift one. The analysis of more than 150 000 galaxies with morphological diversity in a 100h −1 Mpc comoving box size shows that the spin of low-mass, rotationdominated, blue, star-forming galaxies is preferentially aligned with their neighbouring filaments. High-mass, dispersion-dominated, red, quiescent galaxies tend to have a spin perpendicular to nearby filaments. The reorientation of the spin of massive galaxies is provided by galaxy mergers which are significant in the mass build up of high-mass galaxies. We find that the stellar mass transition from alignment to misalignment happens around 3×10 10 M⊙. This is consistent with earlier findings of a dark matter mass transition for the orientation of the spin of halos (5 × 10 11 M⊙ at the same redshift from Codis et al. 2012). With these numerical evidence, we advocate a scenario in which galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the physical and morphological properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as Active Galactic Nuclei is mandatory to quench in situ star formation in massive galaxies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping galaxy spins misaligned with cosmic filaments. It also promotes diversity amongst galaxy properties.

Heating cooling flows with jets

Active galactic nuclei are clearly heating gas in ‘cooling flows’. The effectiveness and spatial distribution of the heating are controversial. We use three-dimensional simulations on adaptive grids to study the impact on a cooling flow of weak, subrelativistic jets. The simulations show cavities and vortex rings as in the observations. The cavities are fast-expanding dynamical objects rather than buoyant bubbles as previously modelled, but shocks still remain extremely hard to detect with X-rays. At late times the cavities turn into overdensities that strongly excite the g modes of a cluster. These modes damp on a long time-scale. Radial mixing is shown to be an important phenomenon, but the jets weaken the metallicity gradient only very near the centre. The central entropy density is modestly increased by the jets. We use a novel algorithm to impose the jets on the simulations.

The Birth of Molecular Clouds: Formation of Atomic Precursors in Colliding Flows

Molecular cloud complexes (MCCs) are highly structured and turbulent. Observational evidence suggests that MCCs are dynamically dominated systems, rather than quasi-equilibrium entities. The observed structure is more likely a consequence of the formation process than something that is imprinted after the formation of the MCC. Converging flows provide a natural mechanism to generate MCC structure. We present a detailed numerical analysis of this scenario. Our study addresses the evolution of an MCC from its birth in colliding atomic hydrogen flows up until the point when H2 may begin to form. A combination of dynamical and thermal instabilities breaks up coherent flows efficiently, seeding the small-scale nonlinear density perturbations necessary for local gravitational collapse and thus allowing (close to) instantaneous star formation. Many observed properties of MCCs come as a natural consequence of this formation scenario. Since converging flows are omnipresent in the ISM, we discuss the general applicability of this mechanism, from local star formation regions to galaxy mergers.

Cooling, Gravity, and Geometry: Flow-driven Massive Core Formation

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.

Connecting the cosmic web to the spin of dark haloes: implications for galaxy formation

We investigate the alignment of the spin of dark matter halos relative (i) to the surrounding large-scale filamentary structure, and (ii) to the tidal tensor eigenvectors using the Horizon 4π dark matter simulation which resolves over 43 million dark matter halos at redshift zero. We detect a clear mass transition: the spin of dark matter halos above a critical massM s 0 ≈ 5(±1) �10 12 M⊙ tends to be perpendicular to the closest large scale filament (with an excess probability up to 12%), and aligned with the intermediate axis of the tidal tensor (with an excess probability of up to 40%), whereas the spin of low-mass halos is more likely to be aligned with the closest filament (with an excess probability up to 15%). Furthermore, this critical mass is redshift-dependent, scaling as M s (z) ≈ M s �(1 +z) −γs with γs = 2.5 ± 0.2. A similar fit for the redshift evolution of the tidal tensor transition mass yields M t ≈ 8(±2) �10 12 M⊙ and γt = 3 ± 0.3. This critical mass also varies weakly with the scale defining filaments. We propose an interpretation of this signal in terms of large-scale cosmic flows. In this picture,most low-mass halos areformed through the winding offlows embedded in misaligned walls; hence they acquire a spin parallel to the axis of the resulting filaments forming at the intersection of these walls. On the other hand, more massive halos are typically the products of later mergers along such filaments, and thus they acquire a spin perpendicular to this direction when their orbital angular momentum is converted into spin. We show that this scenario is consistent with both the measured excess probabilities of alignment w.r.t. the eigen-directions of the tidal tensor, and halo merger histories. On a more qualitative level, it also seems compatible with 3D visualization of the structure of the cosmic web as traced by “smoothed” dark matter simulations or gas tracer particles. Finally, it provides extra support to the disc forming paradigm presented by Pichon et al. (2011) as it extends it by characterizing the geometry of secondary infall at high redshift.

Towards simulating star formation in the interstellar medium

As a first step to a more complete understanding of the local physical processes which determine star formation rates (SFRs) in the interstellar medium (ISM), we have performed controlled numerical experiments consisting of hydrodynamical simulations of a kilo-parsec scale, periodic, highly supersonic and ”turbulent” three-dimensional flow. Using simple but physically motivated recipes for identifying star forming regions, we convert gas into stars which we follow self-consistently as they impact their surroundings through supernovae and stellar winds. We investigate how various processes (turbulence, radiative cooling, self-gravity, and supernovae feedback) structure the ISM, determine its energetics, and consequently affect its SFR. We find that the one-point statistical measurement captured by the probability density function (PDF) is sensitive to the simulated physics. The PDF is consistent with a log-normal distribution for the runs which remove gas for star formation and have radiative cooling, but implement neither supernovae feedback nor self-gravity. In this case, the dispersion, �s, of the log-normal decays with time and scales with p ln(1 + (Mrms/2) 2 ) where Mrms is the root-mean-squared Mach number of the simulation volume, s = ln � , and � is the gas density. With the addition of self-gravity, the log-normal consistently under-predicts the high density end of the PDF which approaches a power law. With supernovae feedback, regardless of whether we consider self-gravity or not, the PDF becomes markedly bimodal with most of the simulation volume occupied by low density gas. Aside from its effect on the density structure of the medium, including self-gravity and/or supernovae feedback changes the dynamics of the medium by halting the decay of the kinetic energy. Since we find that the SFR depends most strongly on the underlying velocity field, the SFR declines in the runs lacking a means to sustain the kinetic energy, and the subsequent high density constrasts. This strong dependence on the gas velocity dispersion is in agreement with Silk’s formula for the SFR (Silk 2001) which also takes the hot gas porosity, and the average gas density as important parameters. Measuring the porosity of the hot gas for the runs with supernovae feedback, we compare Silk’s model for the SFR to our measured SFR and find agreement to better than a factor two.

Self-regulated growth of supermassive black holes by a dual jet/heating AGN feedback mechanism: methods, tests and implications for cosmological simulations

We develop a new sub-grid model for the growth of supermassive Black Holes (BHs) and their associated Active Galactic Nuclei (AGN) feedback in hydrodynamical cosmological simulations. Assuming that BHs are created in the early stages of galaxy formation, they grow by mergers and accretion of gas at a Eddington-limited Bondi accretion rate. However this growth is regulated by AGN feedback which we model using two different modes: a quasar-heating mode when accretion rates onto the BHs are comparable to the Eddington rate, and a radio-jet mode at lower accretion rates. In other words, our feedback model deposits energy as a succession of thermal bursts and jet outflows depending on the properties of the gas surrounding the BHs. We assess the plausibility of such a model by comparing our results to observational measurements of the coevolution of BHs and their host galaxy properties, and check their robustness with respect to numerical resolution. We show that AGN feedback must be a crucial physical ingredient for the formation of massive galaxies as it appears to be the only physical mechanism able to efficiently prevent the accumulation of and/or expel cold gas out of halos/galaxies and significantly suppress star formation. Our model predicts that the relationship between BHs and their host galaxy mass evolves as a function of redshift, because of the vigorous accretion of cold material in the early Universe that drives Eddington-limited accretion onto BHs. Quasar activity is also enhanced at high redshift. However, as structures grow in mass and lose their cold material through star formation and efficient BH feedback ejection, the AGN activity in the low-redshift Universe becomes more and more dominated by the radio mode, which powers jets through the hot circum-galactic medium.

Jet-regulated cooling catastrophe

We present the first implementation of active galactic nuclei (AGN) feedback in the form of momentum-driven jets in an adaptive mesh refinement (AMR) cosmological resimulation of a galaxy cluster. The jets are powered by gas accretion on to supermassive black holes (SMBHs) which also grow by mergers. Throughout its formation, the cluster experiences different dynamical states: both a morphologically perturbed epoch at early times and a relaxed state at late times allowing us to study the different modes of black hole (BH) growth and associated AGN jet feedback. BHs accrete gas efficiently at high redshift (z > 2), significantly pre-heating proto-cluster haloes. Gas-rich mergers at high redshift also fuel strong, episodic jet activity, which transports gas from the proto-cluster core to its outer regions. At later times, while the cluster relaxes, the supply of cold gas on to the BHs is reduced leading to lower jet activity. Although the cluster is still heated by this activity as sound waves propagate from the core to the virial radius, the jets inefficiently redistribute gas outwards and a small cooling flow develops, along with low-pressure cavities similar to those detected in X-ray observations. Overall, our jet implementation of AGN feedback quenches star formation quite efficiently, reducing the stellar content of the central cluster galaxy by a factor of 3 compared to the no-AGN case. It also dramatically alters the shape of the gas density profile, bringing it in close agreement with the β model favoured by observations, producing quite an isothermal galaxy cluster for gigayears in the process. However, it still falls short in matching the lower than universal baryon fractions which seem to be commonplace in observed galaxy clusters.

Rigging dark haloes: why is hierarchical galaxy formation consistent with the inside-out build-up of thin discs?

State-of-the-art hydrodynamical simulations show that gas inflow through the virial sphere of dark matter haloes is focused (i.e. has a preferred inflow direction), consistent (i.e. its orientation is steady in time) and amplified (i.e. the amplitude of its advected specific angular momentum increases with time). We explain this to be a consequence of the dynamics of the cosmic web within the neighbourhood of the halo, which produces steady, angular momentum rich, filamentary inflow of cold gas. On large scales, the dynamics within neighbouring patches drives matter out of the surrounding voids, into walls and filaments before it finally gets accreted on to virialized dark matter haloes. As these walls/filaments constitute the boundaries of asymmetric voids, they acquire a net transverse motion, which explains the angular momentum rich nature of the later infall which comes from further away. We conjecture that this large-scale driven consistency explains why cold flows are so efficient at building up high-redshift thin discs inside out.

Building merger trees from cosmological N-body simulations Towards improving galaxy formation models using subhaloes

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.