Search for: All records

Abstract In certain scenarios, the accreted angular momentum of plasma onto a black hole could be low; however, how the accretion dynamics depend on the angular momentum content of the plasma is still not fully understood. We present three-dimensional, general relativistic magnetohydrodynamic simulations of low angular momentum accretion flows around rapidly spinning black holes (with spina = +0.9). The initial condition is a Fishbone–Moncrief (FM) torus threaded by a large amount of poloidal magnetic flux, where the angular velocity is a fractionfof the standard value. Forf= 0, the accretion flow becomes magnetically arrested and launches relativistic jets but only for a very short duration. After that, free-falling plasma breaks through the magnetic barrier, loading the jet with mass and destroying the jet–disk structure. Meanwhile, magnetic flux is lost via giant, asymmetrical magnetic bubbles that float away from the black hole. The accretion then exits the magnetically arrested state. Forf= 0.1, the dimensionless magnetic flux threading the black hole oscillates quasiperiodically. The jet–disk structure shows concurrent revival and destruction while the gas outflow efficiency at the event horizon changes accordingly. Forf≥ 0.3, we find that the dynamical behavior of the system starts to approach that of a standard accreting FM torus. Our results thus suggest that the accreted angular momentum is an important parameter that governs the maintenance of a magnetically arrested flow and launching of relativistic jets around black holes. 

Abstract The presence of a strong, large-scale magnetic field in an accretion flow leads to the extraction of the rotational energy of the black hole (BH) through the Blandford–Znajek (BZ) process, believed to power relativistic jets in various astrophysical sources. We study rotational energy extraction from a BH surrounded by a highly magnetized thin disk by performing a set of 3D global GRMHD simulations. We find that the saturated flux threading the BH has a weaker dependence on BH spin, compared to highly magnetized hot (geometrically thick) accretion flows. Also, we find that only a fraction (10%–70%) of the extracted BZ power is channeled into the jet, depending on the spin parameter. The remaining energy is potentially used to launch winds or contribute to the radiative output of the disk or corona. Our simulations reveal that the presence of a strong magnetic field enhances the radiative efficiency of the disk, making it more luminous than its weakly magnetized counterpart or the standard disk model. We attribute this excess luminosity primarily to the enhanced magnetic dissipation in the intra-ISCO region. Our findings have implications for understanding X-ray corona formation and BH spin measurements, and interpreting BH transient phenomena. 

We present the largest 3D particle-in-cell shearing-box simulations of turbulence driven by the magnetorotational instability, for the first time employing the realistic proton-to-electron mass ratio. We investigate the energy partition between relativistically hot electrons and subrelativistic ions in turbulent accreting plasma, a regime relevant to collisionless, radiatively inefficient accretion flows around supermassive black holes such as those targeted by the Event Horizon Telescope. We provide a simple empirical formula to describe the measured heating ratio between ions and electrons, which can be used for more accurate global modeling of accretion flows with standard fluid approaches such as general-relativistic magnetohydrodynamics. 

Abstract We present results from a high-cadence multiwavelength observational campaign of the enigmatic changing-look active galactic nucleus 1ES 1927+654 from 2022 May to 2024 April, coincident with an unprecedented radio flare (an increase in flux by a factor of ∼60 over a few months) and the emergence of a spatially resolved jet at 0.1–0.3 pc scales. Companion work has also detected a recurrent quasi-periodic oscillation (QPO) in the 2–10 keV band with an increasing frequency (1–2 mHz) over the same period. During this time, the soft X-rays (0.3–2 keV) monotonically increased by a factor of ∼8, while the UV emission remained nearly steady with <30% variation and the 2–10 keV flux showed variation by a factor ≲2. The weak variation of the 2–10 keV X-ray emission and the stability of the UV emission suggest that the magnetic energy density and accretion rate are relatively unchanged and that the jet could be launched owing to a reconfiguration of the magnetic field (toroidal to poloidal) close to the black hole. Advecting poloidal flux onto the event horizon would trigger the Blandford–Znajek mechanism, leading to the onset of the jet. The concurrent softening of the coronal slope (from Γ = 2.70 ± 0.04 to Γ = 3.27 ± 0.04), the appearance of a QPO, and the low coronal temperature ( $k{T}_e=8_{-3}^{+8}\mathrm{keV}$) during the radio outburst suggest that the poloidal field reconfiguration can significantly impact coronal properties and thus influence jet dynamics. These extraordinary findings in real time are crucial for coronal and jet plasma studies, particularly as our results are independent of coronal geometry. 

