REVIEW 3 major objections 4 minor
A Variable-Slope Smooth-$k$ Filter for Modeling Halo Abundances with Damped and Oscillatory Power Spectra
T0 review · 3 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read A single filter reproduces halo counts for WDM and dark-acoustic-oscillation models
desk verdict The VSMK filter is a genuine analytic extension, but the central claim of a single parameter set is not demonstrated: the parameters are taken from fits to the same simulations, and the scale-inversion premise is asserted rather than proven. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the VSMK window function W_VSMK(k,kM)=[1+(k/kM)^{f(k)}]^{-1}, where f(k) interpolates between the asymptotic slopes beta1 for k/kM << 1 and beta2 for k/kM >> 1. The machinery works through the integrand xi(k,kM)=k^2 P(k) partial(W^2)/partial(kM), whose peak location shifts from k/kM ~ 1 for large halo masses to k/kM << 1 for small masses, mapping filter behavior onto distinct HMF regimes: beta1 sets the small-mass power-law slope dn/dlnM ~ M^{(beta1-3)/3}, while beta2 regulates how strongly DAO oscillations survive in the intermediate-mass HMF.
What would settle it
Compute the exact d sigma^2/dkM integral without the peak approximation for a strong DAO model (h_peak=1) and test whether varying beta2 while holding beta1 fixed leaves the small-mass HMF unchanged, and whether varying beta1 leaves the intermediate-mass oscillation amplitude unchanged; any leakage beyond the claimed few-percent level would falsify the decoupling mechanism.
Extended reading notes
Core claim
Within the extended Press-Schechter formalism, the paper establishes that the asymptotic slope of the halo mass function at small masses is controlled by the low-k/kM behavior of the window function, while the smoothing of dark-acoustic-oscillation features at intermediate masses is controlled by the window at k/kM >~ 1. This scale inversion lets a generalized filter of the form W(k,kM)=[1+(k/kM)^{f(k)}]^{-1}, with f(k)=beta2-(beta2-beta1)[1+(mu k/kM)^delta]^{-1}, decouple the two regimes. Comparing analytic predictions with N-body simulations, the paper shows that a single parameter set reproduces both the WDM small-scale suppression and the DAO oscillations in the halo mass function, with
Load-bearing premise
The scale-inversion mapping — that the small-mass HMF is set by filter behavior at k/kM << 1 and the intermediate-mass DAO imprint by behavior at k/kM >~ 1 — must hold for both WDM and DAO spectra at the chosen mu and delta; if the maximum of the integrand xi(k,kM) does not track these regimes, the two slopes are not truly decoupled.
Editorial extensions
If this is right
- A single VSMK parameter set reproduces halo mass functions for both warm dark matter and DAO-based ETHOS models, removing the need for model-specific SMK recalibration.
- The small-mass slope of the HMF is controlled exclusively by beta1, and the intermediate-mass DAO imprint by beta2, so the two physical regimes can be fitted independently.
- The mass-calibration parameter c continues to shift the HMF along the mass axis without altering its shape or the filter's smoothing of oscillations.
- For ETHOS DAO models at z=5, 8, and 12, the VSMK predictions agree with N-body simulations at the level of Delta <= 0.11 over resolved mass ranges.
- The VSMK filter reduces to the standard SMK filter when beta1=beta2, making it a minimal extension of the existing framework.
Reading between the lines
- If the decoupling holds, future small-scale HMF measurements could separately constrain the free-streaming cutoff (via beta1) and dark-radiation interactions (via beta2), potentially breaking degeneracies that a single-slope filter cannot.
- The same scale-inversion argument may extend to other Press-Schechter observables, such as halo bias or void abundances, whenever the power spectrum is damped and oscillatory.
- Because mu and delta are held fixed across models rather than calibrated, a physical mapping from those parameters to the DAO peak scale or half-mode scale could reduce the filter to its two slope parameters alone.
- The filter is a fitting window rather than a collapse model, so its success suggests that flexible smoothing, not the specific window shape, is what the HMF responds to; simulations reaching lower masses will test whether the beta1 asymptote extrapolates.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces a variable-slope smooth-k (VSMK) filter, Eqs. (3.6)–(3.7), which generalizes the SMK filter by letting the effective logarithmic slope interpolate between β1 at k/kM ≪ 1 and β2 at k/kM ≫ 1. Within the extended Press–Schechter formalism, the paper claims a 'scale-inversion' mechanism: β1 controls the small-mass slope of the halo mass function, dn/dlnM ∝ M^{(β1−3)/3} (Eq. 3.5), while β2 controls the intermediate-mass imprint of dark acoustic oscillations. Using the single parameter set β1 = 4.8, β2 = 3.6, μ = 2.1, δ = 12, c = 3.6, the paper reports agreement with N-body HMFs for a WDM model (Schaeffer & Schneider 2021) and for ETHOS DAO models (Verwohlt et al. 2024; Bohr et al. 2021). The paper concludes that the VSMK filter provides a unified analytic framework for damped and oscillatory power spectra.
Significance. If the claimed decoupling is robust, the VSMK filter is a useful and minimal extension of the SMK filter, with potential applications to non-CDM HMF predictions at small and intermediate masses. The analytic derivation of the small-mass slope in Section 3 and Appendix A is clean and consistent with the SMK limit, and the comparison spans several redshifts and two DAO scenarios. However, the central validation is currently underdetermined: the asymptotic slopes β1 and β2 are taken from the optimal SMK fits to the same simulations used as benchmarks, the transition parameters μ and δ are fixed without a sensitivity analysis, and the scale-inversion premise is only illustrated for two DAO configurations. The manuscript therefore establishes a promising framework, but the advertised universality of the single parameter set is not yet demonstrated.
major comments (3)
- [§3 and Appendix A (Eq. 3.2, Figure 5)] The decoupling claim rests on the assumption that, for kM ≫ khm, the dominant maximum of ξ(k,kM) lies at k/kM < 1/μ, and for kM ≲ khm it lies at k/kM > 1/μ. Figure 5 illustrates this only for m_WDM = 1.61 keV and for kpeak = 35 h/Mpc, hpeak = 1; the paper does not map where the maximum crosses k/kM = 1/μ as a function of kM. Because μ and δ set the transition in units of kM while the DAO peak is at a fixed physical kpeak, the separation of regimes is not guaranteed for other ETHOS parameters. A failure occurs if the dominant ξ-peak lies near the transition region for the mass range where DAO features appear, so that both β1 and β2 contribute simultaneously. Please provide a systematic map of the ξ-peak location as a function of kM/khm and a scan over μ and δ to demonstrate that the one-parameter-per-regime mapping holds.
- [§4 (parameter choice, Table 2, Figure 4)] The 'single parameter set' claim is weakened by the way the parameters are selected. β1 = 4.8 and β2 = 3.6 are the optimal SMK slopes for the WDM and DAO simulations used as benchmarks, and c = 3.6 is likewise the SMK calibration for the DAO model of [41]. The two asymptotic regimes are therefore reproduced partly by construction; the actual test is whether the transition parameters μ and δ are robust and whether the same set works for models not used for calibration. This is further qualified by the text of Figure 4, which states that the Bohr et al. model requires 'minor variations in c and β2.' Please either (a) present a predictive comparison with an uncalibrated model (e.g., hpeak = 0.4 at multiple redshifts without re-fitting β2 or c) or (b) state explicitly which parameters are re-fit and provide a sensitivity analysis of μ and δ.
- [§4 and Data availability statement] The quantitative validation rests on N-body data extracted from published figures, with no published error bars or resolution cuts and no reproducible pipeline. The paper states that 'No new data were generated or analysed' yet presents comparisons to simulation data, and the numerical routines are only 'available from the author upon reasonable request.' This makes it difficult to assess whether the quoted deviations (e.g., ΔVSMK ≤ 0.05 for the DAO case) are robust or sensitive to the extraction procedure. Please make the extracted data and evaluation code available, or at minimum provide a table of the simulation data points and the exact binning/selection criteria used.
minor comments (4)
- [§4, Eq. (4.1) and text below Figure 3] The statement 'Δ3.6(n4.8) → M^{-0.4} − 1 as M → 0' appears incorrect: from Eq. (4.1), Δ3.6(n4.8) = n3.6/n4.8 − 1 ∝ M^{-0.4} − 1, which diverges as M → 0, rather than tending to −1. Please correct the limit or the definition of the deviation.
- [§2.2, Eq. (2.15)] The ETHOS transfer function is difficult to parse as typeset; there are missing parentheses and the term '√h2/4' is ambiguous. Please rewrite the equation with unambiguous notation (e.g., sqrt(h2)/4 and explicit exponentials).
- [Figure 5 caption] The caption labels the two mass scales as 'kM = 30 h Mpc^{-1}' and 'kM = 1000 h Mpc^{-1}' but does not indicate which line style corresponds to which kM; the figure appears to rely on color alone. Please make the legend explicit.
- [§4 and Figure 5] The text uses m_WDM = 0.25 keV for the Schaeffer & Schneider comparison (Figure 3) and m_WDM = 1.61 keV in Figure 5; the relation between these models and the ETHOS parameter mapping (Eq. 2.17) should be stated to avoid confusion.
Circularity Check
Partial circularity: the two VSMK slope parameters are the previously fitted SMK slopes for the same benchmark simulations, and VSMK reduces to those SMK limits by construction in the regimes claimed to be reproduced.
-
fitted input called prediction
[Section 4, Eqs. (3.6)-(3.7), Eq. (3.4), Appendix A]
"WVSMK = [1 + (k/kM)^{f(k)}]^{-1}, f(k)=β2-(β2-β1)[1+(μ k/kM)^δ]^{-1} ... The VSMK filter approaches a SMK filter with β=β1 for k/kM < 1/µ ∼ O(1), and a SMK filter with β=β2 for k/kM > 1/µ as directly follows from Eqs. (3.6) and (3.7). The values of β1 and β2 are chosen to match the optimal SMK slopes reported in [14, 34] and [41], respectively."
The claim 'a single parameter set reproduces both regimes simultaneously' is assembled from two previously fitted slopes: β1 is the optimal SMK slope for the WDM simulation of [14] (4.8), and β2 is the optimal SMK slope for the DAO simulations of [41]/[52] (3.6/3.46). By Eqs. (3.6)-(3.7) and the Appendix quote, VSMK equals SMK(β1) for k/kM < 1/µ and SMK(β2) for k/kM > 1/µ. The WDM small-mass regime is controlled by the k/kM≪1 expansion (Eq. 3.4; HMF slope M^(β1-3)/3), so matching WDM is guaranteed once β1 is the fitted WDM slope; the DAO intermediate imprint is assigned to the k/kM≳1 side where f(k)→β2, so matching DAO is likewise inherited. The remaining parameters µ, δ, c are hand-set transition/calibration parameters, not independent predictions. The comparison is therefore a recombinat
full rationale
The VSMK filter is an explicit ansatz, and the derivation of the small-mass asymptotic slope from the k/kM≪1 expansion is internally consistent. There is no self-citation chain, no imported uniqueness theorem, and no machine-checked claim; the main circularity is the parameter-choice step. β1 and β2 are not predicted or independently fitted here; they are lifted from the optimal SMK fits to the very WDM and DAO simulations that serve as benchmarks. Because the filter is constructed to reduce to SMK(β1) and SMK(β2) in the two regimes that are respectively responsible for the WDM small-mass slope and the DAO intermediate-scale imprint, reproducing both regimes with those values is partly guaranteed by construction. The independent content is limited to the transition parameters µ and δ, which are kept fixed, and the visual/sensitivity check that the two limits coexist without spoiling each other. The asserted scale-inversion mechanism in Section 3 and Appendix A is a physical argument rather than a proof, and the limited number of tested kpeak values (two) and the absence of a µ/δ sensitivity scan affect external validity; these are correctness risks, not additional circularity. Overall, the central 'simultaneous reproduction' claim reduces in part to previously fitted inputs, warranting a partial circularity score of 6.
Assumptions & free parameters
free parameters (5)
- beta_1 =
4.8
- beta_2 =
3.6 (3.46 for the Bohr et al. 2021 case)
- mu =
2.1
- delta =
12
- c =
3.6
assumptions (5)
- domain assumption Press-Schechter/EPS formalism with ellipsoidal collapse and Sheth-Tormen first-crossing distribution (Eqs. 2.2-2.3).
- domain assumption Mass-scale relation M = (4/3) pi rho_bar (c R)^3 for filters whose real-space counterparts have divergent integrals (Eq. 2.11).
- domain assumption The WDM and ETHOS transfer functions (Eqs. 2.13-2.16) faithfully describe the linear power spectra of the simulated models.
- domain assumption For kM >> khm, the variance-derivative integrand is dominated by k/kM << 1, yielding dsigma^2/dkM ~ kM^{-beta_1-1}.
- domain assumption The scale-inversion mapping: small-mass HMF controlled by filter behavior at k/kM << 1, intermediate-mass DAO features by filter behavior at k/kM >~ 1.
Cite this review
Pith. "Pith review of A Variable-Slope Smooth-$k$ Filter for Modeling Halo Abundances with Damped and Oscillatory Power Spectra." pith.science (2026). https://pith.science/paper/TFDS4SFQ
@misc{pith2026260201320,
author = {Pith},
title = {Pith review of: A Variable-Slope Smooth-$k$ Filter for Modeling Halo Abundances with Damped and Oscillatory Power Spectra},
year = {2026},
howpublished = {\url{https://pith.science/paper/TFDS4SFQ}},
note = {Machine review of arXiv:2602.01320}
}
abstract
We introduce a variable-slope smooth-$k$ (VSMK) filter within the Press-Schechter formalism to model halo mass functions derived from damped and oscillatory matter power spectra. While the standard smooth-$k$ approach successfully captures small-scale suppression effects, it intrinsically couples these to oscillatory features at intermediate scales. The VSMK filter generalizes this framework by allowing the effective logarithmic slope of the $k$-space window function to vary smoothly between two asymptotic regimes, thereby decoupling the small-scale suppression of halo abundances from the intermediate-scale oscillatory features characteristic of dark acoustic oscillations (DAO). We compare the analytic predictions obtained with the VSMK filter to $N$-body simulations for warm dark matter and ETHOS-based models with DAO, showing that a single parameter set reproduces both regimes simultaneously. The VSMK filter thus provides a unified and flexible analytic framework for modeling halo abundances in non-cold dark matter scenarios with damped and oscillatory power spectra.
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.