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

Adrien Soudais; Beno\^it Cerutti; Enzo Figueiredo; John Mehlhaff

REVIEW 3 major objections 4 minor 63 references

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

T0 review · 3 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that magnetized spherical accretion onto a Schwarzschild black hole is inherently cyclic, with a three-stage eruption cycle controlled by magnetic reconnection, and that the process reproduces the flaring behavior of Sgr A

desk verdict Kinetic GRPIC evidence for a three-phase flux-eruption cycle in Schwarzschild accretion, with a useful mass-ratio electron-acceleration trend; the analytic scalings are calibrated consistency checks rather than predictions. read the letter →

arxiv 2602.04519 v2 pith:UGVDV4LM submitted 2026-02-04 astro-ph.HE

Enzo Figueiredo , John Mehlhaff , Adrien Soudais , Benoît Cerutti This is my paper · ORCID

classification astro-ph.HE MSC 83C5785A3076W05

keywords blackholeaccretionsphericalmagneticreconnectionparticleaccelerationgeneralrelativisticparticle-in-cellSchwarzschildSgrA*flareskineticsimulations

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper uses global kinetic (particle-in-cell) simulations to show that magnetized spherical accretion onto a non-spinning black hole does not settle into a steady flow. Instead, it cycles: magnetic flux through the horizon grows almost linearly, then a reconnection-regulated phase slows the growth, and finally a large reconnection event expels most of the flux and accelerates particles to nonthermal energies. The cycle is quasi-periodic, lasting about a thousand gravitational timescales, with the eruption itself lasting about a tenth of that. The paper also derives analytic estimates for the saturated horizon flux (phi ~ 60), the exponential flux decay during eruption, and the current-sheet thinning that sets the phase transitions. A sympathetic reader would care because this provides a first-principles kinetic picture of collisionless accretion magnetospheres in the non-spinning case, and suggests that quiescent black holes such as Sgr A* may flare through this self-sustaining cycle.

What carries the argument

The argument rests on three analytic anchors: the flux-transport equation for the horizon magnetic flux, which explains the quasi-linear growth of Phi_H and yields a 1/r equatorial infall profile; a force-balance estimate between magnetic tension and gravity at reconnection X-points, which fixes the saturation ratio Phi_H/sqrt(Mdot) ~ 60; and tearing-instability thresholds of the equatorial current sheet, where the half-opening angle delta/(r-r_H) ~ 0.1 marks the onset of single-X-point reconnection and ~0.05 marks the transition to a plasmoid chain that triggers the eruption. A split-monopole reconnection-layer estimate gives the exponential decay of horizon flux during the eruption, with a

What would settle it

