{"id":"2d6f48de-7d76-4a98-8538-bfd0dcc0895c","arxiv_id":"2504.16161","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Across three population synthesis codes, secondary stars ejected by supernova disruption leave at approximately their pre-supernova orbital velocity, nearly independent of natal kick strength; a correction to the Kiel & Hurley (2009) prescription removes an unphysical kick dependence.","lead":"This paper compares three computer codes that simulate binary star evolution and finds they agree that stars ejected when their companion explodes as a supernova travel at roughly their pre-explosion orbital speed, almost independent of the kick the exploding star receives. The result helps astronomers interpret Gaia measurements of runaway stars and fixes several software bugs that distorted earlier predictions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The abstract's unqualified claim ignores the SN blast-wave impulse on the companion, an effect the three codes do not model and whose importance the paper does not quantify for the disrupted binary population.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing concern: the codes neglect the SN blast-wave impulse, and the central claim depends on that neglect. I agree this is the most serious issue. The KH09 correction, while asserted without a full derivation, is secondary because the paper's main comparison uses Pfahl et al. (2002) and Tauris & Takens (1998), and the correction only affects the historical COSMIC prescription. The single-representative-binary issue is real but partially mitigated by the analytic explanation in Section 3.3, which generalizes the mechanism; the blast-wave neglect is a physical gap that no analytic argument in the paper addresses. The paper itself flags the effect in the introduction, so this is not an unfair demand. The proposed test would settle whether the effect matters: for the representative system it likely does not, but for closer disrupted binaries it could. Thus the reader's CONDITIONAL verdict remains appropriate; the paper should add a qualifier and ideally a quantitative estimate of the blast-wave impulse's contribution.","tokens_in":14325,"tokens_out":7799,"duration_ms":79428,"concrete_test":"Add the analytic ejecta–companion momentum-transfer prescription of Hirai et al. (2018) to one of the codes, or post-process the pre-SN states, for the representative binary and for a grid of post-mass-transfer separations a = 10, 50, 100, 300, 1000 R_sun with the same masses. If the induced companion velocity exceeds ~5 km/s for any separation at which disrupted binaries are expected, the abstract's unqualified claim must be qualified; if not, the concern is quantitatively refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the ejection velocities of secondary stars from disrupted binaries are narrowly distributed about their pre-Supernova orbital velocity. This result is obtained in all three codes, but only because each code treats the SN as instantaneous mass loss plus a natal kick to the compact object, with no momentum transfer from the SN ejecta to the companion. Section 1 explicitly lists the SN blast-wave impulse (Liu et al. 2015; Hirai et al. 2018; Ogata et al. 2021; Wong et al. 2024) and collisions as effects considered in other works, but none of the codes include them. The representative binary is very wide at explosion (P ~ 700 d, a ~ 5 AU), so ejecta–companion interaction is negligible for this system; the agreement in Figure 2 therefore validates the internal consistency of the kick algorithms, not the physical completeness of the claimed narrow distribution. The abstract and Section 5 present the narrow distribution as a property of the binary SN scenario without the qualifier 'in the absence of SN ejecta–companion interaction.' If blast-wave impulses add a velocity component comparable to or larger than the quoted ~5 km/s scatter for a non-negligible subset of disrupted binaries (e.g., systems with smaller post-mass-transfer separations), the central claim as stated is too strong. The paper does not quantify this, and the single representative initial condition cannot rule it out.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses three open-source binary population synthesis codes (COSMIC, COMPAS, binary_c) to evolve a representative massive binary (m1=20 Msun, m2=15 Msun, P=100 d, e=0, Z=0.02) 50,000 times with different natal kicks, applying a common kick magnitude distribution and a common remnant mass prescription. It finds that, across the three codes, the ejection velocities of companion stars from binaries disrupted by the first supernova are tightly clustered around the pre-SN orbital velocity of the companion, with a scatter of only a few km/s despite kicks spanning three orders of magnitude. The paper explains the residual dependence of the ejection velocity on kick direction and magnitude, reports and corrects a missing term in the Kiel & Hurley (2009) derivation, and documents bugs found in all three codes. It concludes that the code implementations are consistent and that companion ejection velocities are primarily set by the pre-SN orbital velocity.","tokens_in":14620,"tokens_out":6046,"duration_ms":57794,"significance":"The paper is a valuable cross-validation of three widely used population synthesis codes and has already led to concrete, publicly documented bug fixes in each of them. It is reproducible in spirit: the data and code are released on Zenodo and GitHub, the test uses 50,000 kick realizations per code, and the input physics is homogenized as far as is practical. The KH09 correction is practically important for interpreting older COSMIC results, and the geometric explanation in Section 3.3 is clear. If the central claim is accepted after appropriate qualification, it simplifies the interpretation of Gaia-era runaway-star velocities from the binary supernova scenario. The main caveat is that the claim, as written, is broader than the physics actually modeled by the three codes.","major_comments":[{"comment":"The central claim that ejection velocities of disrupted companions are narrowly distributed about their pre-SN orbital velocity is stated without the caveat that only instantaneous mass loss plus a natal kick to the compact object is modeled. Section 1 explicitly lists the SN blast-wave impulse on the secondary (Liu et al. 2015; Hirai et al. 2018; Ogata et al. 2021; Wong et al. 2024) and collisions as effects considered in other works, but none of the three codes include those effects. For the particular representative binary the post-mass-transfer separation is large (P ≈ 555–800 d), so ejecta–companion interaction is plausibly negligible for this system; however, the abstract and Section 5 present the narrow distribution as a property of the binary supernova scenario generally. The authors should either qualify the claim to “in the absence of SN ejecta–companion interactions” or quantitatively estimate the blast-wave contribution for the range of post-mass-transfer separations that produce disrupted binaries. As written, the unqualified claim may mislead the interpretation of runaway-star observations.","section":"Abstract; §1; §3.2; §5"},{"comment":"The quantitative support for the population-level claim rests on a single initial binary. The “within ~5 km/s” scatter is measured for one representative system, and the three codes evolve it to a fairly narrow range of pre-SN configurations (v2,preSN between about 9 and 15 km/s). The analytic discussion in §3.3 is general, but the quantitative statement in the abstract is not. To make the central claim robust, the authors should either show that the result holds across a grid of initial masses, periods, and metallicities, especially for systems with smaller post-mass-transfer separations where the compact object remains in the system longer and where blast-wave effects are stronger, or explicitly restrict the claim to the representative system.","section":"§2.2; §3.2; §5"}],"minor_comments":[{"comment":"The caption contains a duplicated word: “each main panel shows shows the difference” should read “each main panel shows the difference”.","section":"Fig. 2 caption"},{"comment":"The sentence “For a compact object that is ahead of the direction of the companion, its will pull the companion towards it” appears to be missing a noun; “its” should be “its gravity” or “it”.","section":"§3.3"},{"comment":"The phrase “These kicks can effect a greater change” should use “affect” rather than “effect”.","section":"§3.3"},{"comment":"There are two distinct 2002 papers by Pfahl et al. in the reference list, but the text cites “Pfahl et al. 2002” without distinguishing which one is intended; the authors should disambiguate these citations.","section":"§1; §2.1"},{"comment":"The variables ν and γ in the corrected KH09 equations are not re-defined in the text, so a reader must consult KH09’s Figure 2 to verify the sign convention; a brief sentence defining these angles would make the correction self-contained.","section":"§4, Eqs. (3) and (4)"},{"comment":"The exact versions of COSMIC and binary_c used for the simulations are not stated as explicitly as the COMPAS version (v3.01.10); giving the exact version numbers would strengthen reproducibility.","section":"§5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid code-comparison paper with concrete, useful bug fixes, and I see no ethical concerns. My main concern is scope: the title and abstract overclaim relative to the physics actually implemented in the three codes, because the SN blast-wave impulse on the companion is acknowledged but not modeled or quantified. This is fixable by qualification or by adding a quantitative estimate, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Tom,\n\nQuick read of Wagg et al. (arXiv:2504.16161). The paper is a careful comparison of secondary-star ejection velocities after supernova disruption across COSMIC, COMPAS, and binary_c. The genuinely new things are: (1) a correction to Kiel & Hurley (2009) that removes an unphysical kick-magnitude dependence in the ejection velocity, and (2) the identification of implementation bugs in all three codes—COSMIC's KH09 error, COMPAS's factor-of-2 overestimate in component velocities, and binary_c's correlated kick sampling. Those fixes are real and valuable, and the authors have made the corrected methods accessible and released data/code.\n\nThe core physical statement—that the companion's ejection velocity tracks its pre-SN orbital velocity rather than the kick magnitude—is not new (Tauris & Takens 1998; Pfahl et al. 2002; Renzo et al. 2019), but the paper confirms it cleanly with 50,000 kick realizations per code and consistent kick/remnant prescriptions. The agreement across three independent implementations is convincing evidence that the two kick algorithms are equivalent and correctly implemented.\n\nThe soft spot is that the abstract and Section 5 state the narrow distribution as a property of the binary supernova scenario without the qualifier 'in the absence of SN ejecta–companion interaction.' The codes treat the SN as instantaneous mass loss plus a natal kick; they do not model the blast-wave impulse on the companion, and the representative binary is wide enough (P ~ 700 d at explosion) that the impulse is negligible for this system. But for closer disrupted binaries—especially those with smaller post-mass-transfer separations—blast-wave effects can add a non-negligible velocity component (Liu et al. 2015; Hirai et al. 2018; Ogata et al. 2021; Wong et al. 2024). The paper cites these works but does not quantify the regime where the conclusion holds. That is a genuine gap in the abstract's claim, though it does not undermine the internal consistency of the code comparison or the bug fixes.\n\nMinor points: the KH09 correction is stated with the missing terms but no full derivation; readers will have to go to KH09 and the code to verify. And the single representative binary limits the generality of the result, though the physical explanation in Section 3.3 is general.\n\nWho's this for? Population synthesis developers and anyone using COSMIC/COMPAS/binary_c results for runaway-star interpretation. If you use those codes, the bug fixes alone justify reading it. The blast-wave caveat matters for applying the result to Gaia observations of runaways.\n\nMy recommendation: send it out for review. The code fixes and cross-code comparison are worth refereeing, and the authors should be pushed on the unqualified claim in the abstract.\n\nBest,\n[You]","headline":"Solid code-comparison paper: the KH09 correction and the three identified bugs are real, the cross-code agreement is credible, but the abstract overstates the result by omitting the SN blast-wave impulse.","tokens_in":15102,"tokens_out":5404,"would_cite":true,"duration_ms":43101,"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":"Ejected companions from supernova-disrupted binaries leave at close to their pre-supernova orbital speed, not at the speed of the natal kick, and this holds across three population-synthesis codes.","keywords":["binary supernova scenario","runaway stars","stellar ejection velocities","natal kicks","population synthesis","post-supernova orbital dynamics","Kiel & Hurley correction","Gaia astrometry"],"falsifier":"Population-level test: measure space motions of a clean sample of Gaia runaway B/O stars whose pre-supernova orbits can be reconstructed from their evolutionary state, such as companions of young neutron stars with known radial-velocity history. If the ejection velocities show a tail toward hundreds of km/s, or scale with the natal kick magnitude inferred from the pulsar, the claim fails. A single well-constrained disrupted binary with measured companion speed more than 5 km/s away from its reconstructed pre-SN orbital velocity would also contradict the prediction.","tokens_in":14181,"feed_emoji":"💫","tokens_out":10398,"duration_ms":88718,"temperature":0.7,"pith_summary":"This paper establishes that when the first supernova disrupts a massive binary, the surviving companion star is flung away at a speed set almost entirely by its pre-supernova orbital velocity, with the natal kick of the newly formed compact object playing only a minor role. The authors test this across three independent population-synthesis codes that use two different kick algorithms, evolving the same representative binary 50,000 times in each code, and find the ejection velocities are centred on the pre-supernova orbital velocity to within about 5 km/s even though kicks span three orders of magnitude. The result matters because runaway and walkaway stars, now catalogued in large numbers by Gaia, are often interpreted as products of this binary-supernova channel; if this paper is right, their measured speeds report the binary's orbital separation at the moment of explosion, not the kick physics. The paper also corrects a missing term in the Kiel & Hurley (2009) kick prescription that had created an unphysical dependence of ejection velocity on kick magnitude, and documents bugs found in all three codes.","feed_headline":"Supernova kicks barely change the companion's ejection speed","feed_subtitle":"Three independent population-synthesis codes agree that runaway stars leave at near their pre-supernova orbital velocity.","key_machinery":"The load-bearing mechanism is the near-instantaneous removal of the primary at core collapse. Because typical natal kicks (hundreds of km/s) are much larger than the binary's orbital speed (roughly 10 km/s here), the compact object leaves the companion's vicinity almost immediately; the companion, no longer held by the centripetal force, coasts away at close to its pre-SN orbital velocity. The two kick algorithms compared—Pfahl et al. (2002), a vector derivation built on the Laplace–Runge–Lenz vector, and Tauris & Takens (1998), a coordinate system aligned with the companion's direction—are shown to be equivalent. The paper's correction to Kiel & Hurley (2009) supplies the missing $V_\\infty \\sin\\nu \\sin\\gamma$ term in the $z$-components of both final velocities, removing an artificial kick-magnitude dependence.","core_discovery":"On the paper's own terms: for a mass-transferring binary that is disrupted at core collapse, the ejection velocity $v_{2,\\mathrm{postSN}}$ of the secondary star is narrowly distributed about its pre-SN orbital velocity $v_{2,\\mathrm{preSN}} = (m_1/(m_1+m_2))v_{\\mathrm{orb}}$, with deviations from the mean within roughly $5\\,\\mathrm{km\\,s^{-1}}$ despite natal kick magnitudes up to $\\sim 1000\\,\\mathrm{km\\,s^{-1}}$. The same holds in COSMIC, COMPAS, and binary_c even though the codes disagree on the pre-SN state: the companion masses are roughly $18\\,M_\\odot$, $28\\,M_\\odot$, and $29\\,M_\\odot$ for COMPAS, COSMIC, and binary_c, with pre-SN orbital velocities $14.6\\,\\mathrm{km\\,s^{-1}}$, $9.2\\,\\mathrm{km\\,s^{-1}}$, and $12.9\\,\\mathrm{km\\,s^{-1}}$ respectively. The authors explain the small residual scatter through the geometry of the kick: strong kicks remove the compact object almost instantaneously, while weaker in-plane kicks can slightly accelerate or decelerate the companion before the compact object leaves, and perpendicular kicks alter the speed least. The paper further claims that its correction to the Kiel & Hurley (2009) equations—adding the missing $V_\\infty$ term to the $z$-components of the final velocities—brings that prescription into agreement with Pfahl et al. (2002) and Tauris & Takens (1998), and that the two algorithms yield identical ejection velocities for identical pre-SN input.","pith_inferences":["If the claim generalises beyond this one representative binary, runaway-star observations become easier to invert: the kick velocity drops out of the ejection-speed relation, leaving the pre-SN orbital separation and mass ratio as the main unknowns.","The same near-instantaneous-removal logic should break down for very tight binaries whose orbital speed is a larger fraction of the kick speed; widening the grid of initial conditions could reveal where the constant-velocity approximation fails.","A natural extension, not pursued here, is to run the corrected kick routines over full galactic populations and compare the predicted runaway velocity distribution directly with Gaia tangential velocities as a function of spectral type.","Because blast-wave impulse and star-compact-object collisions are excluded from all three codes, the narrow velocity prediction is only established for the wide, post-mass-transfer binaries of this study; hydrodynamical tests would show whether close binaries behave differently."],"forward_implications":["Gaia-measured runaway and walkaway speeds become a direct probe of the pre-explosion orbital period, because the ejected companion's velocity is approximately the orbital velocity at disruption.","The natal kick distribution shapes runaway populations mainly by deciding which systems are disrupted and which stay bound, not by setting the speed of the ejected companion.","Results from rapid population synthesis runs that used the original Kiel & Hurley prescription should be revisited; the corrected equations remove a spurious correlation between ejection velocity and kick magnitude.","Because the two independent kick algorithms agree for identical pre-SN inputs, remaining differences between codes in ejection velocities can be attributed to pre-SN evolution rather than to kick implementation.","The documented bugs in COSMIC, COMPAS, and binary_c mean published ejection velocities from the affected versions may need to be re-derived."],"supporting_citations":[{"why":"Establishes the binary supernova scenario: instantaneous mass loss changes the post-SN orbit and can eject the companion.","marker":"Blaauw 1961"},{"why":"Supplies the kick prescription used by COSMIC and COMPAS, with a vector derivation of post-SN orbital parameters.","marker":"Pfahl et al. 2002"},{"why":"Supplies the kick prescription used by binary_c; the paper shows both this and Pfahl et al. give identical ejection velocities.","marker":"Tauris & Takens 1998"},{"why":"The prescription whose missing V_infinity z-term is corrected, removing an unphysical kick-magnitude dependence.","marker":"Kiel & Hurley 2009"},{"why":"Establishes that disruption depends on the ratio of kick speed to pre-SN orbital speed, used to explain the disruption fractions.","marker":"Kalogera 1996"},{"why":"Provides context for walkaway velocities and symmetric mass-loss disruption rates that this result extends.","marker":"Renzo et al. 2019"},{"why":"Supplies the Maxwellian kick dispersion (sigma = 265 km/s) used to draw natal kick magnitudes in all three codes.","marker":"Hobbs et al. 2005"}],"fun_headline_variants":["Runaway stars barely feel supernova kicks","Binary ejection speeds match pre-supernova orbits","Three codes agree: ejection speed ≈ orbital speed","Supernova kicks don't boost companion ejection speed","Companion ejection velocity pinned by pre-SN orbit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The assumption that carries the result is that the supernova blast wave does not directly push the companion, so the companion's final velocity comes only from the sudden loss of the primary's mass and the brief gravitational pull of the fleeing compact object; the paper notes that blast-wave impulse and star-ejecta collisions are effects considered elsewhere but not included in the three codes.","fun_headline_variants_meta":{"raw":{"variants":["Runaway stars barely feel supernova kicks","Binary ejection speeds match pre-supernova orbits","Three codes agree: ejection speed ≈ orbital speed","Supernova kicks don't boost companion ejection speed","Companion ejection velocity pinned by pre-SN orbit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000146,"raw_usage":{"total_tokens":1262,"prompt_tokens":1108,"completion_tokens":154,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":724,"completion_tokens_details":{"reasoning_tokens":82}},"tokens_in":724,"tokens_out":154,"duration_ms":2322,"temperature":1.0,"reasoning_tokens":82,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:10:22.587572+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Population-level test: measure space motions of a clean sample of Gaia runaway B/O stars whose pre-supernova orbits can be reconstructed from their evolutionary state, such as companions of young neutron stars with known radial-velocity history. If the ejection velocities show a tail toward hundreds of km/s, or scale with the natal kick magnitude inferred from the pulsar, the claim fails. A single well-constrained disrupted binary with measured companion speed more than 5 km/s away from its reconstructed pre-SN orbital velocity would also contradict the prediction.","supporting_citations":[],"review_version":1}