{"id":"3a15629b-0a42-42f4-9781-7f7e2995d855","arxiv_id":"1908.08396","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Black hole accretion at 0.6 times the Eddington rate can form a stable puffy disk, with a thin dense core, a thick Keplerian photosphere, and most inflow above the core.","lead":"This paper uses supercomputer simulations to describe a new type of disk of gas spiraling into a black hole, called a puffy disk. It has a thin dense core but a thick, magnetically supported, hot upper layer that changes how such disks shine and look from outside.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The puffy disk's beta~1 stability is inherited from the imposed poloidal-field reservoir; no run tests whether a disk at 0.6 Mdot reaches this state from generic initial data. Thus the new-class claim rests on untested initial-data dependence.","rationale":"The numerical execution appears competent: the authors report inflow equilibrium out to about 20 M, MRI quality factors Qtheta~25 and Qphi~20, and a 15000 GM/c^3 run; the figures are consistent with a stable, turbulent, radiatively supported disk. My concern is not resolution or code validity but physical genericity. The Letter's own Section 2 and Appendix reveal that the beta~1 state is not an emergent outcome: the initial data are constructed to advect poloidal magnetic flux so that magnetic pressure stabilizes the radiation-pressure-dominated disk, and each lower-Mdot run is seeded by scaling down a previous strongly magnetized solution. Thus the simulation proves that a specially prepared initial state evolves to a stable puffy disk, which is a valuable existence result, but the abstract and discussion claim more: a new class of solutions that should replace the slim-disk description at sub-Eddington rates. That claim requires the puffy branch to be an attractor for ordinary accretion flows. No run starts from weak-field, standard conditions, no magnetic-flux parameter study is presented, and the azimuthal wedge plus M1 closure, while standard, add further unquantified uncertainties. The reader's conditionality is therefore appropriate; the decisive missing test is an initial-condition robustness study. If that test confirms convergence to the same puffy state, the concern is resolved; if not, the central claim should be weakened to existence of an initial-data-dependent magnetized disk branch. I mark partial rather than full agreement with the reader because I center the initial-data dependence; the M1-closure issue, while real, is secondary to the stability mechanism in my reading.","tokens_in":10157,"tokens_out":6678,"duration_ms":75735,"concrete_test":"Run the koral setup from this Letter at Mdot = 0.6 Mdot_Edd starting from a standard weakly magnetized Fishbone-Moncrief torus (no quadrupole magnetic reservoir, no rescaled strongly magnetized state), evolving for at least 15000 GM/c^3 with the same grid and coordinates; measure the time-averaged plasma beta and photospheric half-thickness H/r at r~10-20 M. If beta~1 and H~r do not develop spontaneously, the reported puffy solution is an artifact of the chosen initial magnetic-flux advection rather than the generic sub-Eddington state claimed. A complementary control would rerun the published initial conditions with the radial-field advection removed to isolate its role.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper itself flags the key weakness. Section 2 says: \"To ensure thermal stability (Zheng et al. 2011) the disk was made to advect poloidal magnetic field with a significant radial component,\" and the Appendix says only \"reservoirs of mass with quadrupole topology of the magnetic field\" are used, with successively lower-Mdot runs seeded by rescaling a previous \"strongly magnetized stable solution.\" Consequently the beta~1 state that is invoked as the mechanism of thermal stability and as the cause of magnetic puffy support is placed into the simulation via initial data. The central claim that puffy disks are a new class of generic sub-Eddington solutions, offering \"a more realistic description of black hole disks than the slim disk model,\" requires that this state be reached from generic accretion conditions, not merely that it persist when initialized in a specially prepared, magnetically advecting reservoir. The single 0.6 Mdot, Schwarzschild, M1-closure run demonstrates existence of a family of initial data that relaxes to a stable puffy disk; it does not demonstrate that sub-Eddington disks generically occupy this branch. The M1 closure also affects the radiation field, but the initial-field dependence is the more load-bearing gap because it targets the stability mechanism itself.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter reports global, three-dimensional, radiative GRMHD simulations with the koral code of a non-spinning 10 solar-mass black hole at a mass accretion rate of 0.6 Mdot_Edd. The resulting disk has a high-density core with density scale-height h_rho ~ 0.1r, but its photosphere lies at H ~ r, so the disk is simultaneously 'thin' by the density measure and 'thick' by the photospheric measure. The flow remains nearly Keplerian up to the photosphere, the plasma parameter beta = (pgas + prad)/pmag is order unity, and the disk is described as thermally stable despite radiation-pressure dominance. The authors argue that magnetic pressure support and inward advection of radiation define a new class of 'puffy' accretion disks distinct from thin, slim, and thick models, and that the resulting luminosity of 0.36 L_Edd is below the thin-disk expectation because a significant fraction of radiation is captured by the black hole. The Appendix gives numerical details: resolution 320x320x32, a pi/2 azimuthal wedge, MRI quality factors Q_theta ~ 25 and Q_phi ~ 20, inflow equilibrium out to r ~ 20 M over 15000 GM/c^3, and the seeding of lower-Mdot runs by rescaling a previous strongly magnetized stable solution.","tokens_in":10375,"tokens_out":6228,"duration_ms":66355,"significance":"If the claim is robust, this is a significant result: it would establish a sub-Eddington, optically thick disk branch in which the photosphere is decoupled from the density scale-height, magnetic pressure stabilizes a radiation-pressure-dominated state, and radiation advection reduces the effective luminosity. This would challenge the usual identification of sub-Eddington accretion disks with thin disks and has concrete observational consequences for spectra, inclination-dependent images, and timing behavior. The paper's strengths include a genuine numerical experiment with no fitting to data, explicit MRI resolution quality factors, inflow equilibrium, and ray-traced images that are falsifiable predictions. The main caveat is that the claimed new class is not yet shown to be the generic outcome of sub-Eddington accretion rather than a state reached only from specially prepared, strongly magnetized initial data.","major_comments":[{"comment":"The stability mechanism is imposed through the initial conditions rather than shown to emerge. Section 2 states that “to ensure thermal stability … the disk was made to advect poloidal magnetic field with a significant radial component,” and the Appendix states that only reservoirs with quadrupole magnetic-field topology are used, with successively lower-Mdot runs seeded by rescaling a previous “strongly magnetized stable solution,” preserving beta. Because beta ~ 1 is the mechanism invoked for both magnetic support and thermal stability, the paper demonstrates the existence of a family of initial data that relaxes to a puffy disk, but it does not demonstrate that the puffy branch is the generic sub-Eddington solution. Please add at least one 0.6 Mdot_Edd run starting from different initial data, such as a much weaker poloidal field, a purely toroidal field, or beta >> 1 in the initial torus, and show whether it converges to the same puffy state or behaves differently. At minimum, reword the abstract and introduction so the claim is explicitly about a stable state reached from strongly magnetized initial conditions rather than about the generic sub-Eddington solution.","section":"§2 (initial conditions) and Appendix"},{"comment":"Thermal stability is asserted from a single run with no perturbation or convergence test. A single trajectory that reaches inflow equilibrium over 15000 GM/c^3 does not exclude a slow growth of the thermal instability or a long-lived metastable branch. Please provide a quantitative stability diagnostic: time histories of midplane beta, h_rho, H, and luminosity over the full run; a linear-growth or fluctuation analysis of the relevant thermal modes; or a second run at higher resolution or with full 2pi azimuth to show that the equilibrium and its fluctuations are converged. Without such evidence, the word “stable” in the title is stronger than what the simulation alone establishes.","section":"§2, Figs. 3–4"},{"comment":"The radiation field is evolved with the M1 closure (Levermore 1984; Sadowski et al. 2013), which is known to systematically affect the angular distribution of radiation and the diffusion–free-streaming transition. The claimed sub-Eddington luminosity of 0.36 L_Edd and the detailed funnel radiation pattern in Fig. 3 (lower right) could be sensitive to this closure. While M1 is a standard approximation for global radiative GRMHD, the paper should state the expected quantitative uncertainty from M1 and, ideally, test at least one time-averaged snapshot with a different closure or a Monte-Carlo post-processing step to confirm that the inward radiation advection at r < 10 M and the photon capture fraction are not artifacts of the closure.","section":"§2, radiation transport"}],"minor_comments":[{"comment":"“appearence” should be “appearance”.","section":"§1"},{"comment":"The detailed numerical setup is deferred to “Lančová et al. (in preparation)”; for a Letter, at least the initial density and magnetic-field profiles, the radiative boundary conditions, and the exact rescaling procedure used to produce lower-Mdot runs should be stated explicitly.","section":"Appendix"},{"comment":"The lower-right panel caption lists radiation temperature contours “from top to bottom” of 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 MK, but the panel itself shows contours at 6×10^6, 8×10^6, and 10^7 K; please unify these notations.","section":"Fig. 3 caption"},{"comment":"The reference list contains two entries for “Jiang et al. 2019” (arXiv:1904.01674 and ApJ 880, 67); please check whether these are duplicate citations of the same work and, if so, consolidate them.","section":"References"},{"comment":"The use of both “R” and “r” for cylindrical radius is confusing; label the coordinate consistently and state clearly whether the super-Keplerian black contour refers to a time-averaged or instantaneous quantity.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is well executed but the central claim, as stated, overreaches the evidence: the stabilization mechanism is placed in the initial data, and stability is asserted from a single run. I believe the requested control runs and stability diagnostics are feasible within the same numerical framework and should be included. I see no grounds for rejection, provided the authors either supply those tests or carefully restrict the claims to the existence of a stable, strongly magnetized branch."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a solid first report of a new type of radiative GRMHD disk: at 0.6 M_Edd around a Schwarzschild hole, the simulated disk has a thin dense core (h_rho ~0.1r) but a photosphere at H~r, with most accretion happening above the core. That decoupling of density scale-height from photosphere is genuinely new, and the claim that magnetic pressure (beta~1) stabilizes a radiation-pressure-dominated disk addresses a long-standing puzzle. The simulation is well resolved, reaches inflow equilibrium to r~20M, and the ray-traced images are a nice complement.\n\nThe soft spot is the one you've identified: stability is partly built into the initial data. Section 2 says the disk 'was made to advect poloidal magnetic field' specifically to ensure thermal stability, and the appendix describes quadrupole-field mass reservoirs that produce 'strongly magnetized stable solutions.' So the paper demonstrates that a carefully prepared initial state relaxes to a stable puffy disk, not that generic sub-Eddington accretion flows end up there. The single accretion rate, single spin, M1 closure, and wedge domain are reasonable for a Letter but do limit the 'new class of solutions' claim. None of this kills the finding, but it should be framed as an existence proof with a strong hint of generic relevance, not a settled law.\n\nThe paper is honest about these issues, and the authors acknowledge earlier collapsed-disk simulations. The citation pattern is appropriate; they build on Sadowski 2016 and Zheng et al. 2011. My verdict: worth a serious referee. I'd want the initial-data dependence tested in the revision—varied field topologies, maybe an un-magnetized or weaker-field start, and a check that the puffy state is reached rather than imposed. Also, a convergence test would help.\n\nFor a reader in disk theory, this is a useful and stimulating paper. For the field at large, it's a candidate resolution to the thermal stability problem, not yet a confirmed one. I'd accept it for peer review and suggest the authors soften their 'more realistic description' language until the initial-data question is settled.","headline":"A genuinely new disk state in radiative GRMHD, but the stability claim leans on initial data engineered to be magnetically stable.","tokens_in":10996,"tokens_out":2475,"would_cite":true,"duration_ms":24452,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that sub-Eddington black hole accretion disks can settle into a thermally stable “puffy” state: a thin dense core, a thick magnetically supported, optically thick layer, and a photosphere at height comparable to the…","keywords":["black hole accretion disks","radiative GRMHD simulations","puffy disks","radiation pressure","magnetic pressure support","thermal stability","sub-Eddington accretion","M1 closure"],"falsifier":"Run the same radiative GRMHD simulation at 0.6 Eddington rate starting from a radiation-pressure-dominated thin disk with only weak, tangled magnetic fields and no imposed radial poloidal flux: the puffy layer should not form if the claimed state depends on the advected field. Observationally, measure the inclination dependence of the isotropic luminosity and color temperature of a bright sub-Eddington black hole binary; puffy disks predict a bright axial funnel, a dark spot over the hole, and a large spectral color correction, unlike thin or slim disk images.","tokens_in":9952,"feed_emoji":"🌀","tokens_out":9511,"duration_ms":87107,"temperature":0.7,"pith_summary":"This paper reports a new class of black hole accretion disk solutions found in three-dimensional, general-relativistic radiative magnetohydrodynamic simulations. For a non-spinning ten-solar-mass black hole accreting at 0.6 Eddington rate, the disk has a dense core only $h_\\rho\\sim 0.1r$ tall, so it looks thin by density, but its photosphere sits at $H\\sim r$, making it geometrically thick. Most of the inflow, turbulence, and radiation transport happen in this puffy, optically thick, magnetized layer, where magnetic pressure roughly balances gas plus radiation pressure, $\\beta\\sim1$. The authors argue that this configuration is thermally stable where standard radiation-pressure-dominated thin disks are unstable, and that it captures part of its own radiation, so it is dimmer than a thin disk at the same accretion rate.","feed_headline":"At 0.6 Eddington, black hole disks go puffy and stay stable","feed_subtitle":"The disk has a thin dense core but a thick, magnetically supported photosphere that swallows its own radiation.","key_machinery":"The key object is the puffy layer: the region between the density scale-height $h_\\rho$ and the photosphere $H$, which contains most of the mass inflow, turbulence, and advected radiation. The simulation creates this state by advecting poloidal magnetic field with a significant radial component from a quadrupole mass reservoir, so that after magnetorotational instability saturation the plasma parameter $\\beta=(p_{\\rm gas}+p_{\\rm rad})/p_{\\rm mag}\\sim1$. Magnetic pressure then stabilizes the disk against the thermal instability of radiation-pressure-dominated thin disks. Radiation transport is handled with the M1 closure, and the full three-dimensional flow is evolved in a Schwarzschild metric.","core_discovery":"The central claim is that at sub-Eddington rates around 0.6 Eddington, the equilibrium state of a radiative black hole accretion disk is not a canonical thin, slim, or thick disk but a hybrid “puffy” disk. It combines a high-density equatorial core of height $h_\\rho\\sim 0.1r$, so thin by the density scale-height measure, with an extended lower-density optically thick region reaching the photosphere at $H\\sim r$. The whole layer is turbulent and rotates at nearly Keplerian speed up to the photosphere; the accreting fluid is supported in part by magnetic pressure, and much of the radiation is advected inward and swallowed by the black hole rather than escaping vertically. The result is a disk that is thermally stable despite radiation-pressure dominance, with an inner luminosity of about 0.36 Eddington, less than a thin disk at the same accretion rate. The authors conclude that one-dimensional, height-integrated thin and slim disk models miss the meridional flow of matter and radiation that defines this state.","pith_inferences":["If the puffy state is the generic sub-Eddington equilibrium, then classifying disk states by observed Eddington ratio is unreliable: the same intrinsic accretion rate can appear subluminous or super-Eddington depending on viewing angle through the funnel.","The central role of advected poloidal flux suggests puffy disks may be tied to the magnetic flux threading the black hole; varying the initial flux or including black hole spin would test whether the state persists across the parameter space of real sources.","A concrete spectral prediction follows from the paper's structure: the low-density, hot puffy layer should produce a large color correction and a spectrum that is not a sum of thin-disk blackbodies, so joint spectral and inclination fitting of bright X-ray binaries can distinguish puffy from slim states.","The photon stagnation surface implies that radiation advection is not confined to the disk body; coronal models that treat disk and corona as separate thermal components may systematically misattribute the emission of this layer."],"forward_implications":["At accretion rates near 0.6 Eddington, black hole disks can be geometrically thin by density yet geometrically thick by photosphere, so the standard thin/slim/thick trichotomy is incomplete.","A radiation-pressure-dominated disk can be thermally stable when magnetic pressure contributes comparably to gas plus radiation pressure, resolving a long-standing instability that otherwise collapses such disks.","Radiation is advected inward and partly swallowed by the black hole, making the disk's luminosity lower than a thin disk at the same accretion rate.","The observable appearance depends strongly on inclination: a bright funnel near the axis, a dark shadow over the black hole, and obscuration of the near side at large inclinations.","Because most of the inflow occurs in the puffy layer above the dense core, one-dimensional height-integrated disk models miss the dominant accretion flow and radiation transport."],"supporting_citations":[{"why":"Defines the canonical thin α-disk with electron-scattering opacity and radiation-pressure-dominated inner region that the puffy disk is compared against and found dimmer than.","marker":"Shakura & Sunyaev 1973"},{"why":"Defines the slim disk model whose height-integrated advection picture the puffy disk extends beyond with a resolved three-dimensional meridional flow.","marker":"Abramowicz et al. 1988"},{"why":"Supplies the thick-disk funnel morphology and super-Keplerian rotation expectations used to interpret the puffy disk structure.","marker":"Jaroszyński, Abramowicz & Paczyński 1980"},{"why":"Identifies advection of poloidal magnetic field as the stabilizing mechanism that the initial conditions deliberately implement.","marker":"Zheng et al. 2011"},{"why":"Provides the GRMHD simulation setup and the earlier stable magnetized disk results that this work extends to lower accretion rates.","marker":"Ş adowski 2016"},{"why":"Documents the collapse of radiation-pressure-dominated MRI disks without a stabilizing field, the failure mode the puffy disk avoids.","marker":"Mishra et al. 2016"},{"why":"Confirms the thermal instability in shearing-box MRI simulations of radiation-dominated disks, motivating magnetic stabilization.","marker":"Jiang et al. 2019"},{"why":"Introduces the M1 closure for the radiation field used in the simulations.","marker":"Levermore 1984"},{"why":"Describes the GRMHD simulation code in which the simulations are run.","marker":"Ş adowski et al. 2013"},{"why":"Non-radiative simulations finding accretion predominantly above the high-density core, corroborating the puffy inflow pattern.","marker":"Zhu, & Stone 2018"}],"fun_headline_variants":["Thin core, puffy shell: black hole disk reveals new state","Magnetic puffs stabilize disk at 0.6 Eddington","Sub-Eddington disk is optically thick and stable","Black hole disk: thin core, thick puffy envelope","Disk captures its light, stays stable sub-Eddington"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The puffy, stable solution appears because the simulations are started with an advected poloidal magnetic field carrying significant radial flux; if real sub-Eddington disks do not naturally carry that magnetic flux, the state may be an artifact of the initial setup rather than the generic solution.","fun_headline_variants_meta":{"raw":{"variants":["Thin core, puffy shell: black hole disk reveals new state","Magnetic puffs stabilize disk at 0.6 Eddington","Sub-Eddington disk is optically thick and stable","Black hole disk: thin core, thick puffy envelope","Disk captures its light, stays stable sub-Eddington"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000697,"raw_usage":{"total_tokens":3188,"prompt_tokens":1022,"completion_tokens":2166,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":2080}},"tokens_in":638,"tokens_out":2166,"duration_ms":14787,"temperature":1.0,"reasoning_tokens":2080,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:42:36.243682+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same radiative GRMHD simulation at 0.6 Eddington rate starting from a radiation-pressure-dominated thin disk with only weak, tangled magnetic fields and no imposed radial poloidal flux: the puffy layer should not form if the claimed state depends on the advected field. Observationally, measure the inclination dependence of the isotropic luminosity and color temperature of a bright sub-Eddington black hole binary; puffy disks predict a bright axial funnel, a dark spot over the hole, and a large spectral color correction, unlike thin or slim disk images.","supporting_citations":[{"cited_title":"& Sunyaev\\,R.A","cited_arxiv_id":null,"evidence_quote":"Defines the canonical thin α-disk with electron-scattering opacity and radiation-pressure-dominated inner region that the puffy disk is compared against and found dimmer than."},{"cited_title":"1988, ApJ, 332, 646","cited_arxiv_id":null,"evidence_quote":"Defines the slim disk model whose height-integrated advection picture the puffy disk extends beyond with a resolved three-dimensional meridional flow."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies advection of poloidal magnetic field as the stabilizing mechanism that the initial conditions deliberately implement."},{"cited_title":"C., Johnson, L","cited_arxiv_id":null,"evidence_quote":"Documents the collapse of radiation-pressure-dominated MRI disks without a stabilizing field, the failure mode the puffy disk avoids."},{"cited_title":"D., 1984, J","cited_arxiv_id":null,"evidence_quote":"Introduces the M1 closure for the radiation field used in the simulations."},{"cited_title":"M.\\ 2018, , 857, 34","cited_arxiv_id":null,"evidence_quote":"Non-radiative simulations finding accretion predominantly above the high-density core, corroborating the puffy inflow pattern."}],"review_version":1}