REVIEW 3 major objections 4 minor 6 cited by
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.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 01:13 UTC pith:GJCNWYOV
load-bearing objection A credible but model-dependent 3σ stacking excess from southern Seyferts — worth refereeing and citing, but not a discovery on its own. the 3 major comments →
Evidence for neutrino emission from X-ray Bright Seyfert Galaxies in the Southern Hemisphere using Enhanced Starting Track Events with IceCube
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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.
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.
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.
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.
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.
Where Pith is 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.
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.
Axiom & Free-Parameter Ledger
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
axioms (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}
}
read the original 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
Forward citations
Cited by 6 Pith papers
-
On the Apparent Correlation between X-ray and Neutrino Luminosities of Active Galactic Nuclei
The apparent L_nu -- L_X correlation in AGN is an artifact of TS-based selection restricting neutrino flux to a narrow range while distances span four orders of magnitude, making luminosity dominated by the distance term.
-
Particle Acceleration, Coronal Neutrino Production, and the Diffuse Extragalactic Neutrino Background from Supermassive Black Holes
The cosmologically integrated neutrino emission from supermassive black hole coronae in Seyfert galaxies can account for the sub-PeV diffuse extragalactic neutrino flux observed by IceCube.
-
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.
-
Single-source-class interpretation of the diffuse astrophysical neutrino flux
The diffuse astrophysical neutrino flux is interpreted as dominated by a single source class with dominant pγ production for target photon temperatures of 0.1-1 keV.
-
Active Galactic Nuclei as high-energy neutrino sources
AGN (jetted and non-jetted) are viable high-energy neutrino sources, with TXS 0506+056 and NGC 1068 as the only statistically significant associations so far, while gamma-ray blazars contribute only a small fraction o...
-
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...
Pith/arXiv arXiv 2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.