REVIEW 3 major objections 4 minor 3 cited by
Evidence for neutrino emission from X-ray Bright Seyfert Galaxies in the Southern Hemisphere using Enhanced Starting Track Events with IceCube
T0 review · 3 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read This paper reports a 3.0σ excess of neutrino events from a stacking analysis of 13 southern Seyfert galaxies, providing new evidence that Seyferts are sources of high-energy neutrinos.
desk verdict A credible but model-dependent 3σ stacking excess from southern Seyferts — worth refereeing and citing, but not a discovery on its own. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The analysis uses the Enhanced Starting Track Event Selection (ESTES), a dataset of muon-neutrino events with interaction vertices inside the detector, which suppresses atmospheric muons and achieves ~1.4° angular resolution for Southern-sky sources. An unbinned likelihood ratio test weights each candidate source by its predicted neutrino flux from the disk-corona model, which assumes cosmic-ray acceleration in the AGN corona and is calibrated to match NGC 1068. A hybrid background—scrambled data plus Monte Carlo Galactic-plane events—accounts for recently detected Galactic neutrino emission and avoids overestimating background near bright sources.
What would settle it
A concrete check: if the 6.7 excess events do not cluster around the three predicted-brightest sources (Circinus, NGC 7582, ESO 138-1) with the predicted spectral shape, but instead are spread across the sky, the background model is likely wrong. With doubled exposure, the excess should grow toward ~5σ if real; a decline would indicate a fluctuation.
Extended reading notes
Core claim
The central claim is that the collective neutrino signal from 13 southern Seyfert galaxies (excluding Centaurus A) is real, with a best-fit of 6.7 events (single-sided p = 1.3×10⁻³, corresponding to 3.0σ). This is the first stacking evidence from a Southern-Hemisphere starting-track sample and is consistent with the disk-corona model prediction of 4.7 events. The excess is driven mainly by Circinus, NGC 7582, and ESO 138-1, whose best-fit counts correlate with the model predictions.
Load-bearing premise
The analysis assumes the disk-corona model, with parameters fixed by NGC 1068, correctly predicts the relative neutrino fluxes of all 14 Seyferts from their 2–10 keV X-ray luminosities; if this scaling is wrong, the stacking weights are suboptimal and the 3.0σ excess may not indicate true Seyfert emission.
Editorial extensions
If this is right
- If the excess is real, Seyfert galaxies collectively contribute to the observed diffuse neutrino flux, not just NGC 1068.
- The agreement with the disk-corona prediction supports emission from an optically thick corona, explaining the missing gamma-ray counterpart.
- The result motivates dedicated searches for individual Seyferts, which could reach discovery significance with additional exposure or future detectors.
- Weighting sources by X-ray luminosity provides a general method for targeting other AGN populations in neutrino searches.
Reading between the lines
- The 3.0σ significance is below the standard 5σ discovery threshold; if the model's relative flux predictions are partly incorrect, the true source population could be more concentrated (e.g., only Circinus and ESO 138-1 contributing) than the stacking assumes.
- A concrete next test is to add cascade-type events, which have complementary sensitivity in the Southern sky; a combined analysis would either strengthen the excess or reveal inconsistencies.
- The apparent correlation between best-fit neutrino counts and disk-corona predictions for the three brightest sources could be a statistical coincidence; checking whether the excess grows proportionally with exposure is a decisive follow-up.
- The exclusion of Centaurus A, justified by its jet ambiguity, leaves open whether jet-associated emission also produces neutrinos; a separate, jet-model-based stacking could test that.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper searches for high-energy neutrino emission from 14 X-ray bright Seyfert galaxies in the Southern Hemisphere using roughly ten years of IceCube Enhanced Starting Track Events (ESTES). Individual-source searches, assuming either a power law or the disk-corona spectrum, find no source that survives trial corrections; the best individual candidate, the Circinus Galaxy, has a post-trial significance of only 1.8σ. The main result is a stacking analysis in which each source is weighted by the number of events predicted by a disk-corona model normalized to NGC 1068. Excluding Centaurus A, the stacking analysis finds a best-fit excess of 6.7 signal events over background, with a single-sided p-value of 1.3×10⁻³ (3.0σ), which the authors interpret as new independent evidence that Seyfert galaxies contribute to the extragalactic neutrino flux.
Significance. If the result holds, it is a valuable addition to the growing case for Seyfert galaxies as neutrino emitters, extending the NGC 1068 discovery to a population of southern X-ray bright AGNs and demonstrating the power of the ESTES starting-track sample. The paper has real strengths: the likelihood formalism is described in detail, the background model combines data scrambling with Monte Carlo modeling of the Galactic plane, the analysis was validated before unblinding, and the authors show that a more traditional background treatment yields a reduced but still nonzero significance of 2.8σ. The main weakness is that the quoted significance is tied to a specific model weighting whose normalization is calibrated to NGC 1068, and no model-independent collective test is provided for the southern sample.
major comments (3)
- [Secs. 4–5, Table 2] The 3.0σ stacking significance is computed for a signal hypothesis whose source weights are the disk-corona model predictions, with the model normalization fixed to the NGC 1068 flux (Sec. 1). The p-value therefore rejects background for this particular weighted combination, but it does not by itself establish that the selected Seyfert galaxies are neutrino sources if the X-ray-to-neutrino scaling is incorrect. The paper reports no model-independent collective test for the southern sample, unlike the binomial test used in the northern-sky analysis (Abbasi et al. 2025a). I recommend adding such a test (e.g., an unweighted stack or a binomial test of excess directions) or explicitly limiting the abstract's claim to the model-weighted stacking result.
- [Sec. 5] The statement that the best-fit 6.7 events are 'consistent with the model prediction of 4.7 events' is presented as support for the model, but this consistency is weak and partly by construction: the model normalization was calibrated to NGC 1068, and the quoted 68% interval (3.5–10.7) is broad. Please quantify the systematic uncertainty in the predicted event rate arising from the assumed X-ray–luminosity scaling and avoid presenting the 4.7-event prediction as an independent confirmation.
- [Secs. 2 and 5] The source list is selected from the Swift-BAT/BASS catalog using a declination cut and a threshold of 0.1 predicted signal events per source, both based on the same disk-corona model. The stacking p-value is derived for this fixed list, but the paper does not state whether the threshold and source selection were fixed before unblinding or how they affect the trials. If the list or threshold was chosen after inspecting model predictions, a trials correction may be needed; if it was fixed a priori, this should be stated explicitly.
minor comments (4)
- [Table 2] The p-value entries appear garbled: for example, '−0.003' should read '0.003' (and similarly for other rows), and 'ES 138-1' in the table note is a typo for 'ESO 138-1'.
- [Sec. 2] The declination cut 'δ < −5°' and the phrase 'the horizontal region' would benefit from a brief explanation of the geometric/background motivation, since the latter term is not defined in the text.
- [Secs. 5–6] The text says the stacking excess is 'driven by five galaxies–led by three nearby AGNs,' but the individual best-fit counts in Table 4 for the three leading sources (3.6, 1.7, 1.1) already sum to 6.4, and adding the two additional sources gives more than the stacked best fit. Clarify whether Table 4 lists individual-source fits rather than components of the stacking fit.
- [Abstract] The abstract states the result as 'new independent evidence that Seyfert galaxies contribute to the extragalactic flux of high-energy neutrinos.' Given the model-dependent weighting, I suggest adding a brief qualifier such as 'under disk-corona model assumptions' to accurately represent what is tested.
Circularity Check
No significant circularity; the 3.0σ stacking p-value is a background-only rejection and does not reduce to the NGC 1068-calibrated model prediction.
full rationale
The central 3.0σ stacking excess is evaluated with an unbinned likelihood-ratio test against background-only hybrid data+MC datasets constructed by scrambling the data (Sec. 4). The null p-value of 1.3×10−3 is therefore a measured rejection of the background-only hypothesis, not a fitted quantity; it does not reduce to the disk-corona model normalization. The model enters only through the relative signal weights (n_exp) and the assumed spectral shape, and the best-fit n_s=6.7 is a free parameter compared with the model prediction of 4.7 events. That comparison is a consistency check, not the load-bearing evidence. The model parameters were calibrated to NGC 1068 (Sec. 1), so the 4.7-event prediction is not parameter-free, but this is a model-dependence caveat rather than a circular reduction: the southern stacking data are not used to fit the model, and the absolute normalization cancels in the normalized signal PDF. Self-citations to the ESTES selection and prior AGN searches are used as experimental tools and context, not as a uniqueness theorem or as a substitute for the background-only test. The interpretation that the excess comes from Seyfert galaxies is conditional on the disk-corona X-ray-to-neutrino scaling; a wrong scaling could degrade sensitivity but would not make the background rejection circular. No step in the derivation is equivalent by construction to its input.
Assumptions & free parameters
free parameters (2)
- Disk-corona model normalization (proton-to-X-ray luminosity scaling) =
Calibrated to NGC 1068 observed neutrino flux (value not quoted)
- Source selection threshold in predicted signal events =
0.1 events
assumptions (5)
- domain assumption Neutrino flux of a Seyfert is proportional to its intrinsic 2–10 keV X-ray luminosity (disk-corona model)
- domain assumption The neutrino spectrum and efficiency are the same for all candidate Seyferts as for NGC 1068
- domain assumption Galactic plane neutrino background is described by the Fermi-LAT pi0 template with E^-2.7 spectrum and normalization 21.8e-12 TeV cm^-2 s^-1 at 100 TeV
- domain assumption After masking and resampling, RA randomization gives a valid background PDF
- domain assumption ESTES reconstruction and detector simulation correctly model events
Cite this review
Pith. "Pith review of Evidence for neutrino emission from X-ray Bright Seyfert Galaxies in the Southern Hemisphere using Enhanced Starting Track Events with IceCube." pith.science (2026). https://pith.science/paper/GJCNWYOV
@misc{pith2026260210208,
author = {Pith},
title = {Pith review of: Evidence for neutrino emission from X-ray Bright Seyfert Galaxies in the Southern Hemisphere using Enhanced Starting Track Events with IceCube},
year = {2026},
howpublished = {\url{https://pith.science/paper/GJCNWYOV}},
note = {Machine review of arXiv:2602.10208}
}
abstract
IceCube recently reported the observation of TeV neutrinos from the nearby Seyfert galaxy NGC~1068, and the corresponding neutrino flux is significantly higher than the upper limit implied by observations of GeV-TeV gamma rays. This suggests that neutrinos are produced near the supermassive black hole, where the radiation density is high enough to obscure gamma rays. We use a set of muon neutrinos with interaction vertices inside the detector, which have good sensitivity to sources in the Southern sky, from IceCube data recorded between 2011 and 2021. We then search for individual and collective neutrino signals from 14 Seyfert galaxies in the Southern Sky selected from the Swift Burst Alert Telescope (BAT) AGN Spectroscopic Survey. Using the correlations between keV X-rays and TeV neutrinos predicted by disk-corona models, and assuming production characteristics similar to NGC~1068, a collective neutrino signal search reveals an excess of $6.7_{-3.2}^{+4.0}$ events, which is inconsistent with background expectations at the 3$\sigma$ level of significance. In this paper, we present new independent evidence that Seyfert galaxies contribute to the extragalactic flux of high-energy neutrinos.
Figures
Figures from the paper (1 more)
Forward citations
Cited by 3 Pith papers
-
Measurement of the Hubble constant with high-energy neutrinos
Using 12 neutrino-emitting Seyfert galaxies and the Cepheid-plus-TRGB distance to NGC 1068, the authors obtain the first neutrino-based Hubble constant estimate, H0 = 49+40-30 km/s/Mpc, with slope beta = 0.67+0.16-0.25.
-
On the Blueprint of Active Galaxies Producing Neutrinos
Neutrinos from active galaxies are produced in compact X-ray-bright coronae within about ten Schwarzschild radii of the black hole, and such sources may supply the diffuse neutrino flux.
-
Active Galactic Nuclei as high-energy neutrino sources
Review summarizing correlations between astrophysical neutrinos and gamma-ray/radio AGN, specific source associations, proposed production mechanisms, and future prospects.
Reference graph
Works this paper leans on
-
[1]
2013, Science, 342, 1242856, doi: 10.1126/science.1242856 Aartsen, M
Aartsen, M., et al. 2013, Science, 342, 1242856, doi: 10.1126/science.1242856 Aartsen, M. G., et al. 2016, Astropart. Phys., 78, 1, doi: 10.1016/j.astropartphys.2016.01.006 —. 2017, JINST, 12, P03012, doi: 10.1088/1748-0221/12/03/P03012 Aartsen, M. G., Ackermann, M., Adams, J., et al. 2020, PhRvL, 124, 051103, doi: 10.1103/PhysRevLett.124.051103 Aartsen, ...
arXiv 2013
-
[2020]
https://arxiv.org/abs/2005.11208 Bechtol, K., Ahlers, M., Di Mauro, M., Ajello, M., & Vandenbroucke, J. 2017, Astrophys. J., 836, 47, doi: 10.3847/1538-4357/836/1/47 Braun, J., Dumm, J., De Palma, F., et al. 2008, Astropart. Phys., 29, 299, doi: 10.1016/j.astropartphys.2008.02.007 Capanema, A., Esmaili, A., & Murase, K. 2020, Phys. Rev. D, 101, 103012, do...
arXiv 2005
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.