Abstract We present multifrequency (5–345 GHz) and multiresolution radio observations of 1ES 1927+654, widely considered one of the most unusual and extreme changing-look active galactic nuclei (CL-AGNs). The source was first designated a CL-AGN after an optical outburst in late 2017 and has since displayed considerable changes in X-ray emission, including the destruction and rebuilding of the X-ray corona in 2019–2020. Radio observations prior to 2023 show a faint and compact radio source typical of a radio-quiet AGN. Starting in 2023 February, 1ES 1927+654 began exhibiting a radio flare with a steep exponential rise, reaching a peak 60 times previous flux levels, and has maintained this higher level of radio emission for over a year to date. The 5–23 GHz spectrum is broadly similar to gigahertz-peaked radio sources, which are understood to be young radio jets less than ∼1000 yr old. Recent high-resolution Very Long Baseline Array observations at 23.5 GHz now show resolved extensions on either side of the core, with a separation of ∼0.15 pc, consistent with a new and mildly relativistic bipolar outflow. A steady increase in the soft X-ray band (0.3–2 keV) concurrent with the radio may be consistent with jet-driven shocked gas, though further observations are needed to test alternate scenarios. This source joins a growing number of CL-AGNs and tidal disruption events that show late-time radio activity, years after the initial outburst. 

Abstract The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work. 

ABSTRACT Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of 10–100 for ≈60 min. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from particle-in-cell simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles, and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies γe ≲ 106 and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simplified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole, and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard (≲ 100 keV) X-ray energies. 

ABSTRACT Magnetically arrested accretion discs (MADs) around black holes (BHs) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here, we propose instead that it is mainly a dynamically important toroidal field – the result of dynamo action triggered by the significant but still relatively weak vertical field – that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotational instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (standard and normal evolution; strongly magnetized but not ‘arrested’) discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field. 

Abstract Astrophysical jets, launched from the immediate vicinity of accreting black holes, carry away large amounts of power in a form of bulk kinetic energy of jet particles and electromagnetic flux. Here we consider a simple analytical model for relativistic jets at larger distances from their launching sites, assuming a cylindrical axisymmetric geometry with a radial velocity shear, and purely toroidal magnetic field. We argue that as long as the jet plasma is in magnetohydrostatic equilibrium, such outflows tend to be particle dominated, i.e., the ratio of the electromagnetic to particle energy flux, integrated over the jet cross-sectional area, is typically below unity, σ < 1. At the same time, for particular magnetic and radial velocity profiles, magnetic pressure may still dominate over particle pressure for certain ranges of the jet radius, i.e., the local jet plasma parameter β pl < 1, and this may be relevant in the context of particle acceleration and production of high-energy emission in such systems. The jet magnetization parameter can be elevated up to the modest values of σ ≲  ( 10 ) only in the case of extreme gradients or discontinuities in the gaseous pressure, and a significantly suppressed velocity shear. Such configurations, which consist of a narrow, unmagnetized jet spine surrounded by an extended, force-free layer, may require an additional poloidal field component to stabilize them against current-driven oscillations, but even this will not substantially elevate their σ parameter. 

ABSTRACT Turbulent high-energy astrophysical systems often feature asymmetric energy injection: for instance, Alfvén waves propagating from an accretion disc into its corona. Such systems are ‘imbalanced’: the energy fluxes parallel and antiparallel to the large-scale magnetic field are unequal. In the past, numerical studies of imbalanced turbulence have focused on the magnetohydrodynamic regime. In this study, we investigate externally driven imbalanced turbulence in a collision-less, ultrarelativistically hot, magnetized pair plasma using 3D particle-in-cell (PIC) simulations. We find that the injected electromagnetic momentum efficiently converts into plasma momentum, resulting in net motion along the background magnetic field with speeds up to a significant fraction of lightspeed. This discovery has important implications for the launching of accretion disc winds. We also find that although particle acceleration in imbalanced turbulence operates on a slower time-scale than in balanced turbulence, it ultimately produces a power-law energy distribution similar to balanced turbulence. Our results have ramifications for black hole accretion disc coronae, winds, and jets.