Searching for neutrino emission from hard X-ray sources with IceCube

Marcos Santander (for the IceCube Collaboration)

REVIEW 3 major objections 4 minor 16 references

Searching for neutrino emission from hard X-ray sources with IceCube

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read IceCube plans a stacked neutrino search across 828 hard X-ray AGN

desk verdict An honest ICRC proceedings plan for a BASS stacking search, well written but with no data, no sensitivity, and a physical premise that may not transfer from blazars to Seyferts. read the letter →

arxiv 1908.04862 v1 pith:2U6JLEYE submitted 2019-08-13 astro-ph.HE

Marcos Santander (for the IceCube Collaboration) This is my paper · ORCID

classification astro-ph.HE

keywords neutrinoastronomyIceCubeactivegalacticnucleihardX-raysourcesBATAGNSpectroscopicSurveystackinganalysisblazars

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper proposes a search for high-energy neutrinos from the 828 AGN in the BAT AGN Spectroscopic Survey (BASS), using eight years of IceCube muon-neutrino data. Its rationale is that if blazars are neutrino sources, the gamma rays produced alongside neutrinos in hadronic processes cascade in the dense photon fields of the jet and re-emerge as strong hard X-ray to MeV emission, so a hard X-ray-selected catalog should be rich in neutrino emitters. The search stacks all sources in a time-integrated unbinned maximum-likelihood analysis, with five weighting schemes that depend on X-ray flux, luminosity, spectral index, column density, or equal weight. The paper is a plan: it establishes the sample, the method, and the weight definitions, and reports that the catalog's overlap with gamma-ray-selected searches is small, meaning a genuinely new source class is being tested. No neutrino signal or limit is claimed.

What carries the argument

The load-bearing object is the BASS catalog: 828 hard X-ray AGN from the 70-month Swift/BAT survey with measured redshifts, which enters the search as the target list. The analysis machinery is a time-integrated unbinned maximum-likelihood stack, with source probability $S_i = \sum_k \omega_k R_k(\delta_k,\gamma) S(\vec{x}_i,\vec{x}_k,\gamma)/\sum_k \omega_k R_k(\delta_k,\gamma)$, where $\omega_k$ is the per-source weight, $R_k$ the detector acceptance at the source declination, and $\gamma$ a common spectral index. The five weighting schemes — intrinsic 14-195 keV flux, isotropic-equivalent luminosity, inverse X-ray spectral index, total column density $n_H$, and equal weight — are what translate hard X-ray properties into neutrino-search priorities. The hadronic cascade argument, illustrated with TXS 0506+056, supplies the physical reason to expect this translation to work.

What would settle it

Run the identical stacked analysis on the full BASS sample and on a matched control catalog selected by, say, infrared or soft X-ray luminosity; if the control catalogs produce equal or larger test statistics after accounting for trial factors, the hard X-ray cascade premise is not doing the selection work.

Watch

Extended reading notes

Core claim

The central claim is that the BASS hard X-ray AGN catalog is a promising and untested target list for IceCube neutrino searches, and that the proposed stacked analysis will be able to test it. The paper argues that previous AGN searches relied on gamma-ray correlations, while BASS selects sources by 14-195 keV emission from Swift/BAT, a band where hadronic cascades should contribute for neutrino-producing blazars. It defines a source term in an unbinned likelihood that weights each of the 828 AGN by one of five observables and folds in the detector's declination-dependent acceptance, and a test statistic evaluated against scrambled data. The stated near-term deliverable is sensitivities and discovery potentials, followed, in the absence of a signal, by constraints on this AGN class as contributors to the all-sky astrophysical neutrino flux.

Load-bearing premise

The search only improves on gamma-ray-selected searches if neutrino-producing hadronic processes in AGN also make gamma rays that cascade into strong hard X-ray to MeV emission; if that cascade is rare or swamped by other X-ray emission, the BASS sample will not preferentially contain neutrino emitters.

Editorial extensions

If this is right

