More Like this

1. 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
2. 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
3. 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
4. Now You See It, Now You Don’t: Star Formation Truncation Precedes the Loss of Molecular Gas by ∼100 Myr in Massive Poststarburst Galaxies at z ∼ 0.6

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

   Bezanson, Rachel ; Spilker, Justin S. ; Suess, Katherine A. ; Setton, David J. ; Feldmann, Robert ; Greene, Jenny E. ; Kriek, Mariska ; Narayanan, Desika ; Verrico, Margaret ( February 2022 , The Astrophysical Journal)

   We use ALMA observations of CO(2–1) in 13 massive (M_*≳ 10¹¹M_⊙) poststarburst galaxies atz∼ 0.6 to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The poststarburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (Data Release 14) based on their spectral shapes, as part of the Studying QUenching at Intermediate-z Galaxies: Gas, angu$$\overrightarrow{L}\mathrm{ar}$$momentum, and Evolution ($\mathit{SQuIGG}\vec{L}E$) program. Early results showed that two poststarburst galaxies host large H₂reservoirs despite their low inferred star formation rates (SFRs). Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with$M_{{\mathrm{H}}_2}\gtrsim {10}^9$M_⊙. Given their high stellar masses, this mass limit corresponds to an average gas fraction of$〈f_{{\mathrm{H}}_2}\equiv M_{{\mathrm{H}}_2}/M_{*}〉\sim 7\%$or ∼14% using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with theD_n4000 spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched ≲150 Myr prior to observation host detectable CO(2–1) emission, while older poststarburst galaxies are undetected. The large H₂reservoirs and low SFRs in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100–200 Myr, the$\mathit{SQuIGG}\vec{L}E$galaxies require the additional and efficient heating or removal of cold gas to bring their low SFRs in line with standard H₂scaling relations.

    

   more » « less
5. 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)

   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