{"id":"f8f3e994-7809-49a9-abe3-a5542e018126","arxiv_id":"2602.01320","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Introduces a variable-slope smooth-k filter within Press-Schechter theory that uses separate slope parameters for small and intermediate mass regimes, claiming a single parameter set matches N-body halo mass functions for both WDM and DAO models.","lead":"This paper proposes a new mathematical filter for estimating how many dark matter clumps (halos) form when the underlying density fluctuations are damped or oscillatory. It aims to let one set of parameters describe both warm dark matter and models with dark acoustic oscillations, which matters for interpreting future small-scale structure surveys.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Scale-inversion premise is asserted, not demonstrated: for DAO models the separation between β1- and β2-controlled mass regimes depends on the fixed kpeak lying well outside the transition region around k/kM ≈ 1/μ; only two kpeak values are tested, and μ, δ are never varied.","rationale":"The central claim has two parts: (i) the asymptotic slope of the HMF is governed by β1 through Eq. (3.5), which is well supported by the k/kM≪1 expansion; (ii) the intermediate DAO imprint is governed by β2, which is the less secure part. The derivation in Section 3 and Appendix A establishes (i) and a qualitative argument for (ii) based on the maximum of ξ. My check targets the assumption that u_max stays cleanly on one side of the transition for all resolved masses and all ETHOS models. If this fails, the VSMK filter does not truly decouple, and the \"single parameter set\" claim is only a fit to the two tested simulations. I do not accuse the author of anything; the paper explicitly says μ, δ are \"kept fixed\" and Figure 4 permits variations, so the universality statement is affirmatively more limited than the abstract implies. I agree with the reader's weakest assumption and keep the CONDITIONAL verdict: the fix is a concrete sensitivity/robustness test, not a rewrite.","tokens_in":16413,"tokens_out":12496,"duration_ms":118010,"concrete_test":"Recompute the variance derivative (Eqs. 3.1–3.2) for the ETHOS transfer function (Eq. 2.15) over a grid of (hpeak,kpeak), e.g., hpeak=0.2,0.4,0.6,0.8,1.0 and kpeak=10–300 h/Mpc, holding the fixed VSMK parameters (β1=4.8, β2=3.6, μ=2.1, δ=12, c=3.6). For each model, locate u_max(kM)=argmax_u ξ(u kM,kM) and record whether u_max < 1/μ for kM>khm and >1/μ for kM 1/μ. The paper asserts this \"scale inversion\" (Section 3, Appendix A) and illustrates it with Figure 5 for two models (m_WDM=1.61 keV and kpeak=35 h/Mpc, hpeak=1), but does not prove it or map where u_max(kM) crosses 1/μ. Because μ and δ set the transition in units of kM while the DAO peak is at fixed physical kpeak, the separation of regimes is not guaranteed for other ETHOS parameters. A failure is possible when the location of the dominant ξ-peak lies near k/kM≈1/μ for the mass range where DAO features appear; then both β1 and β2 enter, and the claimed one-parameter-per-regime mapping fails. The paper's own Section 4/Figure 4 weakens the universality claim by allowing \"minor variations in c and β2\" for the Bohr et al. model, and neither code/data nor a sensitivity scan over μ, δ is provided (Data availability statement). This is the load-bearing weak point: the decoupling claim is not yet established beyond the two calibrated models.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":16801,"tokens_out":7243,"duration_ms":68187,"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":[{"comment":"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.","section":"§3 and Appendix A (Eq. 3.2, Figure 5)"},{"comment":"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 δ.","section":"§4 (parameter choice, Table 2, Figure 4)"},{"comment":"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.","section":"§4 and Data availability statement"}],"minor_comments":[{"comment":"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.","section":"§4, Eq. (4.1) and text below Figure 3"},{"comment":"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).","section":"§2.2, Eq. (2.15)"},{"comment":"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.","section":"Figure 5 caption"},{"comment":"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.","section":"§4 and Figure 5"}],"recommendation":"major_revision","confidential_remarks":"The scale-inversion premise and the parameter-robustness issue are the main obstacles. The paper is not internally inconsistent, but the current presentation overstates the universality of the single parameter set. With a sensitivity analysis over μ and δ, a demonstration of the ξ-peak crossing condition, and a predictive test on an uncalibrated model, the central claim could be made solid. I do not see a reason for rejection if these are supplied."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Briefly: the VSMK filter is a real, useful extension of the SMK filter, and the paper deserves a serious referee. The asymptotic slope derivation leading to dn/dlnM ~ M^(beta1-3)/3 is clean and consistent with the earlier SMK result. The functional form (Eqs. 3.6–3.7) genuinely decouples the small-scale and intermediate-scale filter behavior, which is the right physics for DAO models. The application to DAO spectra is new, and the paper is clearly written and honest about what it does and does not claim.\n\nThe main soft spot is validation. The \"single parameter set\" is assembled from optimal SMK fits to the very simulations used as benchmarks: beta1=4.8 comes from the WDM fit, beta2=3.6 from the DAO fit. So the two asymptotic regimes are partly guaranteed by construction. The transition parameters mu and delta are fixed with no sensitivity analysis, and the scale-inversion mapping in Section 3 and Appendix A is asserted and illustrated but not proven. The stress-test note is right that for other ETHOS parameters, especially kpeak closer to the transition region, the claimed separation between beta1- and beta2-controlled mass regimes could break down. Only two kpeak values are tested, and Figure 4 relaxes beta2 and c for the Bohr et al. model, which undercuts the universality claim. Simulation points are digitized from published figures with no extraction uncertainties, and the code/data availability statement says only that routines are available upon request — no actual release.\n\nNone of this kills the paper; it means the main claim is not yet independently demonstrated. The VSMK filter is a plausible analytic tool, and the derivation of the small-mass asymptote is solid. What is needed: release the routines, run a sensitivity scan over mu and delta, and test the filter on at least one or two additional DAO parameter sets where kpeak lies closer to the nominal transition scale. Then the universality claim will have some teeth.\n\nWho is it for? Anyone doing analytic HMF modeling in non-CDM scenarios. It deserves a serious referee, but with the expectation of substantial revision on the validation side.","headline":"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.","tokens_in":17337,"tokens_out":2032,"would_cite":true,"duration_ms":21129,"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":"A single filter reproduces halo counts for WDM and dark-acoustic-oscillation models","keywords":["dark matter theory","halo mass function","Press-Schechter formalism","smooth-k filter","warm dark matter","dark acoustic oscillations","ETHOS","window function"],"falsifier":"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.","tokens_in":16256,"feed_emoji":"🌌","tokens_out":5890,"duration_ms":48422,"temperature":0.7,"pith_summary":"The paper argues that the standard smooth-k filter cannot simultaneously describe dark matter models whose power spectra combine a small-scale damping cutoff with dark acoustic oscillations, because its single slope couples the two regimes. It introduces the variable-slope smooth-k (VSMK) filter, whose effective logarithmic slope interpolates between a steep slope at low k/kM and a shallower slope at high k/kM. The paper establishes a scale-inversion property: in the Press-Schechter formalism the small-mass halo mass function is fixed by the filter at k/kM << 1, while the intermediate-mass imprint of dark acoustic oscillations is fixed by the filter at k/kM >~ 1. With one parameter set (beta1=4.8, beta2=3.6, mu=2.1, delta=12, c=3.6), the analytic halo mass function matches N-body simulations for both a warm dark matter model and ETHOS-based DAO models at redshifts 5-12. The result matters because it offers a single analytic description of halo abundances across non-cold dark matter scenarios without model-specific recalibration.","feed_headline":"One filter matches WDM and dark-acoustic-oscillation halo counts","feed_subtitle":"A slope-switching window separates small-mass suppression from intermediate-scale DAO wiggles in one halo-mass-function formula.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["One filter, two dark matter regimes: WDM and DAO","Variable-slope filter untangles halo mass suppression from DAO wiggles","New filter decouples small-scale suppression from DAO oscillations","Single filter fits WDM and dark-acoustic halo counts","Halo mass function: one formula for WDM and DAO"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["One filter, two dark matter regimes: WDM and DAO","Variable-slope filter untangles halo mass suppression from DAO wiggles","New filter decouples small-scale suppression from DAO oscillations","Single filter fits WDM and dark-acoustic halo counts","Halo mass function: one formula for WDM and DAO"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000831,"raw_usage":{"total_tokens":3460,"prompt_tokens":736,"completion_tokens":2724,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":480,"completion_tokens_details":{"reasoning_tokens":2634}},"tokens_in":480,"tokens_out":2724,"duration_ms":19157,"temperature":1.0,"reasoning_tokens":2634,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T05:42:15.572911+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}