skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, September 13 until 2:00 AM ET on Saturday, September 14 due to maintenance. We apologize for the inconvenience.


This content will become publicly available on March 29, 2025

Title: Mapping the Growth of Supermassive Black Holes as a Function of Galaxy Stellar Mass and Redshift
Abstract

The growth of supermassive black holes is strongly linked to their galaxies. It has been shown that the population mean black hole accretion rate (BHAR¯) primarily correlates with the galaxy stellar mass (M) and redshift for the general galaxy population. This work aims to provide the best measurements ofBHAR¯as a function ofMand redshift over ranges of 109.5<M< 1012Mandz< 4. We compile an unprecedentedly large sample with 8000 active galactic nuclei (AGNs) and 1.3 million normal galaxies from nine high-quality survey fields following a wedding cake design. We further develop a semiparametric Bayesian method that can reasonably estimateBHAR¯and the corresponding uncertainties, even for sparsely populated regions in the parameter space.BHAR¯is constrained by X-ray surveys sampling the AGN accretion power and UV-to-infrared multiwavelength surveys sampling the galaxy population. Our results can independently predict the X-ray luminosity function (XLF) from the galaxy stellar mass function (SMF), and the prediction is consistent with the observed XLF. We also try adding external constraints from the observed SMF and XLF. We further measureBHAR¯for star-forming and quiescent galaxies and show that star-formingBHAR¯is generally larger than or at least comparable to the quiescentBHAR¯.

 
more » « less
Award ID(s):
2106990
NSF-PAR ID:
10522735
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
964
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
183
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We investigate the effectiveness of the statistical radio frequency interference (RFI) mitigation technique spectral kurtosis (SK^) in the face of simulated realistic RFI signals.SK^estimates the kurtosis of a collection ofMpower values in a single channel and provides a detection metric that is able to discern between human-made RFI and incoherent astronomical signals of interest. We test the ability ofSK^to flag signals with various representative modulation types, data rates, duty cycles, and carrier frequencies. We flag with various accumulation lengthsMand implement multiscaleSK^, which combines information from adjacent time-frequency bins to mitigate weaknesses in single-scaleSK^. We find that signals with significant sidelobe emission from high data rates are harder to flag, as well as signals with a 50% effective duty cycle and weak signal-to-noise ratios. MultiscaleSK^with at least one extra channel can detect both the center channel and sideband interference, flagging greater than 90% as long as the bin channel width is wider in frequency than the RFI.

     
    more » « less
  2. Abstract

    We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (βp= 0.25) collisionless ion–electron shocks with mass ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ, that scales roughly linearly with the electron temperature jump. Our simulations haveΔϕ/(0.5miush2)0.1–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕin our low-βpshocks. We further focus on twoθBn= 65° shocks: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis the ion skin depth. Within the precursors,ϕhas a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

     
    more » « less
  3. Abstract

    Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low (x¯HI103) or close to unity (x¯HI1). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints onx¯HIatz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits ofx¯HI(z=6.3)<0.79±0.04(1σ),x¯HI(z=6.5)<0.87±0.03(1σ), andx¯HI(z=6.7)<0.940.09+0.06(1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.

     
    more » « less
  4. Abstract

    We measure the CO-to-H2conversion factor (αCO) in 37 galaxies at 2 kpc resolution, using the dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metal ratio. In total, we have ∼790 and ∼610 independent measurements ofαCOfor CO (2–1) and (1–0), respectively. The mean values forαCO (2–1)andαCO (1–0)are9.35.4+4.6and4.22.0+1.9Mpc2(Kkms1)1, respectively. The CO-intensity-weighted mean is 5.69 forαCO (2–1)and 3.33 forαCO (1–0). We examine howαCOscales with several physical quantities, e.g., the star formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength (U¯). Among them,U¯, ΣSFR, and the integrated CO intensity (WCO) have the strongest anticorrelation with spatially resolvedαCO. We provide linear regression results toαCOfor all quantities tested. At galaxy-integrated scales, we observe significant correlations betweenαCOandWCO, metallicity,U¯, and ΣSFR. We also find thatαCOin each galaxy decreases with the stellar mass surface density (Σ) in high-surface-density regions (Σ≥ 100Mpc−2), following the power-law relationsαCO(21)Σ0.5andαCO(10)Σ0.2. The power-law index is insensitive to the assumed dust-to-metal ratio. We interpret the decrease inαCOwith increasing Σas a result of higher velocity dispersion compared to isolated, self-gravitating clouds due to the additional gravitational force from stellar sources, which leads to the reduction inαCO. The decrease inαCOat high Σis important for accurately assessing molecular gas content and star formation efficiency in the centers of galaxies, which bridge “Milky Way–like” to “starburst-like” conversion factors.

     
    more » « less
  5. Abstract

    We present13CO(J= 1 → 0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected12CO emission as a prior, we detect13CO in 41 galaxies via integrated line flux over the entire galaxy and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio12/13I[12CO(J=10)]/I[13CO(J=10)]and the properties of the stars and ionized gas. Higher12/13values are found in interacting galaxies compared to those in noninteracting galaxies. The global12/13slightly increases with infrared colorF60/F100but appears insensitive to other host-galaxy properties such as morphology, stellar mass, or galaxy size. We also present azimuthally averaged12/13profiles for our sample up to a galactocentric radius of 0.4r25(∼6 kpc), taking into account the13CO nondetections by spectral stacking. The radial profiles of12/13are quite flat across our sample. Within galactocentric distances of 0.2r25, the azimuthally averaged12/13increases with the star formation rate. However, Spearman rank correlation tests show the azimuthally averaged12/13does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks,12/13is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on12/13, which further complicates the interpretations of12/13variations.

     
    more » « less