{"id":"ff256505-4f13-4aa9-bc16-2540e9718068","arxiv_id":"2505.22463","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"A new production channel, nu_a + nu_a -> nu_s + nu_s, mediated by a heavy scalar, can generate the observed sterile neutrino dark matter abundance independently of active-sterile mixing.","lead":"This paper proposes a scalar-mediated interaction between ordinary and sterile neutrinos that can produce sterile neutrino dark matter through a new two-to-two reaction, working even at zero mixing angle. It reopens parameter space for sterile neutrino dark matter and shows how future X-ray and neutrino telescopes could test it.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The forward number-changing rate for the benchmark exceeds the Hubble rate at T_i, so the claimed freeze-in behavior and α~1e-3 appear internally inconsistent.","rationale":"I read the paper as a proposal of a novel active-sterile NSI production mechanism, with the central quantitative claim being the freeze-in contours in Figs. 1 and 4 and the benchmark evolution in Fig. 3. The reader's weakest assumption was the two-parameter spectral closure. While that is a legitimate precision concern, the more immediate and severe issue is that the paper's own rate formula, Eq. (10), appears to put the benchmark in the chemical-equilibration regime rather than the freeze-in regime: the forward number-changing rate exceeds H at the initial temperature by a factor of several. If this is correct, the sterile population should thermalize (α→1) and overproduce DM, contradicting the α∼10^-3 shown in Fig. S2 and the entire freeze-in interpretation. The Fig. S1 caption asserts the opposite, so there is an internal inconsistency that must be resolved before the parameter-space claims can be accepted. The proposed concrete test directly evaluates the rate and, if Γ/H>1, would invalidate the central benchmark and the associated contours. I therefore disagree with the reader's assessment that the spectral closure is the weakest link, and recommend rejection unless the rate discrepancy is resolved.","tokens_in":15958,"tokens_out":58252,"duration_ms":604760,"concrete_test":"Compute Eq. (S10) with the full cross section (S16)/(S15) for the benchmark point at T_i=1.5 GeV and T=1 GeV, evaluating Γ_{νaνa→νsνs}(p) for p=T and p=3T, and compare to H. If Γ/H>1, rerun the Boltzmann solution from T_i=1.5 GeV without the α approximation and check whether α approaches 1; if the final Ω_s h^2 changes by more than an order of magnitude, the paper's central parameter region is invalid.","verdict_should_be":"REJECT","load_bearing_attack":"For the benchmark point (m_s=10 keV, sin^2 2θ=10^-12, y_as=3e-4, m_phi=25 GeV), Eq. (10)/(S17) gives Γ_{νaνa→νsνs}(p) = (7π y_as^4/(216 m_phi^4)) p T^4. At the integration start T_i=1.5 GeV with p=T, this gives Γ≈1.6×10^-20 GeV, whereas H≈2.6×10^-21 GeV (g*≈70), so Γ/H≈6 (and Γ/H≈19 for p=3T). Thus the process ν_aν_a→ν_sν_s is faster than Hubble for T≳0.4 GeV. If this rate is correct, the sterile sector should approach chemical equilibrium (α→1) during the early integration, contradicting the freeze-in assumption and the α∼10^-3 shown in Fig. S2 for the same parameters; it would instead overproduce DM. The paper's Fig. S1 caption states that the ν_aν_a↔ν_sν_s rates stay below H over the whole range, which appears inconsistent with Eq. (10). If the rate formula is missing a suppression (e.g., it should contain α T_s^4 for both directions or an additional power of T/m_phi), then the claimed contours in Figs. 1 and 4 and the yield in Fig. 3 are not supported by the stated equations. This is more fundamental than the spectral-closure issue: it questions whether the benchmark is in the freeze-in regime at all.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new production mechanism for sterile-neutrino dark matter, in which a scalar mediator couples active to sterile neutrinos (Eq. 1) and opens the number-changing channel ν_aν_a → ν_sν_s (Fig. 2c). The authors argue that this channel produces sterile neutrinos at a rate independent of the active-sterile mixing angle, so the observed relic abundance Ω_s h² = 0.12 can be obtained even for sin²2θ → 0, evading the Dodelson-Widrow exclusion. They solve integrated Boltzmann equations for the sterile number and energy densities, closing the system with a two-parameter modified Fermi-Dirac ansatz (Eq. 2), and present allowed regions in the (m_s, sin²2θ) and (m_φ, y_as) planes. The paper also discusses X-ray, Lyman-α, and future laboratory sensitivities.","tokens_in":16370,"tokens_out":32874,"duration_ms":337716,"significance":"If the quantitative results are correct, the mechanism is genuinely novel: it opens a previously unexplored production channel for sterile-neutrino DM that does not rely on mixing, a resonance, or a lepton asymmetry. The paper contains a serious effort to confront the parameter space with astrophysical constraints and to identify falsifiable signatures (X-ray and gamma-ray telescopes, beta-decay experiments, DSNB and supernova probes). The analytic collision terms in the heavy-mediator limit are provided in the supplemental material, and the DW limit is recovered as a check. These are strengths. However, the quantitative relic-density contours rest on two assumptions that are not adequately validated: the spectral closure of Eq. (2) and the choice of initial temperature. Both need to be addressed before the central claim can be considered established.","major_comments":[{"comment":"The claim that the final result is insensitive to the initial temperature T_i is not supported. For the benchmark (m_s=10 keV, y_as=3e-4, m_φ=25 GeV), Eq. (10) with p≃T gives Γ_{νaνa→νsνs}/H ≈ 6×10^-3 at T=1.5 GeV, using H≈2.6×10^-18 GeV for g_*=70. This ratio grows as T^3 and reaches unity at T≈6–8 GeV. A standard radiation-dominated universe passes through temperatures above this before cooling to T_i=1.5 GeV, so the zero-abundance initial condition is not the generic one: the sterile sector would have reached chemical equilibrium at higher temperatures, invalidating the freeze-in calculation and leading to a thermal (overproduced) abundance. The paper must either impose an explicit upper bound on the reheating temperature (T_RH ≲ a few GeV for the benchmark) and restrict the parameter space accordingly, or integrate from a sufficiently high temperature and demonstrate that the abundance is not overproduced. As written, the statement of insensitivity to T_i is incorrect and the parameter-space plots in Figs. 1 and 4 implicitly assume a non-standard low reheating temperature that is not stated.","section":"Footnote 2 and Eq. (10)"},{"comment":"The two-moment closure assumes that f_s(p) is always of the form α/(e^{p/T_s}+1), with α and T_s fixed by n_s and ρ_s. All collision terms, including the inverse number-changing rate and the modified DW production term, are evaluated using this ansatz. For the benchmark, the elastic scattering rate ν_aν_s→ν_aν_s is smaller than H over most of the integration range (it becomes comparable to H only near pT∼m_φ², i.e., T≳10 GeV, which is above the chosen T_i), so kinetic equilibrium does not justify the ansatz. The resulting relic-density contours in Figs. 1 and 4 are therefore contingent on an unvalidated spectral shape. Please validate the moment closure against a momentum-resolved solution of Eq. (3) for at least the benchmark point and one additional point, or alternatively demonstrate that the integrated yield is insensitive to the closure by comparing with a different parametrization (e.g., a distribution with a true chemical potential or a spectral-index deformation).","section":"Eq. (2), Eqs. (S5)–(S7), Fig. S2"}],"minor_comments":[{"comment":"I checked the apparent tension between Eq. (10) and the caption of Fig. S1: using the correct Hubble rate H≈2.6×10^-18 GeV at T=1.5 GeV (g_*=70), Eq. (10) gives Γ/H≈6×10^-3 for p=T, so the statement that the number-changing rates stay below H over the plotted range is consistent with Eq. (10). The concern about a rate exceeding H at T_i appears to be based on a numerical error in the Hubble rate.","section":"Eq. (10) vs. Fig. S1"},{"comment":"The second term in the bracket of Eq. (S16) appears to be missing a denominator: as typeset, 2m_φ²(2m_φ²+s) log(...) has dimensions of mass squared rather than being dimensionless. Please check the original LaTeX and ensure the printed formula is correct.","section":"Eq. (S16)"},{"comment":"The active flavor index is suppressed in Eq. (1), but it is not clear whether y_as couples to a single active flavor or to a sum over flavors. This changes the production rate by a factor of the number of flavors and should be stated explicitly, together with the corresponding effect on the relic-density contours.","section":"Eq. (1) and flavor structure"},{"comment":"The text states a Lyman-α lower bound m_s≳8 keV, but Fig. 4 shows viable contours for masses as low as 1 keV. Please indicate clearly in Fig. 4 or its caption which parts of the plotted mass range are excluded by the Lyman-α bound.","section":"Fig. 4 and Lyman-α bound"},{"comment":"The discussion of supernova and DSNB signatures is speculative and would benefit from a brief statement that quantitative analyses are needed before those claims are made.","section":"Potential Impact section"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central idea is worth publishing if the two major points are resolved. The initial-temperature/reheating issue is the more serious one: the present text carries an incorrect claim of insensitivity to T_i, and the parameter space may need to be conditional on a low reheating temperature. The spectral-closure issue is a standard approximation in this literature, but given that the elastic scattering rate is below Hubble for most of the integration, a validation against a momentum-resolved calculation is needed for a quantitative claim. I do not see grounds for rejection, because both issues are fixable within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new ingredient here is the number-changing process ν_a ν_a → ν_s ν_s mediated by a scalar that couples active and sterile neutrinos. That makes production independent of the active-sterile mixing angle, which is a real departure from the DW mechanism and from earlier NSI scenarios that only had active-active or sterile-sterile self-interactions. The paper shows this can open up a parameter region previously excluded by X-ray searches, and it maps out observational prospects. The central logic holds up: the rate calculations in the supplemental material are consistent in the heavy-mediator limit, the DW limit is recovered for small couplings, and the constraints from X-ray, gamma-ray, BBN, and structure formation are folded in sensibly.\n\nThe main soft spot is the two-parameter modified Fermi-Dirac ansatz for the sterile neutrino spectrum, with α and T_s matched to number and energy densities. The paper does not validate this closure against a full QKE solution, and all the rates are evaluated assuming that shape. If freeze-in actually produces a non-thermal spectrum that this ansatz cannot represent, the contours in Figs. 1 and 4 could shift. That is a moderate concern, not a fatal one, but it is the thing I would want to see addressed in revision. Also, the numerical solver is not released, which limits reproducibility.\n\nI checked the stress-test note's claim that the forward rate exceeds Hubble at T_i. It does not hold up. The note uses H ≈ 2.6×10^-21 GeV, but the standard expression H = 1.66√g* T^2 / M_Pl with g* ≈ 70 and T = 1.5 GeV gives H ≈ 2.6×10^-18 GeV. The forward rate from Eq. (10) is about 1.6×10^-20 GeV, so Γ/H ≈ 0.006. The paper's Fig. S1, which shows the rates staying below H, is consistent. The stress-test's arithmetic, not the paper's, is off.\n\nThe paper deserves a serious referee. It is a plausible new production mechanism, the calculations are presented at a level that can be checked, and the phenomenological implications are clearly laid out. Who is it for? DM model builders and sterile neutrino phenomenologists, especially anyone working on X-ray and neutrino observables. I would recommend sending it to peer review, with the request that the authors either provide a QKE sanity check for the spectral closure or at least discuss its limitations, and consider releasing the code.","headline":"Active-sterile NSI producing sterile neutrino DM via the number-changing channel is a genuinely new and plausible mechanism; the main caveat is the unvalidated spectral closure, not the rate-vs-Hubble worry in the stress test.","tokens_in":16878,"tokens_out":3517,"would_cite":true,"duration_ms":36702,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that active-sterile neutrino non-standard interactions can produce sterile-neutrino dark matter at arbitrarily small mixing angles, bypassing the Dodelson-Widrow exclusion.","keywords":["sterile neutrino dark matter","non-standard interactions","freeze-in","Dodelson-Widrow mechanism","number-changing process","scalar mediator","keV-scale dark matter","relic abundance"],"falsifier":"A direct falsifier is to solve the unintegrated Boltzmann equation (3) for $f_s(p)$ without assuming the two-parameter ansatz and recompute the relic density at the benchmark point $m_s=10$ keV, $y_{as}=3\\times10^{-4}$, $m_\\phi=25$ GeV; if the exact momentum distribution changes $\\Omega_s h^2$ by more than the observational uncertainty around $0.12$, the central claim collapses. Less directly, a future X-ray decay-line search that excludes the benchmark region would not falsify the production mechanism but would move the viable contours.","tokens_in":15799,"feed_emoji":"🌌","tokens_out":11699,"duration_ms":104304,"temperature":0.7,"pith_summary":"The paper proposes a previously unexplored production route for sterile-neutrino dark matter: a scalar-mediated non-standard interaction between active and sterile neutrinos. The central number-changing process $\\nu_a\\nu_a \\to \\nu_s\\nu_s$ produces sterile neutrinos at a rate that does not depend on the active-sterile mixing angle, so the observed relic abundance $\\Omega_s h^2 = 0.12$ can be reached even as $\\sin^2 2\\theta \\to 0$. This would bypass the astrophysical exclusion of the standard Dodelson-Widrow mechanism, which requires relatively large mixing angles. If correct, keV-scale sterile neutrinos become viable dark-matter candidates in a broad new parameter region that is testable with future X-ray and gamma-ray telescopes.","feed_headline":"New process makes sterile-neutrino dark matter at zero mixing","feed_subtitle":"A scalar-mediated active-sterile interaction can produce the full dark-matter abundance without fine-tuning.","key_machinery":"The load-bearing mechanism is the number-changing process $\\nu_a\\nu_a \\leftrightarrow \\nu_s\\nu_s$ mediated by the scalar $\\phi$, whose heavy-mediator rate $\\Gamma_{\\nu_a\\nu_a \\to \\nu_s\\nu_s}(p) \\simeq 7\\pi y_{as}^4 p T^4 / (216 m_\\phi^4)$ enters the sterile number- and energy-density equations and supplies the mixing-independent production. To close those equations, the sterile distribution is modeled as a modified Fermi-Dirac form $f_s(p) = \\alpha/(e^{p/T_s}+1)$, with $\\alpha$ and $T_s$ fixed by matching the number and energy densities; $\\alpha$ is the normalization that suppresses the final abundance and $T_s$ controls the spectral shape. The elastic process $\\nu_a\\nu_s \\to \\nu_a\\nu_s$ supplies the thermal potential and damping that modify the residual oscillation production, but it is the number-changing rate that carries the new production path.","core_discovery":"The claim is that the interaction $\\mathcal{L} \\supset y_{as}\\bar{\\nu}_a\\nu_s\\phi + \\mathrm{h.c.}$, with a heavy scalar $\\phi$, opens a number-changing channel $\\nu_a\\nu_a \\leftrightarrow \\nu_s\\nu_s$ whose rate (Eq. 10) is independent of $\\sin^2 2\\theta$. Solving the integrated Boltzmann equations with a modified Fermi-Dirac ansatz for the sterile distribution, the authors find that this channel alone produces the full observed dark-matter abundance for mixing angles that can be arbitrarily small, including exactly zero. The viable region spans sterile masses of roughly $1$--$100$ keV, active-sterile couplings $y_{as}\\sim 10^{-4}$--$10^{-3}$, and mediator masses above about $5$ GeV, while remaining consistent with X-ray, gamma-ray, Lyman-$\\alpha$, BBN and laboratory constraints. No fine-tuned resonance or primordial lepton asymmetry is required.","pith_inferences":["Beyond the paper's sterile-neutrino focus, the mixing-independent freeze-in logic should apply to any number-changing $2\\to 2$ process connecting a thermal bath to an almost empty sector, so similar parameter openings may exist for axion-like particles or other feebly interacting hidden-sector states.","The paper does not solve the full momentum-dependent Boltzmann equation; a dedicated solution that develops spectral features the two-parameter Fermi-Dirac ansatz cannot represent would shift the relic-density contours, and the size of that shift is a quantitative test of the approximation.","In the benchmark window the elastic rate briefly exceeds Hubble while number-changing production stays below it; that kinetic-equilibrium-with-suppressed-chemical-abundance regime could leave a distinctive non-thermal momentum spectrum that future Lyman-$\\alpha$ or 21-cm measurements might distinguish from the standard Dodelson-Widrow spectrum."],"forward_implications":["Sterile-neutrino dark matter can be produced at arbitrarily small active-sterile mixing angles, so the excluded Dodelson-Widrow region is no longer a barrier for keV-scale warm dark matter.","The viable parameter space with $\\Omega_s h^2 = 0.12$ extends to $m_s$ in the $1$--$100$ keV range, $y_{as}\\sim 10^{-4}$--$10^{-3}$ and $m_\\phi\\gtrsim 5$ GeV, where future X-ray missions and beta-decay experiments can search for the decay line and kinematic signatures.","The same scalar interaction gives an energy-dependent opacity for neutrinos passing through dark-matter halos, producing attenuation signatures in astrophysical neutrino spectra.","The $\\nu_a\\nu_a \\to \\nu_s\\nu_s$ process can deplete active neutrinos in core-collapse supernovae, suppressing the diffuse supernova neutrino background relative to standard predictions.","For sub-MeV sterile neutrinos, the process acts as a post-decoupling neutrino cooling channel that mimics self-interacting neutrinos and could ease current cosmological tensions."],"supporting_citations":[{"why":"Defines the Dodelson-Widrow oscillation production mechanism whose excluded parameter space this paper aims to bypass.","marker":"[13]"},{"why":"Supplies the resonant Shi-Fuller alternative requiring a lepton asymmetry, which the new mechanism avoids.","marker":"[23]"},{"why":"Establishes the active-active self-interaction production scenario that still depends on mixing, the main contrast case.","marker":"[32]"},{"why":"Gives the sterile-sterile self-interaction rates and resonance treatment that the collision-term derivations build on.","marker":"[38]"},{"why":"Provides the quantum kinetic equation and Boltzmann formalism used to evolve the sterile neutrino distribution.","marker":"[43]"},{"why":"Sets the observed dark-matter relic abundance that the production contours are matched to.","marker":"[53]"},{"why":"Provides the Lyman-alpha bound on warm dark matter used to set the lower bound on the sterile neutrino mass.","marker":"[68]"},{"why":"Gives the X-ray line constraints that the viable parameter region must evade.","marker":"[14]"},{"why":"Offers a UV-complete model for the effective scalar interaction, supporting the phenomenological setup.","marker":"[42]"}],"fun_headline_variants":["Zero mixing still yields sterile neutrino dark matter","Scalar interaction produces dark matter at zero mixing","New production route for sterile neutrino dark matter","Dark matter from sterile neutrinos without mixing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the sterile-neutrino population, although never in thermal equilibrium, always has the two-parameter modified Fermi-Dirac shape of Eq. (2), with $\\alpha$ and $T_s$ fixed by its number and energy densities; if freeze-in actually produces a differently shaped momentum distribution, the predicted relic abundance and the allowed parameter contours would shift.","fun_headline_variants_meta":{"raw":{"variants":["Zero mixing still yields sterile neutrino dark matter","Scalar interaction produces dark matter at zero mixing","New production route for sterile neutrino dark matter","Dark matter from sterile neutrinos without mixing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000209,"raw_usage":{"total_tokens":1392,"prompt_tokens":912,"completion_tokens":480,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":424}},"tokens_in":528,"tokens_out":480,"duration_ms":5251,"temperature":1.0,"reasoning_tokens":424,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:08:55.647732+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct falsifier is to solve the unintegrated Boltzmann equation (3) for $f_s(p)$ without assuming the two-parameter ansatz and recompute the relic density at the benchmark point $m_s=10$ keV, $y_{as}=3\\times10^{-4}$, $m_\\phi=25$ GeV; if the exact momentum distribution changes $\\Omega_s h^2$ by more than the observational uncertainty around $0.12$, the central claim collapses. Less directly, a future X-ray decay-line search that excludes the benchmark region would not falsify the production mechanism but would move the viable contours.","supporting_citations":[],"review_version":1}