{"id":"a31d07ae-826b-4844-80c7-4b0e5e81e292","arxiv_id":"2509.09012","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"The authors compute D, Js, Jp and Jd factors for 20 dwarf spheroidals from GravSphere fits and give power-law scaling relations in distance, half-light radius, and velocity dispersion.","lead":"This paper computes the geometric factors that determine how bright signals from dark matter decay or annihilation would appear from 20 dwarf galaxies, and provides simple formulas that let other researchers estimate these factors from just three galaxy measurements. It is a practical reference for dark matter searches because the factors they compute are exactly what experiments need to translate a null observation into a bound on particle physics models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Jd values and the Jd scaling relation rest on an unverified Maxwellian velocity distribution; Jd depends on the fourth velocity moment, so a non-Maxwellian DF with the same ρ,β could shift Jd beyond quoted errors.","rationale":"The reader's weakest assumption is that the DM velocity distribution is Maxwellian with anisotropy priors. I agree this is the most load-bearing unverified input. The paper's D and Js factors do not depend on the velocity distribution and are solidly grounded; the comparison with Pace & Strigari and the distance-perturbation check for Sculptor support the D/Js results. The Jp factor is also relatively robust because it depends only on second moments, which are determined by the Jeans equation. The Jd factor, however, explicitly requires the fourth velocity moment, and the Maxwellian form enters through Eq. (9). The paper's own consistency checks vary β but never the shape of the velocity distribution, so the conclusion that Jd is insensitive to the velocity distribution is not supported by the presented evidence. The proposed test would settle this directly. The paper is otherwise careful, with correct line-of-sight integrals and honest discussion of the anisotropy uncertainty, so the reader's CONDITIONAL verdict remains appropriate; my concern does not move the verdict, but it sharpens the condition under which the central Jd scaling claim would hold.","tokens_in":12435,"tokens_out":9573,"duration_ms":106659,"concrete_test":"For one representative galaxy (e.g., Draco), take the best-fit ρ(r) and β(r) from Ref. [1] and construct a non-Maxwellian spherical distribution function (e.g., Eddington inversion for isotropic β, or a two-Gaussian mixture with the same σ_r and σ_t) that reproduces exactly those profiles. Recompute Jd(0.5°) from Eq. (6) using the full velocity integral instead of Eq. (9). If log10 Jd shifts by more than the quoted 1σ error bar in Table 1, the Maxwellian assumption is load-bearing and the scaling relation for Jd is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central novelty is the velocity-dependent Jp and Jd factors and their scaling relations. Jp is robust to the assumed Maxwellian shape because <vrel^2> depends only on the second moments, but Jd (Eq. 6) depends on <vrel^4>, which for a general zero-mean distribution is 2<v^4>+4Σσ_i^4+2<v^2>^2. Equation (9) evaluates this using the Maxwellian f(v,r) of Eq. (1). If the true DM velocity distribution is non-Gaussian, as found in N-body halos, the fourth moment differs even for identical ρ(r) and β(r). The paper's robustness check (Table 2, Fig. 2) varies only βDM (isotropic vs stellar-like), not the functional form of the velocity distribution. The conclusion that anisotropic DM velocities shift Jp/Jd 'well within the error bars' therefore does not cover the Maxwellian assumption itself. Because Jd enters the paper's headline scaling relation Jd ∝ σ^8 d^-2 rh^-1 (Table 5), and because this exponent is motivated by a Gaussian closure, the scaling's validity is load-bearing on an unverified distributional assumption. This concern is explicitly acknowledged indirectly in Section 5 ('we can not measure the DM velocity anisotropy'), but the shape assumption goes unquantified against any non-Maxwellian model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the astrophysical D-factor (decay) and Js-, Jp-, Jd-factors (s-, p-, d-wave annihilation) for 20 dwarf spheroidal galaxies using the DM density profiles, velocity dispersions, and velocity-anisotropy parameters obtained in the companion paper by the same authors (Ref. [1]) from GravSphere Jeans fits. The new elements are the velocity-dependent Jp and Jd factors, which include a possible anisotropy in DM velocities, and empirical power-law scaling relations in line-of-sight stellar velocity dispersion, distance, and half-light radius. The authors compare their D and Js results with Pace & Strigari, report agreement within 1--2 sigma, and propose fixed-exponent scalings D ~ sigma^2 d^-2 rh, Js ~ sigma^4 d^-2 rh^-1, Jp ~ sigma^6 d^-2 rh^-1, and Jd ~ sigma^8 d^-2 rh^-1.","tokens_in":12844,"tokens_out":8397,"duration_ms":99496,"significance":"If the results hold, the paper offers a practical way to estimate indirect-detection geometrical factors for dwarf spheroidals without rerunning full Jeans fits, and it is one of the first systematic treatments of p- and d-wave J-factors with anisotropic DM velocities. The line-of-sight integral reduction in Eqs. (7)--(9) is clean and the explicit tables and reproducibility link to the input profiles are strengths. The central limitation is that Jd depends on the assumed Maxwellian form of the DM velocity distribution, and the robustness checks in the paper only vary the DM velocity anisotropy, not the distributional shape. The paper is therefore a useful step, but the Jd scaling relation should be presented as conditional on the Maxwellian closure unless the shape sensitivity is quantified.","major_comments":[{"comment":"The derivation of Jd assumes the DM velocity distribution is exactly Maxwellian, as in Eq. (1). For fixed density and anisotropy, the d-wave integral depends on the fourth velocity moment: <vrel^4> = 2<v^4> + 4Σσ_i^4 + 2<v^2>^2, so a non-Maxwellian DF with the same ρ(r) and β(r) can change Jd. The robustness check in Table 2 varies only βDM (isotropic, prior, stellar-like), not the functional form of f(v,r). Since N-body halos generally have non-Gaussian velocity distributions, the quoted errors on Jd do not include the dominant shape uncertainty. Please quantify this by repeating the Jd calculation with a family of non-Maxwellian velocity distributions matched to the same ρ, β, σ_r (e.g., a q-Gaussian or a two-component Gaussian with modified kurtosis) and report the spread relative to Table 2. This is essential because the headline Jd ∝ sigma^8 scaling in Table 5 is obtained from the G","section":"Section 3, Eq. (9); Section 5, Table 2"},{"comment":"The quoted 1σ error bars for Jp and Jd are Monte Carlo spreads under the fixed flat priors β∞∈[0,0.56], r0∈[0,3.10] kpc, n∈[0.73,1.36], β0=0. These priors are not derived from data. The alternative-anisotropy check in Table 2 is limited to two scenarios, and for four galaxies (Sextans, NGC 6822, Carina, Fornax) the βDM=β* choice changes Jp/Jd by factors of 2.5--3. The conclusion that the effect is 'well within 1σ' is therefore a statement about the chosen conditional errors, not a bound on the systematic uncertainty in the DM velocity anisotropy. The abstract and conclusions should explicitly state that all Jp/Jd results are conditional on the adopted anisotropy priors and Maxwellian shape, or the analysis should explore a wider prior set.","section":"Section 2, Eq. (2); Table 2"},{"comment":"The paper does not report any goodness-of-fit statistic or residual scatter for the scaling relations. For the fixed 0.5° aperture, the free-exponent fits are not fully consistent with the fixed 'geometrical' exponents: Table 3 gives γ_rh = 0.08±0.27 for D versus the expected 1, and Table 5 gives γ_σLOS = 6.66±0.44 for Jp versus 6 and γ_σLOS = 9.06±0.64 for Jd versus 8. The normalization-only, fixed-exponent fits may therefore have substantial intrinsic scatter. The claim that the fixed-exponent relations 'describe' the computed factors requires a reported RMS scatter in log10(factor) or a chi-square/dof for each fit. Without this, the empirical approximation quality cannot be judged, especially for Jp and Jd.","section":"Section 4, Tables 3--5"}],"minor_comments":[{"comment":"The sentence 'the error bars in the J-factors become smaller' appears in the paragraph discussing the D-factor fits in Table 3; it should likely read 'D-factors'.","section":"Section 4, text near Fig. 3"},{"comment":"Typo: 'it's half-radius' should be 'its half-radius'.","section":"Abstract"},{"comment":"Typo: 'power-low form' should be 'power-law form'.","section":"Eq. (10) paragraph"},{"comment":"When stating θmax = αc/2 ≈ rh/d, note that rh/d is in radians; the text omits the conversion to degrees. Please specify the angular units explicitly.","section":"Section 4, θmax definition"},{"comment":"The definition l± = d cosθ ± sqrt((250 kpc)^2 - (d sinθ)^2) assumes a spherical integration boundary of 250 kpc. This is fine, but it would help to state explicitly that the profile is truncated at 250 kpc consistently with GravSphere.","section":"Section 1, line-of-sight integration range"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about its main limitation in Section 5, but the abstract's 'reliable estimates' wording overstates the status of Jd, whose central scaling relation depends on the unverified Maxwellian shape of the DM velocity distribution. The use of Ref. [1] as input is appropriate, and the paper's contribution is a useful application rather than a breakthrough. I would support publication after the authors add a non-Maxwellian sensitivity test and report scatter for the scaling fits."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuine news in this paper is the computation of p- and d-wave annihilation factors (Jp, Jd) for 20 dSphs, with dark-matter velocity anisotropy included. The line-of-sight integrals are correct, and the scaling relations for D, Js, Jp, Jd are clean and physically motivated. The comparison with Pace & Strigari for D and Js shows consistency within 1-2σ, which is exactly what you'd want from a refined analysis. The authors are also honest about the anisotropy prior: they vary β_DM between isotropic and stellar-like and show the effect stays within the 1σ band for most galaxies, with a few outliers that they check by exclusion.\n\nThe soft spot is the Maxwellian velocity distribution. For Jp, the result depends only on second moments, so the shape of the distribution does not matter. For Jd, the fourth moment enters, and the Maxwellian assumption is doing real work. The robustness check varies only the anisotropy parameter, not the functional form. A non-Maxwellian distribution with the same ρ and β could shift Jd beyond the quoted errors. The paper acknowledges the anisotropy problem but does not quantify the shape problem. That is a genuine limitation, though it does not sink the paper. The scaling relations are fits to computed numbers, not independent predictions, but they are benchmarked against previous work and the fits are stable under galaxy exclusion. That is enough to make them useful.\n\nI would send this to peer review. A referee should push for either a non-Maxwellian test or a clear statement that the Jd errors do not include shape uncertainty. I would cite the J-factor tables and the scaling relations in my own work, and I would bring it to a reading group that cares about the limits of dSph mass reconstruction.","headline":"Useful update: Jp/Jd factors with DM velocity anisotropy for 20 dSphs, plus scaling relations; Jd is conditional on a Maxwellian shape assumption that the paper does not test.","tokens_in":13382,"tokens_out":3365,"would_cite":true,"duration_ms":33807,"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":"Dark matter decay and annihilation rates in dwarf galaxies follow a simple power law in three observables.","keywords":["dark matter","dwarf spheroidal galaxies","D-factor","J-factor","indirect detection","velocity anisotropy","scaling relations","Jeans equation"],"falsifier":"Compare the Maxwellian-based J_p/J_d prediction with the velocity distribution in a realistic N-body simulation of a dwarf-scale halo: if the simulated ⟨v_rel²⟩ differs from 2σ_r²(3−2β) by more than the quoted uncertainties, the velocity-integration step fails. Alternatively, a future dwarf galaxy whose scaling-relation prediction for D or J_s is off by more than 3σ from the Jeans-model integral would falsify the claim that three observables suffice.","tokens_in":12287,"feed_emoji":"🌌","tokens_out":6012,"duration_ms":58448,"temperature":0.7,"pith_summary":"The paper takes dark matter density and velocity profiles for 20 dwarf spheroidal galaxies, obtained by fitting stellar kinematics with a Jeans-equation solver, and integrates them to obtain the geometric factors that control dark matter decay and annihilation signals. It claims that these factors are accurately described by power-law scaling relations in distance, half-light radius, and line-of-sight stellar velocity dispersion, with exponents fixed by elementary Newtonian and dimensional arguments: D ~ σ² d⁻² r_h, J_s ~ σ⁴ d⁻² r_h⁻¹, J_p ~ σ⁶ d⁻² r_h⁻¹, J_d ~ σ⁸ d⁻² r_h⁻¹. The new ingredient is allowing dark matter velocities to be anisotropic, which affects the p-wave and d-wave annihilation factors in a way that remains within the 1σ uncertainties of the isotropic calculation for most galaxies. The practical payoff is that observers can estimate expected dark matter signals for any dwarf galaxy from three easily measured quantities rather than from a full dynamical model.","feed_headline":"Dwarf galaxy dark matter signals follow a 3-variable power law","feed_subtitle":"New fits to 20 dwarf galaxies turn three easy measurements into decay and annihilation signal estimates.","key_machinery":"The central objects are the line-of-sight integrals D = ∫dΩ dl ρ(r), J_s = ∫dΩ dl ρ²(r), and the velocity-weighted generalizations J_p = ∫dΩ dl ρ²(r)⟨v_rel²⟩, J_d = ∫dΩ dl ρ²(r)⟨v_rel⁴⟩. Assuming a Maxwellian velocity distribution with a radially dependent anisotropy β(r), the velocity averages close analytically: J_p = ∫dΩ dl ρ² 2σ_r²(3−2β) and J_d = ∫dΩ dl ρ² 4σ_r⁴(15−20β+8β²). These formulas turn an astrophysical measurement of ρ(r), σ_r(r), β(r) into particle-physics-relevant factors, and the paper evaluates them for 20 dSphs using Jeans-equation fits to stellar observations.","core_discovery":"The central claim is that for a set of 20 dwarf spheroidal galaxies, the phase-space distribution of dark matter reconstructed from stellar kinematics determines the decay D-factor and annihilation J_s-, J_p-, J_d-factors, and that these integrals are captured by empirical scaling relations. When the angular integration is matched to the galaxy half-light radius, the exponents follow the dimensional pattern M(r_h) ∝ σ² r_h and ρ ∝ σ²/r_h², giving log10 D = 16.63 ± 0.01 (σ/5 km/s)² (d/100 kpc)⁻² (r_h/100 pc), log10 J_s = 17.92 ± 0.04 (σ/5 km/s)⁴ (d/100 kpc)⁻² (r_h/100 pc)⁻¹, with analogous expressions J_p = 15.59 ± 0.08 and J_d = 13.68 ± 0.11 with σ powers 6 and 8. The paper also shows that t","pith_inferences":["An untested extension is applying the same scaling relations to ultra-faint dwarfs or dwarf irregulars with sparse data; the functional form may hold, but the fitted normalizations could shift for systems outside the original sample.","The anisotropy dependence of J_p and J_d suggests a way to break the degeneracy: combining stellar proper-motion measurements with realistic simulations of dark matter velocity distributions could replace the flat priors on β∞, r0, n with physically motivated distributions.","Because the scaling relations use only observables, they could be used as a fast Monte Carlo prior in joint analyses of many dwarf galaxies, turning the D/J factors into analytic functions of data with propagated uncertainties."],"forward_implications":["Any dwarf galaxy with measured distance, half-light radius, and line-of-sight velocity dispersion can be assigned D, J_s, J_p, J_d factors through the scaling relations, without redoing the full dynamical analysis.","Bounds on dark matter decay and annihilation from gamma-ray, neutrino, and cosmic-ray observations of dwarf galaxies can be recast using these updated factors; for decay and s-wave annihilation the change relative to previous scaling fits is within 1–2σ.","The first scaling relations for p-wave and d-wave annihilation factors make it possible to constrain models where s-wave annihilation is helicity- or symmetry-suppressed.","For the majority of the 20 galaxies, ignorance of dark matter velocity anisotropy does not dramatically alter J_p and J_d; the main exceptions are Sextans, NGC 6822, Carina, and Fornax, where the choice of anisotropy changes individual factors by a factor 2.5–3."],"fun_headline_variants":["Dark matter signals scale with distance, size, and velocity","Simple scaling laws predict dark matter decay and annihilation","20 dwarf galaxies yield universal dark matter signal formulas","Refined factors for dark matter signals from stellar kinematics","Three measurements estimate dark matter decay and annihilation"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The dark matter velocity distribution is assumed to be Maxwellian with an anisotropy profile of fixed functional form; this velocity distribution is not directly observable, and the derived J_p and J_d factors depend on it.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter signals scale with distance, size, and velocity","Simple scaling laws predict dark matter decay and annihilation","20 dwarf galaxies yield universal dark matter signal formulas","Refined factors for dark matter signals from stellar kinematics","Three measurements estimate dark matter decay and annihilation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000211,"raw_usage":{"total_tokens":1318,"prompt_tokens":875,"completion_tokens":443,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":619,"completion_tokens_details":{"reasoning_tokens":370}},"tokens_in":619,"tokens_out":443,"duration_ms":5247,"temperature":1.0,"reasoning_tokens":370,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T19:48:22.612309+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the Maxwellian-based J_p/J_d prediction with the velocity distribution in a realistic N-body simulation of a dwarf-scale halo: if the simulated ⟨v_rel²⟩ differs from 2σ_r²(3−2β) by more than the quoted uncertainties, the velocity-integration step fails. Alternatively, a future dwarf galaxy whose scaling-relation prediction for D or J_s is off by more than 3σ from the Jeans-model integral would falsify the claim that three observables suffice.","supporting_citations":[],"review_version":1}