Run the same setup with a range of outer boundary radii (say, 0.5x, 1x, 1.7x the nominal r_max) and injection radii while holding resolution fixed: if the cycle period in units of t_g or the saturation phi ~ 60 changes beyond shot-to-shot scatter, the cycle is boundary-related. The paper already reports a 2/3-r_max reduction lowers the maximum Phi_H and stops sigma_H from saturating near 10, so a full scan would settle whether the nominal ~1000 t_g period is intrinsic. Observationally, measure the quasi-period of flares from an isolated stellar-mass black hole (e.g., via X-ray monitoring of a

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that zero-net-angular-momentum accretion onto a Schwarzschild black hole immersed in a vertical magnetic field is intrinsically cyclic. Each cycle has three phases: (i) ideal advection, during which horizon magnetic flux Phi_H and accretion rate grow linearly, (ii) a reconnection-regulated phase during which intermittent reconnection near the horizon slows the flux growth and the ratio Phi_H/sqrt(Mdot) saturates near 60, and (iii) an eruption phase in which a large-scale reconnection event removes roughly three quarters of the horizon flux within ~100 t_g and produces the dominant nonthermal particle acceleration. The phase transitions are a

Load-bearing premise

The simulated cycle depends on a finite simulation box with continuous plasma injection near the outer boundary and a density floor, and on the Bondi radius lying outside the box; if the recurrence time and the saturated horizon flux are set by that finite reservoir rather than by the intrinsic balance between advection and reconnection, the predicted cycle period and Sgr A* luminosity would not transfer to real systems.

Editorial extensions

If this is right

Magnetized spherical accretion onto a Schwarzschild black hole should be quasi-periodically flaring rather than steady, with a cycle period ~10^3 t_g and an eruption duration ~10^2 t_g.
For Sgr A*, the model predicts flare durations of roughly 30 minutes and luminosities near 10^35 erg/s for horizon fields of order 30 G, with most nonthermal emission confined to the eruption phase.
For isolated stellar-mass black holes accreting the interstellar medium, the same mechanism gives hard X-ray flares with durations of milliseconds and cycle periods of tens of milliseconds, potentially contributing to a diffuse Galactic high-energy background.
The electron energy spectrum is expected to harden and extend to higher energies as the proton-to-electron mass ratio increases, so real plasmas should produce electron Lorentz factors up to ~10^4.
Because the eruption removes a substantial fraction of the accumulated horizon flux, the system remains in a quasi-steady state over many cycles, with cycle-averaged efficiencies around 3%.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A direct test of whether the cycle is intrinsic rather than a numerical artifact: if the recurrence time (in units of t_g) or the saturation value phi ~ 60 shifts systematically with the simulation box size or injection radius, the cyclic behavior is set by the finite reservoir. The paper's box-size sensitivity run already hints at this, since shrinking r_max by 2/3 prevents sigma_H from saturatin
The paper's own three-phase picture predicts that the flare duty cycle (eruption duration divided by cycle period) is roughly constant across black hole masses; comparing the observed Sgr A* near-IR duty cycle with future observations of a stellar-mass isolated black hole would test this universality.
If full three-dimensional runs break axisymmetry, eruptions may be azimuthally localized rather than global; the paper's 2D constraint of simultaneous eruption at all azimuths is likely the main reason it cannot reproduce the orbital motion of Sgr A* hotspots.
The analytic force-balance derivation suggests a compact formula: phi_sat ~ 60 depends on the reconnection rate beta_rec and the X-point location; a modest change in beta_rec (e.g., from 0.05 to 0.1) would shift the predicted flare luminosity by a factor of 4, so the luminosity estimate is sensitive to the reconnection rate chosen.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Bottom line: this is a credible first GRPIC demonstration of the flux-eruption cycle in the Schwarzschild limit, run long enough to capture several quasi-steady cycles, and it makes a clean case that reconnection regulates both flux buildup and particle acceleration. The mass-ratio series (pair, 16, 256) is the genuinely new part: the electron tails extend and harden as m_i/m_e grows, and the interpretation via sigma_H(e) ~ (m_i/m_e) sigma_H is plausible. The three-phase decomposition is useful, and the paper is honest about many of its own caveats.

The central simulation claim—cyclic eruptions with an advection phase, a reconnection-regulated phase, and a flaring phase—holds up. I see multiple cycles, a sigma0 scan, and a consistent saturation at Phi_H/sqrt(Mdot) ~ 60. I'm not independently re-running the code, and no code or data are shipped, so I'm trusting the runs as reported. That's a limitation but not a flaw in the science.

The soft spots are real but not fatal. The outer boundary is the main one. The box is finite, with injection at r_inj ~ 29–30 r_g and a density floor, and the Bondi radius lies outside. The reported box-size test that shrinks r_max by 2/3 breaks the saturation, which shows the reservoir matters. The larger-box run is a single test and does not report the recurrence time or isolate the injection radius. So the claim that T_cyc ~ 10^3 t_g and phi ~ 60 are intrinsic is not fully established. The paper does flag the recurrence time as speculative, but the saturation value is used directly in the Sgr A* luminosity estimate, and that deserves a firmer caveat.

Second, the analytic scalings in Section 5 are calibrated to the same simulations. Equation (28) plugs in measured r_X and V_acc and assumes beta_rec = 0.1; Equation (33) is fit to the decay and then used to infer beta_rec ~ 0.05. That is acceptable as a consistency check, but the paper should not present it as a parameter-free prediction. Reframing those as consistency checks would strengthen the paper.

Minor issues: no convergence plots for the spectral indices, and the extrapolation from m_i/m_e = 256 to 1836 is plausible but untested. The Sgr A* applications are speculative but clearly labeled as such, so no complaint there.

This deserves a serious referee. It is a solid simulation paper with a new kinetic result, and the weaknesses are addressable in revision. I'd bring it to a reading group for the mass-ratio result alone.

Referee Report

3 major / 4 minor

Summary. This manuscript presents 2D GRPIC simulations of collisionless spherical accretion onto a Schwarzschild black hole in an initially vertical magnetic field, for pair and electron-ion plasmas. The central result is a robust three-phase accretion-eruption cycle: (I) quasi-linear growth of horizon magnetic flux in an advection-dominated regime; (II) a reconnection-regulated phase in which Φ_H/√Ṁ saturates near ≈60; and (III) a large-scale reconnection event that expels flux, accelerates particles, and resets the cycle. The authors test robustness against σ0, box size, particle number, and mass ratio, and propose analytic scalings for the flux saturation and exponential decay. They apply the model to Sgr A* flares and isolated stellar-mass black holes.

Significance. If correct, the paper provides a first-principles kinetic demonstration that magnetized spherical accretion is cyclic and that reconnection simultaneously regulates flux accumulation and produces nonthermal particles, with direct relevance to Sgr A* flares. The simulations are expensive and long, and the paper reports multiple cycles, a quasi-steady state, and useful numerical tests (including N_PPC up to 50 and a box-size variation). The three-phase picture is likely to be influential even if the analytic scalings are not fully predictive. However, the analytic framework in §5 is calibrated to the simulation, and the boundary-sensitivity of the recurrence time and saturation value is not fully demonstrated.

major comments (3)

[§3.3, §6.1] The box-size convergence test is incomplete for the claims built on T_cyc and φ. Reducing r_max by 2/3 lowers max Φ_H and prevents σ_H from reaching 10, while increasing r_max by 5/3 leaves max Φ_H and σ_H unchanged; however, the recurrence time T_cyc and the saturated value Φ_H/√Ṁ are not reported for the enlarged box. Because r_pml and r_inj move with r_max, this test does not separate the reservoir size from the injection radius. The abstract's 'self-sustaining' statement, the Sgr A* recurrence estimate (~10^3 t_g), and the luminosity estimate (Eq. 34) all assume these quantities are intrinsic rather than set by the finite numerical reservoir. Please report T_cyc and φ for the larger run, or otherwise demonstrate that the cycle period and saturation are converged with respect to the outer boundary and injection shell.
[§5.2, Eqs. (27)-(28)] The derivation of Φ_H/√Ṁ≈60 is not a parameter-free prediction. Equation (27) contains the measured quantities r_X≈3r_g and V_acc≈0.1 and assumes β_rec≈0.1. In addition, the derivation assumes Ṁ(r_H)~Ṁ(r_X)~ρV_acc 4πr_Xδ/α and B_r^up~αΦ_H/(2πr_X^2). The text says this 'fixes the ratio', but the value 60 is calibrated to the simulation. Please reframe this as an order-of-magnitude consistency check and clearly distinguish measured inputs from assumed ones.
[§5.3, Eq. (33)] The exponential decay law is not tested independently: the decay rate αβ_rec/r_X is fitted to the simulated Φ_H and then used to infer β_rec≈0.05. Combined with the assumed β_rec≈0.1 in §5.2, this is an internal calibration exercise. The statement 'matches the behavior in Fig. 3' is therefore not a validation of the model. Please state explicitly that the decay constant is fit, and discuss what would falsify the model (e.g., checking whether the simulated flux loss is indeed governed by the local reconnection rate at r_X).

minor comments (4)

[§4.2, §6.1] The extrapolation from m_i/m_e=256 to 1836 (γ~10^4) is a factor ~7 in mass ratio; the uncertainty in the electron cutoff and the photon energy estimate (Eq. 35) should be stated more prominently.
[§2] The density floor and injection shell are described, but the rate at which plasma is injected (or the effective reservoir mass) is not given; a sentence on how Ṁ0 is set would help reproducibility.
[Fig. 3] The phase-averaged particle distributions in the top panel would benefit from error bars or shaded ranges, given the strong variability noted in the text.
[§5.4] The critical opening angles θ_open=0.1 and 0.05 are motivated by β_rec=0.1, which is itself inferred in §5.3. The phase-transition criterion is therefore partly circular; this should be acknowledged in the text.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the GRPIC simulation is self-contained; Section 5 is an explicitly simulation-calibrated consistency framework, not an independent first-principles prediction.

full rationale

The central claim is the kinetic simulation itself, which is self-contained and not derived from the analytic model. The analytic Section 5 is explicitly framed as "a simplified model to construct a quantitative physical framework based on the simulation results" (Section 1), and each of its apparent predictions uses simulation-measured inputs: the velocity comparison in Section 5.1 plugs in "˙Φ_H≈0.20, measured during phase one"; the saturation value in Section 5.2 uses "r_X ∼3r_g, where the accretion velocity is measured to be V_acc∼0.1" combined with an assumed "β_rec∼0.1"; and Section 5.3 explicitly "fit[s] a decay rate" to the simulation before inferring β_rec≈0.05. These are calibrations or consistency checks, not claims that the framework independently predicts the simulation. The application to Sgr A* is hedged: the recurrence-time identification "is still speculative at this stage; more theoretical work is needed to determine the dependence... on flux supply from larger scales." Self-citations (e.g., Cerutti et al. 2014; Crinquand et al. 2021) support standard plasmoid-reconnection arguments, but the decay law is derived in the text and the plasmoid threshold is also cited to external works, so no load-bearing step reduces to a self-citation. The boundary-sensitivity limitation raised in Section 3.3 is a numerical robustness concern, not a circularity. No step is equivalent to its inputs by construction.

Assumptions & free parameters 7 free parameters · 8 assumptions · 0 invented entities

The simulation is a numerical experiment; the central observation is not circular. However, the analytic model uses measured/assumed inputs (β_rec, r_X, V_acc, fitted decay rate), so its 'predictions' are consistency checks. No new particles or forces are introduced.

free parameters (7)

Reconnection rate β_rec = 0.1 (assumed); 0.05 inferred from fit
Used in Eq. (27)/(28) to reproduce Φ_H/√Ṁ≃60 and in Eq. (33) for the eruption decay; β_rec=0.1 is taken from flat-space expectations, not independently measured in this setup.
X-point radius r_X = ~3 r_g (measured)
Plugged into force-balance Eq. (27) to obtain 60; also used to convert the fitted decay rate to β_rec. It is measured from one run, not predicted.
Accretion velocity V_acc = ~0.1 (measured in phase II)
Input to Eq. (27); measured from the simulation, not derived.
Normalized temperature θ0 = k_B T0/m_i = 1/30
Chosen to satisfy ϑ0 ≤ 2 r_g/r_pml and keep the Bondi radius outside the box; affects the infall speed and Bondi scaling.
Initial magnetization σ0 = 0.05 (reference pair run); 0.3 (mass-ratio runs); explored 0.03–0.3
Chosen as a control parameter; the central dynamics are claimed robust, but the maximum Φ_H depends on σ0.
Outer boundary radius r_pml = ≈30 r_g (r_max≈33 r_g)
The box size determines the available flux reservoir; reducing r_max by 2/3 changes the saturation behavior (§3.3).
Eruption decay rate αβ_rec/r_X = 0.013 t_g^-1 (fitted)
The exponential decay law Eq. (33) is fit to the simulated Φ_H decay; the inferred β_rec=0.05 is then interpreted rather than predicted.

assumptions (8)

standard math 3+1 ideal plasma/field evolution in Kerr-Schild coordinates with FIDO fields; Maxwell equations; Wald solution for initial uniform field.
Basis of the Zeltron GRPIC formulation (Komissarov 2004; Wald 1974).
domain assumption 2D axisymmetry with no toroidal dynamics; eruptions are simultaneous in all azimuths.
The simulations are 2D; the authors note this is in tension with Sgr A* hotspot observations and may impose artificial cutoffs (§6.1).
domain assumption Continuous injection of fresh thermal plasma at the outer shell with density floor n≥n0; absorbing outer layer for waves and particles.
Keeps mass and magnetic flux supply fixed; box-size dependence in §3.3 shows the boundary can affect saturation.
domain assumption Ambient plasma is collisionless; Coulomb collisions are negligible.
Justifies the kinetic PIC approach for radiatively inefficient accretion flows (§1).
domain assumption Resolving the Debye length at the outer boundary with logarithmic radial compression keeps all kinetic scales resolved.
Needed to trust the PIC results; no convergence plots are shown, though N_PPC and box-size tests are described (§2, §3.3).
ad hoc to paper Flux loss at r_X is transferred to the horizon much faster than it reconnects, so dΦ_H/dt = −(αβ_rec/r_X)Φ_H.
Assumed in §5.3 to derive the exponential decay law; not proven, though plausible if the Alfvén transit time is short.
ad hoc to paper A significant part of the accretion flow passes through the equatorial current layer, so Ṁ(r_H)~Ṁ(r_X)~ρV_acc 4πr_Xδ/α.
Needed to replace ρ with Ṁ in the force balance; the authors say it was verified in the simulations but do not quantify it (§5.2).
domain assumption Relativistic collisionless reconnection proceeds at β_rec≈0.1 in flat spacetime, and tearing thresholds are set by plasmoid-chain physics.
Used to interpret the phase transitions and the measured flux decay; drawn from slab-geometry reconnection literature, not from the global curved-spacetime setup.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Emergence of cyclic flux eruptions in kinetic simulations of magnetized spherical accretion onto a Schwarzschild black hole." pith.science (2026). https://pith.science/paper/UGVDV4LM

@misc{pith2026260204519,
  author       = {Pith},
  title        = {Pith review of: Emergence of cyclic flux eruptions in kinetic simulations of magnetized spherical accretion onto a Schwarzschild black hole},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UGVDV4LM}},
  note         = {Machine review of arXiv:2602.04519}
}

