{"id":"8dc86982-4449-41b5-8907-c4553c34e2f4","arxiv_id":"2511.03657","paper_version":2,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A Schwarzschild black hole embedded in a Dehnen dark-matter halo admits homoclinic orbits, and raising halo density, scale radius, or particle energy drives near-horizon motion from regular to chaotic while respecting the surface-gravity chaos bound.","lead":"This paper studies a small star orbiting a black hole that sits inside a dark-matter halo and asks when the orbit becomes chaotic. It finds that the halo's density and size, plus the orbit's energy, can push the motion from regular to chaotic without violating a known bound set by the black hole's surface gravity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Metric–potential mismatch: Eq. (3.5) is not derived from Eq. (2.5); the DM term differs by a factor x/(2M), so the homoclinic tables and chaos bound are tied to a different spacetime.","rationale":"The reader identified the hand-imposed external potentials as the weakest assumption, which is a legitimate concern for the chaos section. However, I find an earlier and more fundamental issue: the effective potential used for the homoclinic orbit does not match the stated metric. This is not a matter of physical modeling choices—it is an internal consistency check. The paper's own Eq. (2.18) and Eq. (3.13) both correspond to f(r)=1-2M/r-32πρ_s r_s^2√(1+r_s/r), whereas the displayed Eq. (2.5) and Eq. (2.11) include an extra 1/r in the halo term. Eq. (3.5) is consistent with the no-1/r form, so the numerical tables and bounds are likely computed for the intended metric, and Eq. (2.5) may be a typographical error. Still, as written, the manuscript contains a mismatch that must be resolved before the quantitative claims can be trusted. The external-potential issue remains secondary: even if the metric is corrected, the chaos analysis is explicitly a toy model with imposed harmonic traps, so the physical claim about 'horizon-induced chaos' is conditional. The homoclinic construction itself is standard and could survive a corrected metric, but the specific values in Tables 1-2 and the L bounds would change. I therefore keep the reader's CONDITIONAL verdict rather than moving to ACCEPT or REJECT.","tokens_in":33034,"tokens_out":29333,"duration_ms":218272,"concrete_test":"Symbolically reduce Eq. (3.3) with f(r) from Eq. (2.5) using x=2M/r and compare with Eq. (3.5); the DM term acquires the extra factors x/(2M) and x^3/M. Then recompute Table 1 rows (ρ_s=0.01, L=3.75) and (ρ_s=0.04, L=4.25) using the correctly reduced V_eff. Also check whether Eq. (2.18) solves f(r)=0 for the metric actually used in Eq. (3.5). If r_un, r_max, or E_homo shift by more than ~1%, or if L_ISCO<L<L_MBO fails, the central homoclinic claim as stated is not supported by the paper's own metric.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative claims—homoclinic bounds L_ISCO<L<L_MBO, Tables 1–2, and the resulting phase-space separatrix—depend on the effective potential Eq. (3.5). But Eq. (3.5) does not follow from the stated metric Eq. (2.5). Starting from Eq. (3.3) with f(r) from Eq. (2.5), the transformation x=2M/r gives f_DM = -(32πρ_s r_s^2/r)√(1+r_s/r) = -(16πρ_s r_s^2 x/M)√(1+x r_s/(2M)). The DM part of f(m^2+L^2/r^2) then contributes, after multiplying by 4M^2/L^2, terms -64πρ_s r_s^2 M m^2 x/L^2 √(...) -16πρ_s r_s^2 x^3/M √(...), not the single printed term -32πρ_s r_s^2 (x^2 + 4M^2 m^2/L^2) √(...). The printed term is instead obtained from f_DM = -32πρ_s r_s^2 √(1+r_s/r) with no 1/r factor, which is also the form used in Eq. (3.13) and which is consistent with the horizon expression Eq. (2.18). Thus either Eq. (2.5) is a typo, or Eqs. (3.5), Tables 1–2, and the homoclinic bounds are computed for a different spacetime than the one introduced as the BH-DM halo metric. Since the homoclinic separatrix is the stated foundation for the chaos analysis, this inconsistency is the most load-bearing issue.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the motion of a massive test particle around a Schwarzschild-like black hole embedded in a Dehnen-(1,4,5/2) dark-matter halo. In the EMRI limit q=10^-5, the authors construct the radial effective potential, locate the ISCO and MBO, and identify a homoclinic orbit for L_ISCO < L < L_MBO that separates bound from plunging geodesics. The second half of the paper places the particle in Painlevé-Gullstrand coordinates and adds externally imposed harmonic potentials a(r) and b(phi). Numerical integration of the resulting equations yields Poincaré sections, orbital plots, and Lyapunov exponents, which the authors use to claim that increasing the halo density, scale radius, or particle energy drives the system from regular to chaotic motion while keeping both the total and radial Lyapunov exponents below the surface-gravity (MSS) bound.","tokens_in":33535,"tokens_out":13478,"duration_ms":109496,"significance":"If established, the homoclinic-orbit result would provide a concrete separatix in a realistic BH-DM halo spacetime and extend earlier near-horizon chaos studies to astrophysical environments. The paper has several strengths: it derives the metric from the Dehnen density profile, uses the conserved energy and angular momentum explicitly, presents numerical tables for homoclinic parameters, and is transparent about the fact that the harmonic potentials are externally imposed. However, the quantitative claims are currently undermined by an apparent inconsistency between the metric and the effective potential, and the chaos section describes a toy model whose connection to the BH-DM halo spacetime is not established. These issues must be resolved before the astrophysical conclusions can be accepted.","major_comments":[{"comment":"The effective potential (3.5) does not follow from the metric (2.5). Starting from Eq. (3.3) with f(r) = 1 - 2M/r - 32πρ_s r_s^2/r sqrt(1+r_s/r), the substitution x=2M/r gives the DM contribution (4M^2/L^2) f_DM (m^2 + L^2/r^2) = -16πρ_s r_s^2 x/M sqrt(1+x r_s/(2M)) (x^2 + 4M^2 m^2/L^2). The printed term in Eq. (3.5), -32πρ_s r_s^2 sqrt(1+x r_s/(2M)) (x^2 + 4M^2 m^2/L^2), is what one would obtain if f_DM had no 1/r factor. Since Eq. (3.5) is used for the ISCO/MBO calculations and for Tables 1-2, all homoclinic parameters and the claimed L_ISCO < L < L_MBO window are not yet shown to hold for the spacetime (2.5). Please correct Eq. (3.5) (or the metric) and recompute; the same inconsistency appears in Eq. (3.13), where the 1/r factor present in Eq. (3.12) is dropped.","section":"Sec. 3, Eq. (3.5)"},{"comment":"The chaos analysis uses hand-imposed harmonic potentials a(r)=1/2 K_r (r-r_c)^2 and b(phi)=1/2 K_phi r_H^2 phi^2, with K_r=100, K_phi=25, and r_c=3.65 chosen rather than derived. The paper explicitly calls these 'an externally imposed constraint rather than an inherent feature' of the spacetime. Therefore the Poincaré sections and Lyapunov exponents in Sec. 4.2 describe a particle in a confining trap, not a particle moving solely in the BH-DM halo geometry. The abstract's claim that increasing halo density/scale radius 'leads to chaos' in an EMRI embedded in a DM halo is not supported unless these potentials are shown to represent a realistic physical interaction. Please either derive the potentials from a physical model or clearly restrict the conclusions to this toy model and adjust the abstract and conclusions accordingly.","section":"Sec. 4.1, Eq. (4.4)"},{"comment":"The energy values used in the chaotic-dynamics section (E=90,95,108,115,118,122,130,132.5, and E=140-156 in Table 3) appear inconsistent with the stated EMRI scaling q=10^-5. In the geodesic section, bound energies are of order 10^-5 with M=1 and m=10^-5 (e.g., Table 1). In Eq. (4.6), with m=10^-5 and external potentials of order unity, E~100 would require momenta that are far larger than the near-horizon EMRI values used earlier. The paper does not state the units of E for the chaos section, and p_phi is computed from Eq. (4.6). This ambiguity affects all Lyapunov exponents and the comparison with kappa in Figs. 16-17. Please specify the normalization explicitly and ensure consistency with q=10^-5, or explain the rescaling used.","section":"Sec. 4.2, Eq. (4.6)"},{"comment":"The horizon radius formula (2.18) is not obviously the solution of f(r)=0 for f given in Eq. (2.5). As written, the structure of Eq. (2.18) differs from what one would obtain from the transcendental equation r - 2M - 32πρ_s r_s^2 sqrt(1+r_s/r)=0, which follows from Eq. (2.5) when the 1/r factor is included. Because r_H enters the surface gravity through Eq. (4.14) and is used in the figures and the chaos bound, the derivation of Eq. (2.18) should be shown explicitly, or the expression should be reconciled with Eq. (2.5).","section":"Sec. 2, Eq. (2.18)"}],"minor_comments":[{"comment":"The entry '10521493' appears to be missing a decimal point; it should likely read 10.521493.","section":"Table 2, row (0.35, 4.25)"},{"comment":"The text refers to 'the fourth and fifth terms' in Eq. (3.5), but Eq. (3.5) as written has four terms. Please revise the wording.","section":"Sec. 3.2, text after Eq. (3.8)"},{"comment":"The color coding of the Poincaré sections is not described in the text. Please add a note stating that different colors represent different randomly chosen initial conditions.","section":"Sec. 4.2.1, Figs. 7-10"},{"comment":"The quantities (λ_T^2 - κ^2) and (λ_r^2 - κ^2) are all negative and of order -0.05, with differences visible only at the level of a few 10^-6. Reporting λ/κ or the ratio (κ-λ)/κ would be more informative and would make the saturation of the bound easier to assess.","section":"Sec. 4.2.4, Figs. 16-17"},{"comment":"Reference [75] is an unpublished companion paper. If it is essential for the gravitational-wave claims in the introduction and conclusions, the connection should be described in more detail or the claims should be softened.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The metric-potential mismatch in Eq. (3.5) is the key risk. If the effective potential cannot be reconciled with Eq. (2.5), the homoclinic tables and the associated separatrix claims would need to be recomputed. The chaos section is also a toy model with externally imposed harmonic traps; the authors should either provide a physical origin for these potentials or substantially weaken the astrophysical conclusions. I do not see this as a rejection if the authors can correct the derivation and reframe the chaos section accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new item here is the homoclinic-orbit construction for a Schwarzschild-like black hole in a Dehnen-(1,4,5/2) halo, plus the claim that halo parameters shift the ISCO/MBO window and the Lyapunov exponents remain under the surface-gravity bound. That is a reasonable extension of existing near-horizon chaos studies, and the paper is honest: it states clearly that the external potentials in the chaos section are imposed, not derived, and it recovers the Schwarzschild limit correctly. The qualitative direction—more halo density or scale radius, or higher energy, tends to break KAM tori—is plausible and consistent with earlier work.\n\nThe soft spots are not minor. Most importantly, the effective potential in Eq. (3.5) does not follow from the metric in Eq. (2.5). If you start from f(r) = 1 - 2M/r - 32πρ_s r_s^2/r sqrt(1+r_s/r) and transform to x = 2M/r, the DM contribution to Veff contains terms proportional to x m^2 and x^3, not the printed (x^2 + 4M^2 m^2/L^2) sqrt(...). The printed form matches a metric missing the 1/r factor in the DM term. Since Tables 1–2 and the homoclinic separatrix are computed from this Veff, the paper's central quantitative results are tied to a different spacetime than the one advertised. This is load-bearing. Either Eq. (2.5) or Eq. (3.5) is a typo, but the authors must reconcile them before the homoclinic bounds can be trusted.\n\nThe chaos section is a separate, structural issue: the harmonic potentials a(r) and b(phi) are chosen, with K_r=100, K_phi=25, r_c=3.65, and the paper itself calls them an externally imposed constraint. That makes the Poincaré sections and Lyapunov exponents descriptive of the toy model, not of the DM halo or the horizon itself. The MSS-bound check is then a numerical observation on that model, not a theorem. No code or data is released, so reproducing the Lyapunov values is unnecessary work.\n\nWho should read this? People working on EMRI waveforms in environmental backgrounds will want to know whether homoclinic orbits survive in a consistent halo metric. That is worth a referee's time, but the paper needs major revision: fix the metric/potential mismatch, recompute or re-derive the tables, and either justify the external potentials physically or soften the horizon-induced-chaos claim.\n\nSend it to peer review, but with a clear request for these fixes.","headline":"A worthwhile but currently inconsistent paper: the homoclinic analysis rests on an effective potential that does not follow from the stated metric, and the chaos section is an openly toy-model study.","tokens_in":34076,"tokens_out":5912,"would_cite":false,"duration_ms":44886,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A Schwarzschild-like black hole embedded in a Dehnen dark-matter halo admits a homoclinic orbit that separates bound from plunging motion, and near-horizon particle motion becomes chaotic as halo density, scale radius, or energy increase—wh","keywords":["extreme-mass-ratio inspiral","dark matter halo","homoclinic orbit","chaos","Lyapunov exponent","surface gravity","Dehnen profile","black hole geodesics"],"falsifier":"Integrate the exact geodesic equations for the metric f(r)=1−2M/r−32π ρ_s r_s^2 sqrt(1+r_s/r)/r without the external harmonic potentials: the unstable equilibrium at r_un should yield a smooth homoclinic loop with zero Lyapunov exponent, and nearby trajectories should either plunge or escape. If, contrary to the paper's expectation, additional bounded chaotic trajectories appear without the traps, the horizon-and-halo origin of the chaos would be supported; if no chaos appears, the positive Lyapunov exponents are artifacts of the harmonic confinement. A more quantitative check: smoothly reduce","tokens_in":32883,"feed_emoji":"🌀","tokens_out":5118,"duration_ms":49329,"temperature":0.7,"pith_summary":"The paper argues that a massive compact object inspiraling around a supermassive black hole embedded in a Dehnen-(1,4,5/2) dark-matter halo can possess a homoclinic orbit—a phase-space separatrix between bound inspiral and plunging into the horizon—when its angular momentum lies between the ISCO and marginally-bound-orbit values. The halo's density and scale radius shift the ISCO and MBO and deform the effective potential. The paper further claims that adding near-horizon harmonic confinement (an external potential) and increasing energy, halo density, or scale radius drives the breakup of KAM tori, producing chaos, while both total and radial Lyapunov exponents stay below the surface gravity of the BH-DM spacetime, so the universal chaos bound is respected. If true, this links dark-matter halo structure to chaotic dynamics in EMRI systems and gives a theoretical handle on gravitational-wave signatures.","feed_headline":"Dark-matter halos can drive near-horizon orbits to chaos","feed_subtitle":"A homoclinic separatrix marks the plunge boundary, and Lyapunov exponents stay below the surface-gravity bound.","key_machinery":"The central object is the homoclinic orbit—the separatrix in the (r, p_r) phase space that emerges from the unstable circular orbit at r_un, defined by the effective potential of the BH-DM metric. The chaotic analysis is carried by the Hamiltonian in Painlevé-Gullstrand coordinates, in which the near-horizon particle is confined by two hand-imposed harmonic potentials, a(r)=½K_r(r−r_c)^2 and b(φ)=½K_φ r_H^2 φ^2; the paper then monitors Poincaré sections and Lyapunov exponents against the surface gravity κ of the combined spacetime.","core_discovery":"For a Schwarzschild-like black hole embedded in a Dehnen-(1,4,5/2) dark-matter halo, the paper establishes, in the extreme-mass-ratio limit q=m/M<<1, the existence of a homoclinic orbit for angular momenta L in (L_ISCO, L_MBO). This orbit, associated with an unstable circular orbit at r_un, forms the boundary between bound and plunging motion. The paper then shows numerically that in Painlevé-Gullstrand coordinates, with the particle confined near the horizon by external harmonic potentials, increasing the halo density ρ_s, scale radius r_s, or particle energy leads to the destruction of invariant tori and the onset of chaos, as measured by Poincaré sections and Lyapunov exponents. Throughou","pith_inferences":["The homoclinic-orbit existence result is geodesic and independent of the external harmonic potentials, so that part of the claim stands on the metric itself; the chaotic regime, however, is conditional on the traps—without them, as the paper's Appendix A shows, the particle would plunge or escape rather than remain confined.","If the hand-imposed harmonic potentials are read as proxies for generic near-horizon perturbations (magnetic fields, accretion disks, or scalar fields), the bounded-chaos result may generalize to those physical environments, but that extension is not established here.","A natural testable extension is to apply the Melnikov method to the homoclinic orbit with the harmonic perturbation, which would yield an analytic criterion for transversal homoclinic intersections and thus an independent prediction for the onset of chaos.","The paper's companion work on gravitational-wave signatures could look for a distinctive dephasing or burst when an inspiraling orbit crosses the homoclinic separatrix; that observable would distinguish chaos-induced modulations from ordinary adiabatic inspiral."],"forward_implications":["If the paper's claim is correct, EMRI systems in galactic nuclei with Dehnen-type halos can exhibit chaotic orbital motion in the near-horizon region, affecting the phase evolution of emitted gravitational waves.","The location of the chaos transition depends on halo density and scale radius, offering a potential way to infer dark-matter halo parameters from observed EMRI waveforms.","The existence of the homoclinic orbit fixes a sharp separatrix between inspiral and plunge, which could improve the accuracy of EMRI waveform templates that currently ignore such boundaries.","The surface-gravity bound on chaos is respected in the BH-DM halo environment, meaning that the universal MSS bound remains valid even when extended dark-matter structures are present.","The results provide a theoretical basis for a companion work connecting horizon-induced chaos to gravitational-wave observables for future space-based detectors."],"fun_headline_variants":["Homoclinic orbits mark chaos onset in dark-matter-halo EMRIs","Halo density and energy push inspirals to near-horizon chaos","Dark-matter halo drives EMRI orbits past homoclinic tipping point","Chaos in black-hole inspirals amplified by dark-matter halo","Near-horizon chaos from dark halo in extreme-mass-ratio systems"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The chaos claim depends entirely on the externally imposed harmonic potentials that confine the particle near the horizon; the paper itself states these are not an inherent feature of the BH-DM spacetime, and Appendix A shows an unconfined particle would plunge or escape.","fun_headline_variants_meta":{"raw":{"variants":["Homoclinic orbits mark chaos onset in dark-matter-halo EMRIs","Halo density and energy push inspirals to near-horizon chaos","Dark-matter halo drives EMRI orbits past homoclinic tipping point","Chaos in black-hole inspirals amplified by dark-matter halo","Near-horizon chaos from dark halo in extreme-mass-ratio systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00019,"raw_usage":{"total_tokens":1210,"prompt_tokens":812,"completion_tokens":398,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":299}},"tokens_in":556,"tokens_out":398,"duration_ms":4455,"temperature":1.0,"reasoning_tokens":299,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T23:54:12.787883+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the exact geodesic equations for the metric f(r)=1−2M/r−32π ρ_s r_s^2 sqrt(1+r_s/r)/r without the external harmonic potentials: the unstable equilibrium at r_un should yield a smooth homoclinic loop with zero Lyapunov exponent, and nearby trajectories should either plunge or escape. If, contrary to the paper's expectation, additional bounded chaotic trajectories appear without the traps, the horizon-and-halo origin of the chaos would be supported; if no chaos appears, the positive Lyapunov exponents are artifacts of the harmonic confinement. A more quantitative check: smoothly reduce","supporting_citations":[],"review_version":1}