Pith. sign in

REVIEW 1 major objections 2 minor 36 references

Seven blazars match the hard X-ray to neutrino luminosity relation calibrated on six AGN.

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 →

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

T0 review · grok-4.3

2026-06-30 19:50 UTC pith:BYRDKTAW

load-bearing objection This extends the prior six-source L_hX-L_nu correlation to seven blazars with consistency checks and distance-controlled tests, but all results stay conditional on the IceCube n_s hat values being real signals. the 1 major comments →

arxiv 2605.15360 v2 pith:BYRDKTAW submitted 2026-05-14 astro-ph.HE

Correlation Between Hard X-Ray and Cosmic Neutrino Sources: From Obscured AGN to Blazars

Emma Kun , Imre Bartos , Claudio Ricci , Santiago del Palacio , Francis Halzen , Julia Becker Tjus , Peter L. Biermann , Anna Franckowiak This is my paper
classification astro-ph.HE
keywords AGNblazarshigh-energy neutrinosIceCubehard X-raysluminosity correlationmultimessenger
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper tests whether seven NuSTAR-observed blazars follow the same L_hX to L_ν correlation reported for six active galactic nuclei that show the strongest individual IceCube signals. It first fits the relation on the six sources with a Bayesian regression that accounts for errors on both axes, obtaining a slope near one and an intrinsic scatter of about 0.6 dex. All seven blazars lie within the posterior predictive distribution of this fit under the assumption that the published IceCube n_s values trace real signals. Two diagnostics that remove the shared distance bias still find a residual luminosity association at roughly 3 sigma. The authors frame the outcome as a conditional consistency check rather than a definitive detection.

Core claim

Calibrating the L_hX-L_ν relation on the six published sources via Bayesian regression yields a slope consistent with β = 1 and intrinsic scatter ∼0.6 dex; all seven new blazars are posterior-predictively consistent with this calibration (χ²₇ = 1.58, p = 0.98) under the working hypothesis that the published IceCube n_s values reflect the signal, while distance-controlled diagnostics (τ|z = 0.69 and flux permutation p = 6.3×10^{-4}) indicate a residual L_hX-L_ν association beyond the distance-induced trend.

What carries the argument

The L_hX-L_ν luminosity relation, calibrated by Bayesian regression with errors on both axes and tested via posterior predictive checks plus distance-controlled rank and permutation diagnostics.

Load-bearing premise

The published IceCube n_s values for both the calibration and test sources represent true neutrino signals rather than background fluctuations.

What would settle it

Adding more hard-X-ray-selected blazars to the sample and finding that their neutrino excesses fall systematically outside the calibrated relation, or that the distance-controlled association vanishes.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The correlation between hard X-ray and neutrino luminosity extends from obscured AGN to blazars.
  • Neutrino production occurs in compact, photon-rich environments near supermassive black holes across these source classes.
  • Both populations lie inside the photohadronic prediction band for the L_hX/L_ν ratio.
  • A full detection-level claim requires either a larger calibration sample or an X-ray-weighted IceCube stacking analysis.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • X-ray catalogs could be used to prioritize targets for neutrino telescopes if the relation holds.
  • Similar photohadronic processes may operate in both obscured AGN and blazars despite their different viewing angles.
  • Future multi-messenger surveys can test whether the slope remains unity when more sources are added.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

1 major / 2 minor

Summary. The manuscript calibrates a Bayesian regression (errors on both axes) of the hard X-ray luminosity L_hX versus neutrino luminosity L_ν relation on six previously published AGN sources, obtaining a slope consistent with β=1 and intrinsic scatter ~0.6 dex. It then tests seven NuSTAR-observed blazars for posterior-predictive consistency with this calibration under the explicit working hypothesis that the published IceCube ˆn_s values represent true signals, reporting χ²_7=1.58 (p=0.98). A null-injection test, a distance-free L_hX/L_ν ratio diagnostic, a redshift-partial rank correlation (τ|z=0.69), and a flux-space permutation test on the joint 13-source sample (p=6.3×10^{-4}) are presented; the results are framed as a conditional consistency check rather than a detection.

Significance. If the working hypothesis that the ˆn_s values trace signal holds, the analysis supplies supporting evidence that the previously reported L_hX-L_ν correlation extends from obscured AGN to blazars and lies within the photohadronic prediction band. Explicit strengths include the upfront statement of the conditioning hypothesis, the null-injection test that quantifies the limited power of the consistency statistic at N_cal=6, and the two distance-controlled diagnostics that address the common d_L² bias. These elements provide a transparent, falsifiable consistency test in multimessenger astrophysics.

