{"id":"08363839-07b0-4f76-8020-aeca3b5f1cbe","arxiv_id":"2603.00226","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Stellar-mass black holes dragging gas from accretion disks out to their Hill radius can power quasi-periodic eruptions, while stellar impactors produce too-asymmetric and tidally unstable flares.","lead":"This paper uses global 3D simulations to compare stars and stellar-mass black holes crashing through the disks around supermassive black holes as a source of quasi-periodic X-ray eruptions. It concludes that black holes, not stars, can power the observed eruptions once gravitational effects out to the Hill radius are included, and proposes a simple scaling law for the interaction size.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The full 10^44–10^48 erg range rests on Eq. (19)'s ad hoc R_sun^2 normalization; a standard drag calculation may shift energies enough to break the claim.","rationale":"The reader's verdict is CONDITIONAL, and I agree with that overall assessment. My primary concern is not exactly the reader's headline weakest assumption (radiative cooling), but the arbitrary energy normalization in Eq. (19) of the semi-analytical model. This matters because the abstract's headline numerical claim — the 10^44–10^48 erg range — is generated by that model, not directly measured in the simulations. The paper's Eq. (6) is a standard dynamical-friction estimate, but the full-range claim in the Abstract/Summary is attributed to the semi-analytical model, so Eq. (19) is load-bearing. The reader did mention 'arbitrary normalization' in their rationale, so I count this as partial agreement. I also note a separate text-level inconsistency: the Introduction gives observed flare energies as 10^40–10^42 erg, while the Abstract claims the full QPE range is 10^44–10^48 erg. This should be reconciled; I did not make it the primary concern because it may be a typo or a different energy definition (e.g., total budget vs. soft-X-ray radiated energy). The proposed check — recomputing the same budgets with standard dynamical friction — would settle whether the ad hoc normalization controls the numerical claim. If the standard-drag numbers still cover the relevant observed energies, the central claim survives and Eq. (19) is demoted to a heuristic; if they do not, the claimed range is unsupported. In either case the verdict stays CONDITIONAL pending this test, so I recommend UNCHANGED.","tokens_in":19262,"tokens_out":11646,"duration_ms":108977,"concrete_test":"Recompute the §5.3 energy budgets for i=18° and i=90° using the standard Chandrasekhar dynamical-friction power F = 4π(G m_sBH)^2 rho_0 lnΛ / v_rel^2, integrated over the disk-crossing path ΔL = 2H/cos(i/2), with the same H=1 R_sun, Σ=10^6 g cm^-2, and m_sBH=100 M_sun. Compare the resulting E_sBH with the quoted 3.52×10^46 erg and 8.81×10^44 erg. If either value changes by more than an order of magnitude, or if the low-to-high-inclination ratio changes by more than a factor of 10, then the claimed 10^44–10^48 erg range is an artifact of the R_sun^2 normalization in Eq. (19) rather than a robust prediction of the model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim — that sBH-disk collisions 'can power the full QPE energy range (10^44–10^48 erg)' — is generated by the semi-analytical model in §5.3, not measured directly in the simulations. The load-bearing step is Eq. (19): dE_sBH/dt = (G m_sBH / R_sun^2) rho_0 v_rel (4/3) pi R_eff^3. This is not derived from the simulations or from standard gravitational drag. The paper adopts R_sun as a 'characteristic normalization scale' because the physical radii lie between ~0.03 and ~10 R_sun. That is a dimensional fudge: an equally plausible normalization (R_B, R_H, or the geometric mean) changes E_sBH by factors of order (R_sun/R_eff)^2 or (R_eff/R_sun)^2, which can reach 10^2–10^4 across the explored inclination range. Since the claimed QPE range spans only four orders of magnitude, the arbitrary normalization — not the hydrodynamics — is what decides whether the full range is covered. The empirical R_eff formula, Eq. (10), is also calibrated to only four runs, but even if R_eff is correct, Eq. (19)'s missing derivation leaves the absolute energy scale unanchored. The absence of radiative transfer (§5.4) is a separate limitation; the immediate numerical claim fails if Eq. (19) is replaced by a standard dynamical-friction estimate.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses global 3D meshless finite-mass (MFM) simulations with GIZMO to compare stellar-mass black hole (sBH) and stellar impactors crossing TDE accretion disks around an SMBH, including the SMBH potential. It reports that stellar impacts produce strongly asymmetric, forward-dominated ejecta and face tidal stability problems, while 100 Msun sBH impacts produce nearly symmetric ejecta with energies that grow at low inclination. The paper introduces an ad hoc effective interaction radius R_eff ~ 0.5 R_B^{1/3} R_H^{2/3}, uses it in a semi-analytical model, and claims that sBH-disk collisions can power the full QPE energy range (10^44-10^48 erg) without invoking intermediate-mass black holes.","tokens_in":19712,"tokens_out":8103,"duration_ms":75450,"significance":"If the quantitative central claim holds, this is a potentially important contribution: the global setup including the SMBH potential is a genuine step beyond local box simulations, and the qualitative distinction between asymmetric stellar ejecta and quasi-symmetric sBH ejecta is directly relevant to the QPE debate. The paper is also commendably transparent that Eq. (10) is ad hoc and that radiative transfer is not included. However, the central quantitative claim rests on a dimensional normalization in Eq. (19) and on an empirical radius calibrated to only four sBH runs; the confirmation in Section 5.3 is therefore partly circular. The strengths are the qualitative simulation findings and the honest discussion of limitations; the weakness is the unanchored absolute energy scale in the semi-analytical model.","major_comments":[{"comment":"The paper states in the Introduction that QPE flare energies are 10^40-10^42 erg, but the abstract, §5.4, and §6 base the central claim on a 'full QPE energy range (10^44-10^48 erg)'. These differ by four orders of magnitude. Since the claimed viability of sBH impacts depends on covering the observed energy range, this inconsistency must be resolved. If the two numbers refer to different quantities (e.g., emitted soft X-ray energy versus total energy budget), define them explicitly.","section":"§1 vs §5.4/§6"},{"comment":"Eq. (19) is not derived from the simulations or from a standard gravitational-drag calculation. The R_sun^2 normalization is introduced as a 'characteristic normalization scale' because the physical scales fall between 0.03 and 10 R_sun. That is a dimensional choice, not a physical result: replacing R_sun^2 by R_eff^2, R_B R_H, or an integrated drag expression changes E_sBH by factors of order (R_eff/R_sun)^2, which across the explored parameter range can shift the claimed 10^44-10^48 erg interval substantially. Because this equation, rather than the simulation data, generates the central quantitative claim, the absolute energy scale is unanchored. Please derive the rate from a controlled calculation (e.g., a calibrated version of Eq. 6) or present the simulation kinetic energies with uncertainties as the primary evidence.","section":"§5.3, Eq. (19)"},{"comment":"The empirical R_eff relation is explicitly ad hoc and is calibrated to only four sBH runs: two inclinations (π/10 and π/2) and three surface densities, all with m_sBH=100 Msun, M_SMBH=10^6 Msun, and H=1 Rsun. The same relation is then used in §5.3 to 'confirm' the simulation trends and to produce the full 10^44-10^48 erg range, so the confirmation is partly circular. No error bars or convergence tests are reported for the ejecta masses and energies in Table 1. Independent runs with different masses, scale heights, and inclinations, and a clear separation between calibration and validation, are needed before extrapolation to the observed QPE population.","section":"§5.2, Eq. (10); Table 1"},{"comment":"The simulations use a polytropic gamma=4/3 disk with no radiative transfer, so shocked gas retains all thermal energy over the hour-long runs. The measured ejecta masses, expansion speeds, and hence the calibration of R_eff may change if radiative cooling is efficient in real TDE disks. The paper acknowledges this in §5.4 but does not quantify the effect on the energy budget. Since the central claim depends on these simulation-derived quantities, the magnitude of the cooling uncertainty should be estimated or bounded.","section":"§3, §5.4"},{"comment":"The simulation energies in Table 1 are estimated as (1/2) M_sh v_rel^2, whereas the semi-analytical model in §5.3 uses the gravitational energy rate in Eq. (19). The 'qualitative reproduction' of the simulation trend therefore compares two different energy definitions. The paper should state which energy channel is meant by 'burst energy' and should show that the semi-analytical model and the simulations measure the same quantity before using the model to extend the simulation results to the full QPE range.","section":"Table 1 vs §5.3"}],"minor_comments":[{"comment":"In §4.2, references to 'Fig. 2(a)' and 'Fig. 2(b)' appear to refer to Fig. 5(a) and Fig. 5(b); please correct the cross-references.","section":"Figures"},{"comment":"There are duplicated or incomplete entries (e.g., Chakraborty et al. 2025 appears twice; multiple Zhou et al. 2025 entries with abbreviated author lists) and some line breaks in the reference list should be cleaned up.","section":"References"},{"comment":"The sentence 'flare energies of 10^40-10^42 erg' should be cross-checked with the later 10^44-10^48 claim; if they denote different quantities, define them explicitly.","section":"§1"},{"comment":"Eq. (15) uses v_k sin i for the residence time while Eq. (13) uses 2 v_k sin(i/2) for the relative velocity; the near-equality at low i is coincidental and the notation should state which velocity enters each step.","section":"§5.3"}],"recommendation":"major_revision","confidential_remarks":"The qualitative simulation results are worth publishing, but the central quantitative claim is built on an ad hoc normalization and a circular calibration, and the observed energy range is stated inconsistently. These issues can be fixed by reframing the claim, deriving or calibrating the energy rate, and adding caveats; I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, this paper is worth reading for the simulations, but don't take the energy-range headline at face value. The genuinely new thing is the global 3D setup: unlike the local patches in Yao et al. and Huang et al., these runs include the SMBH potential, so the Hill sphere is finite and the Bondi-to-Hill transition is actually resolved. The star-vs-sBH comparison is the strongest part. The single star run shows the forward/backward asymmetry cleanly, with forward ejecta dominating by over an order of magnitude, and the sBH runs show much more symmetric ejecta. That qualitative contrast is a real step forward and directly challenges the Linial & Metzger assumption that stellar radii dominate over Bondi radii.\n\nThe soft spot is exactly where you'd expect. Eq. (10), the effective radius, is \"ad hoc\" by the authors' own admission, calibrated to four sBH runs. And Eq. (19) uses R_sun^2 as a normalization because the physical scales happen to be in that range; that is a dimensional fudge, not a derived result. Since the claimed QPE range spans 10^44-10^48 erg, and changing the normalization to something like R_B^2 or R_H^2 shifts the absolute energy by orders of magnitude, the \"full range\" claim is not established by this paper. The semi-analytical model in Sec 5.3 is a fitting exercise, not an independent confirmation. If the authors replaced Eq. (19) with a standard dynamical-friction estimate, the energies could shift enough to break the headline.\n\nI also note the acknowledged limitations: no radiative transfer, fixed H=1 R_sun, polytropic gamma=4/3, no convergence tests, no code/data release. Those are real but not disqualifying for a first global study; they mostly affect the quantitative calibration. The paper is honest about these in Sec 5.4, which I credit. The qualitative conclusions—stellar impacts produce one-sided bursts and tidal stability is a problem, sBH impacts are more symmetric and can be energetically viable at low inclinations—are supported and are likely to survive more detailed modeling.\n\nWho is this for? Anyone working on QPE models, EMRI-disk interactions, or TDE disk physics. It deserves a serious referee. The referee should push for a proper derivation of R_eff or a physically motivated energy deposition formula, error bars from resolution studies, and a radiative-transfer treatment before the \"full QPE energy range\" claim is accepted. I would accept it for review, with major revision expected. Bring it to reading group if you want a good discussion of where semi-analytic calibration crosses the line into fitting.","headline":"Genuinely novel global simulations that make the star/sBH asymmetry case well, but the 'full QPE energy range' claim rests on an ad hoc normalization that the referee should push on.","tokens_in":20183,"tokens_out":2284,"would_cite":true,"duration_ms":23444,"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":"Stellar-mass black holes, not stars, explain quasi-periodic eruptions.","keywords":["quasi-periodic eruptions","stellar-mass black holes","accretion disks","galactic nuclei","hydrodynamic simulations","Bondi-Hoyle accretion","tidal disruption events","X-ray transients"],"falsifier":"Re-run the sBH-disk collision with a radiation-hydrodynamic scheme that allows the shocked gas to cool and radiate during the ~1 hour simulated interaction. If including cooling reduces the ejecta mass or expansion speed by a large factor, the claimed energy range of 10^44–10^48 erg and the R_eff calibration would not survive. Alternatively, a QPE source with a well-measured disk surface density whose flare energy requires an impactor above ~100 solar masses at any inclination would falsify the paper's central claim that no intermediate-mass black hole is needed.","tokens_in":19166,"feed_emoji":"🕳️","tokens_out":6621,"duration_ms":66773,"temperature":0.7,"pith_summary":"Quasi-periodic eruptions (QPEs) are recurring soft X-ray flares in galactic nuclei whose origin is debated. This paper uses global 3D hydrodynamic simulations to compare two proposed impactors crossing a tidal-disruption accretion disk: a solar-mass star and a ~100-solar-mass black hole. It finds that stellar impacts produce strongly one-sided ejecta—the forward outburst dominates by over an order of magnitude—so only one eruption per orbit should be visible, and several observed sources sit dangerously close to the tidal-disruption radius. A stellar-mass black hole, by contrast, gravitationally focuses and heats disk gas across a region extending from its Bondi radius to its Hill radius, producing nearly symmetric two-sided ejecta and, at low orbital inclination, enough energy to match the whole observed range of 10^44–10^48 erg. The paper's central proposal is an effective interaction radius R_eff = 0.5 (R_B R_H^2)^{1/3} that replaces the Bondi-only estimate and makes stellar-mass black holes viable without invoking intermediate-mass black holes.","feed_headline":"Stellar-mass black holes can power all observed QPE energies","feed_subtitle":"Gravitational focusing between the Bondi and Hill radii supplies the missing energy, with no intermediate-mass black hole needed.","key_machinery":"The central object is an ad hoc effective interaction radius R_eff = 0.5 R_B^{1/3} R_H^{2/3}, where R_B is the Bondi radius (the scale at which the intruder's gravity overcomes the gas's thermal and kinetic energy) and R_H is the Hill radius (the scale at which the supermassive black hole's tidal gravity strips gas away). This radius quantifies the region over which a stellar-mass black hole gravitationally focuses, heats, and ejects disk gas during a crossing. It is validated and calibrated by global 3D meshless finite-mass hydrodynamic simulations that resolve the Bondi-to-Hill scales and measure the forward and backward ejecta masses and energies.","core_discovery":"The paper establishes that the standard way of estimating the energy released when a stellar-mass black hole crosses an accretion disk—taking its Bondi radius as the shocked region—underestimates the interaction by a large margin. In global simulations that include the full gravitational potential of the central supermassive black hole, the effective radius over which the intruder gravitationally heats disk gas follows R_eff ≈ 0.5 (R_B^{1/3} R_H^{2/3}), an interpolation between the Bondi and Hill radii with the previously ignored Hill radius dominating. The measured ejecta masses and energies match this scaling, and the paper shows that an sBH of roughly 50–100 solar masses crossing a disk w","pith_inferences":["If the sBH picture is correct, QPEs become an electromagnetic tracer of stellar-mass black holes at galactic centers, potentially the same population that space-based gravitational-wave observatories aim to detect—a connection the paper raises but does not develop.","The R_eff scaling is an empirical fit from adiabatic runs; a testable extension is that faster radiative cooling of the shocked gas should shrink the effective radius and steepen the inclination dependence of flare energy.","The same gravitational-drag enhancement across Bondi-to-Hill scales may apply to other disk-embedded bodies, such as compact objects in AGN disks or planets in protoplanetary disks, where energy-budget estimates currently use only the Bondi or physical radius.","Because the paper fixes the vertical scale height and uses an adiabatic equation of state, a natural next step is to vary H and the cooling time; the claimed 10^44–10^48 erg range could then be refined into a sharper, source-specific prediction."],"forward_implications":["Stellar-mass black holes of roughly 50–100 solar masses on low-inclination orbits become a viable explanation for QPEs such as GSN 069, without requiring intermediate-mass black holes.","The observed alternating long–short recurrence intervals and strong–weak flare contrasts can arise naturally from two mildly asymmetric impacts per orbit of an sBH, whereas the star model predicts only one detectable burst per orbit.","The sBH model predicts eruption durations spanning roughly an order of magnitude through its inclination-dependent interaction radius, matching the observed diversity of QPE durations.","Several short-period QPE sources lie within twice the tidal-disruption radius for a solar-mass star, so the star model faces a severe stability problem; the sBH model has no corresponding lifetime limit.","The formula R_eff ≈ 0.5 (R_B R_H^2)^{1/3}, if correct, boosts the predicted sBH-disk energy budget by orders of magnitude at low inclinations relative to the Bondi-only estimates used in earlier work."],"fun_headline_variants":["Hill radius boosts black hole impact energy for QPEs","Stellar-mass black holes: QPEs without intermediate masses","Bondi-Hill scaling powers QPEs from stellar black holes","Gravitational drag: stellar BHs explain QPE energy budget"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire energy budget rests on treating the disk as an adiabatic polytropic gas with a fixed vertical scale height of one solar radius and no radiative cooling: if the shocked gas cools quickly in real tidal-disruption disks, the heated ejecta would be less massive and slower, and the calibrated effective radius—and with it the computed energies—would shrink.","fun_headline_variants_meta":{"raw":{"variants":["Hill radius boosts black hole impact energy for QPEs","Stellar-mass black holes: QPEs without intermediate masses","Bondi-Hill scaling powers QPEs from stellar black holes","Gravitational drag: stellar BHs explain QPE energy budget"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000721,"raw_usage":{"total_tokens":3131,"prompt_tokens":858,"completion_tokens":2273,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":2201}},"tokens_in":602,"tokens_out":2273,"duration_ms":15725,"temperature":1.0,"reasoning_tokens":2201,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T19:59:33.707605+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the sBH-disk collision with a radiation-hydrodynamic scheme that allows the shocked gas to cool and radiate during the ~1 hour simulated interaction. If including cooling reduces the ejecta mass or expansion speed by a large factor, the claimed energy range of 10^44–10^48 erg and the R_eff calibration would not survive. Alternatively, a QPE source with a well-measured disk surface density whose flare energy requires an impactor above ~100 solar masses at any inclination would falsify the paper's central claim that no intermediate-mass black hole is needed.","supporting_citations":[],"review_version":1}