{"id":"10a161d2-f6a8-46f4-9fc5-01f0a03b14dc","arxiv_id":"2607.07068","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"A Gaia DR3-based census of 105,971 OB stars within 2 kpc maps local Galactic structure and identifies over 4,200 core-collapse supernova or black hole progenitor candidates.","lead":"Astronomers cataloged 105,971 hot, massive OB-type stars within 2 kpc of the Sun using Gaia DR3 data, mapping the Milky Way's spiral arms and star-forming regions. The resulting 3D map reveals the structure of our galactic neighborhood and identifies over 4,200 future supernova or black hole progenitors, finding none will endanger Earth within the next million years.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The claim that more BH progenitors than ccSNe will collapse/explode within 1 Myr rests on SED-fitted masses whose uncertainties (~0.14 dex) far exceed the gaps between adjacent fate thresholds in Table 2 (as small as 1.12 M⊙), making the statistical classification unreliable.","rationale":"The reader correctly identified the progenitor classification as the weakest area and the single-star assumption as a real systematic. I partially agree: the single-star assumption matters, but the more fundamental issue is that the SED-fitted mass precision (~0.14 dex) is insufficient to resolve the fate thresholds in Table 2, some of which are separated by only ~1 M⊙. This affects not just individual classifications but the statistical claim about relative BH vs ccSN numbers. However, the reader's verdict of CONDITIONAL with MODERATE confidence already captures the appropriate level of caution. The paper's primary contributions — the catalog itself, the spatial mapping, and the structural analysis — are well-validated (Section 3 shows reasonable agreement with spectroscopic surveys) and do not depend on the progenitor classification precision. The ccSN/BH progenitor identification is presented as a secondary result with explicit caveats (Appendix F, Section 6.3.3). The authors themselves state they will address this with synthetic modeling in Paper II. Therefore, the verdict should remain CONDITIONAL: the catalog and structural results are sound, but the progenitor classification results — particularly the 'more BH than ccSN' claim — should be treated as preliminary pending proper posterior propagation. No change to the reader's verdict is needed.","tokens_in":36644,"tokens_out":3390,"duration_ms":102581,"concrete_test":"For each of the 4231 progenitor candidates, draw 1000 samples from the SED posterior on M_ZAMS (and fractional age), classify each draw using Table 2 and the Ekström et al. (2012) lifetime models, and compute the fraction of Monte Carlo realizations in which 'more BH progenitors collapse within 1 Myr than ccSNe explode.' If this fraction is below ~90%, the headline claim is not robust to the stated mass uncertainties.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most novel claim — that more direct-collapse BH progenitors will collapse within 1 Myr than ccSNe will explode (Section 6.3.1) — depends on correctly sorting stars into the fate categories defined in Table 2. These categories are separated by ZAMS mass thresholds including a gap of only 1.12 M⊙ between 22.20 and 23.32 M⊙ (failed SN vs neutrino-driven SN). The SED fitting yields log(M/M⊙) with 1σ ≈ 0.14 dex (Section 2.3), corresponding to ~32% uncertainty in linear mass — roughly 7 M⊙ at 22 M⊙. This uncertainty far exceeds the 1.12 M⊙ gap between adjacent fate categories, meaning stars near these boundaries are classified essentially by noise. The authors acknowledge in Appendix F that individual error bars are 'too large to faithfully predict the final fates,' but the statistical claim about relative BH vs ccSN numbers is subject to the same limitation: the classification uses median SED-fitted masses (Section 6.2), collapsing the full posterior to a point estimate. If the mass posteriors were propagated through the Table 2 thresholds, the apparent excess of BH progenitors could dissolve into classification noise. The single-star assumption (the reader's concern) compounds this, but the mass precision relative to threshold spacing is the more fundamental issue — even for truly single stars, the SED masses cannot reliably distinguish a 22 M⊙ star from a 24 M⊙ star, yet the fate model treats these as categorically different outcomes.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This paper presents a catalogue of 105,971 OB-type stars (T_eff > 10,000 K) within 2 kpc of the Sun, identified via astro-photometric SED fitting of Gaia DR3 data combined with 2MASS, IGAPS, and VPHAS+ photometry. The authors validate their catalogue against spectroscopic surveys (APOGEE, LAMOST, Gaia ESP-HS, SHBoost), demonstrate high completeness (>95% across magnitudes), and map the spatial distribution of young stellar populations across the Galactic thin disk. They identify large-scale structures (Sagittarius-Carina arm, Cepheus Spur, Giant Oval Cavity), compare OB star overdensities with young open clusters from Hunt & Reffert (2024), and use a statistical supernova model (Maltsev et al. 2025) to classify 3,998 ccSN progenitor and 233 direct-collapse BH progenitor candidates. The most novel claim is that more BH progenitors are expected to collapse within 1 Myr than ccSNe to explode, interpreted as evidence for a recent massive star formation burst.","tokens_in":37116,"tokens_out":1803,"duration_ms":222122,"significance":"The catalogue represents a substantial expansion over the authors' previous 1 kpc census (Q25), covering four times the area and providing a valuable community resource for spectroscopic follow-up (WEAVE, 4MOST). The cross-validation against four independent spectroscopic surveys is commendable and lends credibility to the derived effective temperatures. The spatial mapping of OB stars and their correlation with young open clusters provides a useful complementary view to existing structure tracers. The application of a recently developed, observationally-constrained supernova model (M25) to a large stellar census is a genuine attempt to move beyond simple mass-threshold progenitor classifications. The identification of specific ccSN/BH progenitor candidates, while subject to the uncertainties discussed below, provides a concrete target list for follow-up. The authors are transparent about limitations in Appendix F and Section 6.3.3.","major_comments":[{"comment":"Section 6.3.1 and the abstract claim that 'more BH progenitors to collapse within the next 1 Myr than ccSN to explode' is indicative of a recent massive star formation burst. This claim is load-bearing for the paper's most novel conclusion, but it rests on classifying stars into fate categories (Table 2) whose thresholds are separated by as little as 1.12 M_sun (between 22.20 and 23.32 M_sun). The SED-fitted log(M/M_sun) has a 1-sigma of ~0.14 dex (Section 2.3), corresponding to ~32% uncertainty in linear mass (~7 M_sun at 22 M_sun). This uncertainty far exceeds the spacing between adjacent fate thresholds. The classification in Section 6.2 uses median SED-fitted ZAMS masses, collapsing the full posterior to a point estimate. If the mass posteriors were propagated through the Table 2 thresholds, the apparent excess of BH progenitors could be an artifact of classification noise near the 8","section":null},{"comment":"Section 6.3.3 acknowledges that *P Cyg is a likely binary donor whose classification changes when binary effects are considered, and the SED fitter assumes all stars are single (Section 3.2.2). Given that a large fraction of massive OB stars are in binary systems, this assumption systematically affects derived ZAMS masses, ages, and therefore both the individual progenitor classifications and the statistical distribution of waiting times. The paper should quantify or at least bound the impact of binarity on the relative BH/ccSN numbers. Without this, the central claim in Section 6.3.1 remains uncertain. A Monte Carlo test injecting a realistic binary fraction and mass-ratio distribution into the SED fitting pipeline, then re-classifying fates, would suffice to demonstrate robustness.","section":null}],"minor_comments":[{"comment":"Section 2.1, Eq. (1): the threshold M_G < 1.5 mag is described as 'more liberal' than the A0V value of 1 mag from Pecaut & Mamajek (2013). It would help to state explicitly what contamination fraction this liberal threshold introduces, or reference the contamination analysis already performed in Q25.","section":null},{"comment":"Section 3.2.1: the comparison with APOGEE shows a tendency to overestimate temperatures, opposite to the other surveys. The authors attribute this to APOGEE's 20,000 K model grid limit. It would strengthen the paper to show a histogram or subsample restricted to T_eff < 20,000 K to confirm that the offset vanishes there.","section":null},{"comment":"Section 4.1, right panel of Fig. 5: the overdensity parameter Delta_Sigma is defined in Appendix D, but the bandwidth choices (h_local = 100 pc, h_mean = 500 pc) are stated only in the appendix. These should be mentioned in the main text caption or Section 4.1 for self-contained reading.","section":null},{"comment":"Table 3: several entries have extremely asymmetric error bars on distance (e.g., chi Oph: 705 +1975/-351 pc). These large upper bounds suggest the SED-fitted distance posterior is poorly constrained for some sources. A note flagging which entries have distance uncertainties exceeding the median would help readers assess the reliability of individual progenitor candidates.","section":null},{"comment":"Section 5: the claim that ~10% of OB stars are found in clustered environments is compared to the ~38% of OB association members that are also cluster members within 1 kpc. The jump from 1 to 2 kpc changes the surface area by a factor of 4, so a direct comparison of these fractions without accounting for volume-dependent completeness may be misleading. A brief clarification would help.","section":null},{"comment":"Figure 12 caption: the text mentions '10 ccSN progenitor candidates we expect to explode within less than 1 Myr (as blue stars)' but Table 3 lists 10 entries above the dashed line, some with tau_cc > 0.8 Myr. The caption should clarify that these are the 10 with the shortest median waiting times, not all with tau_cc < 1 Myr.","section":null},{"comment":"Appendix E: the comparison with alternative spiral arm models is useful but somewhat cursory. A quantitative metric (e.g., cross-correlation between OB star density and model arm positions) would strengthen the comparison beyond visual inspection.","section":null},{"comment":"Section 6.1.2: the minimum M_ZAMS thresholds of 8.55 and 8.85 M_sun for electron-capture and neutrino-driven SNe are derived from specific stellar evolution models (Temaj et al. 2024). A brief note on how sensitive these thresholds are to the choice of model would provide context for the reader.","section":null},{"comment":"Typo in Section 6.3: 'pens within less than 100 years' should be 'happens within less than 100 years'.","section":null},{"comment":"Section 3.2.2: the median offset of 0.06 dex in log(M/M_sun) relative to HR24 is attributed to methodological differences. Given that this offset is comparable to the 1-sigma scatter (~0.14 dex), a brief discussion of whether a simple additive correction would improve agreement would be useful.","section":null}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid piece of work with a well-constructed catalogue, but the most novel claim (more BH than ccSN within 1 Myr) is not yet sufficiently supported. The mass uncertainty issue is the more fundamental concern: even if binarity were fully accounted for, the SED masses cannot distinguish stars on either side of the narrow fate thresholds in Table 2. The authors already acknowledge this in Appendix F for individual sources, but the statistical claim in Section 6.3.1 is subject to the same limitation and needs explicit propagation of uncertainties. If the authors can show via Monte Carlo that the BH excess persists after propagating mass posteriors through the fate thresholds, the claim stands; if not, it should be softened. The catalogue itself and the structural mapping are publishable regardless."},"author_rebuttal":null,"desk_editor":{"model":"glm-5.2","letter":"This paper extends the authors' own 1 kpc OB star census (Q25) out to 2 kpc, yielding 105,971 stars with SED-fitted physical parameters. The catalog itself is the main deliverable — it's well-validated against spectroscopic samples (APOGEE, LAMOST, Gaia ESP-HS), the completeness analysis is thorough, and the spatial maps of local Galactic structure (Sagittarius-Carina arm, Perseus Gap, Cepheus Spur) are a genuine improvement over previous work. The crossmatching with the Hunt & Reffert cluster catalog to quantify clustered vs. field OB star fractions is a nice touch. As a target list for WEAVE and 4MOST, this is useful and timely work. The authors are also commendably honest about limitations throughout — they flag the mass offset relative to HR24, the single-star assumption, and the large error bars on individual progenitor predictions. The stress-test concern about mass precision relative to the fate thresholds in Table 2 is the real issue. The SED fitter yields log(M/M☉) with 1σ ≈ 0.14 dex, which is roughly 7 M☉ at 22 M☉. The gap between the failed-SN and neutrino-driven-SN categories is 1.12 M☉. So stars near these boundaries are being classified by noise. The authors acknowledge this in Appendix F for individual sources, but the statistical claim — more BH progenitors collapsing within 1 Myr than ccSNe exploding — is subject to the same limitation. If the full mass posteriors were propagated through the Table 2 thresholds rather than collapsed to medians, the apparent excess could easily dissolve. The single-star assumption compounds this, but the mass precision problem is more fundamental. The catalog and structure mapping deserve publication. The progenitor classification results should be framed more cautiously — as illustrative rather than statistically robust predictions. The paper is appropriate for a serious referee in MNRAS or similar.","headline":"Solid OB star catalog extending to 2 kpc; the BH-vs-ccSN excess claim is the soft spot","tokens_in":37802,"tokens_out":473,"would_cite":true,"duration_ms":82546,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"105,971 OB stars mapped within 2 kpc reveal local spiral arms and future supernova progenitors","keywords":["OB stars","Gaia DR3","Milky Way spiral arms","core-collapse supernova progenitors","black hole formation","star formation rate","stellar catalogues","SED fitting"],"falsifier":"If spectroscopic follow-up of the 4,200+ progenitor candidates reveals that a substantial fraction have ZAMS masses or ages inconsistent with the SED-fitted values—particularly if binary companions are detected that shift stars across the supernova / direct-collapse boundary—then the reported excess of imminent BH collapses over supernovae could be an artefact of the single-star assumption rather than a real signature of recent star formation.","tokens_in":36936,"feed_emoji":"⭐","tokens_out":1369,"duration_ms":115036,"temperature":0.7,"pith_summary":"The paper constructs the largest and most carefully validated catalogue of OB-type stars (effective temperature above 10,000 K) within 2 kiloparsecs of the Sun, using Gaia DR3 astrometry and photometry combined with multi-band survey data fed into a Bayesian spectral-energy-distribution fitter. From the derived stellar parameters—initial mass, age, distance, and extinction—the authors produce a three-dimensional map of young stellar populations across the local Galactic thin disk, showing where star formation is clustered, how it traces spiral arm segments (notably the Sagittarius-Carina arm), and where large underdensities such as the Perseus Gap / Giant Oval Cavity sit. They then apply a statistical supernova model to classify over 4,200 of these stars as core-collapse supernova or direct-collapse black hole progenitor candidates, estimate waiting times to explosion or collapse, and find that more black hole progenitors are expected to collapse within the next 1 million years than supernovae are expected to explode—despite black hole progenitors being far rarer overall. This excess of imminent collapses is interpreted as evidence for a recent, localized burst of massive star formation in the solar neighbourhood.","feed_headline":"105,971 OB stars mapped, more black-hole collapses than supernovae loom nearby","feed_subtitle":"A Gaia-based census of massive stars within 2 kpc reveals spiral arm structure, confirms clustered star formation, and finds an excess of即将坍","key_machinery":"The central machinery is an astro-photometric Bayesian SED fitter that combines Gaia DR3 parallax and photometry with ground-based optical and near-infrared surveys (2MASS, IGAPS, VPHAS+), constrains extinction using a 3D dust map, and fits stellar atmosphere and evolutionary models to derive initial mass, fractional age, distance, and effective temperature for each star. A statistical core-collapse supernova model (M25) then maps the derived zero-age main-sequence masses onto final fates—successful neutrino-driven supernova with neutron star remnant, failed supernova with black hole remnant, or electron-capture supernova—using metallicity-dependent carbon-oxygen core mass thresholds.","core_discovery":"By fitting physical parameters for 105,971 OB-type stars within 2 kpc and applying a bimodal supernova model, the paper finds that the local Milky Way currently harbours more massive stars on the verge of direct black hole collapse than on the verge of supernova explosion, an imbalance attributed to a recent burst of massive star formation (likely associated with Cyg OB2) rather than to a steady-state rate. The paper also establishes that no OB-type supernova progenitor is expected to explode within 100 parsecs of Earth in the next 1 million years.","pith_inferences":["If the single-star assumption were relaxed, binary mass transfer could strip envelopes and shift progenitors between the supernova and direct-collapse BH channels, potentially altering the reported excess of BH collapses over supernovae. The direction of the effect is not obvious: stripping could push some stars that would have collapsed silently into the explosive channel, or vice versa, dependin","The reported mean waiting time of ~7,800 years between supernova explosions on the 2 kpc scale, combined with the much longer ~15,000-year average inferred from extrapolating the 1 kpc rate over ~400 Myr, suggests that the local Milky Way is currently in a star-formation overdensity phase. If so, the next few thousand years should see a local supernova rate elevated above the long-term Galactic av","The offset between OB star overdensities and young open cluster positions in highly extinguished regions like Cyg OB2 implies that OB star catalogues can probe star formation sites that cluster catalogues miss, making the two tracers complementary rather than redundant for mapping embedded massive star formation."],"forward_implications":["The catalogue of 4,200+ ccSN and BH progenitor candidates provides a prioritized target list for spectroscopic follow-up surveys such as WEAVE-SCIP and 4MOST, which can test the SED-fitted masses and ages against direct spectroscopic measurements.","The finding that more BH collapses than supernovae are expected within 1 Myr, if confirmed, implies that the local core-collapse event rate is time-variable rather than constant, complicating estimates of the Galactic supernova rate that assume steady-state star formation.","The spatial correlation between OB star overdensities and young open clusters, combined with the ~10% membership fraction of OB stars in clusters, constrains the fraction of massive stars born in clustered environments versus the field.","Future Gaia DR4 astrometry (expected to be ~2.25 times more precise) and a potential near-infrared successor (GaiaNIR) would push the catalogue deeper into high-extinction regions toward the Galactic Centre, where current optical surveys are incomplete."],"fun_headline_variants":["More black-hole collapses than supernovae expected in local Milky Way","No OB supernova progenitor within 100 pc set to explode in next 1 Myr","Gaia maps 105,971 OB stars, finds recent massive star-formation burst","Local Milky Way shows excess of imminent BH collapses over supernovae","Nearby OB-star census finds more black-hole collapses than supernovae"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The SED fitting and supernova classification assume every star is a single, non-interacting system. A large fraction of massive OB stars are known to exist in binaries, and binary mass transfer can change a star's final mass, lifetime, and whether it explodes as a supernova or collapses directly into a black hole. The authors acknowledge this limitation explicitly through the case of P Cyg, whose classification flips when binary effects are considered.","fun_headline_variants_meta":{"raw":{"variants":["More black-hole collapses than supernovae expected in local Milky Way","No OB supernova progenitor within 100 pc set to explode in next 1 Myr","Gaia maps 105,971 OB stars, finds recent massive star-formation burst","Local Milky Way shows excess of imminent BH collapses over supernovae","Nearby OB-star census finds more black-hole collapses than supernovae","No supernova within 100 pc in next 1 Myr; BH collapses predominate","Gaia OB census reveals more black-hole than supernova progenitors nearby","105,971 OB stars mapped: recent burst skews collapse outcomes","No supernova within 100 pc expected for next 1 Myr, Gaia census confirms","Excess BH progenitors over supernovae point to recent formation burst"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":1982,"prompt_tokens":669,"completion_tokens":1313,"prompt_tokens_details":null},"tokens_in":669,"tokens_out":1313,"duration_ms":54845,"temperature":1.0,"reasoning_tokens":1041,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T20:36:47.508713+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If spectroscopic follow-up of the 4,200+ progenitor candidates reveals that a substantial fraction have ZAMS masses or ages inconsistent with the SED-fitted values—particularly if binary companions are detected that shift stars across the supernova / direct-collapse boundary—then the reported excess of imminent BH collapses over supernovae could be an artefact of the single-star assumption rather than a real signature of recent star formation.","supporting_citations":[],"review_version":1}