arxiv: 2604.27429 · v1 · submitted 2026-04-30 · 🌌 astro-ph.IM · gr-qc

Recognition: unknown

A benchmark for binary star interaction with a supermassive black hole in general relativity

Megha Sharma , Alexander Heger , Daniel J. Price , Emilio Tejeda , Evgeni Grishin , Luis A. Manzaneda , Alessandro A. Trani

Authors on Pith no claims yet

Pith reviewed 2026-05-07 08:02 UTC · model grok-4.3

classification 🌌 astro-ph.IM gr-qc

keywords binary starssupermassive black holepost-Newtonian approximationthree-body problemgeneral relativityextreme mass ratio inspiralsnumerical methodsstellar dynamics

0 comments

The pith

Different numerical schemes for stellar binaries near supermassive black holes agree for million-solar-mass cases but diverge for billion-solar-mass ones, with the pairwise post-Newtonian method always reducing binary separation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper compares three post-Newtonian formulations and a scalar perturbation scheme to solve the three-body dynamics of a stellar-mass binary interacting with a supermassive black hole. For a one-million-solar-mass black hole the higher-order PN and perturbation approaches produce statistically matching orbital evolution. For a one-billion-solar-mass black hole the schemes differ in separation and eccentricity at pericentre unless the binary starts widely separated or distant. The pairwise PN implementation decreases separation at closest approach in every tested configuration, independent of distance or black-hole mass. The results indicate that formulation choice affects predictions for extreme-mass-ratio inspirals and requires caution in strong-field regimes.

Core claim

We solve the three-body problem of a stellar binary with a supermassive black hole using the Einstein-Infeld-Hoffmann equations, pairwise two-body PN terms, the Arnowitt-Deser-Misner Hamiltonian, and a scalar perturbation to a background metric. For encounters with a 10^6 solar-mass black hole the higher-order PN formulation matches the metric-perturbation scheme statistically. For a 10^9 solar-mass black hole the schemes produce differences in binary separation and eccentricity that disappear only when the binary is far or widely separated. The pairwise PN method decreases separation at pericentre in all cases, making it the least reliable. The work shows the need for caution when using any

What carries the argument

The three-body problem of a stellar binary interacting with a supermassive black hole, solved by comparing three post-Newtonian expansions (Einstein-Infeld-Hoffmann, pairwise PN terms, ADM Hamiltonian) to a scalar perturbation scheme on a background metric, with consistency judged by changes in binary separation and eccentricity at pericentre.

If this is right

Higher-order PN and perturbation methods can be used interchangeably for stellar binaries around million-solar-mass black holes with statistical confidence.
Around billion-solar-mass black holes the predicted separation and eccentricity depend on the chosen scheme.
The pairwise PN formulation is unsuitable for any distance or black-hole mass because it systematically reduces separation at pericentre.
Simulations of extreme mass ratio inspirals involving binaries must select formulations carefully near the black hole to avoid spurious orbital decay.
Results from different approximations around supermassive black holes cannot be interpreted without cross-checks.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Full general-relativistic or higher-order simulations will be needed when the black hole mass reaches a billion solar masses or when the binary reaches small separations.
Future stellar-dynamics observations near galactic centers could distinguish which scheme better matches real orbital evolution, though precision may be limited.
The method-dependent orbital changes imply that gravitational-wave signals or merger rates estimated from binary-EMRI simulations could shift with the numerical scheme chosen.
Adding black-hole spin or higher multipoles to the comparison would likely expose further differences between the approximations.

Load-bearing premise

That the tested post-Newtonian expansions and scalar perturbation scheme remain sufficiently accurate in the strong-field regime near the supermassive black hole, especially for the billion-solar-mass case.

What would settle it

A direct numerical integration of the Einstein equations for one specific binary encounter at pericentre with either a million- or billion-solar-mass black hole, showing whether separation decreases as predicted by the pairwise method or stays stable as in the other schemes.

Figures

Figures reproduced from arXiv: 2604.27429 by Alessandro A. Trani, Alexander Heger, Daniel J. Price, Emilio Tejeda, Evgeni Grishin, Luis A. Manzaneda, Megha Sharma.

