Pith. sign in

REVIEW 3 major objections 5 minor 17 references

Accretion Disk-Outflow/Jet and Hard State ULXs

T0 review · 3 major / 5 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 Highly magnetized advective accretion flows around stellar-mass black holes can produce outflow power in the observed ULX luminosity range, according to GRMHD simulations.

desk verdict GRMHD check of MM19 is real and the ZAMO BZ/BP decomposition is useful, but the abstract overstates mechanical outflow power as ULX luminosity; the body admits it is an upper bound. read the letter →

arxiv 2411.18545 v1 pith:ORM5NE4C submitted 2024-11-27 astro-ph.HE

Mayank Pathak , Banibrata Mukhopadhyay This is my paper · ORCID
classification astro-ph.HE
keywords ultraluminousX-raysourceshardstateGRMHDsimulationsmagnetizedadvectiveaccretionflowsBlandford-ZnajekmechanismBlandford-Paynemagneticallyarresteddisksoutflowpower
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 argues that hard-state ultraluminous X-ray sources (ULXs) can be explained as highly magnetized, advective accretion flows around stellar-mass black holes, without invoking intermediate-mass black holes or modified Eddington limits. Using general relativistic magnetohydrodynamic (GRMHD) simulations of such flows, the authors show that the outflow power they produce falls inside the observed ULX luminosity range, and that the magnetic fields near the black hole reach about $10^7$ gauss. They also decompose the outflow power in the frame of a zero angular momentum observer and identify two components, one tied to black-hole spin (Blandford-Znajek) and one tied to magnetocentrifugal launching from the disk (Blandford-Payne). Because the simulations omit radiative cooling, the computed power is presented as a precursor or upper bound to the X-ray luminosity rather than the luminosity itself.

What carries the argument

The load-bearing object is the outflow power $P(r)=\dot{M}(r)-\dot{E}(r)$, computed from the stress-energy tensor of the GRMHD flow, where $\dot{M}$ is the mass accretion rate and $\dot{E}$ the inward energy flux. It is normalized to physical units by multiplying the dimensionless ratio by $\dot{M}_{\rm phy} c^2$, with $\dot{M}_{\rm phy}=0.05\,\dot{M}_{\rm Edd}$ and a 20-solar-mass black hole. A second object is the ZAMO-frame radial energy flux, which separates the power into a spin-dependent inner peak and an outer disk-launched peak. Together these definitions carry the argument: the first places the simulated flows in the ULX band, and the second attributes the two components to the Blandford-Znajek and Blandford-Payne mechanisms.

What would settle it

Measure the accretion rate of a hard-state ULX with an independently determined black-hole mass; if the rate is substantially below 0.05 Eddington for a ~20 solar-mass black hole, the normalized outflow power falls below the observed ULX band, contradicting the model.

Watch

Extended reading notes

Core claim

The central discovery claimed is that a 20-solar-mass black hole accreting at $0.05$ of the Eddington rate, surrounded by a highly magnetized advective disk, produces time-averaged outflow power in the ULX band. The power is defined as the difference between the mass accretion rate and the inward energy flux, $P(r)=\dot{M}(r)-\dot{E}(r)$, normalized by $\dot{M}_{\rm phy} c^2$. In both the standard (SANE) and magnetically arrested (MAD) initial magnetic configurations, the authors find powers in the ULX range, with MAD flows giving higher power; the required magnetic field at the black hole is about $10^7$ gauss, consistent with earlier steady-state calculations. In the ZAMO frame, the radial power profile shows a peak that grows with black-hole spin and ends near the ergosphere, interpreted as Blandford-Znajek spin extraction, and a second peak around $30$ gravitational radii, interpreted as Blandford-Payne magnetocentrifugal outflow. The authors state that the absence of radiative cooling means the simulated power is only the precursor to, or upper bound on, ULX luminosities.

Load-bearing premise

The load-bearing premise is that a hard-state ULX accretes at about 5% of the Eddington rate onto a 20-solar-mass black hole, so the mechanical outflow power, computed with radiative cooling omitted, can stand in for the X-ray luminosity.

Editorial extensions

If this is right

  • Hard-state ULXs can be modeled without intermediate-mass black holes or a modified Eddington limit.
  • Magnetic field strengths near $10^7$ gauss suffice for ULX-level outflow power, without invoking super-Eddington fields far from the black hole.
  • Magnetically arrested disks give higher outflow power than standard (SANE) disks, making MAD flows a promising engine for the most luminous hard-state ULXs.
  • ZAMO-frame power profiles separate spin-powered (Blandford-Znajek) from disk-launched (Blandford-Payne) outflows, providing a diagnostic for jet-launching mechanisms.
  • Because the simulations exclude radiative cooling, the computed power is an upper bound to the eventual X-ray luminosity, not the luminosity itself.