If hard X-ray selection works, the search can uncover neutrino sources that gamma-ray-selected catalogs miss, since only 7% of BASS sources have a 4FGL counterpart within localization uncertainty.
The five weighting schemes let the analysis determine empirically which X-ray property — flux, luminosity, spectral hardness, or obscuration — best tracks neutrino emission.
With no significant signal, the analysis will set upper limits on the contribution of BASS AGN to the observed all-sky astrophysical neutrino flux.
A detection would extend the population of confirmed neutrino emitters beyond TXS 0506+056 and directly support hadronic cascade models in blazar jets.
The catalog's low redshift (median z ≈ 0.04) means the search is sensitive to the nearest, brightest AGN, where individual source associations are most tractable.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If hard X-ray emission is a valid neutrino tracer, then a population of obscured or Compton-thick AGN — underrepresented in Fermi gamma-ray catalogs — could be a substantial part of the all-sky neutrino flux; the column-density weighting is the first IceCube test of that idea.
The same stacking framework could be applied to other hard X-ray samples, such as INTEGRAL/IBIS or NuSTAR serendipitous surveys, to cross-check that any signal is generic to hard X-ray selection rather than specific to BASS.
A positive correlation with X-ray luminosity would imply a diffuse component of MeV gamma rays from cascades across the AGN population, a prediction that next-generation MeV telescopes could test.
The five single-observable weights could be combined into a joint likelihood or ranking score that uses all four X-ray properties at once, which would likely outperform any single weighting scheme.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

This is a conference proceedings describing a planned IceCube stacking search against 828 hard X-ray AGN from BASS. The paper is exactly what it claims to be: a plan. There is no event sample, no likelihood scan, no limit, no sensitivity estimate. The text says so itself in Section 4. So judge it as a plan, not as a result.

What is genuinely useful: the BASS selection is a new source list for IceCube stacking, and the overlap study in Section 2.1 is a real check. Only 10% of BASS sources have a 4FGL counterpart within 0.2 degrees, so this is not simply a re-run of the Fermi-AGN analyses. The five weighting schemes are clearly defined from BASS quantities, and the unbinned likelihood in Eqs. 3.1–3.3 is written down correctly. The description of the method is clean enough that someone could implement it.

The soft spot is the physical motivation. The paper leans on the hadronic-cascade picture from TXS 0506+056 and Petropoulou et al./Murase et al. to argue that hard X-ray selection preferentially picks neutrino emitters. But BASS is 89% non-beamed AGN, mostly nearby Seyferts whose 14-195 keV emission is usually coronal or disk reprocessed, not jet cascade emission. That does not make the search worthless — it makes the premise an untested assumption. The paper never quantifies the cascade contribution across the actual catalog, so the claim that BASS is a novel source class is plausible but unsupported. Equal weighting does not rescue this, because the issue is at the catalog level, not the weight assignment.

A second, more obvious limitation is the total absence of sensitivity or discovery-potential estimates. For a plan, I can live with that, but it means the paper has no scientific output to scrutinize beyond the internal consistency of the proposal.

Bottom line: this is a fine ICRC proceedings writeup. It is honest, methodologically sound, and the BASS–4FGL overlap check is a small but real contribution. It does not belong in a journal as a standalone paper, because there is no measurement or simulation to review. If the completed analysis with sensitivities and limits comes out, that would be worth a serious referee. For this document, peer review would be an exercise in commenting on a proposal, not on a result.

Referee Report

3 major / 4 minor

Summary. The paper proposes a time-integrated stacked muon-neutrino search with eight years of IceCube data, targeting the 828 AGN in the BAT AGN Spectroscopic Survey (BASS) DR-1. It outlines the unbinned maximum-likelihood formalism (Eqs. 3.1-3.3), the source PDF with per-source weighting, five weighting schemes (flux, luminosity, spectral index, column density, equal), the overlap of BASS with Fermi-4FGL, and the planned test-statistic calibration via coordinate scrambling. The paper explicitly presents this as a plan and states that sensitivities and discovery potentials are future work; no neutrino signal, limit, or validation result is claimed.

Significance. If the proposed search is carried out and the underlying physical premise holds, it would provide a novel test of hadronic neutrino emission using a hard-X-ray-selected AGN sample, complementary to the gamma-ray-selected stacking searches that have so far yielded null results. The technical core is sound: the likelihood, source PDF, and scrambling-based test statistic follow IceCube's standard stacked-search formalism, and the weighting schemes are consistently defined from BASS DR-1 quantities. However, as submitted, the paper contains no measurements, no sensitivity projections, and no validation of the analysis; its current significance is therefore limited to the proposal itself. The motivation linking hard X-ray selection to neutrino emission is plausible but not quantified for the actual BASS catalog, which is dominated by non-beamed Seyferts.

major comments (3)