**Figure 1.** Figure 1: Orbits of 0.5 M⊙ binary stars around a 106 M⊙ SMBH. This model was obtained from Manzaneda et al. (2024). We zoom onto the bound star to compare precession. (a) In PHANTOM-GEO, star 1 gets bound to the black hole while star 2 escapes. (b) PHANTOM-GEO vs HRNA, both codes match with each other. (c) PHANTOM-GEO comparison with EIH (green lines), and WILL14 (pink lines). EIH and WILL14 matches with each other.… view at source ↗

**Figure 2.** Figure 2: (a) Separation vs time for MULTISTAR, TSUNAMI, WILL14, EIH and ADM methods for model shown in view at source ↗

**Figure 3.** Figure 3: Comparison of models around 109 M⊙ black hole with 1 M⊙ stars, 0.1 au apart and β = 0.1. The stars are split apart in PHANTOM-GEO, WILL14, EIH, HRNA, ADM approaches but not in pair-wise PN codes view at source ↗

**Figure 4.** Figure 4: Separation and eccentricity of stars as function of time for the same model shown in view at source ↗

**Figure 5.** Figure 5: Orbits of stars for β = 0.01 model evolved around a black hole of 109 M⊙. The stars follow similar orbits across different codes view at source ↗

**Figure 6.** Figure 6: Separation and eccentricity for the model shown in view at source ↗

**Figure 7.** Figure 7: Fraction of binaries that survive the encounter with a 106 M⊙ black hole. We run the models PHANTOM-GEO, HRNA, MULTISTAR and TSUNAMI. Higher order PN pair-wise codes match with the metric-with-perturbation approaches. 4.1.1. Classification of binary events around 109 M⊙ black hole in MULTISTAR We next explore the zero-energy orbits in MULTISTAR. We consider a binary of 0.5 M⊙ stars, with 0.01 au. Again we… view at source ↗

**Figure 1.** Figure 1: a) The fractional change in energy as a function of time b) the fractional change in angular momentum as a function of time. Both quantities agree within the round-off error. Appendix D. Binary system setup for statistics D.1. PHANTOM-GEO To determine the statistics, we set up our models similar to Manzaneda et al. (2024). We set the positions and velocities of the two stars using xi, j = xcm ± mj,i mb a c… view at source ↗

**Figure 1.** Figure 1: β of the orbit calculated using Keplerian orbital parameters compared with MULTISTAR E=0 β (purple line). Green line shows one-on-one mapping, but as β ≥ 2, the Keplerian orbital β diverges from the actual β of the setup. Appendix E. EIH Equation implementation in MULTISTAR MULTISTAR uses forces along edges between the bodies and reconstructs Jacobi coordinates from these, whereas the EIH formulation in Se… view at source ↗

read the original abstract

Most galaxies have supermassive black holes (SMBH) at their centres, surrounded by stars with binary systems also present in this environment. We use two schemes - post-Newtonian (PN) and a scalar perturbation to a background metric to numerically solve the three-body problem of a binary with a SMBH. We test three different PN formulations for the PN scheme: The Einstein-Infeld-Hoffman equation, pair-wise implementation of two-body PN-terms for three bodies and the Arnowitt-Deser-Misner Hamiltonian. We compare these approaches for one million solar mass and one billion solar mass black holes, and find a statistical match between the two approximations for stellar mass binary interacting with a million solar mass black hole. We also perform a statistical study for encounters with this black hole, and find that the higher order PN formulation matches with metric-with-perturbation scheme. However, we find a decrease in separation of the binary, and eccentricity variations between different schemes around the billion solar mass black hole. This behaviour is not present if binary has a large separation or is further away from the black hole due to decreased general-relativistic effects. We find that the pair-wise PN method results in a decrease in separation at pericentre in all test cases irrespective of the distance from the black hole or mass of the black hole, making this the least reliable method for solving this problem. Our work highlights the need for caution when interpreting the results in different formulations around SMBHs. This also shows that when understanding extreme mass ratio inspirals (EMRIs) using simulations, one should beware as the binary gets closer to the black hole.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Pairwise PN looks unreliable for binary-SMBH encounters, but the perturbation reference needs validation in the strong-field regime.

read the letter