read the original abstract

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.

Figures

Figures reproduced from arXiv: 2602.04519 by the authors.

**Figure 1.** Sketch of the initial GRPIC setup representing zero-angularmomentum accretion onto a central black hole (labeled BH) in the center. On the left, we show in cyan the initial cloud of plasma, at a density n0 and temperature T0, and in pink the injection area of fresh plasma. On the right, the initial uniform magnetic field lines are shown in blue, and the limit of the outer matching layer in red. both with equal num… view at source ↗

**Figure 2.** Top half: Time evolution of: the horizon magnetic flux, ΦH; the horizon accretion rate, M˙ , normalized by M˙ 0 = 4πmin0r 2 injvth; the normalized horizon magnetic flux, ΦH/ √ M˙ ; and the total dissipation rate, E˙, normalized by M˙ 0. Bottom half: Simulation snapshots at t/tg = 2689, 3286, 3596, and 3648. Each snapshot shows the number density n (top left), the electromagnetic energy dissipation J · E (top right),… view at source ↗

**Figure 3.** Top panel: Particle energy distributions at different stages of the eruption cycle. Transparent lines show instantaneous distributions; opaque lines represent distributions time-averaged over each phase. Phases one, two, and three are color-coded, respectively, as green, blue, and red. Bottom panel: Time evolution of ΦH over one eruption cycle. Shaded regions (green, blue, red) indicate the time intervals during wh… view at source ↗

