Pith. sign in

REVIEW 3 major objections 6 minor 27 references

Magnetically arrested advective accretion flows and jets/outflows around stellar mass black holes: Explaining hard state ULXs with GRMHD simulations

T0 review · 3 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read The paper claims that magnetically arrested, advective accretion around a rapidly spinning stellar-mass black hole can explain hard-state ULXs, without super-Eddington accretion or intermediate-mass black holes.

desk verdict The GRMHD runs and the spin/magnetic-field trends are credible, but the claim that they explain hard-state ULX luminosities is not supported: the paper equates mechanical outflow power with X-ray luminosity without any radiation physics, and the CGS normalization is chosen to land in the ULX band. read the letter →

arxiv 2411.18599 v1 pith:RAGHEGYE submitted 2024-11-27 astro-ph.HE

Rohan Raha , Banibrata Mukhopadhyay , Koushik Chatterjee This is my paper · ORCID
classification astro-ph.HE
keywords ultraluminousX-raysourceshardstateGRMHDsimulationsmagneticallyarrestedaccretionadvectiveflowsblackholespinjetsandoutflowsEddingtonmagneticfield
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

The paper tries to resolve a puzzle: observed ultra-luminous X-ray sources (ULXs) in the hard X-ray state shine at $10^{39}$--$10^{40}\,\mathrm{erg\,s^{-1}}$, far above the Eddington limit for a stellar-mass black hole, yet the hard state is normally tied to low, sub-Eddington accretion. Its central claim is that a magnetically arrested, advective, sub-Keplerian accretion flow around a rapidly spinning stellar-mass black hole can produce exactly this luminosity without super-Eddington accretion or an intermediate-mass black hole. The evidence comes from axisymmetric GRMHD simulations of a strongly magnetized, optically thin torus around a rotating black hole, which show outflow efficiencies above unity when the spin and magnetic field are high. Converting those efficiencies with a $20\,M_\odot$ black hole accreting at $0.05\,\dot{M}_{\rm Edd}$ places the mechanical outflow power in the observed ULX band; the paper takes that mechanical power as the source of the ULX luminosity, noting that radiation transport is not yet included. If the claim is right, hard-state ULXs are magnetically dominated flows around stellar-mass black holes, removing the need for heavier black holes or extreme accretion rates.

What carries the argument

The load-bearing object is the magnetically arrested advective accretion flow (MA-AAF): an optically thin, advective, sub-Keplerian disk in which advected poloidal magnetic flux accumulates near the black hole until magnetic pressure balances the ram pressure of the inflowing gas. The energy-conversion mechanism is the twisting of these field lines by frame dragging and disk rotation, which launches helical magnetic waves and jets; the paper quantifies the result through the outflow efficiency $\eta=(\dot{M}-\dot{E})/\dot{M}$, where $\dot{M}$ is the mass accretion rate and $\dot{E}$ is the radial energy flux. This efficiency is then combined with the Eddington magnetic field ceiling $B_{\rm Edd}$ and a Blandford-Znajek-type power law $P\propto \phi_{\rm BH}^{2}a^{2}$ to translate simulated field strengths and spin parameters into CGS outflow powers.

What would settle it

A decisive check is to add radiation transport to the simulations (or post-process the snapshots with a radiation code) and compute the emergent spectrum for the same $20\,M_\odot$, $0.05\,\dot{M}_{\rm Edd}$ case: if the radiative luminosity comes out below $10^{39}\,\mathrm{erg\,s^{-1}}$, or if the magnetic flux required to reach the claimed power exceeds the Eddington magnetic field ceiling, the explanation fails. Similarly, if X-ray polarization or spectral observations of a hard-state ULX show no sign of a magnetically arrested geometry, the model would be ruled out.

Watch

Extended reading notes

Core claim

