Pith. sign in

REVIEW 3 major objections 4 minor 17 references

Shock acceleration in vortex driven magnetic fields of black holes

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

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

Pith's one-line read Shock acceleration in vortex-driven black-hole magnetospheres can push protons to several hundred PeV and cap electrons near 120 GeV.

desk verdict The Fermi-acceleration balance is standard and the energy ceilings are concrete new numbers, but the 400 PeV endpoint rests on an unobserved, self-cited vortex-field bound that the paper takes as a given. read the letter →

arxiv 2607.20209 v1 pith:QM2LO2QS submitted 2026-07-22 astro-ph.HE

Zaza N. Osmanov This is my paper · ORCID
classification astro-ph.HE
keywords Fermiaccelerationshocksupermassiveblackholesvortex-drivenmagneticfieldssynchrotronlossescosmicraysinverseComptonscatteringPeVatrons
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 examines whether first-order Fermi acceleration at relativistic shocks in a supermassive black hole's magnetosphere can produce very-high-energy particles when the magnetic field is vortex-driven. Balancing the shock energy-gain rate against synchrotron and inverse-Compton losses, the author finds that protons around a 10^8-solar-mass black hole can be accelerated to roughly 100 TeV–400 PeV, while electrons are capped near 120 GeV. The reason this matters is that it offers a concrete hadronic route from black-hole magnetospheres to the highest-energy cosmic rays, with a predicted spectral cutoff that observations could test. The calculation's reach is conditional: the field strength B = αB_m is not measured but taken from earlier theoretical proposals, and the lower bound that sets the 400 PeV endpoint comes from a companion work. A sympathetic reader would describe the paper as a direct estimate of what vortex-generated fields would do to shock acceleration, with the magnetic field itself as the main uncertainty.

What carries the argument

The central object is the synchrotron-limited Fermi-I balance at a relativistic shock. The acceleration rate is dE/dt ≈ E/(η r_L/c) ≈ eBc/η, and the average synchrotron power is P_syn ≈ (4e^4B^2E^2)/(9m^4c^7); equating these yields the maximum energy, which scales as (m/m_p)^2 sqrt(M/(α M_⊙)). The vortex-driven field enters through B = αB_m, where B_m ≃ 2.4×10^11 (10^8 M_⊙/M) G is the upper bound obtained by equating magnetic-field energy in a region of size r_g with the black hole's rest energy. Because the synchrotron power grows as B^2 while the acceleration rate grows only as B, the maximum energy falls as α increases; this monotonic decrease is what makes the weaker-field part of the vo

What would settle it

For a 10^8-solar-mass black hole, determine the magnetic-field strength in the shock zone (e.g. via Faraday rotation or synchrotron spectral mapping) and search for the predicted proton cutoff near 100 TeV–400 PeV in neutrino or gamma-ray emission. Finding B ≈ 10^4 G with no cutoff in that band, or a cutoff energy that disagrees with the balance dE/dt ≈ eBc/η versus synchrotron power, would falsify the claim.

Watch

Extended reading notes

Core claim

The paper's central claim is that in the magnetosphere of a supermassive black hole threaded by a vortex-driven field, first-order Fermi shock acceleration is limited almost entirely by synchrotron cooling. Writing the acceleration rate as dE/dt ≈ eBc/η and equating it to the average synchrotron power gives a maximum energy that scales as (m/m_p)^2 sqrt(M/(α M_⊙)), where B = αB_m and B_m is the theoretical upper field set by equating field energy with the black hole's rest energy. For a 10^8-solar-mass black hole, this places protons between roughly 100 TeV and 400 PeV and electrons between roughly 40 MeV and 120 GeV. The paper also argues that inverse Compton scattering, in both Thomson and

Load-bearing premise

The calculation stands or falls on the assumption that the vortex-generated field B = αB_m — with α taking values as small as 10^-6 — actually fills and stays roughly uniform in the shock region; the paper imports this field from earlier theoretical proposals rather than from observation, and the value α=10^-6 that sets the 400 PeV endpoint comes from a companion study.

Editorial extensions

