{"id":"97302821-17e6-44bb-ba8a-9a2eda7797de","arxiv_id":"2411.18751","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"For a Burkert cored M87 halo, the smooth halo, not the black hole density spike, controls the annihilation signal, so M87 constraints on dark matter annihilation weaken by orders of magnitude.","lead":"The authors show that if M87's dark matter halo has a flat, cored center instead of a sharply peaked one, the dark matter annihilation signal from around its supermassive black hole drops enormously. This relaxes earlier claims that M87 rules out common thermal dark matter models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The relaxation claim rests entirely on the adopted 91 kpc Burkert core; the paper does not quantify how well the M87 kinematic data actually discriminate cored from cuspy halos, so the 'no spike' result may not be robust.","rationale":"I read the paper as a conditional demonstration: if M87's halo is Burkert with the [16] parameters, then the smooth-halo J-factor overwhelms the spike and the predicted gamma-ray flux lies below the [6] upper limits. The numerical machinery—the Jeans solution in Eq. (3), the truncated Maxwellian distribution in Eq. (17), and the J/Q integrals in Eqs. (19)–(20)—is internally consistent, and the order-of-magnitude estimates check out. The SIDM section is explicitly a consistency check, and its hand-tuned mass-to-light ratio is a recognized limitation rather than a fatal flaw. The load-bearing step is therefore the halo profile assumption. The paper provides no posterior or evidence for Burkert over NFW from the very data it cites, and because the original M87 constraints change by orders of magnitude between these two assumptions (the paper itself notes Qbar exceeds Jbar by about 1e5 in the NFW-based model of [6]), the central claim of relaxed constraints cannot be fully evaluated without quantifying this model uncertainty. This is exactly the weakest assumption identified by the reader, and it does not change the CONDITIONAL verdict; it sharpens what would be needed to upgrade that verdict to ACCEPT or to justify rejecting the relaxation claim.","tokens_in":15009,"tokens_out":17187,"duration_ms":168458,"concrete_test":"Refit the M87 multiple-tracer kinematics of [14,15] (or use the publicly available likelihood underlying [16]) with a family of dark-matter halo models that includes Burkert (Eq. 1), NFW (Eq. 9), and a generalized profile with a free inner logarithmic slope, using the same stellar anisotropy and mass-to-light nuisance parameters as in [16]. Compute the posterior weight or delta-BIC between the cored and cuspy models. Then draw from the posterior distribution of (rho0, r0, inner slope) and recompute Jbar and Qbar with Eqs. (19)–(20) and the flux comparison of Sec. V. If a non-negligible posterior fraction (e.g., more than 5%) yields fluxes within an order of magnitude of the [6] upper limits, or shows Qbar/Jbar greater than 1, then the central claim of relaxed constraints should be treated as conditional on the cored halo being strongly favored by the data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result—that the smooth halo dominates and M87 annihilation constraints are relaxed—is conditional on the Burkert profile of Eq. (1) being the true dark-matter halo all the way down to r_in = 27.8 pc. The paper adopts the best-fit rho0 and r0 from [16] without propagating uncertainties or giving a model comparison against cuspy profiles. This matters because the J and Q factors in Eqs. (19)–(20) are highly sensitive to the inner density slope: for an NFW cusp the spike slope is gamma_sp = 7/3 and the spike can dominate over the smooth halo. The paper itself notes that for the halo model of [6] the Q factor is about 1e5 times larger than J, which is what makes the NFW-based M87 constraints strong. The kinematic data of [14,15] were fitted in the literature with both cored and cuspy models, and [15] explicitly considered four halo models, including both families. Nothing in this paper establishes that the cored Burkert interpretation is preferred over a cuspy one; it only shows that if the core is present, the spike contribution is negligible. The conclusion is therefore a proof of principle for a cored halo, not a demonstration that M87 no longer constrains dark matter annihilation, unless the census of allowed fits to the [14,15] data is shown to favor a core. The manuscript acknowledges the sensitivity of the constraints to the halo model but does not quantify the prior or posterior weight of the cored solution.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper recalculates the dark matter annihilation J- and Q-factors for M87 using a Burkert cored halo profile taken from De Laurentis & Salucci (2022), combined with an SIDM-motivated density spike of slope 7/4 (and, for comparison, a collisionless cored-halo spike of slope 9/4). A light-mediator dark matter model is adopted, and the velocity-dependent Sommerfeld enhancement is included in the J- and Q-factor computations. The central numerical result is that, for the cored Burkert halo, the smooth halo contribution dominates over the spike contribution by roughly seven orders of magnitude, so the predicted gamma-ray fluxes are several orders of magnitude below the M87 upper limits used by Lacroix et al. (2015). The paper also constructs a self-interacting dark matter halo model including baryonic effects and shows that the conclusion persists.","tokens_in":15346,"tokens_out":37233,"duration_ms":315513,"significance":"If the assumed cored Burkert profile is the correct description of M87's dark matter halo, the paper significantly weakens previously claimed M87 spike constraints on dark matter annihilation and identifies the halo-profile choice as the dominant uncertainty in such constraints. The calculation is transparent, the inputs are clearly stated, and the paper includes useful robustness checks: varying the spike radius from 27.8 pc to 136 pc and the spike slope from 7/4 to 9/4 does not change the conclusion that the smooth halo dominates. The SIDM halo model in Sec. VI is an additional strength, as it shows that a self-consistent cored halo with baryonic contraction still leaves the predicted fluxes far below the observational limits. The main limitation is that the cored profile is adopted rather than established from the kinematic data; the paper does not quantify how strongly the data of Oldham & Auger favor a core over a cusp.","major_comments":[{"comment":"The central claim that the M87 annihilation constraints are relaxed is conditional on the Burkert cored profile of Eq. (1) being the true dark matter halo down to r_in = 27.8 pc. The manuscript adopts the best-fit values rho0 = 6.94e6 M_sun kpc^-3 and r0 = 91.2 kpc from Ref. [16] without reporting the fit quality, the parameter uncertainties, or a model comparison against cuspy halos. Since Ref. [15] explicitly fit four halo models including both cored and cuspy families, the paper should either quantify the support for the cored solution (e.g., by reporting likelihood or AIC-type comparisons from the original fits, or by performing a simple comparison with an NFW profile fit to the same mass-profile data) or explicitly reframe the title and conclusions as a conditional sensitivity study. As it stands, the conclusion that constraints are \"relaxed\" is a proof of principle for a cored halo rather than a demonstration that M87 no longer constrains dark matter annihilation.","section":"Sec. II and Sec. V"}],"minor_comments":[{"comment":"The notation in Eqs. (14)-(16) is inconsistent: Q(Omega) in Eq. (16) is written as a volume integral divided by D^2, which represents the angle-integrated Q (up to the appropriate 4 pi factor), while Eq. (14) adds it to the per-solid-angle J(Omega). Please clarify that the final flux is computed from Jbar and Qbar integrated over the observation window, and define Q(Omega) consistently throughout.","section":"Eq. (20)"},{"comment":"There is a typo in the sentence defining the power-law slope: \"rho_sp propto gamma^{-(3+a)/4}\" should read \"rho_sp propto r^{-(3+a)/4}\".","section":"Sec. IV, after Eq. (20)"},{"comment":"The phrase \"constant light-to-mass ratio\" should be \"constant mass-to-light ratio\".","section":"Sec. II, Eq. (5)"},{"comment":"The text \"fixed by the relict density constraint\" should read \"relic density constraint\".","section":"Sec. V"},{"comment":"The SIDM model fit is described qualitatively (\"after several trials\", \"fits the data well\"). Please provide a quantitative measure of the fit, such as the chi^2 or residual scatter relative to the mass-profile data points in Fig. 4, and specify the range of mass-to-light ratio and cross section explored.","section":"Sec. VI"},{"comment":"Please specify the exact origin of the black \"Upper Limits\" curve: whether it is the observed Fermi-LAT flux upper limit from Ref. [60] as converted in Ref. [6], and over which energy range the limit applies.","section":"Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the numerical machinery appears sound. The main issue is the conditional nature of the cored halo assumption: the authors should either provide quantitative evidence that the cored Burkert profile is a competitive description of the M87 kinematic data, or clearly soften the title and conclusions to a conditional sensitivity study. The robustness checks on the spike radius and slope are useful and should be kept."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a clean, honest conditional calculation. Given a cored Burkert halo for M87 with the [16] best-fit parameters, the spike contributes essentially nothing to annihilation fluxes: Q_bar is about 1e-7 of J_bar, and the projected gamma-ray fluxes stay orders of magnitude below the Fermi-derived limits. The paper applies known J/Q and Sommerfeld machinery rather than inventing new formalism, but the specific result—that a ~90 kpc core moots the SMBH spike for M87—is new and quantitatively useful. The robustness checks against r_sp from 28 to 136 pc and spike slopes 7/4 to 9/4 are the right ones and they hold. The Jeans calculation and density profile handling look internally consistent; I checked the Q_bar scaling and it matches their quoted number.\n\nThe soft spot is exactly what the stress-test note says, and the authors do not hide it: everything rests on the Burkert profile being the true halo down to ~28 pc. They take [16]'s best-fit rho0 and r0 with no uncertainties and no model comparison against cuspy NFW fits. Oldham & Auger themselves fit multiple halo models, including cuspy ones; this paper never quantifies whether the kinematic data actually prefer the core. The paper is careful to frame the result as conditional—\"motivated by recent kinematic studies,\" \"for the cored Burkert profile\"—so it is not misrepresenting the claim. But the title and abstract will be read as \"M87 no longer constrains WIMPs,\" which is only true if the core is real. To make that stronger claim, they would need to propagate halo-fit uncertainties and give the model comparison. The SIDM section is a consistency check with hand-tuned M/L and cross section; it does not resolve the main caveat, though the resulting flux enhancement of only 1.3–2 is a nice reassurance that the conclusion is not sensitive to the exact cored model.\n\nAlso worth noting: no new code or data is released beyond pointing to existing public codes, but everything needed to reproduce the numbers is described or in the cited literature. The citation pattern is fine; the central claim does not depend on self-citations.\n\nWho should read it: anyone working on black hole density spikes, Sommerfeld-enhanced indirect detection, or SIDM halo models of ellipticals. It is a legitimate paper that deserves refereeing. My own verdict would be accept with a request to add uncertainty quantification and a cored-versus-cuspy fit comparison, even a crude one. If those are added, the paper becomes definitive rather than conditional.","headline":"Clean conditional calculation: with M87's 91-kpc Burkert core, the SMBH spike is a non-factor for annihilation fluxes, but the community should still demand a cored-vs-cuspy model comparison before discarding M87 as a WIMP probe.","tokens_in":15896,"tokens_out":5634,"would_cite":true,"duration_ms":52628,"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":"If M87's dark halo is cored, its central density spike contributes almost nothing to dark matter annihilation, relaxing the galaxy's gamma-ray constraints.","keywords":["dark matter annihilation","density spike","M87","Burkert profile","Sommerfeld enhancement","J-factor","gamma-ray constraints","self-interacting dark matter"],"falsifier":"Measure M87's inner dark-matter logarithmic slope $\\gamma$ below roughly $30\\,\\mathrm{pc}$ with kinematic tracers: a detection of $\\rho\\propto r^{-1}$ (cuspy) rather than $\\rho\\simeq$ constant (cored) would raise the spike normalization by orders of magnitude, pushing the projected fluxes toward or above the upper limits derived from [6] and overturning the claimed relaxation.","tokens_in":14753,"feed_emoji":"🌌","tokens_out":8464,"duration_ms":73312,"temperature":0.7,"pith_summary":"This paper tries to establish that M87's black-hole density spike does not dramatically boost dark matter annihilation signals if the galaxy's dark halo is cored rather than cuspy. It adopts the Burkert profile fitted to recent kinematics of stars, globular clusters, and satellite galaxies, and computes velocity-dependent annihilation factors with and without Sommerfeld enhancement. The smooth halo dominates: the $J$-factor is about $8\\times10^{17}\\,\\mathrm{GeV^2\\,cm^{-5}}$, while the spike $Q$-factor is about $6\\times10^{10}\\,\\mathrm{GeV^2\\,cm^{-5}}$, so the spike is negligible. The resulting gamma-ray fluxes are several orders of magnitude below M87 upper limits, meaning the constraints are far weaker than in the widely used NFW spike analysis.","feed_headline":"Cored M87 halo erases the black-hole dark matter spike boost","feed_subtitle":"With a Burkert profile the smooth halo dominates annihilation, so M87 gamma-ray limits no longer bite.","key_machinery":"The argument runs on the ratio of two velocity-weighted line-of-sight integrals: the $J$-factor for the smooth Burkert halo and the $Q$-factor for the spike, which together set the annihilation flux through $d^2\\Psi/dE_\\gamma d\\Omega \\propto J(\\Omega)+Q(\\Omega)$. Both are evaluated with an isotropic Maxwell-Boltzmann velocity distribution truncated at the escape velocity, with Sommerfeld enhancement factors for s- and p-wave annihilation folded into the integrals. The decisive comparison is $\\bar J\\sim 8\\times10^{17}\\,\\mathrm{GeV^2\\,cm^{-5}}$ versus $\\bar Q\\sim 6\\times10^{10}\\,\\mathrm{GeV^2\\,cm^{-5}}$ for the $r^{-7/4}$ spike, a ratio of $10^7$ that makes the spike irrelevant.","core_discovery":"For a cored Burkert halo with central density $\\rho_0=6.94\\times10^6\\,M_\\odot\\,\\mathrm{kpc}^{-3}$, scale radius $r_0=91.2\\,\\mathrm{kpc}$, and black-hole influence radius $r_{\\rm in}=27.8\\,\\mathrm{pc}$, the dark matter spike near M87's supermassive black hole contributes negligibly to annihilation signals. The smooth halo gives $\\bar J\\sim 8\\times10^{17}\\,\\mathrm{GeV^2\\,cm^{-5}}$; the spike, modeled with slope $r^{-7/4}$ for Coulomb-like self-interactions, gives $\\bar Q\\sim 6\\times10^{10}\\,\\mathrm{GeV^2\\,cm^{-5}}$, seven orders of magnitude smaller. The conclusion also holds for a steeper $r^{-9/4}$ collisionless spike and for a larger influence radius of $136\\,\\mathrm{pc}$. For a light-mediator dark matter model with s- and p-wave annihilation and Sommerfeld enhancement, the projected gamma-ray fluxes are several orders of magnitude below the upper limits derived from [6]; a self-consistent SIDM halo that includes the baryonic potential raises the flux by a factor of only $1.3$--$2$ but still leaves it well below the limits.","pith_inferences":["The relaxation is empirical, not model-independent: if future kinematic data inside roughly $30\\,\\mathrm{pc}$ favor a cusp, the old strong M87 bounds would come back with the same calculation.","Since the spike's $Q$-factor is negligible, M87's black-hole spike is not the right target for probing light-mediator annihilation; the observable signal is the extended smooth-halo emission, which requires different angular integration and background treatment.","The same $J/Q$ comparison suggests that other giant ellipticals with kinematically inferred cored halos may have overestimated black-hole spike constraints until their inner slopes are measured.","A direct test is to measure the inner logarithmic density slope of M87 below $\\sim30\\,\\mathrm{pc}$; a cuspy slope near $\\gamma=1$ would flip the paper's central conclusion."],"forward_implications":["M87's gamma-ray upper limits would no longer rule out thermal-relic WIMPs or light-mediator models if the halo is cored.","The spike contributes negligibly, so future M87 analyses should base constraints on the smooth-halo $J$-factor rather than the spike $Q$-factor.","The NFW-based spike model of [6] gives a $Q$-factor $10^5$ times the smooth $J$-factor, so the same observational limits yield vastly different cross-section bounds depending on the assumed halo profile.","An SIDM halo that includes the baryonic potential raises predicted fluxes by only a factor of $1.3$--$2$, leaving them far below the observational upper limits."],"supporting_citations":[{"why":"Supplies the cored Burkert halo parameters $\\rho_0=6.94\\times10^6\\,M_\\odot\\,\\mathrm{kpc}^{-3}$ and $r_0=91.2\\,\\mathrm{kpc}$ that replace the NFW profile.","marker":"[16]"},{"why":"Provides the NFW-based M87 constraints and gamma-ray upper limits against which the predicted fluxes are compared.","marker":"[6]"},{"why":"Gives the self-interacting dark matter spike slope $\\gamma_{\\rm sp}=(3+a)/4$, used to set the $r^{-7/4}$ spike for $a=4$.","marker":"[13]"},{"why":"Supplies the collisionless spike slope and the adiabatic-growth argument that motivates density spikes around black holes.","marker":"[4]"},{"why":"Defines the Burkert density profile used for the smooth halo component.","marker":"[17]"},{"why":"Earlier SIDM spike study with an NFW smooth halo; the present work extends it to a cored halo and finds even weaker constraints.","marker":"[11]"},{"why":"Provides the s- and p-wave Sommerfeld enhancement factors used in the velocity-dependent $J$ and $Q$ integrals.","marker":"[29]"},{"why":"Sets up the light-mediator model, relic-density calculation, and self-interaction benchmark used for the example flux predictions.","marker":"[27]"}],"fun_headline_variants":["Cored M87 halo deflates dark matter spike","M87's Burkert halo mutes black hole dark matter spike","Dark matter annihilation spike vanishes in cored M87","Cored halo soothes M87 dark matter annihilation constraints","M87 halo shape relaxes dark matter limits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that M87's dark halo follows the cored Burkert profile of Eq. (1) with $\\rho_0=6.94\\times10^6\\,M_\\odot\\,\\mathrm{kpc}^{-3}$ and $r_0=91.2\\,\\mathrm{kpc}$ all the way down to the black-hole radius of influence $r_{\\rm in}=27.8\\,\\mathrm{pc}$; if the inner profile is cuspy instead, the spike can dominate and the constraints tighten.","fun_headline_variants_meta":{"raw":{"variants":["Cored M87 halo deflates dark matter spike","M87's Burkert halo mutes black hole dark matter spike","Dark matter annihilation spike vanishes in cored M87","Cored halo soothes M87 dark matter annihilation constraints","M87 halo shape relaxes dark matter limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000221,"raw_usage":{"total_tokens":1456,"prompt_tokens":960,"completion_tokens":496,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":417}},"tokens_in":576,"tokens_out":496,"duration_ms":5609,"temperature":1.0,"reasoning_tokens":417,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:55:00.000348+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure M87's inner dark-matter logarithmic slope $\\gamma$ below roughly $30\\,\\mathrm{pc}$ with kinematic tracers: a detection of $\\rho\\propto r^{-1}$ (cuspy) rather than $\\rho\\simeq$ constant (cored) would raise the spike normalization by orders of magnitude, pushing the projected fluxes toward or above the upper limits derived from [6] and overturning the claimed relaxation.","supporting_citations":[],"review_version":1}