---
title: Self-Force-Based Beyond-GR Waveform Model
url: https://www.emergentmind.com/topics/self-force-based-beyond-gr-waveform-model
type: topic
---

# Self-Force-Based Beyond-GR Waveform Model

Searching arXiv for recent papers on self-force-based beyond-GR waveform modeling and closely related self-force waveform frameworks.
A self-force-based beyond-GR waveform model is a waveform construction in which binary dynamics and radiation are organized in a small-mass-ratio expansion, while departures from general relativity are introduced through modified backgrounds, extra fields, altered regular fields and fluxes, and theory-dependent merger-ringdown structure. Within the literature represented here, the first explicit realization of this program for merger and ringdown is the small-mass-ratio beyond-GR framework of “Black hole mergers beyond general relativity: a self-force approach” [2510.11793]. At the same time, much of the practical architecture derives from earlier GR self-force models that define the reusable pipeline: specify a background spacetime, compute or fit a regularized self-force or equivalent invariant inputs, evolve the worldline or phase-space variables, and generate radiation from the resulting trajectory [1207.2839].

## 1. Definition and boundaries

The topic sits at the intersection of three distinct literatures: GR self-force waveform modeling, beyond-GR perturbation theory, and phenomenological merger deformations. These should not be conflated.

A useful classification is the following.

| Work | Status relative to the topic | Main role |
|---|---|---|
| [2510.11793] | Actual self-force-based beyond-GR waveform model | First-principles merger-ringdown beyond GR in the small-mass-ratio regime |
| [1207.2839] | Not beyond GR | GR reference model for orbit integration, conservative effects, and waveform generation |
| [2112.12265] | Not beyond GR | Second-order self-force inspiral waveform model through 1PA |
| [2506.02189] | Not beyond GR | Offline/online phase-space IMR architecture through plunge and ringdown |
| [2508.00087] | Not beyond GR | Second-order effective-source infrastructure |
| [1703.00865] | Not a waveform model | Frequency-domain PwP source-handling method for scalar self-force |
| [2309.14061] | Beyond-GR, but not self-force-based | Phenomenological merger deformation with flux-consistency closure |

A recurrent misconception is to label any self-force waveform as “beyond-GR.” The Schwarzschild self-force model “Self Force Orbit-Integrated gravitational waveforms for E(I)MRIs” is explicitly a GR construction and is best understood as a GR baseline rather than a beyond-GR model [1207.2839]. Conversely, the merger-ringdown framework of [2510.11793] is genuinely beyond GR, but only in the small-mass-ratio regime and only for a restricted class of effective-field-theory extensions. A separate misconception is to treat phenomenological merger deformations as self-force models; the physically consistent merger parameterization of [2309.14061] is beyond GR in intent but not derived from self-force theory.

## 2. General-relativistic reference architecture

The baseline GR scaffold is particularly clear in the Schwarzschild quasi-circular EMRI/IMRI model of [1207.2839]. The system is a point particle of mass \(\mu\) inspiraling into a Schwarzschild black hole of mass \(M\), with \(\mu \ll M\), and in the explicit numerical example \(\mu=10^{-2}M\). The inspiral runs from \(r_0=10M\) to near the ISCO at \(r=6M\). The worldline is evolved with the local first-order self-force according to
\[
u^{\beta}\nabla_{\beta}u^{\alpha}=\mu^{-1}f^{\alpha}_{\rm SF},
\]
with the orthogonality condition
\[
f^{\rm SF}_{\mu}u^{\mu}=0.
\]
The working force model uses fitted Lorenz-gauge circular-orbit data from Barack and Sago, with \(f^t\) encoding the dissipative sector and \(f^r\) encoding the conservative correction. Turning off the conservative part is implemented by setting \(b_i^A=0\), and in that limit the local self-force inspiral reproduces the energy-balance inspiral [1207.2839].

