pith. sign in

arxiv: 2606.04998 · v1 · pith:FPZYB3V3new · submitted 2026-06-03 · 🌀 gr-qc

Self-force calculations with numerical relativity methods

Nils L. Vu , Nami Nishimura , Thomas Osburn , Jonathan E. Thompson , Lawrence E. Kidder , Samuel D. Upton , Barry Wardell This is my paper

Pith reviewed 2026-06-28 05:15 UTC · model grok-4.3

classification 🌀 gr-qc
keywords self-forceKerr spacetimediscontinuous Galerkinextreme mass-ratio inspiralsscalar chargem-mode decompositionnumerical relativitypoint particle
0
0 comments X

The pith

A numerical method using m-mode decomposition and discontinuous Galerkin discretization computes the scalar self-force in Kerr spacetime with exponential convergence up to spins of 0.998 on orbits as close as the ISCO.

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

The paper introduces a computational scheme that applies m-mode separation of variables together with null slicing in horizon-penetrating coordinates to turn the scalar self-force problem into a set of elliptic PDEs. These equations are discretized with high-order discontinuous Galerkin elements, adaptive mesh refinement, and a parallel multigrid-preconditioned Krylov solver. The resulting scheme produces exponential convergence for the self-force on a scalar point charge even when the charge follows circular equatorial orbits at the ISCO in Kerr spacetimes with spins as high as a=0.998. The approach is implemented inside the open-source SpECTRE code and solves twenty m-modes in parallel within seconds while retaining the structure needed for later extension to gravitational perturbations and generic orbits.

Core claim

Performing an m-mode separation of variables, adding null vtu slicing in horizon-penetrating coordinates, and solving the resulting elliptic PDEs with high-order discontinuous Galerkin discretization yields exponential convergence for the self-force on a scalar point charge in Kerr spacetime, up to spins a=0.998 on circular equatorial orbits as close as the ISCO, despite the non-smooth puncture on the grid.

What carries the argument

m-mode decomposition with null vtu slicing that converts the wave operator into elliptic PDEs, discretized by high-order discontinuous Galerkin finite elements with adaptive refinement and multigrid-Schwarz preconditioning.

If this is right

  • The method supplies the scalar self-force data required as a first step toward second-order gravitational self-force calculations in Kerr.
  • It handles both prograde and retrograde orbits at the ISCO for near-extremal spins.
  • Twenty m-modes can be obtained in parallel on modest resources, opening the route to repeated evaluations during inspiral modeling.
  • The same infrastructure inside SpECTRE can be reused once the scheme is extended to gravitational perturbations.
  • The framework already accommodates more generic orbits once the circular-equatorial restriction is lifted.

Where Pith is reading between the lines

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

  • The same elliptic reduction might be applied to compute self-forces along non-circular or inclined trajectories without changing the core discretization.
  • Embedding the solver inside existing numerical-relativity pipelines could shorten the timeline for producing second-order EMRI waveform templates.
  • The regularity properties demonstrated for the scalar puncture suggest the approach could also treat other singular sources, such as those appearing in black-hole perturbation theory.
  • Because the code is public, independent groups can test the same setup on different background spacetimes or with different matter models.

Load-bearing premise

The m-mode separation together with null slicing produces elliptic PDEs whose solutions remain regular enough for the high-order discontinuous Galerkin scheme to keep exponential convergence around the point-charge puncture.

What would settle it

If the measured convergence rate of the self-force error drops from exponential to algebraic when spin reaches 0.998 or the orbit reaches the ISCO, the claimed robustness of the discretization would be falsified.

Figures

Figures reproduced from arXiv: 2606.04998 by Barry Wardell, Jonathan E. Thompson, Lawrence E. Kidder, Nami Nishimura, Nils L. Vu, Samuel D. Upton, Thomas Osburn.

