skip to main content

This content will become publicly available on May 1, 2023

Title: Coevolution of Brightest Cluster Galaxies and Their Host Clusters in IllustrisTNG
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.
Authors:
; ; ;
Award ID(s):
2107724
Publication Date:
NSF-PAR ID:
10350030
Journal Name:
The Astrophysical Journal
Volume:
931
Issue:
1
Page Range or eLocation-ID:
31
ISSN:
0004-637X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass (M500 ∼ 1013–1015 M⊙) and the redshift range 0 ≲ z ≲ 1. Our joint analysis is based on the combination of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and Sloan Digital Sky Survey multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass–total mass (M500) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass–total mass relation is shallower than unity; $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of M500 agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsicmore »covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution.

    « less
  2. Abstract We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations of the full 2500 deg 2 South Pole Telescope (SPT)-Sunyaev–Zel’dovich cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-IR color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-hosting BCG fraction, which is defined as the fraction of BCGs hosting bright central AGNs over all possible BCGs. Assuming an evolving single-burst stellar population model, we find statistically significant evidence (>99.9%) for a mid-IR excess at high redshift compared to low redshift, suggesting that the fraction of AGN-hosting BCGs increases with redshift over the range of 0 < z < 1.3. The best-fit redshift trend of the AGN-hosting BCG fraction has the form (1 + z ) 4.1±1.0 . These results are consistent with previous studies in galaxy clusters as well as as in field galaxies. One way to explain this result is that member galaxies at high redshift tend to have more cold gas. While BCGs in nearby galaxy clusters grow mostly by dry mergers with cluster members, leading to no increase in AGN activity, BCGsmore »at high redshift could primarily merge with gas-rich satellites, providing fuel for feeding AGNs. If this observed increase in AGN activity is linked to gas-rich mergers rather than ICM cooling, we would expect to see an increase in scatter in the P cav versus L cool relation at z > 1. Last, this work confirms that the runaway cooling phase, as predicted by the classical cooling-flow model, in the Phoenix cluster is extremely rare and most BCGs have low (relative to Eddington) black hole accretion rates.« less
  3. ABSTRACT

    We characterize the properties and evolution of bright central galaxies (BCGs) and the surrounding intracluster light (ICL) in galaxy clusters identified in the Dark Energy Survey and Atacama Cosmology Telescope Survey (DES-ACT) overlapping regions, covering the redshift range 0.20 < z < 0.80. Over this redshift range, we measure no change in the ICL’s stellar content (between 50 and 300 kpc) in clusters with log10(M200m,SZ/M⊙) >14.4. We also measure the stellar mass–halo mass (SMHM) relation for the BCG+ICL system and find that the slope, β, which characterizes the dependence of M200m,SZ on the BCG+ICL stellar mass, increases with radius. The outskirts are more strongly correlated with the halo than the core, which supports that the BCG+ICL system follows a two-phase growth, where recent growth (z < 2) occurs beyond the BCG’s core. Additionally, we compare our observed SMHM relation results to the IllustrisTNG300-1 cosmological hydrodynamic simulations and find moderate qualitative agreement in the amount of diffuse light. However, the SMHM relation’s slope is steeper in TNG300-1 and the intrinsic scatter is lower, likely from the absence of projection effects in TNG300-1. Additionally, we find that the ICL exhibits a colour gradient such that the outskirts are bluer than the core.more »Moreover, for the lower halo mass clusters (log10(M200m,SZ/M⊙) < 14.59), we detect a modest change in the colour gradient’s slope with lookback time, which combined with the absence of stellar mass growth may suggest that lower mass clusters have been involved in growth via tidal stripping more recently than their higher mass counterparts.

    « less
  4. ABSTRACT We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at z = 0 with stellar masses $10^{10.7}\, \mathrm{M}_{\odot } \leqslant M_{\ast } \leqslant 10^{11.9}\, \mathrm{M}_{\odot }$, the total power-law slope has a mean of 〈γ′〉 = 2.011 ± 0.007 and a scatter of $\sigma _{\gamma ^{\prime }} = 0.171$ over the radial range 0.4–4 times the stellar half-mass radius. Several correlations between γ′ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction, and in situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, γ′ is almost constant with redshift below z = 2. Through analysing IllustrisTNG model variations, we show that black hole kinetic winds are crucial to lowering γ′ and matching observed galaxy correlations. The effects of stellar winds on γ′ are subdominant compared to active galactic nucleus (AGN) feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter haloes are well described by steeper than NFW profiles, and they are steeper in the full physics (FP) run than their counterpartsmore »in the dark matter-only (DMO) run. Their inner density slopes anticorrelate (remain constant) with the halo mass in the FP (DMO) run, and anticorrelate with the halo concentration parameter c200 in both the types of runs. The dark matter haloes of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.« less
  5. ABSTRACT

    Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, bv – defined as the ratio of galaxy and DM velocity dispersions – as a function of redshift, host halo mass, and galaxy stellar mass threshold ($M_{\rm \star , sat}$), for massive haloes ($M_{\rm 200c}\gt 10^{13.5} \, {\rm M}_\odot$) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of bv increasing with M200c and decreasing with redshift and $M_{\rm \star , sat}$. The ensemble-estimated theoretical uncertainty in bv is 2–3 per cent, but becomes percent-level when considering only the three highest resolution simulations. We update the mass–richness normalization for an SDSSmore »redMaPPer cluster sample, and find our improved bv estimates reduce the normalization uncertainty from 22 to 8 per cent, demonstrating that dynamical mass estimation is competitive with weak lensing mass estimation. We discuss necessary steps for further improving this precision. Our estimates for $b_v(M_{\rm 200c}, M_{\rm \star , sat}, z)$ are made publicly available.

    « less