Figures from the paper (6 more)

**Figure 5.** Figure 5 [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]

**Figure 4.** Figure 4: Top: map of the azimuthal electric current carried by the electrons (J − ϕ ). Center: map of the azimuthal electric current carried by the ions (J + ϕ ). Bottom: map of the auxiliary azimuthal magnetic field Hϕ. Magnetic field lines are represented in solid black line…

**Figure 6.** Figure 6: Time-averaged energy distributions of the particles during the erupting phase for all electron-ion simulations (including the reference mi = me run). Ions are represented as dashed lines and electrons as solid lines. this case, the electron energy distribution should t…

**Figure 7.** Figure 7: Profiles of V r measured from the simulation presented in Section 3. Solid lines denote time-averages over the intervals defined in [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]

**Figure 8.** Figure 8: Traces of δ(r)/(r−rH) at specific times within (transparent lines), and time-averaged over (solid dashed lines), each phase. When δ/(r − rH) = 0.1, the magnetic field above and below the equatorial current sheet opens at the same angle as the separatrix field lines ne…

**Figure 9.** Figure 9: Conceptual three-phase accretion cycle. Transitions between phases coincide with critical opening angles of the magnetic field lines about the equator. Blue and green field lines are accreted during phases one and two, respectively; black field lines never reach the bl…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

