---
title: 'TEOBResumS-Dal: Unified EOB Waveform for Generic Binaries'
url: https://www.emergentmind.com/topics/teobresums-dal
type: topic
---

# TEOBResumS-Dal: Unified EOB Waveform for Generic Binaries

TEOBResumS-Dal is a public branch of the TEOBResumS effective-one-body waveform framework that extends the EOB description of compact-binary dynamics beyond quasi-circular, comparable-mass black-hole binaries to regimes including large mass ratios, aligned spins, eccentricity, dynamical capture, scattering, and, in later formulations, generic spins, tides, and non-planar motion. Within the EOB picture, the binary is mapped to an effective particle moving in a deformed spacetime, with a conservative Hamiltonian, radiation reaction, and waveform generation. Across the literature, TEOBResumS-Dal and the spelling variant TEOBResumS-Dalí are used for implementations that target EMRI/IMRI-style inspirals, eccentric and unbound black-hole encounters, and eventually a unified model for binary black holes, binary neutron stars, and black-hole–neutron-star systems on arbitrary orbits through merger and including scattering [2310.13578] [2503.14580].

## 1. Definition, naming, and scope

TEOBResumS is an EOB waveform model in which the two-body problem is represented by an effective particle moving in an effective metric, and the waveform is built mode by mode from a factorized, resummed quasicircular baseline. In this framework, the key elements are a Hamiltonian and a radiation reaction force, which incorporate PN, NR, GSF, and test-mass information in a resummed way [2104.07533] [2503.14580].

Within that family, the “Dal” branch is identified as the public implementation used for EMRI/IMRI-style inspirals, with large-mass-ratio capability, aligned spins, eccentricity, and later improvements in the spin-orbit sector. In subsequent work, TEOBResumS-Dalí is presented as an eccentric-orbit iteration of the TEOB family, as a model for dynamical captures and scatterings, and finally as a unified EOB model for generic compact binaries with arbitrary orbits [2406.04108] [2410.05386] [2405.20398] [2503.14580].

The scope attributed to TEOBResumS-Dal broadens substantially across the cited works. One line of development focuses on inspiralling, quasi-circular large-mass-ratio binaries with aligned spins and eccentricity, benchmarked against second-order gravitational self-force information [2310.13578] [2406.04108]. Another line emphasizes noncircular and unbound black-hole dynamics, where the model is used for scatterings, dynamical captures, and waveform phenomenology close to the transition between bound and unbound motion [2405.20398]. A later synthesis describes TEOBResumS-Dalí as a single framework for black holes, neutron stars, and black-hole–neutron-star systems evolving along quasi-circular, eccentric, hyperbolic/scattering, and non-planar orbits through merger, ringdown, or post-merger emission as appropriate [2503.14580].

This suggests that “TEOBResumS-Dal” is best understood not as a single frozen approximant, but as a branch of the TEOBResumS infrastructure whose defining feature is the extension of EOB modeling to generic-orbit and large-mass-ratio problems.

## 2. Dynamical and waveform structure

The EOB structure underlying TEOBResumS-Dal combines a conservative Hamiltonian, dissipative radiation reaction, and multipolar waveform generation. For BBHs, the unified formulation uses a 5PN Hamiltonian with two parameters calibrated to NR quasi-circular simulations. Spin effects include spin-orbit couplings, spin-spin couplings, and spin-induced precession; for spin-spin interactions, the model uses a centrifugal radius \(r_c\), with PN information through NLO order. For BNSs, the Hamiltonian includes gravitoelectric terms, gravitomagnetic terms, and multipolar tidal couplings up to \(\ell=8\), with the reduced tidal deformability written as
\[
\tilde\Lambda=16/13\, (m_1+12m_2)m_1^4\Lambda_1/M + (1\leftrightarrow2).
\]
The default tidal choice in the unified model is the GSF-resummed potential [2503.14580].

