Pith. sign in

REVIEW 4 major objections 6 minor 24 references

Correlation of high energy neutrinos and gamma rays on the direction of Fermi Bubbles

T0 review · 4 major / 6 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 The paper argues that the ten IceCube neutrinos reconstructed inside the Northern Fermi Bubble cannot be connected to the bubble's gamma-ray emission through a single lepto-hadronic model, so they likely have a different origin, while the…

desk verdict An honest, clearly written conference paper, but the Northern Bubble non-association is not established by the background-free event count and hand-picked fits. read the letter →

arxiv 1908.03613 v2 pith:WDDOGKMF submitted 2019-08-09 astro-ph.HE

Paola Alvarez-Hurtado , Nissim Fraija , Antonio Galvan , Antonio Marinelli This is my paper · ORCID
classification astro-ph.HE
keywords FermiBubbleshigh-energyneutrinosIceCubeFermi-LATlepto-hadronicmodelHAWCupperlimitsmulti-messengercorrelationgalacticcenter
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

This paper asks whether the high-energy neutrinos IceCube sees in the direction of the Fermi Bubbles—two giant gamma-ray lobes at the Galactic center—come from the same proton population that makes the bubbles glow in gamma rays. It builds a hybrid lepto-hadronic model and fits it to eight years of Fermi-LAT data, HAWC's upper limits on very-high-energy gamma rays, and ten IceCube neutrino events. For the Northern Bubble, no single choice of model parameters can reproduce both the gamma-ray spectrum and the neutrino flux, leading the authors to conclude those neutrinos probably have a different origin. The Southern Bubble cannot be ruled out because current Cherenkov telescopes do not cover its field of view well.

What carries the argument

The load-bearing object is the hadronic proton spectrum $J_p(E_p) = A_{0\gamma} E_p^{-\alpha_\gamma} \exp[-(E_p/E_{0\gamma})^\beta]$ (Eq. 3.1), built on analytic proton-proton collision spectra, embedded in a hybrid model with a leptonic inverse-Compton component. The method works by fitting this spectrum to gamma-ray data to fix $A_{0\gamma}$, $\alpha_\gamma$, and $E_{0\gamma}$ for two parameter choices per bubble, then comparing the predicted neutrino flux with the IceCube-inferred flux $\phi_\nu E_\nu^2 = n E_\nu/(4\pi A t)$; the inconsistency in the Northern Bubble is the basis for the non-association claim.

What would settle it

Compute a single-zone hadronic model with a different proton spectrum, such as a broken power law or log-parabola, and scan parameters against Fermi-LAT, HAWC upper limits, and IceCube fluxes for the Northern Bubble; if any parameter set fits all three within uncertainties, the non-association conclusion is falsified. Alternatively, a future CTA or SGSO detection of very-high-energy gamma rays from the Southern Bubble with a spectrum matching its neutrino flux would prove association there.

Watch

Extended reading notes

Core claim

The central claim is that the high-energy neutrino events (two tracks and eight showers, above 30 TeV) reconstructed inside the Northern Fermi Bubble's solid angle cannot be connected to the observed Fermi-LAT gamma-ray emission through a single lepto-hadronic model, and therefore those neutrinos likely originate elsewhere. For the Southern Bubble, one of the two model fits can accommodate both the gamma rays and the neutrinos, but without very-high-energy upper limits the association remains unproven rather than excluded.

Load-bearing premise

The conclusion rests on the assumption that the adopted proton spectrum with the two hand-picked parameter sets per bubble covers all plausible hadronic emission models; if some other spectral shape or a multi-component model fit both the gamma-ray and neutrino data simultaneously, the non-association would not follow.

Editorial extensions

If this is right

  • If correct, IceCube's Northern Bubble neutrino events should be treated as background or as tracers of a different source class, not as evidence for hadronic acceleration in the bubble.
  • The Southern Bubble remains a target for next-generation very-high-energy observatories; CTA and SGSO could confirm or reject hadronic emission there.
  • The 1 PeV event IC14 near the Galactic center deserves separate scrutiny; its origin may be the central black hole rather than the bubble.
  • A subdominant hadronic component can coexist with leptonic emission, so future multi-messenger fits should include both components rather than pure hadronic ones.

