Source author record

Julien Devriendt

Julien Devriendt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.HE astro-ph.IM astro-ph astro-ph.SR

Catalog footprint

What is connected

67works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Catalogues of voids as antihalos in the local Universe

A recently-proposed algorithm identifies voids in simulations as the regions associated with halos when the initial overdensity field is negated. We apply this method to the real Universe by running a suite of constrained simulations of the 2M++ volume with initial conditions inferred by the BORG algorithm, along with the corresponding inverted set. Our 101 inverted and uninverted simulations, spanning the BORG posterior, each identify ~150,000 "voids as antihalos" with mass exceeding $4.38\times10^{11} \: \mathrm{M_\odot}$ (100 particles) at $z=0$ in a full-sky sphere of radius 155 Mpc/h around the Milky Way. We calculate the size function, volume filling fraction, ellipticity, central density, specific angular momentum, clustering and stacked density profile of the voids, and cross-correlate them with those produced by VIDE on the same simulations. We make our antihalo and VIDE catalogues publicly available.

preprint2022arXiv

Cosmological simulations of the same spiral galaxy: the impact of baryonic physics

The interplay of star formation and supernova (SN) feedback in galaxy formation is a key element for understanding galaxy evolution. Since these processes occur at small scales, it is necessary to have sub-grid models that recover their evolution and environmental effects at the scales reached by cosmological simulations. We simulate the same spiral galaxy inhabiting a Milky Way (MW) size halo in a cosmological environment changing the sub-grid models for SN feedback and star formation. We test combinations of the Schmidt law and a multi-freefall based star formation with delayed cooling feedback or mechanical feedback. We reach a resolution of 35 pc in a zoom-in box of 36 Mpc. For this, we use the code RAMSES with the implementation of gas turbulence in time and trace the local hydrodynamical features of the star-forming gas. Finally, we compare the galaxies at redshift 0 with global and interstellar medium observations in the MW and local spiral galaxies. The simulations show successful comparisons with observations. Nevertheless, diverse galactic morphologies are obtained from different numerical implementations. We highlight the importance of detailed modelling of the star formation and feedback processes, especially when increasing the resolution of simulations. Future improvements could alleviate the degeneracies exhibited in our simulated galaxies under different sub-grid models.

preprint2022arXiv

Impact of Radiation Feedback on the Formation of Globular Cluster Candidates during Cloud-Cloud Collisions

To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass $M_{\rm GC} \sim 10^5\,M_\odot$ forms within $\approx 3\,{\rm Myr}$ when two GMCs with mass $M_{\rm GMC}=3.6\times10^5\,M_\odot$ collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of $\sim 2$ when the radiation feedback is neglected, and the difference is slightly more pronounced (16\%) when extreme $\rm Lyα$ feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars.

preprint2022arXiv

MgII in the JWST Era: a Probe of Lyman Continuum Escape?

Limited constraints on the evolution of the Lyman Continuum (LyC) escape fraction represent one of the primary uncertainties in the theoretical determination of the reionization history. Due to the intervening intergalactic medium (IGM), the possibility of observing LyC photons directly in the epoch of reionization is highly unlikely. For this reason, multiple indirect probes of LyC escape have been identified, some of which are used to identify low-redshift LyC leakers (e.g. O32), while others are primarily useful at $z>6$ (e.g. [OIII]/[CII] far infrared emission). The flux ratio of the resonant MgII doublet emission at 2796$\dot{\rm A}$ and 2803$\dot{\rm A}$ as well as the MgII optical depth have recently been proposed as ideal diagnostics of LyC leakage that can be employed at $z>6$ with JWST. Using state-of-the-art cosmological radiation hydrodynamics simulations post-processed with CLOUDY and resonant-line radiative transfer, we test whether MgII is indeed a useful probe of LyC leakage. Our simulations indicate that the majority of bright, star-forming galaxies with high LyC escape fractions are expected to be MgII emitters rather than absorbers at $z=6$. However, we find that the MgII doublet flux ratio is a more sensitive indicator of dust rather than neutral hydrogen, limiting its use as a LyC leakage indicator to only galaxies in the optically thin regime. Given its resonant nature, we show that MgII will be an exciting probe of the complex kinematics in high-redshift galaxies in upcoming JWST observations.

preprint2022arXiv

On the Viability of Determining Galaxy Properties from Observations I: Star Formation Rates and Kinematics

We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes we ascertain which physical properties of the galaxy will be recoverable with the HARMONI spectrograph on the European Extremely Large Telescope (ELT). We are able to estimate the galaxy's star formation rate and dynamical mass to a reasonable degree of accuracy, with values within a factor of $1.81$ and $1.38$ of the true value. The kinematic structure of the galaxy is also recovered in mock observations. Furthermore, we are able to recover radial profiles of the velocity dispersion and are therefore able to calculate how the dynamical ratio varies as a function of distance from the galaxy centre. Finally, we show that when calculated on galaxy scales the dynamical ratio does not always provide a reliable measure of a galaxy's stability against gravity or act as an indicator of a minor merger.

preprint2022arXiv

Simulating Jellyfish Galaxies: A Case Study for a Gas-Rich Dwarf Galaxy

We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multi-phase interstellar medium (ISM). We find that the ram-pressure-stripped (RPS) ISM is the dominant source of molecular clumps in the near wake within 10 kpc from the galactic plane, while in-situ formation is the major channel for dense gas in the distant tail of the gas-rich galaxy. Only 20% of the molecular clumps in the near wake originate from the intracluster medium (ICM); however, the fraction reaches 50% in the clumps located at $80\,{\rm kpc}$ from the galactic center since the cooling time of the RPS gas tends to be short due to the ISM--ICM mixing ($\lesssim$ 10 Myr). The tail region exhibits a star formation rate of $0.001-0.01\,{\rm M_{\odot}\,yr^{-1}}$, and most of the tail stars are born in the stripped wake within 10 kpc from the galactic plane. These stars induce bright H$α$ blobs in the tail, while H$α$ tails fainter than $6\times10^{38}\,{\rm erg\,s^{-1}\,kpc^{-2}}$ are mostly formed via collisional radiation and heating due to mixing. We also find that the stripped tails have intermediate X-ray to H$α$ surface brightness ratios (1.5$\lesssim F_{\rm X}/F_{\rm Hα}\lesssim$20), compared to the ISM ($\lesssim$1.5) or pure ICM ($\gg$20). Our results suggest that jellyfish features emerge when the ISM from gas-rich galaxies is stripped by strong ram pressure, mixes with the ICM, and enhances the cooling in the tail.

preprint2022arXiv

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged)

preprint2022arXiv

The Obelisk simulation: galaxies contribute more than AGN to HI reionization of protoclusters

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

preprint2022arXiv

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy.

preprint2021arXiv

Momentum deposition of supernovae with cosmic rays

The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We explore multiple ambient densities, magnetic fields and fractions of supernova energy deposited as CRs ($χ_{\rm CR}$), accounting for cosmic ray anisotropic diffusion and streaming. All our runs have comparable evolutions until the end of the Sedov-Taylor phase. However, our CRMHD simulations experience an additional CR pressure-driven snowplough phase once the CR energy dominates inside the SNR. We present a model for the final momentum deposited by supernovae that captures this new phase: $p_{\rm SNR} = 2.87\times 10^{5} (χ_{\text{CR}} + 1)^{4.82}\left(\frac{n}{\text{cm}^{-3}}\right)^{-0.196} M_{\odot}$ km s$^{-1}$. Assuming a 10% fraction of SN energy in CRs leads to a 50% boost of the final momentum, with our model predicting even higher impacts at lower ambient densities. The anisotropic diffusion of CRs assuming an initially uniform magnetic field leads to extended gas and cosmic ray outflows escaping from the supernova poles. We also study a tangled initial configuration of the magnetic field, resulting instead in a quasi-isotropic diffusion of CRs and earlier momentum deposition. Finally, synthetic synchrotron observations of our simulations using the POLARIS code show that the local magnetic field configuration in the interstellar medium modifies the overall radio emission morphology and polarisation.

preprint2021arXiv

Predicting the Observability of Population III Stars with ELT-HARMONI via the Helium $1640{\rmÅ}$ emission line

Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-processing galaxies from the cosmological, AMR-hydrodynamical simulation NewHorizon with theoretical spectral energy distributions (SED) for Pop. III stars and radiative transfer (i.e. the Yggdrasil Models and CLOUDY look-up tables respectively) we are able to compute the flux of HeII$\lambda1640$ for individual galaxies. From mock 10 hour observations of these galaxies we show that HARMONI will be able to detect Pop. III stars in galaxies up to $z\sim10$ provided Pop. III stars have a top heavy Initial Mass Function (IMF). Furthermore, we find that should Pop. III stars instead have an IMF similar to those of the Pop. I stars, the HeII$\lambda1640$ line would only be observable for galaxies with Pop. III stellar masses in excess of $10^{7}{\,M}_\odot\,{\rm yr}^{-1}$, average stellar age $<1{\rm Myr}$ at $z=4$. Finally, we are able to determine the minimal intrinsic flux required for HARMONI to detect Pop. III stars in a galaxy up to $z=10$.