63 extracted references

[1]

& Kamionkowski, M

Agol, E. & Kamionkowski, M. 2002, MNRAS, 334, 553

2002
[2]

K., Bautz, M

Baganoff, F. K., Bautz, M. W., Brandt, W. N., et al. 2001, Nature, 413, 45

2001
[3]

2018, ApJ, 862, 80

Ball, D., Sironi, L., & Özel, F. 2018, ApJ, 862, 80

2018
[4]

Bisnovatyi-Kogan, G. S. & Ruzmaikin, A. A. 1974, Ap&SS, 28, 45

1974
[5]

Bisnovatyi-Kogan, G. S. & Ruzmaikin, A. A. 1976, Ap&SS, 42, 401

1976
[6]

2019, ARA&A, 57, 467

Blandford, R., Meier, D., & Readhead, A. 2019, ARA&A, 57, 467

2019
[7]

Blandford, R. D. & Znajek, R. L. 1977, MNRAS, 179, 433

1977
[8]

2021, Phys

Bransgrove, A., Ripperda, B., & Philippov, A. 2021, Phys. Rev. Lett., 127, 055101

2021

Show all 63 references

[9]

R., Uzdensky, D

Cerutti, B., Werner, G. R., Uzdensky, D. A., & Begelman, M. C. 2014, ApJ, 782, 104