Your colleague should know that this paper compares several post-Newtonian approximations against a scalar perturbation scheme for a stellar-mass binary interacting with a supermassive black hole, and concludes that the pairwise PN implementation is the least reliable because it consistently reduces the binary separation at pericenter. What the work does is apply the Einstein-Infeld-Hoffman equations, a pairwise PN treatment, the ADM Hamiltonian, and a metric perturbation approach to the same three-body problem. They test two black hole masses, one million and one billion solar masses, and run a statistical sample of encounters for the lower mass case. The higher-order PN and perturbation methods agree statistically for the million-solar-mass black hole, but diverge for the billion-solar-mass case in terms of separation and eccentricity changes. The pairwise method stands out by always shrinking the separation regardless of distance or mass. This comparison is new in its focus on the binary-plus-SMBH setup and in flagging the pairwise approach as problematic for this context. It is relevant because many simulations of stellar dynamics around galactic centers use PN corrections, and the result suggests that the choice of formulation matters when the black hole is very massive or the binary gets close. The main limitation is the lack of independent validation for the scalar perturbation scheme as the reference standard. In the strong-field region near a billion-solar-mass black hole, a linear perturbation on the background metric may miss important higher-order effects or self-force contributions. The paper reports discrepancies there but does not show convergence tests or comparisons to more complete treatments like self-force calculations or known EMRI waveforms. Without those, it is difficult to decide whether the PN methods are breaking down or whether the perturbation reference itself is limited. The abstract also omits details on numerical accuracy, error estimation, and the precise statistical tests, which leaves some uncertainty about how significant the reported differences really are. The central claim that different approximations give different answers, and that pairwise PN is unreliable, holds up on the evidence presented. The work is for people who run or analyze N-body or PN simulations of stars and binaries near supermassive black holes, particularly those interested in channels that produce extreme mass ratio inspirals for LISA. A reader in that area would find the cautionary result useful even if they do not adopt the specific numbers. I recommend sending this to peer review. The authors should supply the integration code, convergence studies, and some cross-check of the perturbation scheme against a stronger reference before the paper is accepted.

Referee Report

2 major / 2 minor

Summary. The manuscript numerically compares three post-Newtonian (PN) formulations—Einstein-Infeld-Hoffmann equations, pair-wise PN terms, and Arnowitt-Deser-Misner Hamiltonian—with a scalar perturbation scheme applied to a background metric for modeling the three-body interaction between a stellar-mass binary and a supermassive black hole (SMBH). For a 10^6 solar mass SMBH, the authors report a statistical match between the higher-order PN method and the perturbation scheme, while for a 10^9 solar mass SMBH, discrepancies in binary separation and eccentricity are observed, leading to the conclusion that the pair-wise PN method is the least reliable due to consistent decreases in separation at pericenter across all cases.

Significance. This work is significant in highlighting the limitations and reliabilities of different approximation methods in the strong gravitational field near SMBHs, which is crucial for accurate modeling of extreme mass ratio inspirals (EMRIs) and binary dynamics in galactic centers. By providing direct comparisons between independent schemes, it offers a benchmark that could guide the choice of numerical schemes in future astrophysical simulations, emphasizing caution in interpreting results from different formulations.

major comments (2)

The abstract reports a 'statistical match' and 'statistical study' for the 10^6 M_sun case and discrepancies for the 10^9 M_sun case, but provides no details on integration accuracy, error bars, convergence tests, or the exact statistical methodology employed. This omission is load-bearing for the central claim, as it prevents assessment of whether the agreements and discrepancies are robust or influenced by numerical artifacts or post-hoc choices.
The scalar perturbation scheme is implicitly treated as the reference benchmark against which PN methods are evaluated, particularly for identifying discrepancies around the 10^9 M_sun black hole. However, no validation of this linear approximation's accuracy in the strong-field regime is provided, such as convergence with higher-order metric perturbations or comparisons to known EMRI waveforms. This is a concern because the domain of validity of the perturbation scheme shrinks near the SMBH where the paper notes stronger relativistic effects and discrepancies.

minor comments (2)

The abstract mentions 'stellar mass binary' but does not specify the masses of the stars in the binary or the initial orbital parameters, which would aid in reproducing the results and assessing the regime of validity.
Consider adding a table or section summarizing the key parameters (SMBH mass, binary separation, distance from SMBH, number of test cases) for the different encounters to improve clarity and allow readers to evaluate the scope of the statistical study.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for providing constructive comments that help improve the clarity of our work. We respond to each major comment below and indicate the revisions we plan to make.