That paper also fixes the operational comparison standard that later beyond-GR work can inherit. It contrasts three evolutions: an energy-balance inspiral, a self-force inspiral with the conservative part turned off, and a full first-order self-force inspiral including the conservative piece. Waveform generation is not a kludge prescription; after orbit construction, the source trajectory is used to evolve the linearized Einstein equations via a sourced Teukolsky equation with hyperboloidal slicing, and the displayed waveform is the real part of \(\psi_4\). The main quantitative result is a cumulative dephasing
\[
\Delta\phi = 8.4 \pm 0.4\ {\rm rad},
\]
roughly \(\frac{4}{3}\) cycles out of \(54.3\) cycles, between the full-self-force and energy-balance waveforms over the inspiral from \(10M\) to near \(6M\). For the overlap study, the threshold window length is
\[
L_{\rm threshold}=816.6\,M,
\]
at which the overlap drops to \(0.96\), interpreted as a \(10\%\) event-rate loss. The same model shows that the gauge-invariant relation \(u^t(\Omega)\) is unchanged within numerical accuracy by the conservative force, while the time taken to traverse that curve is changed [1207.2839].

The later second-order inspiral model of [2112.12265] extends this architecture through first post-adiabatic order. There the inspiral is built from a two-timescale expansion of the Einstein equations through second order in the mass ratio, yielding waveform production in tens of milliseconds. The mode structure is written as
\[
h_{\ell m} = \Big\{\nu h_{\ell m}^{(1)} + \nu^2 \Big[h_{\ell m}^{(1)} +h_{\ell m}^{(2)}-\frac{2x}{3}\dfrac{dh_{\ell m}^{(1)}}{dx}\Big] \Big\}e^{-i m \phi_p},
\]
with the orbital evolution controlled by
\[
\frac{d\Omega}{dt} = \frac{\nu}{M^2}\left[ F_0(x)+\nu F_1(x)\right], \qquad \frac{d \phi_p}{dt} = \Omega.
\]
This provides the GR inspiral-side template for any beyond-GR generalization that can supply analogous invariant functions [2112.12265].

## 3. First-principles beyond-GR construction

The explicitly beyond-GR construction in [2510.11793] is formulated in an EFT language,
\[
{\cal L} = {\cal L}_{\rm EH}[{\sf g}] + {\cal L}_{\rm bGR}[\Psi,{\sf g}],
\]
with the common case of an extra scalar \(\varphi\). The regime is a small or intermediate mass ratio,
\[
\varepsilon := \frac{\mathring m_2}{\mathring m_1} \le 1,
\]
with formal validity at \(\varepsilon \ll 1\). The primary background is Schwarzschild, the motion is quasicircular leading to plunge, and the implementation is restricted to a nonspinning primary. The secondary is skeletonized as a point particle carrying a scalar monopole charge \(q\), with charge-to-mass ratio
\[
\lambda := \frac{q}{\mathring m_2}.
\]

The perturbative content is organized as
\[
h_{\alpha\beta} = \varepsilon h^{(1,0)}_{\alpha\beta} + \varepsilon^2\Bigl(h^{(2,0)}_{\alpha\beta} + \lambda^2 h^{(2,2)}_{\alpha\beta}\Bigr) + O(3),
\]
\[
\varphi = \varepsilon \lambda\, \varphi_{(1,1)} + O(2).
\]
The leading beyond-GR merger corrections therefore appear at order \(\varepsilon^2\lambda^2\). Varying the action yields a coupled metric-scalar system, and the small body obeys
\[
\frac{D^2x^\alpha_p}{d\tau^2} = P^{\alpha\beta}\left[ -\frac{1}{2} \bigl( 2\nabla_{\!\mu} h^{\rm R}_{\nu\beta} - \nabla_{\!\beta} h^{\rm R}_{\mu\nu} \bigr) u^\mu u^\nu + \frac{q}{m_2}\nabla_{\!\beta}\varphi^{\rm R} \right] + O(2),
\]
\[
\frac{dm_2}{d\tau} = - q u^\alpha \nabla_{\!\alpha}\varphi^{\rm R} + O(3),
\]
with
\[
P^{\alpha\beta}=g^{\alpha\beta}+u^\alpha u^\beta.
\]
At the order actually computed, the force decomposes as
\[
f^\alpha = \Bigl(f^\alpha_{(1,0)} + \lambda^2 f^\alpha_{(1,2)}\Bigr)+O(2),
\]
with
\[
f^\alpha_{(1,2)} = P^{\alpha\beta}_{(0)}\nabla_\beta\varphi^{\rm R}_{(1,1)},
\]
and the mass correction is
\[
m_2 = \mathring m_2 \Bigl[ 1 - \lambda^2 \varphi^{\rm R}_{(1,1)} + O(2) \Bigr].
\]