The discovery the paper argues for is that hard-state ULXs can be powered by magnetically arrested advective accretion flows (MA-AAF) around stellar-mass black holes with high spin and strong magnetic fields. In the simulations, poloidal magnetic flux is dragged inward by the accretion flow; when the spin is high ($a=0.998$) and the initial magnetization is strong (plasma-$\beta=0.1$, meaning gas pressure is one-tenth of magnetic pressure), the field forms a magnetic barrier near the horizon, producing a magnetically arrested disk in which matter must diffuse through repeated magnetic barriers. The resulting outflow power grows with both spin and near-horizon magnetic flux, following a Blandford-Znajek-type scaling, and reaches mechanical efficiencies of order unity or higher. For a $20\,M_\odot$ black hole with $\dot{M}=0.05\,\dot{M}_{\rm Edd}$, these efficiencies translate to outflow powers of $10^{39}$--$10^{40}\,\mathrm{erg\,s^{-1}}$, and the simulated magnetic field at the jet base reaches roughly $10^6$ G, consistent with hard-state binaries such as Cygnus X-1. The paper's conclusion is that the combination of high spin and strong magnetic field lets a sub-Eddington stellar-mass accretor radiate at super-Eddington luminosity, explaining the hard-state ULX population without intermediate-mass black holes.

Load-bearing premise

The load-bearing premise is that the mechanical power of the simulated outflows can be equated with the X-ray luminosity observed from ULXs, even though the simulations contain no radiation transport; the conversion also assumes a $20\,M_\odot$ black hole accreting at $5\%$ of Eddington, values chosen to land in the ULX band.

Editorial extensions

If this is right

  • Hard-state ULXs can be understood without super-Eddington accretion rates or intermediate-mass black holes; a roughly $20\,M_\odot$ black hole with high spin and a strong large-scale magnetic field suffices.
  • Outflow power should increase with both black hole spin and near-horizon magnetic flux, so the brightest hard-state ULXs should be the most rapidly spinning and magnetically arrested systems.
  • Magnetic field strengths near the jet base should reach about $10^6$ G for such sources, matching the hard-state X-ray binary Cygnus X-1.
  • The periodic formation and dissipation of magnetic barriers in the magnetically arrested runs provides a natural candidate mechanism for ULX variability.
  • The same magnetically arrested advective flow geometry can be applied to bright hard states of other stellar-mass black hole sources, not only ULXs.

Reading between the lines

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

  • The paper's implicit step is that a large fraction of the mechanical outflow power is eventually radiated as X-rays; if future radiative GRMHD runs show that most of the Poynting flux escapes without thermalizing, the luminosity claim would need to be revised even if the dynamics are correct.
  • Because the key controls are dimensionless (spin and magnetic flux normalized by the accretion rate), the same mechanism should scale up to black holes of other masses, so low-luminosity active galactic nuclei in hard states may be a test bed for this model.
  • A concrete observational test: X-ray polarimetry of a hard-state ULX should show a strong, stable polarization signature if the emission is dominated by ordered magnetic fields near the jet base.
  • The magnetic-barrier cycle seen in the high-spin, strong-field run suggests quasi-periodic flux variations on timescales of thousands of $r_g/c$; checking whether observed ULX variability has such a quasi-periodic component would test the model.
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 / 6 minor

Summary. The paper presents 2D GRMHD simulations of magnetically arrested advective accretion flows around Kerr black holes using the HARMPI code, exploring different black hole spins (a = 0.5, 0.9375, 0.998) and initial plasma beta values (0.1 and 1). The authors compute the outflow power and net efficiency from the stress-energy tensor, find that higher spin and stronger initial magnetic fields produce more powerful outflows, and report magnetic field strengths up to ~10^6 G near the jet base. They then claim that the outflow power reaches levels consistent with hard-state ultraluminous X-ray sources (ULXs) at 10^39-10^40 erg/s, without requiring super-Eddington accretion rates or intermediate-mass black holes, and they attribute the high power to Blandford-Znajek and Blandford-Payne mechanisms.