The waveform sector is multipolar. In the TEOBResumS formulation assessed in the NQC and post-adiabatic study, the modes are written as
\[
h_{\ell m} = h_{\ell m}^{(N,\epsilon)} \hat{h}_{\ell m}\,\hat{h}^{\rm NQC}_{\ell m},
\]
with next-to-quasicircular corrections parameterized by
\[
\hat{h}^{\rm NQC}_{\ell m} = \left(1+a_1^{(\ell,m)}n^{(\ell,m)}_1+a_2^{(\ell,m)}n^{(\ell,m)}_2\right) \exp\!\left[i\left(b_1^{(\ell,m)}n^{(\ell,m)}_3+b_2^{(\ell,m)}n^{(\ell,m)}_4\right)\right].
\]
A central result of that work is the consistent use of NR-informed NQC information in both the waveform and the radiation reaction, replacing iterative determination of \((a_1^{(2,2)},a_2^{(2,2)})\) with direct fits across parameter space [2104.07533].

The same paper emphasizes the post-adiabatic approximation as a computational strategy for accelerating long inspirals. In that implementation, the PA evolution is typically carried out to 8th order on a radial grid with step \(dr=0.1\), from the starting radius \(r_0\) down to \(r=14\), after which the full ODEs are used. Comparing to 611 NR simulations for total mass \(10M_\odot\leq M \leq 200M_\odot\), the model has EOB/NR unfaithfulness well below \(0.01\), with \(78.5\%\) of the cases below \(0.001\) [2104.07533].

For generic-orbit applications, the unified TEOBResumS-Dal formulation uses multipolar waveform amplitudes and phases, multipolar radiation reaction, non-circular corrections for eccentric and hyperbolic motion, Padé resummations for strong-field robustness, and NR-informed waveform completion near merger. For BBHs, the model uses a \(3^{+3}\)PN description of all multipoles up to \(\ell=8\), with the \((2,2)\) mode treated at full 4PN order, and includes spin-dependent horizon absorption effects to LO [2503.14580].

## 3. Large-mass-ratio and GSF-informed developments

A defining application of TEOBResumS-Dal is the large-mass-ratio regime. The EMRI-oriented implementation incorporates aligned spins on both bodies and orbital eccentricity, while recasting conservative first-order gravitational self-force information into resummed EOB potentials and using a \(3^{+19}\)PN-accurate flux at infinity together with an improved horizon-flux prescription [2310.13578].

