{"id":"1b805e2a-cf19-48bf-9357-eda4ecb1f86b","arxiv_id":"2505.04697","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Spin of the central black hole raises the dark matter spike density and the gravitational-wave dephasing it produces, so LISA can detect the dark matter effect more easily.","lead":"This paper studies gravitational waves from a small black hole spiraling into a larger spinning black hole that sits inside a rotating dark matter spike. It finds that the black hole's spin strengthens the dark matter's effect on the waves, improving the chance that the LISA space observatory will detect it.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The spin-enhancement claim rests on an isotropic DF model; the Ferrer spike's rotating, toroidal velocity distribution could change the torque and reverse the prograde/retrograde ordering.","rationale":"The reader's conditional verdict is appropriate. The paper has real strengths: the use of Teukolsky fluxes and FEW for fully relativistic trajectories is a solid pipeline, and the Section IV PN estimate that metric backreaction is negligible for 10^5-10^6 M_sun primaries is a useful, independently checkable check. The qualitative direction that a denser spike produces more dephasing is robust. However, all quantitative detectability claims flow through the dynamical friction force, and the DF model is at odds with the very environment the paper constructs. The Ferrer spike is anisotropic and rotating; Eq. (12) is isotropic, ignores the relative velocity between the secondary and the DM, and evaluates only the equatorial density. This is not a small correction: it directly controls the torque on prograde versus retrograde orbits, which is the mechanism behind the paper's spin-enhancement conclusion. The authors explicitly defer the gravitational Magnus effect, but that deferral leaves the central claim unsecured. A phase-space-resolved DF computation is needed before the quantitative LISA detectability statements can be trusted. The conclusion's reference to a chi = 0.8 result absent from the figures is a minor internal inconsistency, but not the main issue. Reproducibility concerns about missing fit coefficients and code further justify conditionality, but the physics assumption is the load-bearing point.","tokens_in":25801,"tokens_out":4106,"duration_ms":46549,"concrete_test":"Compute the DF force for a circular equatorial orbit at, say, r = 10 M_BH in a chi = 0.7 Ferrer spike using the full phase-space distribution f(E,C,Lz) from Sec. II, via the local integral F_DF = 4 pi G^2 m^2 lnLambda integral d^3 v_DM f(v_DM) (v_sec - v_DM) / |v_sec - v_DM|^3 (with a consistent relativistic generalization). Compare the azimuthal component and the perpendicular components with Eq. (12) for both prograde and retrograde orbits. If the torque differs by more than ~50% from the isotropic formula, recompute the dephasing curves in Fig. 6; if the prograde/retrograde ordering changes, the spin-enhancement claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is that dynamical friction is captured by the isotropic Chandrasekhar formula, Eq. (12), with lnLambda ~ 3, force anti-parallel to the secondary's velocity, and density evaluated on the equatorial plane. The Ferrer et al. spike [51] is toroidal and its phase-space distribution f(E,C,Lz) carries a net rotational velocity. For an equatorial circular orbit, the secondary's velocity and the local DM bulk velocity are not generally aligned: prograde motion is closer to co-rotation while retrograde motion is counter-rotation. Eq. (12) uses the secondary velocity v rather than the relative velocity v - v_DM and ignores force components perpendicular to v. This enters every dephasing, mismatch, and SNR result through Eqs. (11), (14), and (15). A non-aligned drag, or even an aligned drag with relative velocity, changes the torque Ldot_DF and therefore the spin dependence shown in Figs. 6-11. The authors acknowledge the gravitational Magnus effect and defer it to future work, but the headline claim that spin improves detectability and that high-spin prograde orbits are optimal is not secured while this is the only spin-dependent environmental input. The absence of released code and fit coefficients compounds the problem, but the physics assumption is the core issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models extreme mass ratio inspirals in a rotating dark matter spike around a Kerr black hole. It uses the Ferrer et al. relativistic spike profile, fitted by the polynomial scaling relation in Eq. (10), and incorporates dynamical friction through the isotropic Chandrasekhar formula, Eq. (12), into the FEW trajectory and waveform code with Teukolsky fluxes. The authors compute dephasing, mismatch, and SNR for varying spin, halo mass, and scale radius, and also give a 1PN estimate of the metric backreaction from the spike, concluding that the spin of the primary enhances the detectability of DM effects with LISA and that the backreaction is negligible for primary masses ≲ 10^7 Msun.","tokens_in":26025,"tokens_out":6202,"duration_ms":66514,"significance":"If the central claim holds, LISA parameter estimation for EMRIs must account for rotating DM environments, and this work provides one of the first fully relativistic treatments of the spike geometry combined with Teukolsky fluxes. The forward-model structure is a strength: the spike profile and DF formula are external inputs and the dephasing, mismatch, and SNR are outputs, so there is no circularity in the main computation. The paper also gives explicit validity bounds from halo feedback and from the probe-limit approximation, which is commendable. However, the absence of the fit coefficients and code limits reproducibility, and the spin-dependence of the standard DF model is a load-bearing assumption that needs quantitative scrutiny.","major_comments":[{"comment":"The dynamical friction force is modeled with the isotropic Chandrasekhar formula, with the force anti-parallel to the secondary's velocity, ln Lambda ~ 3, and the density evaluated on the equatorial plane. The Ferrer et al. spike [51] is toroidal and its phase-space distribution carries a net rotational velocity, so the relative velocity v - v_DM should enter the drag, and a non-aligned component (the gravitational Magnus effect) may also exist. Because the spin-dependence of the dephasing and mismatch is the central result of the paper, the authors should estimate the magnitude of the correction from using the relative velocity, for example by computing the bulk velocity of the spike, or explicitly soften the claim that the spin of the primary improves detection prospects. The acknowledgment in Section III that the Magnus effect is deferred does not address the aligned relative-velocity correction, which changes the torque even for equatorial circular orbits.","section":"III, Eq. (12)"},{"comment":"The spike density used in all subsequent results is a polynomial fit whose coefficients A_i, B_j and orders n, m are not given, and no code or data release is mentioned. Since the dephasing, mismatch, and SNR results scale with the spike density, the quantitative claims are not reproducible without these coefficients or a supplementary data file. Please include a table of the fit coefficients for each spin value and state the fitting range and error as a function of radius, or make the fitting code available.","section":"II, Eqs. (9)-(10)"},{"comment":"The paper uses both SNR and mismatch to discuss observability, but these measure different things: SNR in Fig. 8 is the loudness of the non-vacuum waveform, not the distinguishability of the DM effect, while mismatch in Figs. 9-11 is the appropriate detectability criterion. The apparent tension that retrograde orbits reach SNR = 20 earlier but have lower mismatch than prograde orbits should be addressed explicitly in the text, so that readers do not infer that SNR supports the 'high-spin prograde optimal' conclusion, which is based on mismatch.","section":"V, Figs. 8-9"}],"minor_comments":[{"comment":"Please define r_pk in the fitting function and clarify whether the Heaviside function H(r - r_mb) creates a discontinuity at the matching point; also report the fitting error as a function of radius rather than only the global 0.1% claim.","section":"II, Eq. (10)"},{"comment":"The abbreviation 'AAK' is used without definition; please spell out 'Augmented Analytic Kludge' at first use and give a reference.","section":"III, Eq. (13) and surrounding text"},{"comment":"The axis label in the left panel for the (10^6 + 50) Msun system appears garbled as '104 7'; please correct the typesetting.","section":"V, Fig. 6"},{"comment":"In the definition of N_cycle, the symbol F is used for the GW frequency but the notation is ambiguous because F appears both in the numerator and in the denominator; please define F_dot = dF/dt and write the integral accordingly.","section":"III, Eq. (16)"},{"comment":"The Poisson equations for the environmental potentials involve the mass current J_H^i; please state explicitly how J_H^i is obtained from the Boyer-Lindquist currents in Eq. (8) via the coordinate transformation, since this is needed to reproduce the metric backreaction estimate.","section":"IV, Eq. (24)"}],"recommendation":"major_revision","confidential_remarks":"The paper's main quantitative claim rests on the standard isotropic Chandrasekhar DF formula, which the authors themselves flag as a simplification. This is a fixable issue if they provide a quantitative estimate of the relative-velocity correction or reframe the claim. The missing fit coefficients and lack of code release are a reproducibility concern that should be addressed before publication. The manuscript fits the journal's scope well, and the forward-model nature of the calculation is a positive feature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is the first to put the Ferrer et al. rotating Kerr spike profile together with fully relativistic Teukolsky fluxes and FastEMRIWaveforms, and it produces concrete spin-dependent predictions for dephasing, mismatch, and SNR. That combination is genuinely new and worth knowing about. The qualitative message—that BH spin changes how strongly a DM spike imprints on EMRI waveforms—is plausible and probably right. The PN estimate of the environmental backreaction on the metric is a thoughtful addition, and the authors are honest about the probe limit and halo-feedback bounds on their mass range.\n\nThe soft spots are real but not fatal. The biggest one is the dynamical friction model. Eq. (12) is the isotropic Chandrasekhar formula with the force anti-parallel to the secondary's velocity, evaluated with the equatorial density. The Ferrer spike is toroidal and has a net rotational velocity; for prograde versus retrograde orbits the relative velocity between the secondary and the DM bulk differs. Using v instead of v - v_DM, and dropping perpendicular force components, can change the torque and therefore the spin ordering in Figs. 6–11. The authors acknowledge the Magnus effect and defer it, but that means the headline claim that high-spin prograde orbits are the best targets is not secured by the present calculation. It is a modeling assumption, not a solved result.\n\nRelatedly, the paper doesn't give the polynomial fit coefficients in Eq. (10) and releases no code or data. For a result whose quantitative claims depend on those fits, that makes independent verification harder than it should be. The fit error is quoted as <0.1%, but coefficients should be in an appendix or ancillary file. There's also a minor internal inconsistency: the conclusion cites a spin-0.8 system for the SNR=20 timing, but the figures show ±0.3 and ±0.7 only. Probably a typo or an unshown run, but it should be fixed.\n\nThe mismatch and SNR statements use standard thresholds (M>0.03, SNR>20) without a full parameter-estimation study. The authors say PE is future work, which is fine, but the detectability language in the abstract and conclusion is a bit stronger than the evidence supports.\n\nAll that said, the paper is a serious, useful step. The central calculation is plausible, the limitations are mostly acknowledged, and the spin-dependent maps will be a reference point for the LISA EMRI literature. It deserves a serious referee; I would send it to review, with the DF modeling and the missing fit coefficients as the main requests for revision.","headline":"First rotating-spike EMRI study with Teukolsky fluxes; useful and mostly sound, but the spin-enhancement claim rests on an isotropic DF model that the toroidal spike calls into question.","tokens_in":26612,"tokens_out":4021,"would_cite":true,"duration_ms":36110,"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":"Spinning black holes make dark matter spikes easier for LISA to detect, and future EMRI parameter estimation must include them.","keywords":["extreme mass ratio inspirals","dark matter spikes","dynamical friction","Kerr black holes","LISA","gravitational wave dephasing","waveform mismatch","signal-to-noise ratio"],"falsifier":"Numerically simulate a small black hole moving through the rotating spike and measure the drag component perpendicular to its velocity. If that perpendicular component is a significant fraction of the parallel drag, the torque balance used in the paper is incomplete and the prograde-versus-retrograde dephasing hierarchy would need revision.","tokens_in":25561,"feed_emoji":"🕳️","tokens_out":8044,"duration_ms":75626,"temperature":0.7,"pith_summary":"The paper argues that when the central black hole spins, the dark matter spike around it is compressed and densified near the equator, and this rotating environment leaves a stronger imprint on extreme mass ratio inspiral waveforms than a non-rotating spike. The central claim is that the spin of the primary black hole improves the chances that LISA will detect dark matter spikes, and that any parameter-estimation study of EMRIs must include the rotating spike or risk biased results. The support comes from dephasing, mismatch, and signal-to-noise calculations for prograde and retrograde orbits around Kerr black holes of $10^5$–$10^6\\,M_\\odot$.","feed_headline":"Spin makes dark matter spikes detectable by LISA","feed_subtitle":"Rotating dark matter spikes change EMRI waveforms enough that LISA can see them, so future parameter studies must include spin.","key_machinery":"The central object is the rotating dark matter spike: a toroidal density distribution obtained by adiabatically growing a Kerr black hole inside a Hernquist halo, with particles conserving their action integrals during the growth. The argument runs through this density profile—parametrized by the fitting formula and scaling relations in Eq. (9) and Eq. (10)—which feeds the Chandrasekhar dynamical friction force of Eq. (12). That force is added to fully relativistic Teukolsky-based gravitational-wave fluxes in an energy and angular-momentum balance evolution (Eqs. (11) and (15)), and the resulting waveforms are compared to vacuum waveforms through dephasing, mismatch, and SNR.","core_discovery":"The paper's central discovery is that rotation changes the dark matter environment itself, not just the orbit: higher spin produces a denser, more centrally concentrated equatorial spike, so the dynamical friction drag on the secondary grows. In inspirals followed all the way to the ISCO, prograde orbits experience larger environmental dephasing because their ISCO is closer to the black hole, letting them accumulate more cycles in the dense inner spike; retrograde orbits, by contrast, reach higher SNR earlier in fixed-time observations. The mismatch between vacuum and dark-matter waveforms exceeds the standard $0.03$ distinguishability threshold for a wider range of halo parameters when spin is included, and this is taken as evidence that LISA detection prospects improve and that rotation must be included in future parameter estimation.","pith_inferences":["If the gravitational Magnus effect—the non-aligned component of dynamical friction—is non-negligible in the toroidal, anisotropic spike, the torque balance and the prograde/retrograde dephasing hierarchy would shift; the authors explicitly leave this to future work.","A direct testable extension is a full Bayesian parameter-estimation study that injects rotating-spike waveforms and recovers with vacuum and non-rotating templates; the expected outcome is a measurable bias in spin and environment parameters.","The prograde/retrograde asymmetry is driven mostly by the ISCO radius, so systems with the highest spins and small secondary masses should show the strongest environmental mismatch; scanning that grid would sharpen the detectability map."],"forward_implications":["LISA parameter-estimation studies that approximate the environment with a static, spherical spike will produce biased estimates of the black hole spin and environment parameters, because rotation changes both the density profile and the orbital evolution.","High-spin, prograde EMRIs are the best targets for detecting dark matter spikes, while retrograde inspirals can become louder sooner in fixed-time observations.","The probe-limit approximation—ignoring the spike's back-reaction on the metric—is safe for primary masses below about $10^7\\,M_\\odot$, so the main dephasing and mismatch results are not contaminated by metric changes for the systems considered.","Including spin lowers the halo mass or increases the scale radius at which a dark matter environment becomes distinguishable from vacuum, so LISA could probe more diffuse halos than non-rotating models suggested."],"supporting_citations":[{"why":"Constructs the fully relativistic rotating dark matter spike density around a Kerr black hole from adiabatic growth; this density is the environmental input to the dynamical friction force.","marker":"[51]"},{"why":"Provides the classical Chandrasekhar dynamical friction formula used to add environmental drag to the inspiral.","marker":"[53]"},{"why":"Supplies the dephasing and mismatch methodology, plus the density scaling relations used to fit spike profiles for waveform studies.","marker":"[15]"},{"why":"Provides the waveform-generation machinery used to evolve inspirals and produce gravitational waveforms for the comparisons.","marker":"[76–79]"},{"why":"Supplies fully relativistic Kerr gravitational-wave fluxes via Teukolsky solutions used in the energy and angular-momentum balance.","marker":"[87]"},{"why":"Supplies the LISA power spectral density used in the mismatch and signal-to-noise ratio computations.","marker":"[91]"},{"why":"Establishes the relativistic spike profile for a non-rotating black hole, which the rotating case generalizes.","marker":"[50]"},{"why":"Establishes the adiabatic spike formation picture on which the rotating spike model builds.","marker":"[49]"}],"fun_headline_variants":["Black hole spin boosts dark matter spike detectability","Spin enhances dark matter spike signals for LISA","Rotating spikes become visible to LISA via spin","Spin-dependent dark matter spikes improve LISA detection"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the drag force is always opposite to the secondary's velocity and has the size given by the standard Chandrasekhar formula with $\\ln\\Lambda\\sim3$, evaluated on the equatorial spike density; if the rotating spike produces a sideways drag or a different effective Coulomb logarithm, the spin-dependent dephasing results would shift.","fun_headline_variants_meta":{"raw":{"variants":["Black hole spin boosts dark matter spike detectability","Spin enhances dark matter spike signals for LISA","Rotating spikes become visible to LISA via spin","Spin-dependent dark matter spikes improve LISA detection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000397,"raw_usage":{"total_tokens":2048,"prompt_tokens":884,"completion_tokens":1164,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":500,"completion_tokens_details":{"reasoning_tokens":1104}},"tokens_in":500,"tokens_out":1164,"duration_ms":8505,"temperature":1.0,"reasoning_tokens":1104,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:23:48.588955+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically simulate a small black hole moving through the rotating spike and measure the drag component perpendicular to its velocity. If that perpendicular component is a significant fraction of the parallel drag, the torque balance used in the paper is incomplete and the prograde-versus-retrograde dephasing hierarchy would need revision.","supporting_citations":[],"review_version":1}