The beyond-GR content is modular in a precise sense. In the current implementation, the leading correction arises from the scalar monopole charge \(q\); theory-dependent curvature couplings such as \(\alpha^{(2)}\) enter only at higher order in the small-\(\varepsilon\) regime. This is why the paper can treat a broad class of EFT extensions while computing a specific leading correction. It also explains the caveat that in dynamical Chern-Simons gravity, where the leading scalar is dipolar rather than monopolar, the theory is indistinguishable from GR at this order [2510.11793].

## 4. Phase-space dynamics and waveform assembly

The beyond-GR merger-ringdown model of [2510.11793] is built by extending the GR post-geodesic plunge formalism into beyond GR. The worldline is parameterized as
\[
x^\mu_p(t,\varepsilon,\lambda) = (t, r_p(t,\varepsilon,\lambda), \theta_p=\pi/2, \phi_p(t,\varepsilon,\lambda)),
\]
and the reduced plunge dynamics is written as
\[
\frac{d\phi_p}{dt} = \Omega_{(0)} + \varepsilon\bigl[\Omega_{(1,0)}+\lambda^2\Omega_{(1,2)}\bigr] + O(2),
\]
\[
\frac{dr_p}{dt} = F_{(0)} + \varepsilon\bigl[F_{(1,0)}+\lambda^2F_{(1,2)}\bigr] + O(2).
\]
The geodesic-order terms are algebraic in \(r_p\), while the first post-geodesic corrections satisfy linear ODEs driven by \(f^\alpha_{(1,k)}\). The plunge solution is not stitched arbitrarily; it is fixed by asymptotic matching to a transition-to-plunge regime using
\[
\Delta r_p := \frac{r_p-6M}{\varepsilon^{2/5}}.
\]
This supplies the boundary data that determine \(\Omega_{(1,2)}\) and \(F_{(1,2)}\) [2510.11793].

Waveform assembly is likewise phase-space based. The asymptotic strain is extracted via
\[
h := \lim_{r\to\infty}\frac{r}{m_1} h_{\bar m\bar m},
\]
and the mode amplitudes are written as
\[
h_{lm} = e^{-im\phi_p} \Bigl\{ \varepsilon H^{(1)}_{lm}(r_p) + \varepsilon^2\bigl[ H^{(2,0)}_{lm}(r_p) + \lambda^2 H^{(2,2)}_{lm}(r_p) \bigr] + O(3) \Bigr\}.
\]
Operationally, waveform generation proceeds by solving for \(r_p(t)\) and \(\phi_p(t)\), evaluating precomputed \(H_{lm}^{(n,k)}(r_p)\), and multiplying by the phase factor \(e^{-im\phi_p}\) [2510.11793].

A closely related GR development makes the organizational principle explicit: the same offline/online phase-space paradigm used for inspiral can be continued through transition, plunge, merger, and ringdown [2506.02189]. In that framework, expensive field solutions are precomputed as functions on orbital phase space, while online waveform generation reduces to integrating a few ODEs and evaluating interpolants. The paper shows that this structure survives even when there is no separation of timescales during plunge. This suggests that the most portable component of self-force waveform modeling is not a particular GR equation, but the phase-space organization of the waveform problem itself [2506.02189].