If this is right

  • If the argument is correct, a 10^8-solar-mass black hole with a vortex field is a PeVatron: protons emerge at 100 TeV–400 PeV, placing supermassive-black-hole magnetospheres among the sources that can feed the observed high-energy cosmic-ray flux.
  • The predicted proton and electron synchrotron cutoffs are sharp and depend on α, so a measured cutoff in gamma-ray or neutrino spectra would directly constrain the magnetic-field fraction in the shock region.
  • Because inverse Compton losses are negligible, most of the electron energy is radiated as synchrotron emission rather than Compton-scattered to very high energies; the model therefore predicts a distinctive spectral shape with little cascade signature.
  • The result implies a non-monotonic role for strong fields: raising the field from α=10^-6 to α=1 lowers the maximum energy, so the strongest-field subregions may be the least efficient at making the highest-energy particles.
  • The same balance scales with black-hole mass and particle mass, so the energy ceiling for lighter black holes or heavier ions can be read off directly from the same formula.

Reading between the lines

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

  • A natural observational test the paper does not spell out: search nearby supermassive-black-hole jets for a hadronic cutoff in the 100 TeV–400 PeV band; an absence of such a cutoff would indicate either that the vortex field does not reach the shock or that the field is not uniform.
  • The assumption of a spatially uniform field is likely too simple; realistic vortex fields decay with distance, so one would expect a range of α values along the shock and a softened, two-component cutoff rather than the single sharp value plotted in the paper's figures.
  • A skeptical reading would stress that the 400 PeV endpoint rests on the theoretical lower bound α=10^-6 from the companion work; the paper's own formula is consistent, but the wide energy range is an output of that bound rather than an independent measurement.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper studies first-order Fermi acceleration at relativistic shocks in supermassive black hole magnetospheres, assuming a 'vortex-driven' magnetic field of the form B = αB_m, with B_m given by Eq. (1) and α a free parameter. The author balances the acceleration rate, Eq. (4), against synchrotron losses, Eq. (5), to derive a maximum particle energy, Eq. (6). For M = 10^8 M_sun and α in [10^-6, 1], the paper claims protons can reach 100 TeV–400 PeV and electrons 40 MeV–120 GeV. Inverse Compton losses are argued to be negligible in both Thomson and Klein–Nishina regimes via Eqs. (7)–(9). The manuscript is a short analytical contribution with no new simulations or observations; its central claim is a prediction conditional on the assumed vortex-field strength and on the shock acceleration model.

Significance. If the vortex-driven magnetic field B = αB_m actually reaches the assumed values in the acceleration zone, the paper identifies a plausible mechanism for producing ultra-high-energy protons (up to hundreds of PeV) in SMBH magnetospheres, which would be relevant for cosmic-ray origin studies. The acceleration-loss balance is standard and the manuscript is transparent enough to allow independent checking, which is a strength. However, the numerical endpoints are dominated entirely by the input parameter α; the α = 10^-6 lower bound is imported from an unpublished, self-authored manuscript (Dvali, Osmanov & Zantedeschi 2025) and is not independently derived or observationally anchored. The central claim is therefore not yet demonstrative, and the paper currently overstates what is established. The core physics is simple and correct in outline, but the manuscript requires correction of a likely algebraic error in Eq. (6) and a more careful framing of the model dependence before the results can be accepted.

