{"id":"d217ba87-aaee-4267-8108-2d767d092670","arxiv_id":"1908.08458","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"No steady neutrino excess is found from 3FHL blazars in eight years of IceCube Northern-sky data; the population contributes at most about 13-17% of the diffuse astrophysical neutrino flux, depending on spectral assumptions.","lead":"An eight-year search of IceCube data finds no steady neutrino signal from 1301 blazars, and places first upper limits on their neutrino output. The limits say these blazars can supply at most 13-17% of the diffuse astrophysical neutrino flux, narrowing the hunt for where cosmic neutrinos come from.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 13–16.7% diffuse-flux bound depends on the equal-strength stacking assumption; an inverse-acceptance weighting test would show whether it is robust.","rationale":"The strongest claim in the paper has two parts: a null result (no significant excess) and a quantitative upper limit on the blazar contribution to the diffuse flux. The null result is supported by the reported p-values and is not affected by model assumptions beyond those of the likelihood. The quantitative limit, however, rests on the equal-strength weighting and common spectral index explicitly introduced in Section 3.2. The reader's weakest_assumption identifies exactly this condition. My concern sharpens it: with W_k=1/M, the signal model is not merely sub-optimal for unequal source strengths; it can be anti-conservative for the interpretation 'total flux from the catalog', because sources with small R_k are down-weighted while high-R_k sources with no signal constrain μ_s downward. The paper's alternative-weighting bands and γ=2/2.19 variants are partial mitigation, but the set of weights is not specified and does not obviously include inverse-acceptance weights. Since the paper is transparent about the equal-strength assumption and presents the diffuse-flux fraction as conditional ('Assuming a global spectral index ...'), this concern does not invalidate the analysis; it argues for a caveat in the abstract/conclusion if the worst-case test shows a large shift. Hence the reader's ACCEPT verdict remains appropriate, and I would not change it without the re-weighting test.","tokens_in":67,"tokens_out":22141,"duration_ms":291036,"concrete_test":"Recompute the 90% C.L. upper limit on the total ν_μ+ν̄_μ energy flux from the 'All blazars' sample in Table 1 (γ=2) using stacking weights W_k∝1/R_k instead of W_k=1/M, and recalculate rmax against the diffuse flux of [10] in the 119 TeV–4.8 PeV band. If rmax exceeds the quoted 16.7% (or reaches 100%), the Section 5 bound is not robust to the equal-strength assumption; if rmax stays below 16.7%, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative statement—Section 5's claim that 3FHL blazars cannot provide more than 13.0–16.7% of the diffuse muon-neutrino energy flux for γ=2—is not a model-independent bound. It is obtained from the stacking likelihood of Eq. (3.3) with W_k=1/M (Section 3.2), i.e., equal flux at Earth for every source, and a single global spectral index. With W_k=1/M, the signal PDF weights each source by its detector acceptance R_k; low-acceptance sources (e.g., near the horizon of the Northern-sky sample) contribute almost nothing to the signal model. If the true neutrino emission is instead concentrated in a few low-acceptance sources, or has very different spectra, the likelihood can be pulled by the many high-acceptance, no-signal sources, and the resulting upper limit on the total catalog flux may not cover the true total flux. The paper partly addresses this by showing bands from 'different weighting schemes' and by quoting γ=2 and γ=2.19, but it does not state which weights are included, and the quoted 13.0–16.7% range could still miss a worst-case weighting such as W_k∝1/R_k. The null result itself is robust; the quantitative diffuse-flux fraction is conditional. A concrete worst-case re-weighting would determine whether the Section 5 conclusion overstates the constraint.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a stacking search for steady high-energy neutrino emission from 1301 blazars in the 3FHL catalog using eight years of IceCube through-going muon data from the Northern hemisphere. The analysis uses the unbinned likelihood ratio of Eqs. (3.1)-(3.3), with the signal PDF built from per-source detector acceptance and an equal physical flux per source (W_k=1/M), and a global spectral index. No significant excess is found in any of the four tested categories (all blazars, HBL, LBL/IBL, FSRQ); the most significant result is the LBL/IBL sample with a post-trial p-value of 0.03. The paper reports 90% C.L. upper limits on the neutrino flux normalization and maximum contributions to the diffuse astrophysical muon-neutrino flux, concluding that 3FHL blazars can account for at most 13.0-16.7% of the diffuse energy flux for a gamma=2 spectrum and 9.7-13.9% for gamma=2.19. It also checks that removing TXS 0506+056 does not change the LBL/IBL result and compares the steady-source limits with the TXS 0506+056 flare.","tokens_in":7363,"tokens_out":5861,"duration_ms":60900,"significance":"If the result holds, it is an important observational constraint: the collective steady neutrino emission of 3FHL blazars is well below the diffuse astrophysical muon-neutrino flux, leaving room for other source classes. The null result itself is robust: it is based on a large event sample, a standard IceCube likelihood-ratio framework, post-trial correction for the flux-threshold scan, and explicit checks of the most interesting source, TXS 0506+056. The paper also provides the first stacking limits on the 3FHL catalog, and the comparison with the TXS 0506+056 flare is a useful consistency check. The main quantitative claim, however, is conditional on the equal-strength stacking assumption W_k=1/M and on a single global spectral index; the paper partially acknowledges this through bands from different weighting schemes, but the schemes are not defined and the model dependence of the diffuse-flux fraction is understated.","major_comments":[{"comment":"The quantitative conclusion in Section 5 that 3FHL blazars 'cannot explain more than 13.0-16.7%' of the diffuse muon-neutrino flux is derived under the equal-flux stacking assumption W_k=1/M stated in Section 3.2. With that choice, the signal PDF in Eq. (3.3) weights each source by its detector acceptance R_k, so low-acceptance sources near the horizon contribute almost nothing; if the true neutrino emission is dominated by a few low-acceptance sources or by sources with very different spectra, the likelihood can be pulled by the many high-acceptance, no-signal sources, and the resulting upper limit may not bound the true total catalog flux. The caption of Fig. 3 refers to 'different weighting schemes' but does not define them, and Table 1 does not state which weights produce the quoted ranges. Please specify the weighting schemes explicitly and add a worst-case weighting such as W_k proportional to 1/R_k, or explicitly qualify the conclusion as applying only to equal-flux populations. Without this, the headline 13.0-16.7% statement overstates the constraint.","section":"Section 3.2, Eq. (3.3), Section 5, Table 1"},{"comment":"The quoted 90% C.L. upper limits and rmax values are presented without any assessment of systematic uncertainties, such as the muon energy scale, angular-resolution uncertainty, or atmospheric neutrino flux normalization. Because the headline result includes flux upper limits quoted to two significant figures, the paper should either include a systematic band or an explicit statement that the limits are statistical only and that a full treatment will appear in a dedicated analysis. This omission does not affect the null detection claim, but it does affect the precision that the reader can assign to the rmax fractions.","section":"Section 4, Table 1"}],"minor_comments":[{"comment":"The footnote states that no trial factor is applied for the four blazar categories because none of the outcomes is significant; this is not a valid statistical principle, since the multiplicity should be accounted for a priori. In this case the conclusion would not change (e.g., a Bonferroni correction gives 4x0.03=0.12), but the p-value label should be clarified.","section":"Section 4, footnote"},{"comment":"The phrase 'world largest' should be 'world's largest'.","section":"Abstract"},{"comment":"The phrase 'pre-trial corrected p-value distributions' appears to be a typo; the distributions shown are pre-trial p-values before the flux-threshold scan correction.","section":"Section 4"},{"comment":"The notation mu_s^j and N^j_tot is introduced without explicit definitions; a short sentence defining these quantities would help the reader, especially because the stacked signal PDF in Eq. (3.3) omits the sample index j.","section":"Eq. (3.1)"},{"comment":"The caption mentions a '1 sigma central band' arising from different weighting schemes, but neither the confidence level of the band nor the weighting schemes are defined; this is related to the major comment and should be specified.","section":"Figure 3 caption"}],"recommendation":"major_revision","confidential_remarks":"This is an ICRC proceedings contribution and should be evaluated as a preliminary conference report. The null search is sound and worth publishing, but the headline diffuse-flux fraction is conditional on the equal-weight stacking assumption and needs explicit qualification or a worst-case reweighting test. I do not doubt the central non-detection; the requested changes are about precision and transparency of the quantitative limits."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Matthias Huber and IceCube report a stacking search for steady neutrino emission from 1301 3FHL blazars in eight years of Northern-sky muon data. The result is a clean null: no category shows a post-trial p-value below 0.03, and the LBL/IBL excess disappears under trial correction. The upper limits on the population flux are the new numbers, and the comparison with TXS 0506+056 is done carefully. The method is not new—it's the standard IceCube unbinned likelihood stacking applied to a new catalog and longer exposure. That's fine for a proceedings; what matters is the constraint.\n\nThe main number—that 3FHL blazars can contribute at most 13.0–16.7% of the diffuse muon-neutrino energy flux for gamma=2—is conditional on the equal-strength stacking weight W_k=1/M. The paper says this, but only in passing, and the conclusion repeats the range without much caveat. More importantly, the paper never tells the reader which 'different weighting schemes' produced the shaded bands in Figure 3, so we can't judge how much the range already covers alternative assumptions. The stress-test suggestion of an inverse-acceptance weighting (W_k proportional to 1/R_k) is a legitimate robustness check; a population dominated by a few low-acceptance sources could evade the stated limit. I'd want to see that before quoting 13–17% as a firm cap.\n\nOther soft spots are minor and typical for proceedings: no released data or code, no explicit detector systematics. The trial-factor footnote is fine—even if you penalize the four categories, 0.03 times 4 is still not significant.\n\nThe citation pattern is appropriate: IceCube analysis papers, Fermi-LAT catalog papers, and the TXS papers. No red flags.\n\nBottom line: this is a workmanlike, useful constraint. Modelers of hadronic blazar emission should have it. It deserves peer review in the sense that a serious editor would send it out rather than desk-reject, but the authors should be asked to state the weighting schemes explicitly and ideally add a worst-case re-weighting. I'd cite it if I were working on blazar neutrino models.","headline":"A clean, standard stacking null for 3FHL blazars with useful limits, though the headline 13–17% diffuse-flux cap leans on an equal-strength assumption that deserves a stated robustness check.","tokens_in":7910,"tokens_out":2541,"would_cite":true,"duration_ms":25643,"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":"No steady neutrino excess from 3FHL blazars; population capped at 13–17% of the diffuse flux.","keywords":["neutrino astronomy","blazars","3FHL catalog","stacking analysis","diffuse neutrino flux","muon neutrinos","steady emission"],"falsifier":"Recompute the same stacking analysis with per-source spectral indices and gamma-ray-flux weights rather than equal weights: a post-trial significant excess in any 3FHL subsample would overturn the null result. A direct eight-year time-averaged flux measurement from any single 3FHL blazar above the population-averaged 90% C.L. flux limit at 100 TeV, at the assumed spectral index, would also contradict the equal-strength steady-emission picture; a short flare would not, because the paper's TXS comparison explicitly shows flaring sources can outshine steady limits.","tokens_in":6864,"feed_emoji":"🔭","tokens_out":9746,"duration_ms":94821,"temperature":0.7,"pith_summary":"Using eight years of northern-sky through-going muon data, this analysis stacks 1301 blazars from the third catalog of hard Fermi-LAT sources (3FHL) to look for steady neutrino emission. No excess over atmospheric background is found for the full sample or for the FSRQ, HBL, and LBL/IBL subsamples. Assuming a common power-law spectrum for all sources, the 90% C.L. upper limits cap the full blazar population at 13.0–16.7% of the diffuse astrophysical muon-neutrino energy flux for $\\gamma = 2$, and at 9.7–13.9% for $\\gamma = 2.19$. If these assumptions hold, steady emission from 3FHL blazars cannot be the dominant explanation of the diffuse astrophysical neutrino flux. The mild excess in LBL/IBL blazars is not driven by TXS 0506+056, whose known neutrino flare stands out as a time-variable event rather than the typical behavior of the population.","feed_headline":"Blazar population can’t supply most cosmic neutrinos","feed_subtitle":"Eight years of muon-track data cap blazar contribution at 13–17% of the diffuse flux.","key_machinery":"The central method is an unbinned likelihood-ratio stacking search. For each candidate source the signal probability density combines a Gaussian spatial term around the source position with an energy term that favors an $E^{-\\gamma}$ signal spectrum over the steep atmospheric background; the stacked population signal is a weighted superposition of these single-source densities, with equal source weights $W_k = 1/M$ in the main result and detector-efficiency weights $R_k$ accounting for declination- and energy-dependent response. A catalog scan orders the sources by integrated gamma-ray flux into ten subsamples, and the test statistic compares the maximum-likelihood signal hypothesis against the background-only hypothesis. This machinery turns the absence of an excess into upper limits on the population's neutrino flux and, by ratio to the measured diffuse flux, into the fractions $r_{\\mathrm{max}}$ quoted for each blazar class.","core_discovery":"The paper reports a stacking search for a steady neutrino signal from 1301 blazars in the 3FHL catalog, using roughly 500,000 through-going muon events accumulated over eight years by the South Pole neutrino telescope. None of the tested samples—all blazars, FSRQs, HBLs, or LBL/IBLs—shows a significant excess; the smallest post-trial p-value, 0.03 for LBL/IBL blazars, is treated as a null result. Converting the null into upper limits, the analysis finds that at 90% C.L. the entire 3FHL blazar population contributes at most 13.0–16.7% of the diffuse astrophysical energy flux between 119 TeV and 4.8 PeV when the population is assumed to share an $E^{-2}$ spectrum, or 9.7–13.9% when the spectrum matches the measured diffuse index $\\gamma=2.19$. Removing TXS 0506+056 barely changes the LBL/IBL result, so the most significant subsample is not dominated by that known flaring blazar. The central claim is that steady neutrino emission from these blazars is below the level required to explain the diffuse astrophysical neutrino flux.","pith_inferences":["One consequence the paper leaves implicit: if only a fraction $f$ of the 1301 blazars actually emit neutrinos, the same data allow that active subset to contribute roughly $(13\\text{--}17\\%)/f$ of the diffuse flux, so the quoted limit constrains the average source rather than every catalog member.","A natural testable extension is to run the same stacking on southern-sky cascade events, extending population coverage to the whole sky and strengthening the constraint on the blazar contribution to the diffuse flux.","The mild LBL/IBL excess at p = 0.03, absent in other classes, may hint that low-to-intermediate synchrotron-peaked blazars behave differently from HBLs; allowing per-source spectral indices in a future search would show whether the hint is physical.","If future gamma-ray catalogs push to fainter fluxes, the cumulative stacking significance pattern could reveal whether neutrino luminosity tracks gamma-ray luminosity or a separate physical property of the jets."],"forward_implications":["If the central result is correct, the 3FHL blazar population cannot produce more than 13.0–16.7% of the diffuse muon-neutrino energy flux under an $E^{-2}$ assumption, so another source class or a differently behaving blazar subpopulation must supply most of the observed astrophysical neutrinos.","The absence of a gamma-ray–neutrino correlation in the flux-threshold scans means gamma-ray brightness alone is not a reliable predictor of steady neutrino output for this catalog.","The FSRQ subsample is constrained to at most 2.9–3.8% of the diffuse flux, making FSRQs particularly unlikely steady neutrino sources.","The agreement between the population limits and the TXS 0506+056 flare studies implies that steady stacking is not sensitive to rare flaring sources, which require flare-triggered searches.","With the assumed diffuse spectral index of $\\gamma=2.19$, the maximal contribution falls to 9.7–13.9%, so the conclusion that these blazars are subdominant does not depend on the choice between the two tested spectra."],"supporting_citations":[{"why":"Supplies the 3FHL catalog of 1301 blazars and their gamma-ray fluxes used for the stacking sample and flux-threshold partitions.","marker":"[5]"},{"why":"Defines the IceCube detector and the eight-year through-going muon data set used in the analysis.","marker":"[6]"},{"why":"Establishes the subdegree median angular resolution of muon events above 10 TeV used in the spatial signal probability density.","marker":"[7]"},{"why":"Provides the diffuse astrophysical muon-neutrino flux spectrum against which the maximum blazar contribution r_max is evaluated.","marker":"[10]"},{"why":"Supplies the HESE-based astrophysical flux estimate used for comparison in the flux-limit figures.","marker":"[11]"},{"why":"Gives the TXS 0506+056 steady and 158-day flaring flux bands that the 3FHL limits are checked against.","marker":"[2]"},{"why":"Documents the 2017 multimessenger detection of TXS 0506+056 that motivates the blazar search.","marker":"[1]"},{"why":"Motivates the association between high-energy neutrinos and high-synchrotron-peaked BL Lacs that informs the source-class choice.","marker":"[4]"}],"fun_headline_variants":["Blazar steady neutrinos capped at 17% of diffuse flux","No steady neutrino excess from 1301 blazars in 8 years","Blazars not main source of cosmic neutrinos, limits show","IceCube finds blazars contribute at most 17% to neutrino sky","Steady blazar neutrino emission below diffuse flux levels"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fixed upper limits assume that every 3FHL blazar produces the same neutrino flux at Earth and that the whole population shares one power-law spectral index; if the true neutrino emission is concentrated in a few sources with different spectra, these population limits do not directly constrain that concentrated emission.","fun_headline_variants_meta":{"raw":{"variants":["Blazar steady neutrinos capped at 17% of diffuse flux","No steady neutrino excess from 1301 blazars in 8 years","Blazars not main source of cosmic neutrinos, limits show","IceCube finds blazars contribute at most 17% to neutrino sky","Steady blazar neutrino emission below diffuse flux levels"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000252,"raw_usage":{"total_tokens":1608,"prompt_tokens":1042,"completion_tokens":566,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":658,"completion_tokens_details":{"reasoning_tokens":474}},"tokens_in":658,"tokens_out":566,"duration_ms":6193,"temperature":1.0,"reasoning_tokens":474,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:39:04.778278+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the same stacking analysis with per-source spectral indices and gamma-ray-flux weights rather than equal weights: a post-trial significant excess in any 3FHL subsample would overturn the null result. A direct eight-year time-averaged flux measurement from any single 3FHL blazar above the population-averaged 90% C.L. flux limit at 100 TeV, at the assumed spectral index, would also contradict the equal-strength steady-emission picture; a short flare would not, because the paper's TXS comparison explicitly shows flaring sources can outshine steady limits.","supporting_citations":[{"cited_title":"Padovani and E","cited_arxiv_id":null,"evidence_quote":"Supplies the 3FHL catalog of 1301 blazars and their gamma-ray fluxes used for the stacking sample and flux-threshold partitions."},{"cited_title":"Ajello et al., Astrophys","cited_arxiv_id":null,"evidence_quote":"Defines the IceCube detector and the eight-year through-going muon data set used in the analysis."},{"cited_title":"Padovani, E","cited_arxiv_id":null,"evidence_quote":"Provides the diffuse astrophysical muon-neutrino flux spectrum against which the maximum blazar contribution r_max is evaluated."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the HESE-based astrophysical flux estimate used for comparison in the flux-limit figures."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the TXS 0506+056 steady and 158-day flaring flux bands that the 3FHL limits are checked against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Motivates the association between high-energy neutrinos and high-synchrotron-peaked BL Lacs that informs the source-class choice."}],"review_version":1}