## 5. Invariant structure and computational infrastructure

Invariant diagnostics play a dual role: they calibrate conservative dynamics and they separate physical content from gauge artifacts. In the Schwarzschild first-order model, the redshift invariant \(u^t(\Omega)\) is unchanged in shape, within numerical accuracy, by inclusion of the conservative force; the change lies in the rate at which the inspiral moves along the curve [1207.2839]. This suggests a two-layer validation strategy: first match invariant relations, then match timing and phase accumulation along them.

The invariant pseudo-Hamiltonian framework of [2507.08081] systematizes that idea at 1PA order. It reformulates the multiscale self-force problem on a six-dimensional phase space and writes each waveform mode as
\[
h_{lm} = \sum_{\bm{k}\in\mathbb{Z}^2} \left[ \epsilon\,\mathring h^{(1)}_{lm\bm{k}}(\mathring\pi_i) +\epsilon^2\,\mathring h^{(2)}_{lm\bm{k}}(\mathring\pi_i,\delta M_A) +O(\epsilon^3) \right] e^{-i k_i \mathring\varphi^i}.
\]
In its final invariant form, the dynamics becomes
\[
\frac{d\mathring\varphi^i}{dt} = \Omega^i_{(0)}(J_j) + \frac{\epsilon}{2}\frac{\partial\left\langle[\mathringH^{\rm sym}_{(1)}]\right\rangle}{\partial J_i} + {\cal O}(\epsilon^2),
\]
\[
\frac{dJ_i}{dt} = -\left\langle\left[\frac{\partial\mathringH^{\rm rad}_{(1)}}{\partial\mathring\varphi^i}\right]\right\rangle -\epsilon^2\left(\left\langle\left[\frac{\partial\mathringH_{(2)}}{\partial\mathring\varphi^i}\right]\right\rangle-\mathring K_i\right) +{\cal O}(\epsilon^3).
\]
A central byproduct is the identification
\[
\left\langle[H_{(1)}]\right\rangle=\mu\langle z_{(1)}\rangle,
\]
which leads to the statement that the on-shell conservative Hamiltonian equals the mechanical energy predicted by the first law of binary black-hole mechanics [2507.08081].

Practical second-order waveform construction requires effective-source technology. The technical infrastructure paper [2508.00087] constructs the multiscale second-order effective source on Schwarzschild and shows that it contains four essential pieces: quadratic coupling of first-order field modes, slow evolution of first-order fields, quadratic products of a first-order puncture field, and the second-order puncture field. Mode by mode, the second-order source is written as
\[
S^{2,\rm eff}_{i\ell m} = -16\pi \bar T^2_{i\ell m} +2\delta^2 R^0_{i\ell m}[\mathring h^1,\mathring h^1] +\frac{4a_{i\ell}}{rf}\left( E^1_{ij\ell m}h^1_{j\ell m} +E^0_{ij\ell m}h^{P2}_{j\ell m} \right),
\]
providing the regularized second-order input that underlies post-adiabatic self-force waveforms [2508.00087].

At the source-handling level, the frequency-domain particle-without-particle method of [1703.00865] is not a waveform model and not beyond GR, but it is structurally relevant. Its master scalar equation,
\[
\left(\square -V^{}_{\ell}(r) \right)\psi^{\ell m}= S^{\ell m}\delta (r-r^{}_{p}(t)),
\]
is turned into homogeneous left/right subdomain problems plus jump conditions by writing
\[
{\cal Q} = {\cal Q}_{-}\Theta^{-}_{p} + {\cal Q}_{+}\Theta^{+}_{p}.
\]
This is a computational template for localized worldline sources in any theory that reduces to wave-type equations with point-particle support [1703.00865].

## 6. Physical content, limitations, and adjacent directions