Significance. If the ULX claim were supported, the paper would offer a novel explanation for hard-state ULXs and would strengthen the case for magnetically arrested accretion flows around stellar-mass black holes. The numerical setup is standard, the outflow-efficiency definition is clear, and the reported spin/magnetic-field trends are plausible and consistent with earlier work. However, the central claim is not established: the simulations contain no radiation physics, and the conversion to CGS units fixes the black hole mass and accretion rate by hand so that the resulting power lands inside the asserted ULX band. The paper's strengths are the careful description of the GRMHD formalism, the time-averaging procedure, and the comparison with a companion BHAC simulation, but these do not compensate for the missing radiation model that the ULX luminosity claim requires.

major comments (3)
  1. [Sections 4 and 5] The paper equates the mechanical outflow power P_out (Equations 19-21) with the observed X-ray luminosities of ULXs, but the GRMHD equations (2)-(15) contain no radiation fields or cooling, as the authors acknowledge in the Introduction ("though no radiation physics included"). The power that reaches large radius in an ideal GRMHD simulation is the total energy flux carried by matter and electromagnetic fields; the fraction that emerges as hard X-rays depends on radiative processes that are not computed. Therefore the statement in Section 5 that "the outflow power achieved in our simulations reaches levels consistent with observed ULXs" is not supported by the simulated quantity.
  2. [Section 4, normalization] The conversion to CGS units fixes MBH = 20 Msun and Mdot = 0.05 Mdot_Edd before evaluating the outflow power. Because P_out = eta * Mdot * c^2 with eta ~ 1.2, this choice yields approximately 1.5e39 erg/s, placing the result inside the claimed 10^39-10^40 erg/s band. The mass accretion rate is not determined by the simulation; it is an input. If Mdot = 0.01 Mdot_Edd were used instead, the power would fall near 3e38 erg/s, below the ULX range. Thus the agreement with the ULX luminosity band is a normalization choice, not a falsifiable prediction, and it cannot support the conclusion that the model "successfully explains" ULX luminosities.
  3. [Section 5] The paper describes the flow as optically thin and advective, which in the standard understanding implies low radiative efficiency. High mechanical outflow power does not translate directly into high luminosity; in advection-dominated flows a large fraction of the dissipated energy can be swallowed by the black hole or converted into kinetic power rather than radiation. Without a radiation model or a radiative efficiency estimate, the claim that "energy stored in strong magnetic fields can generate super-Eddington luminosity" is not demonstrated.
minor comments (6)
  1. [Section 4, first paragraph] The phrase "We have began" should be "We have begun".
  2. [Section 3] The coordinate transformation contains a parameter h that is not defined in the text.
  3. [Section 4, Figure 5 discussion] The Eddington magnetic field is quoted as B_Edd ~ 10^4 G (M/10^9 Msun)^-1/2, but its numerical value for a 20 Msun black hole is not stated; please provide it.
  4. [Section 1] The acronym MA-AAF is introduced but the full phrase "magnetically arrested advective accretion flow" appears only in parentheses; consider spelling it out at first use.
  5. [References] Reference [12] is given as "present volume (2024)" without enough detail; please complete it if permitted.
  6. [Figure 6 caption] The caption appears truncated and does not fully describe the three curves shown in the panel; please complete it.

Circularity Check

2 steps flagged · score 6.0 of 10