preprint2021arXiv

Rivers of Gas I.: Unveiling The Properties of High Redshift Filaments

At high redshift, the cosmic web is widely expected to have a significant impact on the morphologies, dynamics and star formation rates of the galaxies embedded within it, underscoring the need for a comprehensive study of the properties of such a filamentary network. With this goal in mind, we perform an analysis of high-$z$ gas and dark matter (DM) filaments around a Milky Way-like progenitor simulated with the {\sc ramses} adaptive mesh refinement (AMR) code from cosmic scales ($\sim$1 Mpc) down to the virial radius of its DM halo host ($\sim$20 kpc at $z=4$). Radial density profiles of both gas and DM filaments are found to have the same functional form, namely a plummer-like profile modified to take into account the wall within which these filaments are embedded. Measurements of the typical filament core radius $r_0$ from the simulation are consistent with that of isothermal cylinders in hydrostatic equilibrium. Such an analytic model also predicts a redshift evolution for the core radius of filaments in fair agreement with the measured value for DM $(r_0 \propto (1+z)^{-3.18\pm 0.28})$. Gas filament cores grow as $(r_0 \propto (1+z)^{-2.72\pm 0.26})$. In both gas and DM, temperature and vorticity sharply drop at the edge of filaments, providing an excellent way to constrain the outer filament radius. When feedback is included the gas temperature and vorticity fields are strongly perturbed, hindering such a measurement in the vicinity of the galaxy. However, the core radius of the filaments as measured from the gas density field is largely unaffected by feedback, and the median central density is only reduced by about 20%.

preprint2020arXiv

Black hole mergers from dwarf to massive galaxies with the NewHorizon and Horizon-AGN simulations

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. We use the same method to analyse the two simulations, post-processing MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies.

preprint2020arXiv

Detecting the Cosmic Web: Lyα Emission from Simulated Filaments at z=3