Reading between the lines

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

  • The paper's two-fit-per-bubble scan is not exhaustive; a fuller parameter scan over $\alpha_\gamma$, $E_{0\gamma}$, and $\beta$ might find a region where a single model fits both gamma rays and neutrinos, which would overturn the 'cannot' conclusion.
  • If the Northern Bubble non-association holds, the IceCube events likely come from unresolved point sources, such as Sgr A* flares or pulsars, within the bubble region; this could be tested by stacking analyses with improved angular resolution.
  • The same hybrid-model test could be applied to other extended Galactic sources with both gamma-ray and neutrino data, offering a systematic way to identify which structures are truly hadronic emitters.
  • The flux estimate in Eq. (3.2) uses a constant effective area; incorporating IceCube's energy-dependent exposure would yield a more robust neutrino flux and could change the comparison.
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.

Referee Report

4 major / 6 minor

Summary. The paper investigates whether the high-energy neutrinos detected by IceCube in the direction of the Fermi Bubbles can be produced by the same hadronic mechanism that, together with a leptonic inverse-Compton component, accounts for the Fermi-LAT gamma-ray spectrum and the HAWC upper limits. Using 10 IceCube events (8 showers and 2 tracks) inside the bubble region, eight years of Fermi-LAT data, HAWC upper limits, and a proton spectrum with an exponential cutoff, the authors fit two parameter sets per bubble. They conclude that the Northern Bubble neutrinos are not associated with the gamma-ray emission through a single hadronic model, while a Southern Bubble association remains possible. The paper is an ICRC2019 proceedings contribution.

Significance. If firmly established, the non-association claim for the Northern Bubble would be a useful multimessenger constraint on the origin of both the Fermi Bubbles and the IceCube events in that sky region. The paper draws on public data from three observatories and attempts a lepto-hadronic description, which is a reasonable exercise for a proceedings paper. The authors are appropriately cautious about the Southern Bubble and point to CTA and SGSO as future decisive tests. However, the central claim currently rests on a very small event sample and on two ad hoc fits per bubble, and the analysis lacks background subtraction, significance testing, and parameter uncertainties; the result is suggestive rather than conclusive.

major comments (4)
  1. [Section 3, Eq. (3.2)] The neutrino flux used in the comparison is obtained by converting all 10 events in the bubble region into a source flux, with n in Eq. (3.2) set to the total event count. No atmospheric-neutrino background expectation for the Fermi Bubble solid angle is computed, and no significance of the spatial excess is reported. Since the HESE shower events have angular resolutions of order 10-15 degrees and the bubbles cover a large solid angle, the expected number of background events over 2101 days is not negligible. If the observed count is compatible with background, the plotted neutrino flux is not a bubble flux, and the failure of fits N1/N2 to describe it does not support the conclusion of non-association. The authors should provide the background prediction and a Poisson or binomial significance for the excess.
  2. [Section 3, Table 2 and Figs. 2-3] The conclusion that the neutrino and gamma-ray emissions are "not associated through a single model" is based on only two manually chosen parameter sets per bubble. No fitting algorithm, goodness-of-fit measure, or uncertainties on A0γ, αγ, and E0γ are given, and the text describes the fits only qualitatively. Because the claim concerns the whole family of spectra of the form in Eq. (3.1), the authors should scan the parameter space or otherwise demonstrate that the two chosen sets bracket the hadronic models compatible with the gamma-ray data; ideally, they would perform a joint fit to gamma-ray and neutrino data and quote a test statistic. As it stands, the analysis only shows that two particular models do not simultaneously match both datasets.
  3. [Section 3, first paragraph] The timescale argument used to justify the hadronic scenario is internally inconsistent. For the 1 µG case the text quotes t_acc = 3.5 yr, t_pp = 3.5×10^9 yr, and t_esc = 3.2×10^9 yr, which violates the stated condition t_acc < t_pp ≲ t_esc because t_pp exceeds t_esc. This matters because the ordering is used to argue that pp interactions can produce neutrinos before escape. The authors should clarify the magnetic field and confinement assumptions and recompute the timescales, or state which of the quoted values is in error.
  4. [Section 3, Eq. (3.2)] The effective area A is not specified, although the sample combines 8 shower events and 2 track events from different catalogs (HESE and EHE) with different angular resolutions and effective areas. A single effective area cannot correctly convert the mixed event sample into a flux. The authors should provide the effective area as a function of energy and topology used for each event, or restrict the analysis to a homogeneous sample. Without this, the neutrino flux points in Figs. 2-3 are not reproducible and the comparison to the models is not quantitative.