The claimed ULX luminosity match is fixed by hand-chosen MBH=20 Msun and Mdot=0.05 Mdot_Edd, and mechanical outflow power is equated to observed X-ray luminosity despite the absence of radiation physics in the GRMHD runs.

  1. fitted input called prediction [Section 4 ('Throughout, when we report the results in CGS units...') and Section 5 ULX conclusion; see Eqs. (20)-(21).]
    "Throughout, when we report the results in CGS units, we consider a stellar mass black hole of mass MBH = 20 M⊙ with total mass accretion rate, ˙M = 0.05 ˙MEdd, where ˙MEdd = LEdd/ηc2 = 1.39 × 1018(MBH/M) g s−1, considering radiative efficiency η = 0.1."

    The only genuinely simulated output is the dimensionless efficiency eta (Eq. 20). The CGS luminosity is obtained by multiplying eta by the hand-chosen Mdot and MBH: with eta ~ 1.2 (Sec. 4), Pout ~ 1.2 x (0.05 x 1.39e18 x 20 g/s) x c^2 ~ 1.5e39 erg/s, inside the claimed 10^39-10^40 band. The GRMHD equations do not predict Mdot; changing it to 0.01 Mdot_Edd or reducing MBH to 3 Msun moves Pout below the ULX band. The 'match' is therefore a normalization chosen to land in the target range, not a prediction.

  2. other [Section 1 parenthetical '(though no radiation physics included)'; Section 5 'The outflow power achieved...' and 'Future work...']
    "we simulate a sub-Eddington, disk-outflow symbiotic model in the advective regime (hence hard state; though no radiation physics included)"

    Equations (2)-(15) are ideal GRMHD with no radiation transport, so the code outputs mechanical outflow power, not X-ray luminosity. The Summary nevertheless states 'The outflow power achieved in our simulations reaches levels consistent with observed ULXs, with luminosities in the range of 10^39-10^40 ergs/s,' identifying Pout with the observed X-ray band. The paper itself defers 'radiative GRMHD simulations' and 'comparing simulated spectra with observed ULX spectra' to future work, confirming that the conversion from mechanical power to the ULX luminosity is assumed, not derived.

full rationale

The spin/magnetic-field trends and the efficiency eta ~ 1.2 are genuine GRMHD outputs and are independently meaningful. The circularity is concentrated in the ULX claim. The CGS luminosities in Figures 1 and 6 and the Section 5 conclusion are produced by inserting MBH=20 Msun and Mdot=0.05 Mdot_Edd (Section 4) into Pout = eta Mdot c^2 (Eq. 20). Since MBH and Mdot are inputs, not outputs of the simulation, and since the runs contain no radiation physics, the quantitative agreement with the 10^39-10^40 erg/s ULX band is a hand-set normalization rather than a falsifiable prediction. The paper's own admission 'though no radiation physics included' and its deferral of radiative GRMHD and spectral comparison to future work make the identification of mechanical outflow power with X-ray luminosity an assumption. This is partial circularity: one load-bearing quantitative claim reduces to chosen inputs, while the qualitative physics results do not. Score 6 rather than higher because the simulation does contain real dynamical content and external consistency checks (e.g., B ~ 10^6 G at the jet base matching Cyg X-1) that are not themselves circular.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The central ULX luminosity estimate depends on two free scalings, MBH and Mdot, plus an assumed radiative efficiency for the Eddington accretion-rate normalization. The simulations themselves use standard GRMHD in an axisymmetric ideal-MHD setup. The physical content beyond the equations is the magnetically arrested disk scenario and the unstated assumption that mechanical outflow power equals the observable X-ray luminosity.

free parameters (5)
  • Black hole mass MBH = 20 Msun
    Chosen in Section 4 to convert simulation outflow efficiency to physical power in CGS units. It is not constrained by the simulations or by the ULX sample, and it directly sets the luminosity scale.
  • Total mass accretion rate Mdot = 0.05 Mdot_Edd
    Assumed sub-Eddington rate in Section 4. Combined with MBH and the assumed radiative efficiency eta=0.1, it sets Pout = eta * Mdot * c^2 and places the result in the ULX band.
  • Radiative efficiency for Eddington accretion-rate normalization = 0.1
    Used in Section 4 to define Mdot_Edd. The simulated flow is advective and radiatively inefficient by construction, so applying a 10 percent radiative efficiency to normalize the accretion rate is inconsistent with the flow's assumed character.
  • Initial plasma beta values = 0.1 and 1
    Set by hand in Section 3 to control initial magnetic field strength. The results, including whether magnetically arrested disks form, depend strongly on these values.
  • Initial magnetic field normalization = max(pgas)/max(pmag) = 0.1
    Used in Section 3 to set the initial magnetic field strength in the validation run. This is an arbitrary normalization that affects the subsequent evolution.