The standard cosmological model ($Λ$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$α$ detection, motivating the development of modern theoretical predictions. We use hydrodynamical cosmological simulations run with the AREPO code to investigate the potential detectability of large-scale filaments, excluding contributions from the halos embedded in them. We focus on filaments connecting massive ($M_{200c}\sim(1-3)\times10^{12} M_\odot$) halos at z=3, and compare different simulation resolutions, feedback levels, and mock-image pixel sizes. We find increasing simulation resolution does not substantially improve detectability notwithstanding the intrinsic enhancement of internal filament structure. By contrast, for a MUSE integration of 31 hours, including feedback increases the detectable area by a factor of $\simeq$5.5 on average compared with simulations without feedback, implying that even the non-bound components of the filaments have substantial sensitivity to feedback. Degrading the image resolution from the native MUSE scale of (0.2")$^2$ per pixel to (5.3")$^2$ apertures has the strongest effect, increasing the detectable area by a median factor of $\simeq$200 and is most effective when the size of the pixel roughly matches the width of the filament. Finally, we find the majority of Ly$α$ emission is due to electron impact collisional excitations, as opposed to radiative recombination.

preprint2020arXiv

GalICS 2.1: a new semianalytic model for cold accretion, cooling, feedback and their roles in galaxy formation

Dekel & Birnboim06 proposed that the mass-scale that separates late-type and early-type galaxies is linked to the critical halo mass for the propagation of a stable shock and showed that they could reproduce the observed bimodality scale for plausible values of the metallicity of the accreted gas and the shock radius. Here, we take their analysis one step further and present a new semianalytic model that computes the shock radius from first principles. This advancement allows us to compute the critical mass individually for each halo. Separating cold-mode and hot-mode accretion has little effect on the final galaxy masses if feedback does not preferentially couple to the hot gas. We also present an improved model for stellar feedback where 70% of the wind mass is in a cold galactic fountain with a shorter reaccretion timescale at high masses. The latter is the key mechanism that allows us to reproduce the low-mass end of the mass function of galaxies over the entire redshift range 0<z<2.5. Cooling must be mitigated to avoid overpredicting the number density of galaxies with stellar mass greater than a hundred billion Solar masses but is important to form intermediate-mass galaxies. At virial masses greater than three hundred billion Solar masses, cold accretion is more important at high z, where gas is accreted from smaller solid angles, but this is not true at lower masses because high-z filaments have lower metallicities. Our predictions are consistent with the observed metallicity evolution of the intergalactic medium at 0<z<5.

preprint2020arXiv

Horizon-AGN virtual observatory -- 2: Template-free estimates of galaxy properties from colours

Using the Horizon-AGN hydrodynamical simulation and self-organising maps (SOMs), we show how to compress the complex data structure of a cosmological simulation into a 2-d grid which is much easier to analyse. We first verify the tight correlation between the observed 0.3$\!-\!5μ$m broad-band colours of Horizon-AGN galaxies and their high-resolution spectra. The correlation is found to extend to physical properties such as redshift, stellar mass, and star formation rate (SFR). This direct mapping from colour to physical parameter space is shown to work also after including photometric uncertainties that mimic the COSMOS survey. We then label the SOM grid with a simulated calibration sample and estimate redshift and SFR for COSMOS-like galaxies up to $z\sim3$. In comparison to state-of-the-art techniques based on synthetic templates, our method is comparable in performance but less biased at estimating redshifts, and significantly better at predicting SFRs. In particular our "data-driven" approach, in contrast to model libraries, intrinsically allows for the complexity of galaxy formation and can handle sample biases. We advocate that obtaining the calibration for this method should be one of the goals of next-generation galaxy surveys.

preprint2020arXiv

How primordial magnetic fields shrink galaxies

As one of the prime contributors to the interstellar medium energy budget, magnetic fields naturally play a part in shaping the evolution of galaxies. Galactic magnetic fields can originate from strong primordial magnetic fields provided these latter remain below current observational upper limits. To understand how such magnetic fields would affect the global morphological and dynamical properties of galaxies, we use a suite of high-resolution constrained transport magneto-hydrodynamic cosmological zoom simulations where we vary the initial magnetic field strength and configuration along with the prescription for stellar feedback. We find that strong primordial magnetic fields delay the onset of star formation and drain the rotational support of the galaxy, diminishing the radial size of the galactic disk and driving a higher amount of gas towards the centre. This is also reflected in mock UVJ observations by an increase in the light profile concentration of the galaxy. We explore the possible mechanisms behind such a reduction in angular momentum, focusing on magnetic braking. Finally, noticing that the effects of primordial magnetic fields are amplified in the presence of stellar feedback, we briefly discuss whether the changes we measure would also be expected for galactic magnetic fields of non-primordial origin.

preprint2020arXiv

How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $Δz>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

preprint2020arXiv

New Methods for Identifying Lyman Continuum Leakers and Reionization-Epoch Analogues

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, local [SII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SII-BPT diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158μm}$ and [OIII]$_{\rm 88μm}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.

preprint2020arXiv

The impact of AGN feedback on galaxy intrinsic alignments in the Horizon simulations

The intrinsic correlations of galaxy shapes and orientations across the large-scale structure of the Universe are a known contaminant to weak gravitational lensing. They are known to be dependent on galaxy properties, such as their mass and morphologies. The complex interplay between alignments and the physical processes that drive galaxy evolution remains vastly unexplored. We assess the sensitivity of intrinsic alignments (shapes and angular momenta) to Active Galactic Nuclei -AGN- feedback by comparing galaxy alignment in twin runs of the cosmological hydrodynamical Horizon simulation, which do and do not include AGN feedback respectively. We measure intrinsic alignments in three dimensions and in projection at z=0 and z=1. We find that the projected alignment signal of all galaxies with resolved shapes with respect to the density field in the simulation is robust to AGN feedback, thus giving similar predictions for contamination to weak lensing. The relative alignment of galaxy shapes around galaxy positions is however significantly impacted, especially when considering high-mass ellipsoids. Using a sample of galaxy "twins" across simulations, we determine that AGN changes both the galaxy selection and their actual alignments. Finally, we measure the alignments of angular momenta of galaxies with their nearest filament. Overall, these are more significant in the presence of AGN as a result of the higher abundance of massive pressure-supported galaxies.

preprint2020arXiv

Why do extremely massive disc galaxies exist today?

Galaxy merger histories correlate strongly with stellar mass, largely regardless of morphology. Thus, at fixed stellar mass, spheroids and discs share similar assembly histories, both in terms of the frequency of mergers and the distribution of their mass ratios. Since mergers are the principal drivers of disc-to-spheroid morphological transformation, and the most massive galaxies typically have the richest merger histories, it is surprising that discs exist at all at the highest stellar masses (e.g. beyond the knee of the mass function). Using Horizon-AGN, a cosmological hydro-dynamical simulation, we show that extremely massive (M*> 10^11.4 MSun) discs are created via two channels. In the primary channel (accounting for ~70% of these systems and ~8% of massive galaxies) the most recent, significant merger (stellar mass ratio > 1:10) between a massive spheroid and a gas-rich satellite `spins up' the spheroid by creating a new rotational stellar component, leaving a massive disc as the remnant. In the secondary channel (accounting for ~30% of these systems and ~3% of massive galaxies), a system maintains a disc throughout its lifetime, due to an anomalously quiet merger history. Not unexpectedly, the fraction of massive discs is larger at higher redshift, due to the Universe being more gas-rich. The morphological mix of galaxies at the highest stellar masses is, therefore, a strong function of the gas fraction of the Universe. Finally, these massive discs have similar black-hole masses and accretion rates to massive spheroids, providing a natural explanation for why a minority of powerful AGN are surprisingly found in disc galaxies.

preprint2019arXiv

The impact of the connectivity of the cosmic web on the physical properties of galaxies at its nodes

We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall).

preprint2019arXiv

When galaxies align: intrinsic alignments of the progenitors of elliptical galaxies in the Horizon-AGN simulation

Elliptical galaxies today appear aligned with the large-scale structure of the Universe, but it is still an open question when they acquire this alignment. Observational data is currently insufficient to provide constraints on the time evolution of intrinsic alignments, and hence existing models range from assuming that galaxies gain some primordial alignment at formation, to suggesting that they react instantaneously to tidal interactions with the large-scale structure. Using the cosmological hydrodynamical simulation Horizon-AGN, we measure the relative alignments between the major axes of galaxies and eigenvectors of the tidal field as a function of redshift. We focus on constraining the time evolution of the alignment of the main progenitors of massive $z=0$ elliptical galaxies, the main weak lensing contaminant at low redshift. We show that this population, which at $z=0$ has a stellar mass above $10^{10.4}$ M$_\odot$, transitions from having no alignment with the tidal field at $z=3$, to a significant alignment by $z=1$. From $z=0.5$ they preserve their alignment at an approximately constant level until $z=0$. We find a mass-dependence of the alignment signal of elliptical progenitors, whereby ellipticals that are less massive today ($10^{10.4}<M/{\rm M}_\odot<10^{10.7}$) do not become aligned till later redshifts ($z<2$), compared to more massive counterparts. We also present an extended study of progenitor alignments in the parameter space of stellar mass and galaxy dynamics, the impact of shape definition and tidal field smoothing.

preprint2016arXiv

Bursty star formation feedback and cooling outflows

We study how outflows of gas launched from a central galaxy undergoing repeated starbursts propagate through the circumgalactic medium (CGM), using the simulation code RAMSES. We assume that the outflow from the disk can be modelled as a rapidly moving bubble of hot gas at $\mathrm{\sim1\;kpc}$ above disk, then ask what happens as it moves out further into the halo around the galaxy on $\mathrm{\sim 100\;kpc}$ scales. To do this we run 60 two-dimensional simulations scanning over parameters of the outflow. Each of these is repeated with and without radiative cooling, assuming a primordial gas composition to give a lower bound on the importance of cooling. In a large fraction of radiative-cooling cases we are able to form rapidly outflowing cool gas from in situ cooling of the flow. We show that the amount of cool gas formed depends strongly on the 'burstiness' of energy injection; sharper, stronger bursts typically lead to a larger fraction of cool gas forming in the outflow. The abundance ratio of ions in the CGM may therefore change in response to the detailed historical pattern of star formation. For instance, outflows generated by star formation with short, intense bursts contain up to 60 per cent of their gas mass at temperatures $<5 \times 10^4\,\mathrm{K}$; for near-continuous star formation the figure is $\lesssim$ 5 per cent. Further study of cosmological simulations, and of idealised simulations with e.g., metal-cooling, magnetic fields and/or thermal conduction, will help to understand the precise signature of bursty outflows on observed ion abundances.

preprint2016arXiv

Caught in the rhythm: how satellites settle into a plane around their central galaxy

Using the cosmological hydrodynamics simulation Horizon-AGN, we investigate the spatial distribution of satellite galaxies relative to their central counterpart in the redshift range between 0.3 and 0.8. We find that, on average, these satellites tend to be located on the galactic plane of the central object. This effect is detected for central galaxies with a stellar mass larger than 10^10 solar masses and found to be strongest for red passive galaxies, while blue galaxies exhibit a weaker trend. For galaxies with a minor axis parallel to the direction of the nearest filament, we find that the coplanarity is stronger in the vicinity of the central galaxy, and decreases when moving towards the outskirts of the host halo. By contrast, the spatial distribution of satellite galaxies relative to their closest filament follows the opposite trend: their tendency to align with them dominates at large distances from the central galaxy, and fades away in its vicinity. Relying on mock catalogs of galaxies in that redshift range, we show that massive red centrals with a spin perpendicular to their filament also have corotating satellites well aligned with both the galactic plane and the filament. On the other hand, lower-mass blue centrals with a spin parallel to their filament have satellites flowing straight along this filament, and hence orthogonally to their galactic plane. The orbit of these satellites is then progressively bent towards a better alignment with the galactic plane as they penetrate the central region of their host halo. The kinematics previously described are consistent with satellite infall and spin build-up via quasi-polar flows, followed by a re-orientation of the spin of massive red galaxies through mergers.

preprint2016arXiv

Comparing Simulations of AGN Feedback

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

preprint2016arXiv

Modeling Lyman-α Forest Cross-Correlations with LyMAS

We use the Ly-$α$ Mass Association Scheme (LyMAS; Peirani et al. 2014) to predict cross-correlations at $z=2.5$ between dark matter halos and transmitted flux in the Ly-$α$ forest, and compare to cross-correlations measured for quasars and damped Ly-$α$ systems (DLAs) from the Baryon Oscillation Spectroscopic Survey (BOSS) by Font-Ribera et al. (2012, 2013). We calibrate LyMAS using Horizon-AGN hydrodynamical cosmological simulations of a $(100\ h^{-1}\ \mathrm{Mpc})^3$ comoving volume. We apply this calibration to a $(1\ h^{-1}\ \mathrm{Gpc})^3$ simulation realized with $2048^3$ dark matter particles. In the 100 $h^{-1}$ Mpc box, LyMAS reproduces the halo-flux correlations computed from the full hydrodynamic gas distribution very well. In the 1 $h^{-1}$ Gpc box, the amplitude of the large scale cross-correlation tracks the halo bias $b_h$ as expected. We provide empirical fitting functions that describe our numerical results. In the transverse separation bins used for the BOSS analyses, LyMAS cross-correlation predictions follow linear theory accurately down to small scales. Fitting the BOSS measurements requires inclusion of random velocity errors; we find best-fit RMS velocity errors of 399 km s$^{-1}$ and 252 km s$^{-1}$ for quasars and DLAs, respectively. We infer bias-weighted mean halo masses of $M_h/10^{12}\ h^{-1}M_\odot=2.19^{+0.16}_{-0.15}$ and $0.69^{+0.16}_{-0.14}$ for the host halos of quasars and DLAs, with $\sim 0.2$ dex systematic uncertainty associated with redshift evolution, IGM parameters, and selection of data fitting range.

preprint2016arXiv

Redshift and luminosity evolution of the intrinsic alignments of galaxies in Horizon-AGN

Intrinsic galaxy shape and angular momentum alignments can arise in cosmological large-scale structure due to tidal interactions or galaxy formation processes. Cosmological hydrodynamical simulations have recently come of age as a tool to study these alignments and their contamination to weak gravitational lensing. We probe the redshift and luminosity evolution of intrinsic alignments in Horizon-AGN between z=0 and z=3 for galaxies with an r-band absolute magnitude of <-20. Alignments transition from being radial at low redshifts and high luminosities, dominated by the contribution of ellipticals, to being tangential at high redshift and low luminosities, where discs dominate the signal. This cannot be explained by the evolution of the fraction of ellipticals and discs alone: intrinsic evolution in the amplitude of alignments is necessary. The alignment amplitude of elliptical galaxies alone is smaller in amplitude by a factor of ~2, but has similar luminosity and redshift evolution as in current observations and in the nonlinear tidal alignment model at projected separations of > 1 Mpc. Alignments of discs are null in projection and consistent with current low redshift observations. The combination of the two populations yields an overall amplitude a factor of ~4 lower than observed alignments of luminous red galaxies with a steeper luminosity dependence. The restriction on accurate galaxy shapes implies that the galaxy population in the simulation is complete only to an r-band absolute magnitude of <-20. Higher resolution simulations will be necessary to avoid extrapolation of the intrinsic alignment predictions to the range of luminosities probed by future surveys.

preprint2016arXiv

The cosmic evolution of massive black holes in the Horizon-AGN simulation

We analyse the demographics of black holes (BHs) in the large-volume cosmological hydrodynamical simulation Horizon-AGN. This simulation statistically models how much gas is accreted onto BHs, traces the energy deposited into their environment and, consequently, the back-reaction of the ambient medium on BH growth. The synthetic BHs reproduce a variety of observational constraints such as the redshift evolution of the BH mass density and the mass function. Strong self-regulation via AGN feedback, weak supernova feedback, and unresolved internal processes result in a tight BH-galaxy mass correlation. Starting at z~2, tidal stripping creates a small population of BHs over-massive with respect to the halo. The fraction of galaxies hosting a central BH or an AGN increases with stellar mass. The AGN fraction agrees better with multi-wavelength studies, than single-wavelength ones, unless obscuration is taken into account. The most massive halos present BH multiplicity, with additional BHs gained by ongoing or past mergers. In some cases, both a central and an off-centre AGN shine concurrently, producing a dual AGN. This dual AGN population dwindles with decreasing redshift, as found in observations. Specific accretion rate and Eddington ratio distributions are in good agreement with observational estimates. The BH population is dominated in turn by fast, slow, and very slow accretors, with transitions occurring at z=3 and z=2 respectively.

preprint2016arXiv

The Horizon-AGN simulation: morphological diversity of galaxies promoted by AGN feedback

The interplay between cosmic gas accretion on to galaxies and galaxy mergers drives the observed morphological diversity of galaxies. By comparing the state-of-the-art hydrodynamical cosmological simulations Horizon-AGN and Horizon-noAGN, we unambiguously identify the critical role of active galactic nuclei (AGN) in setting up the correct galaxy morphology for the massive end of the population. With AGN feedback, typical kinematic and morpho-metric properties of galaxy populations as well as the galaxy-halo mass relation are in much better agreement with observations. Only AGN feedback allows massive galaxies at the centre of groups and clusters to become ellipticals, while without AGN feedback those galaxies reform discs. It is the merger-enhanced AGN activity that is able to freeze the morphological type of the post-merger remnant by durably quenching its quiescent star formation. Hence morphology is shown to be driven not only by mass but also by the nature of cosmic accretion: at constant galaxy mass, ellipticals are galaxies that are mainly assembled through mergers, while discs are preferentially built from the in situ star formation fed by smooth cosmic gas infall.

preprint2015arXiv

A Detailed Study of Feedback from a Massive Star

We present numerical simulations of a 15 solar mass star in a suite of idealised environments in order to quantify the amount of energy transmitted to the interstellar medium (ISM). We include models of stellar winds, UV photoionisation and the subsequent supernova based on theoretical models and observations of stellar evolution. The system is simulated in 3D using RAMSES-RT, an Adaptive Mesh Refinement Radiation Hydrodynamics code. We find that stellar winds have a negligible impact on the system owing to their relatively low luminosity compared to the other processes. The main impact of photoionisation is to reduce the density of the medium into which the supernova explodes, reducing the rate of radiative cooling of the subsequent supernova. Finally, we present a grid of models quantifying the energy and momentum of the system that can be used to motivate simulations of feedback in the ISM unable to fully resolve the processes discussed in this work.

preprint2015arXiv

Black hole evolution: I. Supernova-regulated black hole growth

The growth of a supermassive black hole (BH) is determined by how much gas the host galaxy is able to feed it, which in turn is controlled by the cosmic environment, through galaxy mergers and accretion of cosmic flows that time how galaxies obtain their gas, but also by internal processes in the galaxy, such as star formation and feedback from stars and the BH itself. In this paper, we study the growth of a 10^12 Msun halo at z=2, which is the progenitor of al group of galaxies at z=0, and of its central BH by means of a high-resolution zoomed cosmological simulation, the Seth simulation. We study the evolution of the BH driven by the accretion of cold gas in the galaxy, and explore the efficiency of the feedback from supernovae (SNe). For a relatively inefficient energy input from SNe, the BH grows at the Eddington rate from early times, and reaches self-regulation once it is massive enough. We find that at early cosmic times z>3.5, efficient feedback from SNe forbids the formation of a settled disc as well as the accumulation of dense cold gas in the vicinity of the BH and starves the central compact object. As the galaxy and its halo accumulate mass, they become able to confine the nuclear inflows provided by major mergers and the BH grows at a sustained near-to-Eddington accretion rate. We argue that this mechanism should be ubiquitous amongst low-mass galaxies, corresponding to galaxies with a stellar mass below <10^9 Msun in our simulations.

preprint2015arXiv

Intrinsic alignment of simulated galaxies in the cosmic web: implications for weak lensing surveys

The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z=1.2 synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation Horizon-AGN. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field `filamentary' direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalog are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to 10 Mpc/h comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate mass galaxies on a similar scale at a level of 10^(-4) which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step toward the use of very realistic catalog of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.

preprint2015arXiv

Intrinsic alignments of galaxies in the Horizon-AGN cosmological hydrodynamical simulation

The intrinsic alignments of galaxies are recognised as a contaminant to weak gravitational lensing measurements. In this work, we study the alignment of galaxy shapes and spins at low redshift ($z\sim 0.5$) in Horizon-AGN, an adaptive-mesh-refinement hydrodynamical cosmological simulation box of 100 Mpc/h a side with AGN feedback implementation. We find that spheroidal galaxies in the simulation show a tendency to be aligned radially towards over-densities in the dark matter density field and other spheroidals. This trend is in agreement with observations, but the amplitude of the signal depends strongly on how shapes are measured and how galaxies are selected in the simulation. Disc galaxies show a tendency to be oriented tangentially around spheroidals in three-dimensions. While this signal seems suppressed in projection, this does not guarantee that disc alignments can be safely ignored in future weak lensing surveys. The shape alignments of luminous galaxies in Horizon-AGN are in agreement with observations and other simulation works, but we find less alignment for lower luminosity populations. We also characterize the systematics of galaxy shapes in the simulation and show that they can be safely neglected when measuring the correlation of the density field and galaxy ellipticities.

preprint2015arXiv

Milking the spherical cow: on aspherical dynamics in spherical coordinates

Galaxies and the dark matter halos that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealised calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are "maximally stable", i.e. that do not evolve at first order when external potentials (which arise from baryons, large scale tidal fields or infalling substructure) are applied. We show that a spherically-symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealised spherically-symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter halos (formed in collisionless cosmological simulations) nearly attain our maximally-stable limit, but their evolution freezes out before reaching it.

preprint2015arXiv

Towards simulating star formation in turbulent high-z galaxies with mechanical supernova feedback

Feedback from supernovae is essential to understanding the self-regulation of star formation in galaxies. However, the efficacy of the process in a cosmological context remains unclear due to excessive radiative losses during the shock propagation. To better understand the impact of SN explosions on the evolution of galaxies, we perform a suite of high-resolution (12 pc), zoom-in cosmological simulations of a Milky Way-like galaxy at z=3 with adaptive mesh refinement. We find that SN explosions can efficiently regulate star formation, leading to the stellar mass and metallicity consistent with the observed mass-metallicity relation and stellar mass-halo mass relation at z~3. This is achieved by making three important changes to the classical feedback scheme: i) the different phases of SN blast waves are modelled directly by injecting radial momentum expected at each stage, ii) the realistic time delay of SNe, commencing at as early as 3 Myr, is required to disperse very dense gas before a runaway collapse sets in at the galaxy centre via mergers of gas clumps, and iii) a non-uniform density distribution of the ISM is taken into account below the computational grid scale for the cell in which SN explodes. The last condition is motivated by the fact that our simulations still do not resolve the detailed structure of a turbulent ISM in which the fast outflows can propagate along low-density channels. The simulated galaxy with the SN feedback model shows strong outflows, which carry approximately ten times larger mass than star formation rate, as well as smoothly rising circular velocity. Other feedback models that do not meet the three conditions form too many stars, producing a peaked rotation curve. Our results suggest that understanding the structure of the turbulent ISM may be crucial to assess the role of SN and other feedback processes in galaxy formation theory.

