Emergence of cyclic flux eruptions in kinetic simulations of magnetized spherical accretion onto a Schwarzschild black hole

The dynamics of black hole magnetospheres critically depend on the black hole spin and on the structure of the accretion flow. In the limit of a Schwarzschild black hole immersed in a zero-net angular momentum flow, accretion is spherical. However, in the presence of a large-scale vertical magnetic field, the classical Bondi accretion model is significantly altered. The frozen-in field is stretched radially as the plasma is pulled inward by gravity. This continues until the restoring force from the magnetic tension suddenly expels the material and resets the field, allowing a new cycle to begin. Although this scenario has been well depicted in previous studies, it remains incomplete as the issues of dissipation and particle acceleration are not yet fully resolved. In this work, we aim to revisit these issues with a first-principles kinetic plasma model. We perform two-dimensional global general relativistic particle-in-cell simulations of magnetized spherical accretion onto a Schwarzschild black hole, for both pair and electron-ion plasmas. The simulations are evolved over long timescales to capture multiple flux eruption events and establish a quasi-steady state. For each accretion cycle, we find that the system goes through three main stages: (i) an ideal advection phase where magnetic flux through the horizon increases quasi-linearly with time; (ii) a reconnection-regulated phase where the net increase of the flux is slowed down by intermittent reconnection events near the horizon; and (iii) a flaring phase when a major, large-scale reconnection event expels the flux, leading to efficient particle acceleration. The emergence of large-amplitude quasi-periodic flux eruptions and concomitant particle acceleration is reminiscent of Sgr A* flaring activity. This phenomenon could also be applicable to quiescent black holes, especially isolated black holes accreting the interstellar medium.

💡 Research Summary

This paper presents the first long‑duration, global general‑relativistic particle‑in‑cell (GR‑PIC) simulations of magnetized spherical accretion onto a non‑spinning (Schwarzschild) black hole, with the goal of elucidating the microphysical processes that drive quasi‑periodic magnetic flux eruptions and associated particle acceleration. The authors consider a zero‑angular‑momentum inflow (Bondi‑type accretion) embedded in an initially uniform vertical magnetic field (the Wald solution for a non‑rotating spacetime). Two plasma compositions are explored: an electron‑positron pair plasma (mass ratio m_i/m_e = 1) and electron‑ion plasmas with reduced mass ratios of 16 and 256, chosen to make the problem computationally tractable while still capturing the essential ion–electron scale separation. The simulations are carried out with the Zeltron GR‑PIC code in spherical Kerr‑Schild coordinates, using a logarithmically spaced radial grid and uniform polar grid, with typical resolutions of N_r = N_θ = 1024–2048 cells. The outer boundary is matched to the initial uniform field, while a thin spherical shell near the outer edge continuously injects fresh plasma to maintain a constant density floor (n ≥ n₀). The initial plasma is cold (k_B T₀/(m_i) = 1/30) and weakly magnetized (σ₀ ≈ 0.05–0.1, β ≳ 1), ensuring that the flow can initially fall freely under gravity.

Over simulation times of several thousand gravitational times (t_g = r_g/c), the system settles into a quasi‑steady, cyclic behavior that repeats many times. Each cycle consists of three distinct phases:

1. Ideal advection phase – The inflowing plasma drags magnetic field lines toward the horizon, causing the magnetic flux threading the black hole (Φ_H) to increase almost linearly with time. The mass accretion rate (Ṁ) also rises, and the ratio Φ_H/√Ṁ grows until it reaches a saturation value (~60 in code units).

2. Reconnection‑regulated phase – Near the equatorial plane a thin current sheet forms as oppositely directed field lines are pushed together by the accumulating flux. Small‑scale, intermittent reconnection events begin to appear, partially dissipating magnetic energy and slowing the growth of Φ_H. Nevertheless, flux continues to accumulate because the gravitational pull still dominates.

3. Flaring (eruption) phase – When the current sheet becomes sufficiently thin, a large‑scale reconnection event is triggered. This event rapidly expels a substantial fraction of the stored flux (Φ_H drops by a factor of ~4 within ~100 t_g), produces a sharp spike in the electromagnetic dissipation rate (J·E), and accelerates particles to high Lorentz factors. The average particle energy ⟨γ⟩ jumps, and a non‑thermal power‑law tail develops in the particle energy distribution.

The authors quantify the energetic efficiency of the process. Averaged over many cycles, the electromagnetic power converted into particle energy is ⟨Ė⟩_cycle ≈ 0.0044 Ṁ₀, giving a cycle‑averaged efficiency η_cycle ≈ 0.03. During the brief flaring interval, the instantaneous efficiency rises to η_flare ≈ 0.2, indicating that roughly 20 % of the accreted rest‑mass energy is instantaneously transferred to particles. These numbers are robust across the different mass‑ratio runs, suggesting that the qualitative picture does not depend sensitively on the exact ion‑to‑electron mass ratio.

Spatially, the simulation reveals a clear dichotomy. Magnetic field lines that connect to the horizon form a highly magnetized funnel (σ ≫ 1) through which plasma streams inward along the field. The rest of the domain is occupied by a dense, weakly magnetized equatorial plasma that drags the vertical field toward the black hole, thereby building up Φ_H. The equatorial current sheet, where reconnection occurs, sits at the interface between these two regions. The large‑scale reconnection that triggers a flare effectively reverses the polarity of the radial field in the equatorial zone, expelling magnetic energy and temporarily halting accretion.

The paper discusses the astrophysical relevance of these results. The quasi‑periodic flux eruptions and associated particle acceleration bear a striking resemblance to the observed near‑infrared and X‑ray flares of Sgr A*, both in timescale (∼10³ t_g) and in the inferred efficiency of converting accretion power into high‑energy particles. Because the simulations deliberately exclude black‑hole spin, the results demonstrate that strong flares can arise purely from magnetically arrested spherical inflow, without invoking the Blandford‑Znajek mechanism. The authors also argue that isolated stellar‑mass black holes accreting from the interstellar medium could exhibit similar episodic flares, potentially detectable as transient radio or X‑ray events.

Limitations are acknowledged. The study is two‑dimensional and axisymmetric, which precludes three‑dimensional instabilities such as kink modes that could modify the reconnection dynamics. The reduced ion‑to‑electron mass ratios, while necessary for computational feasibility, may affect the detailed shape of the particle spectra. Radiative cooling, pair production, and the influence of black‑hole spin are omitted, and future work should incorporate these effects to build a more complete model.

In summary, this work provides the first kinetic‑level demonstration that magnetized spherical accretion onto a Schwarzschild black hole naturally evolves into a cycle of magnetic flux buildup, reconnection‑regulated growth, and violent flux eruption. The large‑scale reconnection events efficiently accelerate particles, offering a plausible microphysical explanation for the flaring activity observed in low‑luminosity galactic nuclei and possibly for transient emission from isolated black holes. The study bridges the gap between fluid GRMHD simulations of magnetically arrested disks and the underlying plasma physics, highlighting the crucial role of collisionless reconnection in black‑hole accretion environments.