{"id":"9d6bbc20-b524-48d3-95f0-815d947aae08","arxiv_id":"2603.17922","paper_version":2,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"GRRMHD puffy-disc simulations of sub-Eddington accretion differ from Shakura–Sunyaev-type models in photosphere thickness, inner-edge location, surface density, and radially varying effective viscosity.","lead":"Radiative GRMHD simulations of mildly sub-Eddington black-hole accretion find a magnetically stabilized “puffy” disc with a thick photosphere around a denser core. Compared with classic analytic models, the inner edge sits closer in, surface density is lower, and effective viscosity rises steeply near the hole.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Reported offsets from analytic models may be numerical artifacts of the GRRMHD setup rather than robust physical differences.","rationale":"The reader already isolated the identical load-bearing premise—that the numerical setup faithfully captures the physical vertical structure and effective viscosity—and correctly flagged the abstract-only status as precluding verification. No stronger internal inconsistency appears in the available text; the concern is therefore confirmatory rather than novel, leaving the UNVERDICTED / LOW-confidence assessment intact.","tokens_in":2094,"tokens_out":406,"duration_ms":9089,"concrete_test":"Re-run the fiducial model at double vertical resolution (and, separately, with a different initial poloidal-field geometry) and recompute the four diagnostic quantities (photospheric H/R, ISCO-relative inner edge, Σ(r), and α_eff(r)). If any quantity shifts by more than ~20 % the offsets are resolution- or setup-dependent and the comparison to analytic models loses force.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on treating the GRRMHD puffy-disc runs as ground truth for photosphere height, truncation radius, surface density, and radial viscosity profile. That premise is insecure because the abstract supplies no resolution study, 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 the scale height, lower the column, and produce a steeply rising effective α near the ISCO. Without those controls the claimed discrepancies with Shakura–Sunyaev-type models cannot be distinguished from numerical bias.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","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.","tokens_in":2270,"tokens_out":850,"duration_ms":12715,"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":[{"comment":"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.","section":"Abstract (comparison framework)"},{"comment":"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.","section":"Abstract (listed differences)"}],"minor_comments":[{"comment":"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.","section":"Abstract"},{"comment":"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.","section":"Abstract"}],"recommendation":"uncertain","confidential_remarks":"Only the abstract was supplied for this review; the full manuscript (figures, resolution studies, tables, methods) was unavailable. A definitive recommendation is therefore impossible. Once the complete paper is in hand the principal risk to check is whether the claimed offsets survive standard numerical-convergence and magnetic-topology tests; if they do, the work is likely suitable for minor or major revision rather than rejection. Scope appears appropriate for an astro-ph.HE journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one thing to know is that this is an abstract-only comparison paper: GRRMHD puffy discs (optically thick, mildly sub-Eddington) versus Shakura–Sunyaev-type models. The claimed offsets are concrete—geometrically thick photosphere, inner edge closer in than analytics, lower surface density, and a radially rising effective α—and they would matter for continuum fitting and warm-corona pictures if they hold. We cannot verify them from what we have.\n\nWhat looks solid on the face of it is the framing. They are not re-deriving puffy discs from scratch; they are using existing simulation machinery associated with this group to run a systematic structural comparison (inner structure, flow properties, viscosity, truncation). Magnetic stabilisation of a denser cooler core under a puffed optically thick layer is a useful organising picture, and the differences they list are falsifiable rather than hand-wavy. Circularity looks low: the numerics are being judged against external analytic models, not fitted to themselves.\n\nThe soft spot is real and load-bearing. The stress-test is right: without resolution studies, magnetic topology scans, radiative-transfer convergence, or boundary checks, we cannot separate physical offsets from numerical bias in a regime where vertical structure is a delicate balance of magnetic support, cooling, and turbulent stress. Under-resolved MRI or photospheric layers can easily inflate scale height, drop column, and produce a steeply rising α near the ISCO. That is not a minor quibble; it is the premise of treating the runs as ground truth. Novelty is incremental—puffy GRRMHD and thin-disc comparisons already exist in this literature—so the value lives entirely in how carefully those diagnostics are controlled.\n\nWho it is for: people who build or use continuum-fitting and warm-corona models for soft-state XRBs, and anyone running or interpreting GRRMHD of sub-Eddington flows. It deserves a serious referee if the full paper ships the resolution and sensitivity tests; otherwise it is a methods note with suggestive plots. I would not cite from the abstract alone. Send it to peer review rather than desk-reject—the claim is sharp enough to be useful once the numerics are stress-tested—but expect the referee to demand those controls before accepting the offsets as physical.","headline":"Abstract-only GRRMHD comparison of puffy discs to thin-disc analytics: interesting offsets, but we cannot yet tell physics from numerics.","tokens_in":2950,"tokens_out":565,"would_cite":false,"duration_ms":5830,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"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.","keywords":["accretion discs","puffy discs","GRRMHD simulations","sub-Eddington accretion","black-hole X-ray binaries","Shakura-Sunyaev model","effective viscosity","radiative magnetohydrodynamics"],"falsifier":"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.","tokens_in":3016,"feed_emoji":"🌀","tokens_out":891,"duration_ms":21519,"temperature":0.7,"pith_summary":"The paper compares general-relativistic radiative magnetohydrodynamic simulations of optically thick, mildly sub-Eddington accretion (puffy discs) with classic analytic models of the Shakura–Sunyaev type used for the soft spectral state of X-ray binaries. Magnetic fields stabilize the flow, producing a stratified structure with a denser cooler core surrounded by a puffed-up optically thick region that resembles a warm corona. The simulations recover some qualitative similarities to analytic models, yet they systematically differ in four concrete respects: the photosphere is geometrically thick, the inner edge lies closer to the black hole than the models assume, the surface density is substantially lower than predicted, and the effective viscosity parameter is not constant but rises steeply in the innermost region. Because analytic models remain the standard framework for interpreting soft-state spectra, these offsets would change how observers map continuum data to accretion rate, disc size, and black-hole spin.","feed_headline":"Puffy black-hole discs break four thin-disc assumptions","feed_subtitle":"Photosphere is thick, inner edge closer in, density lower, viscosity rises steeply","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Puffy discs defy thin-disc models with thick photospheres","Simulated puffy discs put inner edge closer than analytics","Magnetic fields stabilize sub-Eddington puffy accretion discs","Puffy discs show lower density and steeply rising viscosity","Vertically stratified puffy discs break constant-viscosity rule"],"cache_read_input_tokens":128,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Puffy discs defy thin-disc models with thick photospheres","Simulated puffy discs put inner edge closer than analytics","Magnetic fields stabilize sub-Eddington puffy accretion discs","Puffy discs show lower density and steeply rising viscosity","Vertically stratified puffy discs break constant-viscosity rule"]},"model":"grok-4.5","effort":"low","cost_usd":0.004304,"raw_usage":{"total_tokens":1296,"prompt_tokens":820,"num_sources_used":0,"completion_tokens":85,"cost_in_usd_ticks":43040000,"prompt_tokens_details":{"text_tokens":820,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":391,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":820,"tokens_out":85,"duration_ms":3674,"temperature":1.0,"reasoning_tokens":391,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T22:51:01.783350+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"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.","supporting_citations":[],"review_version":1}