The beyond-GR merger-ringdown model of [2510.11793] already produces concrete waveform observables. For the \((2,2)\) mode and the example considered there, the scalar self-force shifts the invariant peak orbital frequency to
\[
m_1\Omega_{\rm peak}
\approx
\frac{1}{3\sqrt6}
+ 0.0029(8)\,\varepsilon^2\lambda^2,
\]
changes the peak amplitude to
\[
|h_{22}|_{\rm max}
\approx
\varepsilon
\left[
1.45 + 5.0(8)\times 10^{-3}\,\varepsilon\lambda^2
\right],
\]
and modifies the fundamental ringdown excitation amplitude according to
\[
A_{220}
\approx
\varepsilon\left\{
2.36+3.65i
+
\bigl[0.33(0)-0.80(7)i\bigr]\varepsilon\lambda^2
\right\}.
\]
Ringdown is beyond GR only through excitation in the current implementation; QNM frequency corrections are not yet computed because the quadratic scalar stress-energy term in the metric perturbation equation is omitted [2510.11793].

The principal limitations are explicit. The present beyond-GR realization is formally valid at small mass ratio, assumes a Schwarzschild primary, treats quasicircular plunge only, computes only a partial leading \(\varepsilon^2\lambda^2\) correction, and retains an effectively GR-like remnant/QNM spectrum. It does not cover Kerr spin, eccentricity, inclination, strong-coupling departures from the EFT hierarchy, or the dCS case at leading order because there the secondary has no scalar monopole [2510.11793].

A nearby but distinct direction is finite-size response of non-black-hole central objects. The compact-star self-force study [2406.02101] shows that, relative to a Schwarzschild black hole at the same orbital frequency, the additional self-force can be represented approximately by a universal frequency-dependent function multiplied by the dynamical tidal deformability. For the conservative radial sector,
\[
F^{\rm r,tide}=\lambda_2^{\rm dyn} c(M\Omega),
\]
and for the dissipative sector,
\[
P^{\rm tide}=(\lambda_2^{\rm dyn})^2a(M\Omega)+\lambda_2^{\rm dyn} b(M\Omega).
\]
In the illustrative EMRI estimate with \((M,\mu)=(10^6,10)M_\odot\), a four-year inspiral, and a supermassive star-like central object, the tide-induced orbital phase shift reaches \(\sim 4.5\times10^3\) rad, while the abstract summarizes the generic scale as \(\mathcal O(10^2)-\mathcal O(10^3)\) rad [2406.02101]. This is not a beyond-GR field-theory model, but it provides a response-theory template for exotic compact objects.

Phenomenological merger deformations remain an adjacent literature rather than part of the self-force line. The model of [2309.14061] introduces a compact beyond-GR merger parameterization on top of IMRPhenomD and enforces physical consistency by feeding the extra radiated energy and angular momentum back into the remnant spin and ringdown frequencies. Its relevance here is methodological: it supplies a flux-consistency closure principle for merger deformations, but it does not derive those deformations from self-force theory [2309.14061].

Taken together, these works define the present meaning of a self-force-based beyond-GR waveform model. In the strict sense, it is presently exemplified by a first-principles small-mass-ratio merger-ringdown framework with scalar-field corrections, asymptotic matching through plunge, and modular phase-space waveform assembly [2510.11793]. In the broader methodological sense, it is supported by a GR ecosystem that already supplies orbit-integrated first-order baselines, second-order post-adiabatic inspirals, effective-source regularization, invariant action-angle dynamics, and offline/online IMR organization [1207.2839], [2112.12265], [2508.00087], [2507.08081], [2506.02189]. The plausible implication is that future beyond-GR extensions will be won less by replacing the overall scaffold than by replacing its physics modules: background spacetime, regular fields, fluxes, waveform amplitudes, and remnant spectra.

Source: https://www.emergentmind.com/topics/self-force-based-beyond-gr-waveform-model