2014
[10]

2023, ApJ, 959, 122

Chernoglazov, A., Hakobyan, H., & Philippov, A. 2023, ApJ, 959, 122

2023
[11]

2021, A&A, 650, A163

Crinquand, B., Cerutti, B., Dubus, G., Parfrey, K., & Philippov, A. 2021, A&A, 650, A163

2021
[12]

2005, MNRAS, 360, L55

Cuadra, J., Nayakshin, S., Springel, V ., & Di Matteo, T. 2005, MNRAS, 360, L55

2005
[13]

2020, MNRAS, 497, 4999

Dexter, J., Tchekhovskoy, A., Jiménez-Rosales, A., et al. 2020, MNRAS, 497, 4999

2020
[14]

2010, ApJ, 725, 450 Event Horizon Telescope Collaboration, Akiyama, K., Alberdi, A., et al

Dodds-Eden, K., Sharma, P., Quataert, E., et al. 2010, ApJ, 725, 450 Event Horizon Telescope Collaboration, Akiyama, K., Alberdi, A., et al. 2024a, ApJ, 964, L26 Event Horizon Telescope Collaboration, Akiyama, K., Alberdi, A., et al. 2024b, ApJ, 964, L26 Event Horizon Telescop...

2010
[15]

2023, Phys

Galishnikova, A., Philippov, A., Quataert, E., et al. 2023, Phys. Rev. Lett., 130, 115201

2023
[16]

2010, Reviews of Modern Physics, 82, 3121

Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Reviews of Modern Physics, 82, 3121

2010
[17]

2003, Nature, 425, 934

Genzel, R., Schödel, R., Ott, T., et al. 2003, Nature, 425, 934

2003
[18]

M., Wright, S

Ghez, A. M., Wright, S. A., Matthews, K., et al. 2004, ApJ, 601, L159 GRA VITY Collaboration, Abuter, R., Amorim, A., et al. 2021, A&A, 654, A22 GRA VITY Collaboration, Abuter, R., Amorim, A., et al. 2018, A&A, 618, L10 GRA VITY Collaboration, Bauböck, M., Dexter, J., et al. 2...

