{"id":"80de4ca4-9496-46c2-b3dc-2802f2602bec","arxiv_id":"1908.03613","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A lepto-hadronic fit to Fermi-LAT and IceCube data finds no single model connects the selected neutrino events to the Northern Fermi Bubble, while the Southern Bubble remains unexcluded.","lead":"This short conference paper asks whether IceCube's highest-energy neutrinos from the Galactic center come from the same source as the Fermi Bubbles' gamma rays. The authors fit a combined lepto-hadronic model to both signals and find the Northern Bubble neutrino events do not match, while the Southern Bubble remains a possible source.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central non-association claim rests on treating 10 IceCube events as a source flux via eq. (3.2) without background subtraction or significance testing; if many are atmospheric, the model comparison cannot support it.","rationale":"The reader identified insufficient model exploration (the hand-picked proton spectra in Table 2) as the weakest assumption. That is a real weakness, but it operates downstream of a more fundamental one: the 10 IceCube events are not demonstrated to be an astrophysical signal from the Fermi Bubbles. The paper provides no atmospheric background estimate and no significance test for the spatial correlation, so the flux inferred by eq. (3.2) may be dominated by background. If that is the case, comparisons in Figures 2 and 3 do not test hadronic association at all. I keep the reader's CONDITIONAL verdict unchanged because the required remedy is well-defined, feasible with public data, and consistent with the paper's own caution about needing further study. The paper is honest and does not overclaim, but the specific limitation that the sample could be background is not addressed in the text.","tokens_in":5420,"tokens_out":8830,"duration_ms":100521,"concrete_test":"Compute the expected number of atmospheric IceCube HESE and EHE events inside the Northern Fermi Bubble mask over 2101 days with the same energy thresholds, topologies, and declination cuts as Table 1, using published IceCube effective areas and standard atmospheric flux models (e.g., Honda et al. or Gaisser et al.). If the background expectation is comparable to or larger than the 10 events (Poisson p approximately greater than 1 percent), eq. (3.2) cannot be used to infer a source flux and the non-association conclusion is unsupported. Conversely, a negligible background expectation would refute this concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's Northern-Bubble conclusion is stated in Section 3 as 'the high-energy neutrinos are not associated with the Fermi Northern Bubble through a single model.' That inference is built on eq. (3.2), which converts the n=10 events in the bubble mask into a neutrino flux phi_nu E^2 = n E_nu / (4 pi A t). This conversion treats every selected event as a signal from the bubbles. No atmospheric-neutrino background expectation is computed, no spatial-excess significance is given, and the effective-area weighting is only sketched. The Fermi Bubbles cover a large solid angle and the HESE showers in Table 1 have median angular resolutions of order 10-15 degrees, so over 2101 days a non-negligible number of atmospheric neutrinos is expected inside that mask. If the 10 events are consistent with background, the 'neutrino flux' plotted against the hadronic models in Figures 2-3 is not a bubble flux, and the failure of fits N1/N2 to also describe it is not evidence against association. A second, related gap is that only two hand-picked parameter sets are tried; eq. (3.1) may not bracket the viable hadronic models. The conclusion may be true, but the evidence presented does not establish it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":5651,"tokens_out":5604,"duration_ms":50469,"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":[{"comment":"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.","section":"Section 3, Eq. (3.2)"},{"comment":"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.","section":"Section 3, Table 2 and Figs. 2-3"},{"comment":"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.","section":"Section 3, first paragraph"},{"comment":"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.","section":"Section 3, Eq. (3.2)"}],"minor_comments":[{"comment":"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.","section":"Section 2, Table 1"},{"comment":"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.","section":"Section 3, after Table 2"},{"comment":"Typos include \"Nothern\" for \"Northern,\" \"Cerenkov\" for \"Cherenkov,\" and \"a few Gauss∼ 2µG\" in Section 3, which should presumably read \"a few µG.\"","section":"Throughout"},{"comment":"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.","section":"Section 3"},{"comment":"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.","section":"Figures 2-3"},{"comment":"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.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings-style contribution with a suggestive but not statistically supported central claim. The main issues (missing background subtraction, lack of a significance test, and the small number of hand-picked fits) are load-bearing for the non-association conclusion but are amenable to revision within the scope of the paper. I recommend major revision rather than rejection, provided the authors can add a background estimate and a parameter scan or otherwise justify that the two fits bracket the viable model space."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThis is a short ICRC conference paper from Alvarez-Hurtado et al. It applies a lepto-hadronic model to ask whether IceCube neutrinos inside the Fermi Bubbles share a common origin with the Fermi-LAT gamma rays. The bottom line: the authors conclude the Northern Bubble is not associated through a single model, and the Southern Bubble remains unproven. If you know the earlier literature, this is not a new result — prior IceCube and HAWC analyses already pointed in that direction. The novelty is the updated event sample (10 events through February 2019) and the explicit two-component modeling.\n\nWhat the paper does well: it is clearly written, uses public data, and is honest about its limits. The hybrid model with a subdominant hadronic component is a reasonable framework, and the authors state plainly that CTA/SGSO data are needed. The discussion of the 1 PeV event IC14 near the Galactic center is a nice touch.\n\nThe soft spots are real and, in one place, load-bearing. The non-association claim is built on equation (3.2), which turns the 10 events in the bubble mask into a neutrino flux with no background subtraction, no significance test, and only a sketched effective area. The sample mixes HESE showers (poor angular resolution) and EHE tracks, and the bubble solid angle is large enough that a few atmospheric neutrinos are expected over 2101 days. If a large fraction of those events are background, the plotted \"neutrino flux\" is not a bubble flux, and the failure of the fits to match it says nothing about association. The two hand-picked fits per bubble (Table 2) are fine for a proceedings plot, but without parameter uncertainties or a scan they cannot exclude a viable hadronic model. There is also a small internal inconsistency in the timescale ordering: for the 1 µG field, they give tpp > tesc, which violates their stated viability condition tacc < tpp ≲ tesc. That is a side issue, since the conclusion rests on the spectral comparison, but it should not have slipped through.\n\nThe paper does not overclaim — the abstract says \"possibly not associated\" — so it deserves credit for caution. But the evidence is too thin to support even that tentative claim. If this were a journal submission, I would not desk-reject it, but I would ask for a proper likelihood analysis with atmospheric background modeling and propagated uncertainties on the spectral parameters.\n\nFor a reading group, it is a good example of how to present a multimessenger search cleanly, and also a good example of how the lack of background modeling weakens a null result. I would probably not cite it in my own work; the earlier papers carry the result.\n\nRecommendation: not strong enough for a full paper, but acceptable as a conference note. If the authors want to make the claim robust, they should add a background estimate and scan the parameter space instead of choosing two fits.\n\nBest.","headline":"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.","tokens_in":6231,"tokens_out":4198,"would_cite":false,"duration_ms":42282,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["Fermi Bubbles","high-energy neutrinos","IceCube","Fermi-LAT","lepto-hadronic model","HAWC upper limits","multi-messenger correlation","galactic center"],"falsifier":"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.","tokens_in":5174,"feed_emoji":"🌌","tokens_out":6679,"duration_ms":64589,"temperature":0.7,"pith_summary":"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.","feed_headline":"Fermi Northern Bubble likely not the source of IceCube neutrinos","feed_subtitle":"A single hadronic model cannot fit gamma-ray and neutrino data there; the Southern Bubble remains an open question.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Identifies the Fermi Bubble gamma-ray emission that defines the target region.","marker":"[1]"},{"why":"Characterizes the Fermi Bubble morphology and spectrum used for the spatial correlation.","marker":"[2]"},{"why":"Supplies the timescale ordering $t_{\\rm acc} < t_{pp} \\lesssim t_{\\rm esc}$ used to justify a hadronic component.","marker":"[7]"},{"why":"Reports the neutrino events near the bubbles and proposes the gamma-neutrino connection being tested.","marker":"[16]"},{"why":"Provides an earlier hadronic-model search for a Fermi Bubble neutrino correlation that this work extends.","marker":"[17]"},{"why":"Supplies the hadronic model framework and the neutrino flux estimate formula used in Eq. (3.2).","marker":"[18]"},{"why":"Provides the analytic proton-proton collision spectra underlying the proton distribution fit.","marker":"[23]"},{"why":"Supplies the very-high-energy upper limits used to constrain the Northern Bubble fit.","marker":"[24]"},{"why":"Gives the energy densities of the seed photon fields used for the inverse-Compton component.","marker":"[25]"}],"fun_headline_variants":["Fermi Northern Bubble unlikely neutrino source","IceCube neutrinos not from Northern Bubble","Southern Bubble still open for neutrino origin","Neutrino correlation favors Southern Fermi Bubble","Hadronic model casts doubt on Northern Bubble neutrinos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Fermi Northern Bubble unlikely neutrino source","IceCube neutrinos not from Northern Bubble","Southern Bubble still open for neutrino origin","Neutrino correlation favors Southern Fermi Bubble","Hadronic model casts doubt on Northern Bubble neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1410,"prompt_tokens":802,"completion_tokens":608,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":418,"completion_tokens_details":{"reasoning_tokens":541}},"tokens_in":418,"tokens_out":608,"duration_ms":6258,"temperature":1.0,"reasoning_tokens":541,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:08:20.358762+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"& Weiner, N., Astrophys","cited_arxiv_id":null,"evidence_quote":"Identifies the Fermi Bubble gamma-ray emission that defines the target region."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Characterizes the Fermi Bubble morphology and spectrum used for the spatial correlation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the timescale ordering $t_{\\rm acc} < t_{pp} \\lesssim t_{\\rm esc}$ used to justify a hadronic component."},{"cited_title":"G., Ackermann, M., Adams, J., Aguilar, J","cited_arxiv_id":null,"evidence_quote":"Reports the neutrino events near the bubbles and proposes the gamma-neutrino connection being tested."},{"cited_title":"& Razzaque, S., Phys","cited_arxiv_id":null,"evidence_quote":"Provides an earlier hadronic-model search for a Fermi Bubble neutrino correlation that this work extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the hadronic model framework and the neutrino flux estimate formula used in Eq. (3.2)."},{"cited_title":"U., Albert, A., Alfaro, R., Alvarez, C., Alvarez, J","cited_arxiv_id":null,"evidence_quote":"Supplies the very-high-energy upper limits used to constrain the Northern Bubble fit."},{"cited_title":"On the Galactic Center Being the Main Source of Galactic Cosmic Rays as Evidenced by Recent Cosmic Ray and Gamma Ray Observations","cited_arxiv_id":"1101.5192","evidence_quote":"Gives the energy densities of the seed photon fields used for the inverse-Compton component."}],"review_version":1}