{"id":"8ca9fd82-2030-4c96-a671-6b87a3c70630","arxiv_id":"2504.21031","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"A limit on warm dark matter from Milky Way satellite halos loses most of its force once non-linear regeneration, stripping, and baryon effects are added, leaving the rotation-curve measurement consistent with the revised estimate.","lead":"This paper argues that a published lower limit on the warm dark matter expansion parameter is too strong because it ignores non-linear structure regeneration, halo stripping, and baryon cold-dark-matter behavior. A generalist reader might care because it reopens whether light warm dark matter is consistent with the observed population of Milky Way satellite galaxies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The revised limit a_hNR ≲ 6e-7 is not produced by the paper's hydrodynamics; it comes from treating a galaxy-count ratio read from the author's prior papers as the WDM halo suppression factor in Eq. (25), with no error budget. The reconciliation is therefore unsupported until that step is validated.","rationale":"Good-faith reading: the paper's qualitative claim—that published WDM limits from galaxy counts ignore non-linear regeneration, stripping, and baryons—is plausible and is supported by cited simulations, including Schneider et al. and Paduroiu et al. The paper also honestly flags its own limitations, notably in Sections 7 and 8. The problem is that the quantitative headline number ahNR ≲ 6e-7 is not an output of the spherical hydrodynamics; it is obtained by substituting a galaxy-count ratio into a halo-mass-function suppression formula with no uncertainty. The reader's concern about 1D spherical initial conditions is related but secondary: even if those integrations are accepted as illustrative, they do not by themselves produce the revised limit. The weakest load-bearing step is the unvalidated use of ndata/nΛCDM in Eq. (25). A full cosmological hydro simulation at the inferred mass is the appropriate test. The verdict should remain conditional: the direction is plausible, but the magnitude is not yet established.","tokens_in":12184,"tokens_out":10286,"duration_ms":106875,"concrete_test":"Run a cosmological hydro simulation of a Milky-Way-mass volume in WDM with mh ≈ 0.4 keV (the value inferred in Section 5), including baryons and a subgrid galaxy-formation model, and compare the predicted number of satellites with Mh > 5.4e8 M⊙ with the observed census used in limit (4). Also compute ndata/nΛCDM from that simulation so Eq. (25) is checked rather than assumed. If the simulated satellite count falls below the observed count at 2σ, the revised bound ahNR ≲ 6e-7 is falsified; if it matches, the reconciliation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is not the spherical collapse itself but the conversion of observed galaxy counts into a WDM mass. In Section 5 the paper estimates ndata/nΛCDM ≈ 0.01 at MPS = 1e9 M⊙ from Figures 1, 2 and 4 of the author's ref [26], then inserts this into Eq. (25), a fitting formula from Schneider et al. for the suppression of the WDM halo mass function, to obtain Mfs ≈ 9e10 M⊙, vhrms(1) ≈ 170 m/s, and ahNR ≈ 6e-7. This is the number that makes the discrepancy 'not significant.' The ratio ndata/nΛCDM is a galaxy-count ratio, not a dark-matter halo ratio; it includes star-formation efficiency, feedback, and survey completeness, none of which are modeled or assigned uncertainties. Using it as the left side of Eq. (25) assumes all baryonic and observational effects cancel in the ratio. If they do not, the inferred Mfs and hence ahNR are biased. The paper's own stripping calculation (Figure 5) does not bridge this gap: it lowers M200 only by a factor 0.29, from 3e10 to 8e9 M⊙, still far above the 5.4e8 M⊙ threshold, and the statement that stripped galaxies extend to zero mass is asserted, not computed. Finally, the baryonic tail calculation invoked for agreement with observations is explicitly deferred in Section 7 ('the calculation is lengthy and will be omitted'), so the n ≈ 0.86 of Eq. (31) is not independently checkable from this paper. The combination of an unvalidated galaxy-count ratio and an omitted calculation leaves the central reconciliation unsupported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper re-examines the warm dark matter limit a_hNR < 6e-8 obtained from the minimum halo mass hosting Milky Way satellites (Eq. 4) and attempts to resolve its discrepancy with the measurement a_hNR = (1.39 +/- 0.24)e-6 from dwarf galaxy rotation curves (Eq. 5). The author argues that the limit neglects nonlinear regeneration of the power spectrum, stripping of halos by neighboring galaxies, and baryons that behave as cold dark matter. Using spherically symmetric hydrodynamic integrations (Section 4), an empirical tail in the transfer function (Eq. 27), and a galaxy-count matching procedure (Eq. 25), the paper obtains a revised estimate a_hNR about 6e-7 and concludes that there is no longer a significant discrepancy. Sections 5 and 6 develop the nonlinear regeneration and stripping arguments, Section 7 discusses baryons, and Section 8 draws the conclusions.","tokens_in":12552,"tokens_out":5768,"duration_ms":54866,"significance":"If the revised calculation were correct, it would matter: it would relax warm dark matter constraints from Milky Way satellites by about an order of magnitude and bring them closer to the dwarf rotation curve measurement, while highlighting nonlinear regeneration, stripping, and baryonic perturbations as previously neglected effects. The author is commendably explicit about the assumptions and limitations, including the spherical symmetry, the single Gaussian initial overdensity, and the omitted numerical integration, so the reader can see exactly where the argument needs support. However, the central quantitative estimate is not established: the galaxy-to-halo conversion in Eq. (25) has no error budget and no baryonic correction, the key tail exponent of Eq. (31) comes from an omitted calculation, and the stripping claim is not computed down to the relevant mass range. The paper is therefore a plausible programmatic argument rather than a completed derivation, and the abstract's central claim is not yet supported by the presented numbers.","major_comments":[{"comment":"The revised limit a_hNR ≈ 6e-7 is obtained by inserting a galaxy-count ratio n_data/n_LCDM ≈ 0.01 at M_PS = 1e9 M_sun, read from Figures 1, 2 and 4 of the author's ref. [26], into the warm dark matter halo suppression formula of Schneider et al. But n_data is a galaxy count, not a dark matter halo count, and no justification is given that star formation efficiency, feedback, and survey completeness cancel in the ratio. Without an error budget and a baryonic correction, the inferred M_fs, v_hrms(1), and a_hNR are unsupported, and this step is load-bearing for the entire reconciliation.","section":"§5, Eq. (25)"},{"comment":"Even taking the revised estimate at face value, the paper's own numbers do not support the claim that there is no longer a significant discrepancy. The revised value a_hNR ≈ 6e-7 is an upper limit, while the measurement (5) is (1.39 ± 0.24)e-6, which lies about 3.3 sigma above the revised limit. The abstract's central claim is therefore contradicted by the quantitative result unless an explicit uncertainty on the revised estimate is provided.","section":"§5 and §8"},{"comment":"The value n ≈ 0.86, which the paper uses to bring the baryonic tail into agreement with the measured range 0.5 ≲ n ≲ 1.1, is stated to come from a numerical integration that is explicitly omitted from the paper: 'the calculation is lengthy and will be omitted'. Because Eq. (31) is used as part of the argument that the tail reconciles predictions with observations, the calculation must be reported, tabulated, or made available in a supplement; otherwise the central reconciliation is not checkable.","section":"§7, Eq. (31)"},{"comment":"The claim that stripped-down galaxies extend all the way to zero mass is asserted rather than computed. The single spherically symmetric stripping simulation in Figure 5 reduces M_200 by a factor 0.29, from 3e10 to 8e9 M_sun, which is still more than an order of magnitude above the threshold M_h < 5.4e8 M_sun of Eq. (6). The figure therefore does not demonstrate that stripped satellites can populate the minimum-mass bin that drives the original limit (4).","section":"§6, Figure 5"},{"comment":"All hydrodynamic simulations are spherically symmetric and initialized with a single Gaussian overdensity at z_i = 65.9. No comparison with three-dimensional simulations is provided, so it remains untested whether these idealized initial conditions faithfully capture the formation, stripping, and survival of Milky Way satellite halos in warm dark matter. If they do not, the conclusion that baryons and stripping erase the discrepancy does not follow; a concrete test would be a comparison with an existing 3D warm dark matter simulation at comparable mass scales.","section":"§4, Figures 1–5"}],"minor_comments":[{"comment":"The text contains several typographical artifacts, including 'disa greement', 'under stand', and 'spe ctrum' in the Abstract and Introduction, which should be corrected.","section":"Abstract and Introduction"},{"comment":"Equation (25) has a mismatched parenthesis: 'n_LCDM (MPS))' should be 'n_LCDM (MPS)'.","section":"Eq. (25)"},{"comment":"The table header spells 'Observble' instead of 'Observable', and the symbol 'greaterorsimilar' appears as raw LaTeX rather than a typeset relation.","section":"Table 1"},{"comment":"The phrase 'galaxies that are stripped by their neighbors' is followed by the word 'coaless' in the text; this appears to be a misspelling of 'coalesce'.","section":"§6"},{"comment":"The axis labels in the figures are difficult to read because of the embedded LaTeX text; the figures should be regenerated with clear typeset labels for the densities and velocities of dark matter and baryons.","section":"Figures 1–5"}],"recommendation":"major_revision","confidential_remarks":"The paper relies heavily on the author's own prior publications (refs. [15, 16, 25, 26, 35]) for the data inputs and for the measured tail exponent, and the numerical derivation behind Eq. (31) is explicitly omitted. Editors may wish to weigh whether the central quantitative claim is sufficiently standalone for this journal, and whether the galaxy-count versus halo-count mismatch in Eq. (25) is acceptable without an independent calibration."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful part of this paper is the argument that the published limit ahNR < 6e-8 from Milky Way satellite counts (Eq. 4) is not a clean linear-theory bound. The author is right that non-linear regeneration, stripping, and baryons acting as cold dark matter are known effects that a robust constraint would need to model, and the paper does a service by forcing that conversation. What is new is the application of these established effects to this specific limit, not the effects themselves. The spherical hydro runs are explicit, and the omission of several key calculations is at least flagged rather than hidden.\n\nThe problem is the number. The revised limit ahNR ≲ 6e-7 is obtained by taking ndata/nΛCDM ≈ 0.01 at MPS = 1e9 from the author's own Figures 1, 2, 4 in ref [26], dropping it into Eq. (25) as if it were a dark-matter halo suppression factor, and reading off Mfs. That is a galaxy-count ratio, not a halo-mass-function suppression. Star-formation efficiency, feedback, survey completeness, and their uncertainties are not modeled. No error bar is attached to the ratio or to the derived ahNR. The measurement being defended, Eq. (5), is the author's own, as is the tail index n in ref [25]; self-citation is not a flaw in itself, but here the data being matched and the parameters that achieve the match come from the same line of work, so the agreement carries little independent weight.\n\nThere are two smaller soft spots. Section 7 explicitly says the baryonic tail calculation 'is lengthy and will be omitted', yet Eq. (31) then quotes n ≈ 0.86 from 'a numerical integration' with no details. And the stripping simulation in Figure 5 only lowers M200 from 3e10 to 8e9 solar masses, which is still 15 times above the 5.4e8 threshold being tested; the claim that stripped galaxies extend to zero mass is asserted, not computed. Together these make the quantitative conclusion essentially uncheckable from the manuscript.\n\nI don't think the author is being dishonest. The paper tells you where it is being rough, and the qualitative direction is consistent with the cited simulation literature. But as written, the 6e-7 number does not follow from the demonstrated work. It is a plausible qualitative reconciliation presented with a false precision.\n\nIf this lands on my desk, I would send it to a referee who knows the WDM small-scale literature, because the question matters and the literature genuinely has a tension here. But I would expect the referee to demand the omitted derivation, an error budget on ndata/nΛCDM, and a demonstration that the spherical stripping scenario is representative before the quantitative claim becomes credible.\n\nYours,","headline":"The paper's qualitative challenge to the satellite-count WDM limit is worth taking seriously, but the revised 6e-7 number is a consistency estimate, not a demonstrated result.","tokens_in":13124,"tokens_out":5251,"would_cite":false,"duration_ms":49432,"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":"This paper argues that the tight published limit on warm dark matter's non-relativistic expansion parameter is invalid because it ignores regeneration, stripping, and baryons, and that with these included the measurement and limit no…","keywords":["warm dark matter","free-streaming cutoff","non-linear power spectrum regeneration","stripped-down galaxies","baryon perturbations","Milky Way satellites","halo mass function","dark matter transfer function"],"falsifier":"A full three-dimensional cosmological hydrodynamic simulation of warm dark matter with $m_h = 0.15$ keV and $v_{\\rm hrms}(1) = 493$ m/s, using the standard free-streaming cutoff without an added empirical tail, should produce a population of stripped-down Milky Way satellites with halo masses below $M_h = 5.4\\times 10^8\\,M_\\odot$ and large cores at $z=0$. If instead the minimum surviving satellite halo mass in such a simulation stays above $5.4\\times 10^8\\,M_\\odot$, the paper's reconciliation fails.","tokens_in":11895,"feed_emoji":"🌌","tokens_out":8229,"duration_ms":71984,"temperature":0.7,"pith_summary":"This paper argues that the published limit $a_{\\rm hNR} < 6\\times 10^{-8}$ from the minimum halo mass hosting Milky Way satellites is not valid as stated, because it leaves out three phenomena: non-linear regeneration of the small-scale power spectrum, stripping of galaxy halos by neighbors, and baryons that act as cold dark matter after decoupling. Including these effects, the author estimates the limit weakens to $a_{\\rm hNR} \\lesssim 6\\times 10^{-7}$, within a factor of roughly 2.3 of the dwarf-galaxy rotation-curve measurement $a_{\\rm hNR} = (1.39 \\pm 0.24)\\times 10^{-6}$. The paper concludes that the measurement is not excluded beyond reasonable doubt and that all galaxy-based warm dark matter limits should incorporate the three effects. If right, the apparent contradiction between the satellite-count limit and rotation-curve measurements disappears.","feed_headline":"Warm dark matter limit may be too strong by tenfold","feed_subtitle":"Including regeneration, halo stripping and baryons erases the gap with rotation-curve measurements.","key_machinery":"The load-bearing object is the warm dark matter transfer function $\\tau^2(k)$, the ratio of the warm to cold dark matter linear power spectra. The paper starts from the standard Gaussian cutoff $\\tau^2(k) = \\exp(-k^2/k_{\\rm fs}^2)$ and then extends it with a power-law tail $\\tau^2(k) = \\exp(-k^n/k_{\\rm fs}^n)$ for $k \\geq k_{\\rm fs}$, with $0.5 \\lesssim n \\lesssim 1.1$ measured and $n \\approx 0.86$ obtained from the baryon calculation. This tail is the device that encodes non-linear regeneration, stripping, and the cold behavior of baryons. The dynamical argument runs through spherically symmetric hydrodynamical integrations starting from a Gaussian overdensity at $z_i = 65.9$; a galaxy forms when the combined warm dark matter and baryon density perturbation reaches 1.69 times the mean density, the standard collapse threshold. The resulting cored isothermal sphere has a core radius $r_c$ and an $r_{200}$ halo mass $M_h$ determined by the constancy of $\\rho_h(r) r^2$, and these quantities give the halo masses compared with the observed satellite minimum.","core_discovery":"On the author's own terms, the central discovery is that the bound $a_{\\rm hNR} < 6\\times 10^{-8}$ is an artifact of a linear-theory extrapolation. The bound rests on the claim that small density perturbations at wavevectors $k > k_{\\rm fs}$ are exponentially suppressed by free-streaming, so a galaxy halo as small as $M_h < 5.4\\times 10^8\\,M_\\odot$ should not exist for warm dark matter with a small particle mass. The paper attempts to show that this is wrong in three ways: the linear spectrum acquires a regenerated non-linear tail; most small galaxies are stripped remnants of larger halos, with a mass distribution extending to zero; and baryons, which decouple cold, seed perturbations below $k_{\\rm fs}$. Spherically symmetric hydrodynamical integrations with $m_h = 0.15$ keV and $v_{\\rm hrms}(1) = 493$ m/s produce halos with $M_h = 3\\times 10^8\\,M_\\odot$, matching the scale of the observed satellites, and a simulated stripped-down galaxy loses enough mass to drop from $M_h = 3\\times 10^{10}\\,M_\\odot$ to $8\\times 10^9\\,M_\\odot$. The conclusion is that the minimum halo mass quoted by the limit is not excluded.","pith_inferences":["This reconciliation implies that warm dark matter with $a_{\\rm hNR}$ near $10^{-6}$ may be nearly indistinguishable from cold dark matter for galaxy formation above $M_{\\rm PS} \\approx 5\\times 10^8\\,M_\\odot$; the discriminating power would then shift to smaller masses or to the internal structure of cores, such as the relation between core radius and circular velocity.","If the baryon-seeded tail is as strong as claimed, the free-streaming cutoff may not be the limiting scale for the first galaxies; instead the baryon perturbation spectrum below $k_{\\rm fs}$ sets the earliest formation epoch, a prediction that could be checked in high-redshift galaxy ultraviolet luminosity functions.","The author's revised limit of $a_{\\rm hNR} \\lesssim 6\\times 10^{-7}$ still sits a factor of about 2.3 below the rotation-curve measurement; this residual gap is not explained, and a quantitative stripped-satellite abundance model could either close it or reopen the discrepancy."],"forward_implications":["The Milky Way satellite bound weakens from $a_{\\rm hNR} < 6\\times 10^{-8}$ to $a_{\\rm hNR} \\lesssim 6\\times 10^{-7}$, a factor of about ten, removing the stated discrepancy with the dwarf-galaxy rotation-curve measurement.","The measured value $a_{\\rm hNR} = (1.39\\pm0.24)\\times 10^{-6}$ is not excluded; because the stripped-satellite mass distribution extends to zero mass, the observed minimum halo mass $M_h < 5.4\\times10^8\\,M_\\odot$ does not constrain warm dark matter as tightly as claimed.","Future galaxy-based limits on warm dark matter must include non-linear regeneration, halo stripping, and baryon perturbations; Lyman-$\\alpha$ limits must additionally model leftover neutral hydrogen clouds in the reionized universe.","Baryon perturbations alone generate a tail with $n \\approx 0.86$ to 0.96, matching the measured range $0.5 \\lesssim n \\lesssim 1.1$, which the paper says is enough to bring predicted stellar mass and ultraviolet luminosity distributions into agreement with observations down to $M_{\\rm PS} \\approx 5\\times10^8\\,M_\\odot$.","A future direct measurement of the linear power spectrum up to $k \\approx 20$ Mpc$^{-1}$ by cosmic microwave background weak lensing can test whether the inferred tail is real or an artifact of galaxy-based proxies."],"supporting_citations":[{"why":"The source of the limit $a_{\\rm hNR} < 6\\times 10^{-8}$ from the minimum halo mass hosting Milky Way satellites, and the target of the paper's critique.","marker":"[2]"},{"why":"Supplies the warm dark matter simulation suppression factor used to translate the satellite constraint into the revised estimate $a_{\\rm hNR} \\lesssim 6\\times 10^{-7}$.","marker":"[27]"},{"why":"Provides the hydrodynamical equations used for the spherical collapse integrations that form the paper's galaxy models.","marker":"[15]"},{"why":"Gives the cored isothermal sphere solution and the cosmological origin of the core radius, used to identify the pivot point and halo mass.","marker":"[16]"},{"why":"Measures the tail exponent range $0.5 \\lesssim n \\lesssim 1.1$ used to bring predicted stellar mass and UV luminosity distributions into line with observations.","marker":"[25]"},{"why":"Compares predicted and observed galaxy stellar mass and ultraviolet luminosity distributions down to $M_{\\rm PS} \\approx 5\\times10^8\\,M_\\odot$, supporting the claim that no discrepancy remains.","marker":"[26]"},{"why":"Calculates the linear power spectrum tail for bosons decoupling while ultra-relativistic, giving the low-momentum enhancement that behaves as cold dark matter.","marker":"[18]"},{"why":"Defines halo mass $M_h$ through $r_{200}$ and abundance matching, setting the comparison basis for the observed satellite minimum mass.","marker":"[21]"}],"fun_headline_variants":["Missing physics erases warm dark matter discrepancy","Warm dark matter bound falls to halo stripping and baryons","Linear-theory limit on warm dark matter is an artifact","Non-linear effects dissolve warm dark matter tension"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the spherically symmetric hydrodynamical integrations, each initialized with a single Gaussian overdensity at redshift $z_i = 65.9$, faithfully represent how real three-dimensional warm dark matter halos form, get stripped by neighbors, and survive to the present; if those idealized runs miss the relevant dynamics, the claim that baryons and stripping erase the discrepancy between the limit and the measurement does not follow.","fun_headline_variants_meta":{"raw":{"variants":["Missing physics erases warm dark matter discrepancy","Warm dark matter bound falls to halo stripping and baryons","Linear-theory limit on warm dark matter is an artifact","Non-linear effects dissolve warm dark matter tension"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000675,"raw_usage":{"total_tokens":3056,"prompt_tokens":914,"completion_tokens":2142,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":530,"completion_tokens_details":{"reasoning_tokens":2080}},"tokens_in":530,"tokens_out":2142,"duration_ms":14652,"temperature":1.0,"reasoning_tokens":2080,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:06:32.201504+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A full three-dimensional cosmological hydrodynamic simulation of warm dark matter with $m_h = 0.15$ keV and $v_{\\rm hrms}(1) = 493$ m/s, using the standard free-streaming cutoff without an added empirical tail, should produce a population of stripped-down Milky Way satellites with halo masses below $M_h = 5.4\\times 10^8\\,M_\\odot$ and large cores at $z=0$. If instead the minimum surviving satellite halo mass in such a simulation stays above $5.4\\times 10^8\\,M_\\odot$, the paper's reconciliation fails.","supporting_citations":[{"cited_title":"Cosmology of Single Species Hidden Dark Matter","cited_arxiv_id":"2305.08943","evidence_quote":"The source of the limit $a_{\\rm hNR} < 6\\times 10^{-8}$ from the minimum halo mass hosting Milky Way satellites, and the target of the paper's critique."},{"cited_title":"(2023) Understanding the Formation of Galaxies with Warm Dark Matter","cited_arxiv_id":null,"evidence_quote":"Provides the hydrodynamical equations used for the spherical collapse integrations that form the paper's galaxy models."},{"cited_title":"(2025) Why Do Galaxies Have Extended Flat Rotat ion Curves? International Journal of Astronomy and Astrophysics , 15, 1-10","cited_arxiv_id":null,"evidence_quote":"Gives the cored isothermal sphere solution and the cosmological origin of the core radius, used to identify the pivot point and halo mass."},{"cited_title":"(2024) Are James Webb Space Telescope Observ ations Consistent with Warm Dark Matter? International Journal of Astronomy and Astrophysics, 14, 45-60","cited_arxiv_id":null,"evidence_quote":"Compares predicted and observed galaxy stellar mass and ultraviolet luminosity distributions down to $M_{\\rm PS} \\approx 5\\times10^8\\,M_\\odot$, supporting the claim that no discrepancy remains."},{"cited_title":"(2008) The dark m at- ter transfer function: free streaming, particle statistics and me mory of gravitational clustering","cited_arxiv_id":null,"evidence_quote":"Calculates the linear power spectrum tail for bosons decoupling while ultra-relativistic, giving the low-momentum enhancement that behaves as cold dark matter."}],"review_version":1}