REVIEW 2 major objections 2 minor 1 cited by
Radiative GRMHD simulations of puffy accretion discs: Numerical versus analytical models of sub-Eddington accretion
T0 review · 2 major / 2 minor · reviewed 2026-07-13 · grok-4.5
Pith's one-line read Simulated puffy accretion discs around stellar-mass black holes share some features with analytic thin-disc models but differ in photosphere thickness, inner-edge location, surface density, and viscosity.
desk verdict Abstract-only GRRMHD comparison of puffy discs to thin-disc analytics: interesting offsets, but we cannot yet tell physics from numerics. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The puffy-disc GRRMHD simulations themselves, which evolve the coupled equations of general-relativistic magnetohydrodynamics and radiation transport to produce self-consistent vertical structure, surface density, and effective viscosity that can be compared directly with analytic Shakura–Sunyaev-type solutions.
What would settle it
A controlled suite of higher-resolution GRRMHD runs with varied initial magnetic geometries and radiative closures that recovers the analytic surface-density profile, a constant alpha viscosity, and the standard ISCO-linked inner edge would falsify the claim that the offsets are physical.
Extended reading notes
Core claim
Optically thick, mildly sub-Eddington GRRMHD accretion flows (puffy discs) are magnetically stabilized and vertically stratified, with a geometrically thick photosphere, an inner edge closer to the black hole than analytic models place it, surface densities far below analytic predictions, and an effective viscosity that rises steeply inward rather than remaining constant.
Load-bearing premise
The numerical setup (resolution, initial magnetic field, radiative transfer treatment, and boundary conditions) faithfully reproduces the physical vertical structure and effective viscosity of real sub-Eddington discs, so the reported offsets from analytic models are physical rather than numerical artifacts.
Editorial extensions
If this is right
- Surface densities used in continuum-fitting models of soft-state X-ray binaries are systematically overestimated relative to the simulated discs.
- The inner edge of a mildly sub-Eddington disc sits closer to the black hole than the analytic truncation radius assumed in most spectral models.
- The effective viscosity parameter cannot be treated as radially constant; it rises steeply in the innermost region.
- The photosphere is geometrically thick, so the observed spectrum is shaped by a warm, optically thick layer rather than a razor-thin surface.
Reading between the lines
- Continuum-fitting spin measurements that assume a thin-disc photosphere and constant alpha may carry a systematic bias once the thicker, lower-density photosphere is taken into account.
- The same offsets should appear in sub-Eddington simulations around supermassive black holes if the underlying microphysics is scale-invariant.
- Time-dependent spectral models that allow a radially rising alpha and a puffed photosphere would be a direct observational test of the simulation results.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript compares general-relativistic radiative magnetohydrodynamic (GRRMHD) simulations of optically thick, mildly sub-Eddington accretion onto a stellar-mass black hole (termed “puffy discs”) with established analytic and semi-analytic thin-disc models in the same regime (Shakura–Sunyaev and related). From the simulations the authors report that magnetic fields stabilize the flow, producing a puffed-up, optically thick region resembling a warm corona around a denser, cooler core. They then compare inner-disc structure, flow properties, effective viscosity, and inner-edge location with standard models, concluding that the simulated discs share some similarities but differ in four principal respects: a geometrically thick photosphere, an inner edge closer to the black hole than analytic models assume, significantly lower surface density, and a radially rising (non-constant) effective viscosity parameter in the innermost region.
Significance. If the reported offsets are robust, the work would supply a concrete, simulation-based calibration of the thin-disc framework that is routinely used to interpret soft-state X-ray binary spectra. Explicit quantification of photosphere height, truncation radius, surface density, and the radial run of α_eff would be directly useful for spectral modelling and for assessing the limitations of purely analytic prescriptions that omit magnetic support and radiative transfer. The use of full GRRMHD in the optically thick, mildly sub-Eddington regime is itself a non-trivial technical contribution.
major comments (2)
- [Abstract (comparison framework)] The central claim treats the GRRMHD puffy-disc runs as ground truth against which Shakura–Sunyaev-type models are judged. That premise is load-bearing, yet the abstract supplies no resolution study, MRI-quality factor, magnetic-field topology scan, radiative-transfer convergence test, or boundary-condition sensitivity check. In the mildly sub-Eddington, optically thick regime the vertical structure is set by a delicate balance among magnetic support, radiative cooling and turbulent stress; modest under-resolution of the MRI wavelength or of the photospheric layers can artificially inflate scale height, lower column density and produce a steeply rising effective α near the ISCO. Without those controls the four listed discrepancies cannot be distinguished from numerical bias.
- [Abstract (listed differences)] The abstract asserts that the photosphere is “geometrically thick,” the inner edge lies “closer to the central black hole than the analytic models assume,” surface density is “significantly lower,” and α_eff “rises steeply in the innermost region.” These quantitative statements are the paper’s main scientific payload, but no numerical values, radial profiles, error bars or comparison tables are given in the available text. A referee cannot assess whether the offsets are statistically significant or merely qualitative impressions until the corresponding figures and tables are examined.
minor comments (2)
- [Abstract] The abstract is clearly written and the four claimed differences are enumerated cleanly; once the full manuscript is available, the same clarity should be preserved in the figure captions and table headings that quantify those differences.
- [Abstract] Terminology “puffy disc” and “warm corona” should be defined operationally (e.g., by optical-depth or density contours) at first use so that later quantitative comparisons remain unambiguous.
Circularity Check
No significant circularity: independent GRRMHD outputs compared to external analytic models
full rationale
This is an abstract-only review of a comparison paper. The claimed results are differences between GRRMHD simulation outputs (photosphere thickness, inner-edge location, surface density, radial viscosity profile) and established external analytic/semi-analytic models (Shakura–Sunyaev and related). Those differences are not forced by construction: the simulations are not fitted to reproduce the analytic quantities being tested, nor is any uniqueness theorem or ansatz imported from the authors’ prior work to define the comparison. The term “puffy disc” is a prior construct of overlapping authors, but that is nomenclature, not a load-bearing self-citation that makes the reported offsets tautological. With only the abstract available there is no equation chain to reduce; the abstract itself presents an empirical numerical-vs-analytic comparison that is self-contained against external benchmarks. Score 0 is therefore the correct, proportionate finding.
Assumptions & free parameters
free parameters (3)
- accretion rate (mildly sub-Eddington)
- initial magnetic field strength/topology
- effective viscosity α (analytic models)
assumptions (3)
- domain assumption General relativistic magnetohydrodynamics with radiative transfer adequately describes optically thick sub-Eddington accretion.
- domain assumption Classic analytic/semi-analytic thin-disc solutions (Shakura–Sunyaev and related) are the correct reference models for the soft spectral state.
- standard math Standard GR and ideal MHD continuum equations hold in the simulated domain.
Cite this review
Pith. "Pith review of Radiative GRMHD simulations of puffy accretion discs: Numerical versus analytical models of sub-Eddington accretion." pith.science (2026). https://pith.science/paper/VQKANI6P
@misc{pith2026260317922,
author = {Pith},
title = {Pith review of: Radiative GRMHD simulations of puffy accretion discs: Numerical versus analytical models of sub-Eddington accretion},
year = {2026},
howpublished = {\url{https://pith.science/paper/VQKANI6P}},
note = {Machine review of arXiv:2603.17922}
}
read the original abstract
A widely accepted picture of an accretion flow in the luminous soft spectral state of X-ray binary systems is a geometrically thin disc structure much like the classic analytic solution of Shakura& Sunyaev. Although the analytic models are troubled by instabilities and miss important aspects of physics, such as magnetic fields, they are successfully used as a framework for interpreting observational data. Here, we compare the results of general relativistic radiative magnetohydrodynamic (GRRMHD) simulations of optically thick, mildly sub-Eddington accretion on a stellar-mass black hole (the puffy disc) with established analytic and semi-analytic accretion models in the same regime. From the simulations, we find that the accretion flow is stabilised by the magnetic field, with a puffed-up, optically thick region resembling a warm corona surrounding a denser and cooler disc core. However, the stratified vertical structure of the disc significantly influences the observational picture of such a system. We analyse the inner disc structure, flow properties, effective viscosity, and inner edge position, and compare them to the predictions of standard models. We find that the simulated discs share some similarities with the models; however, they differ in several important aspects, most notably: the photosphere is geometrically thick, the inner edge is located closer to the central black hole than the analytic models assume, the surface density is significantly lower than analytically predicted, and the effective viscosity parameter is not constant but rises steeply in the innermost region.
Forward citations
Cited by 1 Pith paper
-
The plunging region of thin accretion discs across the black hole spin range
MHD stresses in the plunging region of thin accretion disks increase with prograde black hole spin, and a new flux-freezing model describes the magnetic fields there across the spin range.
Reviewed July 13, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.