Julien Devriendt

Modelling the galaxy bimodality: shutdown above a critical halo mass

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 and Birnboim. After a critical epoch z ∼ 3, only haloes below a critical shock-heating mass M shock ∼ 10 12 Menjoy 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.

galics– I. A hybrid N-body/semi-analytic model of hierarchical galaxy formation

This is the first paper of a series that describes the methods and basic results of the galics model (Galaxies In Cosmological Simulations). galics is a hybrid model for hierarchical galaxy formation studies, combining the outputs of large cosmological N-body simulations with simple, semi-analytic recipes to describe the fate of the baryons within dark matter haloes. The simulations produce a detailed merging tree for the dark matter haloes, including complete knowledge of the statistical properties arising from the gravitational forces. We intend to predict the overall statistical properties of galaxies, with special emphasis on the panchromatic spectral energy distribution emitted by galaxies in the ultraviolet/optical and infrared/submillimetre wavelength ranges. In this paper, we outline the physically motivated assumptions and key free parameters that go into the model, comparing and contrasting with other parallel efforts. We specifically illustrate the success of the model in comparison with several data sets, showing how it is able to predict the galaxy disc sizes, colours, luminosity functions from the ultraviolet to far infrared, the Tully–Fisher and Faber–Jackson relations, and the fundamental plane in the local Universe. We also identify certain areas where the model fails, or where the assumptions needed to succeed are at odds with observations, and pay special attention to understanding the effects of the finite resolution of the simulations on the predictions made. Other papers in this series will take advantage of different data sets available in the literature to extend the study of the limitations and predictive power of galics, with particular emphasis put on high-redshift galaxies.

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.

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.

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.

Formation of structure in molecular clouds: A Case study

Molecular clouds (MCs) are highly structured and turbulent. Colliding gas streams of atomic hydrogen have been suggested as a possible source of MCs, imprinting the filamentary structure as a consequence of dynamical and thermal instabilities. We present a two-dimensional numerical analysis of MC formation via converging H I flows. Even with modest flow speeds and completely uniform inflows, nonlinear density perturbations arise as possible precursors of MCs. Thus, we suggest that MCs are inevitably formed with substantial structure, e.g., strong density and velocity fluctuations, which provide the initial conditions for subsequent gravitational collapse and star formation in a variety of Galactic and extragalactic environments.

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.