{"id":"800f288e-1fab-46fe-a4bb-778a2e468c5f","arxiv_id":"2504.12762","paper_version":4,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"The simplest star-disk collision model for quasiperiodic eruptions fails to simultaneously match the observed luminosity, duration, and temperature for most sources in the sample.","lead":"This paper tests a popular explanation for quasiperiodic eruptions, in which a star repeatedly hits the black hole's accretion disk and the collision produces a hot flare. It matters because the result would narrow the debate over what causes these repeating X-ray bursts, favoring models where a shredded stream of gas hits the disk instead of a whole star.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (12) is not derivable: eliminating mdot between Eqs. (8) and (10) cancels r*, so the red constraints and per-source 'no allowed R_*' regions are unsupported.","rationale":"Algebraic consistency of the central derivation trumps the astrophysical assumptions. The reader's stated weakest assumption (LQ as a proxy for the collision-radius mdot) is a legitimate secondary concern, but the primary failure is more elementary: the one equation that supposedly converts Lp and t_e into an R_* constraint cannot be derived from the paper's own Eqs. (8) and (10). I verified the functional degeneracy: both equations depend on r* and mdot through the single combination r* mdot^{1/2}, so mdot elimination cancels r*. This invalidates the intersection analysis and the per-source verdicts that rely on it. I therefore agree with REJECT, though for a sharper reason than the LQ assumption. The temperature mismatch and Eq. (13)-plus-tidal arguments still raise serious doubts about the model, so the paper remains useful as an observational compilation, but the headline quantitative test needs correction. The reader's verdict should stand unchanged.","tokens_in":15649,"tokens_out":15790,"duration_ms":144620,"concrete_test":"Run a symbolic or numeric elimination of mdot from Eqs. (8) and (10) as printed. Numeric version: for GSN 069 with M6=1, alpha=0.1, Lp=3e42 erg/s, P1=1, solve Eq. (10) for mdot at r*=0.1 and at r*=10, then insert each mdot into Eq. (8). If the two values of t_e are equal (about 0.6 hr for this choice), r* has dropped out and Eq. (12) cannot be a valid constraint. A positive result confirms that the red curves in Figs. 1-3 should be removed and the per-source overlap analysis redone.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (12) is the load-bearing step for the per-source rejection. The paper claims it follows by inferring mdot from Eq. (10) and inserting into Eq. (8), but the algebra cannot produce a constraint on r*. Eq. (8) gives t_e = 0.18 r* alpha^{1/2} mdot^{1/2} (P1/M6)^{2/3} hr, and Eq. (10) gives Lp = 6.1e41 M6 (r*/P1)^{2/3} mdot^{1/3}. Both equations depend on r* and mdot only through the combination q = r* mdot^{1/2}: t_e is linear in q and Lp is proportional to q^{2/3}. Solving Eq. (10) for mdot and substituting into Eq. (8) leaves t_e = 0.18 alpha^{1/2} (P1/M6)^{2/3} (Lp/(6.1e41 M6))^{3/2} P1 hr, with r* canceled exactly. No rearrangement of the same two equations can yield Eq. (12), r* proportional to Lp^{3/4} t_e^{1/2}. Consequently the red lines in Figs. 1-3 are not constraints, and the statements that RX J1301, eRO-QPE1, ASASSN-14ko, and J0230 have no allowed r* because the red and blue curves fail to intersect are unsupported. The blue Eq. (13), based on LQ, is a legitimate curve, and the temperature mismatch from Eq. (14) is a genuine independent problem, but the headline six-of-eight exclusion is not established. In addition, Eq. (9) is internally inconsistent: combining Eq. (10) with Eq. (8) requires E_i at R_diff to scale as r*^{5/3} mdot^{5/6}, not mdot^{-1/6} as printed. The source-specific conclusions must be re-derived before they can be used.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper tests the shock-cooling emission model for quasiperiodic eruptions (QPEs) and QPE-like sources by comparing observed peak luminosity, eruption duration, quiescent luminosity, and peak temperature with analytic model predictions. It derives two constraints on the stellar radius r* (Eqs. 12 and 13), overlays them with tidal-disruption limits for eight sources, and concludes that six of the eight are excluded, with only eRO-QPE3 and eRO-QPE4 requiring r* of order 1 R_sun. The paper also reports a systematic temperature discrepancy between the model prediction (kT_p ~ 10 eV) and observations (~100 eV).","tokens_in":16010,"tokens_out":22321,"duration_ms":190350,"significance":"The paper addresses an active debate on the origin of QPEs and assembles a useful compilation of radiative properties for eight sources. Its cleanest contribution is the peak-temperature prediction, Eq. (14), which is parameter-free in the sense that it does not depend on r*, mdot, alpha, or eta, and the reported factor-of-ten temperature discrepancy is a genuine falsifiable test. If the r*-constraint analysis were valid, the six-of-eight exclusion would be an important discriminator. However, the central algebraic step leading to Eq. (12) is incorrect, so the headline exclusion claim is not currently supported.","major_comments":[{"comment":"Equation (12) is not derivable from Eqs. (8) and (10). Eliminating mdot between these equations cancels r* exactly. From Eq. (10), mdot^{1/2} = [Lp/(6.1e41 M6)]^{3/2} (r*/P1)^{-1}; inserting this into Eq. (8) gives te = 0.18 alpha_{-1}^{1/2} (Lp/6.1e41)^{3/2} P1^{5/3} M6^{-13/6} hr, with no r* dependence. No rearrangement of the same two equations yields an expression for r*, so the red lines in Figs. 1-3 are not constraints on r*, and the statements that RX J1301, eRO-QPE1, ASASSN-14ko, and J0230 have no allowed r* are unsupported. This is the load-bearing step for the paper's main conclusion.","section":"Section 3.6, Eq. (12)"},{"comment":"Equation (9) is internally inconsistent with the preceding equations. Using E_i0 = Msh v_k^2, V0 proportional to R*^2 h, R_diff = vej te, and te from Eq. (8), the adiabatic factor (V0/V_diff)^{1/3} scales as r*^{-1/3} alpha^{-1/2} mdot^{-1/6} P1^{-1/3} M6^{2/3}; multiplying by E_i0 proportional to r*^2 alpha^{-1} mdot P1^{1/3} M6^{5/3} yields E_i(R_diff) proportional to r*^{5/3} alpha^{-3/2} mdot^{5/6} M6^{7/3}. The printed Eq. (9) has mdot^{-1/6} and M6^{1/3}, so the mdot (and M6, alpha) scalings are wrong. Since Eq. (10) is derived from Eq. (9), the luminosity constraint and the subsequent r* analysis should be re-derived.","section":"Section 3.4, Eq. (9)"},{"comment":"The surface-density expression in Eq. (3) appears inconsistent with the stated viscosity law. With h = 1.5 mdot Rg from Eq. (2), nu = alpha sqrt(GMR) (h/R)^2, and R from Eq. (1), one obtains Sigma_d = Mdot/(3 pi nu) proportional to alpha^{-1} P1 mdot^{-1} M6^{-1}, not the printed proportionality to P1 mdot M6. This changes the scaling of Msh in Eq. (4), and hence of te in Eq. (8), Lp in Eq. (10), and the final r* constraints. The authors should either confirm the printed Eq. (3) with a specific disk model or correct it, and then re-derive the subsequent equations.","section":"Section 3.1, Eq. (3)"}],"minor_comments":[{"comment":"The source name 'ASASSIN14-ko' appears to be a typo; the standard name used throughout the rest of the paper is ASASSN-14ko.","section":"Section 1, text"},{"comment":"The units of E_i(R_diff) are written as erg s^{-1}; since this is an internal energy, the units should be erg.","section":"Section 3.4, Eq. (9)"},{"comment":"The sentence 'the shock-cooling emission model from star-disk collisions can work for RX J1301' is inconsistent with the surrounding argument and with the paper's own conclusion in Section 6; this is likely a typo for 'cannot work'.","section":"Section 4.2, paragraph on RX J1301"},{"comment":"The eight predicted-temperature curves and eight observed points are overplotted in a single panel; separate panels with shared axes would make the comparison much easier to read.","section":"Section 4.4, Figure 4"},{"comment":"The retrograde-orbit constraints are derived using the same framework as Eqs. (8)-(13), so they should be re-examined after the errors in Eqs. (3), (9), and (12) are corrected.","section":"Section 5.1, retrograde-orbit analysis"}],"recommendation":"major_revision","confidential_remarks":"The manuscript contains a clear algebraic error in the derivation of Eq. (12) that invalidates the central exclusion claim, and there are additional inconsistencies in the disk-surface-density and internal-energy scalings. I recommend major revision rather than rejection because the framework is repairable: the analysis can be redone by replacing the flawed r* constraint with a consistency test for te, and by correcting the disk scaling laws. The temperature test and the source compilation are valuable and should be retained. No concerns about citation patterns or novelty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper has one solid, independent result and one load-bearing error. The temperature comparison is real; the per-source radius exclusions are not.\n\nWhat the paper does well: it compiles eight QPE/QPE-like sources and uses L_p, t_e, L_Q, and T_p to test the shock-cooling emission model. The sample table is useful, and the parameter-free T_p prediction is a clean check: the model gives kT_p ~ 10 eV while most sources show ~100 eV. That discrepancy is worth knowing about, even though Linial & Metzger (2023) already flagged it.\n\nThe soft spot is fatal to the main claim. Eq. (12) does not follow from Eqs. (8) and (10). Both equations depend on r_* and mdot through the same combination (effectively r_* mdot^{1/2}); eliminating mdot cancels r_* exactly. So the red lines in Figures 1–3 are not constraints on the stellar radius, and the 'no allowed R_*' regions for RX J1301, eRO-QPE1, ASASSN-14ko, and J0230 are unsupported. Eq. (9) is also inconsistent with the printed Eqs. (8) and (10). The blue L_Q constraint (Eq. 13) is legitimate, but alone it does not yield the paper's six-of-eight exclusion.\n\nSo the novel part of the paper—the per-source rejection—does not hold up. The presentation is clear and the authors engage honestly with the literature, but the derivation error is central not peripheral.\n\nRecommendation: send it to peer review anyway, because the T_p discrepancy is a clean challenge to the model and the data compilation has value. But a referee should ask for a full re-derivation of the radius constraints. As it stands, this is a major-revision-or-reject.","headline":"The temperature comparison is genuinely useful, but the paper's central stellar-radius constraint (Eq. 12) is not derivable, so the per-source exclusions are unsupported.","tokens_in":16616,"tokens_out":8810,"would_cite":false,"duration_ms":79755,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that the simplest star-disk collision explanation for quasiperiodic eruptions fails for six of eight sources, because no stellar radius can match the observed luminosity and duration without the star being tidally…","keywords":["quasiperiodic eruptions","star-disk collisions","shock-cooling emission","accretion disks","tidal disruption","X-ray transients","black hole accretion"],"falsifier":"Measure the quiescent disk accretion rate at the collision radius by a route independent of $L_{\\rm Q}$ (for example from the disk continuum normalization or from X-ray variability and reverberation) and compare it with the value inferred from Eq. (11); a significant disagreement would invalidate the $L_{\\rm Q}$-based radius constraint. Alternatively, a sensitive EUV-to-soft-X-ray spectrum of eRO-QPE3 or eRO-QPE4 during an eruption that shows the peak at about 100 eV rather than the predicted about 10 eV would directly falsify the shock-cooling emission model for the sources where the radius test passes.","tokens_in":15369,"feed_emoji":"🌟","tokens_out":16488,"duration_ms":145854,"temperature":0.7,"pith_summary":"Quasiperiodic eruptions are repeated soft X-ray flares from galactic nuclei, and one popular explanation is that an orbiting star periodically plows through the black hole's accretion disk, producing a shock-cooling fireball that radiates the flare. This paper derives what stellar radius $R_\\star$ such a model needs to match the observed peak luminosity, duration, quiescent luminosity, period, and black-hole mass, and it checks that radius against the requirement that the star survive tidal disruption. Across six QPE sources and two QPE-like sources, the test finds that only eRO-QPE3 and eRO-QPE4 admit a roughly solar-radius star; the others either have no self-consistent radius or require a star so large that it would be torn apart. The model also predicts peak temperatures near 10 eV while the observed values cluster near 100 eV. If correct, the result rules out the simplest star-disk shock-cooling picture for most of the sample and redirects attention to variants such as stream-disk collisions or repeated partial tidal stripping.","feed_headline":"No stellar radius fits six of eight quasiperiodic eruptions","feed_subtitle":"Only eRO-QPE3 and eRO-QPE4 survive the radius test; the model's heat is ten times too low.","key_machinery":"The central machinery is the supernova-like shock-cooling diffusion model for the collision ejecta. A star of radius $R_\\star$ sweeps up a disk mass $M_{\\rm sh}\\simeq 4.3\\times10^{-7}\\,r_\\star^2 P_1\\,\\alpha_{-1}^{-1}\\dot{m}M_6\\,M_\\odot$ (Eq. 4); that material is heated to internal energy $E\\simeq M_{\\rm sh}v_{\\rm K}^2$ and expands adiabatically at $v_{\\rm ej}\\simeq(2/3)v_{\\rm K}$. The eruption duration follows from the photon-diffusion condition $\\kappa\\Sigma=c/v_{\\rm ej}$ (Eq. 8), the peak luminosity from the internal energy remaining at the diffusion radius divided by $t_{\\rm e}$ (Eq. 10), and the quiescent luminosity fixes the accretion rate (Eq. 11). Eliminating $\\dot{m}$ two different ways produces the two stellar-radius constraints (Eqs. 12 and 13), while the tidal-disruption limit for a possibly inflated star (Eq. 16) supplies the survival bound. The same scalings give the predicted color temperature (Eq. 14) used for the final comparison.","core_discovery":"The paper treats each eruption as a radiation-dominated shock-cooling fireball: the star sweeps up a disk mass $M_{\\rm sh}\\propto R_\\star^2\\Sigma_d$, the shocked gas expands at roughly $v_{\\rm ej}\\simeq 2v_{\\rm K}/3$, and the flare peaks when the expanding fireball becomes thin enough for photons to escape. Combining this with the $\\alpha$-disk surface density, the quiescent luminosity $L_{\\rm Q}$ (Eq. 11), and the tidal-disruption radius (Eq. 15) yields two independent constraints on the stellar radius $R_\\star$ (Eqs. 12 and 13) plus an upper bound from tidal survival (Eq. 16). Applied to six QPE sources and two QPE-like sources, the paper finds a self-consistent, non-disrupted stellar radius only for eRO-QPE3 ($R_\\star\\sim1\\,R_\\odot$) and eRO-QPE4 ($R_\\star\\sim2\\,R_\\odot$). For GSN 069 and eRO-QPE2 any allowed radius implies repeated partial tidal stripping; for RX J1301, eRO-QPE1, ASASSN-14ko, and J0230 no radius satisfies both luminosity and duration constraints. The predicted peak temperature from Eq. (14) is $kT_{\\rm p}\\sim10$ eV, one order of magnitude below the observed $\\sim100$ eV.","pith_inferences":["The $L_{\\rm Q}$-based constraint is only as good as the assumption that quiescent X-rays trace the same accretion rate that sets the disk surface density at the collision radius; a spectral or variability decomposition of $L_{\\rm Q}$ would show whether the no-solution sources are truly excluded.","The temperature gap suggests that relaxing the full-thermalization assumption, for example with a Comptonized or non-blackbody photosphere, could raise $T_{\\rm p}$ without changing $L_{\\rm p}$ or $t_{\\rm e}$, potentially reviving the model for the two surviving sources.","The same two-constraint overlap test can be applied immediately to any newly discovered QPE: with only $P$, $t_{\\rm e}$, $L_{\\rm p}$, and $L_{\\rm Q}$, one can check whether any stellar radius exists before undertaking timing analysis.","If the stream-disk variant is correct, eruptions should carry a signature of the stream's orbital-energy spread, a collision that lasts several hours rather than being nearly instantaneous, which is a testable discriminator between the two scenarios."],"forward_implications":["Only eRO-QPE3 and eRO-QPE4 can be powered by a whole, roughly solar-mass star colliding with the disk; the other six sources require a different mechanism or repeated partial tidal stripping.","For GSN 069 and eRO-QPE2, any acceptable star-disk collision would be dominated by accretion of tidally stripped stellar material rather than by the shock-cooling emission itself.","The two long-period QPE-like sources (ASASSN-14ko and J0230) are incompatible with the model, because the required radii of $10^3$–$10^4\\,R_\\odot$ are unphysical.","A systematic factor-of-ten shortfall in predicted peak temperature means the fully thermalized blackbody assumption is inadequate even if the radius constraints are relaxed.","If the collider is instead a debris stream shed by the star, the same framework (with $R_\\star^2$ replaced by the stream cross section and the tidal limit dropped) can accommodate four of the six QPE sources, though not the two long-period QPE-like events."],"supporting_citations":[{"why":"Introduces the star-disk collision scenario that the shock-cooling emission model tested here builds on.","marker":"Dai et al. 2010"},{"why":"Supplies the analytic scalings for eruption duration, peak luminosity, and temperature that Eqs. (8), (10), and (14) are checked against.","marker":"Linial & Metzger 2023"},{"why":"Provides the $\\alpha$-disk surface density used to compute the shocked mass per collision (Eq. 3).","marker":"Shakura & Sunyaev 1973"},{"why":"Gives the photon-diffusion criterion $\\kappa\\Sigma=c/v_{\\rm ej}$ that sets the eruption duration (Eq. 8).","marker":"Arnett 1980"},{"why":"Provides the tidal-disruption radius used in the survival constraint (Eq. 15).","marker":"Rees 1988"},{"why":"Supplies the GSN 069 discovery data and the $L_{\\rm p}$, $t_{\\rm e}$, $L_{\\rm Q}$, and $kT_{\\rm p}$ values used in the sample.","marker":"Miniutti et al. 2019"},{"why":"Supplies the eRO-QPE1 and eRO-QPE2 eruption properties used in the stellar-radius test.","marker":"Arcodia et al. 2021"},{"why":"Supplies the eRO-QPE3 and eRO-QPE4 eruption properties and black-hole mass estimates.","marker":"Arcodia et al. 2024b"},{"why":"Supplies the ASASSN-14ko flare durations, luminosities, and temperatures used for the long-period QPE-like test.","marker":"Payne et al. 2021"},{"why":"Supplies the J0230 X-ray eruption data used for the second QPE-like source.","marker":"Evans et al. 2023"}],"fun_headline_variants":["Star-disk collision model fails for 6 of 8 QPEs","Only 2 of 8 QPEs survive star-disk collision test","QPE model's predicted temperature is 10x too low","Most QPEs defy star-disk collision model","eRO-QPE3 and eRO-QPE4 are the only fitting QPEs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole test rests on treating the quiescent X-ray luminosity $L_{\\rm Q}$ as a direct measure of the same mass accretion rate that sets the disk surface density at the collision radius, with radiative efficiency near 0.1, so if $L_{\\rm Q}$ has another origin, or the accretion rate at the collision radius differs from the inner disk value, the $L_{\\rm Q}$-based radius constraint collapses.","fun_headline_variants_meta":{"raw":{"variants":["Star-disk collision model fails for 6 of 8 QPEs","Only 2 of 8 QPEs survive star-disk collision test","QPE model's predicted temperature is 10x too low","Most QPEs defy star-disk collision model","eRO-QPE3 and eRO-QPE4 are the only fitting QPEs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000783,"raw_usage":{"total_tokens":3560,"prompt_tokens":1151,"completion_tokens":2409,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":767,"completion_tokens_details":{"reasoning_tokens":2313}},"tokens_in":767,"tokens_out":2409,"duration_ms":17571,"temperature":1.0,"reasoning_tokens":2313,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:26:20.166627+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the quiescent disk accretion rate at the collision radius by a route independent of $L_{\\rm Q}$ (for example from the disk continuum normalization or from X-ray variability and reverberation) and compare it with the value inferred from Eq. (11); a significant disagreement would invalidate the $L_{\\rm Q}$-based radius constraint. Alternatively, a sensitive EUV-to-soft-X-ray spectrum of eRO-QPE3 or eRO-QPE4 during an eruption that shows the peak at about 100 eV rather than the predicted about 10 eV would directly falsify the shock-cooling emission model for the sources where the radius test passes.","supporting_citations":[],"review_version":1}