{"id":"9ff4c4bc-f078-453b-ae21-1a85d976a9ff","arxiv_id":"2609.26773","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Neutrino-dark matter forward scattering inside AGN black hole spikes alters the flavor composition at Earth, and in the strong-coupling limit the pγ muon-damped prediction falls outside IceCube's 95% flavor contour.","lead":"Neutrinos that climb out of the dense dark matter spike around a supermassive black hole feel a flavor-dependent potential, which can reshuffle which flavor arrives at Earth relative to standard oscillation predictions. If active galactic nuclei dominate the diffuse neutrino flux, these shifts could be large enough to sit outside IceCube's measured flavor contours, making flavor measurements a fresh probe of neutrino-dark matter interactions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"The headline MESE tension depends on pinning neutrino production at the innermost edge (rp ≈ 10 RS ≈ 10^-5 pc) of the cited 10–100 RS source region, where the r^-7/3 spike density is maximal; a factor-of-10 outward shift cuts ρχ by ~215× and raises the required coupling to |εµ|(mχ/GeV)^-1 ~ 2×10^7…","rationale":"Good-faith read: the paper computes complete three-flavor propagation through a r^-7/3 DM spike, checks adiabaticity explicitly (Appendix A), handles antineutrinos with sign-flipped potentials, and is candid that the IceCube comparison is illustrative. The central falsifiable statement is the MESE exclusion of the pγ muon-damped composition. For that statement to hold, three stacked inputs must be true: (i) neutrinos are produced at the spike-density maximum, rp ≈ 10^-5 pc; (ii) the DM background is maximally asymmetric (Δnχ = nχ); (iii) the coupling-weighted density reaches |εµ|(mχ/GeV)^-1 ≥ 10^5. I isolated (i) as the single most load-bearing because the r^-7/3 falloff makes the potential exquisitely sensitive to rp: within the paper's own cited 10–100 RS production region there is a ~215× density swing, so the excluded prediction lives entirely on the innermost edge of the assumed source region, a knife's-edge placement 'for concreteness' rather than derivation. A decisive, cheap test is to re-run the MESE comparison with rp scanned over (and averaged across) that region; if the exclusion collapses for moderate outward shifts, the tension is an input artifact. This is precisely the reader's weakest_assumption, sharpened and quantified, so I agree with that identification and with the CONDITIONAL verdict: the evolution calculation is reported faithfully, but the astrophysical reach rests on unestablished single-point choices, and the asymmetry and coupling-magnitude caveats reinforce the conditionality without compelling a harsher verdict.","tokens_in":15866,"tokens_out":11298,"duration_ms":71874,"concrete_test":"Recompute Fig. 3 for the pγ muon-damped source with |εµ|(mχ/GeV)^-1 = 10^5 fixed, scanning the production radius over the cited 10–100 RS region (rp = 10^-5, 2×10^-5, 5×10^-5, 10^-4 pc), plus an emission-weighted average over that radial range, all other inputs unchanged. If the prediction re-enters the MESE 95% contour for rp ≳ 3×10^-5 pc, or if the weighted-average composition lies inside the contour, the headline tension is an artifact of pinning rp at the spike-density maximum. Secondary check within the same scan: repeat with an asymmetric fraction η = Δnχ/nχ = 0.1 instead of 1; if the excluded region requires εµ∝1/η throughout, the maximal-asymmetry assumption is equally load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The falsifiable claim (pγ muon-damped prediction outside the MESE 95% contour, Figs. 3 and 8) is reached only in the strong-potential regime |εµ|(mχ/GeV)^-1 ≥ 10^5, where the DM potential must vastly exceed Hvac at production. Because the spike density falls as r^-7/3 (Eq. 7), the potential scales as Vµ(rp) ∝ ρχ(rp) ∝ r_p^-7/3. At rp = 10^-5 pc (~10 RS), ρχ ≈ 3×10^18 GeV/cm^3; at 30 RS (3×10^-5 pc) ρχ drops ~5×, and at 100 RS (10^-4 pc) — the outer edge of the very production region the paper itself cites for NGC 1068 (10–100 RS, [45-47]) — ρχ ≈ 1.4×10^16 GeV/cm^3, a ~215× reduction. Within that bracket the coupling needed to retain the MESE-excluding composition rises from ~10^5 to ~2×10^7, and an emission spread over the region yields an energy-dependent intermediate composition (cf. Appendix B, ε=10^2) that the energy-independent MESE comparison does not represent. The choice rp = 10^-5 pc is made 'for concreteness' at the spike-density maximum (r ≈ 9×10^-6 pc, Fig. 1) and is not derived from the production model. Co-load-bearing and equally binary is the maximal-asymmetry assumption Δnχ = nχ (particles only, refs [54,55]): any η = Δnχ/nχ < 1 scales Vµ down proportionally and symmetric DM gives no potential at all, yet no bound on η is imposed. The calculation itself is sound — full three-flavor numerical evolution, explicit adiabaticity diagnostics in Appendix A, honest disclosure that the IceCube comparison is illustrative — so the soft spot is the input selection, not the machinery. The large coupling scale needed (ε ~ 10^5 GF) is also a practical caveat but is disclosed as a limiting-regime illustration.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"The paper proposes that coherent forward scattering of high-energy neutrinos on a dark matter (DM) spike around an AGN's central black hole generates a flavor-dependent matter potential, V_alpha = 2.5e-19 eV (epsilon_alpha Delta n_chi / 1e18 cm^-3) (Eq. 4), which competes with vacuum oscillations at energies around 100 TeV when |epsilon_alpha Delta n_chi| ~ 1e18 cm^-3. Using a Dirac-fermion DM with flavor-dependent vector couplings (Eq. 1) and an adiabatically grown r^-7/3 spike matched to an NFW halo (Eqs. 6-7), the authors solve the full three-flavor radial Schroedinger evolution through the spike (Eq. 9), phase-average the subsequent vacuum propagation (Eq. 10), propagate neutrino and antineutrino components separately (Eq. 11), and compute Earth-frame flavor compositions for p-gamma and p-p pion-decay and muon-damped sources. They compare the results with the IceCube MESE flavor contours under the stated illustrative assumption that AGNs dominate the diffuse flux with common production and spike properties. The central falsifiable claim is that, for a maximally asymmetric (particle-only) DM background and production at r_p = 1e-5 pc, a DM potential acting on the muon flavor and dominating the vacuum terms at production (|epsilon_mu|(m_chi/GeV)^-1 >= 1e5) places the p-gamma muon-damped prediction outside the 95% MESE contour for both signs of epsilon_mu (Fig. 3) and for epsilon_tau < 0 (Fig. 8), suggesting a new probe of neutrino-DM interactions.","tokens_in":16457,"tokens_out":9785,"duration_ms":59832,"significance":"If the central result holds, the proposal is genuinely novel: it uses flavor composition measurements rather than flux attenuation to probe neutrino-DM interactions, and it exploits DM spikes to reach densities where the induced potential competes with vacuum oscillations at ~100 TeV. The paper has real strengths: the three-flavor Hamiltonian and radial-Schroedinger evolution are standard and correctly executed; nonadiabatic transitions are included and diagnosed explicitly through the adiabaticity parameter of Appendix A (Eqs. A2-A3, Fig. 4); the predictions are computed forward at fixed couplings and compared with IceCube data rather than fitted to it, so the claimed exclusion is not circular; and the energy dependence, including SFRD-weighted redshift averaging, is analyzed in Appendix B (Fig. 5). The comparison with IceCube is honestly framed as illustrative, resting on the assumption of AGN domination with common source properties.","major_comments":[{"comment":"The scale of the entire effect is set by the matter potential at the production radius, V_alpha(r_p) proportional to rho_sp(r_p) with rho_sp proportional to r^-7/3 as given in Eq. (7). The paper fixes r_p = 1e-5 pc 'for concreteness' and notes it lies near the spike-density maximum (~9e-6 pc, Fig. 1), yet it simultaneously cites the NGC 1068 production region as 10-100 R_S, i.e., roughly 1e-5 to 1e-4 pc [45-47]. Over that band the spike density falls by about two orders of magnitude (from ~3e18 GeV/cm^3 at 1e-5 pc to ~1e16-5e16 GeV/cm^3 at 1e-4 pc), and because the MESE-excluding composition in Figs. 3 and 8 is only attained at |epsilon_mu|(m_chi/GeV)^-1 >= 1e5, the coupling required to sustain the tension rises by a comparable factor, to ~1e7, at the outer edge of the cited region. The headline exclusion therefore depends on production at the innermost edge of the region quoted from the coronal-models literature, a choice not derived from the production model. The authors should quantify the reach as a function of r_p over the 10-100 R_S band and show that an emission region spread over that band (with the energy-dependent intermediate composition of Appendix B) does not wash out the exclusion.","section":"Section III, Eq. (7) and the 'for concreteness' r_p = 1e-5 pc sentence; Section IV benchmark"},{"comment":"V_alpha is proportional to the net DM number density Delta n_chi = n_chi - n_bar_chi in Eq. (4), and a particle-antiparticle symmetric background produces no potential at all, as the paper correctly notes in Section II. The entire calculation, including the tension claim of Section V, assumes Delta n_chi = n_chi, i.e., a background composed entirely of DM particles (or antiparticles), with no bound imposed on the asymmetry fraction eta = Delta n_chi / n_chi. Since the potential enters linearly in the Hamiltonian of Eq. (3), any eta < 1 scales V_alpha down proportionally, and the flavor shifts in Figs. 2, 3, 7, and 8 shrink smoothly toward the vacuum values as eta approaches zero. This is co-load-bearing with the r_p choice: the paper should state the minimum asymmetry fraction needed at the adopted r_p for the p-gamma muon-damped prediction to leave the 95% MESE contour, and comment on whether asymmetric-DM scenarios of the type cited (refs. 54, 55) can plausibly supply that fraction without violating other constraints.","section":"Section II, Eq. (4); Section III, 'maximally asymmetric background' sentence"},{"comment":"The falsifiable statement of the paper holds only in the strong-potential regime |epsilon_alpha|(m_chi/GeV)^-1 >= 1e5, which is four to five orders of magnitude above the order-unity value at which the DM potential first competes with the vacuum Hamiltonian at production (Eqs. 4-5). The gap is attributable to nonadiabatic transitions destroying the signature at intermediate couplings: Appendix A reports gamma_min ~ 0.1-1 at 1e3-1e4 TeV for |epsilon_mu|(m_chi/GeV)^-1 = 1e2, and Appendix B shows the composition at that coupling moving across the MESE energy range and even turning back for epsilon_mu < 0. As written, the comparison with the energy-independent MESE contour is made only at the upper limit of the coupling range. The paper should present the p-gamma muon-damped flavor composition, and its exclusion of the 95% contour, as a function of |epsilon_alpha|(m_chi/GeV)^-1 over the full range, identifying the coupling window - if any - in which the prediction actually leaves the MESE contour once the energy dependence of the intermediate regime is included.","section":"Section V, Figs. 3 and 8; Appendix B, Fig. 5"}],"minor_comments":[{"comment":"The paper discloses that the IceCube comparison is illustrative and rests on AGN domination with common spike and production properties; given the column density through the inner spike, one additional sentence on the attenuation/regeneration scale for hard scattering (beyond citing ref. 49) would strengthen confidence that neglecting it is quantitatively safe in the adopted density regime.","section":"Section V, MESE comparison"},{"comment":"The envelopes obtained by 'independently varying the oscillation parameters over their marginal 3-sigma ranges' are shown as regions, but the resulting confidence level of those parameter-uncertainty envelopes is not defined; a phrasing such as 'parameter-uncertainty envelope' would be clearer and avoids implying a statistical confidence level.","section":"Captions of Figs. 3, 7, 8"},{"comment":"The Landau entry in ref. [62] appears malformed, with a DOI-like string given in place of a volume and page range; this should be corrected to the standard citation of Landau, Phys. Z. Sowjetunion 2, 46 (1932).","section":"Reference list, ref. [62]"},{"comment":"In the sentence stating that |epsilon_mu|(m_chi/GeV)^-1 >= 1e5 gives compositions 'close to those expected when the DM-induced potential dominates the vacuum terms at production', the order-of-magnitude gap between the naive dominance threshold (order unity) and the adopted 1e5 would benefit from a one-line justification in the main text (the nonadiabaticity argument of Appendix A), so that the strong-coupling choice reads as motivated rather than arbitrary.","section":"Section V, 'Dark matter effects in IceCube flavor measurements'"}],"recommendation":"major_revision","confidential_remarks":"The propagation calculation is sound and honestly presented, but the headline IceCube 'tension' is engineered by two maximally favorable input choices - production at the spike-density maximum (r_p = 1e-5 pc) and a fully particle-asymmetric DM background - compounded by a strong-coupling window (|epsilon|(m_chi/GeV)^-1 >= 1e5) that is four to five orders above the naive threshold. Each of the three major comments asks for a quantitative sensitivity statement rather than a correction, so the paper is fixable within its scope; that is why I recommend major_revision rather than reject. I would also note for the editor that the concurrent works cited as refs. 49 and 58 overlap with the core idea, so the novelty claim should be scoped carefully against those. If the authors convert the 'tension' claim into an exclusion presented as a function of (r_p, eta, |epsilon|), the paper would be considerably more useful and the claim much harder to over-read."},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this is a competent MSW-type calculation applied to a genuinely new environment — the DM spike around an AGN's central black hole — with a full three-flavor numerical treatment that includes nonadiabatic transitions, honest adiabaticity diagnostics, and a forward (not fitted) comparison against the IceCube MESE flavor contours. The headline claim, that the p-gamma muon-damped prediction falls outside the 95% contour for either sign of epsilon_mu, is a real result insofar as the astrophysical inputs are representative.\n\nWhat is new is the environment. Galactic-halo matter effects are prior work (refs 43-44); the r^-7/3 spike boosts the matter potential by ~10^9-10^10 relative to a smooth halo, reaching the interesting regime at ~100 TeV. The machinery is standard and looks correctly executed: the Hamiltonian, radial Schrodinger evolution, phase averaging over the baseline, and the Appendix A adiabaticity parameter all check out. They run full numerics rather than trusting adiabaticity, and they disclose that the IceCube comparison is illustrative. Nothing is fitted — couplings are fixed and predictions compared. The sign-dependence of the muon- and tau-flavor cases is a genuine qualitative feature.\n\nThe soft spots are input selection, not machinery, and the stress-test concern lands. The effect's magnitude pins on the production radius rp = 10^-5 pc, chosen \"for concreteness\" at the spike-density maximum — the inner edge of the 10-100 Rs region the paper itself cites. Shift production to 100 Rs and the r^-7/3 falloff costs ~215x in density, raising the coupling needed to hold the MESE-excluding composition from ~10^5 to ~2x10^7 in (m_chi/GeV)^-1 units; an emission spread over that region gives an energy-dependent intermediate composition (their own Appendix B, epsilon = 10^2) that the energy-independent MESE comparison does not represent. Equally binary is the maximal-asymmetry assumption, particles only: any eta = Delta-n/n < 1 scales the potential down, and symmetric DM yields no potential at all. Neither choice has an uncertainty propagated over it. The large coupling scale epsilon ~ 10^5 is disclosed as limiting-regime illustration, so I do not count it heavily against them. The circularity concern is fair but mild: the effect is enabled by construction, yet the specific outcome — muon-flavor potential putting p-gamma muon-damped outside the 95% contour for both signs — is not something tuned into place.\n\nVerdict matches the reader's: conditionally accept the calculation, hold the conclusion hostage to rp and the asymmetry. This deserves a serious referee; the obvious referee question is exactly the production-radius one, and the authors can either defend it or weaken the claim. I would bring it to reading group and would cite it if I worked in neutrino flavor or neutrino-DM phenomenology.","headline":"Sound and clean three-flavor MSW calculation in a genuinely new environment, but the IceCube tension rides on pinning production at the spike-density maximum and on maximal DM asymmetry.","tokens_in":17058,"tokens_out":2393,"would_cite":true,"duration_ms":16794,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["14.60.Pq","95.35.+d"],"model":"deepseek/deepseek-v4-flash","headline":"Coherent forward scattering on a dark matter spike around an AGN black hole creates a neutrino matter potential that, near 100 TeV, can push predicted flavor ratios outside IceCube's 95% contour.","keywords":["neutrino flavor composition","dark matter spike","coherent forward scattering","active galactic nuclei","MSW-like matter potential","nonstandard neutrino interactions","IceCube MESE analysis","asymmetric dark matter"],"falsifier":"Measure the flavor composition of neutrinos from an individual bright AGN such as NGC 1068 in the 10-1000 TeV band while independently pinning down the production radius from multimessenger data. Under the paper's benchmark, the strong-potential $p\\gamma$ muon-damped composition is approximately (0.02, 0.55, 0.43) for $\\epsilon_\\mu > 0$ or (0.68, 0.08, 0.25) for $\\epsilon_\\mu < 0$, far from the vacuum muon-damped point near (0.18, 0.45, 0.38); a measurement that stays at the vacuum point while the production radius is confirmed near $10^{-5}$ pc would rule the scenario out.","tokens_in":17591,"feed_emoji":"🕳️","tokens_out":15838,"duration_ms":113273,"temperature":0.7,"pith_summary":"The paper sets out to establish that a dark matter spike around an AGN's central black hole can act on neutrinos the way ordinary matter does in the standard MSW (Mikheyev-Smirnov-Wolfenstein) matter effect: coherent forward scattering creates a flavor-dependent potential, $V_\\alpha \\simeq 2.5\\times10^{-19}\\,\\mathrm{eV}\\,(\\epsilon_\\alpha\\,\\Delta n_\\chi/10^{18}\\,\\mathrm{cm}^{-3})$, that alters the flavor composition of the beam as it escapes the source. At neutrino energies near 100 TeV this potential begins to compete with vacuum oscillations once the coupling-weighted net dark matter number density reaches $|\\epsilon_\\alpha\\,\\Delta n_\\chi| \\sim 10^{18}\\,\\mathrm{cm}^{-3}$, a level the adopted NGC 1068 spike can reach at the production radius if the background is maximally asymmetric. Solving the complete three-flavor evolution through the spike, the paper finds Earth flavor compositions for pion-decay and muon-damped sources that depart substantially from vacuum expectations and survive phase averaging over the astronomical baseline. Under the paper's illustrative assumption that AGNs with common production and spike properties dominate the diffuse flux, the strong-potential photohadronic muon-damped prediction falls outside the 95% confidence contour of the IceCube Medium-Energy Starting Events (MESE) flavor measurement, which would make flavor data a new probe of neutrino-dark matter interactions.","feed_headline":"Dark matter spikes push neutrino flavors past IceCube's 95% contour","feed_subtitle":"If AGN spikes are dense and asymmetric, IceCube flavor data near 100 TeV could probe neutrino–dark matter couplings.","key_machinery":"The load-bearing object is the flavor-dependent matter potential $V_\\alpha \\simeq 2.5\\times10^{-19}\\,\\mathrm{eV}\\,(\\epsilon_\\alpha\\,\\Delta n_\\chi/10^{18}\\,\\mathrm{cm}^{-3})$ produced by coherent forward scattering of neutrinos on a dark matter background, the dark-matter analogue of the standard MSW matter potential. Two additional ingredients make it competitive with vacuum oscillations: the adiabatically grown spike profile $\\rho_{\\mathrm{sp}}(r) \\propto r^{-7/3}$ matched to an outer NFW halo, whose density peaks near the adopted production radius $r_p = 10^{-5}$ pc at $\\rho_\\chi(r_p) \\simeq 3\\times10^{18}\\,\\mathrm{GeV}\\,\\mathrm{cm}^{-3}$; and a maximally asymmetric DM background, $\\Delta n_\\chi = \\rho_\\chi/m_\\chi$, since a particle-antiparticle symmetric background gives no potential through the vector current. The calculation solves the full three-flavor radial Schr\\\"odinger equation through this radially varying potential, keeps nonadiabatic transitions, and phase-averages the subsequent vacuum propagation to Earth to obtain the flavor fractions $f_{\\alpha,\\oplus}$.","core_discovery":"The central claim is that coherent forward scattering of high-energy neutrinos on the dark matter spike around an AGN's supermassive black hole generates a flavor-dependent matter potential, $V_\\alpha = 2\\sqrt{2}\\,G_F\\,\\epsilon_\\alpha\\,\\Delta n_\\chi$, where $\\epsilon_\\alpha$ is the flavor-dependent neutrino-DM coupling normalized to the Fermi constant and $\\Delta n_\\chi = n_\\chi - n_{\\bar\\chi}$ is the net particle-antiparticle DM number density. At $E \\sim 100$ TeV this potential competes with the vacuum oscillation terms when $|\\epsilon_\\alpha\\,\\Delta n_\\chi| \\sim 10^{18}\\,\\mathrm{cm}^{-3}$; for a maximally particle- (or antiparticle-) dominated background in the adopted NGC 1068 benchmark that condition is met for $|\\epsilon_\\alpha|\\,(m_\\chi/\\mathrm{GeV})^{-1} \\gtrsim 1$. Evolving the full three-flavor state numerically through the $r^{-7/3}$ spike, including nonadiabatic transitions, and then phase-averaging the vacuum propagation to Earth, the authors find that pion-decay and muon-damped source compositions shift markedly from their vacuum limits, with the direction and size of the shift controlled by the strength, sign, and flavor structure of the interaction. The paper's sharpest quantitative statement is that, when the DM-induced potential acts on the muon flavor and dominates the vacuum terms at production, the photohadronic muon-damped prediction lies outside the 95% C.L. MESE contour for either sign of $\\epsilon_\\mu$, and likewise for $\\epsilon_\\tau < 0$, under the stated assumption that such AGNs dominate the diffuse flux.","pith_inferences":["The two choices the paper fixes for concreteness, production at $r_p = 10^{-5}$ pc near the spike density maximum and a maximally asymmetric DM background, are also the effect's off-switches: if future coronal and photohadronic modeling places AGN neutrino production well outside $10^{-4}$ pc, or if the DM asymmetry is small, the potential drops below the competitive regime and the tension with th","Because the potential is proportional to the net number density $\\Delta n_\\chi$, this probe is inherently a test of asymmetric dark matter: a null flavor measurement would bound the product of the coupling and the asymmetry, not the neutrino-DM coupling itself.","If the effect is real, energy-binned flavor measurements should see the composition swing from the vacuum point toward the strong-potential limit across roughly a decade of energy, and the location of that swing would measure $\\epsilon_\\alpha\\,(m_\\chi/\\mathrm{GeV})^{-1}$.","The note added to the paper records an independent concurrent study of the same mechanism [58], which suggests the spike-induced matter potential is a natural consequence of combining AGN neutrino production with DM spikes rather than an exotic construction."],"forward_implications":["In the strong-potential regime the $p\\gamma$ muon-damped prediction lies outside the 95% MESE contour for both signs of $\\epsilon_\\mu$ and for $\\epsilon_\\tau < 0$, so if such AGNs dominate the diffuse flux the scenario is in tension with current IceCube data.","The matter potential is energy independent while the vacuum terms scale as $E^{-1}$, so in the intermediate coupling regime the flavor composition should vary across the 5 TeV-10 PeV window like an MSW transition, rather than sitting at a fixed point in the flavor triangle.","The sign of $\\epsilon_\\mu$ selects which mass eigenstate the beam is pushed into in the strong-potential limit ($\\nu_3$ for $\\epsilon_\\mu > 0$, $\\nu_1$ for $\\epsilon_\\mu < 0$), so flavor measurements can in principle distinguish a particle-dominated from an antiparticle-dominated DM background.","Interactions acting on different flavors move the predicted composition in different directions, which means flavor data can in principle identify which neutrino flavor couples to the dark matter.","An unspiked Galactic NFW halo would require couplings $10^9$-$10^{10}$ times larger to produce comparable potentials, so the effect is tied to dense spikes around the source black holes rather than to propagation through the Milky Way."],"supporting_citations":[{"why":"Provides the adiabatic-growth dark matter spike profile, $\\rho \\propto r^{-7/3}$, that sets the extreme densities the neutrinos traverse.","marker":"[29]"},{"why":"Provides the adopted NGC 1068 DM-spike benchmark (SMBH mass, NFW matching, spike radius) used in all numerical runs.","marker":"[51]"},{"why":"One of the multimessenger studies motivating neutrino production at 10-100 Schwarzschild radii, which fixes the production radius $r_p$.","marker":"[47]"},{"why":"The IceCube MESE measurement whose 68% and 95% flavor contours the predictions are compared against.","marker":"[10]"},{"why":"Supplies the NuFIT 6.0 oscillation parameters entering the vacuum Hamiltonian and the $3\\sigma$ envelopes of the predictions.","marker":"[50]"},{"why":"Establishes the Wolfenstein matter-potential formalism that the dark matter potential generalizes.","marker":"[27]"},{"why":"Provides the MSW resonance and adiabaticity framework used to diagnose flavor conversion in the spike.","marker":"[28]"},{"why":"Supplies the illustrative IceCube-Gen2 flavor projections used to show future sensitivity.","marker":"[48]"},{"why":"Motivates the maximally asymmetric DM background assumption that makes $\\Delta n_\\chi = \\rho_\\chi/m_\\chi$.","marker":"[54]"},{"why":"Defines the muon-damped source flavor composition $(0:1:0)$ used as a production benchmark.","marker":"[15]"}],"fun_headline_variants":["Neutrino flavor shifts from AGN dark matter spikes","DM spikes alter neutrino flavors, rival IceCube bounds","Dark matter spikes bend neutrino flavors near 100 TeV","AGN dark matter spikes could tip IceCube flavor ratios","Neutrino-DM interactions in spikes test IceCube contours"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The effect's size is set by the matter potential at the production radius, and the paper fixes that radius at $r_p = 10^{-5}$ pc, close to the spike density maximum, while assuming a maximally asymmetric dark matter background (all particles or all antiparticles) so that $\\Delta n_\\chi = \\rho_\\chi/m_\\chi$; if real neutrinos are produced farther out, or the dark matter is nearly particle-antiparticle symmetric, the potential drops and the predicted flavor departures vanish.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino flavor shifts from AGN dark matter spikes","DM spikes alter neutrino flavors, rival IceCube bounds","Dark matter spikes bend neutrino flavors near 100 TeV","AGN dark matter spikes could tip IceCube flavor ratios","Neutrino-DM interactions in spikes test IceCube contours"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":0.00016,"raw_usage":{"cost":0.000159908,"is_byok":false,"cost_details":{"upstream_inference_cost":0.000159908,"upstream_inference_prompt_cost":0.000129108,"upstream_inference_completions_cost":3.08e-05},"total_tokens":1237,"prompt_tokens":1127,"completion_tokens":110,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":256,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":27}},"tokens_in":1127,"tokens_out":110,"cost_source":"openrouter_usage","duration_ms":1612,"image_pages":[],"temperature":1.0,"reasoning_tokens":27,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-09-23T22:56:44.604797+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"falsifier":"Measure the flavor composition of neutrinos from an individual bright AGN such as NGC 1068 in the 10-1000 TeV band while independently pinning down the production radius from multimessenger data. Under the paper's benchmark, the strong-potential $p\\gamma$ muon-damped composition is approximately (0.02, 0.55, 0.43) for $\\epsilon_\\mu > 0$ or (0.68, 0.08, 0.25) for $\\epsilon_\\mu < 0$, far from the vacuum muon-damped point near (0.18, 0.45, 0.38); a measurement that stays at the vacuum point while the production radius is confirmed near $10^{-5}$ pc would rule the scenario out.","supporting_citations":[{"cited_title":"Dark matter annihilation at the galactic center","cited_arxiv_id":"astro-ph/9906391","evidence_quote":"Provides the adiabatic-growth dark matter spike profile, $\\rho \\propto r^{-7/3}$, that sets the extreme densities the neutrinos traverse."},{"cited_title":"Dark matter explanations for the neutrino emission from the Seyfert galaxy NGC 1068","cited_arxiv_id":"2507.16539","evidence_quote":"Provides the adopted NGC 1068 DM-spike benchmark (SMBH mass, NFW matching, spike radius) used in all numerical runs."},{"cited_title":"High-energy neutrinos from the vicinity of the supermassive black hole in NGC 1068","cited_arxiv_id":"2405.20146","evidence_quote":"One of the multimessenger studies motivating neutrino production at 10-100 Schwarzschild radii, which fixes the production radius $r_p$."},{"cited_title":"Coniglione and (on behalf of the KM3NeT collaboration),The km3net neutrino telescope, Journal of Physics: Conference Series632(2015) 012002","cited_arxiv_id":null,"evidence_quote":"The IceCube MESE measurement whose 68% and 95% flavor contours the predictions are compared against."},{"cited_title":"NuFit-6.0: Updated global analysis of three-flavor neutrino oscillations","cited_arxiv_id":"2410.05380","evidence_quote":"Supplies the NuFIT 6.0 oscillation parameters entering the vacuum Hamiltonian and the $3\\sigma$ envelopes of the predictions."},{"cited_title":"Neutrino Oscillations in Matter, Phys. Rev. D17, 2369","cited_arxiv_id":null,"evidence_quote":"Establishes the Wolfenstein matter-potential formalism that the dark matter potential generalizes."},{"cited_title":"Resonance Amplifica- tion of Oscillations in Matter and Spectroscopy of Solar Neutrinos, Sov. J. Nucl. Phys.42, 913","cited_arxiv_id":null,"evidence_quote":"Provides the MSW resonance and adiabaticity framework used to diagnose flavor conversion in the spike."},{"cited_title":"Sensitivity of IceCube-Gen2 to measure flavor composition of Astrophysical neutrinos","cited_arxiv_id":"2308.15220","evidence_quote":"Supplies the illustrative IceCube-Gen2 flavor projections used to show future sensitivity."},{"cited_title":"Review of asymmetric dark matter","cited_arxiv_id":"1305.4939","evidence_quote":"Motivates the maximally asymmetric DM background assumption that makes $\\Delta n_\\chi = \\rho_\\chi/m_\\chi$."},{"cited_title":"Energy dependent neutrino flavor ratios from cosmic accelerators on the Hillas plot","cited_arxiv_id":"1007.0006","evidence_quote":"Defines the muon-damped source flavor composition $(0:1:0)$ used as a production benchmark."}],"review_version":1}