More Like this

1. Mapping the Growth of Supermassive Black Holes as a Function of Galaxy Stellar Mass and Redshift

   https://doi.org/10.3847/1538-4357/ad27cc  

   Zou, Fan; Yu, Zhibo; Brandt, W N; Tak, Hyungsuk; Yang, Guang; Ni, Qingling (March 2024, The Astrophysical Journal)

   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 ( $\overline{\mathrm{BHAR}}$) primarily correlates with the galaxy stellar mass (M_⋆) and redshift for the general galaxy population. This work aims to provide the best measurements of $\overline{\mathrm{BHAR}}$as a function ofM_⋆and redshift over ranges of 10^9.5<M_⋆< 10¹²M_⊙andz< 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 estimate $\overline{\mathrm{BHAR}}$and the corresponding uncertainties, even for sparsely populated regions in the parameter space. $\overline{\mathrm{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 measure $\overline{\mathrm{BHAR}}$for star-forming and quiescent galaxies and show that star-forming $\overline{\mathrm{BHAR}}$is generally larger than or at least comparable to the quiescent $\overline{\mathrm{BHAR}}$. 

   more » « less
2. Length and Velocity Scales in Protoplanetary Disk Turbulence

   https://doi.org/10.3847/1538-4357/ad2c89  

   Sengupta, Debanjan; Cuzzi, Jeffrey N; Umurhan, Orkan M; Lyra, Wladimir (April 2024, The Astrophysical Journal)

   Abstract In the theory of protoplanetary disk turbulence, a widely adopted ansatz, or assumption, is that the turnover frequency of the largest turbulent eddy, Ω_L, is the local Keplerian frequency Ω_K. In terms of the standard dimensionless Shakura–Sunyaevαparameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by $\sqrt{\alpha }H$and $\sqrt{\alpha }c$, in whichHandcare the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as vertical shear instability. Here, we explore the more general case where Ω_L≥ Ω_Kand show that, under these conditions, the characteristic length and velocity scales are respectively $\sqrt{\alpha /R\prime }H$and $\sqrt{\alpha R\prime }c$, where $R\prime \equiv {\mathrm{\Omega }}_L/{\mathrm{\Omega }}_K$is twice the Rossby number. It follows that $\alpha =\tilde{\alpha }/R\prime$, where $\sqrt{\tilde{\alpha }}c$is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and it may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For $R\prime >1$, the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation. 

   more » « less
3. Resolved Measurements of the CO-to-H ₂ Conversion Factor in 37 Nearby Galaxies

   https://doi.org/10.3847/1538-4357/ad23ed  

   Chiang_江, I-Da 宜達; Sandstrom, Karin M; Chastenet, Jérémy; Bolatto, Alberto D; Koch, Eric W; Leroy, Adam K; Sun_孙, Jiayi 嘉懿; Teng, Yu-Hsuan; Williams, Thomas G (March 2024, The Astrophysical Journal)

   Abstract We measure the CO-to-H₂conversion 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)}are ${9.3}_{-5.4}^{+4.6}$and ${4.2}_{-2.0}^{+1.9}{M}_{\odot }{\mathrm{pc}}^{-2}{\left(\mathrm{K}\mathrm{km}{\mathrm{s}}^{-1}\right)}^{-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 ( $\overline{U}$). Among them, $\overline{U}$, Σ_SFR, and the integrated CO intensity (W_CO) 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α_COandW_CO, metallicity, $\overline{U}$, and Σ_SFR. We also find thatα_COin each galaxy decreases with the stellar mass surface density (Σ_⋆) in high-surface-density regions (Σ_⋆≥ 100M_⊙pc⁻²), following the power-law relations ${\alpha }_{\mathrm{CO}(2–1)}\propto {\mathrm{\Sigma }}_{\star }^{-0.5}$and ${\alpha }_{\mathrm{CO}(1–0)}\propto {\mathrm{\Sigma }}_{\star }^{-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
4. Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

   https://doi.org/10.3847/2041-8213/ad1048  

   Cho 조, Hyerin 혜린; Prather, Ben S.; Narayan, Ramesh; Natarajan, Priyamvada; Su, Kung-Yi; Ricarte, Angelo; Chatterjee, Koushik (December 2023, The Astrophysical Journal Letters)

   Abstract Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius ofR_B≈ 2 × 10⁵GM_•/c², whereM_•is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼r⁻¹rather thanr^−3/2and the accretion rate $\dot{M}$is consequently suppressed by over 2 orders of magnitude below the Bondi rate ${\dot{M}}_{\mathrm{B}}$. We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim {10}^{-2}\dot{M}{c}^2\sim 5\times {10}^{-5}{\dot{M}}_{\mathrm{B}}{c}^2$. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback. 

   more » « less
5. Safety First: Stability and Dissipation of Line-tied Force-free Flux Tubes in Magnetized Coronae

   https://doi.org/10.3847/1538-4357/ad3206  

   Rugg, N; Mahlmann, J F; Spitkovsky, A (May 2024, The Astrophysical Journal)

   Abstract Magnetized plasma columns and extended magnetic structures with both footpoints anchored to a surface layer are an important building block of astrophysical dissipation models. Current loops shining in X-rays during the growth of plasma instabilities are observed in the corona of the Sun and are expected to exist in highly magnetized neutron star magnetospheres and accretion disk coronae. For varying twist and system sizes, we investigate the stability of line-tied force-free flux tubes and the dissipation of twist energy during instabilities using linear analysis and time-dependent force-free electrodynamics simulations. Kink modes (m= 1) and efficient magnetic energy dissipation develop for plasma safety factorsq≲ 1, whereqis the inverse of the number of magnetic field line windings per column length. Higher-order fluting modes (m> 1) can distort equilibrium flux tubes forq> 1 but induce significantly less dissipation. In our analysis, the characteristic pitch ${\tilde{\mu }}_0$of flux-tube field lines determines the growth rate ( $\propto {\tilde{\mu }}_0^3$) and minimum wavelength of the kink instability ( $\propto {\tilde{\mu }}_0^{-1}$). We use these scalings to determine a minimum flux tube length for the growth of the kink instability for any given ${\tilde{\mu }}_0$. By drawing analogies to idealized magnetar magnetospheres with varying regimes of boundary shearing rates, we discuss the expected impact of the pitch-dependent growth rates for magnetospheric dissipation in magnetar conditions. 

   more » « less