preprint2014arXiv

Black hole evolution: II. Spinning black holes in a supernova-driven turbulent interstellar medium

Supermassive black holes (BH) accrete gas from their surroundings and coalesce with companions during galaxy mergers, and both processes change the BH mass and spin. By means of high-resolution hydrodynamical simulations of galaxies, either idealised or embedded within the cosmic web, we explore the effects of interstellar gas dynamics and external perturbations on BH spin evolution. All these physical quantities were evolved on-the-fly in a self-consistent manner. We use a `maximal' model to describe the turbulence induced by stellar feedback to highlight its impact on the angular momentum of the gas accreted by the BH. Periods of intense star formation are followed by phases where stellar feedback drives large-scale outflows and hot bubbles. We find that BH accretion is synchronised with star formation, as only when gas is cold and dense do both processes take place. During such periods, gas motion is dominated by consistent rotation. On the other hand, when stellar feedback becomes substantial, turbulent motion randomises gas angular momentum. However BH accretion is strongly suppressed in that case, as cold and dense gas is lacking. In our cosmological simulation, at very early times (z>6), the galactic disc has not yet settled and no preferred direction exists for the angular momentum of the accreted gas, so the BH spin remains low. As the gas settles into a disc (6>z>3), the BH spin then rapidly reaches its maximal value. At lower redshifts (z<3), even when galaxy mergers flip the direction of the angular momentum of the accreted gas, causing it to counter-rotate, the BH spin magnitude only decreases modestly and temporarily. Should this be a typical evolution scenario for BH, it potentially has dramatic consequences regarding their origin and assembly, as accretion on maximally spinning BH embedded in thin Shakura-Sunyaev disc is significantly reduced.

preprint2014arXiv

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 cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. 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 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 and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

preprint2014arXiv

Mergers drive spin swings along the cosmic web

The close relationship between mergers and the reorientation of the spin for galaxies and their host dark haloes is investigated using a cosmological hydrodynamical simulation (Horizon-AGN). Through a statistical analysis of merger trees, we show that spin swings are mainly driven by mergers along the filamentary structure of the cosmic web, and that these events account for the preferred perpendicular orientation of massive galaxies with respect to their nearest filament. By contrast, low-mass galaxies (M_s<10^10 M_sun at redshift 1.5) undergoing very few mergers, if at all, tend to possess a spin well aligned with their filament. Haloes follow the same trend as galaxies but display a greater sensitivity to smooth anisotropic accretion. The relative effect of mergers on spin magnitude is qualitatively different for minor and major mergers: mergers (and diffuse accretion) generally increase the magnitude of the angular momentum, but the most massive major mergers also give rise to a population of objects with less spin left. Without mergers secular accretion builds up the spin of galaxies but not that of haloes. It also (re)aligns galaxies with their filament.

preprint2014arXiv

Swirling around filaments: are large-scale structure vortices spinning up dark halos?

The kinematic analysis of dark matter and hydrodynamical simulations suggests that the vorticity in large-scale structure is mostly confined to, and predominantly aligned with their filaments, with an excess of probability of 20 per cent to have the angle between vorticity and filaments direction lower than 60 degrees relative to random orientations. The cross sections of these filaments are typically partitioned into four quadrants with opposite vorticity sign, arising from multiple flows, originating from neighbouring walls. The spins of halos embedded within these filaments are consistently aligned with this vorticity for any halo mass, with a stronger alignment for the most massive structures up to an excess of probability of 165 per cent. On large scales, adiabatic/cooling hydrodynamical simulations display the same vorticity in the gas as in the dark matter. The global geometry of the flow within the cosmic web is therefore qualitatively consistent with a spin acquisition for smaller halos induced by this large-scale coherence, as argued in Codis et al. (2012). In effect, secondary anisotropic infall (originating from the vortex-rich filament within which these lower-mass halos form) dominates the angular momentum budget of these halos. The transition mass from alignment to orthogonality is related to the size of a given multi-flow region with a given polarity. This transition may be reconciled with the standard tidal torque theory if the latter is augmented so as to account for the larger scale anisotropic environment of walls and filaments.

preprint2012arXiv

Angular momentum transfer to a Milky Way disk at high redshift

An Adaptive Mesh Refinement cosmological resimulation is analyzed in order to test whether filamentary flows of cold gas are responsible for the build-up of angular momentum within a Milky Way like disk at z>=3. A set of algorithms is presented that takes advantage of the high spatial resolution of the simulation (12 pc) to identify: (i) the central gas disk and its plane of orientation; (ii) the complex individual filament trajectories that connect to the disk, and; (iii) the infalling satellites. The results show that two filaments at z>5.5, which later merge to form a single filament at z<4, drive the angular momentum and mass budget of the disk throughout its evolution, whereas luminous satellite mergers make negligible fractional contributions. Combined with the ubiquitous presence of such filaments in all large-scale cosmological simulations that include hydrodynamics, these findings provide strong quantitative evidence that the growth of thin disks in haloes with masses below 10^{12} M_{sun}, which host the vast majority of galaxies, is supported via inflowing streams of cold gas at intermediate and high redshifts.

preprint2012arXiv

Blowing cold flows away: the impact of early AGN activity on the formation of a brightest cluster galaxy progenitor

Supermassive black holes (BH) are powerful sources of energy that are already in place at very early epochs of the Universe (by z=6). Using hydrodynamical simulations of the formation of a massive M_vir=5 10^11 M_sun halo by z=6 (the most massive progenitor of a cluster of M_vir=2 10^15 M_sun at z=0), we evaluate the impact of Active Galactic Nuclei (AGN) on galaxy mass content, BH self-regulation, and gas distribution inside this massive halo. We find that SN feedback has a marginal influence on the stellar structure, and no influence on the mass distribution on large scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge mass content by quenching star formation when the BH is self-regulating, and by depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH proceeds first by a rapid Eddington-limited period fed by direct cold filamentary infall. When the energy delivered by the AGN is sufficiently large to unbind the cold gas of the bulge, the accretion of gas onto the BH is maintained both by smooth gas inflow and clump migration through the galactic disc triggered by merger-induced torques. The feedback from the AGN has also a severe consequence on the baryon mass content within the halo, producing large-scale hot superwinds, able to blow away some of the cold filamentary material from the centre and reduce the baryon fraction by more than 30 per cent within the halo's virial radius. Thus in the very young universe, AGN feedback is likely to be a key process, shaping the properties of the most massive galaxies.

preprint2012arXiv

Connecting the cosmic web to the spin of dark halos: 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 4pi 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 mass tends to be perpendicular to the closest filament, and aligned with the intermediate axis of the tidal tensor, whereas the spin of low-mass halos is more likely to be aligned with the closest filament. Furthermore, this critical mass of 5 10^12 is redshift-dependent and scales as (1+z)^-2.5. We propose an interpretation of this signal in terms of large-scale cosmic flows. In this picture, most low-mass halos are formed through the winding of flows 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.

preprint2012arXiv

Constraining stellar assembly and AGN feedback at the peak epoch of star formation

We study stellar assembly and feedback from active galactic nuclei (AGN) around the epoch of peak star formation (1<z<2), by comparing hydrodynamic simulations to rest-frame UV-optical galaxy colours from the Wide Field Camera 3 (WFC3) Early-Release Science (ERS) Programme. Our Adaptive Mesh Refinement simulations include metal-dependent radiative cooling, star formation, kinetic outflows due to supernova explosions, and feedback from supermassive black holes. Our model assumes that when gas accretes onto black holes, a fraction of the energy is used to form either thermal winds or sub-relativistic momentum-imparting collimated jets, depending on the accretion rate. We find that the predicted rest-frame UV-optical colours of galaxies in the model that includes AGN feedback is in broad agreement with the observed colours of the WFC3 ERS sample at 1<z<2. The predicted number of massive galaxies also matches well with observations in this redshift range. However, the massive galaxies are predicted to show higher levels of residual star formation activity than the observational estimates, suggesting the need for further suppression of star formation without significantly altering the stellar mass function. We discuss possible improvements, involving faster stellar assembly through enhanced star formation during galaxy mergers while star formation at the peak epoch is still modulated by the AGN feedback.

preprint2012arXiv

Feeding compact bulges and supermassive black holes with low angular-momentum cosmic gas at high redshift

We use cosmological hydrodynamical simulations to show that a significant fraction of the gas in high redshift rare massive halos falls nearly radially to their very centre on extremely short timescales. This process results in the formation of very compact bulges with specific angular momentum a factor 5-30$smaller than the average angular momentum of the baryons in the whole halo. Such low angular momentum originates both from segregation and effective cancellation when the gas flows to the centre of the halo along well defined cold filamentary streams. These filaments penetrate deep inside the halo and connect to the bulge from multiple rapidly changing directions. Structures falling in along the filaments (satellite galaxies) or formed by gravitational instabilities triggered by the inflow (star clusters) further reduce the angular momentum of the gas in the bulge. Finally, the fraction of gas radially falling to the centre appears to increase with the mass of the halo; we argue that this is most likely due to an enhanced cancellation of angular momentum in rarer halos which are fed by more isotropically distributed cold streams. Such an increasingly efficient funnelling of low-angular momentum gas to the centre of very massive halos at high redshift may account for the rapid pace at which the most massive supermassive black holes grow to reach observed masses around $10^9$M$_\odot$ at an epoch when the Universe is barely 1 Gyr old.

preprint2012arXiv

Lyman-alpha emission properties of simulated galaxies: interstellar medium structure and inclination effects

[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.

preprint2012arXiv

Satellite Survival in Highly Resolved Milky Way Class Halos

Surprisingly little is known about the origin and evolution of the Milky Way's satellite galaxy companions. UV photoionisation, supernova feedback and interactions with the larger host halo are all thought to play a role in shaping the population of satellites that we observe today, but there is still no consensus as to which of these effects, if any, dominates. In this paper, we revisit the issue by re-simulating a Milky Way class dark matter (DM) halo with unprecedented resolution. Our set of cosmological hydrodynamic Adaptive Mesh Refinement (AMR) simulations, called the Nut suite, allows us to investigate the effect of supernova feedback and UV photoionisation at high redshift with sub-parsec resolution. We subsequently follow the effect of interactions with the Milky Way-like halo using a lower spatial resolution (50pc) version of the simulation down to z=0. This latter produces a population of simulated satellites that we compare to the observed satellites of the Milky Way and M31. We find that supernova feedback reduces star formation in the least massive satellites but enhances it in the more massive ones. Photoionisation appears to play a very minor role in suppressing star and galaxy formation in all progenitors of satellite halos. By far the largest effect on the satellite population is found to be the mass of the host and whether gas cooling is included in the simulation or not. Indeed, inclusion of gas cooling dramatically reduces the number of satellites captured at high redshift which survive down to z=0.

preprint2012arXiv

The Epoch of Disk Settling: z~1 to Now

We present evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the internal kinematics of blue galaxies with stellar masses ranging 8.0 < log M* (M_Sun) < 10.7 over 0.2<z<1.2. DEEP2 provides galaxy spectra and Hubble imaging from which we measure emission-line kinematics and galaxy inclinations, respectively. Our large sample allows us to overcome scatter intrinsic to galaxy properties in order to examine trends in kinematics. We find that at a fixed stellar mass galaxies systematically decrease in disordered motions and increase in rotation velocity and potential well depth with time. Massive galaxies are the most well-ordered at all times examined, with higher rotation velocities and less disordered motions than less massive galaxies. We quantify disordered motions with an integrated gas velocity dispersion corrected for beam smearing (sigma_g). It is unlike the typical pressure-supported velocity dispersion measured for early type galaxies and galaxy bulges. Because both seeing and the width of our spectral slits comprise a significant fraction of the galaxy sizes, sigma_g integrates over velocity gradients on large scales which can correspond to non-ordered gas kinematics. We compile measurements of galaxy kinematics from the literature over 1.2<z<3.8 and do not find any trends with redshift, likely for the most part because these datasets are biased toward the most highly star-forming systems. In summary, over the last ~8 billion years since z=1.2, blue galaxies evolve from disordered to ordered systems as they settle to become the rotation-dominated disk galaxies observed in the Universe today, with the most massive galaxies being the most evolved at any time.

preprint2012arXiv

The Radius of Baryonic Collapse in Disc Galaxy Formation

In the standard picture of disc galaxy formation, baryons and dark matter receive the same tidal torques, and therefore approximately the same initial specific angular momentum. However, observations indicate that disc galaxies typically have only about half as much specific angular momentum as their dark matter haloes. We argue this does not necessarily imply that baryons lose this much specific angular momentum as they form galaxies. It may instead indicate that galaxies are most directly related to the inner regions of their host haloes, as may be expected in a scenario where baryons in the inner parts of haloes collapse first. A limiting case is examined under the idealised assumption of perfect angular momentum conservation. Namely, we determine the density contrast Delta, with respect to the critical density of the Universe, by which dark matter haloes need to be defined in order to have the same average specific angular momentum as the galaxies they host. Under the assumption that galaxies are related to haloes via their characteristic rotation velocities, the necessary Delta is ~600. This Delta corresponds to an average halo radius and mass which are ~60% and ~75%, respectively, of the virial values (i.e., for Delta = 200). We refer to this radius as the radius of baryonic collapse R_BC, since if specific angular momentum is conserved perfectly, baryons would come from within it. It is not likely a simple step function due to the complex gastrophysics involved, therefore we regard it as an effective radius. In summary, the difference between the predicted initial and the observed final specific angular momentum of galaxies, which is conventionally attributed solely to angular momentum loss, can more naturally be explained by a preference for collapse of baryons within R_BC, with possibly some later angular momentum transfer.

preprint2011arXiv

AGN feedback using AMR cosmological simulations

Feedback processes are thought to solve some of the long-standing issues of the numerical modelling of galaxy formation: over-cooling, low angular momentum, massive blue galaxies, extra-galactic enrichment, etc. The accretion of gas onto super-massive black holes in the centre of massive galaxies can release tremendous amounts of energy to the surrounding medium. We show, with cosmological Adaptive Mesh Refinement simulations, how the growth of black holes is regulated by the feedback from Active Galactic Nuclei using a new dual jet/heating mechanism. We discuss how this large amount of feedback is able to modify the cold baryon content of galaxies, and perturb the properties of the hot plasma in their vicinity.

preprint2011arXiv

Are cold flows detectable with metal absorption lines?

[Abridged] Cold gas flowing within the "cosmic web" is believed to be an important source of fuel for star formation at high redshift. However, the presence of such filamentary gas has never been observationally confirmed. In this work, we investigate in detail whether such cold gas is detectable using low-ionisation metal absorption lines, such as CII \lambda1334 as this technique has a proven observational record for detecting gaseous structures. Using a large statistical sample of galaxies from the Mare Nostrum N-body+AMR cosmological simulation, we find that the typical covering fraction of the dense, cold gas in 10^12 Msun haloes at z~2.5 is lower than expected (~5%). In addition, the absorption signal by the interstellar medium of the galaxy itself turns out to be so deep and so broad in velocity space that it completely drowns that of the filamentary gas. A detectable signal might be obtained from a cold filament exactly aligned with the line of sight, but this configuration is so unlikely that it would require surveying an overwhelmingly large number of candidate galaxies to tease it out. Finally, the predicted metallicity of the cold gas in filaments is extremely low (\leq 0.001 Zsun). Should this result persist when higher resolution runs are performed, it would significantly increase the difficulty of detecting filamentary gas inflows using metal lines. However, even if we assume that filaments are enriched to Zsun, the absorption signal that we compute is still weak. We are therefore led to conclude that it is extremely difficult to observationally prove or disprove the presence of cold filaments as the favorite accretion mode of galaxies using low-ionisation metal absorption lines. The Ly-alpha emission route looks more promising but due to the resonant nature of the line, radiative transfer simulations are required to fully characterize the observed signal.

preprint2011arXiv

How AGN feedback and metal cooling shape cluster entropy profiles

Observed clusters of galaxies essentially come in two flavors: non cool core clusters characterized by an isothermal temperature profile and a central entropy floor, and cool-core clusters where temperature and entropy in the central region are increasing with radius. Using cosmological resimulations of a galaxy cluster, we study the evolution of its intracluster medium (ICM) gas properties, and through them we assess the effect of different (sub-grid) modelling of the physical processes at play, namely gas cooling, star formation, feedback from supernovae and active galactic nuclei (AGN). More specifically we show that AGN feedback plays a major role in the pre-heating of the proto-cluster as it prevents a high concentration of mass from collecting in the center of the future galaxy cluster at early times. However, AGN activity during the cluster's later evolution is also required to regulate the mass flow into its core and prevent runaway star formation in the central galaxy. Whereas the energy deposited by supernovae alone is insufficient to prevent an overcooling catastrophe, supernovae are responsible for spreading a large amount of metals at high redshift, enhancing the cooling efficiency of the ICM gas. As the AGN energy release depends on the accretion rate of gas onto its central black hole engine, the AGN responds to this supernova enhanced gas accretion by injecting more energy into the surrounding gas, and as a result increases the amount of early pre-heating. We demonstrate that the interaction between an AGN jet and the ICM gas that regulates the growth of the AGN's BH, can naturally produce cool core clusters if we neglect metals. However, as soon as metals are allowed to contribute to the radiative cooling, only the non cool core solution is produced.

preprint2011arXiv

How Does Feedback Affect Milky Way Satellite Formation?

We use sub-parsec resolution hydrodynamic resimulations of a Milky Way (MW) like galaxy at high redshift to investigate the formation of the MW satellite galaxies. More specifically, we assess the impact of supernova feedback on the dwarf progenitors of these satellite, and the efficiency of a simple instantaneous reionisation scenario in suppressing star formation at the low-mass end of this dwarf distribution. Identifying galaxies in our high redshift simulation and tracking them to z=0 using a dark matter halo merger tree, we compare our results to present-day observations and determine the epoch at which we deem satellite galaxy formation must be completed. We find that only the low-mass end of the population of luminous subhalos of the Milky-Way like galaxy is not complete before redshift 8, and that although supernovae feedback reduces the stellar mass of the low-mass subhalos (log(M/Msolar) < 9), the number of surviving satellites around the Milky-Way like galaxy at z = 0 is the same in the run with or without supernova feedback. If a luminous halo is able to avoid accretion by the Milky-Way progenitor before redshift 3, then it is likely to survive as a MW satellite to redshift 0.

preprint2011arXiv

Most massive halos with Gumbel Statistics

We present an analytical calculation of the extreme value statistics for dark matter halos - that is, the probability distribution of the most massive halo within some region of the universe of specified shape and size. Our calculation makes use of the counts-in-cells formalism for the correlation functions, and the halo bias derived from the Sheth-Tormen mass function. We demonstrate the power of the method on spherical regions, comparing the results to measurements in a large cosmological dark matter simulation and achieving good agreement. Particularly good fits are obtained for the most likely value of the maximum mass and for the high-mass tail of the distribution, relevant in constraining cosmologies by observations of most massive clusters.

preprint2011arXiv

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.

preprint2011arXiv

The angular momentum of baryons and dark matter halos revisited

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense gas filaments, suggesting angular momentum transport by gas differs from that by dark matter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with adaptive mesh refinement, we find that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum, but that it is systematically higher than that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams into the central region of the halo, directly depositing this large amount of angular momentum within a sphere of radius r=0.1rvir. In contrast, dark matter particles pass through the central region unscathed, and a fraction of them ends up populating the outer regions of the halo (r/rvir>0.1), redistributing angular momentum in the process. As a result, large-scale motions of the cosmic web have to be considered as the origin of gas angular momentum rather than its virialised dark matter halo host. This generic result holds for halos of all masses at all redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at the centre of the halos. Despite this injection of angular momentum enriched gas, we predict an amount for stellar discs which is in fair agreement with observations at z=0. This arises because the total specific angular momentum of the baryons remains close to that of dark matter halos. We propose a new scenario where gas efficiently carries the angular momentum generated by large-scale structure motions deep inside dark matter halos, redistributing it only in the vicinity of the disc.

preprint2011arXiv

The environment and redshift dependence of accretion onto dark matter halos and subhalos

A dark-matter-only Horizon Project simulation is used to investigate the environment- and redshift- dependence of accretion onto both halos and subhalos. These objects grow in the simulation via mergers and via accretion of diffuse non-halo material, and we measure the combined signal from these two modes of accretion. It is found that the halo accretion rate varies less strongly with redshift than predicted by the Extended Press-Schechter (EPS) formalism and is dominated by minor-merger and diffuse accretion events at z=0, for all halos. These latter growth mechanisms may be able to drive the radio-mode feedback hypothesised for recent galaxy-formation models, and have both the correct accretion rate and form of cosmological evolution. The low redshift subhalo accretors in the simulation form a mass-selected subsample safely above the mass resolution limit that reside in the outer regions of their host, with ~70% beyond their host's virial radius, where they are probably not being significantly stripped of mass. These subhalos accrete, on average, at higher rates than halos at low redshift and we argue that this is due to their enhanced clustering at small scales. At cluster scales, the mass accretion rate onto halos and subhalos at low redshift is found to be only weakly dependent on environment and we confirm that at z~2 halos accrete independently of their environment at all scales, as reported by other authors. By comparing our results with an observational study of black hole growth, we support previous suggestions that at z>1, dark matter halos and their associated central black holes grew coevally, but show that by the present day, dark matter halos could be accreting at fractional rates that are up to a factor 3-4 higher than their associated black holes.

preprint2011arXiv

The hierarchical build-up of the Tully-Fisher relation

We use the semi-analytic model GalICS to predict the Tully-Fisher relation in the B, I and for the first time, in the K band, and its evolution with redshift, up to z~1. We refined the determination of the disk galaxies rotation velocity, with a dynamical recipe for the rotation curve, rather than a simple conversion from the total mass to maximum velocity. The new recipe takes into account the disk shape factor, and the angular momentum transfer occurring during secular evolution leading to the formation of bulges. This produces model rotation velocities that are lower by ~20-25% for the majority of the spirals. We implemented stellar population models with a complete treatment of the TP-AGB, which leads to a revision of the mass-to-light ratio in the near-IR. I/K band luminosities increase by ~0.3/0.5 mags at redshift z=0 and by ~0.5/1 mags at z=3. With these two new recipes in place, the comparison between the predicted Tully-Fisher relation with a series of datasets in the optical and near-IR, at redshifts between 0 and 1, is used as a diagnostics of the assembly and evolution of spiral galaxies in the model. At 0.4<z<1.2 the match between the new model and data is remarkably good, especially for later-type spirals (Sb/Sc). At z=0 the new model shows a net improvement in comparison with its original version of 2003, and in accord with recent observations in the K band, the model Tully-Fisher also shows a morphological differentiation. However, in all bands the z=0 model Tully-Fisher is too bright. We argue that this behaviour is caused by inadequate star formation histories in the model galaxies at low redshifts. The star-formation rate declines too slowly, due to continuous gas infall that is not efficiently suppressed. An analysis of the model disk scale lengths, at odds with observations, hints to some missing physics in the modeling of disk formation inside dark matter halos.

preprint2011arXiv

The impact of ISM turbulence, clustered star formation and feedback on galaxy mass assembly through cold flows and mergers

Two of the dominant channels for galaxy mass assembly are cold flows (cold gas supplied via the filaments of the cosmic web) and mergers. How these processes combine in a cosmological setting, at both low and high redshift, to produce the whole zoo of galaxies we observe is largely unknown. Indeed there is still much to understand about the detailed physics of each process in isolation. While these formation channels have been studied using hydrodynamical simulations, here we study their impact on gas properties and star formation (SF) with some of the first simulations that capture the multiphase, cloudy nature of the interstellar medium (ISM), by virtue of their high spatial resolution (and corresponding low temperature threshold). In this regime, we examine the competition between cold flows and a supernovae(SNe)-driven outflow in a very high-redshift galaxy (z {\approx} 9) and study the evolution of equal-mass galaxy mergers at low and high redshift, focusing on the induced SF. We find that SNe-driven outflows cannot reduce the cold accretion at z {\approx} 9 and that SF is actually enhanced due to the ensuing metal enrichment. We demonstrate how several recent observational results on galaxy populations (e.g. enhanced HCN/CO ratios in ULIRGs, a separate Kennicutt Schmidt (KS) sequence for starbursts and the population of compact early type galaxies (ETGs) at high redshift) can be explained with mechanisms captured in galaxy merger simulations, provided that the multiphase nature of the ISM is resolved.

preprint2011arXiv

The impact of supernovae driven winds on stream-fed protogalaxies

SNe driven winds are widely thought to be very influential in the high-redshift Universe, shaping the properties of the circum-galactic medium, enriching the IGM with metals and driving the evolution of low-mass galaxies. However, it is not yet fully understood how SNe driven winds interact with their surroundings in a cosmological context, nor is it clear whether they are able to significantly impact the evolution of low-mass galaxies from which they originate by altering the amount of cold material these accrete from the cosmic web. We implement a standard Taylor-Sedov type solution, widely used in the community to depict the combined action of many SN explosions, in a cosmological resimulation of a low mass galaxy at z =9 from the 'Nut' suite. However, in contrast with previous work, we achieve a resolution high enough to capture individual SN remnants in the Taylor-Sedov phase, for which the solution provides an accurate description of the expansion. We report the development of a high-velocity, far-reaching galactic wind produced by the combined action of SNe in the main galaxy and its satellites, which are located in the same or a neighbouring dark matter halo. Despite this, we find that (i) this wind carries out very little mass (the measured outflow is of the order of a tenth of the inflow/star formation rate) and (ii) the cold gas inflow rate remains essentially unchanged from the run without SNe feedback. Moreover, there are epochs during which star formation is enhanced in the feedback run relative to its radiative cooling only counterpart. We attribute this 'positive' feedback to the metal enrichment that is present only in the former. We conclude that at very high redshift, efficient SNe feedback can drive large-scale galactic winds but does not prevent massive cold gas inflow from fuelling galaxies, resulting in long-lived episodes of intense star formation.(abridged)

preprint2010arXiv

Galactic star formation in parsec-scale resolution simulations

The interstellar medium (ISM) in galaxies is multiphase and cloudy, with stars forming in the very dense, cold gas found in Giant Molecular Clouds (GMCs). Simulating the evolution of an entire galaxy, however, is a computational problem which covers many orders of magnitude, so many simulations cannot reach densities high enough or temperatures low enough to resolve this multiphase nature. Therefore, the formation of GMCs is not captured and the resulting gas distribution is smooth, contrary to observations. We investigate how star formation (SF) proceeds in simulated galaxies when we obtain parsec-scale resolution and more successfully capture the multiphase ISM. Both major mergers and the accretion of cold gas via filaments are dominant contributors to a galaxy's total stellar budget and we examine SF at high resolution in both of these contexts.

preprint2010arXiv

Identifying the progenitor set of present-day early-type galaxies: a view from the standard model

We present a comprehensive theoretical study, using a semi-analytical model within the standard LCDM framework, of the photometric properties of the progenitors of present-day early-type galaxies in the redshift range 0<z<1. We explore progenitors of all morphologies and study their characteristics as a function of the luminosity and local environment of the early-type remnant at z=0. In agreement with previous studies, we find that, while larger early-types are generally assembled later, their luminosity-weighted stellar ages are typically older. In dense cluster-like environments, 70% of early-type systems are `in place' by z=1 and evolve without interactions thereafter, while in the field the corresponding value is 30%. Averaging across all environments at z~1, less than 50% of the stellar mass which ends up in early-types today is actually in early-type progenitors at this redshift, in agreement with recent observational work. We develop probabilistic prescriptions which provide a means of including spiral (i.e. non early-type) progenitors at intermediate and high redshifts, based on their luminosity and optical colours. For example, we find that, at intermediate redshifts (z~0.5), large (M_V<-21.5), red (B-V>0.7) spirals have 75-95% chance of being an early-type progenitor, while the corresponding probability for large blue spirals (M_B<-21.5, B-V<0.7) is 50-75%. The prescriptions developed here can be used to address, from the perspective of the standard model, the issue of `progenitor bias', whereby the exclusion of late-type progenitors in observational studies can lead to inaccurate conclusions regarding the evolution of the early-type population over cosmic time. (abridged)

