skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Trinity I: self-consistently modelling the dark matter halo–galaxy–supermassive black hole connection from z  = 0–10
ABSTRACT We present Trinity, a flexible empirical model that self-consistently infers the statistical connection between dark matter haloes, galaxies, and supermassive black holes (SMBHs). Trinity is constrained by galaxy observables from 0 < z < 10 [galaxies’ stellar mass functions, specific and cosmic star formation rates (SFRs), quenched fractions, and UV luminosity functions] and SMBH observables from 0 < z < 6.5 (quasar luminosity functions, quasar probability distribution functions, active black hole mass functions, local SMBH mass–bulge mass relations, and the observed SMBH mass distributions of high-redshift bright quasars). The model includes full treatment of observational systematics [e.g. active galactic nucleus (AGN) obscuration and errors in stellar masses]. From these data, Trinity infers the average SMBH mass, SMBH accretion rate, merger rate, and Eddington ratio distribution as functions of halo mass, galaxy stellar mass, and redshift. Key findings include: (1) the normalization and the slope of the SMBH mass–bulge mass relation increases mildly from z = 0 to z = 10; (2) The best-fitting AGN radiative+kinetic efficiency is ∼0.05–0.06, but can be in the range ∼0.035–0.07 with alternative input assumptions; (3) AGNs show downsizing, i.e. the Eddington ratios of more massive SMBHs start to decrease earlier than those of lower mass objects; (4) The average ratio between average SMBH accretion rate and SFR is ∼10−3 for low-mass galaxies, which are primarily star-forming. This ratio increases to ∼10−1 for the most massive haloes below z ∼ 1, where star formation is quenched but SMBHs continue to accrete.  more » « less
Award ID(s):
1908284
PAR ID:
10492479
Author(s) / Creator(s):
 ; ; ; ; ; ;  ;
