Krzysztof Nalewajko

Magnetized accretion onto spinning black holes (BHs) can accumulate a large magnetic flux across the event horizon and launch a pair of relativistic jets via the Blandford–Znajek mechanism. In the magnetically saturated (arrested) state, excess magnetic flux is ejected from the BH in episodic magnetic flux eruptions, which result in a significant yet temporary reduction of jet power. We analyze results of a high-resolution 3D general-relativistic magnetohydrodynamic numerical simulation of geometrically thick magnetically saturated accretion onto a high-spin Kerr BH for a single cycle of magnetic flux eruption and accumulation. We show that following an eruption, a weakened jet develops a strong helical distortion with a distinct structure of magnetic fields—the poloidal field along the jet core is unaffected by the eruption; while toroidal field lines, ejected from the BH during the eruption and later readvected onto it, form poloidal “bypasses” along the inner jet sheath. Such a distortion may appear in sources fed by geometrically thick accretion flows as an asymmetric superluminal knot, strongly interacting with the jet sheath along an oblique working surface. The jet section repowered by magnetic flux reaccumulated on the BH is tilted by a few degrees, implying significant variations in radiation boost towards observers of BL Lac blazars. The intrinsic structure of the jet spine is consistent with axisymmetric semi-analytical models.

The toroidal magnetic field is a key ingredient of relativistic jets launched by certain accreting astrophysical black holes, and of plasmoids emerging from the tearing instability during magnetic reconnection, which is a candidate dissipation mechanism in jets. Tension of the toroidal field is an anisotropic force that can compress local energy and momentum densities. We investigate this effect in plasmoids produced during relativistic reconnection initiated from a Harris layer by means of kinetic particle-in-cell numerical simulations, varying the system size (including 3D cases), magnetisation, or guide field. We find that: (1) plasmoid cores are dominated by plasma energy density for guide fields up to B_z ∼ B₀; (2) relaxed ‘monster’ plasmoids compress plasma energy density only modestly (by a factor of ∼3 above the initial level for the drifting particle population); (3) energy density compressions by factors ≳10 are achieved during plasmoid mergers, especially with the emergence of secondary plasmoids. This kinetic-scale effect can be combined with a global focusing of the jet Poynting flux along the quasi-cylindrical bunched spine (a proposed jet layer adjacent to the cylindrical core) due to poloidal line bunching (a prolonged effect of tension in the jet toroidal field) to enhance the luminosity of rapid radiation flares from blazars. The case of M87 as a misaligned blazar is discussed.

Toroidal magnetic field is a key ingredient of relativistic jets launched by certain accreting astrophysical black holes, and of plasmoids emerging from the tearing instability during magnetic reconnection, a candidate dissipation mechanism in jets. Tension of toroidal field is an anisotropic force that can compress local energy and momentum densities. We investigate this effect in plasmoids produced during relativistic reconnection initiated from a Harris layer by means of kinetic particle-in-cell (PIC) numerical simulations, varying the system size (including 3D cases), magnetisation, or guide field. We find that: (1) plasmoid cores are dominated by plasma energy density for guide fields up to B_z ~ B_0; (2) relaxed 'monster' plasmoids compress plasma energy density only modestly (by factor ~3 above the initial level for drifting particle population); (3) energy density compressions by factors >~10 are achieved during plasmoid mergers, especially with the emergence of secondary plasmoids. This kinetic-scale effect can be combined with a global focusing of the jet Poynting flux along the quasi-cylindrical bunched spine (a proposed jet layer adjacent to the cylindrical core) due to poloidal line bunching (a prolonged effect of tension of the jet toroidal field) to enhance the luminosity of rapid radiation flares from blazars. The case of M87 as a misaligned blazar is discussed.

Context. Magnetized accretion flow onto a black hole (BH) may lead to the accumulation of poloidal magnetic flux across its horizon, which for high BH spin can power far-reaching relativistic jets. The BH magnetic flux is subject to a saturation mechanism by means of magnetic flux eruptions involving relativistic magnetic reconnection. Such accretion flows have been described as magnetically arrested disks (MAD) or magnetically choked accretion flows (MCAF). Aims. The main goal of this work is to describe the onset of relativistic reconnection and initial development of magnetic flux eruption in accretion flow onto magnetically saturated BHs. Methods. We analyzed the results of 3D general relativistic ideal magnetohydrodynamic (GRMHD) numerical simulations in the Kerr metric, starting from weakly magnetized geometrically thick tori rotating either prograde or retrograde. We integrated large samples of magnetic field lines in order to probe magnetic connectivity with the BH horizon. Results. The boundary between magnetically connected and disconnected domains coincides roughly with enthalpy equipartition. The geometrically constricted innermost part of the disconnected domain develops a rigid structure of magnetic field lines – rotating slowly and insensitive to the BH spin orientation. The typical shape of innermost disconnected lines is a double spiral converging to a sharp inner tip anchored at the single equatorial current layer. The foot-points of magnetic flux eruptions are found to zip around the BH along with other azimuthal patterns. Conclusions. Magnetic flux eruptions from magnetically saturated accreting BHs can be triggered by minor density gaps in the disconnected domain, resulting from the chaotic disconnection of plasma-depleted magnetospheric lines. Accretion flow is effectively channeled along the disconnected lines toward the current layer, and further toward the BH by turbulent cross-field diffusion. Rotation of flux eruption foot-points may contribute to the variability of BH crescent images.

