{"id":"f3fe90a9-942f-4bd1-8c4d-a66b38b30027","arxiv_id":"2411.12007","paper_version":3,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":10,"one_line_summary":"N-body simulations of adiabatic black hole growth in Hernquist haloes yield an empirical spike profile with outer halo depletion and a mass-ratio-dependent spike radius that deviates from analytical predictions.","lead":"Dark matter spikes around black holes have so far been predicted analytically, but this paper runs the first N-body simulations of spike formation in Hernquist haloes and fits a new empirical density profile that depends only on the mass ratio of the black hole to the halo. This matters because dark matter detection forecasts, such as gamma-ray annihilation signals and gravitational wave dephasing from inspirals, depend directly on the assumed spike profile.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fitted spike parameters lie far below the stated convergence radius; the empirical profile and mass-ratio scalings are therefore unsupported by the simulation data.","rationale":"The manuscript's central claim is that the N-body data determine the empirical profile (Eq. 12) with the mass-ratio scalings of Eqs. 14–15. This requires the fitted r_sp and γ_sp to be constrained by converged, physically meaningful particles. That condition fails: using the paper's own resolution criterion (§3.1 and Appendix B), r_conv = 15 r_vir/√N ≈ 0.19 kpc for the 1e4-1e3 run, while Eq. 15 gives r_sp ≈ 0.0037 kpc; for the largest-μ run, r_sp ≈ 0.0086 kpc, still ~20× below r_conv. At r = r_conv the spike term contributes only ~0.5–2% of the density for γ_sp = 7/3–2, so the spike term is not measurable in the converged region. The claimed 1σ errors on α_1...α_6 and the χ²_red ~ 10^-5 values are incompatible with the stated Poisson error model; indeed §3.3 states that 'more than 10^4 particles are present in the least populated bins' for runs with only 1,303 particles, which is internally impossible. The artificial-data test in Appendix C2 does not rescue the claim: its 'zoomed' range starts at 10^-3, which brackets r_sp, whereas the N-body fits start at r_conv ≫ r_sp. The paper honestly flags low-μ γ_sp as uncertain, but the same unresolved range is used for r_sp and β fits at all μ. An additional internal inconsistency is the conclusion's 'depletion up to 20%' versus Eq. 14, which implies ~39% depletion at μ = 0.333. Thus the load-bearing condition is not met: the central empirical profile and its scalings are unsupported by the presented simulation data. The reader's REJECT verdict is appropriate and should remain unchanged.","tokens_in":16434,"tokens_out":13539,"duration_ms":138682,"concrete_test":"Reanalyse the final snapshot of the 1e4-1e3 run using only bins above r_conv = 0.187 kpc. Fit two models: (i) Eq. 12 with r_sp and γ_sp free, and (ii) a pure depleted Hernquist profile (β free, no spike term). Compute Δχ²/dof on the actual binned counts with a proper covariance matrix that accounts for bin correlations. If the spike term does not significantly improve the fit, then r_sp is unconstrained and the claimed Eq. 15 is an artifact of extrapolating below the resolution limit; conversely, a decisive Δχ² would support the profile even if r_sp lies below r_conv.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The fitted spike lies below the declared resolution. For the 1e4-1e3 run (Table 1), ε = 6 r_vir/√N ≈ 0.075 kpc and r_conv = 2.5 ε ≈ 0.187 kpc, while Eq. 15 with μ ≈ 0.074 gives r_sp ≈ 0.0037 kpc — about 50× smaller. Even the largest-μ run (1e4-5e3, μ = 0.333) gives r_sp ≈ 0.0086 kpc, still ~22× below r_conv. At r = r_conv the spike term (r/r_sp)^(1−γ_sp) is only ~0.5–2% of the Hernquist density for γ_sp = 7/3–2, so fits restricted to r > r_conv cannot determine r_sp or γ_sp. The small formal errors quoted in Eqs. 14–15 and the reported χ²_red ~ 10^-5 are therefore not physical constraints. Appendix C2's validation uses a 'zoomed' range 10^-3 ≤ r ≤ 10^0 that brackets the artificial r_sp = 0.05/0.25, unlike the actual N-body fits which start at r_conv ≫ r_sp; it does not test the problematic configuration. The paper itself concedes that low-μ spikes 'manifest below r_conv' and that low-μ γ_sp may be a resolution artifact, yet the same unresolved radius range is used to fit r_sp for all μ. Hence the central empirical profile and its mass-ratio scalings rest on unresolved data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents seven N-body simulations of Hernquist dark matter haloes (N ≈ 1300 particles) with a central black hole grown adiabatically in a modified version of the SWIFT code. From 235 recorded snapshots spanning mass ratios μ = M_BH/M_tot, the authors propose an empirical post-growth density profile (Eq. 12) consisting of the original Hernquist profile multiplied by a depletion factor β plus a power-law spike with slope γ_sp and radius r_sp. They fit β(μ), r_sp(μ), and γ_sp(μ), report scalings in Eqs. (14) and (15), and use the HaloFeedback code to estimate gravitational-wave dephasing for the new profile. The paper claims to be the first fully numerical demonstration of dark matter spike formation.","tokens_in":16725,"tokens_out":11504,"duration_ms":112890,"significance":"If the central claims were supported, this would be a significant contribution: it would offer the first N-body check of adiabatic spike formation in Hernquist haloes, a simple one-parameter empirical profile, and updated predictions for gravitational-wave dephasing. The authors deserve credit for releasing their code, for presenting detailed appendices on the convergence radius and on fitting validation, and for candidly acknowledging some resolution limitations. However, the main claims about the spike radius and slope are extracted from scales far below the stated convergence radius, and the validation appendix does not test the actual fitting configuration. The headline profile and mass-ratio scalings are therefore not established by the presented simulations.","major_comments":[{"comment":"The fitted spike parameters lie far below the resolution limit. For the 1e4-1e3 run, ε = 6 r_vir/√N ≈ 0.075 kpc and, with the adopted δ = 2.5, r_conv = 2.5ε ≈ 0.19 kpc. Equation (15) with μ = 0.074 gives r_sp ≈ 0.0037 kpc, about 50 times smaller than r_conv; even for the highest-μ run (1e4-5e3, μ = 0.333), r_sp ≈ 0.0086 kpc, still about 22 times below r_conv. Since the fitting lower boundary is r_conv, the spike term in Eq. (13) contributes only (r_conv/r_sp)^{1-γ_sp} ≈ 0.5% for γ_sp = 7/3 and about 2% for γ_sp = 2 at the first fitted radius, and it decreases outward. The binned density above r_conv is therefore essentially β times the original Hernquist profile and contains no direct information about r_sp or γ_sp. The two-orders-of-magnitude smaller RMSE reported in Table 2 is an artifact of fitting an extrapolated, unresolved parameter and does not establish a new scaling. This undermines the central claim of the abstract and of Section 4.1.","section":"Sec. 4.1, Table 1, Eq. (15)"},{"comment":"The validation in Appendix C2 does not mimic the actual fitting configuration. The 'zoomed' fit in Tables C1-C4 covers 10^{-3} ≤ r ≤ 10^0 in units of a and brackets the injected artificial spike radii r_sp = 0.05 and 0.25. In the real N-body fits, the lower boundary is r_conv ≈ 10a for the 10^4 M_sun haloes (0.19 kpc versus a = 0.019 kpc), so the actual fitting interval starts at roughly 50 r_sp. The artificial-data tests therefore only demonstrate parameter recovery when the fitted range contains the spike; they do not address the situation in which the spike lies far below the first fitted radius. The paper's own admissions in Section 4.1 that low-μ spikes 'manifest below r_conv' and that the low-μ γ_sp values may be resolution artifacts apply equally to r_sp, so the scaling in Eq. (15) is not supported by the validation presented.","section":"Appendix C2, Sec. 4.1"},{"comment":"The reported bin counts and fit quality are internally inconsistent. Section 3.3 states that more than 10^4 particles are present in the least populated radial bins, but each run has only N = 1303 DM particles; with logarithmic radial binning, the innermost resolved bins can contain at most a few tens to a few hundred particles, not 10^4. The quoted χ²_red values of order 10^{-5} for the secondary fits in Eqs. (14) and (15) and the small 1σ errors therefore cannot follow from the stated Poisson error model applied to single snapshots. Either the error model is mis-described, or the 235 snapshots are being treated as independent even though they are strongly correlated within each of the seven runs. In either case, the formal significance assigned to the empirical scalings is not trustworthy.","section":"Sec. 3.3, Eqs. (14)-(15)"},{"comment":"The gravitational-wave dephasing estimates inherit the resolution problem. Table 3 compares inspirals in the proposed profile with γ_sp = 7/3 and γ_sp = 2, but both cases use an r_sp that is unconstrained by the simulations, and the γ_sp = 2 case uses a slope that the authors themselves attribute to limited fitting range. The qualitative statement that dephasing can be smaller for shallower spikes is reasonable, but the numerical values in Table 3 and the associated conclusions in Section 5 should not be presented as predictions of the simulated profile until the spike parameters are resolved.","section":"Sec. 4.2, Table 3"}],"minor_comments":[{"comment":"Section 3.1 states r_conv = 2ε, while Appendix B concludes that δ = 2.5 should be used; please reconcile this inconsistency.","section":"Sec. 3.1 vs Appendix B"},{"comment":"The caption for panel (c) says the lower sub-panel shows 'values of β divided by the best fit'; it should refer to r_sp.","section":"Fig. 2 caption"},{"comment":"The notation \\(\\tilde{r}_{\\rm sp}\\) in Eq. (12) is not defined until Eq. (13); please state explicitly that \\(\\tilde{r}_{\\rm sp} = r_{\\rm sp}/a\\).","section":"Eqs. (12)-(13)"},{"comment":"The phrase 'fully numerically simulated cold dark matter spikes' in the abstract is stronger than what the resolution analysis in Section 4.1 and Appendix C2 supports; please qualify it.","section":"Abstract and conclusions"},{"comment":"Making the reduced density profiles and the fit catalog publicly available, rather than only upon reasonable request, would substantially improve reproducibility.","section":"Data availability"}],"recommendation":"reject","confidential_remarks":"This is a candid and well-structured first attempt, and the authors are explicit about some limitations. However, the core problem is not a presentation issue: the fitted spike radius is 20-50 times below the stated convergence radius in every run, the artificial-data validation in Appendix C2 tests a different regime, and the reported bin counts are impossible for N = 1303. These issues affect the validity of the headline profile and scalings. A revision would require either a much higher-resolution simulation campaign (the authors themselves estimate ~5×10^10 particles for a reliable γ_sp determination) or a substantial reduction of scope to the depletion parameter β. I therefore recommend rejection rather than major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The novelty is real. As far as I know, this is the first fully N-body attempt to form a dark matter spike by growing a black hole adiabatically in a Hernquist halo. The code is on GitHub, the authors check for adiabaticity, and they are candid that low-mu spikes live below r_conv and that gamma_sp is hard to determine. Good faith effort.\n\nBut the central claim, the empirical profile of Eq. (12) with its mu-dependent r_sp and beta, does not hold. The fitted r_sp is below the convergence radius by factors of 20–50. For the 1e4–1e3 run, r_conv ~ 0.19 kpc while the fitted r_sp is ~ 0.004 kpc. At r = r_conv, the spike term is about half a percent of the Hernquist density, so r_sp and gamma_sp are essentially unconstrained. The Appendix C2 validation uses artificial data where the zoomed range brackets r_sp; the actual fits start at r_conv, which is far outside the spike. The paper even admits low-mu gamma_sp may be a resolution artifact, yet still quotes tight errors on r_sp for all mu. So the main product rests on unresolved data.\n\nThere are also some sloppy inconsistencies. Eq. (3) gives r_vir/a ~ 1e8 for a 1e4 Msun halo, while Table 1 says c ~ 24. The conclusion says depletion up to 20%, but Eq. (14) gives 39% at the highest mu. The text says bins have >1e4 particles when the runs have 1303 particles. And the reported chi^2_red values ~ 1e-5 are implausibly small, suggesting the error model is off by orders of magnitude. None alone kills the paper, but together they erode confidence.\n\nWhat might survive is the depletion parameter beta, since it is fit over r > r_conv where the spike term is negligible. But the paper does not cleanly separate that.\n\nVerdict: this is a genuinely new attempt at an important question and deserves a serious referee. But the central result is not supported by the data. I would reject as is, and ask for either a resolution study with many more particles or a paper restricted to what is actually resolved, like the depletion effect. I would not cite it yet, but it is worth a reading group discussion on methodology.","headline":"First fully numerical attempt at DM spikes, but the fitted spike parameters lie far below the convergence radius, so the empirical profile is unsupported.","tokens_in":17369,"tokens_out":5921,"would_cite":false,"duration_ms":54813,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d"],"model":"deepseek-v4-flash","headline":"Dark matter spikes formed by adiabatic black hole growth follow a new empirical profile, with radius and halo depletion set by a single mass ratio.","keywords":["dark matter spikes","N-body simulations","Hernquist halo","adiabatic black hole growth","empirical density profile","halo depletion","gravitational wave dephasing","intermediate-mass black holes"],"falsifier":"A simulation with enough particles to resolve radii below the predicted spike radius (≈0.002 kpc for the $10^{4}$ M_sun halo) would settle whether the fitted r_sp and γ_sp are physical: if the density in that region does not follow the power-law spike of Eq. (12) with the same fitted parameters, the proposed scalings are numerical artifacts of the limited resolution.","tokens_in":16101,"feed_emoji":"🕳️","tokens_out":10589,"duration_ms":90971,"temperature":0.7,"pith_summary":"The paper reports the first fully numerical N-body simulations of dark matter spike formation by adiabatic growth of a central black hole in Hernquist haloes. It proposes that the final density profile is an empirical formula depending only on the black-hole-to-halo mass ratio μ, with a spike radius that scales differently from the analytical Gondolo-Silk and Merritt predictions and an outer-halo depletion that becomes significant for μ ≳ 0.1. If correct, this changes the expected dark matter densities near intermediate-mass black holes and, consequently, the predicted gravitational-wave dephasing of inspiraling compact objects and the prospects for indirect dark matter detection.","feed_headline":"Dark matter spikes differ from theory in first N-body simulations","feed_subtitle":"A single mass-ratio parameter sets the spike radius and outer-halo depletion, changing predicted gravitational-wave dephasing.","key_machinery":"The central object is the empirical spike profile of Eq. (12), a multiplicative factor on the Hernquist density that combines a depletion term β and a broken power-law spike (r/r_sp)^{1-γ_sp}; the scaling relations of Eqs. (14)-(15) reduce the profile to a function of a single parameter, the mass ratio μ = M_BH/M_tot. The numerical scheme is a modified version of the SWIFT N-body code with a growing point-mass 'DAB' black hole, Hernquist initial conditions drawn from the Eddington distribution function, and the Power et al. (2003) softening with convergence radius r_conv = 2.5ε that defines the resolved fitting range. The gravitational-wave dephasing is computed with the HaloFeedback code of Kavanagh et al. (2020).","core_discovery":"Using the modified SWIFT code with a growing point-mass black hole ('DAB'), the authors simulate seven Hernquist haloes with masses between 3×$10^{3}$ and $10^{5}$ M_sun and black holes grown adiabatically to $10^{3}$–5×$10^{3}$ M_sun, recording 235 snapshots. They fit the final density as ρ(r) = ρ_Hernq(r)[β + (r/r_sp)^{1-γ_sp}] and find that the best-fit depletion and spike radius depend only on μ: β = 1 − 0.998 $μ^{{0.858}}$ and r_sp/a = 0.801 $μ^{{2.29}}$/($μ^{{1.78}}$ + 9.1×$10^{{-4}}$). The spike slope γ_sp is consistent with 7/3 for μ ≳ 0.06 but is not well constrained below; the depletion reaches β ≈ 0.8 at the highest μ. Compared with the 'Modified G&S' spike, the new profile changes the gravitational-wave dephasing of a $10^{3}$ M_sun primary with a solar-mass secondary by up to a factor of two at γ_sp = 7/3, and nearly eliminates the dephasing if the slope is shallower (γ_sp = 2).","pith_inferences":["Editorially: the steep low-μ scaling r_sp ∝ μ^2.29, if real, implies a much stronger dependence of spike size on black hole mass than adiabatic theory's μ^0.5; a simulation with higher resolution at low μ would test whether this is a physical effect or a fitting artifact.","Editorially: the depletion of the outer halo might be observable in the rotation curves or stellar kinematics of dwarf galaxies hosting IMBHs, since β < 1 changes the enclosed mass at radii of order a; a targeted observational search could constrain the profile independently of the simulations.","Editorially: the authors' own convergence analysis shows that γ_sp is the least constrained parameter; a next step would be to run a single simulation with much higher particle number (they estimate ~5×10^10 particles for the 10^4 M_sun halo) to pin down the slope, rather than adding more μ values at fixed resolution."],"forward_implications":["If the profile is correct, the dark matter density around intermediate-mass black holes is lower in the outer spike region than the standard G&S spike for the same μ, and the outer halo is depleted by up to ~20% at μ ~ 0.25.","The spike radius scaling r_sp ~ 0.8 a √μ at high μ and a steeper power at low μ replaces the often-used r_sp = r_h/5, so analyses that assume the Merritt radius will misestimate the spike's normalization and extent.","Gravitational-wave dephasing forecasts for LISA (extreme and intermediate mass-ratio inspirals) should be re-evaluated: the new profile can double the dephasing for a 10^4 M_sun halo at γ_sp = 7/3, or make it nearly vanish if the low-μ slope is actually 2.","The claim that results transfer to NFW haloes at fixed μ (due to the identical inner cusp) means the empirical profile, if confirmed, would apply to cosmologically motivated haloes without modification."],"supporting_citations":[{"why":"Predicts the γ=7/3 spike slope and the analytical spike profile that the N-body results are compared against.","marker":"Gondolo & Silk (1999)"},{"why":"Defines the radius of gravitational influence and the spike radius r_sp = r_h/5 underlying the Modified G&S profile.","marker":"Merritt (2004)"},{"why":"Provides the initial halo density profile and potential used to set up the simulations.","marker":"Hernquist (1990)"},{"why":"Defines the gravitational softening length and convergence radius criterion that sets the resolved radial range.","marker":"Power et al. (2003)"},{"why":"The SWIFT N-body code whose gravity solver and time-stepping were modified to grow a black hole.","marker":"Schaller et al. (2024)"},{"why":"Provides the HaloFeedback code used to compute gravitational-wave dephasing in the simulated profiles.","marker":"Kavanagh et al. (2020)"},{"why":"Supplies the numerical implementation of the Gondolo-Silk formalism used to cross-check the N-body results.","marker":"Bertone et al. (2024)"},{"why":"Gives the mass–concentration relation used to determine the halo scale radius a from the halo mass.","marker":"Correa et al. (2015)"}],"fun_headline_variants":["First N-body dark matter spikes defy theory","Dark matter spike simulations overturn old predictions","One parameter sets spike radius and outer-halo depletion","N-body sims reveal dark spikes scale differently"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fitted spike parameters r_sp and γ_sp are treated as physical even though the spike radius lies far below the convergence radius, so the spike's shape and scalings are inferred by extrapolation rather than directly resolved.","fun_headline_variants_meta":{"raw":{"variants":["First N-body dark matter spikes defy theory","Dark matter spike simulations overturn old predictions","One parameter sets spike radius and outer-halo depletion","N-body sims reveal dark spikes scale differently"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000825,"raw_usage":{"total_tokens":3634,"prompt_tokens":997,"completion_tokens":2637,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":2579}},"tokens_in":613,"tokens_out":2637,"duration_ms":18516,"temperature":1.0,"reasoning_tokens":2579,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:02:40.766731+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A simulation with enough particles to resolve radii below the predicted spike radius (≈0.002 kpc for the $10^{4}$ M_sun halo) would settle whether the fitted r_sp and γ_sp are physical: if the density in that region does not follow the power-law spike of Eq. (12) with the same fitted parameters, the proposed scalings are numerical artifacts of the limited resolution.","supporting_citations":[],"review_version":1}