minor comments (6)
  1. [Section 2, Table 1] The event list in the text includes "IC15" twice, matching both a shower (HESE) and a track (EHE) entry in Table 1. If these are distinct events, the notation should distinguish them; if they are the same event, one entry is erroneous.
  2. [Section 3, after Table 2] The text says "the cut-off energy varies from a minimum of 1.7 TeV to 3 PeV," which contradicts Table 2, where E0γ is listed in PeV and the minimum is 1.7 PeV.
  3. [Throughout] Typos include "Nothern" for "Northern," "Cerenkov" for "Cherenkov," and "a few Gauss∼ 2µG" in Section 3, which should presumably read "a few µG."
  4. [Section 3] The sentence "The maximum energy that protons can be accelerated is obtained through the timescales" is not supported by any calculation or equation; this statement should either be derived or removed.
  5. [Figures 2-3] The figures are not described in enough detail: the neutrino flux points do not appear to show error bars in the text version, and the legends are not visible in the manuscript text. The authors should add error bars and clarify the plotted curves.
  6. [References] Some references are incomplete, e.g., [15] lists only an arXiv identifier without authors or journal, and [24] is missing a DOI. Please complete the bibliography.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the model-comparison test uses independently fitted gamma-ray spectra and observed neutrino counts.

full rationale

The paper's central step is a model-comparison test rather than a derivation whose output is equivalent to its input. Equation (3.2) converts the 10 IceCube events into an observed neutrino-flux estimate from event counts, which is a data-reduction expression, not a fitted parameter or a prediction. The proton spectra in Table 2 are fitted to Fermi-LAT gamma-ray data (and HAWC upper limits for the Northern Bubble), and the corresponding neutrino fluxes are computed through pp interactions using standard gamma-ray/neutrino production machinery. The neutrino flux estimate is not used to adjust the hadronic model parameters, and no equation or parameter is defined in terms of the conclusion. The non-association claim follows from the failure of the fitted models to describe both gamma-ray and neutrino fluxes simultaneously, which is a legitimate falsification test. The self-citations to earlier Fraija et al. works for timescales and astrophysical inputs are not load-bearing for the central conclusion, and no uniqueness theorem or ansatz is smuggled in via citation. The paper is self-contained with respect to external Fermi-LAT, HAWC, and IceCube data, so no circularity is present.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The paper's central comparison rests on fitted proton spectrum parameters (Table 2), an assumed electron index, and standard hadronic and leptonic emission assumptions. It does not introduce new particles or forces. The fitted parameters are adjusted to the gamma-ray data, and the neutrino prediction depends on them, which weakens the test but is not circular.

free parameters (5)
  • Hadronic fit parameters for Northern fit N1 = E0=2 PeV, alpha=2.15, A0=6.5e55 GeV
    Fitted to Fermi-LAT gamma-ray spectrum for the Northern Bubble (Table 2).
  • Hadronic fit parameters for Northern fit N2 = E0=1.7 PeV, alpha=2, A0=3.87e57 GeV
    Fitted to match the IceCube neutrino flux for the Northern Bubble (Table 2).
  • Hadronic fit parameters for Southern fit S1 = E0=2.5 PeV, alpha=2, A0=8.74e56 GeV
    Fitted to the Southern Bubble gamma-ray and neutrino data (Table 2).
  • Hadronic fit parameters for Southern fit S2 = E0=3 PeV, alpha=2.2, A0=2e58 GeV
    Fitted to the Southern Bubble neutrino flux (Table 2).
  • Leptonic electron spectral index = 2.2
    Assumed in Section 3 as dNe/dE ~ E^-2.2, chosen by hand, not fitted.
