{"id":"fc2fc0ff-667d-4f0b-a2a2-b08fbb474a14","arxiv_id":"2505.13701","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A 3D plasma simulation of a thin disc confirms the analytic geodesic model for the plunging region inside the ISCO, finding a 5.3% angular momentum loss and rejecting constant alpha-disc models there.","lead":"This paper tests a theory of what happens to gas after it crosses the last stable orbit around a black hole, using a complex computer simulation of a thin accretion disc. If the theory holds, black hole spin measurements could become more accurate by including light from gas that current models ignore.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermodynamic 'excellent agreement' rests on per-profile epsilon fits; a no-free-parameter prediction with measured entropy and dynamic epsilon is needed to confirm it.","rationale":"Good-faith reading: the paper is a serious numerical test of MB23, with a purpose-built AthenaK simulation, high resolution (cell aspect 1:1:1, 0.025 r_g midplane), MRI quality factors above standard thresholds, and a lower-resolution run in Appendix A. The dynamical comparison (Fig. 3) is strong: the radial velocity profile matches the offset geodesic with u_I taken from the first interior point, not fitted, and the residuals drop toward the horizon. The angular-momentum drop of about 5.3% is directly visible and consistent with prior work. Credit is due for stating the numerical caveats (cooling function, thermal-energy replacement, ad-hoc K model) explicitly. The weak point is the quantitative thermodynamics. The abstract's 'excellent agreement' is evidenced by Fig. 5, but Table 2 shows the comparison is made with epsilon fitted independently for each of Sigma, rho, P, and T, after m is fitted to K. The fitted epsilons disagree by up to a factor of about 3 (0.015 vs 0.044), while the dynamically measured value is 0.041. The paper's averaging argument may explain this, but it is not tested; if accepted, it makes the model hard to falsify with this procedure. The reader's weakest assumption about the K power law is on point, but the more central issue is that the single physical parameter of the model is allowed to float per quantity. A clean, no-free-parameter prediction using the actual measured K(r) and the dynamic epsilon would settle whether the MB23 functional forms are genuinely predictive. This does not change the verdict: CONDITIONAL remains appropriate. If the no-free-parameter test passes, the paper is a strong confirmation; if it fails, the thermodynamic claims are over-stated but the dynamical results and the qualitative non-vanishing thermodynamics still stand, so rejection would be too strong without further evidence.","tokens_in":17662,"tokens_out":10852,"duration_ms":110960,"concrete_test":"Using the simulated K(r) profile directly (interpolated, not the power-law fit) and the dynamic epsilon = 0.041 from Table 2 / Eq. 2, evaluate Eqs. 5-8 for Sigma, rho, P, and T and overlay the resulting curves on the simulated data in Fig. 5. If all four curves lie within the +/- 1 sigma variance bands, the model is predictive and the concern is resolved; if any curve deviates systematically, the 'excellent agreement' is an artifact of the per-profile epsilon and m fitting.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Table 2 shows that the claimed agreement for rho, P, and T (Fig. 5) is obtained after fitting a separate epsilon for each quantity, with m fixed from the K power-law fit. The best-fit epsilons are 0.021 (Sigma), 0.015 (rho), 0.015 (P), and 0.044 (T), whereas the directly measured dynamic epsilon from the radial velocity is 0.041 +/- 0.002. A factor-of-three spread in the single physical parameter of the MB23 model is not small, and the paper's defense that non-linear averaging makes epsilon an 'effective parameter' is plausible but undemonstrated. Because epsilon is fitted per profile, the comparison is a two-parameter fit to each smooth monotonic curve rather than a predictive test of the model. The entropy power-law issue identified by the reader is real but secondary: even a perfect power-law K would leave the per-profile epsilon freedom. The central claim that the MB23 model quantitatively describes plunging thermodynamics is therefore not yet established; it needs a check with no fitted parameters.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a 3D global GRMHD simulation of a thin accretion disc around a Schwarzschild black hole using the AthenaK code, and compares the simulated radial velocity, angular momentum, surface density, density, pressure, and temperature profiles in the plunging region with the analytic MB23 model. The authors find that the radial velocity is well described by an offset geodesic, that the angular momentum drops by about 5.3% across the plunge, and that the thermodynamic profiles match the MB23 solutions provided a non-adiabatic entropy rise is modeled as a power law with a fitted index. The paper argues that the results support the MB23 model and have implications for black hole spin measurements.","tokens_in":17850,"tokens_out":4771,"duration_ms":44439,"significance":"If the match is robust, the paper provides an independent numerical confirmation of the key assumptions of the MB23 plunging-region model: gravity-dominated geodesic dynamics and a small but finite ISCO stress that keeps thermodynamic quantities non-vanishing. This matters for thermal continuum spin measurements that currently truncate the disc at the ISCO. The simulation is a dedicated high-resolution calculation with careful attention to MRI resolution, a lower-resolution cross-check, and direct measurements of the dynamical quantities. The main weakness is that the thermodynamic agreement is obtained after fitting parameters (epsilon per profile, m for the entropy power law), which reduces the predictive power of the comparison.","major_comments":[{"comment":"The claim of excellent agreement for the thermodynamic quantities rests on fitting a separate epsilon for each profile. In Table 2, the best-fit values are 0.021 (Sigma), 0.015 (rho), 0.015 (P), and 0.044 (T), whereas the directly measured epsilon from the radial velocity is 0.041 +/- 0.002. Since epsilon is the same physical parameter in the MB23 model, this factor-of-three spread is not explained by the quoted error bars. The statement in Section 4 that epsilon should be treated as an effective parameter owing to non-linear averaging is plausible but is not demonstrated. Without such a demonstration, the fits to Eqs. (6)-(8) are two-parameter fits to smooth monotonic profiles and do not provide a predictive test of the model. I request a no-free-parameter check, for example using the measured epsilon to predict the Sigma, rho, P, and T profiles without fitting, or a joint fit with a single epsilon, or a synthetic-data demonstration that the averaging effect can produce the observed spread.","section":"Section 4, Table 2, Eqs. (5)-(8)"},{"comment":"The entropy rise is modeled as K = K_I (r/r_I)^(-m) with m fitted to the simulation (m = 2.71 in the high-resolution run, m = 2.08 in the lower-resolution run). This is an ad-hoc prescription. The agreement between the MB23 thermodynamic solutions and the simulation is conditional on this fitted profile. The paper attributes the entropy rise to magnetic reconnection but does not provide a quantitative estimate of the dissipation rate from the simulation that can be compared with the fitted m. I recommend either deriving m from a physical heating model or explicitly framing the thermodynamic test as conditional on the measured entropy profile, and additionally showing the sensitivity of the density, pressure, and temperature fits to the choice of m.","section":"Section 4, Fig. 4, Table 2, entropy power law"},{"comment":"The fitted epsilon values are not stable between the high-resolution and lower-resolution simulations. For example, Sigma gives epsilon = 0.021 in Table 2 but epsilon = 0.051 in Table A2, and rho gives 0.015 versus 0.031. This sensitivity suggests that the per-profile fits are not robust, reinforcing the need for a parameter-free test before the central claim can be accepted.","section":"Section 4, Fig. 5, and Appendix A"}],"minor_comments":[{"comment":"The phrase 'do not to include the plunging fluid' in the abstract and introduction should read 'do not include the plunging fluid'.","section":"Abstract and Section 1"},{"comment":"The paper consistently misspells 'Schwarzschild' as 'Schwarzchild' in a few places, for example in Section 3 and in the Conclusions.","section":"Section 3"},{"comment":"The caption begins 'Theplungingregion' without a space; this is a typographical error.","section":"Figure 7 caption"},{"comment":"The text refers to 'dot-dash lines' for the geodesic solutions, but the offset and pure geodesic models are both plotted as dot-dash lines in different colors; the caption could clarify which color corresponds to which model.","section":"Section 4, Fig. 3 caption"},{"comment":"The data availability statement says numerical results will be shared upon reasonable request; for reproducibility it would be helpful to also release the analysis scripts or key derived profiles alongside the paper.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of MNRAS. The MB23 model was developed in part by two of the co-authors, and the observational consistency checks cited are from the same group, so the thermodynamic comparison is not fully independent. The dynamical results (geodesic velocity match, angular momentum drop) are direct and credible, but the thermodynamic claim needs the requested no-free-parameter test. I recommend inviting a revision with that test rather than accepting the current version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jake and colleagues have done something worthwhile here: they ran a genuinely thin (h/r ~ 0.06) 3D GRMHD disc simulation and used it to test Mummery & Balbus's analytic plunging-region model. The dynamical part of the test is strong and direct. The measured radial 4-velocity tracks the offset geodesic, with the offset taken from the first point inside the ISCO rather than fitted. The 5.3% angular momentum drop is a clean measurement, consistent with earlier simulations and with the observational inference of ~4%. And the order-of-magnitude rise in alpha across the plunge is exactly what you'd expect from flux freezing, so the claim that constant-alpha models are wrong there is well supported. The lower-resolution run backs up the key features.\n\nThe soft spot is the thermodynamics. The model comparison requires two fitted ingredients: m from a power-law fit to the entropy K, and then a separate epsilon fit for each of Sigma, rho, P and T. The best-fit epsilons scatter from 0.015 to 0.044 while the dynamically measured epsilon is 0.041. That's a factor-of-three spread in what is supposed to be a single parameter. The paper's explanation — that nonlinear averaging makes epsilon an effective parameter — is plausible but not demonstrated. So the 'excellent agreement' in Fig. 5 is really a two-parameter fit to each smooth monotonic curve, not an independent prediction of the MB23 model. That doesn't kill the paper, but it should be said plainly. The authors are transparent about it, but the abstract's 'excellent agreement' oversells what is currently a conditional success.\n\nThe other minor issues: the ad hoc cooling function and the cell-averaging of thermal energy in the plunge are standard but still caveats, and the data are not public. None of these undermine the dynamical results.\n\nI'd send this to peer review. The right referee will make them tighten the thermodynamic claim, ideally by testing with epsilon fixed to the measured value or by deriving the entropy rise from a physical model. I would cite the paper for the ISCO stress and the geodesic dynamics.","headline":"Solid GRMHD test of the plunging-region dynamics; the thermodynamic agreement is real but fit-dependent and should not be oversold.","tokens_in":18398,"tokens_out":2856,"would_cite":true,"duration_ms":26259,"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":"A global 3D GRMHD simulation confirms the MB23 analytic geodesic model for the thermodynamics of the plunging region of a thin Schwarzschild disc, provided non-adiabatic heating is included, and measures a finite ISCO stress of δJ ≈ 5.3%.","keywords":["black hole accretion","plunging region","ISCO stress","GRMHD simulation","thin accretion disc","thermal continuum spin measurement","geodesic plunge"],"falsifier":"Run the same physical setup at two or more grid resolutions and with different cooling prescriptions, then compare the fitted entropy index $m$ and angular-momentum drop $\\delta_\\mathcal{J}$; if $m$ changes by more than the quoted uncertainties, the heating law is numerical rather than physical and the thermodynamic confirmation fails.","tokens_in":17404,"feed_emoji":"🕳️","tokens_out":10388,"duration_ms":91720,"temperature":0.7,"pith_summary":"The paper tests whether analytic formulas for the plunging region—the zone between the innermost stable circular orbit (ISCO) of a black hole and its event horizon, where matter falls in almost freely—can describe a realistic thin accretion disc. It reports that they can, provided the flow is not treated as adiabatic: a small radial rise in entropy, modelled as a power law in radius, captures heating attributed to magnetic reconnection. The central finding is a finite but small ISCO stress, corresponding to a ~5.3% drop in angular momentum across the plunge, which barely perturbs the geodesic infall yet keeps density, pressure, and temperature from vanishing near the ISCO. Because black-hole spin measurements from disc spectra commonly truncate the disc at the ISCO, this means the plunging fluid can radiate and should be included.","feed_headline":"Plunging disc gas follows a simple infall model to the horizon","feed_subtitle":"A 5.3% angular-momentum drop keeps gas hot past the ISCO, so spin fits should include it","key_machinery":"The central object is the MB23 offset-geodesic plunging model: a radial 4-velocity $U^r = -c \\sqrt{2 r_g/(3 r_I)} \\left( r_I/r - 1 \\right)^{3/2} - u_I$, where $u_I$ is the small inward speed at the ISCO, written in dimensionless form as $\\epsilon = (u_I/c)\\sqrt{3 r_I/(2 r_g)}$. This closes the equations through mass conservation, vertical hydrostatic equilibrium, and the entropy relation $P = K \\rho^\\gamma$, producing self-similar profiles for the surface density, density, pressure, temperature, and scale height. The paper adds a fitted radial power law $K = K_I (r/r_I)^{-m}$ with $m \\approx 2.71$ to account for non-adiabatic heating, and uses the measured angular-momentum drop $\\delta_\\mathcal{J} \\approx 5.3\\%$ as the diagnostic of the ISCO stress. The ideal-GRMHD simulation supplies the independent numerical data against which these profiles are compared.","core_discovery":"On the paper's own terms, a dedicated 3D general-relativistic magnetohydrodynamic simulation of a thin, weakly magnetised disc around a Schwarzschild black hole shows that the MB23 analytic geodesic-plunge model quantitatively reproduces the plunging-region thermodynamics. The simulated radial 4-velocity follows the offset-geodesic solution, and the angular momentum drops by δJ ≈ 5.3% from ISCO to horizon. The surface density, density, pressure, central temperature, and scale height all agree with the model's self-similar profiles once non-adiabatic heating is represented as a power-law entropy rise K = K_I (r/r_I)^(-m) with fitted index m ≈ 2.71. The authors identify the heating as grid-scale magnetic reconnection in a mid-plane current sheet. The stress is small enough that the plunge remains essentially geodesic, but large enough to dissipate energy near the ISCO and prevent the thermodynamic quantities from vanishing. The paper concludes that constant-α disc models are physically inappropriate inside the plunging region.","pith_inferences":["The power-law entropy description is a pragmatic interpolation; a more physical sub-grid model of magnetic-reconnection heating would let the fitted index $m$ be predicted rather than measured.","If the heating is really grid-scale reconnection driven by flux freezing, the entropy rise should track the magnetic-field amplification set by the plunge, so simulations with different initial field geometries should yield different $m$ values and possibly spin-dependent ISCO stress.","The near-ISCO dissipation implied by the 5.3% stress suggests that the effective inner boundary condition for outer-disc models should be a non-zero stress carrying roughly 5% of the local angular-momentum flux, not the traditional zero-stress condition."],"forward_implications":["Thermal continuum spin measurements that truncate the disc at the ISCO omit a region whose thermodynamic quantities stay non-zero; including the plunging fluid should shift fitted black-hole spins and high-energy spectral tails.","The simulated $\\delta_\\mathcal{J} \\approx 5.3\\%$ matches the roughly 4% inferred from observations of MAXI J1820+070, supporting the reality of a finite ISCO stress.","Constant-$\\alpha$ disc models cannot describe the plunging region, since $\\alpha$ rises by an order of magnitude and the local stress-dissipation coupling fails there.","Because the scaled analytic profiles depend only on $r/r_I$, the same model with $r_I$ set by spin should apply to Kerr black holes as well as Schwarzschild."],"supporting_citations":[{"why":"Supplies the analytic geodesic-plunge thermodynamic model that the simulation is designed to test.","marker":"MB23"},{"why":"Earlier GRMHD test of the same analytic model for thick discs, which this work extends to thin discs.","marker":"Mummery & Stone (2024)"},{"why":"Provides the cooling prescription used to keep the disc thin and reports an earlier finite ISCO stress.","marker":"Penna et al. (2010)"},{"why":"Prior thin-disc GRMHD simulations of the plunging region that the present results are compared with.","marker":"Noble et al. (2010)"},{"why":"Source of the cooling-function approach and an earlier indication of non-zero ISCO stress.","marker":"Shafee et al. (2008)"},{"why":"Observational continuum-fitting result giving $\\delta_\\mathcal{J} \\approx 4\\%$ for MAXI J1820+070, against which the simulated 5.3% is compared.","marker":"Mummery et al. (2024b)"},{"why":"The relativistic thin-disc model with a zero-stress ISCO boundary condition that the paper argues should be replaced.","marker":"Novikov & Thorne (1973)"},{"why":"Argued that magnetic fields can maintain a finite stress at the ISCO, the physical basis for the non-zero stress measured here.","marker":"Krolik (1999)"}],"fun_headline_variants":["Plunge region: 5.3% angular momentum drop keeps gas hot","Plunging disc: 5.3% spin drop, geodesic flow, hot ISCO","Plunge region: magnetic heating, geodesic flow, 5.3% spin loss","GRMHD confirms plunge thermodynamics: 5.3% spin drop","Constant-α discs fail inside plunging region"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The thermodynamic agreement rests on modelling the extra heating as a power-law curve in radius with an index fitted to the simulation; if the real heating profile has a different shape, the reported match for density, pressure, and temperature could be a fitting artifact rather than a confirmation.","fun_headline_variants_meta":{"raw":{"variants":["Plunge region: 5.3% angular momentum drop keeps gas hot","Plunging disc: 5.3% spin drop, geodesic flow, hot ISCO","Plunge region: magnetic heating, geodesic flow, 5.3% spin loss","GRMHD confirms plunge thermodynamics: 5.3% spin drop","Constant-α discs fail inside plunging region"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001211,"raw_usage":{"total_tokens":5014,"prompt_tokens":1005,"completion_tokens":4009,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":621,"completion_tokens_details":{"reasoning_tokens":3907}},"tokens_in":621,"tokens_out":4009,"duration_ms":29789,"temperature":1.0,"reasoning_tokens":3907,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:11:00.942031+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same physical setup at two or more grid resolutions and with different cooling prescriptions, then compare the fitted entropy index $m$ and angular-momentum drop $\\delta_\\mathcal{J}$; if $m$ changes by more than the quoted uncertainties, the heating law is numerical rather than physical and the thermodynamic confirmation fails.","supporting_citations":[{"cited_title":"D., Thorne K","cited_arxiv_id":null,"evidence_quote":"The relativistic thin-disc model with a zero-stress ISCO boundary condition that the paper argues should be replaced."}],"review_version":1}