{"id":"78052ded-efe4-4b33-ab83-0ab1e01353d3","arxiv_id":"2507.03644","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Across thirteen black hole X-ray binaries, type-A QPOs show negative lags and often precede radio flares, type-B QPOs accompany the flares, and inferred soft-intermediate-state jet velocities are moderately relativistic.","lead":"The paper compares X-ray timing and spectral measurements with radio flare timings across thirteen black hole X-ray binaries and reports systematic patterns linking oscillation type, source inclination, and time lag sign. It uses those patterns to argue that the hot corona reshapes before jet ejection and estimates jet speeds near 0.3 to 0.8 times the speed of light.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (10) places eta_acc in the numerator of the accretion-rate formula, but L_x = eta_acc * Mdot * c^2 puts it in the denominator; this biases every beta in Table 3 and may push corrected velocities above c.","rationale":"The qualitative QPO-lag synthesis is plausible and supported by a broad multi-source dataset; the type-A negative-lag precursor pattern and the type-B flare association are valuable working hypotheses, and the paper deserves credit for compiling extensive archival data. However, the abstract's quantitative claim (jets in SIMS have velocities >= 0.3-0.8c) rests on the jet-power equality in Sec. 5.2. The reader flagged the epsilon = Delta_F_nth identification as the weakest assumption; that is legitimate and the authors acknowledge it in Sec. 6.3. My stress-test identifies a more basic problem in Eq. (10), where eta_acc appears in the numerator. The correct expression, obtained from L_x = eta_acc * Mdot * c^2, has eta_acc in the denominator. Because Mdot_Edd as defined already contains a 10% efficiency factor, the published formula underestimates Mdot by roughly 1/eta_acc^2, about 10x for eta_acc = 0.1. Since beta scales as (epsilon*Mdot)^(7/9), the resulting velocities would change by a factor of several and often exceed c, indicating an inconsistency between the radio power estimate, the accretion-rate normalization, and the assumed epsilon values. This is a concrete, checkable algebraic issue rather than a modeling preference, so it should be settled before the velocity numbers are used. The empirical lag patterns and coronal-geometry interpretation are less affected and can stand as a conditional interpretation of the data.","tokens_in":43493,"tokens_out":14659,"duration_ms":150395,"concrete_test":"Independently re-derive Eq. (10) from L_x = 4*pi*D^2*F_x and L_x = eta_acc * Mdot * c^2, obtaining Mdot = (4*pi/c^2)*F_x*D^2/eta_acc. Then recompute Table 3 beta values with this corrected Mdot (recalibrating epsilon for the proper-motion sources if needed). If corrected beta exceeds 1 for several sources, or the 0.3-0.8c range is not recovered, the paper's quantitative jet-velocity claim is not supported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The quantitative jet-velocity result depends on Eq. (10), which reads Mdot = 8.73e-17 * eta_acc * (F_x D^2 / c^2) * (M_BH/M_sun)^-1 * Mdot_Edd, with Mdot_Edd = 1.47e18 (M_BH/M_sun) g/s. After the mass term cancels, the formula gives Mdot proportional to eta_acc * F_x D^2 / c^2. But the stated relation L_x = 4*pi*D^2*F_x = eta_acc * Mdot * c^2 requires Mdot = (4*pi/c^2) * F_x D^2 / eta_acc. Thus Eq. (10) has eta_acc in the wrong place; for eta_acc ~ 0.1 it underestimates Mdot by roughly a factor of 10. Since beta is obtained by equating L_jet proportional to beta^(9/7) with (1/2)*eta_jet*epsilon*Mdot*c^2, beta is proportional to (epsilon*Mdot)^(7/9); a 10x correction in Mdot raises beta by ~6x, pushing many entries in Table 3 above unity. The authors' caveat in Sec. 6.3 that epsilon = Delta_F_nth 'may be overly simplistic' is a legitimate concern, but it is secondary: even granting that identification, the accretion-rate normalization in Eq. (10) must be corrected before any beta can be trusted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a spectro-temporal study of 16 outbursts from 13 black hole X-ray binaries, combining RXTE, HXMT, and AstroSat X-ray data with radio observations from the literature. The authors classify type-A, type-B, type-C, and type-C* QPOs, measure their time lags as a function of source inclination, and associate type-A QPOs with pre-radio-flare epochs and type-B QPOs with the flaring epochs. They argue that the corona evolves from a radially extended to a vertically elongated structure across the type-C-to-type-B transition. In the second half of the paper, they estimate black hole spins by continuum fitting and derive jet velocities by equating a minimum-energy radio-blob power with an assumed fraction of the accretion power, obtaining velocities in the range 0.3-0.8 c during the soft intermediate state.","tokens_in":43822,"tokens_out":6204,"duration_ms":73799,"significance":"If the empirical QPO-lag-inclination pattern holds, the paper would provide a useful phenomenological map connecting QPO subtypes, coronal geometry, and jet ejection, and the assembled lag measurements across many sources would be a valuable reference. The paper is also honest in stating some limitations, particularly in Section 6.3. However, the quantitative jet-velocity claim is not robust: the accretion-rate formula in Equation (10) has a self-contradictory dependence on radiative efficiency, and the epsilon values in Table 3 are either calibrated to the very velocities being predicted or identified with flux changes in a way the authors themselves call 'overly simplistic.' Since the headline '0.3-0.8 c' result depends on this chain, the central quantitative claim is currently not established.","major_comments":[{"comment":"Equation (10) places eta_acc in the numerator of the accretion-rate formula, but the stated relation L_x = eta_acc * Mdot * c^2 requires Mdot to be proportional to F_x D^2 / (eta_acc * c^2). With eta_acc ~ 0.1, the formula underestimates Mdot by roughly an order of magnitude (in addition to a missing 4 pi factor in the flux-to-luminosity conversion). Because Equation (9) gives beta proportional to (epsilon * Mdot)^(7/9), this error changes every beta in Table 3 by a factor of several and can push entries such as the 0.96-0.98 value for XTE J1752-223 above unity. The authors must correct Equation (10) and recompute the table before the velocity claims can be assessed.","section":"Section 5.2, Eq. (10)"},{"comment":"The epsilon values used to predict beta are not determined independently. For H1743-322 (2003 and 2009), XTE J1550-564, MAXI J1535-571, and XTE J1752-223 (F4), epsilon is chosen so that Equation (9) reproduces previously measured proper-motion velocities; the same epsilon convention is then applied to all other sources. Since beta is a monotone function of epsilon in Equation (9), agreement with the known velocities is partially by construction. The additional identification of epsilon with Delta F_nth is acknowledged in Section 6.3 to be 'overly simplistic,' but that caveat applies directly to the entries in Table 3, so the stated velocity range of 0.3-0.8 c is not independently established.","section":"Section 5.2 and Table 3"},{"comment":"Equating the minimum-energy synchrotron power of a radio blob, W_min/t, with the kinetic jet power (1/2) eta_jet epsilon Mdot c^2 is a strong model assumption. W_min is a minimum total energy content of synchrotron-emitting plasma, not necessarily the jet kinetic power, and the adopted values eta_jet = 0.1 and k = 2 or 3 are fixed without a sensitivity analysis. Because beta depends on the 7/9 power of the resulting power ratio, plausible variations in these parameters could shift the headline velocities by tens of percent, and the paper should quantify this uncertainty before presenting 0.3-0.8 c as a precise result.","section":"Section 5.2, Eqs. (6)-(9)"},{"comment":"The central empirical claim that type-A QPOs show negative lags independent of inclination and that type-B QPOs coincide with radio flares rests on manual QPO classification and on a small number of sources, with four high-inclination systems (XTE J1550-564, Swift J1727.8-1613, H1743-322 2003, and GRO J1655-40) explicitly listed as exceptions to the type-C lag-inclination trend. The paper should provide a per-source table of lag signs with uncertainties and a systematic treatment of the exceptions rather than setting them aside, especially because the claimed 'regardless of inclination' property of type-A QPOs is based on only 26 detections across the sample.","section":"Section 4.3 and Figs. 6-7"}],"minor_comments":[{"comment":"The time-lag formula should read delta_t(j) = arg[C(j)] / (2 pi nu_j); as written, delta_t(j) = C(j) / (2 pi nu_j) is a complex quantity, not a real time delay.","section":"Section 3.1, Eq. (5)"},{"comment":"The abstract lists the GX 339-4 outbursts as 2002, 2006, and 2010, while Section 2.1 and Table 1 refer to 2002, 2007, and 2010; the year 2006 appears to be a typo and should be corrected consistently.","section":"Abstract and Section 2.1"},{"comment":"The 4U 1543-47 row lists the outburst year as 2004, but the text and Section 4.2 consistently describe the 2002 outburst; this should be made consistent.","section":"Table 1"},{"comment":"The caption refers to 'Swift J1727.8-0127,' whereas all other parts of the paper use Swift J1727.8-1613; the figure caption should be corrected.","section":"Figure 10 caption"},{"comment":"The high-inclination correlation coefficients are quoted after excluding GRO J1655-40, but the figure does not state why this source is excluded; a brief justification should be given in the text or caption.","section":"Figure 9"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read for you. The QPO-lag part of this paper is worth a careful look; the jet velocity part is not. The genuinely new thing is the systematic multi-source statement that type-A QPOs show negative lags independent of inclination and appear before radio flares in several systems, with type-B QPOs clustering at the flares. That pattern is not, to my knowledge, in the prior literature in this compiled form, and the corona geometry story—radial extended to vertically elongated—is a plausible qualitative interpretation. The authors are also honest that identifying ε with the change in normalized Comptonized flux is crude (Sec. 6.3).\n\nThe load-bearing problem is Eq. (10). They define L_x = 4πD^2F_x = η_acc Mdot c^2, which requires Mdot ∝ 1/η_acc. The printed formula has η_acc in the numerator. The stress-test note caught the right error but understated it: for η_acc ≈ 0.1 the accretion rate is off by roughly η_acc^{-2} ≈ 10^2, not 10. Since β ∝ (ε Mdot)^{7/9}, the corrected velocities scale as η_acc^{-14/9} ≈ 36, so most of Table 3 goes superluminal. This is not a calibration nuance; the velocities as printed are wrong.\n\nSmaller soft spots: QPO typing is manual, with four sources set aside as exceptions rather than modeled; some correlations lack error bars; and the abstract says GX 339-4's 2006 outburst while the text and Table use 2007. These are fixable.\n\nWho is this for? People working on QPO phenomenology and the jet line. The lag–flare sequence deserves serious referee time, but only with the jet velocity section either corrected or demoted to an order-of-magnitude sketch. I would not cite the 0.3–0.8c numbers until the algebra is fixed. Recommendation: send to peer review with a referee who checks the equations, and require a corrected Eq. (10) and re-derived Table 3.","headline":"The QPO lag–flare census is worth a careful look, but the jet velocities are vitiated by an algebra error in Eq. (10).","tokens_in":44373,"tokens_out":6846,"would_cite":false,"duration_ms":69003,"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":"Surveying 16 outbursts of 13 black hole X-ray binaries, the paper reports that type-A QPOs consistently show negative time lags regardless of source inclination, precede radio flares, and signal jet ejection, while type-B QPOs coincide…","keywords":["black hole X-ray binaries","quasi-periodic oscillations","time lags","corona geometry","jet ejection","relativistic jets","radio-X-ray correlation","accretion states"],"falsifier":"Measure the proper motion of the radio ejecta during the soft intermediate state in a source whose velocity is predicted here (for example 4U 1543-47 or XTE J1752-223) and compare it with the paper's Table 3; a resolved speed clearly outside the predicted range would rule out the flux-change identification, while a single outburst with type-A QPOs appearing only after the radio flare peak would undercut the precursor claim.","tokens_in":43271,"feed_emoji":"🕳️","tokens_out":8231,"duration_ms":83942,"temperature":0.7,"pith_summary":"This paper surveys sixteen outbursts of thirteen black hole X-ray binaries, combining X-ray timing and spectra with quasi-simultaneous radio data, to test how the corona changes as jets are launched. It reports a systematic pattern: type-C QPOs have inclination-dependent time lags, type-A QPOs always show negative lags and tend to precede radio flares, and type-B QPOs coincide with the flares. The authors interpret this as the corona shrinking radially and then elongating vertically during the transition from the hard intermediate state to the soft intermediate state, with type-B QPOs tracing a compact or vertically extended corona that resembles a jet base. They further estimate that jets in the soft intermediate state are moderately relativistic, at velocities >0.3–0.8c.","feed_headline":"Type-A QPOs precede radio flares; jets reach 0.8c","feed_subtitle":"A 13-source survey ties QPO lag signs to corona shape and shows intermediate-state jets are moderately relativistic.","key_machinery":"The argument is carried by the time lag of each QPO type—the Fourier phase difference between the 2–6 keV and 6–15 keV lightcurves at the QPO frequency—used as a geometric probe of the corona. A positive lag is read as soft photons being Compton up-scattered in an extended corona before reaching the observer, while a negative lag is read as hard photons being reprocessed back into the disc, an effect that grows when the corona is small or the viewing angle is high. On the jet side, the machinery is the minimum-energy synchrotron formula for a radio-emitting blob, Doppler-corrected and equated through Equation (9) to a fraction of the accretion power, with that fraction set by the change in normalized Comptonized flux between successive X-ray observations around a flare.","core_discovery":"The central claim is that the sequence of QPO types encodes the geometry of the corona leading up to jet ejection. For thirteen sources the authors find that type-C QPOs, seen in harder states, show positive lags for low-inclination systems and negative lags for high-inclination systems, consistent with a large radially extended corona. Type-A QPOs appear near the state transition with negative lags of about 1–10 ms in every source regardless of inclination, and in several outbursts they show up before the radio flare, identifying them as precursors of jet ejection. Type-B QPOs, observed in the soft intermediate state with lower Comptonized flux, coincide with the radio flares and show positive lags in low-inclination sources and mixed lags in high-inclination sources, which the authors read as evidence for a radially compact or vertically elongated corona. Finally, using a minimum-energy jet model normalized to measured proper-motion speeds, the paper estimates jet velocities above 0.3–0.8c during the soft intermediate state, and concludes that the strong radio–X-ray correlation indicates accretion-powered jets.","pith_inferences":["The identification of the Comptonized flux change with the fraction of accretion power carried into the jet could be tested directly in a source with simultaneous X-ray and radio monitoring by separating thermal and Comptonized spectral components at high cadence; if the Comptonized flux drop tracks the radio lightcurve on timescales shorter than the accretion timescale, the outflow interpretation","If the precursor role of type-A QPOs holds up in future outbursts, it suggests that the jet launch itself modifies the corona, and that the delay between type-A and type-B QPOs may measure the vertical growth time of the jet base.","The log–log slope of about 0.33 between jet velocity and X-ray luminosity found here is a candidate scaling for jet-launching models and could be tested against theoretical predictions for radiatively inefficient accretion flows.","An analogous lag-sign ordering might be searched for in X-ray data of active galactic nuclei, where type-C-like low-frequency QPOs are observed, to see whether the geometry sequence is universal."],"forward_implications":["Type-A QPOs can serve as a practical early warning that a transient jet is about to be launched, enabling coordinated multi-wavelength follow-up.","The sign of the type-C QPO lag can be used to infer the inclination of a binary, or the radial extent of its corona, in sources without dynamical inclination measurements.","Type-B QPOs with positive lags identify a vertically elongated corona at the jet base, linking a timing signature directly to jet geometry.","Moderate SIMS jet velocities imply only modest Doppler boosting, changing how intrinsic radio luminosities should be estimated for intermediate-state jets.","The correlation between jet velocity and X-ray luminosity found in this sample suggests that accretion rate, rather than spin, is the main driver of jet speed."],"supporting_citations":[{"why":"Defines the jet line and the HIMS-to-SIMS transition where transient jets are launched; supplies the k=2 vs k=3 ejecta assumptions for the Doppler correction.","marker":"Fender et al. 2004"},{"why":"Previous finding that type-C and type-A QPO lags depend on inclination, which this paper extends to a broader sample and to type-B QPOs.","marker":"van den Eijnden et al. 2017"},{"why":"Gives the minimum-energy synchrotron power formula used to convert radio flux into jet power.","marker":"Longair 2011"},{"why":"Provides the jet radiative efficiency of about 0.1 assumed in the jet-power equation.","marker":"Fender 2001a"},{"why":"Suggested a jet-like corona for type-B QPOs, supporting the vertical-corona interpretation.","marker":"Belloni et al. 2020"},{"why":"Established the radio–X-ray luminosity correlation used here as evidence for accretion-powered jets.","marker":"Corbel et al. 2003"}],"fun_headline_variants":["Type-A QPOs herald jets; type-B track radio flares","Corona evolution traced by QPO lags, tying jets to 0.8c speeds","QPOs decode corona geometry and predict jet velocities up to 0.8c","13 black hole binaries: QPO signatures precede jet ejection and flares","Jet speed up to 0.8c linked to QPO type and corona shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The jet velocities rest on the assumption that the decrease in normalized Comptonized flux between two X-ray observations around a radio flare measures the fraction of accretion power carried into the jet; if that flux change instead reflects a change in the accretion rate, the reported speeds do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Type-A QPOs herald jets; type-B track radio flares","Corona evolution traced by QPO lags, tying jets to 0.8c speeds","QPOs decode corona geometry and predict jet velocities up to 0.8c","13 black hole binaries: QPO signatures precede jet ejection and flares","Jet speed up to 0.8c linked to QPO type and corona shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000479,"raw_usage":{"total_tokens":2550,"prompt_tokens":1301,"completion_tokens":1249,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":917,"completion_tokens_details":{"reasoning_tokens":1143}},"tokens_in":917,"tokens_out":1249,"duration_ms":11635,"temperature":1.0,"reasoning_tokens":1143,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:05:52.160920+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the proper motion of the radio ejecta during the soft intermediate state in a source whose velocity is predicted here (for example 4U 1543-47 or XTE J1752-223) and compare it with the paper's Table 3; a resolved speed clearly outside the predicted range would rule out the flux-change identification, while a single outburst with type-A QPOs appearing only after the radio flare peak would undercut the precursor claim.","supporting_citations":[],"review_version":1}