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)

   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
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. Magnetogram-matching Biot–Savart Law and Decomposition of Vector Magnetograms

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

   Titov, Viacheslav S; Downs, Cooper; Török, Tibor; Linker, Jon A; Prazak, Michael; Qiu, Jiong A (July 2025, The Astrophysical Journal)

   Abstract We generalize a magnetogram-matching Biot–Savart law (BSl) from planar to spherical geometry. For a given coronal current densityJ, this law determines the magnetic field $\tilde{\mathit{B}}$whose radial component vanishes at the surface. The superposition of $\tilde{\mathit{B}}$with a potential field defined by a given surface radial field,B_r, provides the entire configuration whereB_rremains unchanged by the currents. Using this approach, we (1) upgrade our regularized BSls for constructing coronal magnetic flux ropes (MFRs) and (2) propose a new method for decomposing a measured photospheric magnetic field as $\mathit{B}={\mathit{B}}_{\mathrm{pot}}+{\mathit{B}}_T+{\mathit{B}}_{\tilde{S}}$, where the potential,B_pot, toroidal,B_T, and poloidal, ${\mathit{B}}_{\tilde{S}}$, fields are determined byB_r,J_r, and the surface divergence ofB–B_pot, respectively, all derived from magnetic data. OurB_Tis identical to the one in the alternative Gaussian decomposition by P. W. Schuck et al., whileB_potand ${\mathit{B}}_{\tilde{S}}$are different from their poloidal fields ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$, which arepotentialin the infinitesimal proximity to the upper and lower side of the surface, respectively. In contrast, our ${\mathit{B}}_{\tilde{S}}$has no such constraints and, asB_potandB_T, refers to thesameupper side of the surface. In spite of these differences, for a continuousJdistribution across the surface,B_potand ${\mathit{B}}_{\tilde{S}}$are linear combinations of ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$. We demonstrate that, similar to the Gaussian method, our decomposition allows one to identify the footprints and projected surface-location of MFRs in the solar corona, as well as the direction and connectivity of their currents. 

   more » « less
4. (Nearly) Model-independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at z ∼ 5–7 Using Dark Pixel Fractions in Lyα and Lyβ Forests

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

   Jin, Xiangyu; Yang, Jinyi; Fan, Xiaohui; Wang, Feige; Bañados, Eduardo; Bian, Fuyan; Davies, Frederick B.; Eilers, Anna-Christina; Farina, Emanuele Paolo; Hennawi, Joseph F.; et al (January 2023, The Astrophysical Journal)

   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 ( ${\overline{x}}_{\mathrm{H}\mathrm{I}}\lesssim {10}^{-3}$) or close to unity ( ${\overline{x}}_{\mathrm{H}\mathrm{I}}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints on ${\overline{x}}_{\mathrm{H}\mathrm{I}}$atz∼ 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 of ${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.3)<0.79\pm 0.04$(1σ), ${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.5)<0.87\pm 0.03$(1σ), and ${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.7)<{0.94}_{-0.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
5. Efficient modeling of electron kinetics under influence of externally applied electric field in magnetized weakly ionized plasma

   https://doi.org/10.1088/1361-6595/acdaf1  

   Janalizadeh, Reza; Pervez, Zaid; Pasko, Victor P (July 2023, Plasma Sources Science and Technology)

   Abstract We present a theory based on the conventional two-term (i.e. Lorentzian) approximation to the exact solution of the Boltzmann equation in non-magnetized weakly ionized plasma to efficiently obtain the electron rate and transport coefficients in a magnetized plasma for an arbitrary magnitude and direction of applied electric field $\vec{E}$and magnetic field $\vec{B}$. The proposed transcendental method does not require the two-term solution of the Boltzmann equation in magnetized plasma, based on which the transport parameters vary as a function of the reduced electric field $E/N$, reduced electron cyclotron frequency ${\omega }_{\mathrm{c}\mathrm{e}}/N$, and angle $\mathrm{\angle }\vec{E},\vec{B}$between $\vec{E}$and $\vec{B}$vectors, whereNis the density of neutrals. Comparisons between the coefficients derived from BOLSIG+’s solution (obtained via the two-term expansion when $\vec{B}\ne 0$) and coefficients of the presented method are illustrated for air, a mixture of molecular hydrogen (H₂) and helium (He) representing the giant gas planets of the Solar System, and pure carbon dioxide (CO₂). The new approach may be used in the modeling of magnetized plasma encountered in the context of transient luminous events, e.g. sprite streamers in the atmosphere of Earth and Jupiter, in modeling the propagation of lightning’s electromagnetic pulses in Earth’s ionosphere, and in various laboratory and industrial applications of nonthermal plasmas. 

   more » « less