{"id":"8eeb84a7-4ee4-4e7b-aeff-dcad34b7b3d5","arxiv_id":"2509.07570","paper_version":3,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"Bounds on Z'-mediated dark matter are derived from neutrino telescopes, freeze-in, and cosmic strings, but the telescope bounds rest on comparing a Galactic-center flux to faraway sources.","lead":"The paper argues that high-energy neutrinos seen by IceCube, KM3NeT, ANTARES, Baikal-GVD and PAO could come from dark matter annihilating through a new Z' boson, and it maps the allowed model parameters. It combines the annihilation flux with freeze-in relic density and cosmic-string gravitational waves, but the main telescope bound appears to compare a Galactic-center spike flux to extragalactic sources.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (12) computes a Milky Way Galactic-center spike flux, but Fig. 2 applies it to extragalactic sources (NGC 1068, TXS 0506+056, KM3-230213A) with no distance or J-factor rescaling, so the neutrino-telescope bounds are misplaced.","rationale":"The reader's weakest assumption identifies exactly the load-bearing flaw: Eq. (12) computes the Milky Way Galactic-center spike flux, while the datasets in Fig. 2 are extragalactic point sources and diffuse signals, with no J-factors or distance scalings connecting them. This is an internal inconsistency of the calculation, not merely a disagreement with consensus. The paper's own text says the flux is evaluated for the Galactic halo and opening angle θ around the Galactic center, yet the results section claims bounds from NGC 1068, TXS 0506+056, ANTARES, Baikal-GVD, KM3NeT, and PAO. Since every neutrino-telescope bound in the main figure inherits this mismatch, the central claim is unsupported. The freeze-in and gravitational-wave parts of the paper are independent and may be valid, but they do not rescue the neutrino-telescope constraints, and the abstract/conclusion emphasize the neutrino bounds as a key result. I see no way to repair the argument without recomputing source-specific fluxes, so the rejection stands.","tokens_in":18318,"tokens_out":3852,"duration_ms":46296,"concrete_test":"Take the IceCube NGC 1068 dataset. Recomputed the predicted flux by replacing Eq. (12)'s integrand with an NFW+spike DM profile for NGC 1068 at its actual distance (D ≈ 14 Mpc), using its measured SMBH mass and a l.o.s. integral over that galaxy only; compare this source-specific flux with IceCube's observed 1.5–15 TeV excess for the same (g_X, M_Z'). If it is suppressed by many orders of magnitude relative to the Milky-Way GC flux used in the paper, the ICNGC contour in Fig. 2 disappears. As a simpler cross-check, compute the muon-neutrino event counts predicted by Eq. (12) in IceCube's Galactic-center search region and compare with the event counts from the NGC 1068 point-source analysis; if they disagree, the source attribution is the error.","verdict_should_be":"REJECT","load_bearing_attack":"The central bridge from the particle model to data is Eq. (12), which evaluates dφ_ν/dE_ν by integrating ρ_χ^2(r) along lines of sight through the Milky Way: it uses R_sun = 8.5 kpc, R_vir = 200 kpc, R_sp = 0.7 kpc, M_BH = 10^7 M_sun, θ ∈ [0,10°] around the Galactic center, and the Milky Way NFW+spike profile. This is, by construction, the DM-annihilation neutrino flux from the inner Galaxy as seen from Earth. However, the datasets used for the claimed bounds in Fig. 2 are not Galactic-center measurements: NGC 1068 is a Seyfert galaxy at ~14 Mpc, TXS 0506+056 is a blazar at z ≈ 0.336, and KM3-230213A is an ultra-high-energy event with no identified Milky Way counterpart. To compare with any of these, the calculation must replace the Milky Way halo/spike with the DM distribution of that distant source and include a factor ∝ 1/D^2; even a source at 14 Mpc gives a distance suppression of order (kpc/Mpc)^2 ~ 10^-6 relative to a 1 kpc-scale GC region, before J-factor differences. No such source-specific J-factor or distance rescaling appears anywhere in the text. The results section and Fig. 2 simply label contours 'ICNGC', 'ICTXS', 'KM3', 'PAO', etc., while the flux was computed for the Milky Way GC. Thus the central claim that neutrino telescopes exclude 10^-4 < g_X < 10^-3 at M_Z' ~ O(1 TeV) rests on equating a local GC flux with distant extragalactic fluxes. This is not a minor normalization issue: it invalidates every neutrino-telescope exclusion contour in Fig. 2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a minimal U(1)_X extension of the Standard Model with a Dirac dark-matter particle chi that annihilates through a Z' gauge boson into neutrino pairs. It computes the expected neutrino flux from dark-matter spikes around supermassive black holes using Eq. (12), convolves it with a Gaussian detector resolution, and compares it with observed high-energy neutrino events/excesses (IceCube NGC 1068 and TXS 0506+056, ANTARES, Baikal-GVD, KM3NeT KM3-230213A, and Pierre Auger Observatory). It then derives exclusion contours in the g_X-M_Z' plane, adds freeze-in relic-density constraints, and uses cosmic-string gravitational-wave predictions to constrain the U(1)_X breaking scale. The paper claims that neutrino-telescope bounds can be stronger than collider bounds for M_Z' ~ O(1) TeV, with the strongest bounds coming from freeze-in.","tokens_in":18951,"tokens_out":7748,"duration_ms":86999,"significance":"The model-building part is standard and internally consistent: the U(1)_X framework, the seesaw mechanism, the Z'-mediated annihilation cross section, and the spike density profile are all standard tools, and the paper gives enough benchmark parameters to reproduce Eq. (12). The freeze-in and cosmic-string sections, however, are largely adopted from the authors' earlier work [70,71], and the central astronomical comparison is invalid: Eq. (12) computes a Milky Way Galactic-center l.o.s. integral, while Fig. 2 applies it to distant extragalactic sources with no distance or J-factor rescaling. In addition, no likelihood or test statistic is defined, so the claimed exclusion limits are not statistically grounded. If the source-matching issue were corrected, the framework could be of interest, but the main quantitative claims as presented do not follow from the calculation.","major_comments":[{"comment":"Eq. (12) is a Milky Way Galactic-center l.o.s. integral: it uses R_sun = 8.5 kpc, R_vir = 200 kpc, r^2 = l^2 + R_sun^2 - 2 l R_sun cos(theta), and theta in [0,10 deg] around the Galactic center. The exclusion curves labeled ICNGC, ICTXS, KM3, and PAO in Fig. 2 correspond to extragalactic sources/events: NGC 1068 at ~14 Mpc, TXS 0506+056 at z ~ 0.336, and KM3-230213A with no identified Milky Way counterpart. No source-specific J-factor or 1/D^2 suppression is introduced anywhere. For a source at 14 Mpc the distance suppression relative to the kpc-scale Galactic-center region is ~10^-6, and for TXS it is orders of magnitude larger; even the diffuse ANTARES/Baikal-GVD comparisons lack an angular-acceptance/exposure model. The neutrino-telescope bounds in Fig. 2 are therefore misplaced.","section":"Eq. (12), Fig. 2"},{"comment":"No test statistic or likelihood is defined. The text says that bounds are obtained using the 'best-fit values' of each dataset, but it never specifies the event counts, backgrounds, exposure, energy thresholds, or the statistical procedure that converts a predicted flux into an exclusion. Reading contours directly from the 'best-fit' of a single event or a diffuse flux is not a statistically valid way to derive limits. This affects every claimed neutrino-telescope exclusion contour in Fig. 2.","section":"Results and discussions / Fig. 2"},{"comment":"The freeze-in constraints are advertised as providing 'the strongest bounds,' but the cross sections and the constants C in Eq. (21) are taken from the authors' own previous work [70,71] without derivation, and the asymptotic scalings quoted after Eq. (21) are not shown. No independent benchmark of the freeze-in calculation is given. The freeze-in curves in Fig. 2 are therefore not reproducible from the information in this paper, which is a load-bearing issue for the strongest-bound claim.","section":"Freeze-in, Eq. (21)"}],"minor_comments":[{"comment":"The Schwarzschild radius is written as R_S = 2 G M_BH; the factor c^2 should appear (or the paper should explicitly state units with c=1).","section":"Eq. (7)"},{"comment":"The resolution function contains an extra energy-independent factor exp[-(w ln10)^2/2] in the numerator. For a log-normal detector resolution this constant should not appear; the normalization should be 1/(sqrt(2 pi) w ln10 E') and the exponent should be -(log10(E_nu/E'))^2/(2 w^2).","section":"Eq. (16)"},{"comment":"The figure is overcrowded: many labels overlap, and the gray GW contours and dashed neutrino-telescope lines are difficult to distinguish. Separate panels or larger fonts would greatly improve readability.","section":"Fig. 2"},{"comment":"PAO (Pierre Auger Observatory) is a cosmic-ray observatory, not a neutrino telescope. The paper should clarify how PAO data are used—e.g., through cosmogenic neutrino limits—and avoid describing all datasets as neutrino-telescope measurements.","section":"Introduction / Fig. 2"},{"comment":"The phrase 'n_chi to be even and fractional' is confusing, since the numerical choices n_chi = 100 and 10000 are integers. The charge-quantization argument should be stated more precisely.","section":"Framework section"},{"comment":"Reference [68] (arXiv:2507.16539) appears to be directly relevant to dark-matter explanations of the NGC 1068 neutrino excess, but its results are neither compared nor discussed.","section":"References"}],"recommendation":"reject","confidential_remarks":"The paper's main new claim—that neutrino telescopes exclude 10^-4 < g_X < 10^-3 for M_Z' ~ O(1) TeV—rests on the invalid comparison of a Milky Way Galactic-center flux with extragalactic sources. The lack of a likelihood further prevents any quantitative bound from being assessed. The freeze-in and cosmic-string sections largely follow the authors' previous work, so the novel part is the source-mismatched flux comparison. I recommend rejection, though the underlying particle model could be revisited in a future submission with a correct source treatment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: the neutrino-telescope bounds in Fig. 2 are not what they are labeled. Eq. (12) computes the DM-annihilation neutrino flux from the Milky Way's central spike, with R_sun = 8.5 kpc, R_vir = 200 kpc, and a 10-deg opening angle around the Galactic center. The contours labeled ICNGC, ICTXS, and KM3 are attached to extragalactic sources: NGC 1068 at ~14 Mpc, TXS 0506+056 at z ~ 0.34, and KM3-230213A with no identified host. The paper never constructs J-factors for any of those objects and never includes the 1/D^2 suppression. That is not a minor normalization issue; it is a category error. You cannot read a flux off the Milky Way halo and call it a bound on a model using a Seyfert's neutrinos. The ANTARES and Baikal-GVD entries are diffuse measurements, which is a related but different problem: a GC spike signal would be localized within a few degrees, while the observed diffuse astrophysical flux is dominated by extragalactic sources. No likelihood, test statistic, or background model appears anywhere; the bounds come from 'best-fit values' with no error treatment.\n\nWhat is good: the UV-complete U(1)_X model is standard but handled carefully, the annihilation cross section includes the proper kinematic factors, and the detector-energy smearing is a reasonable touch. The freeze-in section has a clean analytic scaling argument, and the cosmic-string GW spectrum follows the canonical formulas. The big exclusion plot compiles a huge set of terrestrial constraints (LEP, CMS, beam-dump, scattering) that is a useful visual summary even if the new telescope lines are bogus. That said, much of this material is recycled from the group's own refs [65,70,71], so the incremental novelty is modest.\n\nMy overall take: the model part is fine and the paper is readable, but the central astrophysical interpretation is broken. The paper could be repaired by recomputing source-specific fluxes for each object, or by reframing it as a Milky Way GC search using, for example, the ANTARES 11-year GC data. As written, the main result is not supported.\n\nFor peer review: I would send it out rather than desk-reject. A competent referee will catch the source mismatch immediately, and the model-building is solid enough that a major revision could yield a useful, if incremental, paper. It deserves a serious referee, but it should not be accepted anywhere near its current form.","headline":"The model-building is careful but the headline neutrino bounds compare a Milky Way Galactic-center spike flux to extragalactic point sources with no J-factor or distance rescaling, which kills the main claim.","tokens_in":19345,"tokens_out":3502,"would_cite":false,"duration_ms":37417,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proposes that the high-energy neutrinos seen by telescopes could come from dark matter annihilating through a Z' boson, and it derives bounds that in some regions beat collider limits.","keywords":["dark matter annihilation","neutrino telescopes","Z' gauge boson","U(1)_X extension","freeze-in","cosmic strings","Galactic-center spike","neutrino flux"],"falsifier":"Compute the J-factor for NGC 1068 itself, using its distance and a plausible dark-matter halo or spike profile, and evaluate the same flux formula; if the predicted neutrino flux is orders of magnitude below the observed TeV excess, then the Galactic-center comparison the bounds rely on is the wrong target. A null gamma-ray search toward the Galactic center at the level implied by the same Z'-mediated annihilation would also rule the model out.","tokens_in":18220,"feed_emoji":"🔭","tokens_out":13364,"duration_ms":141017,"temperature":0.7,"pith_summary":"The paper tries to establish that the energetic neutrinos observed by current neutrino telescopes—including the excess toward the active galaxy NGC 1068 and the blazar TXS 0506+056—can be produced by dark matter annihilation in the dense spike around the Milky Way's central black hole, with the annihilation mediated entirely by a new Z' gauge boson. If true, the same instruments act as dark-matter detectors: the observed flux translates directly into exclusion bounds on the coupling g_X and mediator mass M_Z', and for M_Z' near 1 TeV those bounds are stronger than collider limits. The paper also shows that producing the dark matter through freeze-in, rather than thermal freeze-out, narrows the allowed coupling to a window around 10^-8 to 10^-4 for mediator masses between 1 GeV and 10^4 GeV, and that cosmic-string gravitational waves can probe the symmetry-breaking scale behind the Z'. A careful reader should notice that the flux calculation is performed for the Galactic-center spike and compared directly with events from distant sources.","feed_headline":"Annihilating dark matter via Z-prime can explain cosmic neutrino excess","feed_subtitle":"If right, neutrino telescopes rival or beat colliders for TeV-scale Z-prime searches.","key_machinery":"The load-bearing object is the differential neutrino flux from dark-matter annihilation in the Galactic-center spike, Eq. (12): d phi/dE = <σv>/(8π m_chi^2) (1/3) dN/dE × the integral over solid angle and line of sight of ρ_chi^2(r), with ρ_chi given by an NFW profile that develops a Gondolo–Silk spike around a 10^7 M_sun black hole. This single formula converts the particle-physics parameters (g_X, M_Z', m_chi, n_chi) into an observable count; the paper evaluates the line-of-sight integral over the Milky Way halo, broadens the monochromatic neutrino line with a Gaussian resolution, and compares the result with the observed fluxes. Freeze-in relic-density contours and cosmic-string gravitati","core_discovery":"The paper argues that the high-energy neutrino events observed toward the active galaxy NGC 1068, the blazar TXS 0506+056, and by ANTARES, Baikal-GVD and KM3NeT can be produced by annihilation of a Dirac dark matter particle into neutrino pairs through a new Z' boson, with the emission coming from the dense dark-matter spike around the Milky Way's central black hole. Comparing the predicted monochromatic neutrino flux, smeared by a Gaussian detector resolution, with observed fluxes yields exclusion curves in the [g_X, M_Z'] plane; for M_Z' near 1 TeV the allowed coupling g_X is forced into the range 10^-4 to 10^-3, and these neutrino-telescope bounds can be stronger than collider limits. Req","pith_inferences":["A natural extension not pursued in the paper is to recompute Eq. (12) with each source's own distance and dark-matter distribution; the resulting couplings could differ by orders of magnitude from the Galactic-center comparison.","The monochromatic line prediction suggests a stacking search for line-like features in neutrino-telescope data could test the model without assuming a particular source; a null result would push the interpretation away from observed event energies.","Because the Z' also couples to charged fermions, a cross-check is to compare the neutrino-derived couplings with gamma-ray upper limits from the Galactic center; this is not done in the paper.","The benchmark m_chi = 3 M_Z' forbids Z' decay into dark matter; opening that channel would change both the freeze-in yield and the annihilation spectrum, so the bounds are specific to this kinematic regime."],"forward_implications":["Neutrino telescopes become direct probes of a TeV-scale Z'-mediated dark sector; for M_Z' near 1 TeV the coupling is pinned between about 10^-4 and 10^-3, a region accessible to current and near-future searches.","Collider limits no longer define the frontier for this model: for a wide range of mediator masses the neutrino-derived bounds are comparable or stronger than dilepton and dijet searches.","Freeze-in relic abundance is more restrictive than the telescope flux itself, narrowing g_X to roughly 10^-8–10^-4 for M_Z' between about 1 GeV and 10^4 GeV.","Cosmic-string gravitational-wave observatories can reach the same U(1) breaking scale for M_Z' ≳ 1 TeV, giving an independent cosmological handle on the model."],"supporting_citations":[{"why":"Supplies the observed neutrino excess from NGC 1068 that the DM-annihilation flux is fit against.","marker":"[20]"},{"why":"Multimessenger observation of the flaring blazar TXS 0506+056; supplies the 290 TeV event used for the blazar bound.","marker":"[21]"},{"why":"Gondolo–Silk spike profile; provides the central density enhancement used in Eq. (7) and the annihilation flux integral.","marker":"[22]"},{"why":"ANTARES diffuse cosmic neutrino flux measurement; one of the datasets used to set bounds.","marker":"[24]"},{"why":"Baikal-GVD diffuse neutrino flux measurement; another dataset used for the bounds.","marker":"[27]"},{"why":"KM3NeT ultra-high-energy event; supplies the high-energy data point used to derive bounds at large M_Z'.","marker":"[28]"},{"why":"NFW halo profiles; provide the outer density distribution matched at the spike radius.","marker":"[61, 62]"},{"why":"Dark-matter spike estimates around the Galactic center; supplies the virial radius and line-of-sight cutoff used in the flux integral.","marker":"[66]"},{"why":"Freeze-in production formalism; the relic-abundance calculation that yields the strongest bounds in the g_X–M_Z' plane.","marker":"[5, 6]"},{"why":"Construction of Z'-mediated neutrino and cosmic-string gravitational-wave constraints; supplies the GW spectrum formulas and comparison basis.","marker":"[71]"}],"fun_headline_variants":["Neutrino telescopes rival colliders for Z-prime dark matter","Dark matter via Z-prime may explain cosmic neutrino excess","Cosmic neutrinos traced to dark matter annihilation","Z-prime dark matter strikes neutrino telescopes","Neutrino sky hints at Z-prime dark matter"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The observed neutrino events from distant sources are treated as coming from a dark-matter spike around the Milky Way's own central black hole, with no per-source distance or halo scaling; if they instead come from the sources' own halos or from ordinary astrophysical emission, the exclusion bounds do not apply.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino telescopes rival colliders for Z-prime dark matter","Dark matter via Z-prime may explain cosmic neutrino excess","Cosmic neutrinos traced to dark matter annihilation","Z-prime dark matter strikes neutrino telescopes","Neutrino sky hints at Z-prime dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000565,"raw_usage":{"total_tokens":2483,"prompt_tokens":679,"completion_tokens":1804,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":423,"completion_tokens_details":{"reasoning_tokens":1726}},"tokens_in":423,"tokens_out":1804,"duration_ms":14380,"temperature":1.0,"reasoning_tokens":1726,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T21:56:54.749870+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the J-factor for NGC 1068 itself, using its distance and a plausible dark-matter halo or spike profile, and evaluate the same flux formula; if the predicted neutrino flux is orders of magnitude below the observed TeV excess, then the Galactic-center comparison the bounds rely on is the wrong target. A null gamma-ray search toward the Galactic center at the level implied by the same Z'-mediated annihilation would also rule the model out.","supporting_citations":[{"cited_title":"Hunting for heavy $Z^\\prime$ with IceCube neutrinos and gravitational waves","cited_arxiv_id":"2502.13217","evidence_quote":"Construction of Z'-mediated neutrino and cosmic-string gravitational-wave constraints; supplies the GW spectrum formulas and comparison basis."}],"review_version":1}