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 on to 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 haloes 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.
The recently introduced discrete persistent structure extractor (DisPerSE, Sousbie, Paper I) is implemented on realistic 3D cosmological simulations and observed redshift catalogues; it is found that ...DisPerSE traces very well the observed filaments, walls and voids seen both in simulations and in observations. In either setting, filaments are shown to connect on to haloes, outskirt walls, which circumvent voids, as is topologically required by the Morse theory. Indeed this algorithm returns the optimal critical set while operating directly on the particles. DisPerSE, as illustrated here, assumes nothing about the geometry of the survey or its homogeneity, and yields a natural (topologically motivated) self-consistent criterion for selecting the significance level of the identified structures. It is shown that this extraction is possible even for very sparsely sampled point processes, as a function of the persistence ratio (a measure of the significance of topological connections between critical points). Hence, astrophysicists should be in a position to trace precisely the locus of filaments, walls and voids from such samples and assess the confidence of the post-processed sets as a function of this threshold, which can be expressed relative to the expected amplitude of shot noise. In a cosmic framework, this criterion is shown to level with the friends-of-friends structure finder for the identification of peaks, while it also identifies the connected filaments and walls, and quantitatively recovers the full set of topological invariants (number of holes, etc.) directly from the particles, and at no extra cost as a function of the persistence threshold. This criterion is found to be sufficient even if one particle out of two is noise, when the persistence ratio is set to 3σ or more. The algorithm is also implemented on the SDSS catalogue and used to locate interesting configurations of the filamentary structure. In this context, we carried the identification of an 'optically faint' cluster at the intersection of filaments through the recent observation of its X-ray counterpart by Suzaku.
Accurate cosmology from upcoming weak lensing surveys relies on knowledge of the total matter power spectrum at per cent level at scales k < 10 h Mpc^−1, for which modelling the impact of baryonic ...physics is crucial. We compare measurements of the total matter power spectrum from the Horizon cosmological hydrodynamical simulations: a dark-matter-only run, one with full baryonic physics, and another lacking active galactic nucleus (AGN) feedback. Baryons cause a suppression of power at k ≃ 10 h Mpc^−1 of |${\lt}15{{\ \rm per\ cent}}$| at |$z$| = 0, and an enhancement of a factor of a few at smaller scales due to the more efficient cooling and star formation. The results are sensitive to the presence of the highest mass haloes in the simulation and the distribution of dark matter is also impacted up to a few per cent. The redshift evolution of the effect is non-monotonic throughout |$z$| = 0−5 due to an interplay between AGN feedback and gas pressure, and the growth of structure. We investigate the effectiveness of an analytic ‘baryonic correction model’ in describing our results. We require a different redshift evolution and propose an alternative fitting function with four free parameters that reproduces our results within |$5{{\ \rm per\ cent}}$|. Compared to other simulations, we find the impact of baryonic processes on the total matter power spectrum to be smaller at |$z$| = 0. Correspondingly, our results suggest that AGN feedback is not strong enough in the simulation. Total matter power spectra from the Horizon simulations are made publicly available at https://www.horizon-simulation.org/catalogues.html.
The physics of diffuse gas accretion and the properties of the cold and hot modes of accretion on to proto-galaxies between z= 2 and 5.4 is investigated using the large cosmological simulation ...performed with the ramses code on the MareNostrum supercomputing facility. Galactic winds, chemical enrichment, ultraviolet background heating and radiative cooling are taken into account in this very high resolution simulation. Using accretion-weighted temperature histograms, we have performed two different measurements of the thermal state of the gas accreted towards the central galaxy. The first measurement, performed using accretion-weighted histograms on a spherical surface of radius 0.2Rvir centred on the densest gas structure near the halo centre of mass, is a good indicator of the presence of an accretion shock in the vicinity of the galactic disc. We define the hot shock mass, Mshock, as the typical halo mass separating cold dominated from hot dominated accretion in the vicinity of the galaxy. The second measurement is performed by radially averaging histograms between 0.2Rvir and Rvir, in order to detect radially extended structures such as gas filaments: this is a good proxy for detecting cold streams feeding the central galaxy. We define Mstream as the transition mass separating cold dominated from hot dominated accretion in the outer halo, marking the disappearance of these cold streams. We find a hot shock transition mass of Mshock= 1011.6M⊙ (dark matter), with no significant evolution with redshift. Conversely, we find that Mstream increases sharply with z. Our measurements are in agreement with the analytical predictions of Birnboim & Dekel and Dekel & Birnboim, if we correct their model by assuming low metallicity(≤10−3Z⊙)for the filaments, correspondingly to our measurements. Metal enrichment of the intergalactic medium is therefore a key ingredient in determining the transition mass from cold to hot dominated diffuse gas accretion. We find that the diffuse cold gas supply at the inner halo stops at z= 2 for objects with stellar masses of about 1011.1M⊙, which is close to the quenching mass determined observationally by Bundy et al. However, its evolution with z is not well constrained, making it difficult to rule out or confirm the need for an additional feedback process such as active galactic nuclei.
Abstract
The mass, accretion rate, and formation time of dark matter haloes near protofilaments (identified as saddle points of the potential) are analytically predicted using a conditional version ...of the excursion set approach in its so-called upcrossing approximation. The model predicts that at fixed mass, mass accretion rate and formation time vary with orientation and distance from the saddle, demonstrating that assembly bias is indeed influenced by the tides imposed by the cosmic web. Starved, early-forming haloes of smaller mass lie preferentially along the main axis of filaments, while more massive and younger haloes are found closer to the nodes. Distinct gradients for distinct tracers such as typical mass and accretion rate occur because the saddle condition is anisotropic, and because the statistics of these observables depend on both the conditional means and their covariances. The theory is extended to other critical points of the potential field. The response of the mass function to variations of the matter density field (the so-called large-scale bias) is computed, and its trend with accretion rate is shown to invert along the filament. The signature of this model should correspond at low redshift to an excess of reddened galactic hosts at fixed mass along preferred directions, as recently reported in spectroscopic and photometric surveys and in hydrodynamical simulations. The anisotropy of the cosmic web emerges therefore as a significant ingredient to describe jointly the dynamics and physics of galaxies, e.g. in the context of intrinsic alignments or morphological diversity.
Our statistical understanding of galaxy evolution is fundamentally driven by objects that lie above the surface-brightness limits of current wide-area surveys (μ ∼ 23 mag arcsec−2). While both theory ...and small, deep surveys have hinted at a rich population of low-surface-brightness galaxies (LSBGs) fainter than these limits, their formation remains poorly understood. We use Horizon-AGN, a cosmological hydrodynamical simulation to study how LSBGs, and in particular the population of ultra-diffuse galaxies (UDGs; μ > 24.5 mag arcsec−2), form and evolve over time. For M∗>108M⊙, LSBGs contribute 47, 7, and 6 per cent of the local number, mass, and luminosity densities, respectively (∼85/11/10 per cent for M∗>107M⊙). Today’s LSBGs have similar dark-matter fractions and angular momenta to high-surface-brightness galaxies (HSBGs; μ < 23 mag arcsec−2), but larger effective radii (×2.5 for UDGs) and lower fractions of dense, star-forming gas (more than ×6 less in UDGs than HSBGs). LSBGs originate from the same progenitors as HSBGs at z > 2. However, LSBG progenitors form stars more rapidly at early epochs. The higher resultant rate of supernova-energy injection flattens their gas-density profiles, which, in turn, creates shallower stellar profiles that are more susceptible to tidal processes. After z ∼ 1, tidal perturbations broaden LSBG stellar distributions and heat their cold gas, creating the diffuse, largely gas-poor LSBGs seen today. In clusters, ram-pressure stripping provides an additional mechanism that assists in gas removal in LSBG progenitors. Our results offer insights into the formation of a galaxy population that is central to a complete understanding of galaxy evolution, and that will be a key topic of research using new and forthcoming deep-wide surveys.
Abstract
The variation of galaxy stellar masses and colour types with the distance to projected cosmic filaments are quantified using the precise photometric redshifts of the COSMOS2015 catalogue ...extracted from Cosmological Evolution Survey (COSMOS) field (2 deg2). Realistic mock catalogues are also extracted from the lightcone of the cosmological hydrodynamical simulation Horizon-AGN. They show that the photometric redshift accuracy of the observed catalogue (σz < 0.015 at M* > 1010M⊙ and z < 0.9) is sufficient to provide two-dimensional (2D) filaments that closely match their projected three-dimensional (3D) counterparts. Transverse stellar mass gradients are measured in projected slices of thickness 75 Mpc between 0.5 < z < 0.9, showing that the most massive galaxies are statistically closer to their neighbouring filament. At fixed stellar mass, passive galaxies are also found closer to their filament, while active star-forming galaxies statistically lie further away. The contributions of nodes and local density are removed from these gradients to highlight the specific role played by the geometry of the filaments. We find that the measured signal does persist after this removal, clearly demonstrating that proximity to a filament is not equivalent to proximity to an overdensity. These findings are in agreement with gradients measured in both 2D and 3D in the Horizon-AGN simulation and those observed in the spectroscopic surveys VIPERS and GAMA (which both rely on the identification of 3D filaments). They are consistent with a picture in which the influence of the geometry of the large-scale environment drives anisotropic tides that impact the assembly history of galaxies, and hence their observed properties.
Abstract
The role of the cosmic web in shaping galaxy properties is investigated in the Galaxy And Mass Assembly (GAMA) spectroscopic survey in the redshift range 0.03 ≤ z ≤ 0.25. The stellar mass, ...u − r dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. These trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar mass and large-scale density, the traceless component of the tides from the anisotropic large-scale environment also shapes galactic properties. An extension of excursion theory accounting for filamentary tides provides a qualitative explanation in terms of anisotropic assembly bias: at a given mass, the accretion rate varies with the orientation and distance to filaments. It also explains the absence of type/colour gradients in the data on smaller, non-linear scales.
The intrinsic alignments of galaxies are recognized as a contaminant to weak gravitational lensing measurements. In this work, we study the alignment of galaxy shapes and spins at low redshift (z ∼ ...0.5) in Horizon-AGN, an adaptive-mesh-refinement hydrodynamical cosmological simulation box of 100 h
− 1 Mpc a side with AGN feedback implementation. We find that spheroidal galaxies in the simulation show a tendency to be aligned radially towards overdensities 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.
We present a first principle approach to obtain analytical predictions for spherically averaged cosmic densities in the mildly non-linear regime that go well beyond what is usually achieved by ...standard perturbation theory. A large deviation principle allows us to compute the leading order cumulants of average densities in concentric cells. In this symmetry, the spherical collapse model leads to cumulant generating functions that are robust for finite variances and free of critical points when logarithmic density transformations are implemented. They yield in turn accurate density probability distribution functions (PDFs) from a straightforward saddle-point approximation valid for all density values. Based on this easy-to-implement modification, explicit analytic formulas for the evaluation of the one- and two-cell PDF are provided. The theoretical predictions obtained for the PDFs are accurate to a few per cent compared to the numerical integration, regardless of the density under consideration and in excellent agreement with N-body simulations for a wide range of densities. This formalism should prove valuable for accurately probing the quasi-linear scales of low-redshift surveys for arbitrary primordial power spectra.