preprint2010arXiv

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 onto Super Massive 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 BH growth and associated AGN jet feedback. BHs accrete gas efficiently at high redshift (z>2), significantly pre-heating proto-cluster halos. 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 onto 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 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 beta 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.

preprint2010arXiv

The Origin and Evolution of the Mass-Metallicity Relation using GalICS

The GalICS (Galaxies in Cosmological Simulations) semi-analytical model of hierar- chical galaxy formation is used to investigate the effects of different galactic properties, including star formation rate (SFR) and outflows, on the shape of the mass metallic- ity relation and to predict the relation for galaxies at redshift z=2.27 and z=3.54. Our version of GalICS has the chemical evolution implemented in great detail and is less heavily reliant on approximations such as instantaneous recycling. We vary the model parameters controlling both the efficiency and redshift dependence of the SFR as well as the efficiency of supernova feedback. We find that the factors controlling the SFR influence the relation significantly at all redshifts and require a strong redshift dependence, proportional to 1+z, in order to reproduce the observed relation at the low mass end. Indeed, at any redshift, the predicted relation flattens out at the high mass end resulting in a poorer agreement with observations in this regime. We also find that variation of the parameters associated with outflows has a minimal effect on the relation at high redshift but does serve to alter its shape in the more recent past. We thus conclude that the relation is one between SFR and mass and that outflows are only important in shaping the relation at late times. When the relation is stratified by SFR it is apparent that the predicted galaxies with increasing stellar masses have higher SFRs, supporting the view that galaxy downsizing is the origin of the relation. Attempting to reproduce the observed relation, we vary the parameters controlling the efficiency of star formation and its redshift dependence and compare the predicted relations with Erb et al. (2006) at z=2.27 and Maiolino et al. (2008) at z=3.54 in order to find the best-fitting parameters. (Abridged)