Reading between the lines

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

  • If the mechanical outflow power is converted to radiation with even moderate efficiency, these flows could also account for the kinetic power of ULX jets and winds, not only the X-ray band.
  • The spin dependence of the inner ZAMO peak predicts that higher-spin sources should show stronger BZ-dominated jet power; retrograde-spin simulations could cleanly isolate the BP component.
  • Because the power scales linearly with the adopted accretion rate and black-hole mass, the model makes a quantitative prediction: a ULX with an accretion rate an order of magnitude lower at the same mass would fall out of the ULX band.
  • The 2.5-dimensional axisymmetric setup likely exaggerates the magnetic-barrier oscillations in MAD runs; full 3D simulations may reduce the MAD-versus-SANE power gap.
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 / 5 minor

Summary. The paper uses the BHAC general-relativistic MHD code to simulate axisymmetric, magnetized advective accretion flows around a Kerr black hole with spin a=0.9375, starting from SANE and MAD initial conditions with initial plasma beta=100. It computes mass accretion rates, inward energy fluxes, and net outflow power profiles, normalizes these to physical power using Mdot_phy=0.05 Mdot_Edd and MBH=20 Msun, and reports that the resulting powers lie within the ULX range 3e39-3e41 erg/s. It further reports magnetic field strengths of order 1e7 G, compares the MAD run with an independent HARMPI simulation, and presents ZAMO-frame power profiles attributed to Blandford-Znajek and Blandford-Payne mechanisms. The central claim is that highly magnetized advective accretion flows can explain hard-state ULX luminosities without invoking intermediate-mass black holes or modified Eddington limits.

Significance. If fully established, the result would provide numerical support for the Mondal-Mukhopadhyay (MM19) steady-state model of hard-state ULXs and would strengthen the case that stellar-mass black holes with highly magnetized advective disks can power ULX-like outputs. The paper has several genuine strengths: it uses standard GRMHD methods and publicly available code, it performs a cross-code robustness check with HARMPI, it presents multiple time-averaging definitions of power, and the 1e7 G field strength is an emergent outcome rather than an input. The dimensionless power profiles are also a useful quantity for future comparisons. However, the step from simulated mechanical power to observed ULX X-ray luminosity is not established, and the physical normalization relies on hand-picked values; these issues are load-bearing for the paper's main claim.

major comments (3)
  1. [§3.2, Eq. (6), Fig. 4] The translation from code units to physical power uses hand-picked values Mdot_phy=0.05 Mdot_Edd and MBH=20 Msun. The reported power scales linearly with both quantities, so a lower accretion rate or a smaller black hole mass places the same dimensionless profiles below the ULX band (3e39 erg/s). The paper should present the dimensionless outflow efficiency P/(Mdot c^2) as the primary simulation output and either display the physical power as a function of the assumed Mdot_phy and MBH or clearly state that the ULX-range normalization is an assumption, not a prediction of the simulation.
  2. [§3.2 and §6 (also Abstract)] The simulations do not include radiative cooling, and Section 3.2 explicitly states that the computed outflow power is 'only the precursor to ULX luminosities or the upper bound.' Yet the Abstract claims the systems 'produce high luminosities like ULXs' and Section 6 claims 'high outflow power, well within the observed ULX luminosity range.' The simulated quantity is mechanical power (thermal, kinetic, and Poynting), which can be radiatively inefficient and need not emerge as X-ray luminosity. No specific hard-state ULX is compared with the model, and no radiative efficiency or emission mechanism is supplied. Please qualify the abstract and conclusion or add an explicit conversion from mechanical power to observable X-ray luminosity.
  3. [§5.1, Eq. (7)] The ZAMO-frame energy flux is not defined unambiguously. The contraction T^mu_nu u^nu e^mu has mismatched indices, and the stated radial vector e^mu=(0,1/grr,0,0) is not a unit vector in the Kerr metric unless grr is defined carefully; the unit radial vector should be e^r=1/sqrt(g_rr) or equivalently sqrt(g^rr). In addition, the jet region criterion '-T^mu_nu > 0' is not a well-defined scalar condition; the relevant component and observer frame should be specified. These issues affect the claimed decomposition into Blandford-Znajek and Blandford-Payne contributions.