In that formulation, the conservative dynamics is encoded in EOB potentials such as \(A\), \(D\), and \(Q\), and the phasing is analyzed through
\[
Q_\omega \equiv \frac{\omega^2}{\dot\omega}
\]
with expansion
\[
Q_\omega = \frac{1}{\nu}Q_\omega^{(0)} + Q_\omega^{(1)} + \nu\,Q_\omega^{(2)}+\cdots.
\]
The infinity flux is computed as
\[
F^\infty \equiv \sum_{\ell=2}^{\ell_{\rm max}\sum_{m=-\ell}^{\ell} F_{\ell m}^{\rm Newt}|\hat h_{\ell m}|^2,
\]
and the practical improvement highlighted in the paper is the inclusion of the \(\ell=9\) and \(\ell=10\) modes. For a standard EMRI with \(q=5\times 10^4\) and \(m_2=10\,M_\odot\), the accumulated EOB/2GSF dephasing changes from \(3.1~{\rm rad}\) to \(-1.1~{\rm rad}\) over \((0.0224,0.12)\), and from \(3.0~{\rm rad}\) to \(-0.7~{\rm rad}\) over \([0.045,0.12]\). Over about 1.2382 years of evolution, corresponding to roughly \(1.5\times 10^5\) cycles, the dephasing is \(-0.7~{\rm rad}\), which the paper states is consistent with standard EMRI waveform accuracy requirements [2310.13578].

The 2024 spin-sector upgrade focuses on inspiralling, quasi-circular black-hole binaries with a non-spinning primary and a spinning secondary. The main change is a new gauge choice for the gyro-gravitomagnetic functions
\[
G_S,\qquad G_{S_*},
\]
which enter the spin-orbit Hamiltonian contribution
\[
\hat H_{\rm eff}^{\rm SO} = p_\varphi \left(G_S \hat S + G_{S_*}\hat S_*\right).
\]
The earlier TEOBResumS formulation uses the DJS gauge, in which \(G_S\) and \(G_{S_*}\) depend only on inverse radius and radial momentum. The updated TEOBResumS-Dal implementation adopts the \(\overline{\mathrm{DJS}}\) gauge, allowing dependence on
\[
p^2 = p_r^2 + p_\varphi^2 u^2,
\]
so that one can factor out the complete Kerr spinning-particle expression \(G_{S_*}^{K}\) as
\[
G_{S_*} = G_{S_*}^{K}\,\hat G_{S_*}.
\]
The paper concludes that the new \(\overline{\mathrm{DJS}}\) gauge improves agreement with GSF, especially in the linear-in-spin 1PA contribution: the extracted spin piece \(Q_\omega^{(1,s)}\) is closer to the GSF result, and the integrated spin-driven dephasing is roughly halved compared with the standard DJS choice over the comparison interval \(\omega \in (0.023, 0.09)\). These improved functions are implemented in the public TEOBResumS-Dal code [2406.04108].

The broader GSF motivation was already established in the earlier comparison between TEOBResumS, 2GSF/1PAT1, and NR. That work found that for \(q\le 10\), TEOBResumS is largely indistinguishable from NR while 1PAT1 has significant dephasing; for \(15\lesssim q\lesssim 64\), TEOBResumS and 1PAT1 agree well; and for \(q\gtrsim 64\), 1PAT1 becomes more accurate than TEOBResumS, with EOB–GSF dephasing reaching about \(-1.129~{\rm rad}\) at \(q=128\). The paper identifies a path toward GSF-informed EOB modeling through improved conservative potentials and dissipative fluxes, which is the direction later realized in the TEOBResumS-Dal large-mass-ratio program [2208.01049].

## 4. Eccentricity, noncircular motion, memory, and unbound dynamics

TEOBResumS-Dal is explicitly described as already equipped with eccentricity in the large-mass-ratio branch, and eccentricity remains a central ingredient in later generic-binary formulations [2406.04108] [2503.14580]. In the direct-current memory work, the model becomes the first EOB model to include DC memory contributions. The DC memory is defined as a non-oscillatory, hereditary component of the gravitational-wave signal, with permanent strain offset
\[
\Delta h_{\rm mem} = \lim_{t\to+\infty} h(t)-\lim_{t\to-\infty} h(t).
\]
The paper focuses on the purely non-oscillatory memory contribution in the \(m=0\) spherical-harmonic modes and transforms prior harmonic-coordinate, quasi-Keplerian results into EOB phase-space variables
\[
(r,\varphi,p_r,p_\varphi),
\]
using known 2PN coordinate transformations. The implementation has overall 2.5PN accuracy and preserves the small-eccentricity expansion up to \(e_t^6\) [2410.05386].

For \(m=0\) modes, the model uses
\[
h_{\ell 0} = \frac{8GM\nu}{c^2R}\sqrt{\frac{\pi}{5}\,H_{\ell 0},
\qquad
H_{\ell 0} = H_{\ell 0}^{\rm inst} + H_{\ell 0}^{\rm tail} + H_{\ell 0}^{\rm mem},
\]
with
\[
H_{\ell 0}^{\rm mem} = H_{\ell 0}^{\rm mem,osc} + H_{\ell 0}^{\rm DC},
\]
and the full \(m=0\) waveform assembled in a non-factorized form. A notable implementation detail is the use of \((u,p_{r_*},\dot p_{r_*})\) and the auxiliary combination
\[
Z_p=\sqrt{\dot p_{r_*}^2+p_{r_*}^2u^3},
\qquad
e_{\rm N}=\frac{Z_p}{u^2},
\]
which makes the eccentricity counting cleaner and helps remove spurious higher-order terms beyond the sixth eccentricity order. The paper reports that the DC term produces a clear, cumulative offset in the plus polarization and depends sensitively on the formation eccentricity \(e_i\), with \(e_i=1\) argued to be the most physically motivated default [2410.05386].

For highly noncircular and unbound BBH motion, TEOBResumS-Dal is validated against NR near the transition from scattering to dynamical capture. In that setting, the model uses a waveform factorization
\[
h_{\ell m}=h_{\ell m}^{(N,\epsilon)_{\rm c}\, \hat h_{\ell m}^{(N,\epsilon)_{\rm nc}\, \hat h_{\ell m}^{(\epsilon)_{\rm c}\, \hat h_{\ell m}^{(\epsilon)_{\rm nc},
\]
and for the dominant mode the generic Newtonian noncircular factor is
\[
\hat h_{22}^{(N,0)_{\rm nc} = 1 -\frac{\ddot r}{2r\Omega^2} -\frac{\dot r^2}{2r^2\Omega^2} +\frac{2i\dot r}{r\Omega} +\frac{i\dot\Omega}{2\Omega^2}.
\]
The paper emphasizes that factorizing this generic Newtonian term is crucial for good agreement at apastron and in dynamical capture. In the low-energy regime \(E_0\lesssim 1.02\,M\), mismatches weighted with the zero-detuned, high-power noise spectral density of Advanced LIGO are typically below or around the \(1\%\) level, with only a few spinning binaries slightly above the \(3\%\) threshold [2405.20398].

The unified generic-orbit paper systematizes these capabilities by treating scattering angle as a gauge-invariant diagnostic and by defining
\[
V(r;p_\varphi) \equiv H_{\rm EOB} (r,p_\varphi;p_{r_*}=0)
\]
for interpreting turning points and transitions between scattering and plunge. This places eccentric inspiral, hyperbolic scattering, and dynamical capture within a single EOB dynamical picture [2503.14580].

## 5. Matter effects, precession, and generic compact binaries

A major expansion of TEOBResumS-Dal occurs in BHNS and generic compact-binary modeling. The BHNS extension described in 2025 uses 52 new NR simulations to inform TEOBResumS-Dalí as a multipolar EOB model including precession and eccentricity. The model is designed to handle generic orbits, spin precession, and eccentricity, while incorporating BHNS-specific tidal-disruption physics in the merger and ringdown regime [2507.00113].

That work identifies a BHNS waveform morphology not reducible to BBH behavior with tides added. The merger phenomenology is classified into Type I, Type II, and Type III according to the ratio \(A^{\rm BHNS}_{22}/A^{\rm BBH}_{22}\), with thresholds \(<0.85\), \(>0.99\), and the intermediate range, respectively. The ringdown prescription depends on that classification, and the NQC quantities are extracted at
\[
t^{\rm NQC} \equiv t^{\rm peak} + 2,
\]
except for \((2,1)\), \((3,3)\), and \((4,4)\), where they are evaluated at \(t^{\rm peak}_{22}\). The paper finds that \(\Delta t_{\rm NQC}=4\) works better for BHNS than the BBH-inspired value 1, and reports that the new \((2,2)\) NQC prescriptions improve the merger amplitude by about an order of magnitude relative to the older model [2507.00113].

The same BHNS study highlights the special role of the \((2,0)\) and \((3,0)\) modes. It states that the \((3,0)\) amplitude can be of the same order as \((3,2)\) for \(q\gtrsim 2\), and that the \((2,0)\) can approach the \((2,1)\) contribution in the same region. The authors note that current TEOBResumS-Dalí does not yet include \(m=0\) modes there, identifying them as an important future direction [2507.00113]. This is consistent with the separate DC-memory implementation for \(m=0\) modes in the eccentric-orbit branch [2410.05386].

The unified generic model extends beyond BHNS. It is presented as the first EOB model for the general relativistic dynamics and gravitational radiation of generic compact binaries, explicitly encompassing BBH, BNS, and BHNS systems with tidal interactions, generic spins, multipolar radiation reaction and waveform, and NR-informed merger, ringdown, or post-merger completion. For BNSs, TEOBResumS-Dal includes NRPM or NRPMw post-merger descriptions; for BHNSs, it uses NR-informed remnant black-hole mass and spin models that smoothly deform the BBH case [2503.14580].

Precession in the unified formulation is treated via a hybrid PN-EOB scheme based on the separation of timescales between orbital motion and spin-precession motion. The parameter
\[
\chi_p = \max \left( |\bm{\chi}_{1,\perp}|, \frac{4+3q}{4q^2+3q} |\bm{\chi}_{2,\perp}| \right)
\]
is used in the generic-binary and BHNS descriptions [2503.14580] [2507.00113]. The BHNS validation against the 12-orbit precessing simulation BAM:0223 reports phase differences below \(\sim 0.5\) rad throughout the inspiral and mismatches including all the modes at the one percent level for low inclinations [2507.00113].

## 6. Validation, systematics, and known limitations

The broadest validation claim appears in the generic-binary formulation, which reports comparisons with 1395 high-accuracy NR simulations: 1183 BBHs, 27 BNSs, and 185 BHNSs. For BBHs, \(86.3\%\) of configurations have mismatch below \(0.01\), \(97.9\%\) below \(0.03\), and \(99.2\%\) below \(0.05\). For an eccentric BNS example with \(e_0^{\rm EOB}\sim 0.58\), the mismatch is \(0.9\%\). For a long BHNS simulation with \(q=2.7\), \(\chi_z=0.68\), and \(\Lambda=599\), the mismatch is about \(0.3\%\), while a precessing BHNS with \(\chi_p=0.61\) and \(\Lambda \simeq 494\) yields mismatch \(1.4\%\) [2503.14580].

The low-energy scattering study provides a complementary validation in a regime especially sensitive to noncircular structure. It finds good agreement for waveform phenomenologies, scattering angles, mismatches, and energetics when \(E_0\lesssim 1.02\,M\), but identifies deteriorating performance for higher energies, larger mass ratios, and some spinning direct captures. The paper attributes this to the PN nature of the Hamiltonian and radiation reaction and to the quasi-circular calibration of \(a_6\) and \(c_3\), while also noting that NQC corrections are turned off in that work [2405.20398].

For large-mass-ratio inspirals, the GSF-based benchmark papers supply a different notion of validation. The 2022 TEOBResumS-versus-1PAT1 comparison establishes that current EOB accuracy depends strongly on mass ratio, with a mutual-consistency region around \(q\sim 26\)–32 and growing dephasing at larger \(q\) unless the flux is improved [2208.01049]. The EMRI-oriented Dal implementation responds directly to this by using GSF-informed potentials and a \(3^{+19}\)PN flux with \(\ell\le 10\), obtaining sub-radian dephasing over a year-long standard EMRI [2310.13578]. The 2024 spin-orbit update then reduces the linear-in-spin 1PA discrepancy by changing the gauge choice for \(G_S\) and \(G_{S_*}\) [2406.04108].

A distinct limitation emerges from parameter-estimation systematics. In the 2026 SEOBNRv6EHM study, TEOBResumS-Dal is used as a comparison model for eccentric aligned-spin inference. The paper states that TEOBResumS-Dal can yield biased estimates of eccentricity, masses, and spins in the most challenging high-eccentricity configurations. For SXS:BBH:2527, the optimized mismatch quoted against NR is \(\mathcal{M}_{\rm TEOB} \sim 30\%\), compared with \(\mathcal{M}_{\rm v6} \sim 2\%\) for SEOBNRv6EHM, and the model is not used in the main 26-event LVK analysis because of its higher computational cost for long-duration signals. The same paper also notes that TEOBResumS-Dal does not currently support the backward-in-time waveform integration option used in some SEOBNR analyses [2605.28716].

These results indicate that TEOBResumS-Dal is highly developed and unusually broad in scope, but not uniformly optimal across all inference settings. A plausible implication is that its strongest role is as a physically unified EOB framework—especially for generic-orbit dynamics, GSF-informed large-mass-ratio inspirals, and matter-rich compact binaries—while continued calibration and algorithmic refinement remain necessary in the most demanding high-eccentricity and long-in-band parameter-estimation problems.

Source: https://www.emergentmind.com/topics/teobresums-dal