Publisher / Repository:
MNRAS
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
518
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2123 to 2163
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We present predictions for the high-redshift halo–galaxy–supermassive black hole (SMBH) connection from the Trinity model. Matching a comprehensive compilation of galaxy (0 ≤ z ≤ 13) and SMBH data sets (0 ≤ z ≤ 6.5), Trinity finds: (1) The number of SMBHs with M• > 109 M⊙ in the observable Universe increases by five orders of magnitude from z ∼ 10 to z ∼ 2, and by another factor of ∼3 from z ∼ 2 to z = 0; (2) The M• > 109 and 1010 M⊙ SMBHs at z ∼ 6 live in haloes with ∼(2 − 3) and (3 − 5) × 1012 M⊙; (3) the newly discovered JWST AGN candidates at 7 ≲ z ≲ 11 are overmassive compared to the intrinsic SMBH mass–galaxy mass relation from Trinity, but they are still broadly consistent with Trinity predictions for flux limited AGN samples with Lauer bias. This bias favours the detection for overmassive SMBHs due to higher luminosities at a fixed Eddington ratio. However UHZ1’s M•/M* ratio is still some 1 dex higher than Trinity AGNs, indicating a discrepancy; (4) Trinity underpredicts the number densities of GN-z11 and CEERS_1019 analogues. But given the strong constraints from existing data in Trinity, the extra constraint from GN-z11 and CEERS_1019 does not significantly change trinity model results. (5) z = 6–10 quasar luminosity functions will reduce uncertainties in the trinity prediction of the z = 6–10 SMBH mass–galaxy mass relation by up to ∼0.5 dex. These luminosity functions will be available with future telescopes, such as Roman and Euclid. 
    more » « less
  2. null (Ed.)
    ABSTRACT The James Webb Space Telescope (JWST) is expected to observe galaxies at z > 10 that are presently inaccessible. Here, we use a self-consistent empirical model, the universemachine, to generate mock galaxy catalogues and light-cones over the redshift range z = 0−15. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy–halo relationships, and galaxy–galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at z < 10. We predict that Cycle 1 JWST surveys will very likely detect galaxies with M* > 107 M⊙ and/or M1500 < −17 out to at least z ∼ 13.5. Number density uncertainties at z > 12 expand dramatically, so efforts to detect z > 12 galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits (M* < 107 M⊙ or M1500 > −17). For reionization models, extrapolating luminosity functions with a constant faint-end slope from M1500 = −17 down to M1500 = −12 gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to z ∼ 12. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and light-cones for common fields are available online. 
    more » « less
  3. ABSTRACT Recent observations indicate that galactic outflows are ubiquitous in high-redshift (high-z) galaxies, including normal star-forming galaxies, quasar hosts, and dusty star-forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star-formation histories and galactic outflow properties of galaxies in massive haloes ($$10^{12}\, {\rm M}_{\odot }\ \lt\ M_{\rm vir}\ \lt\ 5\times 10^{12}\, {\rm M}_{\odot }$$) at z ≳ 5.5 in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei. The simulated galaxies resemble z > 4 DSFGs, with star-formation rates of $$\sim\!{1000}\ {\rm M}_{\odot }\, \rm yr^{-1}$$ and molecular gas masses of Mmol ∼ 1010 M⊙. However, the simulated galaxies are characterized by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of ∼0.1, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium. We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of 30 per cent. Finally, compared with FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high-z sample. We argue that both effects originate from the higher total and gas surface densities that characterize high-z massive systems. 
    more » « less
  4. null (Ed.)
    ABSTRACT The galaxy size–stellar mass and central surface density–stellar mass relationships are fundamental observational constraints on galaxy formation models. However, inferring the physical size of a galaxy from observed stellar emission is non-trivial due to various observational effects, such as the mass-to-light ratio variations that can be caused by non-uniform stellar ages, metallicities, and dust attenuation. Consequently, forward-modelling light-based sizes from simulations is desirable. In this work, we use the skirt  dust radiative transfer code to generate synthetic observations of massive galaxies ($$M_{*}\sim 10^{11}\, \rm {M_{\odot }}$$ at z = 2, hosted by haloes of mass $$M_{\rm {halo}}\sim 10^{12.5}\, \rm {M_{\odot }}$$) from high-resolution cosmological zoom-in simulations that form part of the Feedback In Realistic Environments project. The simulations used in this paper include explicit stellar feedback but no active galactic nucleus (AGN) feedback. From each mock observation, we infer the effective radius (Re), as well as the stellar mass surface density within this radius and within $$1\, \rm {kpc}$$ (Σe and Σ1, respectively). We first investigate how well the intrinsic half-mass radius and stellar mass surface density can be inferred from observables. The majority of predicted sizes and surface densities are within a factor of 2 of the intrinsic values. We then compare our predictions to the observed size–mass relationship and the Σ1−M⋆ and Σe−M⋆ relationships. At z ≳ 2, the simulated massive galaxies are in general agreement with observational scaling relations. At z ≲ 2, they evolve to become too compact but still star forming, in the stellar mass and redshift regime where many of them should be quenched. Our results suggest that some additional source of feedback, such as AGN-driven outflows, is necessary in order to decrease the central densities of the simulated massive galaxies to bring them into agreement with observations at z ≲ 2. 
    more » « less
  5. ABSTRACT We measure the rate of environmentally driven star formation quenching in galaxies at z ∼ 1, using eleven massive ($$M\approx 2\times 10^{14}\, \mathrm{M}_\odot$$) galaxy clusters spanning a redshift range 1.0 < z < 1.4 from the GOGREEN sample. We identify three different types of transition galaxies: ‘green valley’ (GV) galaxies identified from their rest-frame (NUV − V) and (V − J) colours; ‘blue quiescent’ (BQ) galaxies, found at the blue end of the quiescent sequence in (U − V) and (V − J) colour; and spectroscopic post-starburst (PSB) galaxies. We measure the abundance of these galaxies as a function of stellar mass and environment. For high-stellar mass galaxies (log M/M⊙ > 10.5) we do not find any significant excess of transition galaxies in clusters, relative to a comparison field sample at the same redshift. It is likely that such galaxies were quenched prior to their accretion in the cluster, in group, filament, or protocluster environments. For lower stellar mass galaxies (9.5 < log M/M⊙ < 10.5) there is a small but significant excess of transition galaxies in clusters, accounting for an additional ∼5–10 per cent of the population compared with the field. We show that our data are consistent with a scenario in which 20–30 per cent of low-mass, star-forming galaxies in clusters are environmentally quenched every Gyr, and that this rate slowly declines from z = 1 to z = 0. While environmental quenching of these galaxies may include a long delay time during which star formation declines slowly, in most cases this must end with a rapid (τ < 1 Gyr) decline in star formation rate. 
    more » « less