minor comments (5)
  1. [Abstract, §1] Please fix typographical issues: 'around107 G' in the Abstract should be 'around 10^7 G', and 'Blanford-Znajek' in Section 1 should be 'Blandford-Znajek'.
  2. [§2] The setup text says 'MBM is the mass of the black hole'; this should read 'MBH'. Also, 'Fishbone Moncrief (FM) tours setup' should be 'torus setup'.
  3. [§3.2, Fig. 4] The caption for Fig. 4 states that 'all definitions' of time averaging are shown, but the text does not specify which line type corresponds to which of the three definitions in the list. Please identify each curve explicitly.
  4. [§3.1, Ref. [13]] The steady-flow radii req=10 rg and req=20 rg are attributed to an in-preparation paper [13]. Since Fig. 3 displays the red dashed lines that define these radii, the values are directly supported by the present paper; citing an unpublished work for them is unnecessary and makes verification harder.
  5. [§5.1] The sentence 'The jet region is considered to be the part of the simulation domain in which -T^mu_nu > 0 [17]' is unclear as written; please specify whether this is a particular component such as -T^r_t, and in which basis or observer frame the condition is applied.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular reduction found: the ULX-range power is an emergent dimensionless simulation ratio multiplied by a physically motivated accretion-rate normalization, and the cited self-results are not load-bearing.

full rationale

Walking the paper's derivation chain: the central object is the code-units net outflow power P(r) = Mdot(r) - Edot(r), with Mdot and Edot evaluated by surface integrals over the GRMHD stress-energy tensor (Eqs. 2-5). This dimensionless ratio is an emergent, time-averaged product of the simulation; it is not fit to ULX data. The dimensional conversion in Eq. (6) multiplies by Mdot_phy c^2 with Mdot_phy = 0.05 Mdot_Edd and MBH = 20 Msun. Those values are physically motivated choices for a stellar-mass advective flow, not constants fitted to the ULX luminosity band, and the simulated profiles are not constructed to land in the shaded region. The 10^7 G field in Sec. 3.3 is likewise an output from the evolved magnetic field, not an input. Self-citations are present but not load-bearing: [13] (in preparation) supplies the req values that are already displayed in Fig. 3, and [7] (MM19) provides the scenario being independently tested rather than the simulation result. The explicit caveat in Sec. 3.2 that radiative cooling is absent, so the power is 'only the precursor to ULX luminosities or the upper bound', is a genuine limitation of the astrophysical inference from mechanical power to X-ray luminosity, but it is not a circularity: the simulation does not assume the ULX luminosity it claims to produce. The dependence of the final power on Mdot_phy and MBH is a sensitivity caveat, not a self-referential reduction. No equation's output is equivalent to its input by construction.

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

The central claim rests on the hand-chosen physical normalization (Mdot_phy, MBH) and on the untested mapping from outflow power to luminosity. The simulation itself introduces no new entities or fitted constants; the 10^7 G field is an emergent quantity, though its relevance depends on the chosen parameters.

free parameters (5)
  • Physical accretion rate scale Mdot_phy = 0.05 Mdot_Edd
    Hand-chosen to represent an advective flow; the resulting physical power scales linearly with this number, so the placement of P(r) in the ULX band depends directly on it (Section 3, after eq. 5).
  • Black hole mass MBH = 20 Msun
    Assumed stellar-mass black hole; the Eddington normalization and Mdot_phy scale with MBH (Section 3).
  • Black hole spin a (main runs) = 0.9375
    Used for the main SANE/MAD power profiles; the BZ component in the ZAMO analysis depends on spin (Sections 2 and 5).
  • Initial plasma beta = 100
    Chosen initial magnetization for both SANE and MAD setups; no parameter survey over beta is presented (Section 2.1).
  • Steady flow radius req = 10 rg (SANE), 20 rg (MAD)
    Read off Fig. 3 and used to restrict the power computation to the 'inflow-outflow equilibrium' region; the values are cited to an in-preparation paper [13].
