{"id":"fec3bd89-9640-4568-9f83-500671dfc7e7","arxiv_id":"1908.06067","paper_version":3,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Dark matter self-interaction cross sections can be accurately parametrized by a two-parameter effective-range formula across many proposed particle models.","lead":"The paper proposes that dark matter self-interactions can be described by just two numbers, the scattering length and effective range, regardless of the underlying particle model. This simple description captures the velocity dependence needed to explain small-scale structure while respecting cluster constraints.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the two-parameter effective-range parametrization is accurate within its disclosed low-energy, short-range domain, and the classical-regime failure is explicitly acknowledged.","rationale":"The reader's weakest assumption identifies the effective-range truncation and the classical-regime failure. I agree this is the most fragile aspect of the argument, but the paper explicitly discloses this limitation and offers extensions, so it does not undermine the appropriately-scoped central claim. The independent numerical validation for Yukawa potentials in the resonant and Born regimes, plus the physically correct deuteron relation in Eq. (18), provide real support for Eq. (8) within its validity domain. The only concrete technical issue found is the typo in Eq. (17), which is secondary and does not affect the main parametrization. Therefore the reader's ACCEPT verdict should stand unchanged.","tokens_in":21002,"tokens_out":21788,"duration_ms":208476,"concrete_test":"Verify Eq. (8) against an independent numerical Schrödinger solver for a Yukawa potential at representative points in Figs. 8 and 9 (e.g., αm/mφ = 1.68 resonance and αm/mφ = 10) over velocities v/α = 0.02-0.2; also verify that for mφ = 3 MeV and v = 1000 km/s the full multi-partial-wave cross section deviates from Eq. (8) as claimed, and re-derive Eq. (17) from the pole condition to confirm the typo.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the standard effective-range expansion k cot δ0 = -1/a + (1/2) r_e k^2, whose validity requires kR << 1 and S-wave dominance. The paper's own Fig. 1 and Sec. II identify the classical regime (mφ ≲ 5 MeV for the benchmark Yukawa parameters) where this condition fails, and Sec. V provides extensions for antiresonances and sharp resonances. Because these limitations are stated rather than hidden, they do not falsify the central claim; they delimit it. The numerical comparisons in Figs. 8 and 9 support the approximation in the resonant and Born regimes where the stated assumptions hold. One minor internal inconsistency is that Eq. (17) as printed gives the S-matrix pole as k = (i a/2)(1 ± sqrt(1 - 2 r_e/a)), whereas the pole condition k cot δ = i k derived from Eq. (4) yields k = (i ± sqrt(2 r_e/a - 1))/r_e. This does not affect Eq. (8), and the deuteron binding-energy check in Eq. (18) uses the correct relation, so it appears to be a typographical error rather than a load-bearing flaw.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes the effective-range expansion (ERE) as a practical, model-independent parametrization of the velocity-dependent S-wave dark matter self-interaction cross section. The central formula, Eq. (8), expresses σ(v) in terms of the scattering length a, the effective range r_e, and the DM mass m, and the authors benchmark this approximation against exact numerical Schrödinger solutions for attractive and repulsive Yukawa potentials (Figs. 1, 8, 9), finding excellent agreement outside the classical regime. They then use the parametrization to derive astrophysical constraints (Figs. 4 and 5), fit the parameters to inferred halo cross sections (Fig. 6), interpret a and r_e in terms of concrete models (contact interaction, light mediator, bound/virtual states, resonances, SIMPs), and propose extensions for antiresonances, sharp resonances, and inelastic scattering (Sec. V). The paper is explicit about the domain of validity of the ERE, identifying the classical regime and antiresonances as places where the two-parameter formula fails.","tokens_in":21257,"tokens_out":12431,"duration_ms":104354,"significance":"If the claims hold, the paper provides a valuable model-independent tool for SIDM simulations and phenomenological studies, analogous to the role of the effective-range expansion in nuclear physics. Its strengths include: (i) the central formula is derived from the standard analytic expansion of k cot δ, not fitted to numerical data; (ii) the approximation is validated against exact quantum-mechanical results for both signs of the Yukawa potential; (iii) the limitations are disclosed rather than hidden, with concrete extensions in Sec. V; (iv) the insight that only |a| and r_e/a can be constrained by velocity-dependent cross sections is stated and demonstrated. The astrophysical fits in Fig. 6 are appropriately hedged as illustrative. The paper does not oversell the parametrization; it explicitly notes where it fails.","major_comments":[],"minor_comments":[{"comment":"Equation (17) as printed is dimensionally inconsistent: the left-hand side k has units of inverse length, while the right-hand side (i a/2)(1 ± sqrt(1 - 2 r_e/a)) has units of length. The correct solution of the pole condition k cot δ = i k is k = (i/r_e)(1 ± sqrt(1 - 2 r_e/a)), and the subsequent deuteron binding-energy check in Eq. (18) and Fig. 10 are consistent with the correct formula. Please correct Eq. (17) and any related text.","section":"Sec. IV.C, Eq. (17)"},{"comment":"The paper clearly acknowledges the classical regime (m_φ ≲ 5 MeV for the benchmark) where the effective-range approximation fails, but since the abstract and introduction emphasize the broad applicability of the two-parameter formula, consider adding an explicit caveat in the abstract that the parametrization applies to short-range interactions with kR ≪ 1.","section":"Sec. II, Fig. 1"},{"comment":"The fitting procedure that produces the benchmark parameter sets S1–S4 is not fully described; please state the likelihood or χ² definition and the velocity distribution parameters (e.g., v0 for each halo class) so that the fit is reproducible.","section":"Sec. III.B, Fig. 6"},{"comment":"In the Born regime, Eq. (15) gives r_e = 4/(m α), which is negative for a repulsive Yukawa potential (α < 0). The text does not comment on this sign, and since Eq. (8) depends only on r_e/a, negative values are physically acceptable; a short remark would avoid confusion.","section":"Sec. IV.B, Eq. (15)"},{"comment":"In the improved antiresonance formula, the parameters a and r_e are not the standard scattering length and effective range, as the text notes; however, using different symbols (e.g., ã and r̃_e) would prevent readers from confusing them with the parameters of Sec. II.","section":"Sec. V.A, Eq. (21)"}],"recommendation":"minor_revision","confidential_remarks":"The paper is well within scope and the central claim is sound. The only substantive issue is the typographical error in Eq. (17), which should be corrected before publication; the remaining points are presentation and clarity improvements. I concur with the reader's positive assessment and believe a minor revision suffices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful, honest paper that makes the effective-range expansion the default two-parameter language for SIDM velocity dependence. It doesn't invent new physics, but it packages textbook scattering theory in a way that should be practical for simulations and phenomenology.\n\nThe genuinely new piece is the framing: earlier papers used effective range in specific models (refs [23-27]), but this one treats it as a model-independent parametrization, shows it interpolates Yukawa, resonant, bound-state and SIMP scenarios, and adds a unitarity-preserving improvement for antiresonances and sharp resonances. That last bit, Eq. (21), is a real addition—the naive hard-sphere-like term violates unitarity, and their factorization into e^{2ikR} times the standard amplitude fixes it.\n\nThe paper earns its keep on the numerics. They benchmark Eq. (8) against exact Schrödinger solutions for attractive and repulsive Yukawa potentials across the Born, resonant, and antiresonance regimes (Figs. 8, 9). The agreement is good in the regime where the assumptions hold, and they are explicit about the failure in the classical long-range regime, where kR ~ 1 and higher partial waves matter. No circular fitting here: the two parameters are extracted from the shape of the exact phase shift, not fitted to the data they claim to predict. The astrophysical fits in Fig. 6 are clearly labeled as illustrative, with proper caveats about tidal stripping.\n\nSoft spots: first, the applicability domain is genuinely narrow—for light mediators of a few MeV, the Yukawa range is comparable to the de Broglie wavelength at cluster velocities, and the two-parameter formula breaks down. That is not hidden; it is stated in Sec. II, but it does limit the \"model-independent\" claim to short-range interactions. Second, there is a minor internal inconsistency in the printed Eq. (17): the pole condition derived from Eq. (4) gives k = (i ± sqrt(2 r_e/a - 1))/r_e, not the expression with a/2 and sqrt(1 - 2r_e/a). Since Eq. (18) uses the correct relation and the main cross section formula is unaffected, this looks like a typo, not a load-bearing error. Third, the constraint that sub-GeV DM is excluded relies on the cluster upper bounds, which are still debated; they acknowledge this.\n\nThe paper is aimed at people doing SIDM simulations or translating bounds into (a, r_e). For that audience it is a solid reference. I'd send it to a competent referee—the typos and the pole formula should be fixed, but the central argument holds up.","headline":"A practical, honest paper that makes the effective-range two-parameter language the default for SIDM velocity dependence; the numerics check out and the limitations are disclosed, so it deserves a real referee.","tokens_in":21794,"tokens_out":2195,"would_cite":true,"duration_ms":20905,"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":"Dark matter's self-scattering reduces to a two-parameter formula.","keywords":["effective range theory","dark matter self-interactions","scattering length","effective range","velocity-dependent cross section","Yukawa potential","resonant dark matter","SIDM"],"falsifier":"Compute the exact S-wave cross section for an attractive Yukawa potential with $m_\\phi=1\\,\\mathrm{MeV}$, $m=3\\,\\mathrm{GeV}$, and $\\alpha=0.04$ at $v=1000\\,\\mathrm{km/s}$, and compare it with Eq. (8) using $a$ and $r_e$ extracted from low velocities; agreement would support the universality claim, while the mismatch in this classical regime would show where the two-parameter formula stops working.","tokens_in":20802,"feed_emoji":"🌌","tokens_out":11522,"duration_ms":93145,"temperature":0.7,"pith_summary":"The paper claims that the velocity-dependent cross section for dark matter self-interactions in halos can be described, without committing to a specific particle model, by two parameters: the S-wave scattering length $a$ and the effective range $r_e$, together with the dark matter mass $m$. This matters because dwarf galaxies seem to require a self-interaction cross section of order $1$--$10\\,\\mathrm{cm^2/g}$ at velocities around $10\\,\\mathrm{km/s}$, while galaxy clusters bound it to below about $1\\,\\mathrm{cm^2/g}$ at velocities near $1000\\,\\mathrm{km/s}$, and translating this velocity dependence into model parameters has usually been done case by case. The paper shows that the same two-parameter formula reproduces the full numerical cross sections for Yukawa forces, for Breit-Wigner resonances, and for scattering through bound states or virtual levels. If right, this gives astrophysical simulations and model comparisons a common, minimal description of self-interacting dark matter.","feed_headline":"Two parameters capture dark matter self-scattering across halo scales","feed_subtitle":"One formula with just two numbers links dwarf galaxies, galaxy clusters, resonances, and bound states.","key_machinery":"The load-bearing object is the effective-range expansion of the phase shift,\n$$\n$k^{{2\\ell+1}}$\\cot\\delta_\\ell(k) \\simeq -\\frac{1}{a_\\$ell^{{2\\ell+1}}$} + \\frac{1}{2} r_{e,\\ell}^{2\\ell-1} $k^{2}$,\n$$\ntruncated at order $k^2$; for S-wave scattering this produces Eq. (4), and written in terms of the relative velocity it becomes Eq. (8). The expansion is justified because $k^{2\\ell+1}\\cot\\delta_\\ell(k)$ is analytic at $k=0$ for finite-range interactions, so low-energy scattering depends only on the scattering length and effective range. The paper also derives a first-order differential equation for the phase shift, which makes the numerical extraction of $a$ and $r_e$ straightforward, and uses the amplitude poles $k_\\pm = (i/a)(1 \\pm \\sqrt{1-2r_e/a})$ to connect those two parameters to bound states, virtual levels, and resonances. For cases where the minimal formula fails, a version with a background phase $e^{2ikR}$ describes antiresonances and sharp resonances.","core_discovery":"The central claim is that the effective-range approximation for non-relativistic scattering works for dark matter self-interactions across the astrophysically relevant velocity range. Concretely, the S-wave cross section is $$\n\\$\\sigma$(v) = 4\\pi $a^{2}$ \\left( \\left(1 - \\frac{1}{8}\\frac{r_e}{a}(m a v)^2\\right)^2 + \\frac{1}{4}(m a v)^2 \\right)^{-1},\n$$ where $v$ is the relative velocity, $m$ the dark matter mass, $a$ the scattering length, and $r_e$ the effective range; only $ma$ and the ratio $r_e/a$ control the velocity dependence. The paper verifies this formula against numerical solutions of the Schrödinger equation for attractive and repulsive Yukawa potentials, shows that it reproduces a Breit-Wigner resonance with an energy-dependent width, and traces its poles in the complex momentum plane to bound states, virtual levels, or resonances. Fitting the formula to semi-analytically extracted halo cross sections gives benchmark parameters such as $a=19.2\\,\\mathrm{fm}$, $r_e=0.01\\,\\mathrm{fm}$, and $m=14.9\\,\\mathrm{GeV}$ for the best-fit case.","pith_inferences":["A measurement of $\\sigma$ at two well-separated velocities (say dwarf and cluster scales) would fix $a$ and $r_e$; a third measurement would test whether nature's self-interactions are truly effective-range-like or require the improved antiresonance formula of Section V.","If future data show a cross section rising toward cluster scales, that would select $r_e/a > 1$ and hence a narrow-resonance interpretation; a flat cross section would point to contact-like models—a diagnostic that does not require building a UV-complete theory.","The classical-regime breakdown for light mediators implies that any model with mediator mass of a few MeV must be treated with the full partial-wave sum once cluster-scale velocities matter; the two-parameter formula is not a universal fit for such models.","The benchmark fits suggest that current data do not strongly constrain $r_e/a$, so upcoming simulations exploring sub-halo dynamics inside the Milky Way may be the first place where the effective range actually becomes measurable."],"forward_implications":["Astrophysical simulations can implement self-interacting dark matter by specifying $(a, r_e, m)$ instead of a full particle model, making different scenarios directly comparable in the same simulation.","For $\\sigma/m \\sim 10\\,\\mathrm{cm^2/g}$ at dwarf scales, dark matter masses below a few GeV are excluded by cluster observations unless the effective-range parameters are tuned to suppress high-velocity scattering.","The velocity dependence of the cross section determines only the ratio $r_e/a$ and the product $ma$; the signs of $a$ and $r_e$ cannot be fixed by $\\sigma(v)$, so the underlying intermediate state (bound, virtual, or resonant) is not observable from the scattering data alone.","The formalism extends to subleading inelastic processes: a complex scattering length relates annihilation and elastic cross sections by $\\sigma_{\\mathrm{an},0}(k) \\simeq \\sigma_0(k) |\\mathrm{Im}\\,a|/(\\mathrm{Re}\\,a)^2$."],"supporting_citations":[{"why":"Introduces the self-interacting dark matter hypothesis that the paper's parametrization is designed to serve.","marker":"[2]"},{"why":"Provides the semi-analytic method and the inferred $\\langle\\sigma v\\rangle/m$ values for clusters, LSB galaxies, and dwarfs that the paper fits with its benchmark parameters.","marker":"[12]"},{"why":"Gives the resonant SIDM scenario whose energy-dependent width $\\Gamma(E)\\propto E^{(2\\ell+1)/2}$ the effective-range formula reproduces.","marker":"[20]"},{"why":"Supplies the proof that low-energy phase shifts depend only on the scattering length and effective range, the theoretical basis of Eq. (8).","marker":"[28]"},{"why":"Supplies the calculations establishing the two-parameter effective-range description of nucleon-nucleon scattering.","marker":"[29]"},{"why":"Provides the light-mediator SIDM formalism and the Hulthén-potential approximation used to compute Yukawa scattering lengths and effective ranges.","marker":"[35]"},{"why":"Textbook reference for amplitude poles, virtual levels, and the background-phase improvement used in Section V.","marker":"[48]"}],"fun_headline_variants":["Two parameters capture all dark matter self-scattering velocity dependence","Effective-range formalism makes dark matter self-interactions model-independent","Two numbers link dwarf galaxies to clusters via dark matter scattering","One formula, two parameters: dark matter self-interactions for all velocities","Scattering length and effective range: dark matter's velocity dependence explained"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The formula assumes the dark-matter force has a short range, so that at every velocity the collision looks point-like and S-wave scattering dominates; for a light mediator whose range rivals the particle's de Broglie wavelength, the approximation breaks down.","fun_headline_variants_meta":{"raw":{"variants":["Two parameters capture all dark matter self-scattering velocity dependence","Effective-range formalism makes dark matter self-interactions model-independent","Two numbers link dwarf galaxies to clusters via dark matter scattering","One formula, two parameters: dark matter self-interactions for all velocities","Scattering length and effective range: dark matter's velocity dependence explained"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000799,"raw_usage":{"total_tokens":3511,"prompt_tokens":935,"completion_tokens":2576,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":2490}},"tokens_in":551,"tokens_out":2576,"duration_ms":17369,"temperature":1.0,"reasoning_tokens":2490,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:57:19.144730+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact S-wave cross section for an attractive Yukawa potential with $m_\\phi=1\\,\\mathrm{MeV}$, $m=3\\,\\mathrm{GeV}$, and $\\alpha=0.04$ at $v=1000\\,\\mathrm{km/s}$, and compare it with Eq. (8) using $a$ and $r_e$ extracted from low velocities; agreement would support the universality claim, while the mismatch in this classical regime would show where the two-parameter formula stops working.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the light-mediator SIDM formalism and the Hulthén-potential approximation used to compute Yukawa scattering lengths and effective ranges."}],"review_version":1}