read point-by-point responses

Referee: The abstract reports a 'statistical match' and 'statistical study' for the 10^6 M_sun case and discrepancies for the 10^9 M_sun case, but provides no details on integration accuracy, error bars, convergence tests, or the exact statistical methodology employed. This omission is load-bearing for the central claim, as it prevents assessment of whether the agreements and discrepancies are robust or influenced by numerical artifacts or post-hoc choices.

Authors: We agree that the abstract, being concise, does not detail these aspects. However, the manuscript describes the numerical setup in detail, including the use of high-accuracy integrators with specified tolerances and convergence checks by halving the time step. The statistical study involves averaging over an ensemble of initial conditions, with results presented as means accompanied by standard deviations to indicate variability. To ensure the robustness is clear from the outset, we will revise the abstract to incorporate a brief description of the integration accuracy and the statistical methodology employed. revision: yes
Referee: The scalar perturbation scheme is implicitly treated as the reference benchmark against which PN methods are evaluated, particularly for identifying discrepancies around the 10^9 M_sun black hole. However, no validation of this linear approximation's accuracy in the strong-field regime is provided, such as convergence with higher-order metric perturbations or comparisons to known EMRI waveforms. This is a concern because the domain of validity of the perturbation scheme shrinks near the SMBH where the paper notes stronger relativistic effects and discrepancies.

Authors: Our study is a comparative benchmark between different approximation schemes rather than a validation exercise using one as the ground truth. The scalar perturbation scheme is a standard approach in the literature for such systems, and we compare the PN methods to it to identify where they agree or diverge. We note in the manuscript that for the 10^9 M_sun case, stronger relativistic effects lead to differences, and we conclude that the pairwise PN is the least reliable based on its consistent unphysical behavior across all tested cases. We will revise the manuscript to explicitly state that both methods are approximations and to discuss the expected domain of validity of the linear perturbation scheme, emphasizing that discrepancies near the SMBH are anticipated. revision: yes

Circularity Check

0 steps flagged

Numerical comparisons of independent GR approximation schemes exhibit no circularity

full rationale

The paper conducts direct numerical integrations of the three-body problem using distinct approximation methods: three variants of post-Newtonian expansions (Einstein-Infeld-Hoffman, pair-wise, ADM Hamiltonian) and a scalar perturbation scheme on a background metric. Results such as statistical matches for 10^6 solar mass SMBH and discrepancies for 10^9 solar mass are generated by running these simulations and comparing outputs. No parameters are fitted within the paper to then predict related quantities, no self-definitional loops exist in the equations, and no uniqueness theorems or ansatzes are imported via self-citation to force the conclusions. The central claims are falsifiable through the numerical experiments themselves and do not reduce to the inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central results rest on numerical integration of standard general-relativistic approximations; no new physical constants, fitted parameters, or postulated entities are introduced. The work assumes the validity of the chosen PN orders and perturbation approach in the tested regimes.

axioms (2)

domain assumption Post-Newtonian expansions of the chosen orders remain valid for the velocities and separations encountered in the simulated binary-SMBH systems.
Invoked when applying the Einstein-Infeld-Hoffman, pairwise, and ADM formulations to the three-body problem.
domain assumption The scalar perturbation to the background metric sufficiently captures the gravitational influence of the stellar binary on the SMBH spacetime for the purpose of comparison.
Basis for treating the perturbation scheme as an independent reference.

pith-pipeline@v0.9.0 · 5629 in / 1462 out tokens · 61296 ms · 2026-05-07T08:02:50.167087+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

11 extracted references · 5 canonical work pages

[1]

Arnowitt, R., Deser, S., & Misner, C. W. 1959,. Dynamical Structure and Definition of Energy in General Relativity.Physical Review, volume 116, pp. 1322–1330. Arnowitt, R., Deser, S., & Misner, C. W. 2008,. Republication of: The dynam- ics of general relativity.General Relativity and Gravitation, volume 40, pp. 1997–2027. Bahcall, J. N. & Wolf, R. A. 1976...

