{"id":"48c50e4d-de30-452a-b077-7011579cf93c","arxiv_id":"2411.12647","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":7,"one_line_summary":"A new public GRMHD simulation library of MAD and SANE accretion flows across five black hole spins confirms SANE spin equilibrium near a*~0.94 and MAD jet-powered spin-down.","lead":"This paper releases a library of 10 general relativistic magnetohydrodynamic simulations of black hole accretion, varying black hole spin and magnetic flux state, and reports how angular momentum, energy, and jet power depend on those parameters. It is a public resource for interpreting Event Horizon Telescope observations of Sgr A* and other low-luminosity black holes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The SANE spin-equilibrium value aeq~0.94 is not actually supported by Eq. 32, whose zero is at a*~0.16; the claim rests on an unquantified interpolation between a*=+0.5 and +0.94.","rationale":"The reader's CONDITIONAL verdict already flags the sparsity of spin sampling and the lack of a quoted uncertainty for Eq. 32, so my concern does not move the verdict. However, my emphasis differs: the reader's weakest assumption is MRI resolution, whereas I find the more immediate logical issue in the paper's own text. Equation 32, as written, has a zero near a*=0.16; if it is meant to support the SANE equilibrium at 0.94, it cannot. If it is meant for MAD models, the text should say so, and the SANE equilibrium should be derived and quantified separately. The public data release and explicit limitation statements are good faith and support the qualitative conclusions, but the precise value aeq~0.94 is an interpolation between two widely spaced spin points without propagated uncertainty. A bootstrap of the time series and one additional SANE run near a*=0.8 would settle whether 0.94 is robust or merely consistent with prior literature. The qualitative claims, including prograde MAD spin-down and outflow efficiency greater than one for a*=+0.94, are supported by the data and by prior work, so a rejection is not warranted.","tokens_in":30466,"tokens_out":10193,"duration_ms":104862,"concrete_test":"Bootstrap the SANE zero crossing: split t=[15,30] x 10^3 tg into at least 30 overlapping sub-intervals, recompute l and e at the horizon for the five SANE runs, and fit a monotone curve to s(a*) for each sub-sample; report the 95% confidence interval for the zero. Independently, add or identify one published SANE simulation at a* between 0.75 and 0.85 and evaluate its s using the same integrals. If the zero moves outside [0.88,0.98] or the bootstrap interval width exceeds 0.1, the aeq~0.94 claim should be revised to a range, and if Eq. 32 is intended as a MAD fit this should be stated explicitly with a separate SANE fit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline quantitative result is the SANE spin-equilibrium value aeq~0.94 in Section 4.6, but the only explicit fit in that section, Eq. 32, has sfit = -3.267a^2 - 8.349a + 1.387, whose positive root is a* approximately 0.157, not 0.94. If the gray curve in Figure 15 is meant to represent the SANE crossing, the equation is inconsistent with the stated equilibrium; if it is instead a MAD fit, then no fit or uncertainty is provided for the SANE crossing. The SANE equilibrium is therefore inferred from the sign change of s between a*=+0.5 and a*=+0.94, with only five spin points and error bars that measure time variability, not systematic or resolution uncertainty. Because s = l - 2a*e, small offsets in l or e are amplified by the factor 2a* near a*=0.94, so a modest systematic error in the measured horizon fluxes could shift the inferred equilibrium by 0.1 or more. The reader's MRI-resolution concern is real and relevant, but the more immediate logical gap is that the paper's own fit does not substantiate the stated aeq, and the sparse spin sampling does not provide a defensible interpolation to 0.94. This does not undermine the qualitative SANE-vs-MAD distinction, which is consistent with prior work, but it does weaken the precision of the central mapping claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a library of ten ideal GRMHD simulations of black hole accretion disks, generated with the KHARMA code, spanning five spins (a* = -15/16, -1/2, 0, +1/2, +15/16) and two magnetization states (SANE and MAD). The authors analyze horizon-penetrating mass, angular momentum, and energy fluxes; inflow equilibrium; time-averaged disk structure; conserved currents; jet power; and black hole spin-up/spin-down. The central claims are that SANE disks closely follow thin-disk angular momentum and energy flux expectations and reach spin equilibrium at a* ~ 0.94, while prograde MAD models experience jet-driven spin-down, with the a* = +0.94 MAD having outflow efficiency greater than 1. The simulation data are publicly released, and the paper emphasizes its use in the EHT analysis of Sgr A*.","tokens_in":30856,"tokens_out":6898,"duration_ms":69709,"significance":"If the quantitative claims hold, this library provides a direct numerical mapping from (a*, magnetization) to horizon fluxes and jet powers, which is valuable for interpreting event-horizon-scale observations. The paper is strong in its transparency: it reports MRI quality factors, inflow equilibrium radii, density floor statistics, and measurement choices, and the code and data are publicly available. The qualitative SANE-versus-MAD distinction is robust and consistent with prior work. However, the precision of the headline SANE spin-equilibrium value is currently not supported by the paper's own fit, and the resolution and measurement-radius choices introduce unquantified systematic uncertainties. These issues are local and fixable, but they affect the central mapping claim.","major_comments":[{"comment":"The quadratic fit s_fit = -3.267 a*^2 - 8.349 a* + 1.387 has a positive root at a* ~ 0.157, not at a* ~ 0.94. The text states that SANE models achieve spin equilibrium at aeq ~ 0.94, but the provided fit does not substantiate this. If Eq. (32) is intended as a fit to the SANE points, it contradicts the claimed crossing; if it is a fit to the combined sample, it does not provide the SANE-specific crossing. The authors should present a SANE-only fit (or a clear interpolation) with quantified uncertainty, or explicitly state that the equilibrium is inferred by interpolation between a* = +0.5 and a* = +0.94 and justify that inference.","section":"§4.6, Eq. (32), Figure 15"},{"comment":"The statement 'l > 0 for all spins indicates a net inflow of angular momentum into the black hole' conflicts with the signed definition of l given in the same section. For retrograde spins (a* < 0), the inflowing gas has negative specific angular momentum, so l should be negative. This inconsistency affects the interpretation of Figure 3 and the spin-up analysis in §4.6, and it should be corrected or clarified.","section":"§4.1.2, first paragraph"},{"comment":"The MRI quality factors in the SANE simulations are Q_theta ~ 5-10 and Q_phi ~ 12-16, below the nominal values of ~10 and ~20, and the paper concedes that the resolution is 'somewhat lower than contemporary studies.' Because the spin-equilibrium result depends on the turbulent stresses that set l and e, the quantitative value aeq ~ 0.94 carries an unquantified resolution uncertainty. The authors should either provide a resolution study of l, e, and s for at least one spin, or explicitly temper the precision of the equilibrium claim.","section":"§3.6, §6"},{"comment":"The inflow equilibrium radii are req ~ 18-32 rg for SANE and ~45-66 rg for MAD, yet the jet powers in Table 3 are evaluated at r = 100 rg. For the MAD +0.94 model with outflow efficiency > 1, the region at 100 rg is beyond the inflow equilibrium radius, so the reported Pjet may not be converged. Please quantify the time variability of Pjet at 100 rg or measure it at a radius within the converged region.","section":"§4.2, Table 2 vs §4.5, Table 3"}],"minor_comments":[{"comment":"The abstract states '30,000 GM/c^3' where the time unit should be written consistently as GM/c^3 (or tg); the use of GM/c^3 vs GM/c^2 is confusing in a few places.","section":"Abstract and text"},{"comment":"The modulation index M3 is defined, but the error estimate quoted as sigma_M3 / sqrt(27) is not fully explained; please clarify how the number of independent 3-hour segments is derived.","section":"§4.1.3"},{"comment":"The columns req(Mdot) and req(tin) are not defined in the table note; they are defined in the text but a brief note would aid readability.","section":"Table 2"},{"comment":"The '1DW' primitive recovery scheme is mentioned without a reference; consider adding Noble et al. (2006) or Mignone & McKinney (2007) at that point.","section":"§3.4"}],"recommendation":"major_revision","confidential_remarks":"This is a useful data-release and benchmarking paper, and the numerical work appears careful. The central quantitative claim (SANE spin equilibrium at a* ~ 0.94) currently conflicts with the presented fit, and the resolution and measurement-radius issues need to be addressed quantitatively. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection. The sign inconsistency in §4.1.2 is likely a simple error but should be resolved before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Ben, here's my read on arXiv:2411.12647. The v3 GRMHD library is a genuinely useful public resource, and the paper is mostly honest about its numerical limitations. But the headline SANE spin-equilibrium number, a_* ~ 0.94, is not supported by the paper's own fit. Equation 32 has a positive root around a_* ~ 0.16, not 0.94. That inconsistency should be fixed before publication.\n\nWhat is actually new: the library itself—ten long GRMHD runs spanning five spins and two magnetization states, publicly released and already used in EHT Sgr A* analysis—and the conserved-current analysis of mass, angular momentum, and energy transport. The paper is careful about floors, ceilings, MRI quality factors, and inflow-equilibrium radii. The authors are transparent that resolution is lower than some contemporary studies and that inflow equilibrium is limited. The qualitative SANE/MAD distinction is robust and consistent with earlier work.\n\nSoft spots. The spin-equilibrium claim is the main one. The abstract and Section 4.6 assert a_eq ~ 0.94, but the quadratic fit they provide crosses zero at a_* ~ 0.16. They might have fitted all ten models together, or meant the fit for MAD only, but as written the equation contradicts the text. The actual case for 0.94 is a sign change between a_* = +0.5 and +0.94, with no uncertainty quoted. Since s = l - 2 a_* e, small offsets in l or e are amplified near a_* = 0.94, so the precision implied by \"0.94\" is not warranted. This is a fixable presentation issue, not a fatal flaw. Second, the MRI resolution is below nominal values (Q_theta ~ 5-10 versus ~10-20). The authors acknowledge it, but it does weaken the quantitative flux values and the spin-equilibrium estimate.\n\nWho this is for: EHT interpreters, GRMHD practitioners, anyone building simulation libraries. It deserves a serious referee. I would send it to a competent referee with a note to check the spin-equilibrium fit and to ask for an explicit statement that 0.94 is a coarse estimate. The paper has real value as a data release and a benchmark.","headline":"The v3 GRMHD library is a genuinely useful public resource, but the headline spin-equilibrium claim is undermined by an internal inconsistency between the stated a_eq ~ 0.94 and the provided quadratic fit.","tokens_in":31358,"tokens_out":3668,"would_cite":false,"duration_ms":34860,"reading_group":"maybe","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 SANE accretion flows spin black holes up to a*≈0.94, while magnetically arrested jets can spin them back down.","keywords":["general relativistic magnetohydrodynamics","black hole accretion","black hole spin equilibrium","magnetically arrested disks","SANE disks","Blandford-Znajek jets","Sagittarius A*","numerical simulations"],"falsifier":"Rerun the SANE $a_* = +0.94$ simulation at double resolution so the MRI quality factors reach the recommended values, recompute the time-averaged spin-up parameter $s$, and check whether the zero crossing stays at $a_* \\approx 0.94$; if it shifts by more than the error bars, the claimed equilibrium is a resolution artifact.","tokens_in":30271,"feed_emoji":"🕳️","tokens_out":7557,"duration_ms":75185,"temperature":0.7,"pith_summary":"This paper tries to establish that a small, carefully chosen library of ten ideal GRMHD simulations can map the two parameters that most control black hole accretion—spin and magnetic flux threading the horizon—onto the fluxes, jet power, and variability that observers can measure. The low-magnetization SANE flows behave like thin disks: their angular momentum and energy fluxes track the standard disk expectation, and they drive the black hole to a spin equilibrium near $a_* \\sim 0.94$. The magnetically arrested MAD flows behave differently: prograde cases launch jets that extract angular momentum so efficiently that the black hole spins down, and the fastest-spinning MAD model has an outflow efficiency greater than one, meaning the jet carries away more energy than the accreting gas supplies. The paper also maps the conserved currents of mass, angular momentum, and energy, and shows that MAD disks are sub-Keplerian, hotter, and more variable than SANE disks. These results matter because this is the simulation library used to interpret horizon-scale millimeter observations of the galactic center and other nearby supermassive black holes.","feed_headline":"SANE disks spin black holes up; MAD jets spin them down","feed_subtitle":"Ten simulations show low-magnetization disks reach spin 0.94, while magnetically arrested jets can reverse the spin.","key_machinery":"The load-bearing machinery is the set of conserved currents built from the Killing vectors of the Kerr metric: the mass current $J_M = \\sqrt{-g}\\,\\rho u^\\mu$, the angular momentum current $J_L = \\sqrt{-g}\\,T^\\mu{}_\\phi$, and the energy current $J_E = -\\sqrt{-g}\\,T^\\mu{}_t$. These are decomposed into electromagnetic and fluid contributions, integrated at the horizon to give specific fluxes $l$ and $e$, and combined into the spin-up parameter $s = l - 2a_*e$, whose sign and zero crossing determine whether a black hole gains or loses spin. The same currents, plotted as time- and azimuth-averaged streamlines, make the Blandford-Znajek extraction pattern visible: inflow in the disk, collimated outflow in the jet.","core_discovery":"The central discovery is that accretion state, not just spin, controls how a black hole evolves. In the SANE models the time-averaged specific angular momentum flux at the horizon matches the thin-disk value at the innermost stable circular orbit, and the spin-up parameter $s = l - 2a_*e$ tracks the thin-disk curve, crossing zero near $a_* \\sim 0.94$; the deviations that do appear in prograde models come from fluid thermodynamic forces. In the prograde MAD models the electromagnetic part of the horizon flux dominates, outward Poynting flux is visible along parabolic jet contours anchored at the horizon, and the $a_* = +0.94$ model reaches an outflow efficiency of about $1.5$, so the black hole loses spin. The same data show MAD disks are sub-Keplerian, roughly an order of magnitude hotter than SANE disks inside $r \\lesssim 10\\,GM/c^2$, and more variable in accretion rate because of intermittent flux eruption events.","pith_inferences":["I would treat the $a_* \\approx 0.94$ equilibrium as provisional until the SANE runs are repeated at resolution high enough to meet the recommended MRI quality factors, since the paper itself reports being below those levels.","A longer-term inference is that black hole spin may settle at an intermediate value if a source alternates between SANE and MAD epochs, because the two states drive spin in opposite directions.","The paper's choice of adiabatic index and omission of electron thermodynamics leaves room for the relative temperature of MAD and SANE disks to shift, which would change the synthetic images built from this library.","The public data release makes a direct test possible: generate synthetic polarized images from the MAD and SANE snapshots and compare predicted image signatures to upcoming higher-resolution observations."],"forward_implications":["If the SANE result holds, any low-magnetization accretion flow will push a black hole toward $a_* \\approx 0.94$, regardless of the spin it started with.","If the MAD result holds, a prograde, magnetically arrested flow can spin the black hole down, so the observed spin of a source encodes its accretion history, not just its initial spin.","An outflow efficiency above one for the highest-spin MAD model means a jet can carry away more mechanical energy than the accreting rest mass supplies, which sets a strict upper limit on the jet power available to accelerate particles.","The measured relation between magnetization, spin, and horizon fluxes gives observers a translation table: a horizon image and light curve can be inverted into a preferred accretion state and spin.","The strong variability difference between MAD and SANE accretion rates implies that light-curve variability is a cheap observable discriminator between the two states."],"supporting_citations":[{"why":"This paper sets the thin-disk spin-up starting point by showing a cold equatorial disk drives the black hole toward maximal spin.","marker":"Bardeen (1970)"},{"why":"This paper adds the photon-capture torque and gives the limiting spin near 0.998 that the SANE thin-disk expectation builds on.","marker":"Thorne (1974)"},{"why":"This paper supplies the earlier GRMHD spin-equilibrium result at a* ≈ 0.94 that the SANE models reproduce.","marker":"Gammie et al. (2004)"},{"why":"This paper provides the mechanism the MAD simulations invoke for outward energy and angular momentum extraction along the jet.","marker":"Blandford & Znajek (1977)"},{"why":"This paper gives the earlier MAD spin-down result that the prograde MAD models are compared against.","marker":"Tchekhovskoy et al. (2012)"},{"why":"This paper supplies the lower MAD spin equilibrium and the fifth-order spin-up fit used as a comparison in this library.","marker":"Narayan et al. (2022)"},{"why":"This paper defines the recommended MRI quality factors used to assess whether the simulations resolve the turbulent transport.","marker":"Hawley et al. (2011)"},{"why":"This paper provides the semi-analytic relativistic viscous flow model that the SANE disk temperature, angular momentum, and scale-height profiles are compared with.","marker":"Popham & Gammie (1998)"}],"fun_headline_variants":["MAD jets spin black holes down; SANE disks spin them up","Accretion state flips black hole spin: SANE up, MAD down","SANE accretion spins black holes to 0.94; MAD jets reverse it","Black hole spin outcome hinges on SANE vs MAD accretion state"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the MRI turbulence that transports angular momentum is resolved well enough that the measured fluxes do not change if the grid is refined; the paper reports MRI quality factors below the recommended values.","fun_headline_variants_meta":{"raw":{"variants":["MAD jets spin black holes down; SANE disks spin them up","Accretion state flips black hole spin: SANE up, MAD down","SANE accretion spins black holes to 0.94; MAD jets reverse it","Black hole spin outcome hinges on SANE vs MAD accretion state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001144,"raw_usage":{"total_tokens":4842,"prompt_tokens":1133,"completion_tokens":3709,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":749,"completion_tokens_details":{"reasoning_tokens":3628}},"tokens_in":749,"tokens_out":3709,"duration_ms":21776,"temperature":1.0,"reasoning_tokens":3628,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:18:05.493179+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the SANE $a_* = +0.94$ simulation at double resolution so the MRI quality factors reach the recommended values, recompute the time-averaged spin-up parameter $s$, and check whether the zero crossing stays at $a_* \\approx 0.94$; if it shifts by more than the error bars, the claimed equilibrium is a resolution artifact.","supporting_citations":[],"review_version":1}