major comments (1)
  1. [Abstract / results section] Abstract and results: The posterior-predictive consistency (χ²_7=1.58, p=0.98) and the distance-controlled diagnostics (τ|z=0.69, permutation p=6.3×10^{-4}) are all conditional on the published ˆn_s values representing signal rather than background; the null-injection test already shows this consistency statistic has limited power at the calibration sample size of six, so the reported p-value does not independently corroborate the signal hypothesis itself.
minor comments (2)
  1. Notation for L_hX and L_ν should be defined at first use and kept consistent with the prior six-source reference.
  2. [Abstract] The abstract is information-dense; a short sentence clarifying that all quantitative claims remain conditional on the signal hypothesis would improve readability for a broad readership.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and for recognizing the transparent, conditional framing of our analysis. We address the single major comment below.

read point-by-point responses

  1. Referee: [Abstract / results section] Abstract and results: The posterior-predictive consistency (χ²_7=1.58, p=0.98) and the distance-controlled diagnostics (τ|z=0.69, permutation p=6.3×10^{-4}) are all conditional on the published ˆn_s values representing signal rather than background; the null-injection test already shows this consistency statistic has limited power at the calibration sample size of six, so the reported p-value does not independently corroborate the signal hypothesis itself.

    Authors: We agree that every reported statistic is conditional on the working hypothesis that the published ˆn_s values trace signal. The manuscript already states this explicitly in the abstract ('under the working hypothesis that the published IceCube ˆn_s values reflect the signal'), in the results section, and in the final interpretation paragraph, where we describe the work as 'a conditional consistency check' and note that a detection-level claim would require either a larger calibration sample or an X-ray-weighted stacking analysis. The null-injection test is presented precisely to quantify the limited power of the χ²_7 statistic at N_cal=6. Because the conditioning and its implications are already foregrounded, we do not believe additional revision is required. revision: no

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation calibrates the L_hX-L_ν relation via Bayesian regression on an external prior sample of six sources, then applies posterior-predictive consistency checks and distance-controlled diagnostics (τ|z and flux permutation) to seven new independent blazars. The reported χ², p-values, and residual-association statistics are computed directly from the new data under an explicitly stated working hypothesis; no equation or step reduces these outputs to the fitted slope/scatter by construction, and the null-injection test confirms the test's limited power without circular dependence. The central claims rest on the new sample and explicit controls rather than self-referential reduction.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim depends on the fitted intrinsic scatter and slope from the six-source Bayesian regression plus the explicit working hypothesis that IceCube ˆn_s values represent signal. No new physical entities are introduced.

free parameters (2)
  • intrinsic scatter = ~0.6 dex
    Fitted via Bayesian regression on the six AGN sources and reported as ~0.6 dex.
  • slope β = 1
    Fitted in the same regression and stated as consistent with 1.
axioms (1)
  • domain assumption Published IceCube ˆn_s values reflect the signal neutrino counts rather than background
    Explicitly invoked as the working hypothesis required for the posterior predictive test and distance-controlled diagnostics.

pith-pipeline@v0.9.1-grok · 5929 in / 1535 out tokens · 43694 ms · 2026-06-30T19:50:13.040888+00:00 · methodology

0 comments
read the original abstract

The origin of high-energy astrophysical neutrinos remains a key open question in multimessenger astrophysics. A correlation between unabsorbed hard X-ray and high-energy neutrino luminosity has been reported in six active galactic nuclei with the highest individual IceCube significances, linking neutrino production to compact, photon-rich environments near supermassive black holes. We study whether the threshold-near IceCube excesses associated with seven NuSTAR-observed blazars are statistically consistent with that relation. Calibrating the $L_\mathrm{hX}$-$L_\nu$ relation on the six published sources via a Bayesian regression with errors on both axes, the slope is consistent with $\beta = 1$ and the intrinsic scatter is $\sim 0.6$\,dex. All seven new blazars are posterior-predictively consistent with this calibration ($\chi^2_7 = 1.58$, $p = 0.98$) under the working hypothesis that the published IceCube $\hat{n}_s$ values reflect the signal. A null-injection test confirms that, at the present calibration sample size, the consistency test does not by itself adjudicate between signal and selected-background origins. A distance-free $L_\mathrm{hX}/L_\nu$ ratio diagnostic places both populations within the photohadronic prediction band, statistically indistinguishable. Two diagnostics that control the common $d_L^{\,2}$ distance bias, a redshift-partial rank correlation ($\tau|z = 0.69$, $\sim\!2.7\,\sigma$) and a flux-space permutation test on the 13-source joint sample ($p = 6.3\times10^{-4}$, $3.23\,\sigma$), indicate a residual $L_\mathrm{hX}$-$L_\nu$ association beyond the distance-induced trend. We interpret these results as a conditional consistency check; a detection-level statement requires either an enlarged calibration set or an X-ray-weighted IceCube stacking likelihood with internal data.