assumptions (5)
  • domain assumption Hadronic proton-proton interactions produce gamma rays and neutrinos in a fixed ratio (Kelner et al. 2006).
    Used to convert the fitted proton spectrum into a neutrino flux prediction; cited to ref [23].
  • domain assumption The ordering t_acc < t_pp < t_esc is required for hadronic acceleration to be viable.
    Stated in Section 3; the paper's own values for B=1 uG give t_pp=3.5e9 yr and t_esc=3.2e9 yr, violating t_pp < t_esc.
  • domain assumption The 10 selected IceCube events are astrophysical neutrinos and are spatially inside the Fermi Bubbles.
    The sample combines HESE, EHE, and alert events with different topologies; no atmospheric background subtraction is described.
  • domain assumption The Fermi-LAT bubble gamma-ray emission and bubble geometry are correctly measured and assigned.
    The analysis relies on 8 years of Fermi-LAT data (Pass 7/8) and on bubble boundaries from Fermi-LAT analyses.
  • domain assumption The proton spectrum is a power law with exponential cutoff of the form in eq. 3.1.
    This functional form is standard in the field but constrains the model flexibility; the conclusion depends on it.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Correlation of high energy neutrinos and gamma rays on the direction of Fermi Bubbles." pith.science (2026). https://pith.science/paper/WDDOGKMF

@misc{pith2026190803613,
  author       = {Pith},
  title        = {Pith review of: Correlation of high energy neutrinos and gamma rays on the direction of Fermi Bubbles},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WDDOGKMF}},
  note         = {Machine review of arXiv:1908.03613}
}
read the original abstract

We study the spatial correlation of astrophysical neutrinos detected by IceCube with the geometry of the two large globular structures located in the center of our Galaxy, known as Fermi Bubbles (FB). Using the Fermi-LAT data collected during 8 years and the upper limits derived by the High Altitude Water Cherenkov (HAWC) gamma-ray observatory, we use a hybrid (lepto-hadronic) model to investigate a possible correlation with the high-energy neutrinos in the direction of the Fermi Bubbles. We find that these events are possibly not associated with the Northern Bubble but do not dismiss a possible correlation with the Southern globular structure. We expect in the coming years to improve the gamma-ray observations through the Cerenkov Telescope Array (CTA) and the Southern Gamma-Ray Survey Observatory (SGSO) observatories to test a possible hadronic emission with the Southern Bubble.

Figures

Figures reproduced from arXiv: 1908.03613 by the authors.

