{"id":"750cc3d2-597b-4dea-9b92-e0c8af35306f","arxiv_id":"1908.10994","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A Fermi-LAT analysis of ten years of data finds a point-like gamma-ray excess at the position of the HH 80-81 protostellar jet (TS about 101, photon index about 3.5) and attributes it most likely to the jet.","lead":"Using ten years of Fermi-LAT data, the authors report a gamma-ray excess toward the protostellar jet system HH 80-81, with a soft spectrum that extends only to about 1 GeV. If the association holds, this would be one of the first detections of gamma-ray emission from a protostellar jet, marking protostars as particle accelerators.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Detection rests on sub-200 MeV photons where the PSF is broad and Galactic diffuse emission dominates; normalization-only systematics do not test spatial model errors.","rationale":"The reader identified the diffuse-background assumption as the weakest point; the stress-test pass reaches the same conclusion and sharpens it. The paper itself provides the strongest evidence for this concern: the TS values for E>200 MeV and E>300 MeV are only 25 and 5, respectively, so the detection is entirely carried by low-energy photons where Fermi-LAT's PSF is large and the Galactic diffuse foreground is strongest. The ±3% normalization rescaling in §2.1 is a useful check but does not address spatial mis-modeling, which is the known dominant systematic for low-latitude, sub-GeV analyses. The source is point-like within the uncertainties, so an unresolved clump of interstellar emission or a faint uncatalogued source inside the 0.28° error circle cannot be excluded by the catalog cross-checks in §3. I also note a minor arithmetic slip in §4: L_gamma/L_jet is about 1.7%, not 0.017%, but this does not affect the central detection claim. The conditional verdict is appropriate: the claim is plausible but depends on a testable systematic that the current paper does not perform. If the proposed independent-template and energy-cut tests preserve the excess, the detection would be much more secure.","tokens_in":7591,"tokens_out":4558,"duration_ms":49265,"concrete_test":"Re-run the §2.1 binned likelihood for the same 10-year dataset using an independent interstellar emission template (e.g., the updated gll_iem_v08 or a GALPROP-based model with free normalization and index) and also with front-converted (PSF0/1) events only, which have better low-energy PSF. Require the TS at the best-fit position to remain >25 above 100 MeV and >10 above 200 MeV, and the best-fit position to remain within 0.3° of IRAS 18162-2048. If either condition fails, the excess is not robust against the dominant low-energy systematic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In §2.1 the authors report TS=101 above 100 MeV, but this drops to 25 above 200 MeV and 5 above 300 MeV. The entire detection significance is therefore supplied by events below ~200 MeV, where the Pass 8 SOURCE-class PSF is broad (>1°) and the low-latitude interstellar emission is intense and spatially complex. The only diffuse systematic tested in §2.1 is a uniform ±3% rescaling of the gll_iem_v07 normalization (TS 76–144). This does not cover the dominant uncertainty: the spatial morphology of the Galactic diffuse emission, including HI/CO gas templates, cosmic-ray gradients, and unresolved sources. Because the source is not resolved (Table 1) and the fitted error circle is 0.28°, a compact residual of mis-modeled interstellar emission or an unrecognized 4FGL-sub-threshold source within the circle could produce the same point-like excess. The cross-checks in §3 rule out catalogued emitters but cannot rule out an uncatalogued or diffuse origin. Absent a test with an independent interstellar emission model or an energy cut that preserves significance, the claim that the gamma-ray excess is 'likely produced in the HH 80-81 jet' is not yet established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports a Fermi-LAT search for gamma-ray emission from the protostellar jet HH 80-81, using ten years of Pass 8 SOURCE-class data centered on IRAS 18162-2048. The analysis finds a point-like excess with TS ≈ 101 above 100 MeV, a power-law photon index Γγ = 3.53 ± 0.11, and a flux Fγ = (5.2 ± 0.4) × 10^-8 photons cm^-2 s^-1 above 100 MeV. The excess is reported as non-extended and non-variable. Catalog cross-checks find no known gamma-ray counterpart within the 0.28° error circle. The authors argue that the jet has sufficient energy and suitable particle populations for the emission to be produced by pp interactions or inverse-Compton scattering. However, the reported significance drops to TS = 25 above 200 MeV and TS = 5 above 300 MeV, so the detection rests almost entirely on sub-200 MeV photons, where the Fermi point-spread function is broad and the Galactic diffuse emission is strong.","tokens_in":7812,"tokens_out":6731,"duration_ms":64675,"significance":"If the detection were robust, this would be an important result: it would be the first gamma-ray detection from a protostellar jet and would support particle acceleration in young stellar outflows. The paper uses standard Fermi-LAT likelihood tools and includes several useful checks: TS maps, an extension test, variability analysis, and catalog cross-matching. The energy-budget argument is simple and transparent. The main weakness is that the central detection significance is not robust to energy cuts and the diffuse-background systematic is only tested via a uniform ±3% normalization rescaling, not via spatial-template variations. The theoretical section also contains an explicit inconsistency between the bremsstrahlung scenario and the radio synchrotron spectrum. These issues are addressable with additional analysis, but they currently prevent the paper from substantiating its title claim as stated.","major_comments":[{"comment":"The detection significance is not robust to energy threshold: TS = 101 above 100 MeV, but only TS = 25 above 200 MeV and TS = 5 above 300 MeV. The claimed detection is therefore entirely driven by photons below ~200 MeV, where the Pass 8 SOURCE-class point-spread function is broad (>1°) and the low-latitude Galactic diffuse emission is bright and spatially complex. The only diffuse-background systematic tested, a uniform ±3% rescaling of the gll_iem_v07 normalization, does not probe the dominant uncertainty: the spatial morphology of the interstellar emission (e.g., gas templates, cosmic-ray gradients, unresolved sources). A compact residual of mis-modeled diffuse emission or an unrecognized 4FGL-sub-threshold source within the 0.28° error circle could mimic or displace this point-like excess. The authors should repeat the analysis with an alternative interstellar emission model or a spatial-template variation, or demonstrate that the excess remains significant above 200 MeV after a more complete treatment of low-energy systematics.","section":"§2.1 and Figure 1"},{"comment":"The disk-extension test gives -2ΔL values of 0.24, 0.86, 0.16, and -1.74 for disk radii 0.1°, 0.3°, 0.5°, and 0.7°, respectively. These values are all well below the TSext > 16 threshold used in Fermi-LAT work, so the test cannot distinguish a point source from a compact extended or diffuse source. Given that the excess is dominated by low-energy photons with a large point-spread function, the conclusion that the emission is point-like is much weaker than the wording in §5 suggests.","section":"§2.2 and Table 1"},{"comment":"The cross-checks against SIMBAD, the Roma blazar catalog, and the ATNF pulsar catalog successfully rule out known cataloged sources, but they cannot exclude an uncataloged source or a compact residual of interstellar emission inside the 0.28° error circle. The nearest 4FGL source is 0.8° away, leaving ample room for a faint unmodeled source that would not appear in 4FGL. The paper should quantify the chance-coincidence probability for such an unresolved source or perform an additional test, such as adding trial sources at various positions within the error circle and examining the resulting TS, to support the identification with the HH 80-81 jet.","section":"§3"},{"comment":"The theoretical discussion does not currently demonstrate that the HH 80-81 jet can produce the observed gamma-rays under a self-consistent particle population. For the relativistic Bremsstrahlung scenario, the authors state that the electron spectral index se ≈ Γγ ≈ 3.5 predicts a synchrotron photon index (se + 1)/2 ≈ 2.3, which they explicitly note is inconsistent with the observed radio photon index of 1.3. For the inverse-Compton scenario, se = 2Γγ − 1 ≈ 6 predicts synchrotron emission peaking near ~0.2 eV and requires an ad hoc cutoff at γc < 3 × 10^5. The paper should either present a multi-wavelength model that simultaneously reproduces the radio and gamma-ray data, or clearly label the theoretical section as a qualitative plausibility argument rather than a quantitative demonstration.","section":"§4"}],"minor_comments":[{"comment":"The text contains a typo: 'pabel b' should read 'panel b'.","section":"§2.1"},{"comment":"The text contains a typo: 'Tabel 1' should read 'Table 1'.","section":"§2.2"},{"comment":"The phrase 'photon index of the observed radio spectrum is 1.3' is ambiguous: it should be stated explicitly whether the quoted value is a photon index Γ (S_ν ∝ ν^{-Γ+1}) or a spectral index α (S_ν ∝ ν^{-α}), because this directly affects the comparison with (se + 1)/2.","section":"§4"},{"comment":"The abstract and conclusions state that the jet properties suffice for producing the observed gamma-rays, but §4 itself notes an inconsistency for the Bremsstrahlung scenario. The summary should be softened to reflect that only a qualitative energy-budget argument is made, not a fully consistent radiative model.","section":"Abstract and §5"},{"comment":"The SED appears to show bins up to 10 GeV, while the text says the spectrum extends only to 1 GeV. Please clarify whether the high-energy points are upper limits, and if so, state the confidence level.","section":"Figure 2 and §2.3"},{"comment":"For the 1-year binned light curve over a 10-year observation, the number of bins and the degrees of freedom should be stated explicitly; the quoted χ²/d.o.f. = 14.02/8 implies roughly 9 bins, which should be consistent with the actual binning.","section":"§2.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a plausible candidate for publication in a letter-style venue if the authors can strengthen the detection claim. The current version's significance rests on sub-200 MeV photons, and the diffuse-systematic treatment is too narrow to support the strong title claim. The theoretical section also contains an admitted inconsistency with the radio data. These issues are within scope to fix: an alternative diffuse model or energy-threshold analysis, a chance-coincidence estimate, and a reframing of the theory section would materially improve the paper. I would not recommend rejection, as the analysis is standard and the scientific question is timely."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a plausible but not yet convincing detection of GeV emission from a protostellar jet. If confirmed, it's a new source class, but the significance lives entirely below 200 MeV, where Fermi's PSF is wide and the Galactic diffuse foreground is hardest to model.\n\nWhat's new: the excess is not in 4FGL, and the prior catalog candidate 3FGL J1819.5-2045c was removed. So the paper is genuinely reporting a new candidate. The analysis is standard and clean: 10-year Pass 8 data, TS map, spectral fit, extension test, variability search, and careful catalog cross-checks. The energy budget argument—jet kinetic luminosity of 3e35 erg/s versus the required 5e33 erg/s in gamma-rays—is fine, and the non-variability is consistent with a steady jet.\n\nThe weak part is the energy and background dependence. TS drops from 101 above 100 MeV to 25 above 200 MeV and 5 above 300 MeV. That's a big red flag. In the sub-200 MeV band, the PSF is degree-scale and the low-latitude interstellar emission is complex. The only diffuse systematic tested is a uniform ±3% normalization rescale, which just changes TS between 76 and 144. That does not address the more important spatial uncertainty—whether the emission model for the Galactic diffuse background is correct in this line of sight. The 0.28° error circle also leaves room for an unresolved or uncatalogued source. The cross-check against known catalogs is thorough but cannot exclude that possibility.\n\nThe theoretical section is plausible but not discriminating. For bremsstrahlung, the predicted synchrotron index doesn't match the radio data; for IC, the predicted synchrotron peak is at 0.2 eV where there is no measurement. pp is fine but generic. So the paper argues energy budget and cooling times are adequate, but does not provide a unique radiation mechanism.\n\nOverall, this is a paper I would send to peer review rather than desk-reject, because the source class would matter and the analysis is competently done. But a referee should insist on a more robust diffuse background treatment—independent interstellar emission models, or at least a spatial template variation—and an energy-stacked significance above 200 MeV. As it stands, the 'detection' is a candidate, not a claim that is nailed down.","headline":"Plausible new gamma-ray excess toward HH 80-81, but the detection is driven by sub-200 MeV photons and the diffuse-background systematics are too weak to nail the association.","tokens_in":8357,"tokens_out":3341,"would_cite":false,"duration_ms":33574,"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":"The HH 80-81 protostellar jet emits gamma rays: a ten-year Fermi-LAT analysis finds a point-like excess with TS > 100 toward the massive protostar IRAS 18162-2048, and the jet's kinetic power can supply the observed luminosity.","keywords":["gamma-ray astronomy","protostellar jets","HH 80-81","Fermi-LAT","non-thermal radiation","particle acceleration","ISM: jets and outflows","massive protostars"],"falsifier":"Re-run the analysis above 200 MeV with an independent diffuse background template built from gas tracers rather than the adopted standard model. The paper itself reports the excess falls to TS $\\approx 25$ above 200 MeV; if that value stays low or vanishes under the alternative background, the claimed detection would be unsupported, whereas a point-like excess that grows with more Fermi data and keeps TS above 100 would confirm the jet origin.","tokens_in":7398,"feed_emoji":"🌠","tokens_out":12642,"duration_ms":123049,"temperature":0.7,"pith_summary":"This paper reports the detection of gamma-ray emission from the HH 80-81 system, a massive protostar whose jet is known from polarized radio observations to contain relativistic electrons. Using ten years of Fermi-LAT data, the authors find a point-like excess with a test statistic above 100, a steep spectrum with photon index $\\Gamma_\\gamma = 3.53 \\pm 0.11$, and no significant variability. They argue that the excess is most plausibly produced in the HH 80-81 jet, rather than by any cataloged gamma-ray source, and that the jet has enough power to supply the observed luminosity. If correct, the result would show that protostellar jets can accelerate particles to gamma-ray-emitting energies, extending the known classes of Galactic gamma-ray sources.","feed_headline":"Fermi detects gamma rays from a protostellar jet","feed_subtitle":"A steady GeV source points to the HH 80-81 jet, powered by a massive baby star.","key_machinery":"The load-bearing object is the HH 80-81 system: a B-type protostar, IRAS 18162-2048, driving a bipolar radio jet whose knots show linearly polarized synchrotron emission, which already demonstrated the presence of relativistic electrons. The analysis machinery is Fermi-LAT binned likelihood: a model with a power-law point source is fitted against a background-only model, with the test statistic $TS = -2\\Delta \\ln L$ measuring significance; the point-like nature is checked with uniform-disk templates and the steadiness with year-binned light curves. The theoretical machinery is a set of cooling-time estimates for pp collisions, relativistic Bremsstrahlung, and inverse-Compton scattering, compared with the jet lifetime, together with the jet kinetic luminosity $L_j = \\frac{1}{2}\\dot{M}_j v_j^2$ used to show that the energy budget closes.","core_discovery":"The paper's central claim is that the HH 80-81 protostellar jet is a source of gamma rays. On ten years of Fermi-LAT Pass 8 data, a point-like excess appears at the position of the massive protostar IRAS 18162-2048 with TS $\\approx 101$ above 100 MeV; it is not extended, has a photon index $\\Gamma_\\gamma = 3.53 \\pm 0.11$, a flux above 100 MeV of $(5.2 \\pm 0.4)\\times 10^{-8}\\,\\mathrm{ph\\,cm^{-2}\\,s^{-1}}$, and is consistent with being steady. After checking the 4FGL catalog, pulsar, blazar, and radio catalogs, no known gamma-ray emitter lies within the 0.28 degree best-fit error circle. The authors conclude that the excess is likely non-thermal emission from the jet and show that both inverse-Compton scattering of relativistic electrons and proton-proton collisions can reproduce the observed spectrum, while the jet's kinetic luminosity exceeds the gamma-ray luminosity by a factor of roughly 60.","pith_inferences":["If this detection holds, the same analysis should be applied to a sample of protostellar jets with polarized radio knots; a handful of such sources would turn a single candidate into a new class of gamma-ray emitters.","The strong drop in significance from below 100 MeV to above 200 MeV suggests the signal is dominated by sub-GeV photons, where the point-spread function is broadest; repeating the analysis with the latest diffuse model and event types is the most direct check.","The paper's own comparison shows the Bremsstrahlung scenario predicts a radio synchrotron spectrum steeper than observed, leaving inverse-Compton and pp as the viable mechanisms; a joint radio-infrared-gamma-ray spectral fit could break this degeneracy.","A testable geometric prediction follows from the jet interpretation: if the emission is inverse-Compton, the gamma-ray centroid should coincide with the inner radio knots and track changes in the infrared seed photon field, so coordinated radio/IR monitoring during continued Fermi observations could discriminate the mechanisms."],"forward_implications":["If the excess is genuinely from the HH 80-81 jet, protostellar jets join supernova remnants and pulsars as Galactic gamma-ray emitters, so star-forming regions contribute to the GeV sky near the Galactic plane.","The steep spectrum and the cutoff near 1 GeV constrain the accelerated particle population: a proton spectral index $s_p \\approx 3.5$ or an electron population with a cutoff below $\\gamma_c \\sim 3\\times 10^5$ can reproduce the data.","Because the gamma-ray emission is steady and point-like at Fermi-LAT resolution, current facilities cannot resolve it; future sub-GeV instruments with sharper point-spread functions could test whether the emission traces the jet axis or a surrounding cloud.","The energy-budget result, with gamma-ray luminosity only about 0.017% of the jet kinetic power, implies the required acceleration efficiency is modest and similar protostellar jets could be detectable with longer exposures."],"supporting_citations":[{"why":"Provides the linearly polarized radio detection of the HH 80-81 jet knots, establishing the relativistic electrons that motivate the gamma-ray search.","marker":"Carrasco-González et al. 2010"},{"why":"Supplies the 4FGL source catalog used to build the background model and to check that no known gamma-ray source lies within the error circle.","marker":"The Fermi-LAT collaboration 2019"},{"why":"Defines the Pass 8 event selection and instrument response functions on which the ten-year Fermi-LAT analysis is based.","marker":"Atwood et al. 2013"},{"why":"Gives the jet physical parameters (density, mass-loss rate, velocity) adopted for the cooling-time and kinetic-power estimates.","marker":"Anglada et al. 2018"},{"why":"Supplies the theoretical framework and seed-photon values for inverse-Compton gamma-ray production in protostellar jets.","marker":"Bosch-Ramon et al. 2010"},{"why":"Provides the pp-interaction gamma-ray yield and the relation between proton and photon spectral indices used to infer ~20 GeV protons.","marker":"Kelner et al. 2006"},{"why":"Provides the radio spectral index and jet morphology used to test the bremsstrahlung/synchrotron prediction and to describe the jet's extent.","marker":"Marti et al. 1993"},{"why":"Supplies the massive-protostar jet phase timescale (~40,000 yr) against which the pp and IC cooling times are compared.","marker":"Guzmán et al. 2012"},{"why":"Identifies IRAS 18162-2048 as a massive B-type protostar, tying the gamma-ray position to the central engine of the system.","marker":"Carrasco-González et al. 2012"}],"fun_headline_variants":["Gamma rays from a protostellar jet: HH 80-81","HH 80-81 jet shines in gamma rays for Fermi","Fermi catches protostellar jet emitting gamma rays","Protostellar jet's gamma-ray emission confirmed","Steady GeV source from HH 80-81 protostellar jet"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim stands on the assumption that the standard Galactic diffuse gamma-ray background model used in the Fermi analysis, probed only by a $\\pm 3\\%$ normalization change, brackets the true background at this low-latitude position; if the true diffuse emission differs by more than that, the point-like excess could shrink, shift, or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Gamma rays from a protostellar jet: HH 80-81","HH 80-81 jet shines in gamma rays for Fermi","Fermi catches protostellar jet emitting gamma rays","Protostellar jet's gamma-ray emission confirmed","Steady GeV source from HH 80-81 protostellar jet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001013,"raw_usage":{"total_tokens":4261,"prompt_tokens":910,"completion_tokens":3351,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":3265}},"tokens_in":526,"tokens_out":3351,"duration_ms":23119,"temperature":1.0,"reasoning_tokens":3265,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:27:25.268580+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the analysis above 200 MeV with an independent diffuse background template built from gas tracers rather than the adopted standard model. The paper itself reports the excess falls to TS $\\approx 25$ above 200 MeV; if that value stays low or vanishes under the alternative background, the claimed detection would be unsupported, whereas a point-like excess that grows with more Fermi data and keeps TS above 100 would confirm the jet origin.","supporting_citations":[{"cited_title":"F., & Carrasco-González, C","cited_arxiv_id":null,"evidence_quote":"Gives the jet physical parameters (density, mass-loss rate, velocity) adopted for the cooling-time and kinetic-power estimates."},{"cited_title":"E., Araudo, A","cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical framework and seed-photon values for inverse-Compton gamma-ray production in protostellar jets."}],"review_version":1}