preprint2009arXiv

Hierarchical models of high redshift galaxies with thermally pulsing asymptotic giant branch stars: comparison with observations

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.

preprint2009arXiv

The impact of TP-AGB stars on hierarchical galaxy formation models

The spectro-photometric properties of galaxies in galaxy formation models are obtained by combining the predicted history of star formation and mass accretion with the physics of stellar evolution through stellar population models. In the recent literature, significant differences have emerged regarding the implementation of the Thermally-Pulsing Asymptotic Giant Branch phase of stellar evolution. The emission in the TP-AGB phase dominates the bolometric and near-IR spectrum of intermediate-age (~1 Gyr) stellar populations, hence it is crucial for the correct modeling of the galaxy luminosities and colours. In this paper for the first time, we incorporate a full prescription of the TP-AGB phase in a semi-analytic model of galaxy formation. We find that the inclusion of the TP-AGB in the model spectra dramatically alters the predicted colour-magnitude relation and its evolution with redshift. When the TP-AGB phase is active, the rest-frame V-K galaxy colours are redder by almost 2 magnitudes in the redshift range z~2-3 and by 1 magnitude at z~1. Very red colours are produced in disk galaxies, so that the V-K colour distributions of disk and spheroids are virtually undistinguishable at low redshifts. We also find that the galaxy K-band emission is more than 1 magnitude higher in the range z~1-3. This may alleviate the difficulties met by the hierarchical clustering scenario in predicting the red galaxy population at high redshifts. The comparison between simulations and observations have to be revisited in the light of our results.

