More Like this

1. 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)

   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
2. Electron Heating in 2D Particle-in-cell Simulations of Quasi-perpendicular Low-beta Shocks

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

   Tran, Aaron ; Sironi, Lorenzo ( April 2024 , The Astrophysical Journal)

   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 ratiom_i/m_e= 200, fast Mach number${\mathcal{M}}_{\mathrm{ms}}=1$–4, and upstream magnetic field angleθ_Bn= 55°–85° from the shock normal$\hat{\mathit{n}}$. 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$\mathrm{\Delta }{\varphi }_{\parallel }/(0.5m_{\mathrm{i}}{u_{\mathrm{sh}}}^2)\sim 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: a${\mathcal{M}}_{\mathrm{s}}=4$(${\mathcal{M}}_{\mathrm{A}}=1.8$) case with a long, 30d_iprecursor of whistler waves along$\hat{\mathit{n}}$, and a${\mathcal{M}}_{\mathrm{s}}=7$(${\mathcal{M}}_{\mathrm{A}}=3.2$) case with a shorter, 5d_iprecursor of whistlers oblique to both$\hat{\mathit{n}}$andB;d_iis 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 the${\mathcal{M}}_{\mathrm{s}}=4$,θ_Bn= 65° case,ϕ_∥shows a weak dependence on the electron plasma-to-cyclotron frequency ratioω_pe/Ω_ce, andϕ_∥decreases by a factor of 2 asm_i/m_eis raised to the true proton–electron value of 1836.

    

   more » « less
3. Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

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

   Grošelj, Daniel ; Sironi, Lorenzo ; Beloborodov, Andrei M. ( July 2022 , The Astrophysical Journal)

   We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizationsσ≲ 10⁻⁴and pair-loading factorsZ_±≲ 10 are studied, whereZ_±is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field whenσexceeds a critical valueσ_Lthat decreases withZ_±. Atσ≲σ_Lthe shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease withZ_±, leading to lowerσ_L. (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as${\overline{E}}_{\mathrm{e}}\propto {\left(Z_{\pm }+1\right)}^{-1}$. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ≈ 5 × 10⁻⁶. The ions then become essentially thermal with mean energy${\overline{E}}_{\mathrm{i}}$, while electrons form a nonthermal tail, extending from$E\sim {\left(Z_{\pm }+1\right)}^{-1}{\overline{E}}_{\mathrm{i}}$to${\overline{E}}_{\mathrm{i}}$. Whenσ= 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows.

    

   more » « less
4. Resolved Measurements of the CO-to-H 2 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)

   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}{(\mathrm{K}\mathrm{km}{\mathrm{s}}^{-1})}^{-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
5. Cosmic-Ray Drag and Damping of Compressive Turbulence

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

   Bustard, Chad ; Oh, S. Peng ( September 2023 , The Astrophysical Journal)

   While it is well known that cosmic rays (CRs) can gain energy from turbulence via second-order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping time$t_{\mathrm{damp}}\sim \rho v^2/{\dot{E}}_{\mathrm{CR}}\propto E_{\mathrm{CR}}^{-1}$becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating at${\dot{E}}_{\mathrm{CR}}\sim \widetilde{ϵ}$, where$\widetilde{ϵ}$is the turbulent energy input rate. In that case, almost all of the energy in large-scale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge–Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering ofE≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large-scale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticE_CR∼E_gscenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large-scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e.,$\widetilde{ϵ}\gg \rho v^3/L$.

    

   more » « less