2004
[19]

2014, Phys

Guo, F., Li, H., Daughton, W., & Liu, Y .-H. 2014, Phys. Rev. Lett., 113, 155005

2014
[20]

Igumenshchev, I. V . & Narayan, R. 2002, ApJ, 566, 137

2002
[21]

C., Lowell, B., et al

Jacquemin-Ide, J., Begelman, M. C., Lowell, B., et al. 2025, arXiv e-prints, arXiv:2510.25842

2025
[22]

& Daughton, W

Ji, H. & Daughton, W. 2011, Physics of Plasmas, 18, 111207

2011
[23]

D., Fish, V

Johnson, M. D., Fish, V . L., Doeleman, S. S., et al. 2015, Science, 350, 1242

2015
[24]

Komissarov, S. S. 2004, MNRAS, 350, 427

2004
[25]

Lam, C. Y . & Lu, J. R. 2023, ApJ, 955, 116

2023
[26]

Y ., Lu, J

Lam, C. Y ., Lu, J. R., Udalski, A., et al. 2022, ApJ, 933, L23

2022
[27]

& Waxman, E

Loeb, A. & Waxman, E. 2007, J. Cosmology Astropart. Phys., 2007, 011

2007
[28]

F., Cowley, S

Loureiro, N. F., Cowley, S. C., Dorland, W. D., Haines, M. G., & Schekochihin, A. A. 2005, Phys. Rev. Lett., 95, 235003

2005
[29]

F., Schekochihin, A

Loureiro, N. F., Schekochihin, A. A., & Cowley, S. C. 2007, Physics of Plasmas, 14, 100703

2007
[30]

Markoff, S., Falcke, H., Yuan, F., & Biermann, P. L. 2001, A&A, 379, L13

2001
[31]

R., Bosch-Ramon, V ., Vieyro, F

Martinez, J. R., Bosch-Ramon, V ., Vieyro, F. L., & del Palacio, S. 2025, A&A, 700, A49

2025
[32]

2018, MNRAS, 475, 1251

Matsumoto, T., Teraki, Y ., & Ioka, K. 2018, MNRAS, 475, 1251

2018
[33]

2025, MNRAS, 541, 3709

Mereghetti, S., Sidoli, L., Ponti, G., & Treves, A. 2025, MNRAS, 541, 3709

2025
[34]

1975, A&A, 44, 59

Meszaros, P. 1975, A&A, 44, 59

1975
[35]

2009, ApJ, 694, L87

Meyer, L., Do, T., Ghez, A., et al. 2009, ApJ, 694, L87

2009
[36]

Narayan, R. & Yi, I. 1994, ApJ, 428, L13

1994
[37]

2020, A&A, 638, A94

Olejak, A., Belczynski, K., Bulik, T., & Sobolewska, M. 2020, A&A, 638, A94

2020
[38]

2019, Phys

Parfrey, K., Philippov, A., & Cerutti, B. 2019, Phys. Rev. Lett., 122, 035101

2019
[39]

2006, ApJ, 643, 1011

Paumard, T., Genzel, R., Martins, F., et al. 2006, ApJ, 643, 1011

2006
[40]

R., et al

Ponti, G., De Marco, B., Morris, M. R., et al. 2015, MNRAS, 454, 1525

2015
[41]

2017, MNRAS, 468, 2447

Ponti, G., George, E., Scaringi, S., et al. 2017, MNRAS, 468, 2447

2017
[42]

2008, A&A, 488, 549

Porquet, D., Grosso, N., Predehl, P., et al. 2008, A&A, 488, 549

2008
[43]

Porth, O., Mizuno, Y ., Younsi, Z., & Fromm, C. M. 2021, MNRAS, 502, 2023

2021
[44]

Quataert, E., Dorland, W., & Hammett, G. W. 2002, ApJ, 577, 524

2002
[45]

M., Quataert, E., & Stone, J

