{"id":"701efdc8-3e02-4017-920d-144e1f6f0532","arxiv_id":"2412.08537","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Upper limits on the scalar-dark-matter neutrino coupling are derived from DSNB attenuation projected for DUNE/Hyper-K and from IceCube data on NGC 1068 and TXS 0506+056, assuming DM spikes around AGN black holes.","lead":"Dark matter that scatters off neutrinos would dim the neutrino light arriving from distant supernovae and from two active galaxies. This paper works out how much dimming current and future detectors could see, and uses that to set new limits on the neutrino-dark matter interaction, with the strongest limits relying on dark matter spikes around black holes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"AGN constraints assume neutrino emission at 4 Schwarzschild radii; moving to plausible jet radii relaxes the σ0 bounds by orders of magnitude.","rationale":"The paper's headline result is the AGN constraint (Fig. 3), and the reader correctly identified the DM spike as a weak point. I find an even more load-bearing geometric assumption: the neutrino source is placed at r = 4Rs so that the line of sight traverses the full spike. This is distinct from the spike's existence; even granting the spike, the emission radius for the blazar TXS 0506+056 is plausibly far outside 4Rs, and the column density scales as r_emit^{-4/3}. A factor of 25 in radius changes the bounds by roughly a factor of 70, so the constraints could easily cease to be competitive if the emission region is at tens or hundreds of Schwarzschild radii. The paper does not quantify or caveat this dependence, unlike the spike existence, which it explicitly acknowledges. The DSNB part of the paper is a defensible sensitivity projection, and the cross-section calculation is standard, so the overall verdict remains CONDITIONAL—but the condition should explicitly include the emission-region geometry, not only the spike profile.","tokens_in":15662,"tokens_out":11119,"duration_ms":116959,"concrete_test":"Recompute the right panel of Fig. 3 for both NGC 1068 and TXS 0506+056 with the lower integration limit in Eq. (4.8) (and Eq. (4.14) for TXS) set to r_emit = 10 Rs, 100 Rs, and 1000 Rs, keeping all other inputs fixed. If the σ0 upper bounds shift upward by more than an order of magnitude, the 'more stringent than DSNB-Xenon1T/SuperK' conclusion depends critically on the 4Rs assumption and should be re-presented as a function of r_emit.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 4 computes the optical depth by integrating the DM column from r = 4Rs to Earth (Eq. 4.8), effectively assuming the neutrino source is at the inner edge of the DM spike. This geometric assumption is as consequential as the spike existence itself but is not flagged in the paper's caveats. For TXS 0506+056, a blazar, the neutrino emission region is commonly modeled at 0.1–10 pc from the central engine; for MBH = 3×10^8 M⊙, that corresponds to 10^3–10^5 Rs. Since the spike density scales as r^{-7/3}, the column density scales as r_emit^{-4/3}; moving r_emit from 4Rs to 100Rs reduces τ by a factor ~70, correspondingly relaxing the upper bound on σ0. Even for NGC 1068, if emission occurs at tens of Rs, the constraints weaken substantially. Thus the central claim that AGN bounds are more stringent than DSNB-Xenon1T/SuperK is not robust without quantifying the emission-radius uncertainty.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a neutrino–scalar dark matter interaction mediated by a fermion, computes the νφ scattering cross section in different kinematic limits, and uses the resulting flux attenuation to constrain the coupling. The DSNB part models the diffuse supernova neutrino background with a Fermi–Dirac spectrum and a star-formation-rate parametrization, computes event rates at Super-Kamiokande, Hyper-Kamiokande, and DUNE, and derives upper bounds on y versus mφ using a Poisson likelihood. The AGN part assumes a DM spike at the centers of NGC 1068 and TXS 0506+056, includes saturation by φφ* annihilation, computes the optical depth for neutrinos emitted at the inner edge of the spike, and uses IceCube event rates to set upper bounds on σ0 versus mφ. The paper concludes that the AGN bounds are more stringent than existing DSNB-Xenon1T and DSNB-SuperK limits, while explicitly cautioning that the AGN results rely on the existence of a DM spike.","tokens_in":15875,"tokens_out":4955,"duration_ms":52920,"significance":"If the constraints are robust, the paper provides useful new limits on neutrino–dark matter interactions in the high-energy regime where σ ∝ Eν, and it combines two complementary sources (MeV DSNB and TeV–PeV AGN neutrinos). The authors make several good choices: they use the νFATE package to treat the cascade attenuation, they provide detailed appendices for the cross sections and the νAr scattering calculation, and they explicitly show the effect of φφ* annihilation on the spike saturation density. The main value of the paper is the quantitative comparison of DSNB and AGN probes within one model; however, as discussed in the major comments, the AGN constraints rest on an unquantified geometric assumption about the neutrino emission radius that is not acknowledged in the paper's caveats.","major_comments":[{"comment":"The optical depth is computed by integrating the DM column from r = 4Rs to the observer, which places the neutrino source at the inner edge of the DM spike. For the adopted r^{-7/3} spike profile, the column density from an emission radius r_emit to Rsp scales as r_emit^{-4/3} (neglecting the Rsp term), so increasing r_emit from 4Rs to 100Rs reduces τ by about a factor of 25^{4/3} ≈ 73 and correspondingly relaxes the upper bound on σ0 by the same factor. For TXS 0506+056, blazar emission regions are commonly modeled at 0.1–10 pc from the central engine, which for MBH = 3×10^8 M⊙ corresponds to roughly 10^3–10^5 Rs; even at the lower end of this range the constraints would weaken by several orders of magnitude. The manuscript's caveat at the end of Section 5 mentions only the existence of the DM spike, not the emission-radius assumption. The authors should either justify 4Rs as the appropriate emission radius or quantify the dependence of all AGN constraints on r_emit, since the claim that the AGN bounds are more stringent than DSNB-Xenon1T/SuperK is not robust to this geometric uncertainty.","section":"Section 4, Eqs. (4.8) and (4.14)"},{"comment":"The AGN constraints are derived by requiring Nsct/N ≥ Q with Q = 0.5 (NGC 1068) and Q = 0.05 (TXS 0506+056), where the text says these values include the IceCube uncertainties. The mapping between the quoted experimental uncertainties and these thresholds is not shown. Because the upper bound on σ0 is set approximately by the condition e^{-τ} ≈ Q (or the equivalent attenuation level), a factor-of-2 change in Q shifts the bounds by a factor of about 2 in σ0. The authors should derive Q from the measured fluxes and their errors, or alternatively show the bounds for a range of Q values, so that the quantitative comparison with previous limits is not tied to an unexplained choice.","section":"Section 4, Eq. (4.11)"}],"minor_comments":[{"comment":"The abstract contains the typo 'Kamionkande'; it should read 'Kamiokande'.","section":"Abstract"},{"comment":"The phrase 'through the upper scattering with cosmic electrons or neutrinos' should presumably be 'through upscattering by cosmic electrons or neutrinos'.","section":"Section 1, paragraph 1"},{"comment":"The sentence 'For TXS 0506+056, , since the energy range...' contains a double comma before 'since'.","section":"Section 4, paragraph after Eq. (4.14)"},{"comment":"The phrase 'As For the ϕϕ∗ annihilation' should read 'As for the ϕϕ∗ annihilation'.","section":"Appendix A, after Eq. (A.10)"},{"comment":"'The contribution for r > Ris negligible' contains a typo; it should be 'r > R is negligible'.","section":"Section 4, paragraph after Eq. (4.14)"},{"comment":"The caption begins with 'T able 1' due to a formatting artifact; it should read 'Table 1'.","section":"Table 1, caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of a hep-ph phenomenology journal and I see no citation or novelty concerns. The central derivations are internally consistent and the DSNB part is in good shape. The main issue is the unquantified emission-radius assumption in the AGN analysis; this is fixable by adding a dedicated study of the r_emit dependence or a convincing argument for 4Rs, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The DSNB/DUNE half of this paper is the cleaner part and deserves take-up: they treat the energy-dependent kinematic regions properly, ship a sensible chi-square analysis with SFR uncertainties, and produce new projected limits from DUNE's nu-Ar channel that are absent from the prior literature. The AGN half is more fragile, but not for circularity reasons; the optical-depth and spike-saturation logic is internally consistent. The main problem is that the AGN optical depth assumes the neutrinos are emitted at the inner edge of the spike, r = 4 Rs, with no discussion of how that choice affects the bounds. For TXS 0506+056, a blazar, the emission region is usually modeled at 0.1–10 pc, i.e. 10^3–10^5 Rs. Since the spike column scales as r_emit^{-4/3}, moving from 4 Rs to 100 Rs already relaxes the sigma0 bound by roughly two orders of magnitude. That is as consequential as the spike existence itself, which the paper does caveat. The stress-test note is correct, and the omission matters: the headline claim that AGN beat DSNB-Xenon1T/SuperK is not robust without quantifying this geometric uncertainty. Other soft spots are minor by comparison: the Q thresholds (0.5 for NGC, 0.05 for TXS) are ad hoc, the nuFATE implementation is not public, and the mF = 10 TeV benchmark deserves more motivation. I do not think the central argument collapses; the DSNB constraints stand on their own, and the AGN constraints are a reasonable first pass with a known, but unquantified, astrophysical lever. This is a paper I would send to a serious referee: it is honest work, the calculations are traceable, and the DUNE sensitivity forecast plus the correlated annihilation treatment are worth publishing after the emission-radius question is addressed. I would cite the DSNB/DUNE part, and I would not cite the AGN bounds without a caveat.","headline":"DSNB/DUNE part is a solid, useful sensitivity study; the AGN bounds are real but ride on an unquantified 4 Rs emission-radius assumption the paper never flags.","tokens_in":16404,"tokens_out":1613,"would_cite":true,"duration_ms":20241,"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":"Neutrino flux attenuation yields new upper bounds on neutrino–dark matter coupling, with AGN sources giving the strongest limits.","keywords":["dark matter","neutrino-dark matter interaction","diffuse supernova neutrino background","active galactic nuclei","DM spike","neutrino attenuation","IceCube","DUNE"],"falsifier":"Measure the dark-matter density in the inner several hundred parsecs of NGC 1068 (e.g., through stellar kinematics or gravitational lensing) and check whether the density at r = 4Rs approaches the assumed spike profile; if the density is lower by an order of magnitude or more, the AGN optical depth and the resulting bounds weaken correspondingly, while the DSNB-DUNE prediction would be unaffected.","tokens_in":15455,"feed_emoji":"🔭","tokens_out":3969,"duration_ms":43860,"temperature":0.7,"pith_summary":"The paper argues that a light scalar dark matter particle φ that couples to neutrinos through a heavy fermion mediator F imprints itself by attenuating neutrino fluxes from cosmic sources. Using two very different source classes—the MeV-energy Diffuse Supernova Neutrino Background and TeV-PeV neutrinos from AGNs NGC 1068 and TXS 0506+056—it converts attenuation into upper bounds on the coupling y and on the cross-section parameter σ0. The DSNB analysis, with Super-K data and projected DUNE/Hyper-K exposures, sets limits in the low-mass region; the AGN analysis, incorporating a dark-matter spike around each supermassive black hole, yields constraints that the paper finds more stringent than earlier DSNB-Xenon1T and DSNB-SuperK limits for mφ between 1e-6 and 1 GeV. The AGN bounds are explicitly conditional on the existence of a DM spike at the center of each AGN, which the paper flags as its main caveat.","feed_headline":"AGN neutrinos set strongest bounds on neutrino-dark matter coupling","feed_subtitle":"Using NGC 1068 and TXS 0506+056, limits on σ0 beat DSNB-based constraints across 1e-6 to 1 GeV.","key_machinery":"The central object is the νφ scattering cross section, which has distinct energy dependences in different kinematic regions: σ ∝ Eν² for Eν ≪ mφ ≈ mF, an energy-independent form when mF = mφ ≫ Eν, σ ∝ Eν⁻¹ for Eν ≫ mφ, and σ ∝ Eν in the heavy-mediator limit mF² ≫ Eν mφ ≫ mφ². The argument also relies on the transmittance $e^{{-τ}}$ from the optical depth integral, the DM spike density profile ρ_sp with a saturation density set by φφ* annihilation, and the IceCube event-count comparison through the ratio N_sct/N ≥ Q for each AGN source.","core_discovery":"The central claim is that neutrino–scalar dark matter scattering, mediated by a fermion F, attenuates astrophysical neutrino fluxes in an energy-dependent way, and that measuring this attenuation at existing and upcoming detectors sets meaningful upper limits on the coupling. For the DSNB, the cross section behaves as σνφ ∝ Eν² at low energy or becomes energy independent in the degenerate-mass limit, and the resulting flux suppression is computed including cosmological redshift and then compared with Super-K observations and projected DUNE/Hyper-K event rates. For AGNs, the paper chooses the kinematic region mF² ≫ Eν mφ ≫ mφ², where σνφ ∝ Eν, and computes the optical depth through a DM spike profile around each supermassive black hole, including saturation of the spike density from φφ* annihilation. The resulting constraints on σ0 for NGC 1068 and TXS 0506+056 are stronger than the DSNB-based limits over a broad mass range, with σ0 running from about 7 × 10⁻³⁸ cm² to 3 × 10⁻²⁷ cm² for mφ between 1e-6 and 1 GeV.","pith_inferences":["A more robust version of the AGN bound would marginalize over the spike parameters (MBH, tBH, rh) rather than fixing them to representative values; the paper's central caveat is that the bounds weaken substantially if the spike is absent.","Because the νφ cross section grows linearly with neutrino energy, the same analysis could be extended to map the dark-matter column density as a function of redshift by stacking neutrino sources at different distances.","The coupling of spike saturation to φφ* annihilation suggests that neutrino attenuation and dark-matter self-annihilation could be constrained jointly with future multimessenger observations, separating the two effects through the spectral shape of the attenuation.","If no DM spike forms at the center of AGNs, the DSNB-based limits become the more robust channel; the paper's comparison of AGN and DSNB bounds therefore depends on which astrophysical environment is better understood."],"forward_implications":["DUNE, Hyper-Kamiokande, and Super-K can probe the neutrino–dark matter coupling y in the low-mass region through a detectable suppression of the DSNB flux in the open energy window between 10.8 and 26.4 MeV.","The AGN constraints from NGC 1068 and TXS 0506+056 are more stringent than the earlier DSNB-Xenon1T and DSNB-SuperK bounds for mφ between about 1e-6 and 1 GeV.","For TXS 0506+056, the high neutrino energy makes the Eν-linear cross section very effective, so its constraint remains strong even with a much lower allowed event-fraction Q = 0.05.","Dark-matter self-annihilation saturates the spike density for heavier φ, weakening the AGN bounds and making the σ0 upper bound rise steeply for mφ above about 1e-4 GeV."],"supporting_citations":[{"why":"Supplies the DM spike density profile and the adopted NGC 1068 parameters (MBH, tBH, rh) used for the optical depth calculation.","marker":"[28]"},{"why":"Provides the spike profile form with the saturation-density modification from dark-matter annihilation at AGN centers.","marker":"[27]"},{"why":"Gives IceCube's measured NGC 1068 neutrino flux and best-fit spectral index used to normalize the event-count comparison.","marker":"[23]"},{"why":"Gives IceCube's neutrino observations of TXS 0506+056 that define the second AGN energy range and flux normalization.","marker":"[24]"},{"why":"Defines the transmittance e^{-τ} and the optical-depth integral used to compute neutrino attenuation along the line of sight.","marker":"[32]"},{"why":"Supplies the νFATE package that solves the cascade equation for neutrino flux attenuation with the paper's implemented νφ cross section.","marker":"[34]"},{"why":"Provides the Super-K DSNB observed events and background estimates used in the chi-square analysis for present data.","marker":"[19]"},{"why":"Provides DUNE background event estimates and the liquid-argon exposure assumptions used in the ν-Ar event calculation.","marker":"[40]"}],"fun_headline_variants":["AGN neutrinos tighten dark matter coupling limits","Neutrino attenuation sets new dark matter bounds","DSNB and AGN neutrinos probe dark matter scattering","AGN neutrino flux yields stronger dark matter limits","Neutrino dark matter interaction constrained by AGN flux"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The AGN constraints assume that a dark-matter spike with a specific density profile and fixed parameters (black-hole mass, age, and influence radius) exists at the center of NGC 1068 and TXS 0506+056, and if the spike does not form or has different parameters, the resulting bounds change substantially.","fun_headline_variants_meta":{"raw":{"variants":["AGN neutrinos tighten dark matter coupling limits","Neutrino attenuation sets new dark matter bounds","DSNB and AGN neutrinos probe dark matter scattering","AGN neutrino flux yields stronger dark matter limits","Neutrino dark matter interaction constrained by AGN flux"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000291,"raw_usage":{"total_tokens":1768,"prompt_tokens":1084,"completion_tokens":684,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":700,"completion_tokens_details":{"reasoning_tokens":609}},"tokens_in":700,"tokens_out":684,"duration_ms":7344,"temperature":1.0,"reasoning_tokens":609,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:44:35.911003+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the dark-matter density in the inner several hundred parsecs of NGC 1068 (e.g., through stellar kinematics or gravitational lensing) and check whether the density at r = 4Rs approaches the assumed spike profile; if the density is lower by an order of magnitude or more, the AGN optical depth and the resulting bounds weaken correspondingly, while the DSNB-DUNE prediction would be unaffected.","supporting_citations":[],"review_version":1}