{"id":"fd41af16-6b13-4231-8f32-213188ec890d","arxiv_id":"2511.12570","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A dark-matter spike built from the full orbital motion of its particles has ~50% more energy density near the black hole and produces metric deviations ~2.5 times larger than mass-only models.","lead":"This paper calculates how a dense 'spike' of dark matter around a black hole warps spacetime, including the motion of the dark matter particles rather than treating them as a motionless fluid. It finds the moving particles add about 50% more energy density near the spike and produce spacetime distortions roughly 2.5 times larger than earlier, simpler models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Factor-2.5 backreaction claim conflates EMT content with choice of truncated Einstein equations; no same-method control.","rationale":"The reader's weakest assumption concerns the unquantified residual of the dropped θθ/conservation equation. That is a valid self-consistency concern, but it may be numerically benign because the violation is first-order in the metric deviation (a relative effect of ~10^-7). The more direct threat to the paper's headline number is that the 2.5 factor is computed against a different solution scheme (Ref. [47] solves the full system with pr=0). The paper's own two-equation truncation could systematically inflate the deviation relative to a fully constrained solve, independent of the physics. A same-method control would isolate the EMT contribution. This is a concrete, feasible check that does not require solving the full self-consistent problem. The EMT derivation (Sec. II) and the energy-condition checks are solid; the issue is confined to the backreaction comparison. I therefore agree with the reader's CONDITIONAL verdict: the paper should provide either the control calculation or a quantitative demonstration that the truncation residual is negligible. My concern does not change the verdict, so I set UNCHANGED.","tokens_in":16104,"tokens_out":13065,"duration_ms":109145,"concrete_test":"Using the same spectral method and boundary conditions as in Sec. III.B, recompute the metric with the source T replaced by the mass-density-only version (ρ_M from Eq. (35), p_r=p_t=0). Compare the ratios |g−g_Sch|_full/|g−g_Sch|_mass and |m−m_Sch|_full/|m−m_Sch|_mass near the ρ_E peak. If the ratios are not ≈2.5, the claimed factor is contaminated by the different truncation. Also evaluate the residual of Eq. (47) with the full-T metric; if it exceeds a few percent of the local |g−g_Sch|, the truncated solution is not a good approximation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that the full dynamical EMT produces metric deviations ~2.5 times larger than a mass-density-only treatment—rests on a comparison between two inequivalent calculations. The authors solve only the tt and rr Einstein equations (43)-(44) with their full T (ρ_E, p_r≠0), while the alternative method (Ref. [47]) solves the full system (43),(44),(47) with p_r=0. Thus the 2.5 factor conflates the change in the EMT with the change in which subset of Einstein equations is enforced. The dropped θθ/conservation equation (47) is knowingly not satisfied in the perturbed metric, and its residual is never quantified; if that residual is comparable to the deviations (which are only O(10^-7) in g), the reported ratio is an artifact of the truncation rather than a physical effect of kinetic energy. The text attributes the enhancement primarily to kinetic energy (Sec. IV), but no same-method control—e.g., solving the same two equations with the mass-density-only T—is presented. Without that control, one cannot distinguish the physical contribution of the full EMT from the bias introduced by the underdetermined system.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs the energy-momentum tensor of a dark-matter spike formed by adiabatic growth of a Schwarzschild black hole inside a Hernquist halo, using the Einstein-cluster distribution-function formalism. From the Hernquist DF and conserved actions in the Schwarzschild background, the authors derive integral expressions for the mass current and the diagonal EMT components (Eqs. (35)-(38)), including a kinetic contribution to the energy density and non-vanishing radial pressure. They report that near the spike the kinetic energy enhances the energy density by roughly 50% relative to the rest-mass density, that the EMT satisfies the standard energy conditions, and that the backreacted static spherically symmetric metric, obtained by solving only the tt and rr Einstein equations (43)-(44) with this EMT as a fixed source, deviates from Schwarzschild by up to about 2.5 times more than the mass-density-only treatment of Ref. [47]. The paper explicitly acknowledges that the conservation equation (47) is not satisfied in the perturbed metric and calls the construction semi-consistent.","tokens_in":16371,"tokens_out":3919,"duration_ms":37806,"significance":"If established, the result would be a useful quantitative step beyond mass-density-only treatments of DM-spike backreaction, showing that orbital kinetic energy and pressure can matter for the metric. The main strength is the explicit phase-space construction: the integrals (35)-(38) are concrete, depend only on macroscopic Hernquist parameters, and the paper checks the energy conditions numerically. There is no obvious circularity: the kinetic enhancement and the metric deviation are outputs of the orbital-invariant calculation, not fitted targets. However, the central comparison is currently not a controlled one, and the unquantified violation of the conservation equation leaves the reported metric-source correspondence vulnerable.","major_comments":[{"comment":"The paper's central claim that the full dynamical EMT produces metric deviations roughly 2.5 times larger than mass-density-only treatments is not established by the comparison shown. The curves labelled EC from Ref. [47] are obtained by solving the full system (43),(44),(47) with zero radial pressure, whereas the present solution uses the full EMT but solves only (43) and (44), explicitly dropping (47). The factor 2.5 therefore conflates two independent changes: the content of the source and the subset of Einstein equations enforced. A same-method control, e.g., solving exactly (43)-(44) with the mass-density-only source built from Eq. (35), is not presented. Without that control, the attribution of the enhancement to kinetic energy is underdetermined.","section":"Sec. IV, Fig. 5; Sec. III.A, Eqs. (43)-(44) and (47)"},{"comment":"The paper states that after backreaction is included, 'this statistical construction no longer guarantees exact conservation, hence (47) cannot be simultaneously fulfilled in the perturbed metric.' The size of the conservation residual is never estimated. Since the reported deviations are small (g deviations of order 10^-7-10^-6 and m deviations of order 10^-5 near the spike), the residual in ∇_μ T^{μν} could be comparable to or larger than the claimed effect. The authors should quantify the residual in the solved metric, or otherwise show that it is negligible. As it stands, the metric obtained from the truncated system cannot be unambiguously identified as the metric sourced by the claimed T^{μν}.","section":"Sec. III.A, Eq. (47)"},{"comment":"The choice to solve only (43) and (44) is motivated by the fact that (46)/(47) is automatically satisfied in the fixed Schwarzschild background. That motivation is insufficient: consistency in the background does not imply that the residual remains small once g and m are changed by O(10^-6)-O(10^-5). A linearized estimate of the conservation residual in the perturbed metric, or an iterative scheme that recomputes the EMT in the updated metric, is needed to support the statement that this truncation 'captures the dominant physical effects'.","section":"Sec. III.A"}],"minor_comments":[{"comment":"Typo 'symmertry' should be 'symmetry'.","section":"Eq. (40) and surrounding text"},{"comment":"The right panels plot ratios of deviations against the alternative method. The caption and axis labels should clarify that the comparison is between two different solution schemes, not two sources in the same scheme; otherwise the reader may misread the factor 2.5 as a pure EMT effect.","section":"Fig. 5"},{"comment":"Reference [35] has a formatting error: 'arXiv:2508.20238 [gr-qc]]' has a misplaced bracket. Several other references would benefit from journal/volume updates where available.","section":"References"},{"comment":"The adiabatic-invariant relation is central but the numerical evaluation of ϵ(E,L^2) is not described. A brief statement of the root-finding procedure and accuracy would improve reproducibility.","section":"Sec. II.B, Eq. (22)"}],"recommendation":"major_revision","confidential_remarks":"The paper has a solid and interesting EMT construction, but the headline factor-2.5 comparison is not controlled and the conservation residual is unquantified. Both issues appear fixable within the manuscript's scope: a same-method control and a residual estimate, or a fully consistent iterative scheme, would determine whether the enhancement is physical or an artifact of the truncated equations. I therefore recommend major revision rather than rejection. No concerns about novelty disclosure or citation pattern were identified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper's real contribution is the explicit full energy-momentum tensor for a dark-matter spike in the Einstein cluster picture: kinetic energy density, radial pressure, and tangential pressure from the Hernquist distribution function after adiabatic growth. That part is done cleanly. The peak energy density comes out about 50% above the rest-mass density, and the tensor satisfies the standard energy conditions. These are genuine, parameter-free outputs from the orbital integrals, not fit to anything. A nice property is that the EMT is independent of the DM particle mass, depending only on macroscopic halo parameters and the BH mass. I would cite this paper for the EMT.\n\nThe backreaction section is where it gets shaky. The claim that the metric deviation is 2.5 times larger than previous mass-density-only treatments is not supported as stated. The authors solve only the tt and rr Einstein equations using a T fixed in a Schwarzschild background, and compare that to a full-system solution with p_r=0. That comparison conflates the change in the EMT with the change in which subset of Einstein equations is enforced. They also knowingly drop the conservation equation (47) and never quantify its residual in the perturbed metric. Since the deviations are only 1e-7 to 1e-5 in geometric units, an unquantified residual of comparable relative size could change the ratio. A same-method control—solving the same two equations with a mass-density-only T—is missing, so the paper's attribution of the factor 2.5 primarily to kinetic energy is not established.\n\nYour stress-test note is on target. The paper is honest about the limitation and frames the metric result as coming from a semi-consistent iterative procedure, but the abstract's wording invites the reader to take the factor 2.5 at face value. The EMT part stands; the backreaction number should be treated as tentative until someone either iterates the DF in the perturbed metric or at least quantifies the conservation residual and runs the ρ_M-only control.\n\nNo code or data are included, so reproducibility is moderate. The citation pattern looks normal; the key prior work is credited.\n\nWho gets value: people modeling DM spikes around BHs, especially the gravitational-wave and lensing communities. It deserves peer review, with the backreaction section needing a clear revision before the factor 2.5 becomes a citable result.","headline":"Solid EMT derivation; the factor-2.5 backreaction comparison is not apples-to-apples and should be treated as provisional.","tokens_in":16867,"tokens_out":6603,"would_cite":true,"duration_ms":55108,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.20.-q","04.40.-b","95.35.+d"],"model":"deepseek-v4-flash","headline":"Accounting for the orbital kinetic energy and radial pressure of a dark-matter spike makes the black-hole spacetime deviate from Schwarzschild about 2.5 times more than rest-mass-only treatments find.","keywords":["dark matter spike","black hole backreaction","Einstein cluster","energy-momentum tensor","anisotropic pressure","adiabatic growth","Schwarzschild metric","cuspy halo profile"],"falsifier":"Evaluate the conservation residual C^r and the tangential Einstein equation using the paper's numerical g(r), m(r), and T^μ_ν; a residual comparable to the claimed deviations would show the metric is an artifact of the truncated system. An iterative re-computation of the distribution function and T in the backreacted metric, repeated to convergence, would also test whether the factor-2.5 enhancement persists.","tokens_in":15948,"feed_emoji":"🕳️","tokens_out":11774,"duration_ms":95721,"temperature":0.7,"pith_summary":"The paper sets out to show that a dark-matter spike around a black hole gravitates through its full dynamical content—the kinetic energy of bound orbits and the anisotropic pressure they generate—not just through its rest-mass density. Within a cluster model of collisionless particles, the authors compute the spike's complete energy-momentum tensor after adiabatic black-hole growth, using a cuspy halo profile and Milky Way parameters, and find that near the spike the kinetic term adds roughly 50% to the energy density relative to the rest-mass term. Feeding this tensor into the Einstein equations as a fixed source, the metric functions g(r) and m(r) deviate from Schwarzschild by about 2.5 times more than in earlier mass-density-only treatments near the spike's peak. If this is right, earlier spike-backreaction models have understated how strongly the spike curves spacetime. The computation is a step toward a self-consistent relativistic description of black holes embedded in dark-matter halos.","feed_headline":"Dark matter spike bends the black-hole metric 2.5x more than thought","feed_subtitle":"Kinetic energy adds ~50% to the spike's density and bends spacetime around the Milky Way's black hole.","key_machinery":"The central object is the Einstein-cluster energy-momentum tensor of the spike, T^μ_ν = ∫ u^μ u_ν [M f^(4)(x,p)] √(-g) d^4p, built from the phase-space distribution f(E,L²)=f^(halo)(ε(E,L²)), with the orbital energy mapped through the adiabatic-invariant relation between the halo orbit and a Schwarzschild geodesic. It yields explicit integrals for the energy density, radial pressure, and tangential pressure. The backreaction metric is obtained by solving the tt and rr Einstein equations for g(r) and m(r), with this tensor as a fixed source and Schwarzschild boundary conditions at the spike's inner edge.","core_discovery":"Within the Einstein-cluster framework, the paper derives the complete diagonal energy-momentum tensor of the spike—energy density, radial pressure, tangential pressure—from a distribution function conserved through adiabatic black-hole growth, evaluated in a fixed Schwarzschild background. Near the spike, the full energy density is about 1.5 times the rest-mass density; the radial pressure is two orders of magnitude smaller but nonzero, becoming comparable to the tangential pressure at large radius; all standard energy conditions hold. Solving the tt and rr Einstein equations with this tensor as a fixed source yields metric functions whose deviations from Schwarzschild exceed those of a mass","pith_inferences":["The paper attributes the enhancement mainly to kinetic energy but does not separately quantify how much comes from the kinetic term versus the non-zero radial pressure; a run that switches on each term independently would isolate the two effects.","Because the conservation equation is knowingly not satisfied in the perturbed metric and no residual is reported, the solved (g,m) may not be the exact metric sourced by the claimed tensor; computing the residual of the unenforced equations is a direct check the paper leaves open.","Applying the same construction to other halo profiles (a different cuspy model or a cored profile) would test whether the ~50% kinetic share and the ~2.5 factor are robust or specific to the chosen cuspy profile.","If a fully self-consistent calculation confirms the stronger deviation, dark-matter spikes would imprint larger-than-expected signatures in gravitational-wave phase evolution and in photon orbits around the Galactic center black hole."],"forward_implications":["Near the peak of the energy density, the backreacted metric functions deviate from Schwarzschild by about 2.5 times more than a mass-density-only treatment, so earlier estimates understate the spike's spacetime effect.","The kinetic contribution makes up about half the peak energy density, so the spike's total gravitational source is not well approximated by its rest-mass density alone.","The radial pressure is nonzero in the spike and becomes comparable to the tangential pressure at large radius; setting it to zero is not a globally valid simplification.","The derived energy-momentum tensor satisfies the null, weak, strong, and dominant energy conditions, so it is physically admissible as a source.","The mass-function deviation m(r)-m_Schwarzschild is larger than the g(r) deviation, indicating the main backreaction is a redistribution of the effective mass profile seen by orbiting particles."],"fun_headline_variants":["Black hole metric bends 2.5x more with full dark matter spike","Dark matter spike kinetic energy boosts density 50%, warping black hole","Full spike tensor in Einstein cluster yields 2.5x metric deviation","Kinetic term adds 50% to spike density, bends black hole metric 2.5x"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The computation solves only the tt and rr Einstein equations, using a source evaluated in the fixed Schwarzschild background, and the paper acknowledges that the conservation equation is not satisfied in the perturbed metric; if the conservation residual is not negligible, the reported metric deviations are not the metric sourced by the claimed tensor.","fun_headline_variants_meta":{"raw":{"variants":["Black hole metric bends 2.5x more with full dark matter spike","Dark matter spike kinetic energy boosts density 50%, warping black hole","Full spike tensor in Einstein cluster yields 2.5x metric deviation","Kinetic term adds 50% to spike density, bends black hole metric 2.5x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000806,"raw_usage":{"total_tokens":3366,"prompt_tokens":723,"completion_tokens":2643,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":467,"completion_tokens_details":{"reasoning_tokens":2566}},"tokens_in":467,"tokens_out":2643,"duration_ms":18360,"temperature":1.0,"reasoning_tokens":2566,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T21:59:48.843994+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the conservation residual C^r and the tangential Einstein equation using the paper's numerical g(r), m(r), and T^μ_ν; a residual comparable to the claimed deviations would show the metric is an artifact of the truncated system. An iterative re-computation of the distribution function and T in the backreacted metric, repeated to convergence, would also test whether the factor-2.5 enhancement persists.","supporting_citations":[],"review_version":1}