assumptions (4)
  • domain assumption The Kerr metric and ideal GRMHD equations accurately describe the accretion flow
    The BHAC code solves these equations in MKS coordinates; no explicit validation against 3D or radiative runs is given in this paper.
  • domain assumption Axisymmetry (2.5D) captures the relevant outflow physics
    The azimuthal dimension has only one cell and the authors note 2D effects exaggerate MAD barrier formation, so quantitative power may differ in 3D (Section 2).
  • ad hoc to paper The jet region is identified by -T^t_r > 0
    Used to compute the ZAMO-frame outflow power; follows Penna et al. (2013) [17] but is a modelling choice for this analysis (Section 5.1).
  • ad hoc to paper Mechanical outflow power is an upper bound to the observable luminosity
    No radiative cooling is included, so the paper treats the computed power as a 'precursor to ULX luminosities or the upper bound' (Section 3.2). This assumption bridges the simulation output and the abstract's luminosity claim.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Accretion Disk-Outflow/Jet and Hard State ULXs." pith.science (2026). https://pith.science/paper/ORM5NE4C

@misc{pith2026241118545,
  author       = {Pith},
  title        = {Pith review of: Accretion Disk-Outflow/Jet and Hard State ULXs},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ORM5NE4C}},
  note         = {Machine review of arXiv:2411.18545}
}
abstract

Ultraluminous X-ray sources (ULXs) have been objects of great interest for the past few decades due to their unusually high luminosities and spectral properties. A few of these sources exhibit super-Eddington luminosities assuming them to be centering around stellar mass objects, even in their hard state. It has been shown via numerical steady state calculations that ULXs in hard state can be interpreted as highly magnetised advective accretion sources around stellar mass black holes. We use general relativistic magnetohydrodynamic (GRMHD) framework to simulate highly magnetised advective accretion flows around a black hole and show that such systems can indeed produce high luminosities like ULXs. We also verify that the magnetic fields required for such high emissions is around $10^7$ G, in accordance with previous numerical steady state calculations. We further present power profiles for zero angular momentum observer (ZAMO) frame. These profiles show interesting features which can be interpreted as effects of emission due to the Blandford-Znajek and Blandford-Payne mechanisms.

Figures

Figures reproduced from arXiv: 2411.18545 by the authors.

Figure 1
Figure 1. Density contours for SANE with magnetic field streamlines. The top [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Density contours for MAD with magnetic field streamlines. The top [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Time averaged accretion rate profiles. Red dashed line indicates [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Power profiles for the all definitions of time average. Shaded region [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]
Figure 5
Figure 5. Figure 5: (a) Density contour overplotted with magnetic field streamlines at [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]
Figure 6
Figure 6. Figure 6: Magnetic field profiles for SANE and MAD simulations. [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 7
Figure 7. Figure 7: Magnetic field components for SANE and MAD simulations. [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]
Figure 8
Figure 8. Figure 8: Power profiles for ZAMO observer for different [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

  1. [1]

    Ebisawa, P

    K. Ebisawa, P. Życki, A. Kubota, T. Mizuno, K.-y. Watarai, ApJ, 597, 780, 2003

    work page 2003

  2. [2]

    M. C. Begelman, ApJ, 568, L97, 2002

    work page 2002

  3. [3]

    Herold, Phys

    H. Herold, Phys. Rev. D, 19, 2868, 1979

    work page 1979

  4. [4]

    Bachetti et al., Nature, 514, 202, 2014

    M. Bachetti et al., Nature, 514, 202, 2014

    work page 2014

  5. [5]

    Fürst et al., ApJ, 831, L14, 2016

    F. Fürst et al., ApJ, 831, L14, 2016

    work page 2016

  6. [6]

    G. L. Israel et al., Science, 355, 817, 2017

    work page 2017

  7. [7]

    Mondal, B

    T. Mondal, B. Mukhopadhyay, MNRAS, 482, L24, 2019

    work page 2019

  8. [8]

    Porth, H

    O. Porth, H. Olivares, Y. Mizuno et al., Comput. Astrophys. 4, 1, 2017

    work page 2017

Show all 17 references
  1. [9]

    L. G. Fishbone, V. Moncrief, ApJ, 207, 962, 1976

    1976

  2. [10]

    Chatterjee, R

    K. Chatterjee, R. Narayan, ApJ, 941, 30, 2022

    2022

  3. [11]

    R. Raha, K. Chatterjee, B. Mukhopadhyay, present volume, 2024

    2024

  4. [12]

    Narayan, A

    R. Narayan, A. Chael, K. Chatterjee et al., MNRAS, 511, 3795, 2022

    2022

  5. [13]

    Pathak, B

    M. Pathak, B. Mukhopadhyay, et al., in preparation

  6. [14]

    Tchekhovskoy, R

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

    2011

  7. [15]

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

    1977

  8. [16]

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

    1982

  9. [17]

    Penna, R

    R. Penna, R. Narayan, A, Sądowski, 436, 3741, 2013

    2013

Pith tools

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