Magnetized accretion flow onto a black hole (BH) may lead to accumulation of poloidal magnetic flux across its horizon, which for high BH spin can power far-reaching relativistic jets. The BH magnetic flux is subject to a saturation mechanism by means of magnetic flux eruptions involving relativistic magnetic reconnection. Such accretion flows have been described as magnetically arrested (MAD) or choked (MCAF). The main goal of this work is to describe the onset of relativistic reconnection and initial development of magnetic flux eruption in accretion flow onto magnetically saturated BH. We analyze the results of 3D ideal GRMHD numerical simulations in the Kerr metric, starting from weakly magnetized geometrically thick tori rotating either prograde or retrograde. We integrate large samples of magnetic field lines in order to probe magnetic connectivity with the BH horizon. The boundary between magnetically connected and disconnected domains coincides roughly with enthalpy equipartition. The geometrically constricted innermost part of the disconnected domain develops a rigid structure of magnetic field lines - rotating slowly and insensitive to the BH spin orientation. The typical shape of innermost disconnected lines is a double spiral converging to a sharp inner tip anchored at the single equatorial current layer. The footpoints of magnetic flux eruptions are found to zip around the BH along with other azimuthal patterns. Magnetic flux eruptions from magnetically saturated accreting BHs can be triggered by minor density gaps in the disconnected domain, resulting from chaotic disconnection of plasma-depleted magnetospheric lines. Accretion flow is effectively channeled along the disconnected lines towards the current layer, and further towards the BH by turbulent cross-field diffusion. Rotation of flux eruption footpoints may contribute to the variability of BH crescent images.

Black holes (BH) acquire relativistic magnetospheres by accreting magnetized gas. Once they collect significant magnetic flux across the horizon, aided by the spin they can drive powerful relativistic jets by the Blandford-Znajek mechanism. Large enough BH magnetic flux backreacts on the accretion flow, which has been described in terms of arresting or choking. Magnetic flux eruptions have been identified as the mechanism of BH magnetic flux saturation. These eruptions can potentially dissipate a large fraction of magnetic energy in the BH magnetosphere by means of relativistic magnetic reconnection, accelerating particles and producing flares of non-thermal radiation. We analyze the results of 3D general-relativistic ideal magnetohydrodynamic (GRMHD) numerical simulations of accretion flows onto magnetically saturated Kerr BHs, focusing on the initiation of magnetic flux eruptions.

We present and analyze a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a nonrotating, stellar-mass black hole. The simulations are initialized using the Shakura–Sunyaev model with a mass accretion rate of M ̇ = 3 L Edd / c 2 (corresponding to L = 0.17L _Edd). Our previous work demonstrated that such disks are thermally unstable when accretion is driven by an α-viscosity. In the present work, we test the hypothesis that strong magnetic fields can both drive accretion through magnetorotational instability and restore stability to such disks. We test four initial magnetic field configurations: (1) a zero-net-flux case with a single, radially extended set of magnetic field loops (dipole), (2) a zero-net-flux case with two radially extended sets of magnetic field loops of opposite polarity stacked vertically (quadrupole), (3) a zero-net-flux case with multiple radially concentric rings of alternating polarity (multiloop), and (4) a net-flux, vertical magnetic field configuration (vertical). In all cases, the fields are initially weak, with a gas-to-magnetic pressure ratio ≳100. Based on the results of these simulations, we find that the dipole and multiloop configurations remain thermally unstable like their α-viscosity counterpart, in our case collapsing vertically on the local thermal timescale and never fully recovering. The vertical case, on the other hand, stabilizes and remains so for the duration of our tests (many thermal timescales). The quadrupole case is intermediate, showing signs of both stability and instability. The key stabilizing factor is the ability of specific field configurations to build up and sustain strong, P _mag ≳ 0.5P _tot, toroidal fields near the midplane of the disk. We discuss the reasons why certain configurations are able to do this effectively and others are not. We then compare our stable simulations to the standard Shakura–Sunyaev disk.

We propose a simple analytical jet model of magnetic jets, in which radially averaged profiles of main physical quantities are obtained based on conservation laws and some results of published general relativistic magnetohydrodynamic jet simulations. We take into account conversion of the magnetic energy flux to bulk acceleration in jets formed around rotating black holes assuming the mass continuity equation and constant jet power, which leads to the Bernoulli equation. For assumed profiles of the bulk Lorentz factor and the radius, this gives us the profile of the toroidal magnetic field component along the jet. We then consider the case where the poloidal field component is connected to a rotating black hole surrounded by an accretion disc. Our formalism then recovers the standard formula for the power extracted from a rotating black hole. We find that the poloidal field strength dominates over the toroidal one in the comoving frame up to large distances, which means that jets should be more stable to current-driven kink modes. The resulting magnetic field profiles can then be used to calculate the jet synchrotron emission.

Relativistic magnetized jets, such as those from AGN, GRBs, and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and nonthermal radiation. Here, we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron–positron pair plasma. Generalizing recent treatments by Alves et al. and Davelaar et al., we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases, a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.