[Section 1 and Figure 1] The central motivation assumes that hard X-ray (14-195 keV) selection, as embodied by BASS, preferentially identifies neutrino emitters because hadronic cascades produce strong hard-X-ray/MeV emission. This premise is illustrated with a flaring-blazar model (TXS 0506+056), but BASS itself is roughly 89% non-beamed AGN, mostly Seyferts with median redshift z~0.04, whose hard X-ray emission is generally attributed to coronal or disk reprocessing rather than jet hadronic cascades. The manuscript never quantifies the hadronic-cascade contribution to the BAT band across the catalog, nor does it provide evidence that sources lacking such a component would be neutrino-quiet. Because this premise sits at catalog selection rather than weight assignment, the equal-weight scheme does not hedge the risk. Please add a quantitative estimate or at least a detailed, catalog-specific discussion of the expected cascade contribution, or otherwise the claim that this search is well-motivated over gamma-ray-selected stacks is unsupported.
[Section 4] The paper explicitly defers sensitivity and discovery-potential estimates to 'next steps'. Without pseudo-experiments, background-only scrambling studies, or signal-injection tests, the manuscript provides no quantitative evidence that the proposed analysis can actually constrain or detect neutrino emission from this source class. The feasibility of the plan cannot be assessed from the current text. At minimum, the paper should include expected sensitivity curves (e.g., the 90% CL upper limit on the stacked flux normalization as a function of spectral index) and discovery potential for representative assumptions about the BASS source population, using the described likelihood and selection.
[Section 3.3] The five weighting schemes are described, but the manuscript does not specify how the five resulting test-statistic values will be handled for multiple testing. Since the analysis will maximize the likelihood ratio for each weighting scheme (and for the continuous spectral index), the final significance must account for the look-elsewhere effect over the five schemes. The paper should state a concrete procedure, such as calibrating the final p-value by repeating the entire five-scheme maximization on scrambled data, or using a single pre-defined primary weighting scheme.

minor comments (4)

[Section 3.1, text below Eq. (3.1)] There is a typo: 'atmoshpheric' should be 'atmospheric'.
[Section 3.3, Figure 5] Figure 5 shows cumulative distributions of luminosity, flux, spectral index, and column density, but the text refers to them as the 'weights'. Please clarify how these distributions translate into the actual ω_k values (e.g., whether the weights are the raw quantities or normalized versions).
[Section 2.1] The statement that the overlap with previous Fermi-based studies is 'not significant' is based on the 7% source overlap and a qualitative comparison of flux distributions in Figure 4. A quantitative statement (e.g., a Kolmogorov-Smirnov test on the flux distributions) would make this claim more robust.
[References] Reference [10] is incomplete as written ('These proceedings 916 (2019)'); it should include the author, title, and full proceedings information.

Circularity Check

0 steps flagged · score 0.0 of 10

The paper is a plan for a future stacking search; it fits no parameters, reports no signal, and makes no derivation whose output reduces to its input.

full rationale

No significant circularity was found. The paper explicitly presents plans for a time-integrated stacked neutrino search using 8 years of IceCube muon-neutrino data against 828 BASS hard X-ray AGN; it does not claim a detection, a limit, or a derived physical quantity. The likelihood in Eq. 3.1 is the standard IceCube unbinned likelihood, the source PDF in Eq. 3.2 combines externally chosen weights with detector acceptance, and the test statistic in Eq. 3.3 is evaluated against a scrambled background distribution. The five weighting schemes are defined directly from BASS catalog observables (flux, luminosity, spectral index, column density, and equal weights) and none is fitted to the neutrino data or to the search outcome. The physical motivation that hadronic cascades produce hard X-ray to MeV emission, supported by the TXS 0506+056 modeling in Figure 1, is indeed an assumption and could be challenged on astrophysical grounds for mostly non-beamed Seyfert sources, but an unquantified or debatable premise is not circular reasoning. The overlap estimate with Fermi 4FGL is a computed catalog comparison, not a self-referential derivation. Citations to IceCube collaboration papers are used for standard detector description, muon-neutrino performance, and previous search results; the cited TXS hadronic modeling is from an independent group (Gao et al.) and is not used to forbid alternatives. No equation is shown to equal its own input, no fitted parameter is renamed as a prediction, and no self-citation chain carries the central argument. The central claim is only that the analysis is feasible and motivated, and that claim is self-contained and consistent with the standard stacked-search formalism.

Assumptions & free parameters 1 free parameters · 4 assumptions · 0 invented entities

The paper reports no fitted values; the only analysis parameter, the spectral index gamma, is to be optimized in the future search. The plan relies on the hadronic cascade assumption, the accuracy of the BASS catalog, standard IceCube event properties, and a standard cosmology. No new particles, forces, or entities are introduced.