assumptions (5)
  • domain assumption Ideal MHD approximation with infinite conductivity, u_mu F^mu nu = 0.
    Invoked in Section 2, equation (9). No resistivity, ambipolar diffusion, or other non-ideal effects are included.
  • domain assumption The accreting plasma is optically thin, advective, and non-radiative; no radiation feedback or cooling is included.
    Stated in the abstract and Section 5. This is essential for the hard-state identification but also prevents the simulation from predicting X-ray luminosity.
  • domain assumption 2D axisymmetry with N_phi=1 is sufficient to capture magnetically arrested disk and outflow behavior.
    Section 3 sets the grid to 256x256x1. Axisymmetric simulations cannot capture non-axisymmetric instabilities that can limit magnetic flux accumulation and modify jet power.
  • domain assumption A Fishbone-Moncrief torus is an appropriate initial condition for a sub-Keplerian advective accretion flow.
    Section 3 initializes the disk from a Fishbone-Moncrief torus. Real accretion flows may not start from such an equilibrium configuration.
  • standard math Standard GRMHD equations in Kerr-Schild coordinates describe the accretion flow.
    Section 2 uses conservation of particle number, energy-momentum, and Maxwell's equations. These are standard background equations for the code.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Magnetically arrested advective accretion flows and jets/outflows around stellar mass black holes: Explaining hard state ULXs with GRMHD simulations." pith.science (2026). https://pith.science/paper/RAGHEGYE

@misc{pith2026241118599,
  author       = {Pith},
  title        = {Pith review of: Magnetically arrested advective accretion flows and jets/outflows around stellar mass black holes: Explaining hard state ULXs with GRMHD simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RAGHEGYE}},
  note         = {Machine review of arXiv:2411.18599}
}
abstract

An optically thin advective accretion disk is crucial for explaining the hard state of black hole sources. Using general relativistic magnetohydrodynamic (GRMHD) simulations, we investigate how a large-scale, strong magnetic field influences accretion and outflows/jets, depending on the field geometry, magnetic field strength, and the spin parameter of the black hole. We simulate a sub-Eddington, advective disk-outflow system in the presence of a strong magnetic field, which likely remains in the hard state. The model simulations based on HARMPI successfully explain ultra-luminous X-ray sources (ULXs) in the hard state, typically observed with luminosities ranging from $10^{39}$ - $10^{40}$ ergs s$^{-1}$. Our simulations generally describe the bright, hard state of stellar-mass black hole sources without requiring a super-Eddington accretion rate. This work explores the characteristics of ULXs without invoking intermediate-mass black holes. The observed high luminosity is attributed to the energy stored in the strong magnetic fields, which can generate super-Eddington luminosity. The combined energy of the matter and magnetic field leads to such significant luminosity.

Figures

Figures reproduced from arXiv: 2411.18599 by the authors.

