{"id":"2250e989-d51d-4d41-aec9-e4b4143f72c2","arxiv_id":"2502.09426","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"This paper identifies 12 massive prestellar core candidates in 14 protoclusters and estimates prestellar lifetimes of 50 to 240 thousand years, about 10 to 30 times the free-fall time.","lead":"Using ALMA observations of 14 massive protoclusters, this study finds 30 massive cores with no detected outflow, including 12 with masses above 16 solar masses, the best candidates yet for the quiet precursors of high-mass stars. It uses the ratio of these cores to outflow-driving cores to estimate how long the prestellar phase lasts, about 50 to 240 thousand years, suggesting that collapse is delayed by turbulence, magnetic fields, or rotation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Prestellar classification rests entirely on outflow non-detection; the paper's own scenario list admits geometries where a high-mass protostar produces no detectable CO/SiO outflow, which would inflate the lifetime ratios.","rationale":"The reader's weakest assumption matches mine: if outflow non-detection does not imply the absence of a highly accreting high-mass protostar, then the sample is contaminated and the lifetime ratios are overestimated. The paper's careful treatment of projection effects (Sect. 4.2) is commendable, but it does not cover the distinct physical scenarios listed in Sect. 7.8, where a protostar could accrete without ejecting detectable CO/SiO inside the core. The paper honestly flags these possibilities, but the headline conclusion about 10-30 free-fall times depends on their being rare. A methanol maser survey is a clean, decisive test because Class II masers trace high-mass protostars directly and would not be suppressed by the same mechanisms that hide outflows. This concern reinforces the CONDITIONAL verdict without warranting rejection or a change in the reader's assessment.","tokens_in":58537,"tokens_out":4373,"duration_ms":42127,"concrete_test":"Observe the 12 HM PSC candidates with the JVLA in the 6.7 GHz CH3OH maser transition (or ATCA at 12.2 GHz) to an rms of ~1 mJy. Class II methanol masers are a virtually unambiguous tracer of high-mass protostellar accretion, independent of outflow brightness or orientation. If any candidate shows maser emission, reclassify it as protostellar and recompute the prestellar lifetimes and free-fall ratios; if all 12 are maser-free, the prestellar interpretation receives independent confirmation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that high-mass prestellar cores live 10 to 30 free-fall times depends on the 30 PSC candidates (12 above 16 M_sun) being genuinely prestellar, i.e., containing no highly accreting high-mass protostar. The On-Off CO/SiO detection method is the sole discriminator. The false-negative calculation in Sect. 4.2 (1.75% for canonical outflows) covers only geometric projection for gas already entrained above the 6 km/s threshold; it does not address the physically distinct alternatives the paper itself lists in Sect. 7.8. In particular, item 5 (cores mostly empty, dust concentrated in smaller fragments such that jets escape without interacting with the observed core) and item 3 (dead zones suppressing magneto-centrifugal ejection) allow a high-mass protostar to accrete without driving a detectable outflow inside the 2700-au beam. If even 3 of the 12 robust candidates were misclassified, the prestellar-to-protostellar ratio would drop from ~0.44 to ~0.25, reducing the derived lifetimes and free-fall ratios by nearly a factor of two and weakening the conclusion that collapse is slowed by turbulence, magnetic fields, or rotation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript uses the ALMA-IMF survey of 14 massive protoclusters to search for high-mass prestellar cores. From the 580 gravitationally bound cores of Louvet et al. (2024), the authors isolate 141 cores with mass >8 Msun at an assumed 20 K dust temperature and classify each as protostellar or prestellar using a new automated On-Off aperture method that searches for CO(2-1) and SiO(5-4) high-velocity line wings relative to a DCN-based VLSR, complemented by outflow lobe maps. After rejecting cores whose apparent excess emission can be attributed to neighbouring outflows, 42 candidates without outflow remain; adopting temperature-dependent masses between 20 K and PPMAP values reduces this to 30 likely prestellar cores with M>8 Msun, including 12 above 16 Msun. The authors combine prestellar/protostellar count ratios with an adopted 300 kyr high-mass protostellar lifetime, apply corrections for missing UCHII protostars, and consider two scenarios for envelope-mass evolution, deriving prestellar lifetimes of 120-240 kyr for 8-16 Msun cores and 50-100 kyr for 30-55 Msun cores, corresponding to 10-30 free-fall times. They interpret these large ratios as evidence that collapse is slowed by turbulence, magnetic fields, or rotation.","tokens_in":58804,"tokens_out":9375,"duration_ms":90730,"significance":"If the classification holds, this is the first statistically meaningful sample of high-mass prestellar core candidates from a single survey, and the inferred lifetime relative to free-fall is an important observational constraint on high-mass star formation. The paper's strengths are substantial: the Out-Off outflow method is described in enough detail to be reproduced, per-candidate spectra and outflow maps are shown in Appendix E, the noise is estimated per velocity channel (Sect. 3.2), and the method is validated against visual classifications with 95% agreement in W43 (Sect. 4.3). The lifetime derivation is also not circular: it uses count ratios and an externally adopted protostellar lifetime, and the free-fall comparison is a post-hoc benchmark rather than an input. The main risk to significance is that the prestellar classification relies entirely on outflow non-detection, and the manuscript's own Sect. 7.8 lists physical scenarios in which a highly accreting high-mass protostar could hide inside one of the candidates. A quantitative treatment of these contamination channels is needed before the lifetime and 10-30 free-fall claims can be regarded as secure.","major_comments":[{"comment":"The central lifetime and free-fall-ratio claims rest on outflow non-detection being a reliable proxy for the absence of a highly accreting high-mass protostar. Section 4.2 quantifies only the geometric projection false-negative rate (1.75-7.35% for entrained gas already above the 6 km/s threshold), but Section 7.8 itself lists physically distinct failure modes: dead zones suppressing magneto-centrifugal ejection (item 3), cores that are mostly empty with dust concentrated in smaller fragments so that jets escape without interacting with the observed core (item 5), and long accretion outages (item 4). These scenarios would place a strongly accreting object inside a core with no detectable CO/SiO outflow inside the 2700-au beam, and they are not folded into the false-negative budget. The impact is load-bearing: if only a few of the 12 candidates above 16 Msun were misclassified, the prestellar-to-protostellar ratio would drop from roughly 0.44 to about 0.25, reducing the derived lifetimes and free-fall ratios by nearly a factor of two. The manuscript should either place quantitative upper bounds on these contamination channels using independent tracers or literature limits on outflow duty cycle and accretion-ejection coupling, or it should explicitly recast the lifetime and 10-30 free-fall conclusions as upper limits conditional on the assumed outflow-accretion link. As written, the Sect. 7.8 caveat undermines the load-bearing inference rather than resolving it.","section":"Sect. 7.8 (with Sect. 4.2)"},{"comment":"The sample of 30 likely prestellar cores with M>8 Msun and the robust set of 12 above 16 Msun is threshold-sensitive to the assumed dust temperature. Several Table 3 entries have mass ranges that cross the selection cuts: G008.67 #4 (7.8 Msun, [6.6-8.9]), G333.60 #8 (6.7, [5.4-8.0]), G338.93 #16 (6.5, [4.8-8.2]), W43-MM2 #9 (7.8, [6.7-8.8]), W51-E #28 (7.5, [6.7-8.2]), and W51-IRS2 #22 and #24 (7.6-7.7 with lower limits near 6) can fall below 8 Msun at the PPMAP end; W43-MM2 #12 (16.0, [14.0-18.1]) and W51-IRS2 #8/#9 (17.4, [14.0-20.9] and 16.5, [13.0-19.9]) can fall below 16 Msun. The text partly acknowledges this through the [9-14] range in Sect. 6 and through the black/brown temperature-extreme points in Figs. 8-9, but the lifetime estimates in Table 5 and the abstract are quoted as fixed numbers. The authors should state explicitly which temperature choice defines the 12 robust candidates and should propagate the threshold migration into the quoted lifetime ranges, or present the lifetimes as conditional on the adopted temperature set.","section":"Table 3 and Sect. 4.5/5.1"}],"minor_comments":[{"comment":"The text contains a typo: 'Bonnort-Ebert spheres' should read 'Bonnor-Ebert spheres'.","section":"Sect. 1"},{"comment":"The solar-mass symbol appears inconsistently as 'Md', 'M d', and 'M_sun' in the abstract, Table 3, and body text; please unify the notation (preferably as M_sun).","section":"Throughout"},{"comment":"The abstract quotes '10 to 30 free-fall times,' but Sect. 7.5 and Fig. 10 present values up to about 40 for the constant-envelope (clump-fed) scenario; please align the quoted range with the figures or explicitly state which scenario the abstract range refers to.","section":"Abstract and Sect. 7.5"},{"comment":"The caption says the figure displays the 12 most massive PSC candidates, but G333.60 #1A and #1B are combined in a single panel; state this explicitly so the reader knows that 12 objects are shown in 11 panels.","section":"Fig. 5 caption"},{"comment":"The meaning of the 'Mean masses' and the adopted central mass should be defined in the table header or footnote; the current text refers only to Sect. 4.5, and the relation between the adopted mass and the lower/upper limits is not immediately clear from the table itself.","section":"Table 3"},{"comment":"The calculation behind the 1.75% and 7.35% geometric false-negative fractions would be easier to verify if the assumed distributions of opening angle, inclination, and maximum outflow velocity were given in a short appendix or explicit formula.","section":"Sect. 4.2"}],"recommendation":"major_revision","confidential_remarks":"I recommend major revision rather than rejection. The candidate list, per-core outflow spectra, and automated method are valuable and largely reproducible, and the central inference is defensible in principle. The blocking issue is interpretational: the prestellar classification rests on outflow non-detection, and the paper's own alternative-scenario list in Sect. 7.8 shows that the false-negative budget is not closed. A quantitative contamination analysis, or a clearly qualified statement of the lifetimes as upper limits, would put the manuscript within reach of acceptance. There is no novelty or scope concern; the dependence on companion ALMA-IMF papers for masses and noise is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: this is the paper that finally gives high-mass prestellar cores a real sample. Thirty candidates above 8 Msun, twelve above 16 Msun, drawn homogeneously from 14 protoclusters. That is genuinely new—previous work had at most a handful of candidates, often with conflicting mass calibrations. The classification is carefully documented: every candidate has CO and SiO spectra, outflow maps, and DCN VLSR fits in the appendices, and the automated On-Off method is validated against the Nony et al. visual classifications with 95% agreement in W43. The 12 robust candidates, mostly in crowded cluster centers, are plausible targets for follow-up. I consider the sample construction and the rarity conclusion solid.\n\nThe soft spot is the one the paper itself identifies in Sect. 7.8: prestellar status is defined by absence of a CO/SiO outflow. That works if every high-mass protostar drives a detectable outflow in its core, which is the standard assumption, but the paper lists physically distinct ways it could fail—dead zones suppressing magneto-centrifugal ejection, or cores where the dust is concentrated in small fragments so jets escape without entraining enough gas. The false-negative estimate in Sect. 4.2 covers only geometric projection for gas already above the 6 km/s threshold; it does not quantify these alternative channels. If even three of the twelve are misclassified, the prestellar-to-protostellar ratio drops from about 0.44 to 0.25, cutting the derived lifetimes and the free-fall ratios by nearly a factor of two. That would weaken but not eliminate the central conclusion: even with a decreasing protostellar lifetime, the authors show ratios of 7-10, still well above the free-fall expectation.\n\nThe lifetime numbers themselves depend on an adopted 300 kyr protostellar lifetime and on which mass-evolution scenario you pick (constant vs decreasing envelope). That is not circular—the prestellar lifetimes come from count ratios times an external timescale—but it means the absolute numbers carry an uncertainty the paper does not fully propagate. The dust temperature issue is handled honestly with mass ranges that straddle the 8 and 16 Msun thresholds.\n\nBottom line: this is a careful, honest, and useful paper, for anyone working on the initial conditions of high-mass star formation—observers who want a catalog of targets for JWST or high-resolution ALMA follow-up, and theorists who need empirical lifetimes to compare against simulations. The catalog and the 12 candidates are the main contribution; the lifetime ratios are suggestive rather than definitive. It deserves a serious referee, who should push on the misclassification sensitivity and ask for a quantitative treatment of the Sect. 7.8 scenarios. I would cite it for the sample and the method.","headline":"First real sample of high-mass prestellar core candidates; the lifetime ratios are suggestive but rest on an assumption the authors themselves flag.","tokens_in":59506,"tokens_out":3116,"would_cite":true,"duration_ms":29698,"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":"The paper identifies 30 outflow-free cores above 8 $M_\\odot$ as massive prestellar candidates — 12 above 16 $M_\\odot$ — and estimates their lifetime at 50–240 kyr, or 10–30 free-fall times.","keywords":["high-mass star formation","prestellar cores","protostellar cores","molecular outflows","core lifetimes","free-fall time","ALMA-IMF survey","galactic protoclusters"],"falsifier":"A sub-arcsecond search for embedded heating in the 12 cores above 16 $M_\\odot$ — mid-infrared continuum or hot-core molecular emission at high spatial resolution — would settle the central claim: several detections of compact warm dust or hot-core chemistry would mean the outflow-free cores are not all starless, and the reported prestellar lifetimes would drop toward the free-fall times. The temporal version is equally concrete: re-observing the same fields after a few years and finding newly appeared outflow wings from any of the 12 candidates would falsify their prestellar status directly.","tokens_in":58362,"feed_emoji":"⭐","tokens_out":15598,"duration_ms":122423,"temperature":0.7,"pith_summary":"The paper sets out to find the direct precursors of high-mass stars: compact, massive dust cores that have not yet ignited a protostar inside them. Using the ALMA-IMF survey's roughly 580 cores in 14 massive protoclusters, it classifies the 141 most massive cores by whether they drive CO or SiO outflows, the expected signature of active accretion. Cores with no outflow in either tracer are treated as prestellar, and the paper reports 30 such cores above 8 $M_\\odot$, including 12 above 16 $M_\\odot$ that it proposes as the best current candidates for high-mass star precursors. Counting these against outflow-driving protostellar cores and assuming a 300 kyr protostellar lifetime, it derives prestellar lifetimes of roughly 120–240 kyr for the 8–16 $M_\\odot$ range and 50–100 kyr for the 30–55 $M_\\odot$ range. Because these timescales are 10–30 times longer than the cores' free-fall times, the paper concludes that collapse is slowed by non-thermal support — turbulence, magnetic fields, or rotation — turning the high-mass prestellar phase into a measurable stage that formation models must reproduce.","feed_headline":"Found: 12 starless cores massive enough to birth high-mass stars","feed_subtitle":"Their 50–240 kyr lifetimes run 10–30× the free-fall time, implying turbulence or fields stall collapse.","key_machinery":"The load-bearing tool is the On-Off spectrum method: for each of the 141 cores above 8 $M_\\odot$, a core-averaged (On) spectrum is taken inside the continuum ellipse and an annular (Off) background between 2.5 and 3.5 times the ellipse FWHM is subtracted, isolating high-velocity line-wing emission that belongs to the core itself rather than its surroundings. Emission beyond $\\pm 6$ km s$^{-1}$ from the core's DCN(3–2) velocity, reaching 5$\\sigma$ in either CO or SiO, marks a protostellar outflow; per-channel noise maps built from 150 random On-Off positions set the significance, and moment-zero maps of blue- and red-shifted lobes provide a spatial cross-check that rejects contamination by neighbouring outflows. The same counts feed the lifetime estimate: the prestellar-to-protostellar number ratio is multiplied by an assumed protostellar lifetime of 300 $\\pm$ 100 kyr, the protostar count is corrected upward for expected ionising (UCHII) phases using the protostellar-evolution models described in Appendix D, and mass bins are built under two scenarios for envelope-mass evolution, a constant-mass clump-fed view and a decreasing-mass core-fed view.","core_discovery":"The central claim is that a homogeneous, statistically useful sample of high-mass prestellar cores can be pulled out of the ALMA-IMF survey by systematic outflow detection: 30 cores above 8 $M_\\odot$ show no CO(2–1) or SiO(5–4) high-velocity emission, and 12 of them exceed 16 $M_\\odot$, with masses up to about 54 $M_\\odot$, mostly in the crowded central clumps of the protoclusters. The paper defines a high-mass prestellar core as a compact core (deconvolved size below ~5000 au) massive enough to form a high-mass star that is not yet hosting any highly accreting protostar, and it argues that any strongly accreting high-mass protostar must drive a detectable outflow inside its own core. Comparing the 30 starless cores with 52 outflow-driving protostellar cores above 8 $M_\\odot$, and adopting a 300 kyr protostellar lifetime with corrections for ionising UCHII protostars and for two envelope-evolution scenarios, the paper obtains prestellar lifetimes of 120–240 kyr (8–16 $M_\\odot$) and 50–100 kyr (30–55 $M_\\odot$). The resulting ratios of lifetime to free-fall time, 10 to 30, are the paper's main quantitative result: they imply that the collapse of massive cores is delayed by non-thermal support of turbulent, magnetic, or rotational origin.","pith_inferences":["Because the paper's own projection analysis allows up to ~7% of outflow detections to be missed, roughly one of the 12 high-mass candidates could still be protostellar; any such misclassification shortens the reported lifetimes in proportion, and this contamination could be measured directly by the high-resolution mid-infrared follow-up the paper proposes.","The 10–30 free-fall-time ratio is a quantitative target for magnetized, turbulent collapse simulations: models should reproduce massive cores that persist for many free-fall times without forming any embedded protostar, and can be checked against the reported DCN line widths of 1–2 km s$^{-1}$ with virial parameters near unity.","Whether the clump-fed or core-fed scenario applies could be settled within the same sample by comparing envelope masses of young versus evolved protostellar cores: replenished envelopes would favour the clump-fed view and the longer lifetimes, while steadily depleted envelopes favour the core-fed view and the shorter values.","The On-Off aperture method is resolution-portable, so applying the same outflow-based classification to surveys at other physical scales or evolutionary stages could map how the prestellar-to-protostellar ratio changes across the high-mass regime."],"forward_implications":["The 12 cores above 16 $M_\\odot$, with masses up to ~54 $M_\\odot$, provide a concrete target list for studying the initial conditions of high-mass star formation in detail with follow-up line and continuum observations.","High-mass prestellar cores are genuinely rare — about 2% of all ALMA-IMF cores — but they exist up to ~50 $M_\\odot$, and they concentrate in the crowded central clumps where high-mass stars are expected to form.","The prestellar phase lasts 50–240 kyr, roughly an order of magnitude shorter than the ~1.2 Myr lifetime of low-mass prestellar cores, so high-mass reservoirs must accumulate gas quickly, consistent with inflow and converging-flow pictures rather than slow quasi-static contraction.","Ratios of lifetime to free-fall time of 10–30 require non-thermal support; the paper notes that ambipolar-diffusion times are about 30 free-fall times, in the right range to explain the longest-lived cores.","The lifetime estimates bracket the true value between the two mass-evolution scenarios: 240 → 120 kyr for 8–16 $M_\\odot$ and 100 → 50 kyr for 30–55 $M_\\odot$ for constant versus decreasing protostellar envelopes."],"supporting_citations":[{"why":"Supplies the catalog of ~580 gravitationally bound cores at a common 2700 au resolution that the outflow classification and lifetime statistics are built on.","marker":"Louvet et al. (2024)"},{"why":"Provides the On-Off differential-spectrum technique for outflow identification and the baseline 300 kyr protostellar lifetime adopted for the lifetime ratios.","marker":"Duarte-Cabral et al. (2013)"},{"why":"Establishes the observational link between protostellar accretion and CO outflow activity that grounds the outflow-based definition of prestellar cores.","marker":"Bontemps et al. (1996)"},{"why":"Defines the ALMA-IMF survey, its targets, and data products from which the 14 protocluster cores and line cubes are drawn.","marker":"Motte et al. (2022)"},{"why":"Earlier pre- and protostellar classification of the W43 cores with the same CO transition, used to validate the automated method at 95% agreement.","marker":"Nony et al. (2023)"},{"why":"Protostellar evolution models used to estimate how long protostars spend as ionising UCHII regions, correcting the protostar counts per mass bin.","marker":"Hosokawa & Omukai (2009)"},{"why":"Supplies the low-mass prestellar lifetime of ~1.2 Myr and lifetime-to-free-fall ratio (~8) that frame the comparison for the high-mass values.","marker":"Könyves et al. (2015)"},{"why":"Presents the line-cube reductions and the DCN(3–2) velocity fitting procedure used to set each core's V_LSR.","marker":"Cunningham et al. (2023)"},{"why":"Provides the PPMAP dust temperatures that set the adopted mass ranges for prestellar and protostellar cores.","marker":"Dell'Ova et al. (2024)"},{"why":"Statistical lifetime estimates of massive YSOs that constrain the adopted 300 kyr protostellar lifetime.","marker":"Mottram et al. (2011)"}],"fun_headline_variants":["30 starless high-mass cores; lifetimes 10–30× free-fall time","ALMA-IMF finds 12 massive starless cores, lifetimes 50–240 kyr","12 hefty prestellar cores, >16 Msun, persist 10–30 free-fall times","High-mass starless cores: 30 found, collapse stalled 10–30×","Most massive starless cores yet: 12 in crowded protocluster centers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire prestellar classification rests on the premise, which the paper itself flags in Sect. 7.8, that any massive, actively accreting protostar inside a core would always produce a detectable CO or SiO outflow, so a core without one must be starless; if the accretion-ejection link can silently fail — through a dead zone, a nearly face-on outflow, or an interrupted accretion episode — some of the 12 candidates would actually hide young protostars and the derived lifetimes would shrink.","fun_headline_variants_meta":{"raw":{"variants":["30 starless high-mass cores; lifetimes 10–30× free-fall time","ALMA-IMF finds 12 massive starless cores, lifetimes 50–240 kyr","12 hefty prestellar cores, >16 Msun, persist 10–30 free-fall times","High-mass starless cores: 30 found, collapse stalled 10–30×","Most massive starless cores yet: 12 in crowded protocluster centers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000335,"raw_usage":{"total_tokens":1990,"prompt_tokens":1209,"completion_tokens":781,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":825,"completion_tokens_details":{"reasoning_tokens":668}},"tokens_in":825,"tokens_out":781,"duration_ms":7364,"temperature":1.0,"reasoning_tokens":668,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T21:32:57.143294+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A sub-arcsecond search for embedded heating in the 12 cores above 16 $M_\\odot$ — mid-infrared continuum or hot-core molecular emission at high spatial resolution — would settle the central claim: several detections of compact warm dust or hot-core chemistry would mean the outflow-free cores are not all starless, and the reported prestellar lifetimes would drop toward the free-fall times. The temporal version is equally concrete: re-observing the same fields after a few years and finding newly appeared outflow wings from any of the 12 candidates would falsify their prestellar status directly.","supporting_citations":[{"cited_title":"2013, , 558, A125","cited_arxiv_id":null,"evidence_quote":"Provides the On-Off differential-spectrum technique for outflow identification and the baseline 300 kyr protostellar lifetime adopted for the lifetime ratios."},{"cited_title":"2023, , 678, A194","cited_arxiv_id":null,"evidence_quote":"Presents the line-cube reductions and the DCN(3–2) velocity fitting procedure used to set each core's V_LSR."},{"cited_title":"2024, , 687, A217","cited_arxiv_id":null,"evidence_quote":"Provides the PPMAP dust temperatures that set the adopted mass ranges for prestellar and protostellar cores."},{"cited_title":"C., Hoare , M","cited_arxiv_id":null,"evidence_quote":"Statistical lifetime estimates of massive YSOs that constrain the adopted 300 kyr protostellar lifetime."}],"review_version":1}