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: Concordance between Observations and Simulations in the Evolution of the Mass Relation between Supermassive Black Holes and Their Host Galaxies
Abstract We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at 0.2 < z < 1.7 using well-matched observations and multiple state-of-the-art simulations (e.g., MassiveBlackII, Horizon-AGN, Illustris, TNG, and a semianalytic model). The observed sample consists of 646 uniformly selected Sloan Digital Sky Survey quasars (0.2 < z < 0.8) and 32 broad-line active galactic nuclei (AGNs; 1.2 < z < 1.7) with imaging from Hyper Suprime-Cam (HSC) for the former and Hubble Space Telescope (HST) for the latter. We first add realistic observational uncertainties to the simulation data and then construct a simulated sample in the same manner as the observations. Over the full redshift range, our analysis demonstrates that all simulations predict a level of intrinsic scatter of the scaling relations comparable to the observations that appear to agree with the dispersion of the local relation. Regarding the mean relation, Horizon-AGN and TNG are in closest agreement with the observations at low and high redshift ( z ∼ 0.2 and 1.5, respectively), while the other simulations show subtle differences within the uncertainties. For insight into the physics involved, the scatter of the scaling relation, seen in the SAM, is reduced by a factor of two and closer to the observations after adopting a new feedback model that considers the geometry of the AGN outflow. The consistency in the dispersion with redshift in our analysis supports the importance of both quasar- and radio-mode feedback prescriptions in the simulations. Finally, we highlight the importance of increasing the sensitivity (e.g., using the James Webb Space Telescope), thereby pushing to lower masses and minimizing biases due to selection effects.  more » « less
Award ID(s):
1909933
PAR ID:
10351054
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
933
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
132
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We use IllustrisTNG simulations to explore the dynamic scaling relation between massive clusters and their—central—brightest cluster galaxies (BCGs). The IllustrisTNG-300 simulation we use includes 280 massive clusters from the z = 0 snapshot with M 200 > 10 14 M ⊙ , enabling a robust statistical analysis. We derive the line-of-sight velocity dispersion of the stellar particles of the BCGs ( σ *,BCG ), analogous to the observed BCG stellar velocity dispersion. We also compute the subhalo velocity dispersion to measure the cluster velocity dispersion ( σ cl ). Both σ *,BCG and σ cl are proportional to the cluster halo mass, but the slopes differ slightly. Thus, like the observed relation, σ *,BCG / σ cl declines as a function of σ cl , but the scatter is large. We explore the redshift evolution of the σ *,BCG − σ cl scaling relation for z ≲ 1 in a way that can be compared directly with observations. The scaling relation has a similar slope at high redshift, but the scatter increases because of the large scatter in σ *,BCG . The simulations imply that high-redshift BCGs are dynamically more complex than their low-redshift counterparts. 
    more » « less
  2. ABSTRACT The James Webb Space Telescope will have the power to characterize high-redshift quasars at z ≥ 6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e. with bolometric luminosity $$L_{\rm bol}\geqslant 10^{46}\, \rm erg/s$$) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with $$L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$$ opens a path to understand the build-up of more normal BHs at z ≥ 6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z ≥ 4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the MBH − M⋆ scaling relation at z = 0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z ≥ 4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with $$L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$$ observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z ∼ 6, with the same diversity of subgrid physics, in order to gain statistics on the most extreme objects at high redshift. 
    more » « less
  3. We present a Hubble Space Telescope (HST) weak gravitational lensing study of nine distant and massive galaxy clusters with redshifts 1.0 ≲  z  ≲ 1.7 ( z median  = 1.4) and Sunyaev Zel’dovich (SZ) detection significance ξ  > 6.0 from the South Pole Telescope Sunyaev Zel’dovich (SPT-SZ) survey. We measured weak lensing galaxy shapes in HST/ACS F 606 W and F 814 W images and used additional observations from HST/WFC3 in F 110 W and VLT/FORS2 in U HIGH to preferentially select background galaxies at z  ≳ 1.8, achieving a high purity. We combined recent redshift estimates from the CANDELS/3D-HST and HUDF fields to infer an improved estimate of the source redshift distribution. We measured weak lensing masses by fitting the tangential reduced shear profiles with spherical Navarro-Frenk-White (NFW) models. We obtained the largest lensing mass in our sample for the cluster SPT-CL J2040−4451, thereby confirming earlier results that suggest a high lensing mass of this cluster compared to X-ray and SZ mass measurements. Combining our weak lensing mass constraints with results obtained by previous studies for lower redshift clusters, we extended the calibration of the scaling relation between the unbiased SZ detection significance ζ and the cluster mass for the SPT-SZ survey out to higher redshifts. We found that the mass scale inferred from our highest redshift bin (1.2 <  z  < 1.7) is consistent with an extrapolation of constraints derived from lower redshifts, albeit with large statistical uncertainties. Thus, our results show a similar tendency as found in previous studies, where the cluster mass scale derived from the weak lensing data is lower than the mass scale expected in a Planckν ΛCDM (i.e. ν Λ cold dark matter) cosmology given the SPT-SZ cluster number counts. 
    more » « less
  4. We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 <  z <  0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posteriorP(M200|ζ̂, λ̂,z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂,z) observables. This method enables simultaneous constraints on richnessλ-mass and tSZE detection significanceζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of theλ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity. 
    more » « less
  5. Abstract We analyze the cooling and feedback properties of 48 galaxy clusters at redshifts 0.4 < z < 1.3 selected from the South Pole Telescope (SPT) catalogs to evolve like the progenitors of massive and well-studied systems at z ∼ 0. We estimate the radio power at the brightest cluster galaxy (BCG) location of each cluster from an analysis of Australia Telescope Compact Array data. Assuming that the scaling relation between the radio power and active galactic nucleus (AGN) cavity power P cav observed at low redshift does not evolve with redshift, we use these measurements in order to estimate the expected AGN cavity power in the core of each system. We estimate the X-ray luminosity within the cooling radius L cool of each cluster from a joint analysis of the available Chandra X-ray and SPT Sunyaev–Zel’dovich (SZ) data. This allows us to characterize the redshift evolution of the P cav / L cool ratio. When combined with low-redshift results, these constraints enable investigations of the properties of the feedback–cooling cycle across 9 Gyr of cluster growth. We model the redshift evolution of this ratio measured for cool-core clusters by a log-normal distribution Log -  ( α + β z , σ 2 ) and constrain the slope of the mean evolution to β = −0.05 ± 0.47. This analysis improves the constraints on the slope of this relation by a factor of two. We find no evidence of redshift evolution of the feedback–cooling equilibrium in these clusters, which suggests that the onset of radio-mode feedback took place at an early stage of cluster formation. High values of P cav / L cool are found at the BCG location of noncool-core clusters, which might suggest that the timescales of the AGN feedback cycle and the cool core–noncool core transition are different. This work demonstrates that the joint analysis of radio, SZ, and X-ray data solidifies the investigation of AGN feedback at high redshifts. 
    more » « less