free parameters (1)

Spectral index gamma (proposed analysis) = to be determined by maximization of the test statistic in Eq. 3.3
The likelihood in Eq. 3.1 assumes a power-law neutrino spectrum with index gamma; this is a free parameter of the planned search, not fixed in the paper.

assumptions (4)

domain assumption Hadronic neutrino production in AGN is accompanied by gamma rays that cascade in the source and re-emerge as hard X-ray to MeV emission.
Used in Section 1 to justify selecting BASS hard X-ray sources; Figure 1 illustrates the hadronic model for TXS 0506+056.
domain assumption The BASS DR-1 catalog is a reliable, complete all-sky AGN sample in the 14-195 keV band.
The entire source list of 828 AGN is taken from BASS (Ref. [15]); the analysis assumes the catalog positions, fluxes, luminosities, spectral indices, and column densities are accurate.
domain assumption IceCube muon-neutrino events have an angular resolution of about 1 degree above 10 TeV and an atmospheric background spectrum with index ~3.7.
Stated in Section 3.1 and used in the unbinned likelihood (Eq. 3.1) to separate signal from background.
domain assumption Standard flat Lambda-CDM cosmology with Omega_L=0.7, Omega_M=0.3, H0=70 km/s/Mpc is used for luminosity distances.
Adopted in Section 3.3 for the luminosity weighting scheme.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Searching for neutrino emission from hard X-ray sources with IceCube." pith.science (2026). https://pith.science/paper/2U6JLEYE

@misc{pith2026190804862,
  author       = {Pith},
  title        = {Pith review of: Searching for neutrino emission from hard X-ray sources with IceCube},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2U6JLEYE}},
  note         = {Machine review of arXiv:1908.04862}
}

read the original abstract

The IceCube neutrino observatory, a cubic-kilometer particle detector at the South Pole, first announced the discovery of an astrophysical flux of high-energy neutrinos in the TeV-PeV range in 2013, followed in 2017 by the detection of a high-energy neutrino event in temporal and directional correlation with the flaring gamma-ray blazar TXS 0506+056. This observation, combined with archival neutrino detections in 2014-2015, has provided compelling evidence for the detection of the first high-energy astrophysical neutrino source. A promising way of detecting additional sources is to correlate neutrino detections with sources where a hadronic electromagnetic signature is observed. If blazars are a significant source of neutrinos, the high-energy gamma rays produced in pionic decays in coincidence with the neutrinos may cascade in the strong photons fields present in blazar jets, leading to strong emission in the hard X-ray to MeV gamma-ray energy range. We here present plans for a search for neutrino emission from a large sample of hard X-ray sources from the BAT AGN Spectroscopic Survey (BASS).

Figures

Figures reproduced from arXiv: 1908.04862 by the authors.

**Figure 1.** Modeling of the spectral energy distribution of the blazar TXS 0506+056 during the 2017 flare potentially associated with a high-energy neutrino for purely leptonic (a) and purely hadronic (b) models. The hadronic models tend to overproduce flux in the hard X-ray to MeV range due to hadronic pair production and cascading effects (from [7]). The angular distribution of the astrophysical neutrinos discovered by IceCu… view at source ↗

**Figure 2.** Skymap in equatorial coordinates of the BASS sources that will be used in the proposed search. The dashed line indicates the Galactic plane, with a black marker showing the Galactic Center position. Marker colors show the intrinsic (i.e. deabsorbed) energy flux of the source in the 14-195 keV band in logarithmic scale. Most BASS sources are nearby, with a median of z ∼ 0.04 and 90% of the sources located within a re… view at source ↗

**Figure 3.** Left: Redshift distribution of the BASS catalog sources. Right: Hard X-ray luminosity (isotropicequivalent) for the BASS sources as a function of redshift. Different markers indicate AGN classes, while the colors indicate the derived intrinsic luminosity in the 14-195 keV band. The highest luminosities are mostly associated with blazars thanks to the contribution of Doppler boosting. AGN with a potential gamma-ray … view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: Gamma-ray energy fluxes above 100 MeV for all 4FGL sources (light blue) and those included in BASS (purple). The BASS gamma-ray sources have a slightly higher energy flux. the Cherenkov light emitted by relativistic charged particles produced in high-energy neutrino in…

