{"id":"0d3b7e69-1716-4dfd-a8f7-6ec7ad446a22","arxiv_id":"2502.01261","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A 15-year Fermi-LAT analysis detects gamma rays from the HH 80-81 protostellar jet with a harder power-law spectrum and finds both leptonic and hadronic models can explain the emission.","lead":"Using 15 years of Fermi-LAT data, the authors detect gamma-ray emission from 300 MeV to 100 GeV near the HH 80-81 protostellar jet and argue the jet is the most likely source. The result supports the idea that massive protostellar jets can accelerate particles to relativistic energies, though the hadronic or leptonic origin remains undecided.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The association of the Fermi-LAT excess with HH 80-81 is not tested against the J1819 hypothesis: the localization is computed from a model that already places the source at HH 80-81, and the competing compact radio source is never fitted as a gamma-ray point source.","rationale":"The reader's conditional verdict is appropriate, and the reader's weakest assumption matches the concern identified here: the gamma-ray excess is attributed to HH 80-81 rather than to J1819 or to residual diffuse emission. The paper is careful to hedge the association as 'most probable' and explicitly acknowledges the ambiguity in the radiative models. My stress-test sharpens the concern rather than changing the verdict. The radiative fitting is honest about the leptonic/hadronic degeneracy, and the energy-injection timescales are presented as compatible, not definitive. The main soft spot is that the source association is decided by a localization that already assumes the source is at HH 80-81, and the only other plausible counterpart, J1819, is never fitted as a gamma-ray source. A simple likelihood comparison with a point source at J1819 would settle whether the data actually prefer HH 80-81. Until such a test is done, the central claim remains conditional on an assumption that is currently untested. Therefore the reader's CONDITIONAL verdict should stand unchanged.","tokens_in":19846,"tokens_out":5241,"duration_ms":56081,"concrete_test":"Re-run the binned likelihood analysis with the exact same ROI, energy range, and diffuse model, but place the point source at the J1819 coordinates (18h19m36.9s, -20d36m32s) instead of at IRAS 18162-2048, with no source at the IRAS position. Compare the global log-likelihood and TS of this model with the IRAS-centered model, and also fit both sources simultaneously with free normalizations. If the J1819-only model gives comparable or higher likelihood, or if the two-source fit does not clearly prefer the IRAS source, then the HH 80-81 association is not established. As a control for the diffuse background, repeat the comparison with the Galactic diffuse normalization scaled by +/-3%; if the TS difference between the two source placements is not robust to this variation, the association is even weaker.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the 15-year Fermi-LAT excess above 300 MeV is most plausibly powered by the HH 80-81 protostellar jet. The load-bearing step is the source identification in Sect. 3. The authors remove the catalog source 4FGL J1818.5-2036, insert a point source at the IRAS 18162-2048 coordinates, and then localize the excess with that model. They report that IRAS 18162-2048 falls inside the 1-sigma contour while J1819 falls outside the 2-sigma contour. This is not a discriminating test of the alternative hypothesis: with the source already fixed at IRAS coordinates, the localization is biased toward those coordinates, and the LAT point-spread function is several degrees at the energies where most of the 1445 predicted counts lie. J1819 is separated by only about 0.2 degrees from IRAS 18162-2048, so the two hypotheses cannot be separated by positional overlap alone. The non-variability argument is also weak: many flat-spectrum radio quasars and AGN do not show significant Fermi variability on 15-year timescales, and the light curve shown in Fig. 4 is dominated by large error bars and upper limits. In parallel, the paper states in Sect. 2.1 that a roughly 1% variation in the Galactic diffuse template corresponds to the measured source flux, yet no diffuse-model variation is propagated into the detection significance or the spectral index. The two most dangerous alternatives, J1819 as a gamma-ray source and a diffuse-template artifact, are therefore excluded by assumption rather than by a likelihood comparison.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes 15 years of Fermi-LAT data toward the HH 80-81 protostellar jet, reports a gamma-ray excess above 300 MeV, and uses positional, variability, and multiwavelength arguments to associate it with HH 80-81. It then fits the gamma-ray SED with leptonic (relativistic Bremsstrahlung) and hadronic (π0 decay) models, finding both compatible and deriving total particle energies and injection times that are consistent with the jet lifetime.","tokens_in":20225,"tokens_out":5764,"duration_ms":53114,"significance":"If the association with HH 80-81 holds, this would be an important piece of evidence that protostellar jets accelerate particles to relativistic energies and may contribute to Galactic cosmic rays. The paper uses the full 15-year Fermi-LAT dataset and openly documented analysis tools (fermipy, naima, GAMERA). The radiative modeling is clearly described, the authors acknowledge the leptonic/hadronic degeneracy, and they give quantitative energy and timescale arguments. However, the central source identification has not been properly confronted with the principal alternative, the compact radio source J1819, and the detection claim relies on a diffuse-model assumption that is not tested for systematics.","major_comments":[{"comment":"The exclusion of all data below 300 MeV is justified by the statement that a roughly 1% variation in the Galactic diffuse template corresponds to the measured source flux, yet no systematic uncertainty from the diffuse model is propagated into the detection significance or the spectral index. This is load-bearing because the claimed 5σ excess and the hard spectrum (Γ=2.62±0.12) are derived with a fixed gll_iem_v07 template. Please repeat the likelihood analysis with alternative diffuse-model normalizations or different interstellar emission templates and report how the TS and spectral index change. If the excess does not survive at >5σ under reasonable diffuse-model variations, the central detection claim is not robust.","section":"Sect. 2.1 (Table 1, Fig. 2)"},{"comment":"The association with HH 80-81 is not actually tested against the J1819 hypothesis. The localization is computed after replacing 4FGL J1818.5-2036 with a point source fixed at the IRAS coordinates, so the resulting best-fit position is biased toward those coordinates. To make the identification discriminating, the authors should fit a point source at the J1819 coordinates and compare the maximum likelihood (or TS) with the IRAS-centered model, and also show the localization result when the test source is initialized at J1819. Given that the LAT PSF is several degrees for the bulk of the 1445 predicted counts and the separation is only about 0.2 deg, the statement that J1819 is outside the 2σ region cannot by itself exclude this alternative. The non-variability argument is also weak: the light curve in Fig. 4 is dominated by non-detections and upper limits, and many flat-spectrum radio quasars do not show variability in Fermi data on 15-year timescales.","section":"Sect. 3 (Fig. 5, Table 3)"}],"minor_comments":[{"comment":"There are numerous typographical and rendering errors (e.g., 'di ffusive' in the abstract, 'spacial' in Sect. 2.1, 'loosing' for 'losing' in Sect. 4, and square symbols that replace minus signs, degree signs, and other LaTeX constructs in the provided text). These should be corrected before publication.","section":"Abstract and throughout"},{"comment":"The cross-reference to 'Table 3' appears before that table is introduced; Table 3 is in Sect. 3. Use a forward reference or restructure.","section":"Sect. 2.2"},{"comment":"The placeholder 'GGD catalog (?)' should be replaced with a proper catalog name or citation.","section":"Sect. 3"},{"comment":"The Araudo et al. (2008) reference is incomplete: 'arXiv preprint arXiv: . . . [arXiv:arXiv:0806.2306v1]' needs a full bibliographic entry.","section":"References"},{"comment":"The phrase 'the detected spectrum can be explained by both hadronic and leptonic particle components' is a fitted consequence rather than a test, as the paper itself acknowledges the degeneracy in Sect. 4. Consider rephrasing to indicate that the data are consistent with both models without implying a unique physical conclusion.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has two load-bearing weaknesses that the authors can address with additional analysis: a diffuse-template systematic test and a direct likelihood comparison of the J1819 hypothesis versus the IRAS-centered hypothesis. The paper's scientific value is moderate and the subject is topical for A&A. The character-substitution artifacts in the source file also need editorial cleanup."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper gives the deepest look yet at the Fermi-LAT excess near HH 80-81, with 15 years of data and a spectrum that extends above 700 MeV. The new spectral index of 2.62±0.12 is meaningfully harder than Yan et al.'s 3.53, and the authors are careful to check point-source morphology, spectral curvature, and variability. The radiative fits are honest: both hadronic and leptonic models can reproduce the SED, and they do not claim to break the degeneracy. The likelihood gain (~42) from replacing the 4FGL source with a point source at the IRAS position is a real piece of evidence favoring the protostellar jet over the catalog position.\n\nThe soft spots are in the source identification. The 300 MeV energy cutoff is justified by noting that a ~1% change in the Galactic diffuse template equals the source flux, but no systematic error from the template is propagated through the detection significance or the spectral index. The index could shift outside the quoted uncertainty. More seriously, the localization is performed with the source already fixed at the IRAS coordinates, so the statement that J1819 lies outside the 2σ contour is not a fair test: the model is biased toward IRAS. The authors never fit J1819 as a candidate gamma-ray point source, so the two hypotheses are not directly compared. The non-variability argument is also weak, since many AGN are steady in gamma-rays over 15 years and the light curve is mostly upper limits.\n\nThese are not fatal flaws. The paper is clearly written, uses standard tools, and the conclusion is appropriately hedged as 'most probable.' But the central claim is conditional on excluding the J1819 and diffuse-residual alternatives by assumption rather than by likelihood comparison. I would send this to a referee, but the referee should ask for a dedicated J1819 fit and a diffuse-model systematics study. If those come back clean, the case for a protostellar jet accelerator is considerably stronger.\n\nWorth discussing in a reading group as a case study in source association pitfalls.\n\nBest.","headline":"A useful, honest Fermi-LAT study of HH 80-81, but the source association rests on a circular localization and an unpropagated diffuse-systematic, so the central claim remains conditional.","tokens_in":20797,"tokens_out":4145,"would_cite":true,"duration_ms":37690,"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":"After 15 years of Fermi-LAT data, the paper argues that the gamma-ray excess above 300 MeV in this region most probably comes from the HH 80-81 protostellar jet, not from the neighboring radio source J1819, and that the spectrum can be…","keywords":["cosmic rays","gamma rays","massive young stellar objects","IRAS 18162-2048","protostellar jets","HH 80-81"],"falsifier":"Track the >300 MeV excess with an instrument of better angular resolution than Fermi-LAT, such as one of the next-generation imaging atmospheric Cherenkov arrays: if the centroid is resolved onto J1819, or if the flux varies on month-to-year timescales in step with J1819's radio flaring, the paper's association fails. A matched hard-X-ray observation of J1819 that detects emission at or above the level expected for an AGN would also remove the paper's strongest alternative candidate.","tokens_in":19668,"feed_emoji":"🌠","tokens_out":12072,"duration_ms":99685,"temperature":0.7,"pith_summary":"This paper sets out to determine whether the HH 80-81 protostellar jet, a 10-parsec outflow driven by a ~20-solar-mass protostar, can accelerate particles to relativistic energies and emit gamma rays. Working with 15 years of Fermi-LAT data from 300 MeV to 100 GeV, the authors find a ~5-sigma point-like excess at a position compatible within 1 sigma with IRAS 18162-2048, the star that powers the jet. They go through every plausible counterpart in the error region and conclude that HH 80-81, not the compact radio source J1819 nor residual Galactic diffuse emission, is the most probable source of the excess. Radiative fits show that both a hadronic population (proton-proton collisions producing pion decay) and a leptonic population (relativistic bremsstrahlung) can reproduce the measured spectrum, with total particle energies and injection times that fit inside the jet's ~40,000-year lifetime. If the association holds, the jet is a genuine gamma-ray emitter and a local accelerator of relativistic particles.","feed_headline":"5σ gamma-ray excess points to protostellar jet HH 80-81","feed_subtitle":"15 years of Fermi-LAT data place the excess on the jet; protons or electrons could power it within its lifetime.","key_machinery":"The physical mechanism assumed is diffusive shock acceleration, in which charged particles gain energy by crossing the shocks formed where the jet rams into dense interstellar material; the paper leans on the existing detection of polarized non-thermal radio emission from HH 80-81 as evidence that this acceleration is happening. The modeling machinery is a power-law-with-exponential-cutoff particle distribution, $\\phi(E) = \\phi_0 (E / 1\\ \\mathrm{GeV})^{-\\Gamma} \\exp(-E/E_{\\mathrm{cutoff}})$, with the gamma-ray photon index carried over to the particle index and with cutoff energies set by cooling: about $3$ TeV for electrons (synchrotron) and about $12$ TeV for protons (pion decay at $100\\ \\mathrm{cm}^{-3}$). Radiative fitting then compares relativistic bremsstrahlung and proton-proton pion decay against the Fermi-LAT spectrum, and the decisive diagnostic is the implied total energy and injection time relative to the jet's known lifetime.","core_discovery":"On its own terms, the paper's central claim is that the gamma-ray excess detected above 300 MeV is real and belongs to HH 80-81. After swapping the catalog source 4FGL J1818.5-2036 for a point-like source at the protostar's coordinates, the likelihood improves substantially, and the source has a test statistic $\\mathrm{TS}\\approx 29$ above 300 MeV (about $5\\sigma$) and a power-law spectrum with photon index $\\Gamma = 2.62 \\pm 0.12$, harder than the previous detection. The source is stable over 15 years, pointing away from a blazar interpretation, and the compact radio source J1819 sits outside the $2\\sigma$ containment while the protostar sits inside the $1\\sigma$ region. Both a leptonic model (relativistic bremsstrahlung) and a hadronic model ($\\pi^0$ decay from proton-proton collisions) fit the spectrum, with inverse Compton scattering ruled out energetically; the required particle energies are $(1.4 \\pm 0.3)\\times 10^{46}$ erg and $(2.9 \\pm 0.5)\\times 10^{47}$ erg for a density of $100\\ \\mathrm{cm}^{-3}$, and the corresponding injection times, about $10^3$ yr and $2\\times 10^4$ yr, are compatible with the jet lifetime of about $4\\times 10^4$ yr. The paper therefore concludes that HH 80-81 is the most probable counterpart, while leaving the hadronic-versus-leptonic question open.","pith_inferences":["If the excess is hadronic, the implied proton energy of roughly $3\\times10^{47}$ erg accumulated over about $2\\times10^4$ years makes HH 80-81 a non-negligible local source of cosmic rays; extrapolating to the population of massive protostellar jets would require knowing their duty cycle, which this paper does not address.","The modeling degeneracy could be broken by the predicted cutoffs: electrons should run out of energy near $3$ TeV and protons near $12$ TeV, so a future imaging atmospheric Cherenkov telescope detection above a few hundred GeV with a measured cutoff would discriminate the two populations.","The paper's spatial coincidence between the gamma-ray excess and the L291 molecular gas suggests a hadronic contribution, because pion decay scales with target density; mapping the excess with better angular resolution could test whether the emission follows the dense gas rather than the jet axis.","A direct extension would be to use the gamma-ray flux as a calorimeter: combining the measured total energy with an assumed cosmic-ray injection spectrum yields an estimate of the jet's mechanical feedback on its host cloud, a step the paper does not take."],"forward_implications":["The gamma-ray excess above 300 MeV is a genuine point-like detection at about 5 sigma, so HH 80-81 can be studied as a gamma-ray-emitting protostellar jet and not merely as a catalog coincidence.","Because both hadronic and leptonic models fit, the particle accelerator at the termination shock can supply either relativistic electrons or protons, and the injected energy is consistent with a few percent of the jet's kinetic power.","The 100-300 MeV band is too contaminated by Galactic diffuse emission for reliable spectral measurements; future analyses of this source should start at 300 MeV or higher.","Inverse Compton scattering is effectively excluded as the dominant gamma-ray mechanism under the assumed IR photon density, narrowing the radiative channels to bremsstrahlung and pion decay.","A future gamma-ray flare coincident with a radio or infrared flare from the jet would settle the association between young stellar objects and gamma-ray emission."],"supporting_citations":[{"why":"Supplies the only detection of linear polarization consistent with non-thermal synchrotron along a protostellar jet, the physical evidence that HH 80-81 accelerates particles.","marker":"Carrasco-González et al. (2010)"},{"why":"Provides the earlier 100 MeV-1 GeV gamma-ray detection that this work extends and revises with 15 years of data and a different low-energy treatment.","marker":"Yan et al. (2022)"},{"why":"Supplies the X-ray non-thermal component and the magnetic-field value used in the radiative fits.","marker":"Rodríguez-Kamenetzky et al. (2019)"},{"why":"Gives the radiative cooling-timescale equations used to compute cutoff energies and injection times for the jet.","marker":"Bosch-Ramon et al. (2010)"},{"why":"Provides the acceleration-efficiency estimate and the prediction that protons reach higher cutoff energies than electrons.","marker":"Araudo et al. (2021)"},{"why":"Supplies the jet lifetime of about $4\\times10^4$ yr against which the particle injection times are compared.","marker":"Qiu et al. (2019)"},{"why":"Supplies the adopted 1400 pc distance that sets the physical scale and luminosity of HH 80-81.","marker":"Zucker et al. (2020)"},{"why":"Supplies the ambient density of $100\\ \\mathrm{cm}^{-3}$ and outflow velocities used in the hadronic and leptonic models.","marker":"Bally & Reipurth (2023)"}],"fun_headline_variants":["5σ gamma-ray excess points to protostellar jet HH 80-81","Protostellar jet HH 80-81 likely gamma-ray source","Fermi-LAT links HH 80-81 jet to gamma-ray excess","HH 80-81 jet candidate for gamma-ray production"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the association of the gamma-ray excess with HH 80-81 itself, which rests on positional overlap within one sigma, the absence of variability, and the lack of X-ray detection of the alternate candidate J1819; if the excess is actually powered by J1819 or is an artifact of the Galactic diffuse template, the conclusion that the protostellar jet accelerates relativistic particles does not follow.","fun_headline_variants_meta":{"raw":{"variants":["5σ gamma-ray excess points to protostellar jet HH 80-81","Protostellar jet HH 80-81 likely gamma-ray source","Fermi-LAT links HH 80-81 jet to gamma-ray excess","HH 80-81 jet candidate for gamma-ray production"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001469,"raw_usage":{"total_tokens":6062,"prompt_tokens":1254,"completion_tokens":4808,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":870,"completion_tokens_details":{"reasoning_tokens":4729}},"tokens_in":870,"tokens_out":4808,"duration_ms":37227,"temperature":1.0,"reasoning_tokens":4729,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T15:52:01.574672+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track the >300 MeV excess with an instrument of better angular resolution than Fermi-LAT, such as one of the next-generation imaging atmospheric Cherenkov arrays: if the centroid is resolved onto J1819, or if the flux varies on month-to-year timescales in step with J1819's radio flaring, the paper's association fails. A matched hard-X-ray observation of J1819 that detects emission at or above the level expected for an AGN would also remove the paper's strongest alternative candidate.","supporting_citations":[{"cited_title":"Detection of Gamma-Rays from the Protostellar Jet in the HH 80-81 System","cited_arxiv_id":"1908.10994","evidence_quote":"Provides the earlier 100 MeV-1 GeV gamma-ray detection that this work extends and revises with 15 years of data and a different low-energy treatment."},{"cited_title":"E., Araudo, A","cited_arxiv_id":null,"evidence_quote":"Gives the radiative cooling-timescale equations used to compute cutoff energies and injection times for the jet."},{"cited_title":"T., Padovani, M., & Marcowith, A","cited_arxiv_id":null,"evidence_quote":"Provides the acceleration-efficiency estimate and the prediction that protons reach higher cutoff energies than electrons."},{"cited_title":"2019, ApJ, 871, 141","cited_arxiv_id":null,"evidence_quote":"Supplies the jet lifetime of about $4\\times10^4$ yr against which the particle injection times are compared."},{"cited_title":"S., Schlafly, E","cited_arxiv_id":null,"evidence_quote":"Supplies the adopted 1400 pc distance that sets the physical scale and luminosity of HH 80-81."},{"cited_title":"& Reipurth, B","cited_arxiv_id":null,"evidence_quote":"Supplies the ambient density of $100\\ \\mathrm{cm}^{-3}$ and outflow velocities used in the hadronic and leptonic models."}],"review_version":1}