{"id":"8d03fad2-40cf-4747-82f2-b0d47da55adf","arxiv_id":"2608.05534","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"No in-phase timing modulation is found in GSN 069's quasi-periodic eruptions, ruling out a companion supermassive black hole with projected light-travel amplitude A0 >= 254 M at 95% credibility.","lead":"This paper analyzes X-ray flare timing from the galaxy GSN 069 and finds no signal of a supermassive black hole binary, contrary to a recent claim. It also maps where a companion black hole would have torn apart the star powering the flares, tightening the excluded zone.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central null result is only as secure as an unpublished cycle-number assignment; a one-cycle slip in the O-C branches can create the exact in-phase signal Eq. (13) excludes.","rationale":"The paper makes a clear, falsifiable null claim. The timing model itself is straightforward: the light-travel-time modulation is a sinusoid added to arrival times, and the Bayesian fit is internally consistent. Credit is due for the ZLK survival diagnostic: the TDE boundary is reported as nearly insensitive to initial eccentricity, and the SA equations are given in Appendix B. However, none of that bears on the decisive data-reduction step. The prior for A0 is broad, but that is not the issue; the issue is that the input O-C series is not independently reproducible. The manuscript signals this itself: footnote 1 concedes that the competing detection is 'likely' a false alarm and defers to a same-group note (Zhou et al. 2026) rather than presenting the cycle-number validation here. That is a missing support explicitly located in the text. The reader's weakest-assumption diagnosis (cycle numbering) matches mine, and their CONDITIONAL verdict is the right one. I would not change it: the scientific claim is plausible but not yet settled until the preprocessing is either published or shown to be robust to the ±1-cycle ambiguity.","tokens_in":14835,"tokens_out":6180,"duration_ms":63925,"concrete_test":"Release the arrival times with cycle numbers, then recompute the fit under all alternate assignments obtained by shifting one contiguous segment of the even/odd branches by ±1 cycle—the class of mismatches that produces spurious in-phase structure. The test fails if any self-consistent assignment yields log B > -1 or a posterior for A0 inconsistent with zero at 90% credibility. Ideally, an independent group should blind-reduce the original X-ray light curve to O-C points without consulting Zhou et al. (2026), then re-run the same model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Eq. 13, A0 < 254 M•, log B = -3.31) is a statement about the O-C timing residuals of GSN 069. The rejection of the competing in-phase detection (Miniutti et al. 2025) is delegated to footnote 1 and to the same-group preprint Zhou et al. (2026), whose analysis is not reproduced or summarized. In an alternating long/short QPE light curve, assigning integer cycle numbers is a discrete preprocessing step: shifting one branch by one cycle changes which intervals are labeled even/odd and can convert the genuine anti-phase apsidal-precession signal into an apparent in-phase component, or erase a real one. No arrival-time table, cycle-number list, or code is provided, so the reader cannot tell whether Eq. (13) is an artifact of this assignment. The Bayes factor compares two hypotheses only after this data-dependent preprocessing, so the 'strong support for the base hypothesis' is conditional on an assumption that is not part of either model. If an independent re-reduction permits a different but self-consistent cycle assignment under which the in-phase amplitude is nonzero, the null claim is overturned. The survival maps in Sec. 4 are a separate diagnostic and do not rescue Eq. (13).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes two ways to search for a wide SMBH companion in QPE sources. Section 3 adds the light-travel-time modulation from the host SMBH's motion about the SMBHB center of mass to the QPE timing model of Zhou et al. (2025b). For GSN 069 the authors report A0 < 254 M• at 95% credibility and log B = -3.31 favoring the base hypothesis, hence no SMBHB. Section 4 computes ZLK survival maps for a solar-type stellar EMRI and shows that a TDE boundary in the (i_in,ini, eps_GR) plane is insensitive to initial eccentricity. The paper concludes that QPE timing and stellar survival are complementary probes of hidden SMBH companions.","tokens_in":15133,"tokens_out":8267,"duration_ms":73550,"significance":"If the null result is robust, the paper makes a useful contribution by quantifying the sensitivity of QPE timing to SMBH companions and by providing a survival-based exclusion map. The derivation of Eq. (8) is clean and parameter-free for a circular outer orbit, and the posterior for A0 is internally consistent with zero. The SA evolution equations in Appendix B are given in detail. However, the headline constraint is not yet reproducible: it depends on a discrete cycle-number assignment dismissed in a one-sentence footnote citing an unpublished same-group note, and no arrival-time table or code is provided. The significance of the direct probe is therefore conditional on an unverifiable preprocessing step.","major_comments":[{"comment":"The central null result is conditional on a cycle-number assignment that is neither tabulated nor reproduced. The manuscript dismisses the competing in-phase detection by Miniutti et al. (2025) in footnote 1 by citing Zhou et al. (2026), an unpublished same-group note. In an alternating long/short QPE sequence, shifting one branch by one cycle changes which intervals are labelled even/odd and can create or destroy exactly the in-phase component that Eq. (13) constrains. Because no arrival-time table, cycle-number list, or analysis code is provided, a reader cannot determine whether A0 < 254 M• and log B = -3.31 are artifacts of the preprocessing. The authors should either include the cycle-number table and demonstrate that the result is stable under all self-consistent integer assignments, or soften the claim to be conditional on the adopted assignment.","section":"Sec. 3, footnote 1, Eqs. (13)-(14)"},{"comment":"The reported log Bayes factor is not accompanied by the likelihood model, evidence integrals, or a definition of the log base (natural vs base 10), and it is computed under the same assumed cycle numbering. Given that the competing detection is dismissed solely by the cycle-number assignment, Eq. (14) cannot be interpreted as a model-independent preference for the base hypothesis. Please specify the exact likelihood, prior volume, and numerical evidence (with uncertainties) and provide the data products needed to recompute it.","section":"Sec. 3, Eq. (14)"},{"comment":"The indirect survival constraint is not translated into the astrophysical parameter space of the putative SMBHB in GSN 069. The scan is presented in (i_in,ini, eps_GR) after marginalizing over a grid of M_e and a_out, but the paper does not state which (M_e, a_out, inclination) combinations are excluded for the observed source, nor does it connect the assumed a_in=300 M• and M•=4e5 M_sun to the posterior constraints from Section 3. The statement that 'the region enclosed by the TDE contour is excluded' holds only under the unproven assumption that the SMO is a surviving solar-type star; without that assumption, the map does not constrain the companion. Please add a projection of the excluded region onto (M_e, a_out) and a discussion of the stellar-survival assumption.","section":"Sec. 4, Fig. 3"}],"minor_comments":[{"comment":"Please fix missing spaces in compound terms such as 'in-phasemodulation' and 'anti-phasemodulation'.","section":"Throughout"},{"comment":"Define the convention for labeling even/odd eruptions and which recurrence interval belongs to each branch; the in-phase/anti-phase classification depends on this convention.","section":"Fig. 1"},{"comment":"State explicitly whether the Bayes factor is a natural logarithm and describe the numerical evidence estimate and its uncertainty.","section":"Eq. (14)"},{"comment":"Define the parameters a, q_r,ini, q_z,ini, q_phi,ini, theta_min, and sigma_sys in the main text or in a parameter table; the corner plot is not self-contained.","section":"Fig. 4"},{"comment":"Add a data availability statement specifying where the GSN 069 arrival times, cycle-number list, and analysis code are available; without these, the central result cannot be checked.","section":"Data availability"},{"comment":"Please standardize author names in citations: 'Sniegowska' has a stray leading apostrophe and 'Huang Xiaoshan' appears in text while the reference list uses 'Huang, X.'.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The manuscript's central claim rests on footnote 1 and an unpublished same-group note (Zhou et al. 2026), which is a citation-pattern concern. I recommend requesting that the authors either include the cycle-number assignment and arrival times in the paper or present a stability analysis over all self-consistent assignments before publication. The paper is otherwise within scope for astro-ph.HE, and the direct-probe formalism is worth preserving."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper's core claim—no SMBHB-induced in-phase timing modulation in GSN 069, A0 < 254 M•—is a useful null result, but it is only as secure as the eruption cycle-number assignment, and that assignment is not shown. The authors dismiss the competing in-phase detection (Miniutti et al. 2025) in a footnote, citing their own unpublished note (Zhou et al. 2026). Without the arrival-time table, the cycle-number list, or code, a referee cannot check whether a one-cycle slip in labeling even/odd eruptions creates or erases the in-phase signal. That is not a manufactured flaw; in an alternating long/short light curve, the even/odd labels are a discrete preprocessing step, and the model comparison is conditional on them. The paper should either include the cycle-number assignment, show the O-C residuals, or fold the competing claim into the analysis rather than footnote it away.\n\nWhat is genuinely new and solid: the separation of the SMBHB light-travel-time effect as an in-phase branch, the Bayesian exclusion for GSN 069, and the ZLK survival maps using single-averaged octupole plus 1PN terms. The survival diagnostic is independent of the cycle-number issue. The TDE boundary is nearly unchanged between e_in,ini = 10^-2 and 10^-5, which is a nice robustness check. The tidal apsidal precession comparison in Appendix C is also careful—it shows tidal effects are negligible at the relevant eccentricities.\n\nThe soft spots beyond the cycle-number issue are minor. The survival map uses a fiducial a_in = 300 M• and M• = 4e5 M⊙, but these are motivated by GSN 069. The direct probe's priors are the same as Zhou et al. 2025b, and the posteriors look internally consistent. The Bayes factor log B = -3.31 is reported with uncertainty, which is good practice.\n\nWho is this for? Anyone working on QPE timing models, SMBHB searches in galactic nuclei, or hierarchical triple dynamics applied to EMRI/TDE systems. It deserves a serious referee: the methodology is new, the null is testable, and the survival maps are a standalone contribution. The right outcome is likely 'revise' rather than 'reject'—add the cycle-number details, provide data/code or at least the O-C residuals, and engage directly with the Miniutti et al. detection instead of delegating to an unpublished note.\n\nMy bottom line: send it to review, but treat the A0 < 254 M• exclusion as provisional until the cycle-number assignment is independently checkable.","headline":"A solid null-result paper with a genuinely new survival diagnostic; the central timing exclusion depends on an unpublished cycle-number assignment that should be made public before the claim is treated as settled.","tokens_in":15735,"tokens_out":2905,"would_cite":true,"duration_ms":26712,"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":"Timing of GSN 069's X-ray bursts rules out a black hole binary.","keywords":["quasi-periodic eruptions","GSN 069","supermassive black hole binary","extreme-mass-ratio inspiral","light-travel-time timing","von Zeipel-Lidov-Kozai mechanism","X-ray timing","accretion disk"],"falsifier":"Re-fit the GSN 069 eruption timings with an independently determined cycle-number assignment, or with a systematic shift of the even/odd labels; if a significant in-phase modulation appears (for example $\\log B>0$ favoring the binary model), the no-companion conclusion is an artifact of the numbering. An independent detection of a sub-parsec binary in GSN 069, for instance through radio imaging or gravitational-wave observations, would also falsify the exclusion.","tokens_in":14602,"feed_emoji":"🕳️","tokens_out":9339,"duration_ms":73142,"temperature":0.7,"pith_summary":"This paper asks whether the quasi-periodic X-ray eruptions from the galactic nucleus GSN 069 carry any timing imprint of a supermassive black hole binary. In the standard picture the eruptions are produced by a stellar-mass object repeatedly crossing an accretion disk around a central black hole, so the eruption clock traces the inner orbit. The authors derive that a distant companion black hole would make the host wobble around the binary's center of mass, adding a light-travel-time delay that appears as an in-phase modulation of even and odd eruption intervals. Fitting this binary-modulation model to GSN 069's observed timings, they find the modulation amplitude is consistent with zero, at 95% credibility $A_0<254\\,M_\\bullet$, and the log Bayes factor $\\log B = -3.31\\pm0.25$ strongly favors the single-black-hole model over the companion model. They also show that a surviving stellar orbiter would be tidally disrupted for a range of companion configurations, giving an independent constraint. If correct, GSN 069 shows no evidence for a supermassive black hole binary at the sensitivity of its eruption timing.","feed_headline":"Timing of GSN 069's X-ray bursts rules out a black hole binary","feed_subtitle":"No binary wobble appears in the eruption clock, and star survival independently shrinks the allowed companion space.","key_machinery":"The load-bearing object is the extended QPE timing model: a Bayesian fit to the observed eruption arrival times in which the base hypothesis (single SMBH, apsidal precession producing anti-phase even/odd branches) is compared with a binary-modulation hypothesis that adds the term $\\delta t_{\\rm binary} = -A_0\\cos(\\omega t + \\Delta\\Phi_0)$ from the host's motion about the binary center of mass. The three new parameters $A_0$, $\\omega$, and $\\Delta\\Phi_0$ carry the companion signature, and the Bayes factor between the two models is the quantitative verdict. The secondary machinery is the hierarchical-triple evolution of Liu and Lai (single-averaged, octupole order, 1PN precession), used to map the maximum inner eccentricity over initial inclination and the relativistic-precession parameter $\\varepsilon_{\\rm GR}$; the boundary $a_{\\rm in}(1-e_{\\rm in,max})=r_{\\rm TDE}$ marks where a surviving star would be tidally disrupted.","core_discovery":"On the paper's own terms, the central claim is a null result: the QPE timing of GSN 069 contains no detectable in-phase modulation of the even and odd recurrence branches, the signature an external supermassive black hole would imprint through the host's motion about the binary center of mass. The posteriors place the modulation amplitude at $A_0<254\\,M_\\bullet$ (95% credibility), and the log Bayes factor $\\log B = -3.31\\pm0.25$ is read as strong support for the base hypothesis with no companion. The same fit recovers the familiar anti-phase apsidal-precession modulation and a central mass $\\log_{10}(M_\\bullet/M_\\odot)=5.6\\pm0.1$, consistent with earlier measurements. For the indirect probe, survival maps computed with single-averaged octupole equations and first-post-Newtonian precession show that a companion in part of the excluded parameter space would drive a solar-type stellar orbiter into tidal disruption, while the TDE boundary itself is nearly insensitive to the assumed initial eccentricity.","pith_inferences":["The authors' own caveat about cycle-number assignment is the main internal risk: if an independent count of eruption cycles changes the even/odd labeling, the in-phase amplitude could shift, so the null result is only as secure as that labeling.","The two diagnostics probe different regions of companion parameter space, so combining them for a single source brackets the allowed companion mass and separation more tightly than either alone.","A longer timing baseline for GSN 069, or application to sources with more eruptions, could push the $A_0$ upper limit well below $254\\,M_\\bullet$ and either reveal a weak signal or sharpen the exclusion.","If any future QPE source shows a genuine in-phase modulation, the same framework would yield a direct measurement of companion mass and separation from $A_0$ and $\\omega$."],"forward_implications":["GSN 069's eruption timing does not support the claimed in-phase binary signature, so the supermassive-black-hole-binary interpretation of this source is disfavored unless the underlying cycle-number assignment is revised.","The same light-travel-time diagnostic can be applied to other QPE sources with long timing baselines, turning a null-result pipeline into a search for hidden black hole companions in nearby galactic nuclei.","The stellar-survival map provides an independent exclusion region whose tidal-disruption boundary is robust to the assumed initial inner eccentricity, so surviving QPE orbiters can constrain companion mass and separation.","The recovered central mass of about $4\\times10^5\\,M_\\odot$ under the binary model confirms that QPE timing remains a viable dynamical measurement of the central black hole even when a companion is allowed."],"supporting_citations":[{"why":"Provides the GSN 069 discovery data and the alternating long-short recurrence pattern that the timing model fits; the empirical basis of the analysis.","marker":"Miniutti et al. 2019"},{"why":"Supplies the base QPE timing model, parameter priors, and the Bayesian framework that the binary-modulation hypothesis extends.","marker":"Zhou et al. 2025b"},{"why":"Argues that previously claimed in-phase detections are false alarms from mismatched cycle numbering, the key assumption behind the null result.","marker":"Zhou et al. 2026"},{"why":"The earlier claim of an in-phase modulation that the paper tests and rejects; defines the positive signal the data are checked against.","marker":"Miniutti, G. et al. 2025"},{"why":"Gives the single-averaged hierarchical-triple evolution equations used to compute the eccentricity excitation and survival maps.","marker":"Liu & Lai 2018"},{"why":"Supplies the octupole-order and 1PN precession treatment of ZLK eccentricity excitation that sets the tidal-disruption boundary.","marker":"Liu et al. 2015"},{"why":"Defines the relativistic-precession parameter $\\varepsilon_{\\rm GR}$ used to organize the survival scans.","marker":"Huang, X. et al. 2026"}],"fun_headline_variants":["GSN 069's X-ray burst timing finds no binary black hole","No supermassive binary in GSN 069 from QPE timing","GSN 069's eruption clock excludes a black hole binary","QPE timing and stellar survival rule out binary in GSN 069","No sign of companion black hole in GSN 069's QPEs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The null result depends on the eruption cycle numbers used in the timing analysis being assigned correctly; a different assignment can spuriously create or erase the in-phase modulation that would signal a companion.","fun_headline_variants_meta":{"raw":{"variants":["GSN 069's X-ray burst timing finds no binary black hole","No supermassive binary in GSN 069 from QPE timing","GSN 069's eruption clock excludes a black hole binary","QPE timing and stellar survival rule out binary in GSN 069","No sign of companion black hole in GSN 069's QPEs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000787,"raw_usage":{"total_tokens":3525,"prompt_tokens":1050,"completion_tokens":2475,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":666,"completion_tokens_details":{"reasoning_tokens":2380}},"tokens_in":666,"tokens_out":2475,"duration_ms":15199,"temperature":1.0,"reasoning_tokens":2380,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T11:19:03.742264+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the GSN 069 eruption timings with an independently determined cycle-number assignment, or with a systematic shift of the even/odd labels; if a significant in-phase modulation appears (for example $\\log B>0$ favoring the binary model), the no-companion conclusion is an artifact of the numbering. An independent detection of a sub-parsec binary in GSN 069, for instance through radio imaging or gravitational-wave observations, would also falsify the exclusion.","supporting_citations":[],"review_version":1}