{"id":"92a193d6-d5b6-47fe-ab2a-5dc69d356132","arxiv_id":"2502.03328","paper_version":3,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper derives a range of sterile neutrino masses and transition magnetic moments that would reproduce the LORRI cosmic optical background excess, but the excess is only a 1.5-sigma signal.","lead":"This paper asks whether decaying 'sterile' neutrinos, a proposed dark matter particle, could produce the unexplained optical glow that NASA's New Horizons spacecraft measured. The authors calculate which magnetic strengths would fit the glow, but that glow is only a weak, uncertain signal.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8) yields Iλ per unit wavelength, but Section IV.A sets it equal to the band-integrated LORRI residual 2.99 nW/m2/sr without integrating over 0.4–0.9 μm; this units mismatch shifts the derived dNNγ range.","rationale":"Good-faith read: the paper's goal is to translate a claimed unresolved COB intensity into a constraint on the sterile-to-sterile transition magnetic moment using an EFT one-loop amplitude. The EFT setup and the analytic expressions for the decay width (Eqs. 2–3) are standard; the UV discussion in Section V is speculative but clearly labelled. I checked the redshift integral leading to Eq. (8): using the correct delta function δ(hν(1+z)−Eγ) in Eq. (6) and integrating Eq. (5) reproduces the scaling of Eq. (8), so the derivation of the spectrum itself is not the problem. The load-bearing weakness is the comparison of a differential quantity to an integrated measurement. The LORRI value 2.99 nW/m2/sr is the integral over the instrumental bandpass, not a measurement at one wavelength. Equating Iλ to this number implicitly assumes the band is 1 μm wide and flat and that the DM-decay spectrum is constant across it. Both assumptions fail: the decay spectrum has a finite onset at λe = hc/Δm and falls with λobs through H(z)(1+z)^2, and the LORRI response is nontrivial. The resulting error is larger when Δm is near 2 eV, where most of the LORRI band lies blueward of λe. This directly changes the central output of the paper, not just its interpretation. The reader's weakest assumption about the residual's astronomical origin is also real: at 1.5σ, I299 may vanish under better foreground modelling, and I806 has a different central value. But the band-integration error is independent of whether the residual is real; it makes the quoted dNNγ range unreliable even if the residual is entirely sterile-neutrino decay. The paper also labels its result an upper bound while actually fitting the central value; a genuine upper limit would be larger. I therefore support the reader's REJECT, with the concrete test above as the sharpest single check.","tokens_in":12828,"tokens_out":22831,"duration_ms":200265,"concrete_test":"Recompute the model intensity by integrating Eq. (8) over the LORRI bandpass: for a grid of (m1, Δm), evaluate I_model(m1, Δm, Γ) = ∫_{0.4 μm}^{0.9 μm} Iλ(λ) R(λ) dλ using the LORRI response R(λ) from Ref. [47], then solve for Γ such that I_model equals 2.99 nW/m2/sr (and, for the upper limit, 2.99 + 1.64×2.03 nW/m2/sr). Convert to dNNγ via Eq. (3) and plot against the paper's Fig. 2. If the dNNγ values shift by more than ~50% for any (m1, Δm), the band-integration issue is decisive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (8) is the mean specific intensity per observed wavelength Iλ (nW/m2/sr/μm). The LORRI residual I299 = 2.99 ± 2.03 nW/m2/sr from Ref. [12] is the intensity integrated over the 0.4–0.9 μm band with the LORRI response. Section IV.A, however, says 'we first fix the obtained intensity to I299', i.e., it sets Iλ(λpiv) = I299. This assumes an effective bandwidth of 1 μm and ignores the shape of the bandpass and of the predicted spectrum. The error is not a constant factor: for Δm ≈ 2 eV, λe = hc/Δm ≈ 0.62 μm, so the model predicts zero flux at λpiv = 0.608 μm; the band-integrated intensity is the integral of the declining Iλ over 0.62–0.9 μm and is roughly an order of magnitude smaller than Iλ at the line onset. For larger Δm, λe moves bluer and the bandpass samples higher redshifts, again changing the normalization. Consequently, the required Γ and dNNγ quoted in the abstract and Figs. 2–3 are under-estimated by a model-dependent factor (roughly 2–3 in dNNγ). Additionally, the analysis fits the central value rather than deriving an upper limit; a 95% upper bound would use I299 + 1.64σ. Both issues must be fixed before the claimed bounds can be taken seriously.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies radiative decay of a keV-scale sterile neutrino pair through a sterile-to-sterile transition magnetic moment, N1 -> N2 + gamma, and computes the mean specific intensity of the resulting cosmic optical background. It compares this intensity with the LORRI anomalous COB measurements I806 = 8.06 +/- 1.92 nW/m2/sr and I299 = 2.99 +/- 2.03 nW/m2/sr, solves for the required decay width and transition magnetic moment dNNgamma, and claims upper bounds on dNNgamma for m1 = 1-40 keV and Delta m = 2-10 eV. The paper also compares with X-ray bounds and sketches UV completions that produce small mass splittings.","tokens_in":13231,"tokens_out":10297,"duration_ms":96815,"significance":"The topic is timely: the LORRI COB excess is an interesting anomaly, and an effective-field-theory treatment of sterile-neutrino radiative decay is a reasonable framework. If the claimed bounds were correct, they would provide a new astrophysical constraint on a poorly constrained sterile-to-sterile transition magnetic moment. The paper is clearly organized, uses standard formulas for the decay width and line intensity, and makes a useful comparison to existing X-ray bounds. However, the central quantitative comparison is not performed correctly: a per-unit-wavelength intensity is equated to a band-integrated intensity, and the procedure yields required parameter values rather than upper bounds. The paper does not provide machine-checked derivations, reproducible code, or a detailed bandpass treatment, and the numerical results therefore cannot be accepted as stated.","major_comments":[{"comment":"The quantity computed in Eq. (8) is the mean specific intensity I_lambda in nW/m2/sr/um, but the analysis sets this quantity equal to the LORRI residuals I806 and I299, which are band-integrated intensities in nW/m2/sr over 0.4-0.9 um. Setting I_lambda(lambda_piv) = I299 assumes an effective bandwidth of 1 um and ignores the LORRI bandpass. For Delta m near 2 eV, lambda_e = hc/Delta m is about 0.62 um, so the predicted line begins just redward of the pivot wavelength 0.608 um; the band-integrated flux is the integral of the declining I_lambda over 0.62-0.9 um and is much smaller than I_lambda at the line onset. For larger Delta m the mismatch changes in a model-dependent way. The required Gamma and dNNgamma are therefore underestimated by a factor that must be recomputed with the actual LORRI response; this is the central quantitative comparison of the paper.","section":"Section IV.A and Eq. (8)"},{"comment":"The procedure does not yield an upper bound. The text says 'we first fix the obtained intensity to I806/I299' and solves for Gamma and dNNgamma; the resulting curves are the values required to reproduce the excess, not limits. An upper bound would follow from requiring the predicted intensity to be smaller than the observed excess, or smaller than a chosen upper limit such as I299 + 1.64 sigma, and would be an inequality on dNNgamma. As they stand, the quoted ranges are either required values or, if interpreted as the minimum needed to explain all of the excess, lower bounds on dNNgamma. The abstract and title claim of 'upper bounds' is therefore not supported by the calculation.","section":"Title, Abstract, and Section IV.A"},{"comment":"The uncertainty on the excess is not used. The analysis matches the central values 2.99 +/- 2.03 nW/m2/sr and 8.06 +/- 1.92 nW/m2/sr rather than deriving a limit from the upper end. Since the 2.99 value is only about 1.5 sigma above zero, a 95% upper limit on a BSM contribution would be roughly 6.3 nW/m2/sr, which would shift the dNNgamma bands by about a factor sqrt(6.3/2.99) in dNNgamma before the bandpass correction. The plotted bands appear to reflect only the Delta m variation and not the intensity uncertainty.","section":"Section IV.A and Figs. 2-3"},{"comment":"I806 from Ref. [10] and I299 from Ref. [12] are not independent measurements. Ref. [12] is a reanalysis of the same LORRI data with a revised diffuse Galactic light estimate and supersedes the earlier value. Treating both as separate anomalous intensities and deriving two bands inflates the parameter range and gives the impression of a robust constraint where there is one current excess measurement. The analysis should be based on the current best estimate I299 and should treat the older I806 only as a check of systematic variation, not as an equal-weight constraint.","section":"Section I and Section IV.A"},{"comment":"The calculation sets Omega_chi,0 equal to the total dark matter density, implicitly assuming that the keV sterile neutrinos constitute all of the dark matter. No dark-matter fraction f is introduced. Since the predicted intensity scales linearly with f, every derived dNNgamma value scales as f^{-1/2} when only a fraction f of the dark matter is in the decaying state. The claimed bounds are therefore conditional on f = 1 and are not conservative for f < 1. The paper should state this assumption and propagate f through the constraint.","section":"Section III, Eq. (8)"}],"minor_comments":[{"comment":"The abstract quotes dNNgamma in the range 3e-13 to 1e-9 eV^-1 for the I299 analysis, while the conclusion quotes 3e-12 to 1e-9 eV^-1 for sterile neutrino masses 1-40 keV; the discrepancy should be resolved and the mass range and measurement should be specified in both places.","section":"Abstract and Conclusion"},{"comment":"The Froggatt-Nielsen numerical example is internally inconsistent: with M ~ TeV and epsilon ~ 1e-6, epsilon M is of order MeV, not keV. The parameters should be adjusted or the scaling relations stated more carefully.","section":"Section V, Fig. 4"},{"comment":"There are typographical errors, including 'intesnity' in the introduction and inconsistent uses of 'Fig. (2)' versus 'Fig. 2'; a careful proofread would improve the presentation.","section":"Section I"},{"comment":"The caption refers to 'intensity I = 2.99 +/- 2.03 nW/m2/sr' while the vertical axis of the figure is I_lambda; distinguishing the band-integrated intensity from the specific intensity would prevent the units mismatch from being obscured.","section":"Fig. 1(b) caption"}],"recommendation":"reject","confidential_remarks":"The manuscript has been on arXiv since February 2025 and has been updated, but the central issues remain. I would not recommend major revision unless the authors are willing to redo the intensity comparison using the actual LORRI bandpass, propagate uncertainties, and reframe the result either as a required parameter range or as a genuine upper limit. The use of the superseded 8.06 nW/m2/sr value as an equal-weight constraint is also problematic. The UV-origin section is largely extraneous to the main claim and contains a numerical inconsistency."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plainly: this paper has one genuinely new idea—using the LORRI COB excess to constrain the sterile-to-sterile transition magnetic moment for quasi-degenerate keV sterile neutrinos—but the central claim as written doesn't survive. The authors set the specific intensity Iλ from Eq. (8) equal to the band-integrated LORRI residual 2.99 ± 2.03 nW/m2/sr without integrating over the 0.4–0.9 μm band. For Δm = 2 eV the rest wavelength is 0.62 μm, so the predicted Iλ at the pivot 0.608 μm is exactly zero; the band-integrated intensity is the integral of a declining spectrum and can be an order of magnitude smaller than Iλ at line onset. This shifts the derived dNNγ by a model-dependent factor, roughly 2–3.\n\nSecond, the 'upper bound' language is wrong. The calculation solves for the decay width that reproduces the central value of the excess. That is a required parameter range, not an upper bound. A genuine 95% upper limit would use I299 + 1.64σ ≈ 6.3 nW/m2/sr, which weakens the bound. The excess itself is only 1.5σ above zero, so the whole exercise is contingent on the anomaly being real.\n\nWhat's good: the EFT setup is standard, the decay width formula is correct, and the comparison with X-ray bounds from Natwariya & Nayak is useful. The UV discussion of eV-level splittings via Froggatt–Nielsen or clockwork is clearly speculative but honestly labeled. The paper is transparent about the assumptions, and the algebra from width to intensity is internally consistent.\n\nFor all that, the flaws are load-bearing: the quoted range 3×10^-13 – 10^-9 eV^-1 is not an upper bound, and the units mismatch invalidates the normalization. The paper is salvageable—a proper band integration, propagated errors, and honest reframing as a required range for an optional explanation would fix it. As is, I wouldn't cite the numbers.\n\nThe audience is astroparticle physicists working on decaying dark matter and the COB. They'll find the formalism and the X-ray comparison useful, but they should not quote the derived dNNγ. I'd send it to peer review—a good referee can force the revision—but I'd expect major changes before publication.","headline":"New application of the LORRI COB excess to the sterile-to-sterile transition magnetic moment, but the claimed bound is really a fit to a marginal central value and the intensity comparison mixes units.","tokens_in":13732,"tokens_out":3890,"would_cite":false,"duration_ms":36359,"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":"Sterile neutrinos can explain the cosmic optical excess","keywords":["cosmic optical background","sterile neutrino","transition magnetic moment","dark matter decay","LORRI","keV sterile neutrino","low-energy effective field theory","line intensity mapping"],"falsifier":"Re-analyse the LORRI fields with a revised diffuse-galactic-light or scattered-starlight model: if the residual falls below roughly $1\\,\\mathrm{nW/m^2/sr}$, the inferred $d_{NN\\gamma}$ range no longer has an excess to explain. Alternatively, a dedicated search for the predicted quasi-monochromatic line at $2$–$10\\,\\mathrm{eV}$ in the COB spectrum that finds no line at the required decay width would rule out the sterile-neutrino explanation.","tokens_in":12636,"feed_emoji":"🔭","tokens_out":11450,"duration_ms":87676,"temperature":0.7,"pith_summary":"The paper argues that the unexplained optical glow measured by LORRI on NASA's New Horizons spacecraft—an anomalous intensity of $2.99 \\pm 2.03\\,\\mathrm{nW/m^2/sr}$ after subtracting galaxy counts and foreground light—can be produced by radiatively decaying sterile neutrinos with masses of order keV. If that is right, the same dark-matter candidate that many models invoke would also be the source of a puzzling astrophysical signal, and the measurement would translate directly into bounds on the sterile-to-sterile transition magnetic moment, a poorly constrained particle-physics parameter. Working in a low-energy effective field theory, the authors derive that values of $d_{NN\\gamma}$ between $3\\times10^{-13}\\,\\mathrm{eV^{-1}}$ and $10^{-9}\\,\\mathrm{eV^{-1}}$ are needed to account for the reported excess. The central result is this inferred parameter range, offered as an upper bound on the transition magnetic moment for $\\mathcal{O}(\\mathrm{keV})$-mass sterile neutrinos.","feed_headline":"Sterile neutrinos can explain the cosmic optical excess","feed_subtitle":"A decaying keV-scale dark-matter candidate with a magnetic dipole moment matches LORRI's 2.99 nW/m2/sr excess.","key_machinery":"The central object is the radiative decay of a heavy sterile neutrino through the effective dipole operator $\\mathcal{O}_{NN\\gamma} = (\\bar{N}^c_R i\\sigma_{\\mu\\nu} N_R) F^{\\mu\\nu}$ in the low-energy effective field theory. The decay width is $\\Gamma_{m_1} = (2|d_{NN\\gamma}|^2/\\pi)\\, m_1^3 (2-\\delta)^3 \\delta^3$ with $\\delta = \\Delta m/m_1$, and the mean specific intensity from decaying dark matter is $I_\\lambda = (c/4\\pi)(\\Omega_{\\chi,0}\\rho_c c^2 \\Gamma_{m_1})/(\\lambda_{\\mathrm{obs}}(1+z)H(z))\\,(\\Delta m/m_1)$. Matching this $I_\\lambda$ to LORRI's measured excess, with $\\Delta m$ chosen so the emitted photons fall in the $0.4$–$0.9\\,\\mu\\mathrm{m}$ band after redshift, is what converts an astrophysical brightness into a bound on the coupling $d_{NN\\gamma}$.","core_discovery":"Within a low-energy effective field theory, a keV-scale sterile neutrino $m_1$ can decay radiatively to a slightly lighter sterile neutrino $m_2$ plus a photon through a sterile-to-sterile transition magnetic moment $d_{NN\\gamma}$, with mass splitting $\\Delta m = m_1 - m_2$. The paper shows that if such sterile neutrinos are all of the dark matter, the integrated light from these decays reproduces LORRI's anomalous residual for $m_1$ in the 1–40 keV range, $\\Delta m$ in the 2–10 eV range, and decay widths around $10^{-22}$–$10^{-21}\\,\\mathrm{s^{-1}}$. The corresponding transition magnetic moment is $3\\times10^{-13}\\,\\mathrm{eV^{-1}}$ to $10^{-9}\\,\\mathrm{eV^{-1}}$ for the $2.99\\pm2.03\\,\\mathrm{nW/m^2/sr}$ excess and roughly an order of magnitude larger for the earlier $8.06\\pm1.92\\,\\mathrm{nW/m^2/sr}$ excess. These are presented as upper bounds on $d_{NN\\gamma}$ at keV masses, lying below existing X-ray line limits; sterile-to-active and active-to-active decay channels are found unable to explain the excess.","pith_inferences":["Because the $2.99\\pm2.03\\,\\mathrm{nW/m^2/sr}$ residual is only about $1.5\\sigma$ from zero, a foreground model that shaves off part of the residual would push the required $d_{NN\\gamma}$ upward or eliminate the need for sterile-neutrino decays altogether.","The model predicts a quasi-monochromatic photon line at energy $\\Delta m \\sim 2$–$10\\,\\mathrm{eV}$; high-resolution COB spectroscopy looking for a line rather than a broadband excess would directly test the sterile-neutrino origin.","The same effective-field-theory machinery applies to the cosmic ultraviolet, X-ray, and infrared backgrounds; agreement across bands would strengthen the case, disagreement would point to a different source of the excess.","If competing explanations such as decaying axions or unmodelled intra-halo light account for part of the excess, the quoted numbers should be read as upper limits on $d_{NN\\gamma}$ rather than a required value."],"forward_implications":["The required decay widths ($\\sim 10^{-22}$–$10^{-21}\\,\\mathrm{s^{-1}}$) are far below the inverse age of the universe, so sterile neutrinos remain a viable cold dark matter candidate while producing the excess.","Sterile-to-active and active-to-active decay channels cannot reproduce the excess, singling out the sterile-to-sterile transition magnetic moment within this framework.","The inferred $d_{NN\\gamma}$ values lie below existing X-ray upper bounds, so the scenario is not currently excluded.","Future surveys such as SPHEREx, GALEX, and ULTRASAT can probe the same radiative-decay signal in adjacent wavelength bands and sharpen or overturn the inferred range.","A verified excess would make the COB measurement a new, largely model-independent probe of keV-scale neutrino magnetic moments, complementing collider searches for long-lived particles."],"supporting_citations":[{"why":"Postman et al. 2024 re-analysis; supplies the anomalous rest intensity 2.99 ± 2.03 nW/m2/sr that the paper explains.","marker":"[12]"},{"why":"Lauer et al. 2022; provides the earlier LORRI excess 8.06 ± 1.92 nW/m2/sr used as the second target I806.","marker":"[10]"},{"why":"Creque-Sarbinowski & Kamionkowski 2018; gives the line-intensity mapping formula for the mean specific intensity from decaying dark matter.","marker":"[5]"},{"why":"Atre et al., Bondarenko et al., and Balantekin & Vassh; supply the radiative-decay width expression for heavy sterile neutrinos used in Eq. (2).","marker":"[30-32]"},{"why":"Natwariya & Nayak 2022; provides the X-ray upper bounds on decay width versus mass that delimit the allowed dNNγ region.","marker":"[43]"},{"why":"Beltrán et al. 2024; supplies the low-energy effective field theory operators and the UV completion for small mass splittings.","marker":"[25]"},{"why":"Madau 1995 and Inoue et al. 2014; justify the 2–10 eV photon-energy window via intergalactic absorption at higher redshift.","marker":"[45, 46]"}],"fun_headline_variants":["Sterile neutrino decay matches cosmic background excess","keV sterile neutrinos explain anomalous COB light","Magnetic moment decay fits LORRI's extra glow","Dark matter's sterile neutrinos shine in excess","Sterile neutrino magnetic moment bound from COB"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that the $2.99\\pm2.03\\,\\mathrm{nW/m^2/sr}$ residual in the cosmic optical background is a genuine excess—about $1.5\\sigma$ above zero—and that keV sterile neutrinos comprising all of the dark matter are its only source.","fun_headline_variants_meta":{"raw":{"variants":["Sterile neutrino decay matches cosmic background excess","keV sterile neutrinos explain anomalous COB light","Magnetic moment decay fits LORRI's extra glow","Dark matter's sterile neutrinos shine in excess","Sterile neutrino magnetic moment bound from COB"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000171,"raw_usage":{"total_tokens":1377,"prompt_tokens":1155,"completion_tokens":222,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":771,"completion_tokens_details":{"reasoning_tokens":149}},"tokens_in":771,"tokens_out":222,"duration_ms":3014,"temperature":1.0,"reasoning_tokens":149,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T05:08:26.839857+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-analyse the LORRI fields with a revised diffuse-galactic-light or scattered-starlight model: if the residual falls below roughly $1\\,\\mathrm{nW/m^2/sr}$, the inferred $d_{NN\\gamma}$ range no longer has an excess to explain. Alternatively, a dedicated search for the predicted quasi-monochromatic line at $2$–$10\\,\\mathrm{eV}$ in the COB spectrum that finds no line at the required decay width would rule out the sterile-neutrino explanation.","supporting_citations":[{"cited_title":"Measurement of the Cosmic Optical Background using the Long Range Reconnaissance Imager on New Horizons","cited_arxiv_id":"1704.02989","evidence_quote":"Lauer et al. 2022; provides the earlier LORRI excess 8.06 ± 1.92 nW/m2/sr used as the second target I806."},{"cited_title":"Barducci, E","cited_arxiv_id":null,"evidence_quote":"Beltrán et al. 2024; supplies the low-energy effective field theory operators and the UV completion for small mass splittings."}],"review_version":1}