preprint2009arXiv

The Skeleton: Connecting Large Scale Structures to Galaxy Formation

We report on two quantitative, morphological estimators of the filamentary structure of the Cosmic Web, the so-called global and local skeletons. The first, based on a global study of the matter density gradient flow, allows us to study the connectivity between a density peak and its surroundings, with direct relevance to the anisotropic accretion via cold flows on galactic halos. From the second, based on a local constraint equation involving the derivatives of the field, we can derive predictions for powerful statistics, such as the differential length and the relative saddle to extrema counts of the Cosmic web as a function of density threshold (with application to percolation of structures and connectivity), as well as a theoretical framework to study their cosmic evolution through the onset of gravity-induced non-linearities.

Julien Devriendt

What is connected

Connect this record

See the researcher in context

Building this map preview

67 published item(s)

Catalogues of voids as antihalos in the local Universe

Cosmological simulations of the same spiral galaxy: the impact of baryonic physics

Impact of Radiation Feedback on the Formation of Globular Cluster Candidates during Cloud-Cloud Collisions

MgII in the JWST Era: a Probe of Lyman Continuum Escape?

On the Viability of Determining Galaxy Properties from Observations I: Star Formation Rates and Kinematics

Simulating Jellyfish Galaxies: A Case Study for a Gas-Rich Dwarf Galaxy

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