major comments (3)
  1. [§2, Eq. (6)] The η dependence and the numerical coefficient in Eq. (6) are inconsistent with the stated balance. Setting dE/dt = eBc/η equal to the synchrotron power P_syn = 4e^4B^2E^2/(9m^4c^7) yields E_max = (3/2) m^2 c^4 / [e^{3/2} sqrt(η B)], i.e. E_max ∝ 1/√η, not √η as printed. Numerically, for protons with M = 10^8 M_sun and B = αB_m, one obtains E_max ≈ (410 TeV)/√η × sqrt(M / (α 10^8 M_sun)), not 590√η. The printed formula overestimates the endpoint by a factor of about η (e.g., a factor ~3 for η = 10) and also disagrees with the abstract's 'several hundred PeV' range. The figures and conclusions were presumably made with a corrected formula, but the manuscript as written contains a load-bearing algebraic error that must be fixed and the numerical results re-stated consistently.
  2. [§2 after Eq. (1); §2, B = αB_m] The central numerical claim depends entirely on the assumption that the magnetic field in the shock region takes the vortex-generated value B = αB_m with α as low as 10^-6. This lower bound is taken from Dvali, Osmanov & Zantedeschi (2025), a submitted manuscript co-authored by the present author. No independent derivation or observational constraint is provided in this paper, and the result scales as α^{-1/2}. An order-of-magnitude error in α changes the energy endpoint by a factor ~√10; if the vortex field does not thread the actual acceleration region, the entire predicted energy range is unsupported. The manuscript should either (i) provide a self-contained derivation or a testable observable consequence of the α = 10^-6 bound, or (ii) clearly frame the results as conditional on the vortex-field hypothesis and avoid the categorical 'demonstrated' language in the abstract and conclusi
  3. [§2, Eqs. (7)–(9)] The inverse Compton analysis is incomplete. The Thomson-regime check leading to the conclusion that the Thomson regime 'violates' the applicability condition is not shown in detail; the condition itself is written as 'E_IC kT / (m^2 c^4) << 1', which appears dimensionally inconsistent and should involve the electron mass and the photon energy in the electron rest frame. The Klein–Nishina ratio τ_KN/τ_acc in Eq. (9) is stated without derivation. In particular, the accretion temperature T in Eq. (8) depends on u∞ (the ambient sound speed) and n∞; only n∞ = 1 cm^-3 is specified, while u∞ is left undefined. Since the claim that synchrotron is the sole limiting process for electrons relies on this comparison, the derivation should be completed and u∞ chosen (or the result shown to be insensitive to it). Until then, the electron energy bounds are not fully established.
minor comments (4)
  1. [§2, Eq. (6)] The notation '2e√ηαeB' is ambiguous. Presumably it means 2 e sqrt(η α e B) or equivalently 2 e^{3/2} sqrt(η α B), but as printed it could be misread. Use parentheses or a clearer radical sign.
  2. [§2, after Eq. (1)] The phrase 'can be as much as 10^6 times smaller' is awkward; it should read 'can be as small as 10^-6 times the upper limit'.
  3. [§2, Fig. 2] For electrons at α = 10^-6, η = 1, Eq. (6) (even with the corrected 1/√η form) gives ~175 GeV, while Fig. 2 shows a maximum of ~120 GeV. The figure axes and the text should be checked for consistency.
  4. [References] The in-text citation '(Shapiro 1983)' should be '(Shapiro & Teukolsky 1983)'. Also, some self-citations (Osmanov 2025, 2026) appear to be to papers in press; please confirm that all references are complete and appropriately attributed.

Circularity Check

1 steps flagged · score 4.0 of 10

The claimed proton/electron energy range rests largely on a load-bearing, self-authored submitted paper for the vortex-field lower bound α=10^-6; the radiative-balance algebra itself is independent.

  1. self citation load bearing [Section 2, text after Eq. (6) (arXiv:2607.20209)]
    "Theoretical analysis shows that the minimum magnetic field strength can be as much as 10 6 (α= 10 −6) times smaller than the aforementioned theoretical upper limit (Dvali et al. 2025)."

    The headline result—'protons can be accelerated to energies of several hundred PeV' for M=10^8 M_sun—is obtained from Eq. (6), where E_max ∝ (M/(α 10^8 M_sun))^{1/2}; the value α=10^-6 is the input that produces the several-hundred-PeV endpoint. That lower bound is not derived or tested in this paper; it is imported from Dvali, Osmanov & Zantedeschi 2025, a submitted paper co-authored by the present author. The manuscript supplies no independent derivation, observational check, or error estimate for this α value, so the central quantitative claim depends on a load-bearing self-citation chain rather than on an externally verified result.

full rationale

