{"id":"ba3f9b3b-b234-4ea7-8df7-b6216921e203","arxiv_id":"2501.12767","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Gamma-ray emission up to about 1 TeV is detected from four HII regions hosting young massive star clusters, supporting stellar winds as a contributor to Galactic cosmic rays.","lead":"Astronomers analyzed Fermi satellite data toward four star-forming clouds and found gamma-ray glow reaching energies up to about one trillion electron volts. The result suggests winds from very young, massive stars can accelerate cosmic rays, and future ground-based telescopes should be able to test it.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1 TeV claim rests on residual TS maps after removing 4FGL sources; without a likelihood test adding an extended WISE template, those residuals could be diffuse-model or source-confusion artifacts.","rationale":"The paper's logic has two stages. Section 2 establishes a statistical correlation between unidentified 4FGL sources and WISE HII regions; that result is quantitative (127 vs ~50, sigma=12) and is a useful independent hint. However, a catalog correlation does not by itself identify which HII regions are gamma-ray emitters; it only supplies candidates. The physical claim in the abstract and conclusions comes from Section 3, where four selected regions are examined in Fermi-LAT data. The only evidence shown is residual TS maps obtained after the 4FGL sources near the HII regions are removed from the background model. The interpretation of those maps as genuine extended emission assumes that the 4FGL catalog plus the standard diffuse templates constitute a complete and accurate model of the sky in each field. In the inner Galactic plane this assumption is known to be fragile: the interstellar emission model is uncertain at the tens-of-percent level, and source confusion means that unresolved or mis-associated sources are common. If the true sky in one of these fields contains a component not in the model—a molecular cloud illuminated by cosmic rays, an unresolved pulsar, or a variation of the diffuse background—the residual map will reproduce that component and can mimic spatial correlation with the WISE dust emission. The paper does not perform the decisive control: a likelihood-ratio test that adds an extended source with the WISE morphology to the full model and quantifies the improvement. It also does not give uncertainties for the highest-energy flux points, so the 'up to 1 TeV' statement cannot be checked. This is the single most load-bearing concern because it attacks the primary evidence itself: if the residual maps are artifacts, the correlation remains interesting but the four-region detection and the pre-supernova wind-acceleration conclusion collapse. The other issues raised by the reader—lack of code/data, unspecified selection function—concern reproducibility rather than correctness and would not, by themselves, invalidate the claim if the residual analysis were sound. I therefore agree with the reader's weakest assumption and do not propose changing the CONDITIONAL verdict.","tokens_in":4118,"tokens_out":10034,"duration_ms":112020,"concrete_test":"Using Fermitools with the same data selection and diffuse models as Peron et al. (2024), re-fit each of the four fields twice: (a) keep all 4FGL sources in the model, (b) add an extended source whose spatial template is the WISE 22-um image of the HII region. Report TS_ext = 2*Delta(logL) for (b) vs (a), the best-fit extension (if left free), and the TS of the highest-energy bin (E>100 GeV) with a trials factor for the four regions and energy bins. If TS_ext < 25 in any field, or the best-fit extension is inconsistent with the WISE radius, or the highest-energy bin is consistent with an upper limit at 95% confidence, the residual maps of Fig. 2 do not establish the claimed association or the 1 TeV emission.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 constructs the detection by optimizing a background model of 4FGL point sources plus diffuse templates, then removing the 4FGL sources listed in Table 1 and presenting the leftover TS maps in Fig. 2. The central claim ('emission up to 1 TeV' from the four WISE regions) is load-bearing on the null hypothesis that the remaining sky is fully described by that model. All four regions are at |b| <= 3.5 deg and l ~ 18-52 deg, i.e. the crowded inner Galaxy, exactly where the interstellar emission model has known systematics and where 4FGL is incomplete or confused. Any mismodeling of the diffuse background or of a non-removed point source will appear in the residual map and can peak at the HII position, especially after the 'potentially associated' sources are deleted. The manuscript reports neither the TS for adding an extended source with the WISE template to the full model nor the per-bin significance of the highest-energy points; the SED was extracted with the extension fixed to WISE, so Fig. 3 cannot by itself validate the morphology. The discussion even hedges: 'further study would be needed to confirm whether this type of sources could be seen at higher energies.' Thus the 1 TeV part of the claim is not independently testable from this text.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a Fermi-LAT analysis of four WISE HII regions associated with young massive star clusters, claiming gamma-ray emission up to 1 TeV from G018.426+01.922, G028.746+03.458, G040.554+02.443, and G051.978+00.542. The detection is based on residual TS maps obtained after removing nearby 4FGL sources from the background model, and the spectral energy distributions are extracted with the source extension fixed to the WISE size. The paper also summarizes a statistical correlation between unidentified Fermi-LAT sources and WISE HII regions, reporting a global significance of Σ=12 and Σ>35 in the inner Galaxy. The authors conclude that these young systems are pre-supernova cosmic-ray accelerators and that future CTAO/ASTRI observations can test the emission at TeV energies.","tokens_in":4317,"tokens_out":4060,"duration_ms":39316,"significance":"If the claimed 1 TeV emission is confirmed, it would provide strong evidence that collective stellar winds in very young massive star clusters accelerate particles before any supernova has occurred, with direct implications for the contribution of star clusters to Galactic cosmic rays. The statistical correlation between HII regions and unidentified Fermi sources is a useful quantitative step toward resolving source confusion. The paper also benefits from explicitly stating the young ages (< 3 Myr) of the targets, which cleanly separates the stellar-wind component from supernova remnants, and the Monte Carlo approach to evaluate the correlation is appropriate. However, the current manuscript does not yet present the statistical evidence needed to fully support the central 1 TeV claim, as detailed in the major comments.","major_comments":[{"comment":"The claim of emission up to 1 TeV rests entirely on residual TS maps after subtracting 4FGL sources and a diffuse model. The manuscript never reports the test statistic for adding an extended source with the WISE template to the full background model, nor the per-energy-bin significance of the highest-energy spectral points. In the crowded inner Galaxy (|b| ≤ 3.5°, l ~ 18–52°), the interstellar emission model is known to have systematics and 4FGL is incomplete or confused, so the residual TS peaks could be artifacts of mismodeled diffuse emission or unremoved point sources. The paper itself hedges in Section 4 ('Further study would be needed to confirm whether this type of sources could be seen at higher energies'), which underscores that the present evidence is not conclusive. A likelihood-ratio test comparing the full model with and without the WISE-extended source, together with energy-resolved TS maps, is required to support the detection claim.","section":"Section 3, Fig. 2"},{"comment":"The spectral energy distribution is extracted with the source extension fixed to the WISE size, so the SED cannot independently validate the morphological association. The claimed morphological match between gamma-ray residuals and infrared emission is only qualitative (visual comparison in Fig. 2). The paper should provide a quantitative morphology test, e.g., TS_ext comparing a point-like source against a source extended with the WISE shape, and report the best-fit position and extension with uncertainties. Without this, the statement that 'the morphology of the Fermi-LAT emission matches well with the infrared emission' is not tested against alternative hypotheses.","section":"Section 3, Fig. 3"},{"comment":"The Monte Carlo significance calculation for the correlation between unidentified Fermi sources and WISE HII regions is described too briefly to assess the validity of the reported Σ values, especially for the longitude-sliced analysis. It is unclear whether the 1000 simulated catalogs with 'random extraction in l and b' preserve the actual source density or the longitude-dependent distribution when slices of 60° are considered. If the simulation draws uniformly over the whole sky, the high Σ > 35 in the inner Galaxy could reflect the higher real source density there rather than a physical association. The manuscript should specify the simulation procedure for the slice test, or explicitly refer to the detailed derivation in Peron et al. (2024) with the necessary information.","section":"Section 2, Fig. 1"}],"minor_comments":[{"comment":"The caption describes 'light-blue circles' for the HII region extensions, but the text in Section 3 refers to 'red circles' for the 4FGL source localization uncertainties; the color coding in the figure and text should be made consistent.","section":"Section 3, Fig. 2 caption"},{"comment":"The column header contains a typo: '4FGL souces' should be '4FGL sources'. Also, the table does not give uncertainties on the extension and 22 μm flux, which would be helpful for judging the precision of the association.","section":"Table 1"},{"comment":"The abstract uses 'we will present' for the statistical investigation, but the results are already presented in Section 2; the tense should be adjusted to reflect the actual content.","section":"Abstract"},{"comment":"There are two 'Peron et al. 2024' entries with different coauthor lists (Nature Astronomy and ApJ 972, L22); the in-text citations do not distinguish them. The authors should label them 2024a and 2024b or use the author abbreviations consistently.","section":"References"},{"comment":"The Fermi-LAT analysis details (version of the fermitools, event selection, energy range, and the specific diffuse models used) are not stated. For a proceedings paper this may be acceptable, but a reference to the detailed analysis in Peron et al. (2024) would improve reproducibility.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the residual TS maps is legitimate and directly affects the central claim of the paper. The manuscript currently does not provide a likelihood test that would distinguish genuine extended emission from diffuse-model or source-confusion artifacts. I agree with the conditional verdict: the paper needs a major revision with additional statistical tests. The authors should also clarify the simulation procedure for the longitude-sliced correlation, as the current description is too ambiguous to evaluate the reported Σ values."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Giada,\n\nQuick read: this is a conference proceedings, not a full paper. Most of the content is a summary of Peron et al. 2024 (Nature Astronomy and ApJ 972, L22), which the authors cite openly. The genuinely new bit is the follow-up analysis of four WISE HII regions with Fermi-LAT, presented as residual TS maps and SEDs up to 1 TeV. The statistical correlation between WISE HII regions and unidentified 4FGL sources, with Sigma ~ 12 and >35 in the inner Galaxy, is a strong result and is well explained here.\n\nWhat the paper does well: the selection criterion (young clusters, <3 Myr, so no SN contamination) is physically motivated, and the correlation method using Monte Carlo against random catalogs is sound. The residual maps do show emission that spatially matches the infrared tracers, which is suggestive. For a proceedings, the description is clear and the references point to the detailed analyses.\n\nSoft spots, in proportion: the 1 TeV claim is load-bearing but rests on residual TS maps after removing 4FGL sources and the diffuse model. The text does not report a likelihood test that adds an extended source with the WISE template, nor per-bin significances for the highest-energy points. All four regions sit at |b| <= 3.5 deg in the inner Galaxy, where diffuse mismodeling and source confusion are known issues. So the emission could in principle be an artifact of the leftover background. That said, this is a proceedings—the full details are in the referenced papers, and the strain does not necessarily mean the result is wrong. It means this text alone cannot independently confirm the 1 TeV part of the claim. The selection function for the follow-up targets is also not quantified, but again, the proceedings format is compressed.\n\nWho this is for: someone wanting a quick overview of the star-cluster/HII-region gamma-ray connection, or a citable pointer to the team's broader work. The paper itself does not deserve a serious referee as a standalone research article—its value is as a summary—but the underlying science does merit attention. If the authors submit this as a full paper, it needs the missing likelihood analysis and a fuller background treatment.\n\nMy call: if this crossed my desk as a proceedings contribution, I'd let it through without fuss. If it were submitted as a regular paper, I'd send it to review but with the expectation that the 1 TeV claim be backed by an extended-template test. For your own work, cite the ApJ paper, not this proceedings.\n\nBest,\n[Your name]","headline":"A compact proceedings paper that mostly repackages already-published results; the four HII-region follow-ups and the strong correlation are real but need the cited ApJ papers for full support.","tokens_in":4892,"tokens_out":1413,"would_cite":false,"duration_ms":16916,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Young massive star clusters embedded in HII regions can be detected in gamma rays before any supernova explodes, and Fermi-LAT data reveal four such regions emitting up to about 1 TeV.","keywords":["star clusters","HII regions","gamma-ray astronomy","Fermi-LAT","cosmic-ray acceleration","stellar winds","unidentified gamma-ray sources","WISE catalog"],"falsifier":"Recomputing the residual maps with an updated Fermi source catalog or a revised interstellar emission model and finding that the four excesses disappear, or assigning any of them to a known pulsar or supernova remnant, would overturn the claimed association.","tokens_in":3884,"feed_emoji":"🌌","tokens_out":7163,"duration_ms":62623,"temperature":0.7,"pith_summary":"The paper sets out to show that very young massive star clusters, still wrapped in the gas clouds where they formed, are gamma-ray emitters in their own right. Because the clusters are younger than about 3 million years, no supernova has yet exploded in them, so any gamma-ray emission must come from particles accelerated by the collective winds of the massive stars themselves. Analyzing Fermi-LAT data, the authors find a highly significant positional correlation between unidentified gamma-ray sources and infrared-selected HII regions, with 127 matches versus about 50 expected by chance, and targeted analysis of four regions reveals emission up to roughly 1 TeV whose shape follows the infrared emission. If this reading is right, stellar winds alone can accelerate cosmic rays in dense environments, and young clusters contribute to the Galactic cosmic-ray budget before they ever produce supernova remnants.","feed_headline":"Four HII regions emit gamma rays up to 1 TeV before any supernova","feed_subtitle":"Fermi-LAT residual emission matches infrared images, tying unidentified sources to stellar-wind acceleration.","key_machinery":"The analysis rests on two linked tools. The first is a catalog cross-correlation: unidentified sources from the Fermi-LAT 4FGL catalog are matched to the WISE catalog of HII regions, an infrared 22-micron sample that traces hot dust around massive young stars, and the significance of the number of matches is assessed against 1000 randomized catalogs through the statistic Sigma = (N_real - <N_sim>)/sigma_sim. The second is a Fermi-LAT residual analysis: the background model includes Galactic and extragalactic diffuse emission and 4FGL point sources; after optimization, the sources overlapping each candidate HII region are removed, and the leftover test-statistic maps are compared with WISE infrared contours. The physical mechanism carrying the argument is hadronic gamma-ray production, in which protons accelerated by stellar winds collide with the dense gas of the HII region, at densities of $10^{2}$ to $10^{3}$ $cm^{-3}$, and produce gamma rays.","core_discovery":"On the paper's own terms, the discovery is that the four WISE HII regions G018.426+01.922, G028.746+03.458, G040.554+02.443, and G051.978+00.542 are genuine GeV gamma-ray sources: after the standard Fermi-LAT background, made of Galactic and extragalactic diffuse emission plus 4FGL point sources, is optimized and the catalog sources overlapping each region are removed, residual test-statistic maps show emission up to about 1 TeV that is spatially matched to the 22-micron infrared morphology of the HII regions. Because these clusters are younger than the supernova timescale, the authors attribute the emission to hadronic interactions of protons accelerated by stellar winds in gas with density of order $10^{2}$ to $10^{3}$ $cm^{-3}$. The accompanying catalog-wide correlation, with significance Sigma = 12 overall and greater than 35 in the inner Galaxy, is used to argue that the four targeted objects are not isolated curiosities but representatives of a population of pre-supernova cosmic-ray accelerators.","pith_inferences":["If the correlation generalizes, a sizable share of the currently unidentified Fermi-LAT sources in the inner Galaxy could be young stellar clusters, which would shift the inferred source population of the 4FGL catalog.","A natural test is to rank all matched HII regions by 22-micron flux and gas density and predict which should be the next brightest gamma-ray sources; the present paper selects by infrared brightness and size, so this ranking is not yet tested.","The same age-based selection could be applied to older clusters above about 3 million years to isolate when supernova remnants start to dominate; if emission persists or brightens past that age, the pre-supernova attribution would need revision.","Because the HII-region gas provides a target for hadronic collisions, neutrino production at comparable flux is expected from the same proton population, so a neutrino observation from one of the four directions would independently test the hadronic origin."],"forward_implications":["If the four detections hold, young massive clusters are confirmed as pre-supernova particle accelerators, and the observed emission traces wind-accelerated protons interacting with dense gas.","The statistical correlation, with Sigma about 12 overall and above 35 in the inner Galaxy, implies that many unidentified Fermi sources coincident with HII regions are likely the same kind of stellar-wind accelerators, not chance alignments.","The measured GeV fluxes fall within the sensitivity of the planned ASTRI and CTAO arrays, so the next generation of Cherenkov telescopes can test whether the emission extends to TeV energies and probe the parent proton spectrum.","With acceleration efficiency around one percent, as found in the earlier Vela pilot study, the total power from such clusters may be large enough to contribute a meaningful fraction of the Galactic cosmic rays, supplementing supernova remnants at the high-energy end."],"supporting_citations":[{"why":"Supplies the WISE 22-micron HII region catalog used for the correlation analysis and for selecting the four targets.","marker":"Anderson et al. 2014"},{"why":"Supplies the 4FGL catalog of Fermi-LAT sources, used both for the statistical matching and for the background model from which candidate sources are removed.","marker":"Abdollahi et al. 2022"},{"why":"The earlier Vela molecular cloud pilot study that first detected HII-region gamma-ray emission and established the order-one-percent acceleration efficiency that motivates this extension.","marker":"Peron et al. 2024"},{"why":"Provides the extended-source sensitivity curves for ASTRI and CTAO used to compare the measured spectra with future observability.","marker":"Celli & Peron 2024"},{"why":"Gives the theoretical model of particle acceleration by stellar winds in young clusters that predicts gamma-ray emission.","marker":"Morlino et al. 2021"},{"why":"Characterizes the powerful winds of young massive clusters, supporting the acceleration environment assumed by the paper.","marker":"Celli et al. 2024"}],"fun_headline_variants":["HII regions shine in gamma rays before supernovae","Gamma rays from star clusters without supernovae","Fermi-LAT spots 4 HII regions as gamma-ray sources","Pre-supernova cosmic-ray accelerators found in HII regions","Young star clusters emit gamma rays up to 1 TeV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the 4FGL catalog plus the standard Galactic and extragalactic diffuse models completely describe everything in each field except the HII region, so the leftover gamma-ray signal is real emission rather than an artifact of an incomplete background model.","fun_headline_variants_meta":{"raw":{"variants":["HII regions shine in gamma rays before supernovae","Gamma rays from star clusters without supernovae","Fermi-LAT spots 4 HII regions as gamma-ray sources","Pre-supernova cosmic-ray accelerators found in HII regions","Young star clusters emit gamma rays up to 1 TeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000687,"raw_usage":{"total_tokens":3143,"prompt_tokens":1005,"completion_tokens":2138,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":621,"completion_tokens_details":{"reasoning_tokens":2052}},"tokens_in":621,"tokens_out":2138,"duration_ms":13568,"temperature":1.0,"reasoning_tokens":2052,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:48:26.251142+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recomputing the residual maps with an updated Fermi source catalog or a revised interstellar emission model and finding that the four excesses disappear, or assigning any of them to a known pulsar or supernova remnant, would overturn the claimed association.","supporting_citations":[{"cited_title":"D., Bania, T","cited_arxiv_id":null,"evidence_quote":"Supplies the WISE 22-micron HII region catalog used for the correlation analysis and for selecting the four targets."},{"cited_title":"2022, The Astrophysical Journal Supplement Series, 260","cited_arxiv_id":null,"evidence_quote":"Supplies the 4FGL catalog of Fermi-LAT sources, used both for the statistical matching and for the background model from which candidate sources are removed."},{"cited_title":"2021, Monthly Notices of the Royal Astronomical Society, 504, 6096","cited_arxiv_id":null,"evidence_quote":"Gives the theoretical model of particle acceleration by stellar winds in young clusters that predicts gamma-ray emission."}],"review_version":1}