{"id":"fa72d7ba-e630-480c-8df4-fc62e2898db7","arxiv_id":"2506.12642","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Number-changing dark matter self-interactions can significantly deplete dark matter spikes around supermassive black holes, while the 2 to 1 semi-annihilation generally preserves the spike structure.","lead":"Dark matter that can change the number of its own particles could destroy the dense dark matter spikes around supermassive black holes. The authors calculate when this happens and show it changes the expected signals from these structures.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The all-escape assumption for n→m final products in Eq. (3.10) likely overstates dissolution; recapture by 2→2 scattering raises effective plateau densities, weakening the claim of significant spike depletion for n≥3.","rationale":"The paper's principal novelty is the systematic treatment of n→m processes in DM spikes, and the headline claim is that for n≥3, freeze-out-favored cross-sections significantly deplete the spike. The depletion is driven by the dissolution term in Eq. (3.10), whose plateau solution Eq. (3.12) is presented as a universal density cap. The physical justification of that cap is the assertion that all final-state DM particles are relativistic and escape. If that assertion is wrong, the cap is too low. The paper itself contains the ingredients to see the problem: ξ(r) in Eq. (3.8) is defined as the probability that a final product interacts before radius r, and for the large σ2→2/mχ values in benchmarks D and H, ξ≈1 inside the dissolution radius. A scattered boosted particle does not necessarily escape; it can lose enough energy to become bound, so the net loss per event is n−m rather than n. This is not a disagreement with community consensus—it is an internal tension between the dissolution treatment and the self-heating treatment. The paper's single sentence about the capture rate being small refers to the boosted flux, not to the density profile, so the tension is unresolved. The internal inconsistency with benchmarks E–H (where 3→2 is subdominant) makes the situation worse: the abstract's claim is not demonstrated for the displayed representative parameters. I therefore concur with the reader's CONDITIONAL verdict. The proposed optical-depth test directly checks whether the all-escape assumption holds in the parameter region where the depletion claim is made. No ad hominem is intended; the analytic framework is a useful step, but the loss term must be coupled to capture before the headline claim can be accepted without qualification.","tokens_in":14774,"tokens_out":13578,"duration_ms":160901,"concrete_test":"Compute the scattering optical depth τ(r)=∫_r^{r_esc} dr' [ρ_b(r')/mχ] σ2→2 for a boosted final particle produced at the dissolution radius R_diss, using the density profile from Sec. 3.1–3.2 and the benchmark parameters E–H and the freeze-out benchmark mχ=10 MeV, α3→2=1. If τ≫1, the all-escape assumption of Sec. 3.3 fails; rerun Eq. (3.10) with the net loss n−m and recompute the J-factors in Fig. 3. If the resulting J3 differs by more than 10% from the paper's values, the abstract's depletion claim is overestimated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.3 (Eq. 3.10) removes all n initial particles per n→m event, justified by the statement in Sec. 3.3 that the m final particles are 'boosted to relativistic velocities' and 'able to escape from their orbit at r.' This is the load-bearing mechanism for the abstract's n≥3 depletion claim, since the dissolution plateau (Eq. 3.12) is what caps the spike density. But in the strongly self-interacting regime considered (σ2→2/mχ up to 5×10^-2 cm^2/g), the optical depth for a boosted final particle to scatter before reaching the escape radius is large. The self-heating section already defines ξ(r) (Eq. 3.8) as the probability that a final particle interacts before radius r; for large σ2→2/mχ, ξ≈1 throughout the region where dissolution matters. If a final particle scatters and is recaptured, the net number loss is n−m, not n: for 3→2 this raises the plateau density of Eq. (3.14) by sqrt(3); for 4→2, by 2^{1/3}. The paper's statement that 'the capture rate is much smaller than 1' concerns the boosted flux, not the density profile, so it does not resolve this. In addition, the benchmarks E–H (mχ=100 MeV) already show 3→2 subdominant, so the claim relies on lower masses; the recapture correction is largest exactly where σ2→2 is large enough to matter.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies how number-changing dark-matter self-interactions (n→m processes) affect the density profile of dark-matter spikes around supermassive black holes. It combines three effects: isothermal core formation from 2→2 scattering, self-heating from boosted final-state particles, and density dissolution from 2→0 and n→m processes. The authors derive analytic plateau densities, present benchmark profiles for 2→1 and 3→2 processes at mχ=100 MeV, and compute generalized J-factors. The central claim is that for n≥3 processes with cross-sections favored by thermal freeze-out the spike is significantly depleted, while semi-annihilation 2→1 largely preserves the spike.","tokens_in":15114,"tokens_out":6661,"duration_ms":75898,"significance":"The framework is timely and useful: if the depletion is real, it caps the central spike density and suppresses J-factors for annihilation and boosted-DM searches. The rate equations and plateau-density formulas are straightforward and correctly derived for the stated model, and the J-factor comparison with the NFW baseline is a clear phenomenological output. No parameters are fitted to the target result; the benchmark cross-sections are taken from prior SIMP studies. However, the manuscript's central claim is not currently supported by its own benchmarks and rests on an unquantified all-escape assumption for n→m final states.","major_comments":[{"comment":"Equation (3.10) removes n particles per n→m event, justified by the sentence following it that final-state particles are relativistic and escape the orbit. However, in the strongly self-interacting regime adopted in benchmarks E-H (σ2→2/mχ up to 5×10^-2 cm^2/g), the optical depth ξ(r) defined in Eq. (3.8) is of order unity throughout the region where dissolution is relevant, so most final-state particles scatter before escaping. If a scattered particle is recaptured, the net number loss is n−m rather than n; for 3→2 this raises the plateau density in Eq. (3.14) by sqrt(3), and for 4→2 it raises Eq. (3.15) by 2^(1/3). The statement in Sec. 3.2 that 'the capture rate is much smaller than 1' concerns the boosted flux observed at Earth, not the density profile, so it does not resolve the inconsistency. Please quantify the recapture probability self-consistently and either modify Eq. (3.10) or justify the all-escape limit.","section":"Sec. 3.3, Eq. (3.10)"},{"comment":"The abstract claims 'for n≥3, the spike is significantly depleted for n→m cross-sections favored by DM production via thermal freeze-out.' However, Fig. 2 and the text following it state that for the benchmark 3→2 rates (E-H), neither self-heating nor dissolution significantly alters the density profile, and that for general n≥3 the cross-sections required to make these effects dominant are no longer perturbative. This is an internal contradiction between the headline claim and the paper's own benchmarks. Please revise the abstract and conclusions to match the benchmark results, or demonstrate explicitly that the benchmarks are not representative of the freeze-out-favored parameter region.","section":"Abstract vs. Sec. 3.3 and Fig. 2"},{"comment":"The manuscript states that 'We have checked the capture rate is much smaller than 1 in our interesting parameter region,' but no calculation, equation, or figure is provided for this check. Because this statement is used to separate the effect on the boosted flux from the effect on the density profile, and because the density-profile effect is the load-bearing part of the n≥3 depletion claim, the check should be written out explicitly or replaced by a proper treatment of recapture in the dissolution equation.","section":"Sec. 3.2, self-heating discussion"}],"minor_comments":[{"comment":"The header for the 2→1 column should read [cm^3 s^-1] rather than [cm^3 s]; the current notation is dimensionally inconsistent.","section":"Table 1"},{"comment":"The labels R2→1 and R3→2 appearing in the figures are not defined in the captions; please define them as the self-heating core radii for the respective processes.","section":"Figs. 1 and 2 captions"},{"comment":"The sentence 'We have verified that this behavior remains for general n→m processes with n≥3' is an unsupported assertion in the text; if it is to be retained, please provide the underlying calculation or a supplementary figure showing a representative 4→2 case.","section":"Sec. 3.3, final paragraph"},{"comment":"The heat time-scale estimate uses 'the typical radius' without a precise specification; please state exactly how that radius is chosen from the isothermal/self-heating profile.","section":"Sec. 3.2, Eq. (3.9)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper builds a clean framework for n→m number-changing processes in DM spikes — dissolution plateaus, self-heating, isothermal cores — and the rate equations and plateau densities (Eqs. 3.10–3.12) are correct for the stated model. The genuinely new piece is applying 3→2 and 4→2 processes to spikes and computing the generalized J-factors; the 2→1 semi-annihilation result (the spike structure survives) is a clean, useful finding. The calculations are self-contained and the citation pattern is honest, with prior work by Chu-Garcia-Cely, Kamada-Kim, Shapiro-Paschalidis, and Hochberg et al. properly credited.\n\nThe first soft spot is the mismatch between the abstract and the benchmarks. The abstract says n≥3 processes significantly deplete the spike at freeze-out-favored cross-sections. The E–H benchmarks at mχ=100 MeV show the opposite, and Sec. 3.3 goes further, claiming the subdominance holds for all n≥3 with perturbative couplings. That is not the same statement as the abstract, and the 'verified for general n≥3' claim is not backed by an actual argument in the text. The depletion regime actually lives at the canonical reference point, mχ ~ 10 MeV with α3→2 ~ 1, where ρ_pl ≈ 1.9×10^13 GeV/cm³ cuts below the spike and would suppress J_3 substantially. The authors should add that regime to the benchmark coverage or narrow the abstract; as written, the paper argues against its own headline.\n\nSecond, the all-escape assumption in Eq. (3.10). Counting all n initial particles as lost because the m final particles are relativistic is too strong when σ2→2 is large: the boosted particles can scatter, deposit energy, and be recaptured before escaping, making the net loss closer to n−m. The stress-test version of this overstates the correction — capture plausibly needs many scatterings, not one, and the plateau density shifts by only √3 or 2^{1/3} — but the tension is real. Self-heating assumes these same final particles deposit energy in the spike, while dissolution assumes they leave, and the paper never reconciles the two. A paragraph of honest discussion would fix it.\n\nWho this is for: people computing indirect-detection signals or boosted-DM fluxes from galactic centers, and SIMP phenomenologists. The J-factor formulas are a useful reference even where the depletion claim is contested. Send it to review — it deserves a serious referee, but expect a required revision that aligns the abstract with Sec. 3.3 and the benchmark coverage.","headline":"Competent spike/SIMP framework with correct rate equations, but the abstract's n≥3 depletion claim contradicts the paper's own benchmarks and needs reconciling.","tokens_in":15697,"tokens_out":16519,"would_cite":true,"duration_ms":197287,"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":"Number-changing dark matter self-interactions that convert three or more particles into fewer ones can deplete the dense spikes around supermassive black holes, capping the central density far below the standard power-law profile and…","keywords":["dark matter spike","self-interacting dark matter","number-changing processes","3 to 2 annihilation","semi-annihilation","SIMP dark matter","J-factor","supermassive black hole"],"falsifier":"Compute the optical depth for a boosted final-state particle from an $n\\to m$ reaction to undergo a $2\\to 2$ scattering before it leaves the spike; if a non-negligible fraction is recaptured for the benchmark cross-sections, the dissolution equation overcounts particle loss and the predicted depletion and $J$-factors would need to be revised.","tokens_in":14541,"feed_emoji":"🕳️","tokens_out":6594,"duration_ms":71009,"temperature":0.7,"pith_summary":"The paper asks whether dark matter self-interactions that change particle number can be tested in the ultra-dense dark matter spikes expected around supermassive black holes. It argues that processes converting $n$ initial dark matter particles into $m$ final ones ($n>m$), which are generic in self-interacting dark matter models, can proceed at observable rates there. For $n\\geq 3$, cross-sections favored by thermal freeze-out production of the dark matter significantly deplete the spike, while the $2\\to 1$ semi-annihilation process generally preserves its structure. These density changes affect the $J$-factors that control predicted signals, so the spike's particle physics must be included in phenomenological forecasts.","feed_headline":"3-to-2 dark matter reactions can deplete galaxy-center spikes","feed_subtitle":"Freeze-out-sized n-to-m cross-sections cap spike density far below the standard profile, reshaping signal forecasts.","key_machinery":"The central objects are the plateau density $\\rho_{\\mathrm{pl}}$ from the dissolution equation (the density below which a number-changing process cannot deplete the spike within the halo age) and the three radius scales $R_c$ (isothermal core), $R_{n\\to m}$ (self-heating core), and $R_{\\mathrm{diss}}$ (dissolution radius) that partition the spike profile. The argument operates by evolving the local dark matter number density $n_\\chi(r,t)$ with $\\dot{n}_\\chi = -\\langle \\sigma_{2\\to 0}v\\rangle n_\\chi^2 - (n/n!)\\langle \\sigma_{n\\to m}v^{n-1}\\rangle n_\\chi^n$, assuming all $n$ initial particles are lost because the $m$ final-state particles are relativistic and escape. The radius where the resulting density saturates at $\\rho_{\\mathrm{pl}}$ gives the depletion boundary; comparing $R_{\\mathrm{diss}}$, $R_c$, and $R_{n\\to m}$ determines which effect dominates the observed profile. Cross-sections are parametrized as $\\sigma_{n\\to m}v^{n-1} \\equiv \\alpha_{n\\to m}^n / m_\\chi^{3n-4}$, following the freeze-out literature, so the benchmark values correspond to couplings that also set the relic abundance.","core_discovery":"The paper's central claim is that the fate of a dark matter spike is governed by a competition among four effects: isothermal core formation from $2\\to 2$ self-scattering, self-heating by boosted final-state particles from $n\\to m$ reactions, dissolution of the central density by number-changing processes, and (when present) $2\\to 0$ annihilation. For representative spike parameters and $n\\geq 3$ processes such as $3\\to 2$, the dissolution plateau density $\\rho_{\\mathrm{pl}} = m_\\chi \\big( (N-2)!\\,/\\,\\langle \\sigma_{N\\to M} v^{N-1}\\rangle\\, t_{\\mathrm{age}} \\big)^{1/(N-1)}$ bounds the spike density at the level needed for freeze-out, substantially flattening the inner profile. For the $2\\to 1$ semi-annihilation the plateau is not restrictive, so the spike shape is instead set by core formation and self-heating. The paper concludes that these effects significantly modify the $J$-factors for photon, neutrino, and boosted dark matter signals relative to naive NFW-based expectations.","pith_inferences":["If recapture of boosted final-state particles proves efficient, the effective particle loss per $n\\to m$ event drops from $n$ to $n-m$, weakening the dissolution depletion; the paper's plateau-density bounds would then be upper limits rather than typical densities.","The same dissolution logic should apply to other high-density dark matter environments, such as the centers of some dwarf galaxies or halos around smaller black holes, where the plateau density could be tested without relying on the Milky Way spike's uncertain stellar-heating history.","Combining the plateau-density cap with stellar-heating constraints suggests that the observable spike signal may be dominated by the outer spike region, making line-like semi-annihilation signatures more promising than searches for $n\\geq 3$ process signals.","A numerical simulation tracking boosted particles as they propagate through the spike, rather than the single-scattering efficiency $\\xi(r)$, would settle whether self-heating and dissolution act in the same direction or partially compensate."],"forward_implications":["For $n\\geq 3$ processes with freeze-out-favored cross-sections, the central spike density is capped near the plateau density, so annihilation fluxes and boosted dark matter fluxes from the inner spike are markedly lower than collisionless-spike predictions.","For $2\\to 1$ semi-annihilation, the spike structure is preserved in general, so semi-annihilating dark matter can still produce strong boosted-dark-matter signals from galactic centers.","The $J_2$ factor is enhanced relative to the NFW expectation only when $\\sigma_{2\\to 2}/m_\\chi \\lesssim 10^{-4}\\,\\mathrm{cm^2\\,g^{-1}}$ and $\\langle\\sigma_{2\\to 1}v\\rangle \\lesssim 10^{-22}\\,\\mathrm{cm^3\\,s^{-1}}$; otherwise the spike yields less signal than the naive profile.","When $\\sigma_{2\\to 2}/m_\\chi \\gtrsim 10^{-4}\\,\\mathrm{cm^2\\,g^{-1}}$, self-heating and core formation lower $J_3$ below the NFW expectation even for perturbative $3\\to 2$ couplings.","Phenomenological studies that use dark matter spikes should include these density modifications rather than assuming the bare power-law spike."],"supporting_citations":[{"why":"Supplies the adiabatic spike power-law profile that the paper modifies.","marker":"[53]"},{"why":"Supplies the cross-section parametrization and the freeze-out reference couplings used to set benchmark rates.","marker":"[27]"},{"why":"Supplies the isothermal core density and velocity profiles for self-interacting dark matter cusps around massive black holes.","marker":"[63]"},{"why":"Introduces the self-heating mechanism and the efficiency parameter used to estimate heat injection.","marker":"[44]"},{"why":"Gives the self-heating time-scale formula used to determine the self-heating core radius.","marker":"[45]"},{"why":"Gives the one-scattering-per-halo-age criterion that defines the isothermal core radius.","marker":"[61]"}],"fun_headline_variants":["Dark matter spikes flattened by 3-to-2 reactions","Spike density capped by number-changing dark matter","Self-interacting dark matter erases central spikes","n-to-m processes reshape dark matter spike profiles","Semi-annihilation preserves spikes, 3-to-2 depletes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that in every $n\\to m$ process all $n$ initial dark matter particles are lost from the local density because the $m$ final-state particles are relativistic and escape the spike, with no account of the fraction that is recaptured by strong $2\\to 2$ self-scattering before escaping.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter spikes flattened by 3-to-2 reactions","Spike density capped by number-changing dark matter","Self-interacting dark matter erases central spikes","n-to-m processes reshape dark matter spike profiles","Semi-annihilation preserves spikes, 3-to-2 depletes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1380,"prompt_tokens":944,"completion_tokens":436,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":356}},"tokens_in":560,"tokens_out":436,"duration_ms":5481,"temperature":1.0,"reasoning_tokens":356,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:46:24.389771+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the optical depth for a boosted final-state particle from an $n\\to m$ reaction to undergo a $2\\to 2$ scattering before it leaves the spike; if a non-negligible fraction is recaptured for the benchmark cross-sections, the dissolution equation overcounts particle loss and the predicted depletion and $J$-factors would need to be revised.","supporting_citations":[{"cited_title":"Core formation from self-heating dark matter","cited_arxiv_id":"1803.09762","evidence_quote":"Introduces the self-heating mechanism and the efficiency parameter used to estimate heat injection."},{"cited_title":"Escalating core formation with dark matter self-heating","cited_arxiv_id":"1911.09717","evidence_quote":"Gives the self-heating time-scale formula used to determine the self-heating core radius."}],"review_version":1}