{"id":"072ecde7-a817-481e-ae8e-fc394b885e9d","arxiv_id":"2411.16847","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A Gaia-based kinematic comparison shows black holes have mixed natal kicks: several are much faster than 90% of their local stars, while others match the local population, supporting both strong-kick and weak-kick formation channels.","lead":"This paper compares the motions of 12 black holes in binary systems with the motions of nearby stars, using Gaia data, and finds that some black holes were kicked hard at birth while others were not. The result matters because natal kicks determine which black hole binaries survive, merge, and produce gravitational waves.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The percentile classification rests on an age-unmatched comparison sample: local red giants are assumed to represent each BH's birth population, so the 6 'no-kick' systems could hide real kicks; the paper itself flags this for Cyg X-1 but does not test it.","rationale":"The reader's CONDITIONAL verdict is appropriate. The most load-bearing assumption is the one the paper itself labels implicit in Section 3.1: local reference samples dominated by old red giants represent the birth population of each BH binary. This assumption directly controls the velocity percentile and the kick/no-kick split in Table 2 and underlies the abstract's '6 BHs have kinematics typical of their local populations' claim. The paper's caveats admit the age mismatch for Cyg X-1 but do not quantify or correct it. I also note that the printed Equation (4) is inconsistent with Table 2: the table values equal Vpec,local - V68%, not V_total - V68%, where V_total includes the full azimuthal velocity. This looks like a typo-level error rather than the central flaw, but it should be corrected. The absence of propagated uncertainties on velocity percentiles and minimum kicks is a related weakness: MAXI J1820+070 flips from no-kick in Table 2 to a small-kick classification in Table 3 between two defensible metrics, so the 6/6 split is fragile. None of this overturns the qualitative mixed-origin conclusion, because V404 Cyg and VFTS 243 provide independent weak-kick evidence, and the strong-kick systems exceed the 90th percentile of even an old-giant baseline. However, the quantitative half/no-kick count should be treated as conditional on age matching. The verdict should remain CONDITIONAL, with the authors asked to add an age-matched comparison test or soften the '6 typical' wording.","tokens_in":27518,"tokens_out":6555,"duration_ms":69937,"concrete_test":"Recompute Table 2 with age- and chemistry-matched reference samples. For each of the 12 targets, select comparison stars from Gaia DR3 that plausibly share the binary's birth population: for young HMXBs such as Cyg X-1, use young disk stars identified from isochrones, [alpha/Fe], or OB association membership; for old LMXBs, use red giants selected to match the system's estimated age. Recompute V68, V90, V95 and velocity percentiles exactly as in Section 3.2, then count how many systems cross the 68% and 90% thresholds relative to Table 2. If the no-kick count drops from 6 to 4 or fewer, or if any of the six 'typical' systems, especially Cyg X-1 or MAXI J1820+070, moves above the 68th percentile, the headline 'half no kick' is not robust to the age-mixing assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the construction of the no-kick baseline. Section 3.1 builds each reference sample from Gaia DR3 stars within 500 pc and D +/- Delta D, requiring only RV error < 5 km/s and tangential-velocity error < 30 km/s, and the paper admits these are mostly red giants. The percentile radii V68/V90/V95 and the resulting kick/no-kick split in Table 2 are then interpreted as natal-kick evidence. But a BH binary and a random local red giant do not necessarily share a birth population: red giants have been dynamically heated over Gyrs, and their present-day velocity dispersion is larger than that of a younger or thin-disk birth population. For a young system such as Cyg X-1, the paper itself notes in Section 4.2.10 that most comparison stars are older and that an ideal comparison sample would be age-matched; no age or chemistry selection is applied. The bias is not conservative: comparing a kicked BH to an overly hot, old population shifts it toward the 'no-kick' side of the 68% threshold. Thus the claim that 6 BHs have kinematics typical of their local populations, disfavoring kicks of about 50 km/s or more, depends directly on an untested assumption. The Section 4.4 disk-crossing test changes Vpec by about 10 km/s but does not address age matching; it reuses the same comparison-sample construction. Because the mixed-origin conclusion also has independent support from V404 Cyg and VFTS 243, this is a robustness concern rather than a refutation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper estimates natal kicks for 12 Galactic black-hole binaries by comparing their Gaia-based space velocities and Galactic orbits with local comparison samples drawn from Gaia DR3 (500 pc projected radius, distance-matched, with quality cuts on parallax and velocity errors; the samples are mostly red giants). Using Toomre diagrams, the authors define 68/90/95% velocity thresholds for each local population and classify each BH as showing or not showing evidence of a kick, reporting 'minimum kicks' as the excess over the 68% threshold. They find that 6 of 12 systems lie beyond the 68th percentile, including 4 beyond the 90th (Swift J1727.8-162, GRO J1655-40, GRS 1124-684, Gaia BH1), while 6 lie within it (MAXI J1820+070, MAXI J1305-704, V404 Cyg, 3A 0620-003, Cyg X-1, Gaia BH2). Combining these classifications with independent dynamical constraints on V404 Cyg (kick ≲ 5 km/s from triple stability) and VFTS 243 (≲ 10 km/s from orbital circularity), they argue that the data are consistent with a mixed origin: some BHs form with strong kicks and some with negligible kicks, while cautioning that binary selection biases favor weak kicks and that a bimodal distribution is not directly proven. An appendix compares the kinematics of BH and NS binaries.","tokens_in":98,"tokens_out":20500,"duration_ms":389773,"significance":"If the central claims hold, the paper makes a useful methodological contribution: it replaces the cold thin-disk birth assumption of earlier kick studies (e.g., Atri et al. 2019; Zhao et al. 2023) with empirically measured local velocity dispersions, yielding more conservative minimum kicks and a population-level dichotomy rather than a unimodal distribution. Strengths worth emphasizing: the no-kick baseline is measured from Gaia DR3 percentiles rather than fitted; the paper is transparent about assumptions (age-match caveat in §3.1 and §4.2.10, arbitrariness of the percentile choice in §3.2, birth-location test in §4.4); it quantifies the binary selection bias toward weak kicks (§4.4); and it checks each system against dedicated literature studies (§4.2). The strong-kick conclusion is robust to the comparison-sample construction (Swift J1727.8-162 and GRO J1655-40 exceed the 95th percentile), and the existence of at least some very-weak-kick BHs rests on independent dynamical evidence (V404 Cyg, VFTS 243), so the mixed-origin conclusion survives even if individual kinematic classifications shift.","major_comments":[{"comment":"The reported minimum kicks cannot be reproduced from the written equations. Equation (4) defines V as the full Galactocentric 3D speed sqrt(V_R^2 + V_z^2 + V_phi^2), which for a thin-disk system is dominated by the azimuthal component V_phi ~ 230 km/s, so max(0, V - V68%) would be roughly 100-200 km/s for every system and could never be zero; yet six rows of Table 2 report V_kick,min = 0. Every row instead satisfies V_kick,min = V_pec,local - V68% within rounding (e.g., GRO J1655-40: 138 - 76 = 62; GRS 1124-684: 115 - 67 = 48; XTE J1118+480: 123 - 87 = 36; Gaia BH1: 79.4 - 50.0 = 29; and the zero rows correspond to V_pec,local < V68%). The intended estimator is evidently the target's excess beyond the 68% Toomre radius in the (V_phi, V_perp) plane, which is well defined, but it is not what Eqs. (4)-(5) state. Please correct the equations or the table, and align the numbers quoted in Sections 4.2.4, 4.2.6, 4.2.7, and 4.2.11 with the corrected definition.","section":"§3.2, Eqs. (4)–(5), Table 2"},{"comment":"The no-kick baseline rests on an untested age-matching assumption. The comparison samples are dominated by red giants (old, dynamically heated), and the paper states only that it assumes the BH binaries and reference samples probe similar stellar populations; for Cyg X-1 it explicitly acknowledges that most comparison stars are older than the system, yet no age or chemistry selection, and no test of the resulting bias, is provided. Because old giants have a larger V68% than a young thin-disk birth population, the comparison is not conservative for young systems: Cyg X-1's peculiar velocity of 24.2 km/s (18th percentile) could be consistent with a sizable kick if judged against an age-matched sample, and the §4.2.10 statement that the present-day peculiar velocity is a 'reasonable estimate' of the kick is in tension with the concern just raised. The §4.4 disk-crossing test addresses birth location, not birth population. Please quantify the age-mixing bias (e.g., by stratifying the comparison sample by color or isochrone age, or by constructing a young-star comparison sample for Cyg X-1) and revise the claim that six BHs have kinematics that disfavor kicks of ≳50 km/s accordingly.","section":"§3.1, §4.2.10, §4.4"},{"comment":"The kick/no-kick classification is quoted without propagated uncertainties. The 10,000-sample Monte Carlo described in §3.1 is used only to assign error bars to V_pec,local and V_pec,circ; the velocity percentiles, the V68%/V90%/V95% thresholds, the V_kick,min values, and the counts of systems beyond each contour are all reported as single numbers, even though the peculiar-velocity errors are large enough to move systems across thresholds (e.g., MAXI J1305-704 at 71+41/-28 km/s reaches its V68% = 110 km/s only in the upper error tail, and V404 Cyg at the 60th percentile is eight points inside the 68% contour). The binomial probabilities in §4.1 treat the counts as exact. Please propagate the Monte Carlo samples through the percentile-rank calculation and report, for example, the probability that each system is classified as kicked, and recast the 'at least 4' and 'half' statements with these uncertainties.","section":"§3.1, Table 2, §4.1"}],"minor_comments":[{"comment":"The quoted probabilities (0.098, 0.021, 0.099) are exact-count binomial probabilities, but the appropriate statistics for 'observing six systems beyond the 68% contour' (etc.) are the tail probabilities, which are approximately 0.18 (≥6 of 12 at p=0.32), 0.025 (≥4 of 12 at p=0.10), and 0.12 (≥2 of 12 at p=0.05); the 68% count is therefore not significant at the 90% level, and the text should either state tail probabilities or explicitly say 'exactly'.","section":"Abstract and §4.1"},{"comment":"The statement that the four systems above the 90th percentile 'were born with kicks of ≳100 km/s' overstates what the adopted thresholds imply: the relevant V90% values for the four systems are 142, 108, 98, and 77 km/s (Gaia BH1), so the ≳100 km/s inference is only directly supported for Swift J1727.8-162 and GRO J1655-40; please tie the statement to the individual V90% thresholds.","section":"Abstract"},{"comment":"'Blauuw kick' is a typo for 'Blaauw kick' (Blaauw 1961).","section":"Section 1"},{"comment":"Two distinct Miller-Jones et al. (2009) papers (ApJL 706, L230 and MNRAS 394, 1440) are cited without year-letter disambiguation; please introduce 2009a/2009b in the text and reference list.","section":"References and §4.2.8"},{"comment":"The analysis depends entirely on Gaia DR3 queries and Monte Carlo sampling, but no code, catalog of the comparison samples, or query parameters are provided; a reproducibility statement or data release would substantially strengthen the paper.","section":"Throughout"},{"comment":"'A minimum natal kick of > 14 km/s' should read 'of 14 km/s' (or '≥14 km/s') to match the V'_kick,min = 14 entry in Table 3.","section":"§4.2.2"},{"comment":"The orbits are integrated for only 1 Gyr although §4.4 describes LMXBs as several Gyr old; a sentence justifying the integration time as illustrative of the present-day orbital structure would prevent confusion.","section":"Figure 4 and §4.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is close to publishable after a revision that (i) aligns Eqs. (4)-(5) with Table 2 (the table is almost certainly the correct version, as I verified V_kick,min = V_pec,local - V68% row by row), (ii) adds a test of the age-matching assumption, and (iii) propagates uncertainties into the classifications. The central mixed-origin claim is substantially supported by independent literature constraints (V404 Cyg, VFTS 243) and is robust to the comparison-sample construction on the strong-kick side, which lowers the risk of this revision. For the editor's awareness, the binomial significance of the 6/12 count beyond the 68% contour is modest (tail probability ≈0.18), so the paper's strongest quantitative statements are the >90th/95th percentile counts and the independent weak-kick cases rather than the 68% count; the authors should be encouraged to hedge that specific claim. The citation pattern is appropriate, and the paper builds on rather than hides prior work, including the authors' own Gaia BH1/BH2 papers."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely useful: instead of treating peculiar velocity as a direct kick proxy, it compares each BH binary to the actual Gaia DR3 velocity distribution of stars in its vicinity. That's a real improvement over Atri, Zhao, and earlier work, and it produces a cleaner qualitative picture: at least a few BHs are kinematically hotter than 90-95% of their local populations (Swift J1727.8-162, GRO J1655-40), while a few others sit comfortably inside the local dispersion. The authors also deserve credit for flagging their own main assumption—that the comparison stars and the BH binary trace the same birth population—and for the careful discussion of individual systems against the prior literature.\n\nThe soft spots are real but fixable. First, Eq. (4) defines V as the full 3D Galactocentric speed (including V_phi), yet the tabulated Vkick,min values match Vpec,local - V68%, not V - V68%. That's a clear internal inconsistency; the text and table cannot both be right. The authors clearly intended the peculiar velocity, but it needs to be stated and used consistently. Second, the comparison samples are dominated by old red giants. For a young system like Cyg X-1, comparing to a hot, dynamically heated population will push it toward the \"no kick\" side—the paper admits this in Section 4.2.10 but does not test it. The disk-crossing test in Section 4.4 changes the sample location, not the age mix. So the claim that six BHs disfavor kicks of ≳50 km/s is weaker than it looks. The strong-kick end is more robust, since those systems are beyond 90-95% even of an older, hotter population. Third, there are no propagated uncertainties on the V68% thresholds or on the resulting kick/no-kick classification; the percentiles are treated as exact.\n\nThat said, the central qualitative conclusion—that some BHs form with strong kicks and some with essentially none—does not rest solely on the kinematic comparison. V404 Cyg's triple and VFTS 243's circular orbit provide independent, model-robust evidence for very weak kicks. So even if the age-matching problem were severe, the mixed-origin claim would survive in some form. This is not a paper with a load-bearing flaw; it's a paper with a presentation error and a robustness gap.\n\nI'd send it to peer review. The method is worth publishing, and the flaws are correctable. A serious referee should ask for the equation fix, a quantitative test of the age-matching assumption (e.g., a young-star subsample or a kinematic model), and error bars on the classifications. It's the kind of paper that will be cited for the method and the revised individual kick estimates, so it deserves careful refereeing rather than a desk rejection.","headline":"A useful and mostly sensible re-analysis of BH natal kicks whose qualitative mixed-origin conclusion is probably right, but the quantitative thresholds have an equation/table mismatch and a real age-matching weakness that need fixing before publication.","tokens_in":28388,"tokens_out":2498,"would_cite":true,"duration_ms":26504,"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":"Black holes are born with a mix of strong and weak natal kicks, and the observed split matches a mixed-origin scenario of direct collapse and supernova kicks.","keywords":["black hole natal kicks","stellar kinematics","Gaia DR3","X-ray binaries","Toomre diagrams","Galactic orbits","direct collapse black holes","supernova kicks"],"falsifier":"Measure the ages (by asteroseismology, spectroscopy, or Gaia colors) of the local comparison stars around the six systems classified as \"no kick\" and compare them with the ages of the binaries' companions; if the comparison stars turn out to be systematically much older, then a black hole born in a younger, colder population would look unkicked next to elderly heated giants, and the no-kick classifications would collapse. If age-matched samples reproduce the same split, that would confirm the result.","tokens_in":27297,"feed_emoji":"🕳️","tokens_out":6784,"duration_ms":59056,"temperature":0.7,"pith_summary":"Using Gaia DR3 astrometry and radial velocities, this paper tries to establish whether stellar-mass black holes receive a \"natal kick\" when they form, and how big such kicks are. For 12 black-hole binaries with well-measured distances, the authors compare each black hole's present-day space velocity and Galactic orbit with the velocity dispersion of its local stellar population, rather than treating its raw peculiar velocity as the kick. They find that half the systems are kinematically hotter than at least 68% of their local stars, four systems are hotter than 90% of local stars (implying kicks $\\gtrsim 100$ km s$^{-1}$), and six systems look kinematically ordinary (disfavoring kicks $\\gtrsim 50$ km s$^{-1}$). Two systems with independent orbital evidence, V404 Cyg and VFTS 243, require kicks $\\lesssim 10$ km s$^{-1}$. The conclusion is that black holes do not all form gently: some are born with strong kicks and some with almost none, a split that matters for how black-hole binaries form and merge.","feed_headline":"Half of black hole binaries show natal kicks, six show none","feed_subtitle":"Gaia DR3 data put four BHs beyond 90% of local stars, pointing to two birth channels.","key_machinery":"The argument is carried by the Toomre diagram: a plot of the non-azimuthal velocity $\\sqrt{V_R^2 + V_z^2}$ against the azimuthal velocity $V_\\phi$ for the black-hole binary and for several hundred nearby Gaia DR3 stars. Around the comparison sample's median azimuthal velocity the authors draw semicircles that enclose 68%, 90%, and 95% of the local stars; these contours represent how hot the local stellar population is, i.e. how much of a black hole's motion could be ordinary dynamical heating rather than a kick. The minimum kick is the distance from the black hole's velocity to the 68% contour, $V_{\\rm kick,min} = \\max(0,\\, V - V_{68\\%})$, with $V=\\sqrt{V_R^2+V_z^2+V_\\phi^2}$. Galactic orbits integrated backward in time for 1 Gyr serve as a visual cross-check: systems classified as kicked have puffier, more eccentric orbits, while unkicked systems stay near circular orbits close to the disk midplane.","core_discovery":"The central claim is that a black hole's raw present-day velocity is a biased measure of its natal kick, because most black-hole binaries are old and their orbits have been dynamically heated by the Galaxy. Once a local comparison sample from Gaia DR3 is used to define 68%, 90%, and 95% velocity-dispersion contours on a Toomre diagram, the minimum kick is $V_{\\rm kick,min} = \\max(0,\\, V - V_{68\\%})$. By this metric, half of the 12 disk systems show at least weak kick evidence, one-third exceed the 90th percentile of their local population, and two systems (Swift J1727.8-162 and GRO J1655-40) exceed the 95th percentile, implying birth kicks $\\gtrsim 100$ km s$^{-1}$. The other half are consistent with their local populations, so kicks $\\gtrsim 50$ km s$^{-1}$ are disfavored for them. The paper therefore argues that it would be wrong to conclude that most black holes form with weak kicks merely because two well-studied systems did, and that the data are consistent with a mixed origin: direct collapse with weak kicks for some black holes, supernova kicks for others.","pith_inferences":["Because binaries that survive strong kicks are preferentially destroyed, the strong-kick cases in this sample likely underrepresent the true fraction of strongly kicked black holes in the field, so the field population may be even more kick-heavy than the 50% found here.","Applying the same local-dispersion comparison to neutron-star binaries would give a uniform, apples-to-apples measurement of neutron-star versus black-hole kicks; the paper's appendix suggests neutron-star binaries are more dispersed, but not dramatically.","The main caveat is testable: age-dating the comparison stars would let future work replace the implicit same-population assumption with an explicit, age-matched baseline.","With more astrometric black-hole binaries from future Gaia releases, the percentile method could map the kick distribution's shape and distinguish a true bimodal distribution from a broad unimodal one."],"forward_implications":["Minimum kicks inferred for most systems are smaller than the raw peculiar velocities quoted in earlier work, so treating present-day velocity alone as the kick overestimates black-hole kicks in old populations.","At least four black holes (Swift J1727.8-162, GRO J1655-40, GRS 1124-684, and Gaia BH1) are kinematically hotter than 90% of their local stars, consistent with birth kicks of roughly $\\gtrsim 100$ km s$^{-1}$.","Six black holes in the sample have kinematics consistent with their local populations, disfavoring kicks $\\gtrsim 50$ km s$^{-1}$; V404 Cyg and VFTS 243 independently require kicks of only $\\lesssim 5$-$10$ km s$^{-1}$.","Binary population synthesis and gravitational-wave merger-rate models that assume all black holes form with negligible kicks are missing a real population of strongly kicked black holes.","The sample as a whole is consistent with a mixed-origin scenario, some black holes forming by direct collapse with weak kicks and others in supernovae with strong kicks, though the uncertainties are too large to prove the kick distribution is truly bimodal."],"supporting_citations":[{"why":"Supplies the DR3 proper motions and radial velocities used for both the black-hole binaries and the local comparison samples.","marker":"Gaia Collaboration et al. (2023)"},{"why":"Provides the geometric distances with direction-dependent priors used for comparison-sample stars and some targets.","marker":"Bailer-Jones et al. (2021)"},{"why":"Establishes the potential-kick-velocity method that this paper compares against and which ignores the local velocity dispersion.","marker":"Atri et al. (2019)"},{"why":"Supplies the recent population-wide black-hole kick inference and the neutron-star binary catalog used for comparison.","marker":"Zhao et al. (2023)"},{"why":"Measured Swift J1727.8-162's large peculiar velocity, the strongest single-system kick case in the sample.","marker":"Mata Sanchez et al. (2024)"},{"why":"Shows V404 Cyg is a hierarchical triple whose survival requires a natal kick of only about 5 km s$^{-1}$.","marker":"Burdge et al. (2024)"},{"why":"Uses VFTS 243's near-circular, non-tidally-synchronized orbit to constrain its natal kick to about 10 km s$^{-1}$.","marker":"Vigna-Gómez et al. (2024)"},{"why":"Provides the discovery and population-synthesis natal-kick estimate for Gaia BH1.","marker":"El-Badry et al. (2023a)"},{"why":"Gives the BlackCAT catalog of dynamically confirmed black-hole X-ray transients from which the starting sample is drawn.","marker":"Corral-Santana et al. (2016)"}],"fun_headline_variants":["Half of black holes show natal kicks, six do not (Gaia DR3)","Black holes: strong natal kicks for some, none for others","Gaia data: some black holes get supernova kicks, others none","Mixed black hole births: strong kicks for a third, weak for rest","Natal kick census: half of black holes show evidence, six do not"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the stars found today within 500 pc of each black hole represent the stellar population in which the black hole's progenitor was born, so that the spread of their velocities is the right baseline for \"no kick.\"","fun_headline_variants_meta":{"raw":{"variants":["Half of black holes show natal kicks, six do not (Gaia DR3)","Black holes: strong natal kicks for some, none for others","Gaia data: some black holes get supernova kicks, others none","Mixed black hole births: strong kicks for a third, weak for rest","Natal kick census: half of black holes show evidence, six do not"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000364,"raw_usage":{"total_tokens":2055,"prompt_tokens":1133,"completion_tokens":922,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":749,"completion_tokens_details":{"reasoning_tokens":825}},"tokens_in":749,"tokens_out":922,"duration_ms":8422,"temperature":1.0,"reasoning_tokens":825,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:48:25.550840+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the ages (by asteroseismology, spectroscopy, or Gaia colors) of the local comparison stars around the six systems classified as \"no kick\" and compare them with the ages of the binaries' companions; if the comparison stars turn out to be systematically much older, then a black hole born in a younger, colder population would look unkicked next to elderly heated giants, and the no-kick classifications would collapse. If age-matched samples reproduce the same split, that would confirm the result.","supporting_citations":[{"cited_title":"The black hole low mass X-ray binary V404 Cygni is part of a wide hierarchical triple, and formed without a kick","cited_arxiv_id":"2404.03719","evidence_quote":"Shows V404 Cyg is a hierarchical triple whose survival requires a natal kick of only about 5 km s$^{-1}$."}],"review_version":1}