{"id":"335f775c-70ab-4c93-82fe-4decab5df3e9","arxiv_id":"1908.06429","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Follow-up HAWC observations of SS 433 with a new energy estimator confirm the TeV emission from the jet lobes e1 and w1.","lead":"This HAWC follow-up paper reports new TeV gamma-ray measurements of the jet interaction regions of SS 433, the microquasar, using 1,039 days of data and a new energy estimator. The measured fluxes agree with the 2018 Nature detection, confirming the lobes as point sources and supporting a leptonic origin.","discovery_kind":"replication","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The follow-up 'confirmation' is largely a re-analysis of the same 1,017 days; for w1 the new flux is 3.5e-16 versus the published 2.1e-16, and systematics for the energy-estimator results are explicitly absent.","rationale":"After reading the paper, I find the central claim plausible but under-supported. The reader's weakest_assumption correctly identifies the missing systematics and the untested comparability of the energy-estimator analysis with the published fractional-hit-bin analysis. I partly agree, and would add two sharpenings. First, the reader's strongest_claim misquotes Table 1: it lists dN/dE_w1 = 3.5e-16 as the 'published value from Abeysekara et al. 2018', whereas Table 1 lists the published w1 value as 2.1e-16 and the new energy-estimator value as 3.5e-16. This matters because the w1 agreement is the less clean of the two lobes. Second, the new dataset contains 1,017 of the same 1,039 days as ref. [4], so the apparent agreement primarily validates the new energy estimator against the old one on the same events, rather than establishing source persistence on an independent dataset. The paper itself flags the ongoing status and missing systematics, so the appropriate verdict remains CONDITIONAL rather than fully ACCEPT; there is no basis for REJECT because the claims are modest and the missing uncertainties are explicitly disclosed. The concrete test above would settle whether the on-array selection and energy-estimator systematics bias the flux normalization and whether the w1 offset is real.","tokens_in":5630,"tokens_out":6561,"duration_ms":67669,"concrete_test":"Apply the new energy estimator to the same event set used in ref. [4], including off-array events, with the same 1,017-day exposure, and fit the SS 433 region with the same multi-source model. If the resulting e1 and w1 normalizations reproduce the published 2.4 and 2.1 x 10^-16 TeV^-1 cm^-2 s^-1 values within statistical errors, the on-array-only selection is not biasing the flux comparison. If w1 shifts by roughly 1.4e-16 when the event selection is changed, the discrepancy is a selection or energy-estimator artifact. In the same pass, estimate the energy-estimator systematics by varying hadronness/containment cuts and the energy reconstruction, and report those uncertainties before claiming confirmation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the new on-array-only energy-estimator analysis reproduces the published SS 433 lobe fluxes. This requires the new analysis to be a valid, comparable measurement, and three linked conditions are insecure. First, Section 4 and the Table 1 caption state that the energy-estimator results are currently missing systematic uncertainties, so the quoted errors are purely statistical and no energy-scale or effective-area systematic is assessed. Second, the event selection differs: the new dataset contains only on-array events, while the published analysis also used off-array events; if the on-array cut has energy-dependent efficiency, the inferred normalization at 20 TeV shifts. Third, the agreement for w1 is weaker than stated: the new value is 3.5e-16 versus the published 2.1e-16 TeV^-1 cm^-2 s^-1, an offset of 1.4e-16 that is at most about 1 sigma before adding missing systematics. Additionally, 1,017 of the 1,039 days are the same data as ref. [4], so the agreement is largely a reprocessing of the same events, not new evidence that the emission persists on an independent dataset.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This ICRC 2019 proceedings paper reports a follow-up HAWC measurement of the TeV gamma-ray emission from the jet interaction regions e1 and w1 of the microquasar SS 433. Using 1,039 days of HAWC data and a new energy estimator that uses only on-array events, the authors perform a simultaneous maximum-likelihood fit of the two lobes and the nearby extended source MGRO J1908+06, with power-law spectra of fixed index 2.0 and a pivot energy of 20 TeV. The fitted flux normalizations at 20 TeV are dN/dE_e1 = (2.5^{+1.1}_{-0.8}) x 10^{-16} TeV^{-1} cm^{-2} s^{-1} and dN/dE_w1 = (3.5^{+1.2}_{-0.9}) x 10^{-16} TeV^{-1} cm^{-2} s^{-1}, compared with the published values of 2.4^{+0.6+1.3}_{-0.5-1.3} and 2.1^{+0.6+1.2}_{-0.5-1.2} from Abeysekara et al. 2018. The paper argues that this agreement supports the earlier discovery and the leptonic interpretation of the emission. Residual significance maps and pixel-significance histograms are shown to demonstrate that the fitted source model leaves a residual map consistent with background.","tokens_in":5921,"tokens_out":3024,"duration_ms":29521,"significance":"If the new measurement is taken at face value, it provides a useful cross-check of the only known VHE gamma-ray detection from a microquasar's jet termination regions, and it supports the conclusion that the lobes are point-like at HAWC angular resolution. The paper makes good use of a new, higher-resolution energy estimator that was previously validated on the Crab Nebula, and the residual maps and significance histograms provide a clear visual check that the multi-source fit does not leave obvious unmodeled excess or oversubtraction. However, the significance of the result is limited because the bulk of the data are the same as the earlier publication, the new energy-estimator results are presented without systematic uncertainties, and the w1 flux differs from the published value by more than the wording 'good agreement' implies. These issues must be addressed before the confirmation claim can be considered quantitatively established.","major_comments":[{"comment":"The Table 1 caption and Section 4 state that the energy-estimator results are currently missing systematic uncertainties. Since the central claim is that the new flux normalizations are in agreement with the published values, and the published values include large systematic uncertainties (e.g., +1.3/-1.3 on e1), a comparison without any estimate of the energy-scale or effective-area systematic uncertainty is incomplete. The authors should either provide a systematic uncertainty for the energy-estimator fluxes or explicitly state which range of agreement can be claimed before systematics are included.","section":"Section 4, Table 1"},{"comment":"The 'follow-up' dataset is not substantially independent of the earlier detection: 1,017 of the 1,039 days are the same data used in reference [4], with only 22 additional days. Thus the agreement between the two analyses is primarily a reprocessing of the same events with a different energy estimator and event selection, not new evidence that the emission persists on an independent dataset. The language in the abstract and Section 6 ('confirms', 'follow-up measurements') should be tempered, and the overlap of the datasets should be stated explicitly when the comparison is discussed.","section":"Section 3, Section 6"},{"comment":"The agreement for the w1 lobe is weaker than stated. The new flux is 3.5^{+1.2}_{-0.9} x 10^{-16} TeV^{-1} cm^{-2} s^{-1}, while the published value is 2.1^{+0.6+1.2}_{-0.5-1.2} x 10^{-16} TeV^{-1} cm^{-2} s^{-1}. The central values differ by 1.4 x 10^{-16}, which is about 1 sigma of the new statistical uncertainty and more than 2 sigma of the published statistical uncertainty. Before the missing systematic uncertainties are added, calling this 'good agreement' is an overstatement; a quantitative compatibility statement (e.g., a chi-square or p-value) is needed.","section":"Table 1, Section 5"},{"comment":"The comparison between the new and published results changes two analysis ingredients at once: the energy estimator and the event selection (on-array only versus on-array plus off-array). If the on-array selection has an energy-dependent efficiency, the flux normalization at 20 TeV could shift. The authors should either apply the new energy estimator to the same event sample used in [4], or apply the old fractional-hit-bin analysis to the on-array-only sample, to isolate the effect of each change on the fitted normalizations.","section":"Section 3"}],"minor_comments":[{"comment":"The luminosity unit in the abstract, '10^40 erg s^-2', should be 'erg s^-1' (or erg/s).","section":"Abstract"},{"comment":"There is a typo in the author list line: 'F or a complete author list' should read 'For a complete author list'.","section":"Author list"},{"comment":"The phrase 'more than 99.9% originating from cosmic rays' would read better as 'more than 99.9% of which originate from cosmic rays'.","section":"Section 2"},{"comment":"The reference list entries are not in a consistent style; for example, reference [1] mixes a 2004 publication date with a pagination that includes the journal name out of order. A uniform journal formatting would improve readability.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a short ICRC proceedings paper, so the level of technical detail expected is lower than for a full journal article. Even so, the central claim of agreement with the published SS 433 fluxes is currently not fully supported because the new results lack systematics, the w1 agreement is marginal, and the dataset is dominated by the same events as the earlier publication. The authors should be encouraged to provide the missing systematics or an explicit compatibility test and to moderate the confirmation language. The reliance on references [4] and [8] for the analysis chain is reasonable for a conference proceeding and is not a novelty concern in itself."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Unofficial take: this is exactly what it looks like — a conference proceeding that re-analyzes the same HAWC data with a newer energy estimator and confirms the earlier Nature detection. The genuinely new piece is the energy-estimator analysis, validated on the Crab in ref [8], applied to SS 433 with 1,039 days of on-array events. The paper is transparent that it is a follow-up and that the analysis chain goes back to [4] and [8]. It does not oversell the result.\n\nWhat it does well: the residual maps and significance histograms are the right checks, and they support the claim that the simultaneous fit subtracts MGRO J1908+06 and the two lobes cleanly. The e1 flux is consistent with the published value. That is real confirmation value, even if it is not an independent detection.\n\nSoft spots, in order. The energy-estimator numbers in Table 1 have no systematic uncertainties; the text says so, but that means the quoted errors are statistical only, so the \"good agreement\" claim for w1 is weaker than it looks. New w1 is 3.5 vs the published 2.1 in units of 1e-16 TeV^-1 cm^-2 s^-1 — about one sigma before adding energy-scale or effective-area systematics. Second, the new dataset uses only on-array events while the published analysis also used off-array events, so the two measurements are neither event-independent nor selection-identical. And 1,017 of the 1,039 days are the same data, so this is a reprocessing confirmation, not evidence that the emission persists on a fresh time baseline. Third, the leptonic-origin remark leans on point-source morphology; that argument was already in the Nature paper, and a point-source limit alone does not sharply discriminate leptonic from hadronic. The paper itself hedges this appropriately.\n\nOverall, the central claim holds up: the lobes are seen again in TeV gamma rays with consistent flux. The missing systematics and the selection difference should be fixed before anyone quotes these as standalone measurements.\n\nWho this is for: someone tracking HAWC analysis methods or SS 433/W50 gets a compact confirmation, but should still cite the Nature paper for the discovery. For peer review: it is a legitimate measurement worth refereeing as a proceedings or short letter, but I would not expect a main journal to treat it as new physics.","headline":"A short, honest ICRC proceeding that re-processes HAWC data with a new energy estimator and confirms the published SS 433 lobe fluxes, though it adds no new physics and its systematics are still missing.","tokens_in":6443,"tokens_out":2140,"would_cite":false,"duration_ms":24072,"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":"Follow-up HAWC measurements with 1,039 days of data and a new energy estimator reproduce the published TeV jet-lobe fluxes of SS 433, strengthening the case that the emission is real and persistent.","keywords":["SS 433","HAWC","TeV gamma-ray astronomy","microquasar jets","jet interaction regions","W 50","leptonic emission","very high energy gamma rays"],"falsifier":"Compute the energy-estimator systematic uncertainties for the e1 and w1 fits and check whether the flux normalizations at 20 TeV remain within the published values; a shift exceeding the statistical errors would overturn the confirmation claim.","tokens_in":5409,"feed_emoji":"🔭","tokens_out":6638,"duration_ms":57997,"temperature":0.7,"pith_summary":"This paper sets out to confirm the 2018 HAWC discovery that the microquasar SS 433 emits multi-TeV gamma rays from the two jet interaction regions, e1 and w1, where the jets collide with the W 50 supernova remnant about 40 parsecs from the binary. Using 1,039 days of HAWC data and a new energy estimator that uses only well-contained on-array events, the authors simultaneously fit the two lobes and the nearby extended source MGRO J1908+06. They obtain flux normalizations at 20 TeV of $2.5^{+1.1}_{-0.8} \\times 10^{-16}$ for e1 and $3.5^{+1.2}_{-0.9} \\times 10^{-16}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ for w1, in good agreement with the published values. A sympathetic reader would care because an independent data set and energy reconstruction that reproduce the original signal would confirm that the TeV emission is a persistent feature rather than a statistical fluctuation or an artifact of one analysis method.","feed_headline":"New HAWC data confirm the TeV jet lobes of SS 433","feed_subtitle":"With 1,039 days of data and a new energy estimator, the east and west jet lobe fluxes match the original discovery.","key_machinery":"The analysis is carried by a simultaneous maximum-likelihood fit of three sources: point-source models for the two SS 433 jet lobes and an electron diffusion morphology for MGRO J1908+06, the bright extended TeV source that contaminates the field. A semi-circular region of interest excludes the part of MGRO J1908+06 nearest the Galactic plane, avoiding the need to model diffuse Galactic emission. The new ingredient is the HAWC energy estimator, which uses on-array events (showers whose cores land on the main array) and is designed to be more reliable than the fractional-hit-bin method above about 10 TeV, where most photomultiplier tubes are hit. The fit uses a power-law spectrum with the index fixed at 2.0 and pivot energy 20 TeV.","core_discovery":"The central claim is that the TeV gamma-ray emission from the SS 433 jet termination regions is real and persistent, and that it remains consistent with compact, point-like sources at the X-ray-defined positions e1 and w1. The new measurement with the energy estimator yields fluxes that agree with the earlier fractional-hit-bin analysis: at 20 TeV, $dN/dE_{e1} = (2.5^{+1.1}_{-0.8}) \\times 10^{-16}$ and $dN/dE_{w1} = (3.5^{+1.2}_{-0.9}) \\times 10^{-16}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$. Because the lobes appear as point sources and do not show the wide angular spread expected if the gamma rays came from pion decay of spreading protons, the paper argues the measurement supports a leptonic production mechanism, though a hadronic-only model is not completely ruled out.","pith_inferences":["Editorial inference: If the missing systematic uncertainties on the energy estimator shift the flux normalizations by more than their statistical errors, the apparent agreement with the published values could weaken, so the confirmation is not fully closed until those systematics are reported.","Editorial inference: A direct comparison of the on-array-only and off-array event samples could reveal whether the energy-estimator selection introduces an energy-dependent exposure that partially mimics the published flux.","Editorial inference: The same simultaneous-fitting approach with an energy estimator could be applied to other compact TeV sources near the Galactic plane to test whether the point-source, leptonic interpretation for SS 433 is specific to W 50 or more general among microquasars."],"forward_implications":["The e1 and w1 lobes remain detected at TeV energies with an independent event selection and energy reconstruction, so the original discovery is not tied to the fractional-hit-bin method alone.","The agreement of the flux normalizations means the spectral energy distribution of the lobes is consistent with a power law of index 2.0 around 20 TeV.","The point-source appearance of the lobes favors in-situ acceleration of electrons at the jet termination regions and disfavors gamma rays from a widely spread hadronic proton population.","With additional cumulative data, the lobe spectra can be measured in more detail, including the systematic uncertainties currently missing from the energy-estimator results.","The residual maps after subtracting the fitted sources are consistent with background, validating the multi-source modeling of the crowded SS 433 and MGRO J1908+06 field."],"supporting_citations":[{"why":"Supplies the published HAWC discovery fluxes and the fractional-hit-bin analysis that the follow-up reproduces.","marker":"[4]"},{"why":"Describes the new HAWC energy estimator used to reconstruct the follow-up photon energies.","marker":"[8]"},{"why":"Provides the high-level analysis framework used for the maximum-likelihood fits.","marker":"[9]"},{"why":"Supplies the maximum-likelihood fitting framework used to fit the lobes and MGRO J1908+06 simultaneously.","marker":"[10]"},{"why":"Provides the electron diffusion morphology model used for MGRO J1908+06.","marker":"[11]"},{"why":"Gives the ROSAT X-ray contours that locate the jet interaction regions and support the positional identification.","marker":"[12]"},{"why":"Identifies the X-ray jet interaction regions e1, e2, e3, w1, and w2 that define where the TeV emission is expected.","marker":"[2]"}],"fun_headline_variants":["HAWC follow-up pins TeV emission to SS 433's jet lobes","SS 433's TeV jets: point-like, persistent, leptonic","HAWC data solidify TeV origin in SS 433's jet impact zones","Leptonic TeV light from SS 433's jet termination spots"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The measurement stands on the assumption that the new energy estimator and the on-array-only event selection reconstruct energies and exposure with no unaccounted energy-dependent bias large enough to change the fitted fluxes.","fun_headline_variants_meta":{"raw":{"variants":["HAWC follow-up pins TeV emission to SS 433's jet lobes","SS 433's TeV jets: point-like, persistent, leptonic","HAWC data solidify TeV origin in SS 433's jet impact zones","Leptonic TeV light from SS 433's jet termination spots"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000853,"raw_usage":{"total_tokens":3721,"prompt_tokens":976,"completion_tokens":2745,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":2663}},"tokens_in":592,"tokens_out":2745,"duration_ms":19845,"temperature":1.0,"reasoning_tokens":2663,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:45:05.349135+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the energy-estimator systematic uncertainties for the e1 and w1 fits and check whether the flux normalizations at 20 TeV remain within the published values; a shift exceeding the statistical errors would overturn the confirmation claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the published HAWC discovery fluxes and the fractional-hit-bin analysis that the follow-up reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the high-level analysis framework used for the maximum-likelihood fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the electron diffusion morphology model used for MGRO J1908+06."},{"cited_title":"(1996) ROSAT observations of the W 50/SS 433 system","cited_arxiv_id":null,"evidence_quote":"Gives the ROSAT X-ray contours that locate the jet interaction regions and support the positional identification."},{"cited_title":"and Ogelman, H","cited_arxiv_id":null,"evidence_quote":"Identifies the X-ray jet interaction regions e1, e2, e3, w1, and w2 that define where the TeV emission is expected."}],"review_version":1}