Figure 1
Figure 1. a) Neutrinos events collected during 2101 days. b) Location of 10 neutrino events spatially correlated with the Fermi Bubbles [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Lepto-hadronic spectrum for the Northern bubble. electrons are accelerated with protons and are described by a simple power-law dNe/dE ∼ E −2.2 e . We assume that the high-energy photons are generated by the IC scatterings from CMB, IR and SL photons. The dominant supply of low-energy photons is transferred through CMB photons with an energy density of 4 × 10−13 erg cm−3 . Other seed photon field also contributes to… view at source ↗
Figure 3
Figure 3. Lepto-hadronic spectrum for the Southern bubble [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 23 canonical work pages

  1. [1]

    & Weiner, N., Astrophys

    Dobler, G., Finkbeiner, P., Cholis, I., Slatyer, T. & Weiner, N., Astrophys. J. 717, 825–842 (2010)

    work page 2010

  2. [2]

    Su, M., Slatyer, T. R. & Finkbeiner, D. P., Astrophys. J.724, 1044 (2010)

    work page 2010

  3. [3]

    P., Astrophys

    Finkbeiner D. P., Astrophys. J. 614, 186–193 (2004)

    work page 2004

  4. [4]

    Snowden, S., Egger, R., Freyberget, M. et al. Astrophys. J. 485, 125 (1997)

    work page 1997

  5. [5]

    Carretti, E., Crocker, R., Staveley-Smith, L. et al. Nature 493, 66–69 (2013)

    work page 2013

  6. [6]

    & Sun, M, Astrophys

    Mou, G., Yuan, F., Gan, Z. & Sun, M, Astrophys. J. 811, 37 (2015)

    work page 2015

  7. [7]

    Crocker, R. M. & Aharonian, Phys. Rev. Lett. 106, 101102 (2011)

    work page 2011

  8. [8]

    Lacki & Brian, C., MNRAS: Letters 444, L39–L43 (2014)

    work page 2014

Show all 24 references
  1. [9]

    O., Dogiel, V

    Cheng, K., Chernyshov, D. O., Dogiel, V . A., Ko, C. M. & Ip, W.,Astrophys. J. Letts. 731, L17 (2011)

    2011

  2. [10]

    & Ko, C.,Astrophys

    Cheng, K., Chernyshov, D., Dogiel, V . & Ko, C.,Astrophys. J. 790, 23 (2014)

    2014

  3. [11]

    M., Perez, M

    Fraija, N., González, M. M., Perez, M. & Marinelli, A., Astrophys. J. 753, 40 (2012)

    2012

  4. [12]

    Fraija, N., MNRAS 437, 2187–2200 (2014)

    2014

  5. [13]

    & Aguilar-Ruiz, E., Astroparticle Physics 89, 14–22 (2017)

    Fraija, N., Marinelli, A., Galván-Gámez, A. & Aguilar-Ruiz, E., Astroparticle Physics 89, 14–22 (2017)

    2017

  6. [14]

    & Marinelli, A., Astrophys

    Fraija, N. & Marinelli, A., Astrophys. J. 830, 81 (2016)

    2016

  7. [15]

    & Zweibel, E., [ arXiv:1802.05636]

    Yang, K., Ruszkowski, M. & Zweibel, E., [ arXiv:1802.05636]

    arXiv

  8. [16]

    G., Ackermann, M., Adams, J., Aguilar, J

    IceCube Collaboration, Aartsen, M. G., Ackermann, M., Adams, J., Aguilar, J. A., Ahlers, Ahrens, M. et. al. Phys. Rev. Lett. 113, 101101 (2014)

    2014

  9. [17]

    & Razzaque, S., Phys

    Lunardini, C. & Razzaque, S., Phys. Rev. Lett.108, 221102 (2012)

    2012

  10. [18]

    Razzaque, S., Phys. Rev. D 88, 081302 (2013)

    2013

  11. [19]

    & de Diego, J

    Fraija, N., Araya, M., Galvan-Gamez, A. & de Diego, J. A., (2018), [ arXiv:1811.01108]

    2018 arXiv

  12. [20]

    & de Diego, J

    Fraija, N., Aguilar-Ruiz, E., Galván-Gámez, A., Marinelli, A. & de Diego, J. A. , Astrophys. J. 481, 4461–4471 (2018)

    2018

  13. [21]

    Fraija, N., Astrophys. J. 783, 44 (2014)

    2014

  14. [22]

    Shukurov, A., Rodrigues, L. F. , Bushby, P. J., Hollins, J. & Rachen, J. P., A&A623, A113 (2019).1809.03595

    2019 arXiv

  15. [24]

    U., Albert, A., Alfaro, R., Alvarez, C., Alvarez, J

    Abeysekara, A. U., Albert, A., Alfaro, R., Alvarez, C., Alvarez, J. D., Arceo, R. et. al. Astrophys. J. 842, 85 (2017)

    2017

  16. [25]

    Guo, Y ., 32nd ICRC, 271 (2011), [arXiv:1101.5192]. 7

    2011 arXiv

Pith tools

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