The core derivation is not circular by construction: Eq. (6) follows from a standard synchrotron-loss balance against Fermi-I acceleration, and no data are fitted in this paper. If the vortex field B=αB_m is granted, the resulting energy limits are straightforward algebra. However, the numerical interval claimed in the abstract and conclusions ('100 TeV to 400 PeV', 'several hundred PeV') depends critically on the α=10^-6 lower bound, which is taken from a submitted paper co-authored by the present author and is not independently substantiated here. Because E_max scales as α^{-1/2}, an order-of-magnitude uncertainty in α changes the predicted energy by a factor ~√10, making the self-citation load-bearing for the headline result. This warrants a moderate circularity score under the self-citation-load-bearing rule, but not a high score, since the remaining radiative and acceleration physics is independent and the paper does not redefine the field in terms of the predicted energies.

Assumptions & free parameters 4 free parameters · 6 assumptions · 0 invented entities

The calculation imports a non-observed magnetic field (B=αB_m) and a standard-but-unvalidated acceleration time-scale; it contributes no new entities. The main free inputs are α and η, plus an unspecified u∞ for the inverse-Compton side check.

free parameters (4)
  • α (magnetic field scaling) = adopted range 10^-6 to 1
    B=αB_m; the claimed proton energy range 100 TeV–400 PeV is driven entirely by scanning this parameter over the vortex-field uncertainty; no derivation of α is given.
  • η (acceleration efficiency) = 1, 3, 10 in figures
    t_acc = η r_L/c; chosen by hand to show curves; no physical model selects η, and as printed Eq. (6) has an incorrect η dependence.
  • u∞ (ambient sound speed) = unspecified
    Enters Eq. (8) via the Bondi accretion rate and therefore sets the numerical prefactor in Eq. (9); no value is given anywhere, so the IC-negligibility claim cannot be checked.
  • n∞ (ambient density) = 1 cm^-3
    Chosen as a 'natural astrophysical parameter' for the IC/KN side calculation; not used in the final synchrotron-limited energy.
assumptions (6)
  • domain assumption Vortex-driven magnetic field exists with B=αB_m, B_m=c^4/(M G^{3/2}) ≈ 2.4e11 (1e8 M_sun/M) Gauss
    Taken wholesale from Dvali et al. 2021 and Dvali et al. 2025; not derived or tested here. Eq. (1) and §2.
  • domain assumption Plasma flows near the BH are supersonic and form relativistic shocks; Alfvén speed ≈ c because magnetic energy dominates plasma rest energy
    §2: 'Plasma flows ... supersonic' and the discussion of Eq. (2); the printed comparison nm_p c^2 << B is dimensionally sloppy.
  • domain assumption Fermi-I shock acceleration time-scale is t_acc = η r_L/c with η≥1
    §2, Eqs. (3)–(4); standard Bohm-limit estimate, but no diffusion or shock-compression model is given.
  • standard math Synchrotron power Eq. (5) and the Thomson/KN inverse-Compton cooling formulas are valid in this regime
    Eq. (5), Eq. (7), and the KN expression after Eq. (8); standard results from Rybicki & Lightman and Blumenthal & Gould.
  • domain assumption Maximum energy is found by equating dE/dt with radiative losses; escape and spatial confinement are not limiting
    §2 equilibrium balance; ignores Hillas/escape constraints which could further lower E_max if the Larmor radius exceeds the acceleration region.
  • domain assumption Bondi accretion rate and thin-disk temperature formula Eq. (8) describe the accretion flow
    Eq. (8) after Carroll; u∞ is never specified; used only to argue that inverse Compton cooling is negligible.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Shock acceleration in vortex driven magnetic fields of black holes." pith.science (2026). https://pith.science/paper/QM2LO2QS

@misc{pith2026260720209,
  author       = {Pith},
  title        = {Pith review of: Shock acceleration in vortex driven magnetic fields of black holes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QM2LO2QS}},
  note         = {Machine review of arXiv:2607.20209}
}
abstract

The aim of this work is to investigate particle acceleration by shock waves in the magnetospheres of supermassive black holes, where the magnetic field is vortex-driven. For this purpose, we investigated the first-order Fermi acceleration process at relativistic shocks in the magnetospheres of supermassive black holes. In our analysis, we considered synchrotron radiation and inverse Compton scattering as the principal energy-loss mechanisms limiting the maximum attainable particle energy. We investigated the acceleration of both protons and electrons and demonstrated that, for a typical supermassive black hole with a mass of $10^8$ solar masses, protons can be accelerated to energies of several hundred PeV, whereas electrons can attain energies of up to approximately 120 GeV.