The Obelisk simulation: galaxies contribute more than AGN to HI reionization of protoclusters

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Momentum deposition of supernovae with cosmic rays

Predicting the Observability of Population III Stars with ELT-HARMONI via the Helium $1640{\rmÅ}$ emission line

Rivers of Gas I.: Unveiling The Properties of High Redshift Filaments

Black hole mergers from dwarf to massive galaxies with the NewHorizon and Horizon-AGN simulations

Detecting the Cosmic Web: Lyα Emission from Simulated Filaments at z=3

GalICS 2.1: a new semianalytic model for cold accretion, cooling, feedback and their roles in galaxy formation

Horizon-AGN virtual observatory -- 2: Template-free estimates of galaxy properties from colours

How primordial magnetic fields shrink galaxies

How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

New Methods for Identifying Lyman Continuum Leakers and Reionization-Epoch Analogues

The impact of AGN feedback on galaxy intrinsic alignments in the Horizon simulations

Why do extremely massive disc galaxies exist today?

The impact of the connectivity of the cosmic web on the physical properties of galaxies at its nodes

When galaxies align: intrinsic alignments of the progenitors of elliptical galaxies in the Horizon-AGN simulation

Bursty star formation feedback and cooling outflows

Caught in the rhythm: how satellites settle into a plane around their central galaxy

Comparing Simulations of AGN Feedback

Modeling Lyman-α Forest Cross-Correlations with LyMAS

Redshift and luminosity evolution of the intrinsic alignments of galaxies in Horizon-AGN

The cosmic evolution of massive black holes in the Horizon-AGN simulation

The Horizon-AGN simulation: morphological diversity of galaxies promoted by AGN feedback

A Detailed Study of Feedback from a Massive Star

Black hole evolution: I. Supernova-regulated black hole growth

Intrinsic alignment of simulated galaxies in the cosmic web: implications for weak lensing surveys

Intrinsic alignments of galaxies in the Horizon-AGN cosmological hydrodynamical simulation

Milking the spherical cow: on aspherical dynamics in spherical coordinates

Towards simulating star formation in turbulent high-z galaxies with mechanical supernova feedback

Black hole evolution: II. Spinning black holes in a supernova-driven turbulent interstellar medium

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

Mergers drive spin swings along the cosmic web

Swirling around filaments: are large-scale structure vortices spinning up dark halos?

Angular momentum transfer to a Milky Way disk at high redshift

Blowing cold flows away: the impact of early AGN activity on the formation of a brightest cluster galaxy progenitor

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

Constraining stellar assembly and AGN feedback at the peak epoch of star formation

Feeding compact bulges and supermassive black holes with low angular-momentum cosmic gas at high redshift

Lyman-alpha emission properties of simulated galaxies: interstellar medium structure and inclination effects

Satellite Survival in Highly Resolved Milky Way Class Halos

The Epoch of Disk Settling: z~1 to Now

The Radius of Baryonic Collapse in Disc Galaxy Formation

AGN feedback using AMR cosmological simulations

Are cold flows detectable with metal absorption lines?

How AGN feedback and metal cooling shape cluster entropy profiles

How Does Feedback Affect Milky Way Satellite Formation?

Most massive halos with Gumbel Statistics

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

The angular momentum of baryons and dark matter halos revisited

The environment and redshift dependence of accretion onto dark matter halos and subhalos

The hierarchical build-up of the Tully-Fisher relation

The impact of ISM turbulence, clustered star formation and feedback on galaxy mass assembly through cold flows and mergers

The impact of supernovae driven winds on stream-fed protogalaxies

Galactic star formation in parsec-scale resolution simulations

Identifying the progenitor set of present-day early-type galaxies: a view from the standard model

Jet-regulated cooling catastrophe

The Origin and Evolution of the Mass-Metallicity Relation using GalICS

Hierarchical models of high redshift galaxies with thermally pulsing asymptotic giant branch stars: comparison with observations

The impact of TP-AGB stars on hierarchical galaxy formation models

The Skeleton: Connecting Large Scale Structures to Galaxy Formation