Figure 1
Figure 1. Figure 1: FIG. 1. Computational domain for [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Preconditioning of the iterative linear solver. Each GMRES iteration is preconditioned with one multigrid V-cycle, [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Convergence of the DG method for modes up to [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: gives an overview of the m-mode solutions we obtain with our method. All m-modes up to m = 20 are shown on a log scale, demonstrating the decay in power over mode number. Note that regular zero crossings in the solution of the regularized field within the worldtube display as sharp drops on the chosen log scale, particularly for the m = 0 mode. Also visible are features in the solution within the worldtube… view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Convergence of the preconditioned linear solver, ex [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Error budget in the extracted self-force. Upper panels show the value of the self-force and lower panels the relative [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
read the original abstract

To model gravitational waveforms from extreme mass-ratio inspirals (EMRIs) for the upcoming LISA space mission, gravitational self-force calculations are needed to second order in perturbation theory. However, to date these calculations have only been attempted for the simplest case of circular orbits in Schwarzschild spacetime. In this work, we present a new computational method aimed at performing generic second-order self-force calculations in Kerr spacetime using methods from the adjacent field of numerical relativity. We perform an $m$-mode separation of variables, add null ("$vtu$") slicing in horizon-penetrating coordinates, and solve the resulting elliptic PDEs using high-order discontinuous Galerkin discretization, adaptive mesh-refinement, and an iterative Krylov-type linear solver with parallelizable multigrid-Schwarz preconditioning. We find that our method achieves exponential convergence for the self-force on a scalar point charge in Kerr spacetime up to spins of $a=0.998$ (Thorne limit) on circular equatorial orbits as close as the ISCO (prograde and retrograde), despite the non-smooth puncture on the grid. We solve for 20 $m$-modes in parallel in a few seconds and retain the flexibility to extend the method to gravitational self-force and more generic orbits in the future. The code to perform these calculations is publicly available in the open-source numerical relativity code SpECTRE.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a new numerical method for computing the scalar self-force on a point charge in Kerr spacetime. It combines m-mode separation of variables, null (vtu) slicing in horizon-penetrating coordinates, high-order discontinuous Galerkin discretization with adaptive mesh refinement, and a parallelizable multigrid-Schwarz preconditioned Krylov solver. The central claim is that this approach achieves exponential convergence for the self-force on circular equatorial orbits up to a=0.998 at the ISCO (prograde and retrograde), despite the non-smooth puncture, with the code released publicly in SpECTRE; the framework is positioned as extensible to second-order gravitational self-force and generic orbits.

Significance. If the exponential convergence holds with quantitative substantiation, the work would represent a meaningful advance toward second-order self-force calculations in Kerr, directly relevant to LISA EMRI waveform modeling. The public availability of the SpECTRE implementation and the parallel m-mode treatment are concrete strengths that could facilitate community adoption and extension beyond the scalar circular case.

major comments (2)
  1. [Abstract] Abstract: the claim that exponential convergence is achieved for a=0.998 at the ISCO is load-bearing for the paper's central result, yet the abstract (and by extension the manuscript) provides no quantitative error tables, convergence plots versus polynomial degree, or direct comparisons to known Schwarzschild results; without these, it is impossible to verify that the DG method retains exponential rates rather than degrading to algebraic convergence near the puncture.
  2. [Method] Method section on m-mode elliptic system: the regularity assumption underlying the exponential DG convergence—that the solutions of the m-mode PDEs remain sufficiently smooth away from the particle after null slicing and puncture subtraction—requires explicit demonstration for a=0.998, where ergosphere and near-horizon effects could reduce regularity; standard DG theory predicts only algebraic rates if the remainder is not C^infty, and this step is not secured by the provided description.
minor comments (2)
  1. [Abstract] The abstract states that 20 m-modes are solved in parallel in a few seconds; a minor clarification on wall-clock scaling with spin and orbit radius would help readers assess practicality.
  2. [Introduction] Notation for the vtu slicing and horizon-penetrating coordinates should be defined at first use to avoid ambiguity for readers outside numerical relativity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We address the two major comments point by point below, providing clarifications from the full manuscript and indicating revisions we will make to improve substantiation of the exponential convergence results.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that exponential convergence is achieved for a=0.998 at the ISCO is load-bearing for the paper's central result, yet the abstract (and by extension the manuscript) provides no quantitative error tables, convergence plots versus polynomial degree, or direct comparisons to known Schwarzschild results; without these, it is impossible to verify that the DG method retains exponential rates rather than degrading to algebraic convergence near the puncture.

    Authors: The manuscript body contains the requested quantitative evidence: convergence plots versus polynomial degree appear in Section 5.1 (Figures 7–9), including explicit exponential fits for a=0.998 at the ISCO; error tables are given in Table 3; and direct comparisons to published Schwarzschild scalar self-force values are discussed in Section 4.2. We agree the abstract itself lacks these details and will revise it to include a quantitative statement such as “achieving exponential convergence with relative errors below 10^{-9} at polynomial degree 6.” We will also add a parenthetical reference to the relevant figures and table. revision: partial

  2. Referee: [Method] Method section on m-mode elliptic system: the regularity assumption underlying the exponential DG convergence—that the solutions of the m-mode PDEs remain sufficiently smooth away from the particle after null slicing and puncture subtraction—requires explicit demonstration for a=0.998, where ergosphere and near-horizon effects could reduce regularity; standard DG theory predicts only algebraic rates if the remainder is not C^infty, and this step is not secured by the provided description.

    Authors: We acknowledge that the current Methods description does not explicitly analyze the smoothness of the m-mode remainder for a=0.998. The numerical results in Section 5 nevertheless demonstrate clear exponential (not algebraic) convergence at this spin, consistent with the expected DG rate. We will add a short paragraph to the m-mode elliptic system subsection explaining the regularity properties preserved by the vtu slicing and puncture subtraction, supported by references to existing analyses of Kerr regularity. A complete mathematical proof of C^infty smoothness lies beyond the scope of this numerical-methods paper. revision: partial

Circularity Check

0 steps flagged

Numerical method for elliptic PDEs with reported convergence exhibits no circularity

full rationale

The paper derives elliptic PDEs via m-mode separation and null slicing, then solves them numerically with high-order DG, AMR, and iterative solvers. Reported exponential convergence is an output of the discretization applied to the PDE system, not a quantity fitted or defined in terms of itself. No load-bearing self-citations, no parameters renamed as predictions, and no ansatze or uniqueness theorems that reduce the central claim to prior author work. The method and results are self-contained against external benchmarks (numerical convergence tests).

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The contribution is a numerical implementation; no new physical free parameters, ad-hoc axioms, or invented entities are introduced beyond the standard self-force perturbation framework in Kerr spacetime.

axioms (1)
  • domain assumption The scalar self-force can be obtained from the solution of the wave equation with a point source in Kerr spacetime using m-mode decomposition.
    The entire computational strategy rests on this separation being valid and on the resulting elliptic problems being well-posed.

pith-pipeline@v0.9.1-grok · 5792 in / 1395 out tokens · 35448 ms · 2026-06-28T05:15:49.918047+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

96 extracted references · 4 canonical work pages

  1. [1]

    There is no known spheroidal equivalent of a tensor- spherical harmonic basis in which the metric per- turbation equations separate. The direct metric perturbation problem can be formulated as a sys- tem of equations in which all l-modes are coupled together, which is an obstacle that would affect the field equations at both first and second order

  2. [2]

    puncture

    The poor convergence arising from the non- smoothness of the second-order source is exacer- bated by the increased mode coupling. Although it is possible to formulate the source construction as a Clebsch-Gordon mode-coupling problem (see, for example, Ref. [ 28] for a demonstration in the case of the Teukolsky equation), this relies on a projection onto s...

  3. [3]

    We attempted an associated numerical implementation using our iterative discretized elliptic PDE solver (see Sec

    Hyperboloidal andvtuslicing The formulation to this point, based on t-slicing and ˜Ψm, is nearly equivalent to that of Osburn and Nishimura [45] (other than using z rather than θ). We attempted an associated numerical implementation using our iterative discretized elliptic PDE solver (see Sec. IV), but that early version of our code required a prohibitive...

  4. [4]

    This brings several advantages including more efficient discretizations (by focusing resolution where it is most needed) and a simplified boundary treatment

    Horizon-penetrating coordinates and compactification The elimination of asymptotic oscillations via hyper- boloidal slicing enables compactification of the radial co- ordinate. This brings several advantages including more efficient discretizations (by focusing resolution where it is most needed) and a simplified boundary treatment. We also found that com...

  5. [5]

    Equations in first-order flux form Our final form of the field equation enables the use of DG methods by expressing the equations in first- order form with a divergence producing the principal part (see Sec. IV). To manipulate the field equations into this form, it is convenient to introduce the 2D coordinates qi ≡(r, z) and to define the 2D flux vector f...

  6. [6]

    The vanishing of F r m at the horizon means that we do not have to explicitly impose boundary conditions there

    Boundary conditions Our formulation of the field equation avoids coordinate singularities at the horizon and reduces boundary condi- tions to simple regularity conditions. The vanishing of F r m at the horizon means that we do not have to explicitly impose boundary conditions there. Similarly, factoring out the sin|m|θ behavior in Eq. (2.15) causes F z m ...

  7. [7]

    interior

    Note that the coefficients Aijk are functions of the particle position and four-velocity only, and that for our case of circular equatorial orbits we have that Aijk = 0 whenever eitherjorkis odd. B. Mode decomposition Next we want to calculate the m-mode punctures Ψ P m from ΦP according to the φ integral in Eq. (2.16) (note that the same logic defines Ψ ...

  8. [8]

    Domain decomposition Figure 1 gives an overview of our DG domain decompo- sition. We split the two-dimensional rectangular domain with coordinates qi = (r, z) (or qi = (r, z2) for equatorial symmetry) in 5 × 3 irreducibleblocks(or 5 × 2 for equato- rial symmetry) to define the different regions of our prob- lem: in the radial direction we have a horizon-p...

  9. [9]

    Specifically, our h refinement strategy is to split all elements containing the puncture in half in both dimensions in each refinement step

    Mesh refinement We split blocks into many smallerelements, Ω k, by repeatedly dividing the logical cube in half along either of the two dimensions ξ and η (h refinement). Specifically, our h refinement strategy is to split all elements containing the puncture in half in both dimensions in each refinement step. Since we place the particle at the center of ...

  10. [10]

    int” denotes the quantities on the interior side of the element, and “ext

    DG residuals Finally, we discretize the PDEs, △mΨm = Sm. To this end, we define nodal DG residuals on the grid by projecting the equations onto the same polynomial basis used to expand fields in each element, ψp(ξ). Project- ing the first-order flux divergence form of the equations, Eq. (2.40), and following Ref. [ 65], gives the DG residuals in strong fo...

  11. [11]

    (3.14), and between the v, t, and u domains, Eqs

    Imposing jump conditions With the DG numerical flux we also impose the jump conditions across the worldtube boundary, Eq. (3.14), and between the v, t, and u domains, Eqs. (2.36) and (2.37). To this end, we modify the exterior quantities Ψ ext m and (niF i m)ext once they were received by an element. To impose worldtube jumps, Eq. (3.14), we modify Ψext m...

  12. [12]

    Note that on boundaries where the perpendic- ular flux vanishes, niAi m = 0, neither of the two boundary terms contributes, as Ψ ext m is ignored and ( niF i m)ext = 0

    Imposing boundary conditions On external boundaries we impose boundary conditions through the DG numerical flux by a choice of Ψ ext m and (niF i m)ext. Note that on boundaries where the perpendic- ular flux vanishes, niAi m = 0, neither of the two boundary terms contributes, as Ψ ext m is ignored and ( niF i m)ext = 0. This occurs on the horizon boundary...

  13. [13]

    (4.3) with Gauss quadrature on the Legendre-Gauss collocation points of each element [75, 76]

    Matrix representation To obtain a discrete representation of the DG resid- uals, we evaluate the integrals in Eq. (4.3) with Gauss quadrature on the Legendre-Gauss collocation points of each element [75, 76]. This step reduces Eq. (4.3) to −M Di · F i m +M·β mΨm +γ i mM Di ·Ψ m +M L DT i · Aij mnj(Ψ∗ m −Ψ m) −M L· (niF i m)∗ −n iF i m =M·S m, (4.10) where...

    work page doi:10.3030/101200625

  14. [14]

    Bakeret al., The Laser Interferometer Space An- tenna: Unveiling the Millihertz Gravitational Wave Sky, arXiv:1907.06482 [astro-ph.IM] (2019)

    J. Bakeret al., The Laser Interferometer Space An- tenna: Unveiling the Millihertz Gravitational Wave Sky, arXiv:1907.06482 [astro-ph.IM] (2019)

    Pith/arXiv arXiv 1907

  15. [15]

    Colpiet al.(LISA), LISA Definition Study Report, arXiv:2402.07571 [astro-ph.CO] (2024)

    M. Colpiet al.(LISA), LISA Definition Study Report, arXiv:2402.07571 [astro-ph.CO] (2024)

    Pith/arXiv arXiv 2024

  16. [16]

    Afshordiet al.(LISA Consortium Waveform Work- ing Group), Waveform modelling for the Laser Interfer- ometer Space Antenna, Living Rev

    N. Afshordiet al.(LISA Consortium Waveform Work- ing Group), Waveform modelling for the Laser Interfer- ometer Space Antenna, Living Rev. Rel.28, 9 (2025), arXiv:2311.01300 [gr-qc]

    Pith/arXiv arXiv 2025

  17. [17]

    Abacet al.(ET), The Science of the Einstein Telescope, JCAP2026(03), 081, arXiv:2503.12263 [gr-qc]

    A. Abacet al.(ET), The Science of the Einstein Telescope, JCAP2026(03), 081, arXiv:2503.12263 [gr-qc]

    Pith/arXiv arXiv

  18. [18]

    Evanset al., A Horizon Study for Cosmic Explorer: Science, Observatories, and Community, arXiv:2109.09882 [astro-ph.IM] (2021)

    M. Evanset al., A Horizon Study for Cosmic Explorer: Science, Observatories, and Community, arXiv:2109.09882 [astro-ph.IM] (2021)

    Pith/arXiv arXiv 2021

  19. [19]

    C. P. L. Berry, S. A. Hughes, C. F. Sopuerta, A. J. K. Chua, A. Heffernan, K. Holley-Bockelmann, D. P. Mi- 17 haylov, M. C. Miller, and A. Sesana, The unique po- tential of extreme mass-ratio inspirals for gravitational- wave astronomy, Bull. Am. Astron. Soc.51, 42 (2019), arXiv:1903.03686 [astro-ph.HE]

    arXiv 2019

  20. [20]

    B. P. Abbottet al.(The LIGO Scientific Collaboration and Virgo Collaboration), GWTC-1: A Gravitational- Wave Transient Catalog of Compact Binary Mergers Ob- served by LIGO and Virgo during the First and Sec- ond Observing Runs, Phys. Rev. X9, 031040 (2019), arXiv:1811.12907 [astro-ph.HE]

    Pith/arXiv arXiv 2019

  21. [21]

    R. Abbottet al.(The LIGO Scientific Collaboration and Virgo Collaboration), GWTC-2: Compact Binary Coales- cences Observed by LIGO and Virgo During the First Half of the Third Observing Run, Phys. Rev. X11, 021053 (2021), arXiv:2010.14527 [gr-qc]

    Pith/arXiv arXiv 2021

  22. [22]

    R. Abbottet al.(The LIGO Scientific Collabora- tion, Virgo Collaboration, and KAGRA Collaboration), GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run, Phys. Rev. X13, 041039 (2023), arXiv:2111.03606 [gr-qc]

    Pith/arXiv arXiv 2023

  23. [23]

    M. A. Scheelet al., The SXS collaboration’s third catalog of binary black hole simulations, Class. Quant. Grav.42, 195017 (2025), arXiv:2505.13378 [gr-qc]

    arXiv 2025

  24. [24]

    Healy and C

    J. Healy and C. O. Lousto, Fourth RIT binary black hole simulations catalog: Extension to eccentric orbits, Phys. Rev. D105, 124010 (2022), arXiv:2202.00018 [gr-qc]

    arXiv 2022

  25. [25]

    K. Jani, J. Healy, J. A. Clark, L. London, P. Laguna, and D. Shoemaker, Georgia Tech Catalog of Gravita- tional Waveforms, Class. Quant. Grav.33, 204001 (2016), arXiv:1605.03204 [gr-qc]

    Pith/arXiv arXiv 2016

  26. [26]

    Fergusonet al., Second MAYA catalog of binary black hole numerical relativity waveforms, Phys

    D. Fergusonet al., Second MAYA catalog of binary black hole numerical relativity waveforms, Phys. Rev. D112, 044043 (2025), arXiv:2309.00262 [gr-qc]

    arXiv 2025

  27. [27]

    Hamilton, E

    E. Hamilton, E. Fauchon-Jones, M. Hannam, C. Hoy, C. Kalaghatgi, L. London, J. Thompson, D. Yeeles, S. Ghosh, S. Khan, P. Kolitsidou, and A. Vano-Vinuales, Precessing binary-black-hole numerical relativity cata- logue (complete data release), 10.5281/zenodo.7677297 (2023), dOI: 10.5281/zenodo.7677297

    work page doi:10.5281/zenodo.7677297 2023

  28. [28]

    Y. Mino, M. Sasaki, and T. Tanaka, Gravitational radia- tion reaction to a particle motion, Phys. Rev. D55, 3457 (1997), arXiv:gr-qc/9606018

    Pith/arXiv arXiv 1997

  29. [29]

    T. C. Quinn and R. M. Wald, Axiomatic approach to electromagnetic and gravitational radiation reaction of particles in curved spacetime, Phys. Rev. D56, 3381 (1997), arXiv:gr-qc/9610053

    Pith/arXiv arXiv 1997

  30. [30]

    Poisson, A

    E. Poisson, A. Pound, and I. Vega, The motion of point particles in curved spacetime, Living Rev. Rel.14, 7 (2011), arXiv:1102.0529 [gr-qc]

    Pith/arXiv arXiv 2011

  31. [31]

    Barack and A

    L. Barack and A. Pound, Self-force and radiation reaction in general relativity, Rept. Prog. Phys.82, 016904 (2019), arXiv:1805.10385 [gr-qc]

    Pith/arXiv arXiv 2019

  32. [32]

    Pound and B

    A. Pound and B. Wardell, Black hole perturbation the- ory and gravitational self-force, inHandbook of Gravi- tational Wave Astronomy, edited by C. Bambi, S. Kat- sanevas, and K. D. Kokkotas (Springer, Singapore, 2022) arXiv:2101.04592 [gr-qc]

    arXiv 2022

  33. [33]

    Pound, E

    A. Pound, E. Poisson, and B. G. Nickel, Limitations of the adiabatic approximation to the gravitational self-force, Phys. Rev. D72, 124001 (2005), arXiv:gr-qc/0509122

    Pith/arXiv arXiv 2005

  34. [34]

    Burke, G

    O. Burke, G. A. Piovano, N. Warburton, P. Lynch, L. Speri, C. Kavanagh, B. Wardell, A. Pound, L. Durkan, and J. Miller, Assessing the importance of first postadia- batic terms for small-mass-ratio binaries, Phys. Rev. D 109, 124048 (2024), arXiv:2310.08927 [gr-qc]

    arXiv 2024

  35. [35]

    Pound, B

    A. Pound, B. Wardell, N. Warburton, and J. Miller, Second-Order Self-Force Calculation of Gravitational Binding Energy in Compact Binaries, Phys. Rev. Lett. 124, 021101 (2020), arXiv:1908.07419 [gr-qc]

    arXiv 2020

  36. [36]

    Wardell, A

    B. Wardell, A. Pound, N. Warburton, J. Miller, L. Durkan, and A. Le Tiec, Gravitational Waveforms for Compact Binaries from Second-Order Self-Force Theory, Phys. Rev. Lett.130, 241402 (2023), arXiv:2112.12265 [gr-qc]

    arXiv 2023

  37. [37]

    Albertini, A

    A. Albertini, A. Nagar, A. Pound, N. Warburton, B. Wardell, L. Durkan, and J. Miller, Comparing second- order gravitational self-force, numerical relativity, and effective one body waveforms from inspiralling, quasicir- cular, and nonspinning black hole binaries, Phys. Rev. D 106, 084061 (2022), arXiv:2208.01049 [gr-qc]

    arXiv 2022

  38. [38]

    van de Meent, A

    M. van de Meent, A. Buonanno, D. P. Mihaylov, S. Os- sokine, L. Pompili, N. Warburton, A. Pound, B. Wardell, L. Durkan, and J. Miller, Enhancing the SEOBNRv5 effective-one-body waveform model with second-order gravitational self-force fluxes, Phys. Rev. D108, 124038 (2023), arXiv:2303.18026 [gr-qc]

    arXiv 2023

  39. [39]

    Spiers, A

    A. Spiers, A. Pound, and B. Wardell, Second-order per- turbations of the Schwarzschild spacetime: Practical, co- variant, and gauge-invariant formalisms, Phys. Rev. D 110, 064030 (2024), arXiv:2306.17847 [gr-qc]

    arXiv 2024

  40. [40]

    Miller, B

    J. Miller, B. Wardell, and A. Pound, Second-order per- turbation theory: the problem of infinite mode coupling, Phys. Rev. D94, 104018 (2016), arXiv:1608.06783 [gr-qc]

    Pith/arXiv arXiv 2016

  41. [41]

    Spiers, Analytically separating the source of the Teukolsky equation, Phys

    A. Spiers, Analytically separating the source of the Teukolsky equation, Phys. Rev. D109, 104059 (2024), arXiv:2402.00604 [gr-qc]

    arXiv 2024

  42. [42]

    S. R. Green, S. Hollands, and P. Zimmerman, Teukolsky formalism for nonlinear Kerr perturbations, Class. Quant. Grav.37, 075001 (2020), arXiv:1908.09095 [gr-qc]

    arXiv 2020

  43. [43]

    Toomani, P

    V. Toomani, P. Zimmerman, A. Spiers, S. Hollands, A. Pound, and S. R. Green, New metric reconstruc- tion scheme for gravitational self-force calculations, Class. Quant. Grav.39, 015019 (2022), arXiv:2108.04273 [gr-qc]

    arXiv 2022

  44. [44]

    Bourg, B

    P. Bourg, B. Leather, M. Casals, A. Pound, and B. Wardell, Implementation of a Green-Hollands- Zimmerman-Teukolsky puncture scheme for gravitational self-force calculations, Phys. Rev. D110, 044007 (2024), arXiv:2403.12634 [gr-qc]

    arXiv 2024

  45. [45]

    S. R. Dolan, C. Kavanagh, and B. Wardell, Gravitational Perturbations of Rotating Black Holes in Lorenz Gauge, Phys. Rev. Lett.128, 151101 (2022), arXiv:2108.06344 [gr-qc]

    arXiv 2022

  46. [46]

    S. R. Dolan, L. Durkan, C. Kavanagh, and B. Wardell, Metric perturbations of Kerr spacetime in Lorenz gauge: circular equatorial orbits, Class. Quant. Grav.41, 155011 (2024), arXiv:2306.16459 [gr-qc]

    arXiv 2024

  47. [47]

    Wardell, C

    B. Wardell, C. Kavanagh, and S. R. Dolan, Sourced metric perturbations of Kerr spacetime in Lorenz gauge, Class. Quant. Grav.42, 205007 (2025), arXiv:2406.12510 [gr-qc]

    arXiv 2025

  48. [48]

    Barack and D

    L. Barack and D. A. Golbourn, Scalar-field perturbations from a particle orbiting a black hole using numerical evolution in 2+1 dimensions, Phys. Rev. D76, 044020 (2007), arXiv:0705.3620 [gr-qc]

    Pith/arXiv arXiv 2007

  49. [49]

    Barack, D

    L. Barack, D. A. Golbourn, and N. Sago, m-Mode Regu- larization Scheme for the Self Force in Kerr Spacetime, Phys. Rev. D76, 124036 (2007), arXiv:0709.4588 [gr-qc]

    Pith/arXiv arXiv 2007

  50. [50]

    S. R. Dolan and L. Barack, Self force via m-mode regu- larization and 2+1D evolution: Foundations and a scalar- 18 field implementation on Schwarzschild, Phys. Rev. D83, 024019 (2011), arXiv:1010.5255 [gr-qc]

    Pith/arXiv arXiv 2011

  51. [51]

    S. R. Dolan, L. Barack, and B. Wardell, Self force via m-mode regularization and 2 + 1D evolution: II. Scalar- field implementation on Kerr spacetime, Phys. Rev. D84, 084001 (2011), arXiv:1107.0012 [gr-qc]

    Pith/arXiv arXiv 2011

  52. [52]

    S. R. Dolan and L. Barack, Self-force via m-mode regu- larization and 2+1D evolution: III. Gravitational field on Schwarzschild spacetime, Phys. Rev. D87, 084066 (2013), arXiv:1211.4586 [gr-qc]

    Pith/arXiv arXiv 2013

  53. [53]

    Thornburg and B

    J. Thornburg and B. Wardell, Scalar self-force for highly eccentric equatorial orbits in Kerr spacetime, Phys. Rev. D95, 084043 (2017), arXiv:1610.09319 [gr-qc]

    Pith/arXiv arXiv 2017

  54. [54]

    Miller and A

    J. Miller and A. Pound, Two-timescale evolution of extreme-mass-ratio inspirals: waveform generation scheme for quasicircular orbits in Schwarzschild spacetime, Phys. Rev. D103, 064048 (2021), arXiv:2006.11263 [gr-qc]

    arXiv 2021

  55. [55]

    Miller, B

    J. Miller, B. Leather, A. Pound, and N. Warburton, World- tube puncture scheme for first- and second-order self-force calculations in the Fourier domain, Phys. Rev. D109, 104010 (2024), arXiv:2401.00455 [gr-qc]

    arXiv 2024

  56. [56]

    Mathews and A

    J. Mathews and A. Pound, Postadiabatic waveform- generation framework for asymmetric precessing bina- ries, Phys. Rev. D112, 104078 (2025), arXiv:2501.01413 [gr-qc]

    arXiv 2025

  57. [57]

    Lewis, T

    J. Lewis, T. Kakehi, A. Pound, and T. Tanaka, Postadia- batic dynamics and waveform generation in self-force the- ory: An invariant pseudo-Hamiltonian framework, Phys. Rev. D113, 064046 (2026), arXiv:2507.08081 [gr-qc]

    arXiv 2026

  58. [58]

    Osburn and N

    T. Osburn and N. Nishimura, New self-force method via elliptic partial differential equations for Kerr inspiral mod- els, Phys. Rev. D106, 044056 (2022), arXiv:2206.07031 [gr-qc]

    arXiv 2022

  59. [59]

    Panosso Macedo, P

    R. Panosso Macedo, P. Bourg, A. Pound, and S. D. Upton, Multidomain spectral method for self-force calculations, Phys. Rev. D110, 084008 (2024), arXiv:2404.10083 [gr- qc]

    arXiv 2024

  60. [60]

    Bourg, A

    P. Bourg, A. Pound, S. D. Upton, and R. Panosso Macedo, Simple, efficient method of calculating the Detweiler- Whiting singular field to very high order, Phys. Rev. D 110, 084007 (2024), arXiv:2404.10082 [gr-qc]

    arXiv 2024

  61. [61]

    28th Capra Meeting on Radiation Reaction in General Relativity

    P. Bourg, The SS-piece of the second-order source (2025), talk at the “28th Capra Meeting on Radiation Reaction in General Relativity”

    2025

  62. [62]

    28th Capra Meeting on Radiation Reaction in General Relativity

    C. Dyson, Second-order source generation for m-modes in Kerr: Regular-regular and singular-regular contributions (2025), talk at the “28th Capra Meeting on Radiation Reaction in General Relativity”

    2025

  63. [63]

    I. Vega, P. Diener, W. Tichy, and S. L. Detweiler, Self- force with (3+1) codes: A Primer for numerical relativists, Phys. Rev. D80, 084021 (2009), arXiv:0908.2138 [gr-qc]

    Pith/arXiv arXiv 2009

  64. [64]

    Dhesi, H

    M. Dhesi, H. R. R¨ uter, A. Pound, L. Barack, and H. P. Pfeiffer, Worldtube excision method for intermediate- mass-ratio inspirals: Scalar-field toy model, Phys. Rev. D 104, 124002 (2021), arXiv:2109.03531 [gr-qc]

    arXiv 2021

  65. [65]

    N. A. Wittek, M. Dhesi, L. Barack, H. P. Pfeiffer, A. Pound, H. R. R¨ uter, M. S. Bonilla, N. Deppe, L. E. Kid- der, P. Kumar, M. A. Scheel, W. Throwe, and N. L. Vu, Worldtube excision method for intermediate-mass-ratio inspirals: Scalar-field model in 3 + 1 dimensions, Phys. Rev. D108, 024041 (2023), arXiv:2304.05329 [gr-qc]

    arXiv 2023

  66. [66]

    N. A. Wittek, A. Pound, H. P. Pfeiffer, and L. Barack, Worldtube excision method for intermediate-mass-ratio in- spirals: Self-consistent evolution in a scalar-charge model, Phys. Rev. D110, 084023 (2024), arXiv:2403.08864 [gr- qc]

    arXiv 2024

  67. [67]

    N. A. Wittek, L. Barack, H. P. Pfeiffer, A. Pound, N. Deppe, L. E. Kidder, A. Macedo, K. C. Nelli, W. Throwe, and N. L. Vu, Relieving Scale Disparity in Binary Black Hole Simulations, Phys. Rev. Lett.134, 251402 (2025), arXiv:2410.22290 [gr-qc]

    arXiv 2025

  68. [68]

    N. L. Vuet al., A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code, Phys. Rev. D105, 084027 (2022), arXiv:2111.06767 [gr- qc]

    arXiv 2022

  69. [69]

    Gourgoulhon, P

    E. Gourgoulhon, P. Grandcl´ ement, J.-A. Marck, J. Novak, and K. Taniguchi, LORENE, http://www.lorene.obspm. fr(2016)

    2016

  70. [70]

    Ansorg, B

    M. Ansorg, B. Bruegmann, and W. Tichy, A single-domain spectral method for black hole puncture data, Phys. Rev. D70, 064011 (2004), arXiv:gr-qc/0404056

    Pith/arXiv arXiv 2004

  71. [71]

    H. P. Pfeiffer, L. E. Kidder, M. A. Scheel, and S. A. Teukolsky, A multidomain spectral method for solving elliptic equations, Comput. Phys. Commun.152, 253 (2003), arXiv:gr-qc/0202096

    Pith/arXiv arXiv 2003

  72. [72]

    L. J. Papenfort, S. D. Tootle, P. Grandcl´ ement, E. R. Most, and L. Rezzolla, New public code for initial data of unequal-mass, spinning compact-object binaries, Phys. Rev. D104, 024057 (2021), arXiv:2103.09911 [gr-qc]

    arXiv 2021

  73. [73]

    Rashti, F

    A. Rashti, F. M. Fabbri, B. Br¨ ugmann, S. V. Chaurasia, T. Dietrich, M. Ujevic, and W. Tichy, New pseudospectral code for the construction of initial data, Phys. Rev. D 105, 104027 (2022), arXiv:2109.14511 [gr-qc]

    arXiv 2022

  74. [74]

    Ury¯ u and A

    K. Ury¯ u and A. Tsokaros, New code for equilibriums and quasiequilibrium initial data of compact objects, Phys. Rev. D85, 064014 (2012), arXiv:1108.3065 [gr-qc]

    Pith/arXiv arXiv 2012

  75. [75]

    Assumpcao, L

    T. Assumpcao, L. R. Werneck, T. P. Jacques, and Z. B. Etienne, NRPyElliptic: A fast hyperbolic relaxation el- liptic solver for numerical relativity, I: Conformally flat, binary puncture initial data, Phys. Rev. D105, 104037 (2022), arXiv:2111.02424 [gr-qc]

    arXiv 2022

  76. [76]

    Brandt and B

    S. Brandt and B. Br¨ ugmann, A simple construction of initial data for multiple black holes, Phys. Rev. Lett.78, 3606 (1997), arXiv:gr-qc/9703066

    Pith/arXiv arXiv 1997

  77. [77]

    Deppe, W

    N. Deppe, W. Throwe, L. E. Kidder, N. L. Vu, K. C. Nelli, C. Armaza, M. S. Bonilla, F. H´ ebert, Y. Kim, P. Kumar, G. Lovelace, A. Macedo, J. Moxon, E. O’Shea, H. P. Pfeiffer, M. A. Scheel, S. A. Teukolsky, N. A. Wittek, et al., SpECTRE v2024.09.29 , 10.5281/zenodo.13858965 (2024)

    work page doi:10.5281/zenodo.13858965 2024

  78. [78]

    N. L. Vu, Discontinuous Galerkin scheme for elliptic equa- tions on extremely stretched grids, Phys. Rev. D110, 084062 (2024), arXiv:2405.06120 [gr-qc]

    arXiv 2024

  79. [79]

    S. L. Detweiler and B. F. Whiting, Self-force via a Green’s function decomposition, Phys. Rev. D67, 024025 (2003), arXiv:gr-qc/0202086

    Pith/arXiv arXiv 2003

  80. [80]

    Vega and S

    I. Vega and S. L. Detweiler, Regularization of fields for self- force problems in curved spacetime: Foundations and a time-domain application, Phys. Rev. D77, 084008 (2008), arXiv:0712.4405 [gr-qc]

    Pith/arXiv arXiv 2008

Showing first 80 references.