Figures

Figures reproduced from arXiv: 2607.20209 by the authors.

Figure 1
Figure 1. — Maximum attainable energy of protons versus α. The set of parameters is: m = mp, M = 108M⊙. 10 -6 10 -4 10 -2 10 0 , 10 7 10 8 10 9 10 10 10 11 E max (eV) [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. — Maximum attainable energy of electrons [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 1 linked inside Pith

  1. [1]

    WorldvScientific Publishing Co

    Aharonian, F.A., 2004, Very High Energy Cosmic Gamma Radiation - A Crucial Window on the Extreme Universe. WorldvScientific Publishing Co. Pte. Ltd. (2004) Anile A. M., 1989, Relativistic fluids and

    2004

  2. [2]

    Press, Cam- bridge Bell A.R., 1978, MNRAS, 182, 147 Bell A.R., 1978, MNRAS, 182, 443

    Magneto-Fluids, Cambridge Univ. Press, Cam- bridge Bell A.R., 1978, MNRAS, 182, 147 Bell A.R., 1978, MNRAS, 182, 443

    1978

  3. [3]

    & Woltjer, L., 1990, Active Galactic Nuclei (Springer-Verlag Berlin Heidelberg GmbH)

    Blandford, R.D., Netzer, H. & Woltjer, L., 1990, Active Galactic Nuclei (Springer-Verlag Berlin Heidelberg GmbH)

    1990

  4. [4]

    & Gould, R.J., 1970, Rev

    Blumenthal, G.R. & Gould, R.J., 1970, Rev. Mod. Phys., 42, 237-271

    1970

  5. [5]

    & Ostlie, Dale A.,An intro- duction to modern astrophysics and cosmology, Pearson (2010) Dvali G., Kuhnel, F

    Carroll, Bradley W. & Ostlie, Dale A.,An intro- duction to modern astrophysics and cosmology, Pearson (2010) Dvali G., Kuhnel, F. & Zantedeschi, M., 2021, Phys. Rev. L., 12, 061302

    2010

  6. [6]

    & Zantedeschi, M., 2025, (submitted to JCAP) arXiv: 2502.15510

    Dvali, G., Osmanov, Z.N. & Zantedeschi, M., 2025, (submitted to JCAP) arXiv: 2502.15510

    arXiv 2025

  7. [7]

    Fermi, E., 1949, Phys. Rev. 75, 1169 Kim H.B. & Kim J., 2013, Int. J. Mod. Phys. D, 22, 1350045

    1949

  8. [8]

    Kotera, K., & Olinto, A.V., 2011, ARA&A49, 119

    2011

Show all 17 references
  1. [9]

    & Osmanov, Z.N., 2025, MNRAS, 541, 3445

    Nikuradze, N. & Osmanov, Z.N., 2025, MNRAS, 541, 3445

    2025

  2. [10]

    Osmanov, Z.N., 2026, IJMPD, 35, 2650026

    2026

  3. [11]

    Osmanov, Z.N., 2025, ApSS, 370, 116

    2025

  4. [12]

    Osmanov, Z.N., 2021, Galaxies, 9, 6

    2021

  5. [13]

    2023, Galaxies, 4, 226 4

    Osmanov, Z.N., Mahajan, S.M. 2023, Galaxies, 4, 226 4

    2023

  6. [14]

    Phys., 99, 30

    Osmanov, Z., Mahajan, S.; Machabeli, G., 2018, Astropart. Phys., 99, 30

    2018

  7. [15]

    & Duffy, P., 2007, ApSS, 309, 119-125

    Rieger, F.M., Bosch-Ramon, V. & Duffy, P., 2007, ApSS, 309, 119-125

    2007

  8. [16]

    L., & Teukolsky, S

    Shapiro, S. L., & Teukolsky, S. A. 1983, Black

    1983

  9. [17]

    & Lightman, A.P

    Rybicki, G.B. & Lightman, A.P. 2007, Radiative Processes in Astrophysics, Wiley, 2007 This 2-column preprint was prepared with the AAS LATEX macros v5.2. 5

    2007

Pith tools

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