{"id":"77a79010-b096-4c23-961d-c2ab2d736353","arxiv_id":"2505.10036","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Repeated tidal resonances from an inspiralling EMRI secondary shift a dark matter particle's semi-major axis by less than one percent, leaving the DM spike nearly intact in the non-overlapping-orbit regime.","lead":"Dark matter particles orbiting a supermassive black hole are repeatedly tugged by the gravity of a small black hole spiraling inward, and this paper calculates how much those repeated tugs push the dark matter around. It finds the cumulative effect is under one percent in the cases studied, meaning this resonant scattering channel is too weak to reshape the dark matter spike that LISA will probe.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Sec. V 'regardless of initial orbits' claim is not supported in the orbital-overlap regime, because the Teukolsky reconstruction (Eq. 72) and all numerical runs stop before overlap; the authors themselves say the effect may be enhanced there.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the central conclusion is established only for DM orbits that remain entirely inside the EMRI's radial range, while the paper's own Sec. V states that the resonance effect may be enhanced when the orbits overlap. This is not a hidden internal inconsistency; it is an explicit scope limitation, but it directly undercuts the unqualified phrasing of the abstract and Sec. V ('regardless of the initial orbits'). A DM particle with a larger initial semilatus rectum than the pDM = 25M used in Tables I-II would encounter the overlap region during the EMRI inspiral, and the Teukolsky reconstruction used to compute the interaction Hamiltonian is not valid there. The exponential decay of |H| with resonance order, the paper's central explanatory mechanism, is demonstrated only on the non-overlap branch, so extrapolating it to the overlap branch is unsupported. The fast-resonance phase-maximization prescription is conservative and therefore not a concern for the negative result; the analytic derivation and the non-overlap numerics are serious and credible. Because the authors themselves flag the overlap limitation, the appropriate verdict remains CONDITIONAL: the paper should be accepted only with the explicit restriction that the small-effect conclusion applies to the non-overlapping regime until the overlap contribution is computed or bounded.","tokens_in":15145,"tokens_out":14568,"duration_ms":161967,"concrete_test":"Extend the Sec. III A reconstruction by adding the r' < r branch of the Teukolsky Green's function (the counterpart of Eq. (72), with the homogeneous solutions assigned to source and field point swapped) and rerun the Fig. 4 evolution with initial DM semilatus recta pDM = 30M, 40M, and 50M, keeping eDM = 0.2, xDM = 0.7 and the same EMRI parameters. Continue the integration past the point where r_p,EMRI < r_a,DM instead of stopping. If the cumulative fractional changes in pDM and eDM from the overlapping segment stay below about 1%, the Sec. V claim survives; if they reach several percent or more, the general 'regardless of initial orbits' conclusion must be restricted to the non-overlapping regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing gap is the exclusion of the orbital-overlap regime from the central claim. The interaction Hamiltonian in Sec. III A is derived from the Teukolsky Green's function branch (Eq. 72) that is valid only when the source radius r' exceeds the field-point radius r, and every numerical run in Sec. IV is stopped as soon as the two radial ranges overlap (Fig. 4 note; Sec. V). The authors state explicitly in Sec. V that in this excluded regime 'the DM particle and the EMRI secondary will get closer to each other in this case, and thus the effect of resonance will be enhanced.' Consequently, the Sec. V conclusion that the cumulative effect 'remains very small regardless of the initial orbits of the DM particles' is not established for exactly the orbital configurations in which the resonance interaction is expected to be strongest. Tables I-II and Fig. 6 only sample pDM = 25M with eDM ≤ 0.4, keeping the DM apocenter below the EMRI pericenter through the entire inspiral. DM particles with larger pDM (for example, pDM around 50M for the Fig. 4 EMRI parameters) would enter the overlap region before the EMRI reaches the relevant part of its inspiral; these are the particles missing from the claimed result. The exponential-decay argument based on Fig. 5 is computed on the non-overlap branch and cannot be assumed to carry over, because the Green's function, and hence the Fourier coefficients H_{n,N}, differ in the overlap branch. Thus the headline claim is conditional on a regime that may contain the dominant effect.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript investigates whether an extreme-mass-ratio inspiral (EMRI) in the galactic center can deplete a dark-matter (DM) spike through tidal resonances between the secondary and individual DM particles on bound Kerr geodesics. The setup is a restricted three-body problem: the EMRI evolves by gravitational-wave backreaction (quadrupole formula), while each DM particle is treated as a test particle whose orbital elements change discretely at each resonance satisfying n_a ω_a = N_a Ω_a. The interaction Hamiltonian is derived from first principles using Teukolsky metric reconstruction (Eqs. (59)-(88)), Fourier-decomposed in the action-angle variables of both orbits, and the cumulative effect of O(10^3) resonances per EMRI event is computed numerically with Black Hole Perturbation Toolkit modes. For a benchmark run (central mass 10^6 M_sun, η=10^-4, spin 0.9, EMRI from p=126 to p=72, DM at p=25), the DM semi-major axis changes by about 0.3%, and a scan over DM eccentricity and inclination (Tables I-II, Fig. 6) gives changes below about 1% in semi-major axis and up to about 4.5% in eccentricity. The authors attribute the smallness to an exponential decay of the resonance Hamiltonian with resonance order K (Fig. 5). All numerical runs are stopped before the DM and EMRI radial ranges overlap, and the authors state in Sec. V that the overlap regime may enhance the effect.","tokens_in":15420,"tokens_out":18875,"duration_ms":187939,"significance":"This is a potentially useful negative result for the DM-spike/EMRI literature: within the non-overlapping-orbit regime, the tidal-resonance feedback channel does not significantly modify DM orbits, so the spike survives this specific channel and EMRI-phase signatures of a spike are not erased by it. The paper's strengths are concrete: the interaction Hamiltonian is derived rather than fitted; the exponential suppression of high-order resonances is extracted from numerical evaluation of the Fourier coefficients rather than assumed as an ansatz; the authors are transparent about the six-orders-of-magnitude scatter at fixed K and about the excluded overlap regime; and the numerical pipeline (public Teukolsky toolkit, datasets promised at an archived repository) is reproducible in principle. The significance is conditional on the domain of validity: the claimed universality 'regardless of initial orbits' currently covers only orbits that never cross the EMRI's radial range, a subset for which the interaction is weakest; quantifying the excluded regime is the main open task.","major_comments":[{"comment":"The paper's headline claim, that the cumulative resonant effect on DM orbits 'remains very small regardless of the initial orbits of the DM particles' (Sec. V), is established only for the non-overlap regime, which is precisely the regime of weakest interaction. The reconstructed metric (81) is built from the retarded Green's function branch (72) that holds for source radius r' larger than the field-point radius r, and every numerical run is terminated as soon as the DM and EMRI radial ranges overlap (note to Fig. 4; Sec. IV). Tables I and II and Fig. 6 fix pDM = 25M, so for the Fig. 4 EMRI (pEMRI from 126M down to 72M) the DM apocenter remains below the EMRI pericenter throughout and the excluded regime is never probed. For example, a DM particle with pDM = 50M and eDM = 0.4 has apocenter ≈ 83M and would begin overlapping the EMRI's radial range once pEMRI drops below ≈ 100M, i.e., during the studied inspiral window; these are the configurations for which the authors themselves state that the resonance effect 'will be enhanced' (Sec. V). As it stands, the conclusion is conditional on the overlap-free assumption, and the spike-survival statement in the abstract should be qualified accordingly, or supplemented by a quantitative estimate of the overlap-regime contribution (for instance, using the complementary branch of the Green's function or a local treatment of near-coincidence encounters).","section":"Sec. V; Eq. (81); Fig. 4 note"},{"comment":"The claim that the smallness holds 'regardless of the initial orbits' is an extrapolation from a narrow parameter slice. The numerical scans hold pDM fixed at 25M, vary eDM only over {0.2, 0.3, 0.4} and xDM over the values shown in Tables I-II and Fig. 6, and do not vary the EMRI mass ratio or spin; the sentence in Sec. IV that 'systematic parameter variation studies reveal negligible quantitative differences' refers to EMRI parameters and is not substantiated by any figure or table in the manuscript. The authors should state explicitly the region of parameter space over which the conclusion is claimed (e.g., the set of (pDM, eDM, xDM) such that the DM radial range never intersects the EMRI's during the relevant inspiral), and ideally add at least one run with a larger pDM or a different mass ratio to demonstrate that the exponential-decay mechanism persists.","section":"Sec. IV; Tables I-II; Fig. 6"},{"comment":"The truncation at Kmax = 15 is asserted to be convergent without supporting evidence. Appendix A states that 'convergence of the results has been verified,' but no figure, table, or quantitative criterion is given (e.g., the change in the cumulative Δa/a or Δe when Kmax is increased from 12 to 15, or the largest fractional contribution from the highest-K retained resonances). This matters because Figs. 8-9 show that the decay of the interaction Hamiltonian with K becomes slower as the inclination increases, and the high-inclination points in Fig. 6 (xDM down to 0.1) are part of the headline result; an undiagnosed truncation bias in exactly those cases could affect the stated smallness of the cumulative effect. A short convergence table or a statement of the error criterion is needed.","section":"Appendix A; Sec. IV"}],"minor_comments":[{"comment":"The prescription that selects 'the phase to maximize the contribution from the resonance' for fast crossings should be labeled as an upper-envelope estimate for that channel, not a phase-averaged value; since the resulting cumulative change is still small, the conclusion is conservative under this choice, but the manuscript should say so explicitly so that Fig. 4 and Tables I-II are not read as typical trajectories.","section":"Sec. IV (fast-crossing phase)"},{"comment":"In the derivation following Eqs. (47)-(53), the jump is Δj_a = n_a Θ with Θ = 8√H/(π√k) and k = G_ab n_a n_b, so the denominator in Eq. (54) should be √k (or a defined norm |n|_G), not the Euclidean |n|; as printed, the equation is inconsistent with its own derivation.","section":"Eq. (54)"},{"comment":"The 'exponential decay' of the Hamiltonian with K is inferred from the envelope of a scatter of points spanning six orders of magnitude at fixed K; quoting an e-folding scale per unit K and the K range over which the envelope is exponential would make the claim quantitative and testable, and would directly support the Kmax = 15 truncation.","section":"Fig. 5"},{"comment":"The motivation for neglecting direct scattering is derived for η ≲ 10^-6, while the numerical study uses η = 10^-4; for η = 10^-4 the crossover radius in Eq. (8) is about 478M, so the sentence following Eq. (8) should state explicitly that the conclusion applies to the LISA-band radii (≲ O(100)M) for the adopted mass ratio.","section":"Sec. II A"},{"comment":"First sentence of the concluding paragraph: 'the effect on the overall effect on the DM orbital evolution' is redundant and should be reworded.","section":"Sec. V"},{"comment":"The statement that 'the results for other cases are similar' is presented without supporting evidence; a supplementary panel or a quantitative statement on the robustness of the resonance count would help the reader assess the O(10^3) estimate.","section":"Sec. II C (Fig. 2 discussion)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on the authors' own earlier work (Refs. [24], [28], [31]) for the action-angle formalism and the interaction Hamiltonian; an editor may wish to confirm with the authors that the Teukolsky-reconstruction derivation in Sec. III adds to Ref. [21] (Silva-Hirata) rather than duplicating it, and that Ref. [24] contains the leading-order version of the jump formula used here. Substantively, my recommendation hinges on scope: the abstract's unqualified claim exceeds the computed domain, but the authors' own disclosure of the overlap limitation indicates awareness of the gap; a solid revision would add either a quantitative estimate of the overlap regime or an explicit scope restriction in the abstract, together with the missing convergence evidence for Kmax = 15."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Kakehi et al. compute, for the first time, the cumulative effect of multiple tidal resonances on a dark-matter particle orbiting a Kerr black hole as an EMRI inspirals. The machinery is serious: they derive the averaged interaction Hamiltonian from metric reconstruction via the Teukolsky formalism, use the Black Hole Perturbation Toolkit for mode functions, and evolve the DM orbit through thousands of resonances. The main finding—that the resonance Hamiltonian decays exponentially with resonance order K, so only a handful of low-order resonances matter—looks robust in the regime they can treat, and the resulting total change in semi-major axis stays below about 1% for their sampled parameters. That is a useful boundary condition for DM-spike and EMRI modeling. The single-resonance Hamiltonian overlaps with Silva and Hirata (2022), but the cumulative multi-resonance application to DM spike evolution is new.\n\nThe paper is honest about its main limitation: the Teukolsky Green's function branch they use (Eq. 72) is valid only when the EMRI secondary lies outside the DM particle's radial range, and every numerical run is stopped as soon as the two orbits overlap. The Sec. V sentence saying the effect 'remains very small regardless of the initial orbits' is therefore too strong. The authors themselves note in the very next paragraph that in the overlapping regime the interaction should be enhanced, because the perturber and the particle get closer. So the correct statement is: within the non-overlap regime, the cumulative resonance effect is small. Whether overlap resonances can disrupt the spike remains open.\n\nOther weaknesses are minor. The truncation at K_max=15 is justified only by an assertion of convergence; a convergence plot or error estimate would help. The fast-crossing phase is chosen to maximize the effect rather than sampled from a distribution—if anything this makes the small-effect conclusion conservative, but it would be nice to see the phase dependence. The datasets are promised but not yet archived.\n\nThe derivation is standard in structure and the numerical pipeline plausible; neither the reader's skepticism nor mine found a load-bearing contradiction. This is a conditional result with an explicit boundary, and it deserves a serious referee. I would send it to review, with a request to soften the 'regardless' claim and to state clearly which regimes are covered.","headline":"Solid derivation and numerical pipeline showing cumulative tidal resonances barely affect DM orbits before the EMRI and DM radial ranges overlap, but the 'regardless of initial orbits' claim overreaches into the overlap regime where the effect may be enhanced.","tokens_in":15984,"tokens_out":3449,"would_cite":true,"duration_ms":32623,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Tidal resonances from an inspiraling compact object barely move dark-matter particles near a supermassive black hole.","keywords":["dark matter spike","extreme mass ratio inspiral","tidal resonance","Kerr geodesics","Teukolsky equation","halo feedback","gravitational wave backreaction","restricted three-body problem"],"falsifier":"Run the same cumulative-resonance evolution with a metric-reconstruction scheme that is valid when the EMRI secondary and the DM particle orbits overlap, and check whether the total semi-major-axis change exceeds the ~1% level for the same fiducial parameters; if it does, the paper's conclusion that the DM spike survives this channel fails.","tokens_in":14899,"feed_emoji":"🕳️","tokens_out":7411,"duration_ms":70016,"temperature":0.7,"pith_summary":"Dark matter around a supermassive black hole is thought to form a dense spike, and an extreme-mass-ratio inspiral (EMRI) plunging through it could shake that spike through tidal resonances: whenever a combination of the DM particle's orbital frequencies matches one of the EMRI's, the small gravitational pull can accumulate into a lasting change of the DM orbit. The paper derives the full interaction Hamiltonian for these resonances in Kerr spacetime and evolves DM particles through the thousands of resonances they encounter as the EMRI spirals inward. It finds that, for orbits that do not overlap the EMRI's radial range, the cumulative effect is tiny—typically under one percent change in semi-major axis and a few percent in eccentricity—because the resonance amplitude decays exponentially with the resonance order, so only a few low-order resonances matter. A sympathetic reader would take this as evidence that this particular backreaction channel does not destroy the DM spike, leaving the spike available as a gravitational-wave probe of dark matter.","feed_headline":"Tidal resonances leave dark-matter spikes nearly untouched","feed_subtitle":"Thousands of resonant kicks shift a typical dark-matter orbit by under one percent around a supermassive black hole.","key_machinery":"The central object is the averaged interaction Hamiltonian between the EMRI secondary and a DM test particle, expressed in action-angle variables of Kerr geodesics: ⟨H_int⟩(χ) = Σ_{(n,N)∈Res} H_{n,N} $e^{{iχ}}$, with χ = n_a q^a + N_a Q^a and resonance condition n_a ω̂^a + N_a Ω̂^a = 0. Around a single resonance the resonant angle obeys a pendulum equation d²χ/dτ² = -α sin χ, so the orbital jump is Δj_a = (8√H/(π|n|)) n_a, proportional to √H. The amplitudes H_{n,N} are obtained from metric reconstruction of the Teukolsky equation, and the load-bearing numerical observation is that they decay exponentially with K = Σ(|n_a|+|N_a|), which makes higher-order resonances harmless and concentrates the effect in a few low-order events.","core_discovery":"On its own terms, the paper's claim is that the repeated tidal resonances between an EMRI secondary and a DM particle—each event shifting the particle's actions by an amount proportional to the square root of the interaction Hamiltonian—do not accumulate into a significant orbital change. Numerically, for a central black hole of mass $10^{6}$ solar masses, mass ratio $10^{-4}$, spin 0.9, and a range of DM eccentricities and inclinations, the semi-major axis changes by at most roughly 0.8% and the eccentricity decreases by a few percent over the inspiral, with 2612 resonances crossed in the reference case. The reason is that the Fourier amplitude of the interaction Hamiltonian decays exponentially as the resonance order K increases; only a handful of low-order resonances produce noticeable jumps, and the rest are exponentially suppressed. The paper therefore concludes that, in the regime where the DM orbit and the EMRI orbit do not radially overlap, the DM spike survives this resonant-scattering channel.","pith_inferences":["The paper stops where the DM and EMRI orbits begin to overlap; the authors expect the resonance effect to be enhanced there. If a matter-aware metric reconstruction confirms that expectation, the 'spike survives' conclusion would need to be narrowed to non-crossing orbits only.","The exponential decay of H with K mirrors the general fact that Fourier coefficients of a smooth function decay exponentially, so a similar smallness may hold for other environmental backreaction channels that scan through resonances.","A testable extension is to compute the same cumulative kicks with the EMRI's mass ratio pushed toward η ~ 10^-3 or with higher spin, where the decay onset in K shifts and the total effect could grow above the percent level.","The scatter of six orders of magnitude in H at fixed K is left open; if that scatter correlates with orbital phase or inclination, resonance networks could be modeled statistically rather than event-by-event."],"forward_implications":["For an EMRI with mass ratio 10^-4 around a 10^6-solar-mass black hole, a DM particle in the LISA band crosses on the order of 10^3 resonances, yet its semi-major axis changes by less than a percent.","The DM spike is not dissolved by tidal resonances as long as the DM particle's orbit stays outside the EMRI's radial range, so a spike can persist to be probed by EMRI gravitational-wave phase shifts.","Naive estimates that sum equal-sized resonance kicks overestimate the effect; exponential suppression of high-order Fourier harmonics makes the total kick small.","The small number of effective resonances means numerical studies can truncate the resonance sum at modest order K with verified convergence (K_max = 15 here)."],"supporting_citations":[{"why":"Establishes that a small perturbation around a Kerr geodesic causes net orbital drift only at tidal resonances, motivating the resonance search.","marker":"[20]"},{"why":"Provides the first explicit expression for orbital-element evolution across a resonance in Kerr, which this paper re-derives in a more formal way.","marker":"[21]"},{"why":"Supplies the metric reconstruction procedure from the Teukolsky variable used to derive the interaction Hamiltonian.","marker":"[32]"},{"why":"Gives the quadrupole energy-loss formula used to evolve the EMRI orbit by gravitational-wave backreaction.","marker":"[23]"},{"why":"Introduces the action-angle variables for Kerr geodesics in which the resonance dynamics is formulated.","marker":"[27]"},{"why":"Provides the relation between the change in action variables and resonance integers used to derive the resonant jump.","marker":"[29]"},{"why":"Supplies the Teukolsky mode functions used in the numerical evaluation of the interaction Hamiltonian.","marker":"[35]"},{"why":"Maps the conserved constants to the orbital elements (p,e,x) used in the numerical survey.","marker":"[36]"}],"fun_headline_variants":["Tidal resonances barely dent dark-matter spikes","DM spikes survive EMRI tidal resonances","Resonant scattering leaves DM spikes unscathed","Exponential suppression keeps DM spikes intact","EMRI kicks alter dark-matter orbits by <1%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation assumes the DM particle's orbit never crosses the EMRI secondary's radial range, and the numerics are stopped at the first overlap; in the overlapping regime, where the two bodies get closest, the resonance effects could be larger and are not computed.","fun_headline_variants_meta":{"raw":{"variants":["Tidal resonances barely dent dark-matter spikes","DM spikes survive EMRI tidal resonances","Resonant scattering leaves DM spikes unscathed","Exponential suppression keeps DM spikes intact","EMRI kicks alter dark-matter orbits by <1%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000704,"raw_usage":{"total_tokens":3149,"prompt_tokens":893,"completion_tokens":2256,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":509,"completion_tokens_details":{"reasoning_tokens":2184}},"tokens_in":509,"tokens_out":2256,"duration_ms":15505,"temperature":1.0,"reasoning_tokens":2184,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:17:06.662899+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same cumulative-resonance evolution with a metric-reconstruction scheme that is valid when the EMRI secondary and the DM particle orbits overlap, and check whether the total semi-major-axis change exceeds the ~1% level for the same fiducial parameters; if it does, the paper's conclusion that the DM spike survives this channel fails.","supporting_citations":[],"review_version":1}