More Like this

1. Jets with a Twist: The Emergence of FR0 Jets in a 3D GRMHD Simulation of Zero-angular-momentum Black Hole Accretion

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

   Lalakos, Aretaios; Tchekhovskoy, Alexander; Bromberg, Omer; Gottlieb, Ore; Jacquemin-Ide, Jonatan; Liska, Matthew; Zhang, Haocheng (March 2024, The Astrophysical Journal)

   Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius,R_B= 10³R_g. The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideR_Band suppress BH accretion to 1.5% of the Bondi rate, ${\dot{M}}_{\mathrm{B}}$. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of ${\dot{M}}_{\mathrm{B}}$reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% ${\dot{M}}_{\mathrm{B}}$. Because the BAD and RAD states tangle up the jets and destroy them well insideR_B, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies. 

   more » « less
2. 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
3. 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
4. Acceleration and Spectral Redistribution of Cosmic Rays in Radio-jet Shear Flows

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

   Webb, G. M.; Xu, Y.; Biermann, P. L.; Al-Nussirat, S.; Mostafavi, P.; Li, G.; Barghouty, A. F.; Zank, G. P. (November 2023, The Astrophysical Journal)

   Abstract A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf₀(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrum $f_0(\infty ,p\prime )$with the Green’s function $G(r,p;p\prime )$, which describes the monoenergetic spectrum solution in which $f_0\to \delta (p-p\prime )$asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution for $G(r,p;p\prime )$. In this paper, we explore for the first time, solutions for more general and realistic forms for $f_0(\infty ,p\prime )$. The flow velocityu=u(r)e_zis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering time $\tau {(r,p)={\tau }_0\left(p/p_0\right)}^{\alpha }$in the shear flow region 0 <r<r₂, and $\tau {(r,p)={\tau }_0\left(p/p_0\right)}^{\alpha }{\left(r/r_2\right)}^s$, wheres> 0 in the regionr>r₂is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distribution ${\psi }_p(r,p;p\prime )$that particles observed at (r,p) originated fromr→ ∞ with momentum $p\prime$. The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described. 

   more » « less
5. Characterizing the fundamental bending vibration of a linear polyatomic molecule for symmetry violation searches

   https://doi.org/10.1088/1367-2630/ace471  

   Jadbabaie, Arian; Takahashi, Yuiki; Pilgram, Nickolas H; Conn, Chandler J; Zeng, Yi; Zhang, Chi; Hutzler, Nicholas R (July 2023, New Journal of Physics)

   Abstract Polyatomic molecules have been identified as sensitive probes of charge-parity violating and parity violating physics beyond the Standard Model (BSM). For example, many linear triatomic molecules are both laser-coolable and have parity doublets in the ground electronic $\tilde{X}{}^2{\mathrm{\Sigma }}^{+}\left(010\right)$state arising from the bending vibration, both features that can greatly aid BSM searches. Understanding the $\tilde{X}{}^2{\mathrm{\Sigma }}^{+}\left(010\right)$state is a crucial prerequisite to precision measurements with linear polyatomic molecules. Here, we characterize the fundamental bending vibration of ${}^{174}$YbOH using high-resolution optical spectroscopy on the nominally forbidden $\tilde{X}{}^2{\mathrm{\Sigma }}^{+}\left(010\right)$ $\to \tilde{A}{}^2{\mathrm{\Pi }}_{1/2}\left(000\right)$transition at 588 nm. We assign 39 transitions originating from the lowest rotational levels of the $\tilde{X}{}^2{\mathrm{\Sigma }}^{+}\left(010\right)$state, and accurately model the state’s structure with an effective Hamiltonian using best-fit parameters. Additionally, we perform Stark and Zeeman spectroscopy on the $\tilde{X}{}^2{\mathrm{\Sigma }}^{+}\left(010\right)$state and fit the molecule-frame dipole moment to $D_{\mathrm{m}\mathrm{o}\mathrm{l}}=2.16\left(1\right)$Dand the effective electrong-factor to $g_S=2.07\left(2\right)$. Further, we use an empirical model to explain observed anomalous line intensities in terms of interference from spin–orbit and vibronic perturbations in the excited $\tilde{A}{}^2{\mathrm{\Pi }}_{1/2}\left(000\right)$state. Our work is an essential step toward searches for BSM physics in YbOH and other linear polyatomic molecules. 

   more » « less