{"id":"e83215a4-4491-411b-beff-d2c7328f7c83","arxiv_id":"2507.22506","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"GRMHD simulations of low-angular-momentum accretion onto rapidly spinning black holes show quasi-steady standing shocks in the weakly magnetized SANE regime, with shock position shifting outward for higher spin.","lead":"This paper uses general-relativistic magnetohydrodynamic simulations to study gas with low angular momentum falling onto a spinning black hole. It reports stable standing shocks near the black hole, a feature predicted analytically for decades but never before seen in these magnetized simulations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'first GRMHD evidence' claim rests on 2.5D runs in a regime where the authors state MRI growth stalls; the 3D validation is at α=10, so α=50 standing shocks are not shown to survive full turbulence.","rationale":"The reader's weakest assumption identifies the same load-bearing issue: the main results rely on 2.5D axisymmetric simulations in a late-time window where MRI-driven turbulence is not sustained, and the only 3D test is at a different, lower-α parameter point. That is the critical condition for the paper's central claim of 'first GRMHD evidence' of standing shocks. If the α=50 shocks are artifacts of 2D axisymmetry or of the density-confinement parameter, the strong novelty claim is not established; if they survive a full 3D run with sustained turbulence, the paper's conclusion is substantially supported. I agree with the reader's conditional assessment: the result is plausible and important if correct, but the evidence is incomplete. A targeted 3D α=50 run with turbulence diagnostics is the single check that would settle the concern. I therefore recommend no change to the reader's verdict.","tokens_in":14298,"tokens_out":4748,"duration_ms":57492,"concrete_test":"Run a full 3D GRMHD simulation of the α=50, a*=0.9 case with a resolution comparable to the 2.5D runs (e.g., at least 512×128×64, or a short resolution study from 256×80×64 upward) and evolve past 15,000 tg. During t=12,000–15,000 tg, measure equatorial radial velocity and pressure profiles, locate the discontinuities, and compute Maxwell and Reynolds stresses to confirm MRI-driven turbulence is sustained. If the two quasi-steady jumps near r≈4–5 rg and r≈40–50 rg persist with small radial variability and the flow remains turbulent, the central claim is strengthened; if the shocks disappear or oscillate chaotically in 3D, the 2.5D α=50 result is not evidence for standing shocks in full GRMHD.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the quasi-steady discontinuities are genuine standing shocks in turbulent, weakly magnetized GRMHD accretion. The paper's Section 3 states that after t ≳ 10,000 tg, 'MRI can not continuously grow the magnetic fields due to the 2D nature of the simulations', yet the shock analysis in Section 4.1 is time-averaged over 12,000–15,000 tg, precisely inside this decayed-MRI window. For α=50, where the shocks are strongest, the flow may be effectively laminar, and the features could be centrifugal pressure-balance structures produced by the ad hoc α=50 density confinement rather than evidence that shock formation is robust in turbulent GRMHD accretion. The 3D validation in Section 5 uses α=10, not α=30 or 50, at a lower resolution of 256×80×64 compared with the 2.5D runs, and reports only an oscillating single shock at r≈37±5 rg with no turbulence-stress diagnostics. Thus, even if the α=10 3D run shows a shock, it does not establish the α=50 2.5D result that carries the headline claim. Additionally, no Rankine-Hugoniot or Mach-number verification is provided; the radial velocity and pressure jumps in Figure 6 are consistent with shocks but are not proof. The claim is plausible, but the current evidence does not yet support it at the advertised strength.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 2.5D axisymmetric and one 3D GRMHD simulations, performed with BHAC, of low-angular-momentum accretion onto Kerr black holes. A sub-Keplerian parameter F=0.4 and an inclination parameter α (5, 10, 30, 50) are used, with black hole spins a*=0.8, 0.9, and 0.94. The central claim is that quasi-steady standing shocks form in the weakly magnetized SANE regime—specifically a two-shock structure around r≈4–5 rg and r≈40–50 rg for α=50—and that this is the first GRMHD demonstration of such standing shocks. The outer shock location is reported to move outward with increasing spin, and a 3D run at α=10 is presented as validation, showing an oscillating shock at r≈37±5 rg.","tokens_in":14521,"tokens_out":2726,"duration_ms":32231,"significance":"If the claim is correct, the paper would be a significant advance: it would show that standing shocks, previously obtained in idealized 1.5D/2D hydrodynamic and semi-analytic models, can survive in a genuinely multi-dimensional, magnetized GRMHD flow, with direct consequences for interpreting Sgr A* flares and QPOs in X-ray binaries. The paper has notable strengths: it uses the public BHAC code rather than a purpose-built solver, no free parameters are fitted to the claimed shock location, the SANE regime is checked via normalized magnetic flux, and the parameter study spans several spins and inclination parameters. The 3D validation, although limited, is an honest attempt to go beyond axisymmetry. However, the advertised strength of the claim currently exceeds what the presented diagnostics establish.","major_comments":[{"comment":"Section 3 states that after t≳10,000 tg the MRI cannot continuously grow the magnetic fields due to the 2D nature of the simulations, yet the shock analysis in Section 4 is time-averaged over t=12,000–15,000 tg and Figure 7 uses t=14,000–15,000 tg. This means the shocks are identified inside a window where the flow is explicitly described as having decaying MRI, so the repeated characterization of the flow as 'turbulent, weakly magnetized' in Section 4.1 is not supported for the analyzed epoch. Please provide turbulence diagnostics (e.g., Maxwell stress, magnetic energy evolution, α-viscosity parameter) over the averaging window, or identify the shocks in an earlier epoch where MRI is active, or explicitly restrict the claim to a quasi-steady but non-turbulent regime.","section":"Section 3 and Section 4"},{"comment":"The 3D validation is performed for α=10 at resolution 256×80×64, while the headline two-shock result is obtained in 2.5D for α=50 (and partially α=30) at an effective resolution of 1024×512. The α=10 3D run therefore does not validate the α=50 case, and the statement 'we expect to see a similar shock structure for higher values' is an extrapolation. To support the central claim, either run a 3D simulation with α=50 (even at reduced resolution), or provide a concrete physical argument, backed by diagnostics, for why the α=10 3D result transfers to α=50.","section":"Section 5"},{"comment":"The identification of features as shocks rests entirely on visual discontinuities in the radial velocity and pressure along the equatorial plane at a single time (Figure 6) and in ten slices over 14,000–15,000 tg (Figure 7). No Rankine-Hugoniot jump conditions, Mach-number or magnetosonic-Mach-number checks, compression-ratio tests, or comparison with the sonic-point location are provided. Without such verification, the discontinuities could be, for example, centrifugal pressure-balance layers or numerical artifacts. Please add a quantitative shock test, at minimum an upstream Mach number greater than unity and consistency with the expected density/temperature jumps, and ideally a resolution study showing the jump sharpens or converges with resolution.","section":"Section 4.1, Figure 6 and Figure 7"}],"minor_comments":[{"comment":"The paper interchangeably uses '2D', '2.5D', and 'axisymmetric (2.5D)'; please define the terminology once and use it consistently, since 2.5D usually implies axisymmetry with phi derivatives retained, whereas 2D can be ambiguous.","section":"Throughout"},{"comment":"The text says the 3D snapshots are at ts=8000, 9000, and 10,000, while the following paragraph describes t=10,000–12,000 tg and Figure 9 mentions t=10000 and 12000; the time ranges are inconsistent and should be reconciled.","section":"Section 5 and Figure 9"},{"comment":"The text uses a*=0.94, while Figure 3's caption and panels label a*=0.9375; please use one value throughout or explain the difference.","section":"Section 3 and Figure 3"},{"comment":"There is a typo and spacing issue in 'aroundr = 40−50 rg' and the phrase 'by (Nakayama 1994)' should read 'by Nakayama (1994)'.","section":"Section 4.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of ApJ and addresses a timely topic, but the central 'first GRMHD evidence' claim currently rests on 2.5D runs analyzed after MRI has stalled, with a 3D validation at a different parameter point. The required additional diagnostics (turbulence indicators, jump-condition checks, and preferably a 3D high-α run) are substantial but not unreasonable for a revision. I would advise the editor that the 'first' claim should be softened unless the verification is substantially strengthened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Low-angular-momentum accretion is a well-trodden place for semi-analytic shock solutions, and GRHD simulations have seen them, but GRMHD has been genuinely unclear. If standing shocks really form in turbulent, weakly magnetized flows, that is a real bridge between an old paradigm and contemporary simulation. What the paper adds: BHAC runs in 2.5D across several spins, a new inclination parameter alpha that funnels material toward the equator, and one 3D run. The shocks show up as sharp radial velocity and pressure jumps on 1D equatorial slices, and Figure 7 shows the feature sitting at the same radius over ten time slices. That is a decent case for persistence rather than a transient blob. The authors are honest about the central caveat: they state that after about 10,000 tg, the 2D nature of the runs means MRI cannot continuously grow the field. The quantitative shock analysis is averaged over 12,000-15,000 tg, exactly inside that window. So for the alpha=50 cases where the shocks look cleanest, the flow may be effectively laminar, and the structure could be a centrifugal feature imposed by the alpha=50 density confinement rather than proof that shocks survive turbulence. The 3D validation is at alpha=10, not 30 or 50, at 256x80x64, and shows a single oscillating shock near r~37 rg. That is an existence proof at low alpha, but it does not test the parameter branch carrying the headline claim. There is no convergence study, no Rankine-Hugoniot check, no Mach number or turbulence-stress diagnostic. So the paper is not as strong as its abstract advertises; the claim should be softened to evidence from 2.5D GRMHD with partial 3D confirmation at lower inclination. That said, the result is plausible, the code is standard, the setup is clearly described, and the literature is well covered. This deserves a serious referee; a good revision would add an alpha=50 3D run or justify the cost, a convergence test, and a jump-condition check. I would bring it to the reading group; it is a good prompt for discussing what counts as first evidence.","headline":"Plausible first GRMHD standing shocks, but the headline claim rides on 2.5D runs taken after MRI decayed, and the 3D validation skips the exact parameters used for the strongest shocks.","tokens_in":15121,"tokens_out":2167,"would_cite":false,"duration_ms":27738,"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 low-angular-momentum, weakly magnetized accretion onto a rotating black hole produces two quasi-steady standing shocks, the first time standing shocks have been reproduced in full GRMHD simulations.","keywords":["accretion","black holes","GRMHD simulations","standing shocks","low angular momentum","SANE regime","quasi-periodic oscillations"],"falsifier":"Run a high-resolution 3D GRMHD simulation of the alpha = 50, a* = 0.94 case for at least 15,000 tg and check the equatorial radial profiles of pressure and radial velocity: if the two sharp jumps at about 4-5 rg and 40-50 rg do not form or do not persist over several thousand tg, the claim of standing shocks in GRMHD SANE low-angular-momentum flows would be contradicted.","tokens_in":1832,"feed_emoji":"🕳️","tokens_out":2438,"duration_ms":70231,"temperature":0.7,"pith_summary":"The paper aims to show that standing shocks, abrupt jumps in pressure and velocity where infalling gas slows from supersonic to subsonic, are genuine features of magnetized, turbulent black hole accretion, not artifacts of idealized hydrodynamic or semi-analytic treatments. Using 2.5D axisymmetric GRMHD simulations of low-angular-momentum flows in the weakly magnetized SANE regime, the authors report two quasi-steady shocks at roughly 4 to 5 gravitational radii and 40 to 50 gravitational radii when the inflow is strongly concentrated toward the equatorial plane, with the outer shock moving outward as the black hole spin increases. A 3D run at a lower inclination parameter shows an oscillating shock at about 37 gravitational radii, lending support to the 2.5D picture. If correct, this connects several decades of shock theory to state-of-the-art GRMHD simulations and offers a concrete physical mechanism for observed variability such as quasi-periodic oscillations and flares.","feed_headline":"First GRMHD evidence of standing shocks in black hole accretion","feed_subtitle":"New simulations show two quasi-steady shocks in low-angular-momentum flows onto spinning black holes.","key_machinery":"The central control parameter is the inclination parameter $\\alpha$, which concentrates the Fishbone-Moncrief torus density toward the equatorial plane via rho proportional to exp[-($\\alpha$ cos $\\theta$)^2], combined with a reduced angular momentum fraction F = 0.4 and a weak poloidal magnetic field. These choices establish a pressure balance between centrifugal force (amplified by frame dragging near a spinning black hole) and the ram pressure of the infalling gas; when $\\alpha$ is large this balance supports two standing shocks, at roughly 4-5 rg and 40-50 rg, with increasing spin pushing the outer shock outward.","core_discovery":"The authors report the first GRMHD evidence of standing shocks in low-angular-momentum accretion flows. In the weakly magnetized SANE regime, when the initial torus is confined near the equatorial plane (inclination parameter alpha = 50), two standing shocks form at r approximately 4-5 rg and r approximately 40-50 rg, and the outer shock shifts outward for higher black hole spin because frame dragging strengthens centrifugal support. The shock locations remain quasi-steady between 14,000 and 15,000 tg, and a 3D validation run at alpha = 10 shows a quasi-steady, gently oscillating shock at about 37 rg. This demonstrates that such structures are not limited to idealized hydrodynamics but arise naturally in turbulent, weakly magnetized GRMHD flows around a rotating black hole.","pith_inferences":["The 2.5D runs may underestimate how full 3D turbulence disrupts the shocks; a 3D run at alpha = 50 and high spin would be the decisive test of whether the strong two-shock structure survives.","The finding suggests that earlier GRMHD studies that saw no standing shocks may have used initial conditions that were too close to Keplerian or insufficiently concentrated toward the equatorial plane, with alpha as the key control.","If standing shocks require the weakly magnetized SANE regime, shock presence could serve as an observational discriminant between SANE and magnetically arrested (MAD) accretion states.","The connection between shock radius and spin implies that correlated timing and spectral observations of quasi-periodic oscillations might be used to infer both the accretion geometry and the black hole spin in a single source."],"forward_implications":["Standing shocks are not artifacts of 1.5D idealized flows: they can persist in turbulent magnetized GRMHD simulations, so shock-based analytic models gain a firmer footing.","The outer shock location depends on black hole spin through frame dragging, so measuring shock positions in observations could offer a spin diagnostic.","The hot, dense post-shock region acts as a post-shock corona that can up-scatter soft photons, linking these simulations to the hard X-ray tails of black hole X-ray binaries.","Quasi-steady radial oscillation of the shock front (seen in the 3D run) provides a plausible source of low-frequency quasi-periodic oscillations and flares such as those observed in Sgr A* and X-ray binaries.","If confirmed in full 3D at high alpha, the shock structures give a concrete target for general relativistic radiative transfer calculations that would predict direct observational signatures."],"supporting_citations":[{"why":"Supplies the BHAC code that solves the GRMHD equations used in all simulations.","marker":"Porth et al. (2017)"},{"why":"Documents the BHAC code and its GRMHD setup, load-bearing for the numerical method.","marker":"Olivares et al. (2019)"},{"why":"Previous GRMHD study of low-angular-momentum accretion whose setup is closely followed and which did not report standing shocks, defining the baseline the paper extends.","marker":"Dihingia & Mizuno (2024)"},{"why":"Recent GRMHD simulations of low-angular-momentum flow that saw only local density jumps, not standing shocks, providing the direct counterexample the paper addresses.","marker":"Olivares et al. (2023)"},{"why":"Established the analytic theory of standing shocks in low-angular-momentum accretion flows that the simulations are testing.","marker":"Chakrabarti (1989a)"},{"why":"Provides the initial torus solution that is modified by the sub-Keplerian angular momentum fraction and inclination parameter.","marker":"Fishbone & Moncrief (1976)"},{"why":"Identifies the magnetorotational instability that drives the turbulence in the simulated accretion flow.","marker":"Balbus & Hawley (1991)"},{"why":"Defines the SANE regime of weakly magnetized accretion that the simulations are designed to occupy.","marker":"Narayan et al. (2012)"}],"fun_headline_variants":["GRMHD simulations show first standing shocks in black hole accretion","Standing shocks found in GRMHD black hole accretion for first time","GRMHD simulation reveals quasi-steady standing shocks near black hole","Low angular momentum black hole accretion yields standing shocks in GRMHD"],"cache_read_input_tokens":17152,"weakest_assumption_plain":"The main results come from 2.5D axisymmetric simulations that the authors acknowledge cannot sustain MRI growth of magnetic fields after about 10,000 gravitational times, and the single 3D validation was run at the lower inclination parameter alpha = 10 rather than the alpha = 50 cases where the two-shock structure is strongest; if full 3D turbulence at high alpha destroys or prevents the shocks, the central claim fails.","fun_headline_variants_meta":{"raw":{"variants":["GRMHD simulations show first standing shocks in black hole accretion","Standing shocks found in GRMHD black hole accretion for first time","GRMHD simulation reveals quasi-steady standing shocks near black hole","Low angular momentum black hole accretion yields standing shocks in GRMHD"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000996,"raw_usage":{"total_tokens":4189,"prompt_tokens":885,"completion_tokens":3304,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":501,"completion_tokens_details":{"reasoning_tokens":3229}},"tokens_in":501,"tokens_out":3304,"duration_ms":27084,"temperature":1.0,"reasoning_tokens":3229,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:35:30.888302+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a high-resolution 3D GRMHD simulation of the alpha = 50, a* = 0.94 case for at least 15,000 tg and check the equatorial radial profiles of pressure and radial velocity: if the two sharp jumps at about 4-5 rg and 40-50 rg do not form or do not persist over several thousand tg, the claim of standing shocks in GRMHD SANE low-angular-momentum flows would be contradicted.","supporting_citations":[],"review_version":1}