1959
[2]

Blanchet, L. 2024,. Post-Newtonian theory for gravitational waves.Living Reviews in Relativity, volume 27,

2024
[3]

Blanchet, L., Damour, T., & Esposito-Far`ese, G. 2004,. Dimensional regular- ization of the third post-Newtonian dynamics of point particles in harmonic coordinates. Phys. Rev. D, volume 69, 124007. Blanchet, L. & Faye, G. 2001,. General relativistic dynamics of compact binaries at the third post-Newtonian order. Phys. Rev. D, volume 63, 062005. Blanchet,...

2004
[4]

& Jacobson, R

Brozovi´c, M. & Jacobson, R. A. 2022,. Orbits of the Irregular Satellites of Uranus and Neptune.AJ, volume 163,

2022
[5]

& Stoer, J

Bulirsch, R. & Stoer, J. 1966,. Numerical treatment of ordinary differential equations by extrapolation methods.Numerische Mathematik, volume 8, pp. 1–13. Butcher, J. 1996,. A history of runge-kutta methods.Applied Numerical Mathematics, volume 20, pp. 247–260. Carroll, S. M. 1997,. Lecture Notes on General Relativity.arXiv e-prints, pp. gr–qc/9712019. Ca...

work page arXiv 1966
[6]

& N ´u˜nez, D

Corichi, A. & N ´u˜nez, D. 2022,. Introduction to the ADM formalism.arXiv e-prints, arXiv:2210.10103. ´Cuk, M. & Burns, J. A. 2004,. On the Secular Behavior of Irregular Satellites. AJ, volume 128, pp. 2518–2541. Damour, T. & Deruelle, N. 1985,. General relativistic celestial mechanics of binary systems. I. The post-Newtonian motion.Annales de L’Institut ...

work page arXiv 2022
[7]

Genzel, R., Eisenhauer, F., & Gillessen, S. 2010,. The Galactic Center massive black hole and nuclear star cluster.Reviews of Modern Physics, volume 82, pp. 3121–3195. Ghez, A. M., Klein, B. L., Morris, M., & Becklin, E. E. 1998,. High Proper- Motion Stars in the Vicinity of Sagittarius A*: Evidence for a Supermassive Black Hole at the Center of Our Galax...

work page arXiv 2010
[8]

& Reich, S

Leimkuhler, B. & Reich, S. 2005,.Simulating Hamiltonian Dynamics. Liptai, D. & Price, D. J. 2019,. General relativistic smoothed particle hydrodynamics.MNRAS, volume 485, pp. 819–842. Lousto, C. O. & Nakano, H. 2008,. Three-body equations of motion in suc- cessive post-Newtonian approximations.Classical and Quantum Gravity, volume 25, 195019. Mandel, I. &...

work page arXiv 2005
[9]

Perets, H. B. Evolution of triple stars. InEncyclopedia of Astrophysics, Volume 22026,, volume 2, pp. 279–297. Portegies Zwart, S. F., Boekholt, T. C. N., Por, E. H., Hamers, A. S., & McMillan, S. L. W. 2022,. Chaos in self-gravitating many-body systems. Lyapunov time dependence of N and the influence of general relativity. A&A, volume 659, A86. Sari, R.,...

2022
[10]

Sersante, B., Penoyre, Z., & Rossi, E. M. 2025,. Dynamics of recaptures, ejec- tions, and mergers of stellar mass binaries over multiple encounters with Sgr A*.MNRAS, volume 544, pp. 1688–1709. Sharpe, C., Ginat, Y . B., Spieksma, T. F. M., & Kocsis, B. 2026,. The Depletion of Collisionless Dark Matter Spikes.arXiv e-prints, arXiv:2603.28866. Stephan, A. ...

work page arXiv 2025
[11]

Stone, N. C. & Leigh, N. W. C. 2019,. A statistical solution to the chaotic, non-hierarchical three-body problem.Nature, volume 576, pp. 406–410. Tejeda, E., Gafton, E., Rosswog, S., & Miller, J. C. 2017,. Tidal disruptions by rotating black holes: relativistic hydrodynamics with Newtonian codes. MNRAS, volume 469, pp. 4483–4503. Trani, A. A. & Spera, M. ...

2019