{"id":"17ccf6c4-3772-4755-aed0-8624e37814ee","arxiv_id":"2501.16908","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Sterile neutrinos suppress the matter power spectrum and halo mass and velocity functions while increasing halo pairwise infall velocities, with effects up to tens of percent for the most massive cases studied.","lead":"This paper runs large cosmological simulations to predict how light sterile neutrinos would change the clustering of matter and dark matter halos. It provides fitting formulas that future surveys could use to search for these particles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central numerical predictions and fitting formulae rest on the linear-response treatment of sterile neutrino clustering (Eqs. 2.8/A.10) at scales where CDM is nonlinear; a full particle-based comparison is needed before the quantitative claims are trusted.","rationale":"The reader's weakest-assumption analysis and my own stress test converge on the same load-bearing point: the linear-response treatment of sterile-neutrino overdensities is the single assumption on which all quantitative results depend. The paper is otherwise careful: it uses refitted cosmological parameters, checks resolution and cosmic variance, and presents fitting formulae derived from four model points. But none of those checks validates the validity of Eq. (2.8)/(A.10) in the strongly nonlinear regime. The fact that the paper explicitly acknowledges the assumption in the Conclusion makes it a named limitation rather than a hidden flaw, and it is a standard method for active neutrinos. However, sterile neutrinos with m_phy = 1–2 eV and ΔN_eff = 0.2–0.4 have larger free-streaming lengths than typical active-neutrino cases, so the linear response should be particularly accurate; conversely, the model with meff = 0.8 eV contributes a non-negligible fraction of the total matter density (f_ν ~ 0.06), and the claimed ~40% suppression is large enough that a few-percent error in the neutrino clustering would change the headline numbers. The proposed test—one particle-based run for model C2—is expensive but decisive, and it directly targets the only untested link in the argument. If the test passes, the conditional acceptance is fully justified; if it fails, the fitting formulae would need recalibration. I therefore keep the reader's CONDITIONAL verdict and recommend no change.","tokens_in":16038,"tokens_out":3666,"duration_ms":37363,"concrete_test":"Run one high-resolution particle-based sterile-neutrino simulation for the most extreme model, C2 (m_phy = 2 eV, ΔN_eff = 0.4), using the same box size (1000 h^-1 Mpc), same refitted cosmology (Table 2), and same seed. Represent sterile neutrinos as tracer particles with Fermi–Dirac thermal velocities and the same total Ω_νs, with at least comparable effective resolution to the CDM particles; alternatively, use an independent implementation of the grid-based method with a nonlinear neutrino solver. Compare z = 0 total-matter P(k), halo-halo pairwise velocity v_hh(r), and HMF against the corresponding linear-response results. If the differences in R(k) = ΔP/P0 exceed ~2–3% at k > 1 h/Mpc, or if v_hh and HMF deviations shift by more than the quoted statistical errors, the fitting formulae (Eqs. 4.4 and 4.13) and the claimed enhancement/suppression amplitudes require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's core quantitative claims—up to ~40% suppression of P(k) and ξ(r), the sign and magnitude of pairwise-velocity deviations, and 40–50% reductions in halo mass and velocity functions—all pass through Eq. (2.8)/(A.10), the linearized Vlasov evolution of δ_νs sourced by the fully nonlinear CDM–baryon field. The method is validated only by resolution convergence (Fig. 1) and a two-seed cosmic-variance check (Fig. 11), neither of which tests the linear-response approximation itself. At k > 1 h/Mpc, where the paper claims the m_phy–ΔN_eff degeneracy is broken and where R̄ is averaged over [0.7, 2.5] h/Mpc, CDM fluctuations are strongly nonlinear; the sterile-neutrino overdensity is assumed to remain linear and to respond through the free-streaming kernel Φ(k(s−s′)). If true sterile-neutrino clustering on these scales deviates from that linear solution—for example, because nonlinear potential wells induce extra neutrino overdensities, or because the linear response misestimates the potential correction in Eq. (2.7)—then the suppression amplitude, the pairwise-velocity enhancement, and all three fitting formulae are systematically biased. The Conclusion itself flags the assumption ('the sterile neutrino over-density responds to the nonlinear CDM plus baryon density according to the linear theory'), so the concern is not hidden; it is simply untested against a full particle-based sterile-neutrino simulation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the impact of eV-scale sterile neutrinos on cosmological large-scale structure using N-body simulations that incorporate both active and sterile neutrinos through a grid-based linear response approximation (LRA). The authors run simulations for four sterile-neutrino models defined by (m_phy, ΔN_eff), with cosmological parameters refitted to Planck+BAO data, and compare against a no-sterile-neutrino fiducial run. They report that sterile neutrinos suppress the total matter power spectrum and the two-point correlation function, reduce the halo mass and maximum-circular-velocity functions, and increase the magnitude of the halo-halo pairwise velocity. They also present three fitting formulae for the averaged fractional deviations of the power spectrum, halo pairwise velocity, and pairwise velocity dispersion as functions of m_phy and ΔN_eff.","tokens_in":16439,"tokens_out":10361,"duration_ms":98852,"significance":"If the quantitative results are reliable, the paper provides a useful set of predictions and fitting formulae for using LSS observables to constrain eV-scale sterile neutrinos, and it highlights a physically interesting breaking of the m_phy–ΔN_eff degeneracy at small scales. The work is strengthened by the use of refitted cosmological parameters, consistency checks with two random seeds, and a simulation method that is standard in the massive-neutrino literature. However, the fitting formulae are calibrated on only four simulation models with three free parameters, and the linear-response approximation is not directly validated in the high-k regime where the main quantitative claims are made; these issues limit the strength of the quantitative conclusions.","major_comments":[{"comment":"The fitting formulae for \\bar R, \\bar R^v_hh, and \\bar R^σ_hh are calibrated on only four simulation models (T1–T4 or B1–C2) while each formula has three free coefficients. With only one degree of freedom, the quoted parameter uncertainties are not statistically meaningful, and the expressions are effectively interpolations through four points rather than validated fitting functions. The authors should either run additional models spanning the (m_phy, ΔN_eff) plane, or explicitly present Eqs. (4.4) and (4.13) as four-point interpolations and remove or heavily qualify the reported coefficient errors. This matters because the fitting formulae are a central deliverable of the paper.","section":"§4.1, Eq. (4.4); §4.2, Eq. (4.13) and Table 3"},{"comment":"The quantitative predictions at k ∈ [0.7, 2.5] h/Mpc and the halo statistics in §4 all pass through the linear-response evolution of the sterile-neutrino overdensity, Eq. (2.8)/(A.10), sourced by the fully nonlinear CDM-baryon field. The validation presented in Fig. 1 and Fig. 11 tests only resolution convergence and seed dependence, not the validity of the linear-response approximation itself in the k range where the paper claims degeneracy breaking and reports up to ~40% suppression. I request either a direct comparison against a particle-based sterile-neutrino simulation in this regime, or a quantitative estimate of the LRA error from published tests at k > 1 h/Mpc. Without this, the accuracy of the headline suppression amplitudes and of the fitting formulae is not established.","section":"§2.3, Eq. (2.8); Appendix A.1; §3.1, Fig. 1"},{"comment":"The sign reversal between the particle-particle pairwise velocity, whose magnitude decreases with m_phy and ΔN_eff, and the halo-halo pairwise velocity, whose magnitude increases for r > 4 h^-1 Mpc, is a key qualitative result. The paper does not explain this reversal. If it arises from halo bias or from the specific halo mass range [10^13, 10^14] h^-1 M_sun, this should be demonstrated, because it directly affects the interpretation of the fitting formula in Eq. (4.13).","section":"§4.2, Figs. 6 and 7"}],"minor_comments":[{"comment":"The lower panel is labelled 'P(k)/P0(k)' but the plotted quantity is the fractional deviation ΔP/P0; please relabel the axis to avoid confusion.","section":"Fig. 3"},{"comment":"The caption reads '∆Nphy = 0.4'; this should be 'ΔN_eff = 0.4'.","section":"Fig. 11 caption"},{"comment":"The fitting formulae are presented as functions of m_phy and ΔN_eff, but the simulations also vary the other cosmological parameters via the Planck+BAO refit. The text should state explicitly that Eqs. (4.4) and (4.13) describe the combined effect of sterile neutrinos plus the accompanying refit, not the isolated free-streaming effect at fixed background cosmology.","section":"§4.1 and §4.2"},{"comment":"The redshifts at which the pairwise velocity, halo mass function, and velocity function results are evaluated are not stated in the figure captions; please specify that the results are at z = 0 (or state the relevant redshift in each caption).","section":"§4.2 and §4.3"},{"comment":"A brief sentence in §3.1 noting the expected accuracy of the linear-response approximation at the k values used, with a citation to the relevant validation studies, would help the reader assess the robustness of the results without requiring a new simulation.","section":"§2.3"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful and honest incremental paper. It extends the linear-response grid method used for active neutrinos to eV-scale sterile neutrinos, and gives quantitative predictions for P(k), 2PCF, pairwise velocities, and halo mass/velocity functions, including fitting formulae. The genuinely new pieces are the sterile-specific templates and the demonstration that m_phy and DeltaN_eff are not degenerate even at fixed meff. The simulations use standard tools and the Planck+BAO refit for each model is a nice touch; the two-seed check shows the P(k) results are stable.\n\nThe main caveat is exactly the one the conclusion flags: sterile neutrino overdensity is evolved with the linear response equation (2.8/A.10) while CDM goes nonlinear. Validation is resolution convergence and a two-seed run, not a particle-based sterile neutrino simulation. I don't think this is a fatal flaw—the LRA has been tested in the active neutrino context and the physics is similar—but the k>1 h/Mpc predictions and fitting formulae are method-dependent until checked against a full particle simulation. If I were refereeing, I would ask for that comparison or a clear statement of the regime of validity.\n\nMinor issues: the Rbar fit uses four models with three free parameters, so it is not a strong test of the functional form; the halo sample requires at least 200 particles while the stated mass range [1e13,1e14] includes halos below the implied completeness limit; and several of the modeled cases (meff=0.4, 0.8 eV) sit above existing LSS bounds, though the paper is upfront that these are illustrative. The self-citations are to the group's earlier active-neutrino work and are method reuse, not a problem.\n\nBottom line: the central conclusions are well supported. This deserves a serious referee. Send it to review and ask for an LRA validation test before acceptance. I'd cite it if I worked on sterile neutrino constraints.","headline":"Useful, incremental sterile-neutrino LSS paper with honest fitting formulae; the main caveat is that the linear-response treatment is untested against particle simulations on nonlinear scales.","tokens_in":16922,"tokens_out":3834,"would_cite":true,"duration_ms":36500,"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":"Light eV-scale sterile neutrinos would leave a two-sided imprint on the cosmic web: suppressed matter clustering and halo abundances, but enhanced halo pairwise infall velocities, quantified by new fitting formulae.","keywords":["sterile neutrinos","large-scale structure","matter power spectrum","halo mass function","halo velocity function","pairwise velocity","N-body simulations","linear response approximation"],"falsifier":"Run the same five cosmologies with a full particle-based sterile neutrino population at matched box size and resolution, and compare the matter power spectrum, halo pairwise velocity, and halo mass and velocity functions at $k \\gtrsim 1\\,h\\,{\\rm Mpc}^{-1}$; a disagreement larger than the quoted cosmic variance would show the linear-response predictions are not accurate there.","tokens_in":15886,"feed_emoji":"🌌","tokens_out":8846,"duration_ms":78000,"temperature":0.7,"pith_summary":"The paper sets out to show that a relic population of eV-scale sterile neutrinos, parameterised by a physical mass $m_{\\rm phy}$ and an extra radiation abundance $\\Delta N_{\\rm eff}$, leaves measurable imprints on the cosmic large-scale structure. Using N-body simulations in which sterile neutrinos are evolved as a linearly responding density field, it finds that the total matter power spectrum and the two-point correlation function are suppressed by up to roughly 40 percent, and the halo mass and circular-velocity functions by 40–50 percent at their high ends. At the same time, the infall speed of halo pairs is enhanced by up to about 15 percent, while the pairwise velocity dispersion drops by a few percent. The paper also provides fitting formulae that turn these fractional deviations into functions of $m_{\\rm phy}$ and $\\Delta N_{\\rm eff}$, so observations could in principle reconstruct the sterile neutrino parameters.","feed_headline":"Sterile neutrinos damp cosmic structure up to 40 percent","feed_subtitle":"eV-scale sterile neutrinos would cut matter clustering and halo counts while making halo pairs fall faster—a signature surveys can hunt.","key_machinery":"The load-bearing machinery is the linear-response approximation for the sterile neutrino overdensity, implemented on a grid inside a modified N-body code. In this scheme the sterile neutrino phase-space distribution is split into an unperturbed Fermi-Dirac part and a linear perturbation, and the linearised Vlasov equation (Eq. 2.8 / A.10) evolves $\\tilde{\\delta}_{\\nu_s}(s,k)$ using the free-streaming kernel $\\Phi(q)$, which is the Fourier transform of the normalised momentum distribution. This overdensity is then folded into the total matter density field through Eq. (2.7), so the gravitational force on CDM particles includes the neutrino smoothing without requiring neutrino particles. The paper's fitting formulae are quadratic expansions in $m_{\\rm phy}$ and $\\Delta N_{\\rm eff}$ around the no-sterile fiducial model, with coefficients fixed by the simulations; that expansion is the object that converts simulation outputs into a usable observable prediction.","core_discovery":"The central discovery, on the paper's own terms, is that light sterile neutrinos act as a free-streaming hot component whose gravitational back-reaction reshapes structure in a characteristic two-sided way. The dilution of the total matter overdensity by the weakly clustering sterile neutrinos suppresses power at $k \\gtrsim 0.1\\,h\\,{\\rm Mpc}^{-1}$ and lowers the abundance of massive halos, yet the same slower growth makes infalling halo pairs fall faster at separations above a few Mpc, because the surrounding matter distribution is less clumped and the pairwise streaming relation responds to the changed correlation function. The paper claims that these effects are robust across its resolution tests and that the parameter degeneracy between $m_{\\rm phy}$ and $\\Delta N_{\\rm eff}$ is broken at $k \\gtrsim 1\\,h\\,{\\rm Mpc}^{-1}$ in the power spectrum and at $r \\lesssim 4\\,h^{-1}\\,{\\rm Mpc}$ in the correlation function. If correct, the fitting formulae it presents let a measurement of $\\bar{R}$, $\\bar{R}^v_{hh}$, or $\\bar{R}^{\\sigma}_{hh}$ be translated directly into constraints on sterile neutrino mass and thermalisation.","pith_inferences":["By extension, the same grid-based pipeline applies to any decoupled relic with a Fermi-Dirac-like momentum distribution and small overdensity, such as thermally produced eV-scale QCD axions; the paper notes this, and the natural extension is that the quadratic fitting-formula structure would carry over with rescaled coefficients.","Since the power spectrum breaks the $m_{\\rm phy}$–$\\Delta N_{\\rm eff}$ degeneracy at small scales while the pairwise velocity is nearly degenerate in $m_{\\rm eff}$, combining the two observables in a joint analysis should constrain the sterile neutrino parameters more tightly than either channel alone.","A testable extension is to measure the predicted roughly 15 percent halo pairwise velocity enhancement at $r\\sim6\\text{–}20\\,h^{-1}\\,{\\rm Mpc}$ with kinematic Sunyaev-Zeldovich or redshift-space streaming data, which current and near-future surveys have the pair counts to attempt."],"forward_implications":["The predicted suppression of the matter power spectrum and two-point correlation function reaches roughly 40 percent for $m_{\\rm phy}=2\\,{\\rm eV}$, $\\Delta N_{\\rm eff}=0.4$, and the $m_{\\rm phy}$–$\\Delta N_{\\rm eff}$ degeneracy is broken at $k\\gtrsim1\\,h\\,{\\rm Mpc}^{-1}$.","Halo mass and cumulative circular-velocity functions are suppressed by 40–50 percent at their high-mass and high-speed ends for the same parameters, so cluster counts and velocity-selected samples become sensitive probes.","For halos in the mass range $[10^{13},10^{14}]\\,M_\\odot\\,h^{-1}$, the halo pairwise infall velocity increases by up to about 15 percent while its dispersion decreases by about 2 percent, giving a velocity-space signature opposite in sign to the density suppression.","Eq. (4.4) and Eq. (4.13) with the tabulated coefficients provide direct fits for the averaged fractional deviations, meaning that a measured deviation can be mapped back to a region in the $m_{\\rm phy}$–$\\Delta N_{\\rm eff}$ plane.","Because the background cosmology is refitted for each sterile neutrino model, the quoted impacts are those that would survive a joint CMB+BAO calibration; they are not artefacts of holding other parameters fixed."],"supporting_citations":[{"why":"Supplies the linear-response scheme that computes neutrino overdensities from the nonlinear CDM field.","marker":"[23]"},{"why":"Provides the grid-based neutrino N-body method that this paper extends to sterile neutrinos.","marker":"[48]"},{"why":"Gives the early-universe 3+1 active-sterile mixing calculation behind the $\\Delta N_{\\rm eff}$-parametrised sterile neutrino distribution.","marker":"[47]"},{"why":"Furnishes the CMB data and fiducial cosmological parameters used both for the fixed-cosmology tests and the MCMC refits.","marker":"[11]"},{"why":"Supplies the BAO data included in the parameter refits for each sterile neutrino model.","marker":"[26]"},{"why":"Computes the linear CDM and neutrino power spectra used to set initial conditions at $z=99$.","marker":"[51]"},{"why":"Performs the MCMC parameter estimation that yields the refitted cosmologies in Table 2.","marker":"[52]"},{"why":"Identifies the halos whose mass and velocity functions and pairwise velocities are analysed.","marker":"[53]"}],"fun_headline_variants":["Sterile neutrinos damp power but speed halo infall","Light sterile neutrinos leave twin cosmic fingerprints","Hot sterile neutrinos: less clump, faster fall","How eV sterile neutrinos reshape structure growth","Sterile neutrinos slow clustering, boost pair velocities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the sterile neutrino overdensity follows the linear evolution equation even on scales where CDM is strongly nonlinear, so its clustering can be computed from the nonlinear CDM field rather than from a full particle treatment; the paper tests this by resolution convergence but not against a particle-based sterile neutrino simulation.","fun_headline_variants_meta":{"raw":{"variants":["Sterile neutrinos damp power but speed halo infall","Light sterile neutrinos leave twin cosmic fingerprints","Hot sterile neutrinos: less clump, faster fall","How eV sterile neutrinos reshape structure growth","Sterile neutrinos slow clustering, boost pair velocities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000216,"raw_usage":{"total_tokens":1429,"prompt_tokens":937,"completion_tokens":492,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":418}},"tokens_in":553,"tokens_out":492,"duration_ms":5142,"temperature":1.0,"reasoning_tokens":418,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T05:42:46.106085+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same five cosmologies with a full particle-based sterile neutrino population at matched box size and resolution, and compare the matter power spectrum, halo pairwise velocity, and halo mass and velocity functions at $k \\gtrsim 1\\,h\\,{\\rm Mpc}^{-1}$; a disagreement larger than the quoted cosmic variance would show the linear-response predictions are not accurate there.","supporting_citations":[],"review_version":1}