Ressler, S. M., Quataert, E., & Stone, J. M. 2018, MNRAS, 478, 3544

2018
[46]

M., Quataert, E., White, C

Ressler, S. M., Quataert, E., White, C. J., & Blaes, O. 2021, MNRAS, 504, 6076

2021
[47]

2022, ApJ, 924, L32

Ripperda, B., Liska, M., Chatterjee, K., et al. 2022, ApJ, 924, L32

2022
[48]

C., Anderson, J., Casertano, S., et al

Sahu, K. C., Anderson, J., Casertano, S., et al. 2022, ApJ, 933, 83

2022
[49]

Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337

1973
[50]

Shapiro, S. L. & Teukolsky, S. A. 1983, Black holes, white dwarfs and neutron stars. The physics of compact objects

1983
[51]

& Tanuma, S

Shibata, K. & Tanuma, S. 2001, Earth, Planets and Space, 53, 473

2001
[52]

A., & Giannios, D

Sironi, L., Uzdensky, D. A., & Giannios, D. 2025, ARA&A, 63, 127

2025
[53]

Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS, 418, L79

2011
[54]

A., Loureiro, N

Uzdensky, D. A., Loureiro, N. F., & Schekochihin, A. A. 2010, Phys. Rev. Lett., 105, 235002 von Fellenberg, S. D., Roychowdhury, T., Michail, J. M., et al. 2025, ApJ, 979, L20 V os, J., Cerutti, B., Mo´scibrodzka, M., & Parfrey, K. 2025, Phys. Rev. Lett., 135, 015201

2010
[55]

Wald, R. M. 1974, Phys. Rev. D, 10, 1680

1974
[56]

Werner, G. R. & Uzdensky, D. A. 2024, ApJ, 964, L21

2024
[57]

R., Uzdensky, D

Werner, G. R., Uzdensky, D. A., Begelman, M. C., Cerutti, B., & Nalewajko, K. 2018, MNRAS, 473, 4840

2018
[58]

R., Uzdensky, D

Werner, G. R., Uzdensky, D. A., Cerutti, B., Nalewajko, K., & Begelman, M. C. 2016, ApJ, 816, L8

2016
[59]

2022, A&A, 665, L6

Wielgus, M., Moscibrodzka, M., V os, J., et al. 2022, A&A, 665, L6

2022
[60]

2018, ApJ, 863, 15

Witzel, G., Martinez, G., Hora, J., et al. 2018, ApJ, 863, 15

2018
[61]

Zelenyi, L. M. & Krasnoselskikh, V . V . 1979, Soviet Ast., 23, 460

1979
[62]

& Hoshino, M

Zenitani, S. & Hoshino, M. 2008, ApJ, 677, 530

2008
[63]

2021, ApJ, 922, 261 Article number, page 12 of 12

Zhang, H., Sironi, L., & Giannios, D. 2021, ApJ, 922, 261 Article number, page 12 of 12

2021

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.

[1] [1]

& Kamionkowski, M

Agol, E. & Kamionkowski, M. 2002, MNRAS, 334, 553

2002

[2] [2]

K., Bautz, M

Baganoff, F. K., Bautz, M. W., Brandt, W. N., et al. 2001, Nature, 413, 45

2001

[3] [3]

2018, ApJ, 862, 80

Ball, D., Sironi, L., & Özel, F. 2018, ApJ, 862, 80

2018

[4] [4]

Bisnovatyi-Kogan, G. S. & Ruzmaikin, A. A. 1974, Ap&SS, 28, 45

1974

[5] [5]

Bisnovatyi-Kogan, G. S. & Ruzmaikin, A. A. 1976, Ap&SS, 42, 401

1976

[6] [6]

2019, ARA&A, 57, 467

Blandford, R., Meier, D., & Readhead, A. 2019, ARA&A, 57, 467

2019

[7] [7]

Blandford, R. D. & Znajek, R. L. 1977, MNRAS, 179, 433

1977

[8] [8]

2021, Phys

Bransgrove, A., Ripperda, B., & Philippov, A. 2021, Phys. Rev. Lett., 127, 055101

2021