**Figure 5.** Figure 5: Luminosity (top left), intrinsic flux (top right), spectral index (bottom left) and column density (bottom right) weights based on the BASS catalog values. References [1] S. R. Kelner, F. A. Aharonian, and V. V. Bugayov, Phys. Rev. D74 (2006) 034018. [Erratum: Phys. Re…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 12 canonical work pages

[1]

S. R. Kelner, F. A. Aharonian, and V . V . Bugayov,Phys. Rev. D74 (2006) 034018. [Erratum: Phys. Rev.D79,039901(2009)]

work page 2006
[2]

S. R. Kelner and F. A. Aharonian, Phys. Rev. D78 (2008) 034013. [Erratum: Phys. Rev.D82,099901(2010)]

work page 2008
[3]

Cooray, Royal Society Open Science 3 (Mar, 2016) 150555

A. Cooray, Royal Society Open Science 3 (Mar, 2016) 150555

work page 2016
[4]

Petropoulou, S

M. Petropoulou, S. Dimitrakoudis, P. Padovani, A. Mastichiadis, and E. Resconi, Mon. Not. Roy. Astron. Soc. 448 (2015) 2412–2429

work page 2015
[5]

Mon. Not. Roy. Astron. Soc

M. Petropoulou and A. Mastichiadis, "Mon. Not. Roy. Astron. Soc." 447 (Feb, 2015) 36–48

work page 2015
[6]

Murase, F

K. Murase, F. Oikonomou, and M. Petropoulou, Astrophys. J. 865 (2018) 124

work page 2018
[7]

S. Gao, A. Fedynitch, W. Winter, and M. Pohl, Nat. Astron. 3 (2019) 88–92

work page 2019
[8]

IceCube Collaboration, M. G. Aartsen et al., Science 342 (2013) 1242856

2013

Show all 16 references

[9]

IceCube Collaboration, M. G. Aartsen et al., Astrophys. J. 835 (2017) 45

2017
[10]

Huber, These proceedings 916 (2019)

IceCube Collaboration, M. Huber, These proceedings 916 (2019). 6 Searching for neutrino emission from hard X-ray sources with IceCube

2019
[11]

IceCube, Fermi-LA T, MAGIC, AGILE, ASAS-SN, HA WC, H.E.S.S., INTEGRAL, Kanata, Kiso, Kapteyn, Liverpool Telescope, Subaru, Swift NuSTAR, VERITAS, VLA/17B-403 Collaboration, M. G. Aartsen et al., Science 361 (2018) eaat1378

2018
[12]

IceCube Collaboration, M. G. Aartsen et al., Science 361 (2018) 147–151

2018
[13]

Rodrigues, S

X. Rodrigues, S. Gao, A. Fedynitch, A. Palladino, and W. Winter, Astrophys. J. 874 (2019) L29

2019
[14]

S. D. Barthelmy et al., Space Sci. Rev. 120 (2005) 143

2005
[15]

Koss et al., Astrophys

M. Koss et al., Astrophys. J. 850 (2017) 74

2017
[16]

IceCube Collaboration, M. G. Aartsen et al., JINST 12 (2017) P03012. 7

2017

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

S. R. Kelner, F. A. Aharonian, and V . V . Bugayov,Phys. Rev. D74 (2006) 034018. [Erratum: Phys. Rev.D79,039901(2009)]

work page 2006

[2] [2]

S. R. Kelner and F. A. Aharonian, Phys. Rev. D78 (2008) 034013. [Erratum: Phys. Rev.D82,099901(2010)]

work page 2008

[3] [3]

Cooray, Royal Society Open Science 3 (Mar, 2016) 150555

A. Cooray, Royal Society Open Science 3 (Mar, 2016) 150555

work page 2016

[4] [4]

Petropoulou, S

M. Petropoulou, S. Dimitrakoudis, P. Padovani, A. Mastichiadis, and E. Resconi, Mon. Not. Roy. Astron. Soc. 448 (2015) 2412–2429

work page 2015

[5] [5]

Mon. Not. Roy. Astron. Soc

M. Petropoulou and A. Mastichiadis, "Mon. Not. Roy. Astron. Soc." 447 (Feb, 2015) 36–48

work page 2015

[6] [6]

Murase, F

K. Murase, F. Oikonomou, and M. Petropoulou, Astrophys. J. 865 (2018) 124

work page 2018

[7] [7]

S. Gao, A. Fedynitch, W. Winter, and M. Pohl, Nat. Astron. 3 (2019) 88–92

work page 2019

[8] [8]

IceCube Collaboration, M. G. Aartsen et al., Science 342 (2013) 1242856

2013