Figures

Figures reproduced from arXiv: 2605.15360 by Anna Franckowiak, Claudio Ricci, Emma Kun, Francis Halzen, Imre Bartos, Julia Becker Tjus, Peter L. Biermann, Santiago del Palacio.

Figure 1
Figure 1. Figure 1: Posterior-predictive consistency test. The blue line shows the plug-in posterior-predictive median of the K2024 relation calibrated on the six published sources only (shown as red filled circles), with shaded 68% (dark) and 95% (light) bands set by ±σint,med and ±1.96 σint,med around the median line. Open blue squares show the seven new blazars analyzed in this work. The dashed line marks LhX = Lν. All sev… view at source ↗
Figure 2
Figure 2. Figure 2: Null calibration of the posterior-predictive consistency test. Left: Distribution of χ 2 7 under a background-only null in which ˆns values are drawn from Poisson(λbg) truncated to ˆns > 0, for four background expectations λbg ∈ {0.3, 0.5, 1.0, 2.0}. The black dashed curve is the χ 2 7 distribution that would obtain if the test were correctly calibrated against the K2024 relation (median ∼ 7). The observed… view at source ↗
Figure 3
Figure 3. Figure 3: Distance-free log(LhX/Lν) ratio diagnostic. Red filled circles: six K2024 calibration sources; open blue squares: seven new blazars from this work. Vertical dashed lines show inverse-variance-weighted means with shaded 68% bands. The grey band indicates the photohadronic prediction range κ ∈ [0.1, 1]. Population-weighted κ values are 0.67 (calibration) and 0.43 (new sample); the two means differ by −0.19±0… view at source ↗
Figure 4
Figure 4. Figure 4: Flux-space permutation test as a finite-sample diagnostic. Left: Observed joint sample with the weighted forward regression and its ±σ w scat band. Middle: Null distribution from a luminosity-space permutation that breaks the d 2 L coupling in the null but preserves it in the observation, yielding an inflated nominal significance of 5.07 σ. Right: Null distribution from the flux-space permutation adopted h… view at source ↗
Figure 5
Figure 5. Figure 5: Monthly Fermi-LAT 0.1–100 GeV photon flux (top) and NuSTAR 15–55 keV energy flux (bottom) of 3C 454.3 as a function of time. Filled circles show Fermi-LAT detections (TS ≥ 4); downward triangles indicate upper limits (TS < 4). The dashed line marks the quiescent mean flux over 2017–2022. Orange shaded bands and dotted vertical lines indicate the epochs of NuSTAR observations. 14 http://www.astropy.org [PI… view at source ↗
Figure 6
Figure 6. Figure 6: Fermi-LAT and NuSTAR light curves of, from top to bottom: PKS 1441+25 (FSRQ, z = 0.939; the neutrino spectral index was fitted freely); MITG J201534+3710 (FSRQ, z = 0.858); and S3 0458−02 (FSRQ, z = 2.286). Symbols and lines as in [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Fermi-LAT and NuSTAR light curves of, from top to bottom: OJ 287 (BL Lac, z = 0.306); S2 0109+22 (IBL, z = 0.265; the neutrino spectral index was fitted freely); and Ton 599 (FSRQ, z = 0.725). Symbols and lines as in [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

36 extracted references · 35 canonical work pages · 4 internal anchors

  1. [1]

    G., Abbasi, R., Abdou, Y., et al

    Aartsen, M. G., Abbasi, R., Abdou, Y., et al. 2013, PhRvL, 111, 021103, doi: 10.1103/PhysRevLett.111.021103

    work page doi:10.1103/physrevlett.111.021103 2013

  2. [2]

    , keywords =

    Aartsen, M. G., Ackermann, M., Adams, J., et al. 2020, PhRvL, 124, 051103, doi: 10.1103/PhysRevLett.124.051103

    work page doi:10.1103/physrevlett.124.051103 2020

  3. [3]

    Abbasi, et al

    Abbasi, R., Ackermann, M., Adams, J., et al. 2022, Science, 378, 538, doi: 10.1126/science.abg3395

    work page doi:10.1126/science.abg3395 2022

  4. [4]

    2024, arXiv e-prints, arXiv:2406.07601, doi: 10.48550/arXiv.2406.07601

    Abbasi, R., Ackermann, M., Adams, J., et al. 2024, arXiv e-prints, arXiv:2406.07601, doi: 10.48550/arXiv.2406.07601

    work page doi:10.48550/arxiv.2406.07601 2024

  5. [5]

    2023, ApJS, 265, 31, doi: 10.3847/1538-4365/acbb6a

    Abdollahi, S., Ajello, M., Baldini, L., et al. 2023, ApJS, 265, 31, doi: 10.3847/1538-4365/acbb6a

    work page doi:10.3847/1538-4365/acbb6a 2023

  6. [6]

    B., et al

    Ackermann, M., Ajello, M., Atwood, W. B., et al. 2015, ApJ, 810, 14, doi: 10.1088/0004-637X/810/1/14

    work page doi:10.1088/0004-637x/810/1/14 2015

  7. [7]

    L., Ansoldi, S., Antonelli, L

    Ahnen, M. L., Ansoldi, S., Antonelli, L. A., et al. 2015, ApJL, 815, L23, doi: 10.1088/2041-8205/815/2/L23

    work page doi:10.1088/2041-8205/815/2/l23 2015

  8. [8]

    L., Biermann, P

    Allen, M. L., Biermann, P. L., Chieffi, A., et al. 2024, Astroparticle Physics, 161, 102976, doi: 10.1016/j.astropartphys.2024.102976

    work page doi:10.1016/j.astropartphys.2024.102976 2024

  9. [9]

    Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et a...

    work page doi:10.1051/0004-6361/201322068 1996

  10. [10]

    J., Cowie, L

    Barger, A. J., Cowie, L. L., Mushotzky, R. F., & Richards, E. A. 2001, AJ, 121, 662, doi: 10.1086/318742

    work page doi:10.1086/318742 2001

  11. [11]

    C., Blandford, R

    Begelman, M. C., Blandford, R. D., & Rees, M. J. 1984, Reviews of Modern Physics, 56, 255, doi: 10.1103/RevModPhys.56.255

    work page doi:10.1103/revmodphys.56.255 1984

  12. [12]

    2022, ApJ, 939, 43, doi: 10.3847/1538-4357/ac9588

    Eichmann, B., Oikonomou, F., Salvatore, S., Dettmar, R.-J., & Becker Tjus, J. 2022, ApJ, 939, 43, doi: 10.3847/1538-4357/ac9588

    work page doi:10.3847/1538-4357/ac9588 2022

  13. [13]

    1991, ApJL, 380, L51, doi: 10.1086/186171

    Haardt, F., & Maraschi, L. 1991, ApJL, 380, L51, doi: 10.1086/186171

    work page doi:10.1086/186171 1991

  14. [14]

    2022, International Journal of Modern Physics D, 31, 2230003, doi: 10.1142/S0218271822300038

    Halzen, F. 2022, International Journal of Modern Physics D, 31, 2230003, doi: 10.1142/S0218271822300038

    work page doi:10.1142/s0218271822300038 2022

  15. [15]

    Gruppuso, F

    Hewett, P. C., & Wild, V. 2010, MNRAS, 405, 2302, doi: 10.1111/j.1365-2966.2010.16648.x HI4PI Collaboration, Ben Bekhti, N., Fl¨ oer, L., et al. 2016, A&A, 594, A116, doi: 10.1051/0004-6361/201629178

    work page doi:10.1111/j.1365-2966.2010.16648.x 2010

  16. [16]

    Data analysis recipes: Fitting a model to data

    Hogg, D. W., Bovy, J., & Lang, D. 2010, arXiv e-prints, arXiv:1008.4686, doi: 10.48550/arXiv.1008.4686 IceCube Collaboration. 2022, Dataset, doi: 10.21234/03fq-rh11

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1008.4686 2010

  17. [17]

    , keywords =

    Inoue, Y., Khangulyan, D., & Doi, A. 2020, ApJL, 891, L33, doi: 10.3847/2041-8213/ab7661

    work page doi:10.3847/2041-8213/ab7661 2020

  18. [18]

    2022, ApJ, 938, 87, doi: 10.3847/1538-4357/ac8794

    Kawamuro, T., Ricci, C., Imanishi, M., et al. 2022, ApJ, 938, 87, doi: 10.3847/1538-4357/ac8794

    work page doi:10.3847/1538-4357/ac8794 2022

  19. [19]

    I., Kovalev, Y

    Kellermann, K. I., Kovalev, Y. Y., Lister, M. L., et al. 2007, Ap&SS, 311, 231, doi: 10.1007/s10509-007-9622-5

    work page doi:10.1007/s10509-007-9622-5 2007

  20. [20]

    Kelly, B. C. 2007, ApJ, 665, 1489, doi: 10.1086/519947

    work page doi:10.1086/519947 2007

  21. [21]

    J., Ricci, C., Trakhtenbrot, B., et al

    Koss, M. J., Ricci, C., Trakhtenbrot, B., et al. 2022, ApJS, 261, 2, doi: 10.3847/1538-4365/ac6c05

    work page doi:10.3847/1538-4365/ac6c05 2022

  22. [22]

    , keywords =

    Kovalev, Y. Y., Pushkarev, A. B., G´ omez, J. L., et al. 2025, A&A, 700, L12, doi: 10.1051/0004-6361/202555400

    work page doi:10.1051/0004-6361/202555400 2025

  23. [23]

    2023, A&A, 679, A46, doi: 10.1051/0004-6361/202346710

    Kun, E., Bartos, I., Becker Tjus, J., et al. 2023, A&A, 679, A46, doi: 10.1051/0004-6361/202346710

    work page doi:10.1051/0004-6361/202346710 2023

  24. [24]

    , keywords =

    Kun, E., Bartos, I., Becker Tjus, J., et al. 2021, ApJL, 911, L18, doi: 10.3847/2041-8213/abf1ec

    work page doi:10.3847/2041-8213/abf1ec 2021

  25. [25]

    B., et al

    Kun, E., Bartos, I., Tjus, J. B., et al. 2024, PhRvD, 110, 123014, doi: 10.1103/PhysRevD.110.123014

    work page doi:10.1103/physrevd.110.123014 2024

  26. [26]

    On the Apparent Correlation between X-ray and Neutrino Luminosities of Active Galactic Nuclei

    Luo, J.-J., Lu, M.-X., & Liang, Y.-F. 2026, On the Apparent Correlation between X-ray and Neutrino Luminosities of Active Galactic Nuclei, https://arxiv.org/abs/2605.13588

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  27. [27]

    S., French, H

    Miller, J. S., French, H. B., & Hawley, S. A. 1978, in BL Lac Objects, ed. A. M. Wolfe, 176–187

    1978

  28. [28]

    S., & M´ esz´ aros, P

    Murase, K., Kimura, S. S., & M´ esz´ aros, P. 2020, PhRvL, 125, 011101, doi: 10.1103/PhysRevLett.125.011101

    work page doi:10.1103/physrevlett.125.011101 2020

  29. [29]

    Neronov, A., Savchenko, D., & Semikoz, D. V. 2024, PhRvL, 132, 101002, doi: 10.1103/PhysRevLett.132.101002

    work page doi:10.1103/physrevlett.132.101002 2024

  30. [30]

    The 105 month Swift-BAT all-sky hard X-ray survey

    Oh, K., Koss, M., Markwardt, C. B., et al. 2018, ApJS, 235, 4, doi: 10.3847/1538-4365/aaa7fd

    work page internal anchor Pith review doi:10.3847/1538-4365/aaa7fd 2018

  31. [31]

    M., Assef, R

    Padovani, P., Alexander, D. M., Assef, R. J., et al. 2017, A&A Rv, 25, 2, doi: 10.1007/s00159-017-0102-9

    work page doi:10.1007/s00159-017-0102-9 2017

  32. [32]

    2016, MNRAS, 458, 2836, doi: 10.1093/mnras/stw472

    Paiano, S., Landoni, M., Falomo, R., Scarpa, R., & Treves, A. 2016, MNRAS, 458, 2836, doi: 10.1093/mnras/stw472

    work page doi:10.1093/mnras/stw472 2016

  33. [33]

    V., Burenin, R

    Plavin, A. V., Burenin, R. A., Kovalev, Y. Y., et al. 2024, JCAP, 2024, 133, doi: 10.1088/1475-7516/2024/05/133

    work page doi:10.1088/1475-7516/2024/05/133 2024

  34. [34]

    , keywords =

    Ricci, C., Chang, C.-S., Kawamuro, T., et al. 2023, ApJL, 952, L28, doi: 10.3847/2041-8213/acda27

    work page doi:10.3847/2041-8213/acda27 2023

  35. [35]

    2024, MNRAS, 528, 5990, doi: 10.1093/mnras/stae273

    Israyelyan, D. 2024, MNRAS, 528, 5990, doi: 10.1093/mnras/stae273

    work page doi:10.1093/mnras/stae273 2024

  36. [36]

    On the Absorption of X-rays in the Interstellar Medium

    Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914, doi: 10.1086/317016

    work page internal anchor Pith review doi:10.1086/317016 2000