Figure 1
Figure 1. (a) The density contour at the flow vertical plane in log-scale at [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Density contour in log-scale at time t = 20, 000rg/c with the white lines indicating magnetic field lines for the case a = 0.998, initial plasma￾β = 0.1 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Density contour at the flow vertical plane in log-scale at time [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Density contour at the flow vertical plane in log-scale at time [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: Magnetic field in Gauss for (i) a = 0.998, initial plasma-β = 0.1, (ii) a = 0.5, initial plasma-β = 0.1, and (iii) a = 0.998, initial plasma-β = 1 [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: Outflow Power in erg/s: (a) Power for (i) [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [21]

    Chatterjee, R

    K. Chatterjee, R. Narayan, Astrophys. J., 941, 30 (2022)

    work page 2022

  2. [1]

    N. I. Shakura and R. A. Sunyaev, Astron. and Astrophys., 24, 337 (1973)

    work page 1973

  3. [2]

    Chakrabarti, L.G

    S. Chakrabarti, L.G. Titarchuk, Astrophys. J., 455, 623 (1995)

    work page 1995

  4. [3]

    Narayan, I

    R. Narayan, I. Yi, Astrophys. J., 452, 710 (1995)

    work page 1995

  5. [4]

    A. R. Das, B. Mukhopadhyay, Astrophys. J., 955, 86 (2023)

    work page 2023

  6. [5]

    R. P. Fender, T. M. Belloni, E. Gallo, MNRAS, 355, 1105 (2004)

    work page 2004

  7. [6]

    Mondal and B

    T. Mondal and B. Mukhopadhyay, MNRAS, 495, 350 (2020)

    work page 2020

  8. [7]

    S. R. Datta, T. Mondal, and B. Mukhopadhyay, MNRAS, 513, 204 (2022)

    work page 2022

Show all 27 references
  1. [8]

    Mondal, B

    T. Mondal, B. Mukhopadhyay, Mon. Not. R. Astron. Soc. Lett., 482, L24 (2019)

    2019

  2. [9]

    Mondal, B

    T. Mondal, B. Mukhopadhyay, MNRAS, 476, 2396 (2018)

    2018

  3. [10]

    Mondal, B

    T. Mondal, B. Mukhopadhyay, MNRAS, 486, 3465 (2019)

    2019

  4. [11]

    Mukhopadhyay, K

    B. Mukhopadhyay, K. Chatterjee, Astrophys. J., 807, 43 (2015)

    2015

  5. [12]

    Pathak, and B

    M. Pathak, and B. Mukhopadhyay, present volume (2024)

    2024

  6. [13]

    H. Feng, P. Kaaret, Astrophys. J., 696, 1712 (2009)

    2009

  7. [14]

    Soria, G

    R. Soria, G. Risaliti, M. Elvis, G. Fabbiano, S. Bianchi, Z. Kuncic, Astrophys. J., 695, 1614 (2009)

    2009

  8. [15]

    Soria, D

    R. Soria, D. S. Wong, MNRAS, 372, 1531 (2006)

    2006

  9. [16]

    H. Feng, F. Rao, P. Kaaret, Astrophys. J., 710, L137 (2010)

    2010

  10. [17]

    C. F. Gammie, J. C. McKinney, and G. T ´oth, Astrophys. J., 589, 444 (2003)

    2003

  11. [18]

    van Leer, J

    B. van Leer, J. Comput. Phys., 32, 101 (1979)

    1979

  12. [19]

    Harten, J

    A. Harten, J. Comput. Phys., 49, 357 (1983)

    1983

  13. [20]

    Fishbone, V

    L.G. Fishbone, V. Moncrief, Astrophys. J., 207, 962 (1976)

    1976

  14. [22]

    Porth, H

    O. Porth, H. Olivares, Y. Mizuno, Z. Younsi, L. Rezzolla, M. Moscibrodzka, H. Falcke, M. Kramer, CompAC, 4, 1 (2017)

    2017

  15. [23]

    R. Raha, B. Mukhopadhyay, K. Chatterjee, S.M. Gopika, Astronomy Reports, 67, S189 (2023)

    2023

  16. [24]

    G. S. Bisnovatyi-Kogan, A. A. Ruzmaikin, ApSS, 28, 45 (1974)

    1974

  17. [25]

    D.Blandford, R

    R. D.Blandford, R. L. Znajek, MNRAS, 179, 433 (1977)

    1977

  18. [26]

    D.Blandford, D

    R. D.Blandford, D. G. Payne, MNRAS, 199, 883 (1982)

    1982

  19. [27]

    Zdziarski, P

    A.A. Zdziarski, P. Pjanka, M. Sikora, L. Stawarz, MNRAS, 442, 3243 (2014)

    2014

Pith tools

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