pith. sign in

arxiv: 2605.10874 · v2 · pith:NGRZKTVHnew · submitted 2026-05-11 · 🌀 gr-qc

Cusp Formation in Merging Black Hole Horizons

Pith reviewed 2026-05-20 21:54 UTC · model grok-4.3

classification 🌀 gr-qc
keywords black hole mergersapparent horizonscuspsmultipole momentsnumerical relativityhead-on collisionquasi-local horizons
0
0 comments X

The pith

Cusps form on black hole horizons during mergers and connect the initial black holes to the final remnant.

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

The paper investigates the formation of cusps on quasi-local horizons in the head-on collision of two non-spinning black holes using numerical methods. These cusps are shown to play a central role in linking the properties of the two separate initial black holes to those of the merged remnant black hole. The study details the evolution of the horizon mass and higher mass multipole moments specifically at the cusp locations. A phenomenological model is suggested to describe this behavior, offering an alternative to gravitational wave analysis for understanding merger outcomes.

Core claim

Cusps forming in otherwise smoothly evolving horizons play a central role in connecting the two initially separate black holes with the final remnant. For the head-on collision of two non-spinning black holes, the mass and higher mass multipole moments behave in specific ways at the cusp, which can be captured by a phenomenological model.

What carries the argument

Quasi-local black hole horizons and the cusps that form on them during the merger process.

Load-bearing premise

The numerical evolution accurately captures the quasi-local horizon geometry and the formation of cusps without significant gauge or resolution artifacts altering the multipole behavior.

What would settle it

If higher resolution simulations or different gauge choices reveal substantially different multipole moment values or behaviors at the cusp points, the reported results would be challenged.

Figures

Figures reproduced from arXiv: 2605.10874 by Badri Krishnan, Daniel Pook-Kolb, Shilpa Kastha, Stamatis Vretinaris.

Figure 1
Figure 1. Figure 1: FIG. 1: MOTS structure during the head-on collision. The second row shows close-ups of the configuration in the [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Mean curvature tr( [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Ricci scalar [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Maximum value of the Ricci scalar [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Integral of ∆ [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Difference of the Ricci scalar at the south pole of [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Ricci scalar [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Ricci scalar [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Fitted values for [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Fitted values for the scaling parameter [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Multiple moments of [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Even derivatives of the Ricci scalar [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: Ricci scalar [PITH_FULL_IMAGE:figures/full_fig_p018_13.png] view at source ↗
read the original abstract

An important question in binary black hole mergers is to connect properties of the remnant black hole to those of the two initial black holes. These properties include not only the final mass and spin of the remnant, but also higher multipoles and answers to other questions such as, for a given initial configuration, which quasi-normal modes of the final black hole are excited, and what are the amplitudes of these modes? Such questions have thus far been primarily addressed through a study of the emitted gravitational wave signal. In this paper we consider a different alternative, namely using quasi-local black hole horizons themselves to establish the link between the initial and final states. Recent work has elucidated the behavior of black hole horizons in a merger. Cusps forming in such otherwise smoothly evolving horizons have been shown to play a central role in connecting the two initially separate black holes with the final remnant. In the present work, we will discuss from a numerical perspective how such cusps form in detail for the head-on collision of two non-spinning black holes. We show how the mass and higher mass multipole moments behave at the cusp and suggest a phenomenological model.

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 / 0 minor

Summary. The paper presents a numerical study of the head-on collision of two non-spinning black holes, focusing on the formation of cusps in the evolving apparent horizons. It claims that these cusps play a central role in connecting the initial separate horizons to the final remnant black hole, and it examines the time evolution of the horizon mass and higher multipole moments at the cusp, proposing a phenomenological model to describe the behavior.

Significance. If the numerical results hold under rigorous validation, the work would offer a quasi-local horizon-based alternative to gravitational-wave analyses for linking initial and final black-hole properties, including potential insights into multipole evolution and quasi-normal mode excitation during mergers.

major comments (2)
  1. [Numerical results and cusp analysis] The manuscript provides no convergence tests, error bars, or explicit validation of the horizon finder and multipole extraction near cusp formation (as implied by the numerical perspective in the abstract). This is load-bearing for the central claim, since gauge drift or under-resolved curvature at the narrow bridge between horizons could alter the reported multipole time series by amounts comparable to the signal.
  2. [Phenomenological model] The phenomenological model is fitted directly to the extracted numerical multipoles without reported uncertainties, robustness checks against different gauges, or comparison to independent evolutions. This makes it unclear whether the model captures physical behavior or propagates numerical artifacts into the claimed connection between initial and remnant states.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting important issues regarding numerical validation and the phenomenological model. We address each major comment below and indicate the changes planned for the revised version.

read point-by-point responses
  1. Referee: [Numerical results and cusp analysis] The manuscript provides no convergence tests, error bars, or explicit validation of the horizon finder and multipole extraction near cusp formation (as implied by the numerical perspective in the abstract). This is load-bearing for the central claim, since gauge drift or under-resolved curvature at the narrow bridge between horizons could alter the reported multipole time series by amounts comparable to the signal.

    Authors: We agree that the absence of explicit convergence tests and error estimates weakens the presentation of the numerical results, especially near the cusp. The current manuscript emphasizes the qualitative formation process and multipole evolution but does not quantify numerical uncertainties. In the revision we will add convergence tests at multiple grid resolutions for both the horizon mass and the extracted multipole moments, report error bars based on the differences between resolutions, and include additional discussion of the horizon finder’s behavior in the bridge region to address possible gauge and resolution effects. revision: yes

  2. Referee: [Phenomenological model] The phenomenological model is fitted directly to the extracted numerical multipoles without reported uncertainties, robustness checks against different gauges, or comparison to independent evolutions. This makes it unclear whether the model captures physical behavior or propagates numerical artifacts into the claimed connection between initial and remnant states.

    Authors: We acknowledge that the model was presented without accompanying uncertainties or systematic checks. We will revise the relevant section to report fit uncertainties derived from the numerical data, add a brief robustness analysis with respect to small gauge variations, and clarify the model’s intended scope as a phenomenological description rather than a definitive physical law. Direct comparison with independent evolutions is not available for this specific head-on configuration; we will note this limitation explicitly and suggest it as a topic for future work. revision: partial

Circularity Check

0 steps flagged

No significant circularity in numerical horizon analysis

full rationale

The paper performs a numerical simulation of head-on non-spinning black hole mergers to describe cusp formation on apparent horizons and the associated evolution of mass and higher multipole moments, then proposes a phenomenological model based on those observed behaviors. No load-bearing step reduces a claimed result or prediction to a fitted parameter or self-citation by construction; the central claims rest on the output of the external numerical evolution code rather than re-deriving quantities from the same data via the paper's own equations. The work is self-contained against the simulation benchmarks and does not invoke uniqueness theorems or ansatze that collapse to prior author work.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The central claim rests on the assumption that the numerical relativity simulation faithfully represents quasi-local horizon geometry; no free parameters, axioms, or invented entities are explicitly introduced in the abstract.

pith-pipeline@v0.9.0 · 5735 in / 1125 out tokens · 39123 ms · 2026-05-20T21:54:22.557812+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

77 extracted references · 77 canonical work pages · 42 internal anchors

  1. [1]

    punctures

    for two non-spinning and uncharged black holes. Thus, we consider Euclidean spaceR 3 with two points (the “punctures”)x 1,2 excluded. Letδ ij be the Euclidean metric onR 3, anddthe Euclidean distance between the punctures:d=|x 1 −x 2|. The extrinsic curvature is trivial:K ij = 0, and the 3-metric is conformally flat: hij = Φ4δij. The conformal factor whic...

  2. [2]

    Thus, since these are non-spinning black holes, att touch, the mass ofS inner is obtained from the masses ofS 1 andS 2 as Min = q M 2 1 +M 2 2 .(23) This value ofM in can then be used as an initial value as we followS inner, initially backwards in time, through tbifurcate and then eventually towards the final black hole. It follows thatM in < M 1 +M 2 and...

  3. [3]

    For highly asymmet- ric systems, i.e

    (this is the same limit as obtained in [66]). For highly asymmet- ric systems, i.e. smallη, we will haveM in ≈M 1 +M 2. When we connect this to the final remnant mass (by trac- ingS inner backwardsin time till we reacht bifurcate, and then forward in time till we reach the final state of the remnant) the final mass will be larger than the above value ofM ...

  4. [4]

    The behavior ofS 1 andS 2: the mass isM 1 +M 2 which increases monotonically for−∞< t < t touch, as the two horizons absorb infalling radiation

  5. [5]

    Att touch: A discontinuous decrease in the mass when we go fromM 1 +M 2 →M in

  6. [6]

    north-pole

    Approach to equilibrium:M in generally increases till we reach the final equilibrium state (this in- crease is however not always monotonic [33]). The discontinuous jump in (2) above is larger (in abso- lute value) than the mass increases in (1) and (3), so that the final remnant mass is less than the sum of the initial masses. This is consistent with the...

  7. [7]

    THE MASS MUL TIPOLE MOMENTS A T THE MERGER In this section we present analytical results regarding the mass multipole moments ofS inner at the merger, i.e

    Thus: ζ⋆ = 2R2 1 R2 in −1 = 2 M 2 1 M 2 in −1 = M 2 1 −M 2 2 M 2 1 +M 2 2 .(29) V. THE MASS MUL TIPOLE MOMENTS A T THE MERGER In this section we present analytical results regarding the mass multipole moments ofS inner at the merger, i.e. att=t touch. A. The mass dipole moment at the merger The mass dipole momentM 1 (or alternatively the geo- metric momen...

  8. [8]

    Black hole fusion in the extreme mass ratio limit

    R. Emparan, M. Martinez, and M. Zilhao, Black hole fusion in the extreme mass ratio limit, Phys. Rev. D97, 044004 (2018), arXiv:1708.08868 [gr-qc]

  9. [9]

    The numerical relativity breakthrough for binary black holes

    U. Sperhake, The numerical relativity breakthrough for binary black holes, Class. Quant. Grav.32, 124011 (2015), arXiv:1411.3997 [gr-qc]

  10. [10]

    Black-hole binaries, gravitational waves, and numerical relativity

    J. Centrella, J. G. Baker, B. J. Kelly, and J. R. van Meter, Black-hole binaries, gravitational waves, and nu- merical relativity, Rev. Mod. Phys.82, 3069 (2010), arXiv:1010.5260 [gr-qc]

  11. [11]

    Massive Black Hole Binary Evolution

    D. Merritt and M. Milosavljevi´ c, Massive Black Hole Bi- nary Evolution, Living Reviews in Relativity8, 8 (2005), arXiv:astro-ph/0410364 [astro-ph]

  12. [12]

    C. L. Rodriguez, M. Zevin, P. Amaro-Seoane, S. Chat- terjee, K. Kremer, F. A. Rasio, and C. S. Ye, Black holes: The next generation—repeated mergers in dense star clusters and their gravitational-wave properties, Phys. Rev. D100, 043027 (2019)

  13. [13]

    Dreyer, B

    O. Dreyer, B. J. Kelly, B. Krishnan, L. S. Finn, D. Gar- rison, and R. Lopez-Aleman, Black hole spectroscopy: Testing general relativity through gravitational wave ob- servations, Class. Quant. Grav.21, 787 (2004), arXiv:gr- qc/0309007

  14. [14]

    Spectroscopy of Kerr black holes with Earth- and space-based interferometers

    E. Berti, A. Sesana, E. Barausse, V. Cardoso, and K. Bel- czynski, Spectroscopy of Kerr black holes with Earth- and space-based interferometers, Phys. Rev. Lett.117, 101102 (2016), arXiv:1605.09286 [gr-qc]

  15. [15]

    Modeling the remnant mass, spin, and recoil from unequal-mass, precessing black-hole binaries: The Intermediate Mass Ratio Regime

    Y. Zlochower and C. O. Lousto, Modeling the rem- nant mass, spin, and recoil from unequal-mass, precess- ing black-hole binaries: The Intermediate Mass Ratio Regime, Phys. Rev. D92, 024022 (2015), [Erratum: Phys.Rev.D 94, 029901 (2016)], arXiv:1503.07536 [gr-qc]

  16. [16]

    Remnant mass, spin, and recoil from spin aligned black-hole binaries

    J. Healy, C. O. Lousto, and Y. Zlochower, Remnant mass, spin, and recoil from spin aligned black-hole binaries, Phys. Rev. D90, 104004 (2014), arXiv:1406.7295 [gr-qc]

  17. [17]

    C. O. Lousto and Y. Zlochower, Black hole binary rem- nant mass and spin: A new phenomenological formula, Phys. Rev. D89, 104052 (2014), arXiv:1312.5775 [gr-qc]

  18. [18]

    The final spin from binary black holes in quasi-circular orbits

    F. Hofmann, E. Barausse, and L. Rezzolla, The final spin from binary black holes in quasi-circular orbits, Astro- phys. J. Lett.825, L19 (2016), arXiv:1605.01938 [gr-qc]

  19. [19]

    The final mass and spin of black hole mergers

    W. Tichy and P. Marronetti, The Final mass and spin of black hole mergers, Phys. Rev. D78, 081501 (2008), arXiv:0807.2985 [gr-qc]

  20. [20]

    Predicting the direction of the final spin from the coalescence of two black holes

    E. Barausse and L. Rezzolla, Predicting the direction of the final spin from the coalescence of two black holes, Astrophys. J. Lett.704, L40 (2009), arXiv:0904.2577 [gr- qc]

  21. [21]

    Modelling the final state from binary black-hole coalescences

    L. Rezzolla, Modelling the final state from binary black- hole coalescences, Class. Quant. Grav.26, 094023 (2009), arXiv:0812.2325 [gr-qc]. 18 4.5 5.0 5.5 6.0 t/M 1.25 1.50 1.75 2.00 2.25 R ttouch 4.5 5.0 5.5 6.0 t/M −5 −4 −3 −2 ∂2 s R ttouch 4.5 5.0 5.5 6.0 t/M 20 40 60 ∂4 s R ttouch 4.5 5.0 5.5 6.0 t/M −1.5 −1.0 −0.5 ∂6 s R ×103 ttouch FIG. 12: Even der...

  22. [22]

    The final spin from the coalescence of aligned-spin black-hole binaries

    L. Rezzolla, P. Diener, E. N. Dorband, D. Pollney, C. Reisswig, E. Schnetter, and J. Seiler, The Final spin from the coalescence of aligned-spin black-hole binaries, Astrophys. J. Lett.674, L29 (2008), arXiv:0710.3345 [gr- qc]

  23. [23]

    High-accuracy mass, spin, and recoil predictions of generic black-hole merger remnants

    V. Varma, D. Gerosa, L. C. Stein, F. H´ ebert, and H. Zhang, High-accuracy mass, spin, and recoil predic- tions of generic black-hole merger remnants, Phys. Rev. Lett.122, 011101 (2019), arXiv:1809.09125 [gr-qc]

  24. [24]

    The spin expansion for binary black hole merger: new predictions and future directions

    L. Boyle and M. Kesden, The spin expansion for binary black hole merger: new predictions and future directions, Phys. Rev. D78, 024017 (2008), arXiv:0712.2819 [astro- ph]

  25. [25]

    Binary black hole merger: symmetry and the spin expansion

    L. Boyle, M. Kesden, and S. Nissanke, Binary black hole merger: Symmetry and the spin expansion, Phys. Rev. Lett.100, 151101 (2008), arXiv:0709.0299 [gr-qc]

  26. [26]

    Cao and C.-Z

    Z.-J. Cao and C.-Z. Liu, Spin expansion for head-on col- lision of two black holes, Int. J. Mod. Phys. D20, 43 (2011)

  27. [27]

    Final spin of a coalescing black-hole binary: an Effective-One-Body approach

    T. Damour and A. Nagar, Final spin of a coalesc- ing black-hole binary: An Effective-one-body approach, Phys. Rev. D76, 044003 (2007), arXiv:0704.3550 [gr-qc]

  28. [28]

    Merger states and final states of black hole coalescences: a numerical-relativity-assisted effective-one-body approach

    T. Damour, A. Nagar, and L. Villain, Merger states and final states of black hole coalescences: a numerical- relativity-assisted effective-one-body approach, Phys. Rev. D89, 024031 (2014), arXiv:1307.2868 [gr-qc]. 19

  29. [29]

    Isolated and dynamical horizons and their applications

    A. Ashtekar and B. Krishnan, Isolated and dynamical horizons and their applications, Living Rev. Rel.7, 10 (2004), arXiv:gr-qc/0407042

  30. [30]

    Fundamental properties and applications of quasi-local black hole horizons

    B. Krishnan, Fundamental properties and applications of quasi-local black hole horizons, Class. Quant. Grav.25, 114005 (2008), arXiv:0712.1575 [gr-qc]

  31. [31]

    Dynamical Horizons and their Properties

    A. Ashtekar and B. Krishnan, Dynamical horizons and their properties, Phys. Rev.D68, 104030 (2003), arXiv:gr-qc/0308033

  32. [32]

    Ashtekar and B

    A. Ashtekar and B. Krishnan, Quasi-Local Black Hole Horizons: Recent Advances (2025), arXiv:2502.11825 [gr- qc]

  33. [33]

    Black hole boundaries

    I. Booth, Black hole boundaries, Can. J. Phys.83, 1073 (2005), arXiv:gr-qc/0508107

  34. [34]

    T. A. Apostolatos, C. Cutler, G. J. Sussman, and K. S. Thorne, Spin induced orbital precession and its modula- tion of the gravitational wave forms from merging bina- ries, Phys. Rev. D49, 6274 (1994)

  35. [35]

    Introduction to Isolated Horizons in Numerical Relativity

    O. Dreyer, B. Krishnan, D. Shoemaker, and E. Schnetter, Introduction to isolated horizons in numerical relativity, Phys. Rev. D67, 024018 (2003), arXiv:gr-qc/0206008

  36. [36]

    D. A. B. Iozzoet al., Comparing Remnant Proper- ties from Horizon Data and Asymptotic Data in Nu- merical Relativity, Phys. Rev. D103, 124029 (2021), arXiv:2104.07052 [gr-qc]

  37. [37]

    The existence and stability of marginally trapped surfaces

    D. Pook-Kolb, O. Birnholtz, B. Krishnan, and E. Schnet- ter, Existence and stability of marginally trapped sur- faces in black-hole spacetimes, Phys. Rev. D99, 064005 (2019), arXiv:1811.10405 [gr-qc]

  38. [38]

    Pook-Kolb, O

    D. Pook-Kolb, O. Birnholtz, B. Krishnan, and E. Schnet- ter, Self-intersecting marginally outer trapped surfaces, Phys. Rev. D100, 084044 (2019)

  39. [39]

    kinematic fitting

    D. Pook-Kolb, O. Birnholtz, B. Krishnan, and E. Schnet- ter, MOTS Finder version 1.1 (2019), 10.5281/zen- odo.2591105

  40. [40]

    Pook-Kolb, O

    D. Pook-Kolb, O. Birnholtz, B. Krishnan, and E. Schnet- ter, Interior of a binary black hole merger, Phys. Rev. Lett.123, 171102 (2019)

  41. [41]

    Pook-Kolb, O

    D. Pook-Kolb, O. Birnholtz, J. L. Jaramillo, B. Krishnan, and E. Schnetter, Horizons in a binary black hole merger I: Geometry and area increase (2020), arXiv:2006.03939 [gr-qc]

  42. [42]

    Pook-Kolb, O

    D. Pook-Kolb, O. Birnholtz, J. L. Jaramillo, B. Kr- ishnan, and E. Schnetter, Horizons in a binary black hole merger II: Fluxes, multipole moments and stability (2020), arXiv:2006.03940 [gr-qc]

  43. [43]

    Booth, R

    I. Booth, R. A. Hennigar, and D. Pook-Kolb, Ulti- mate fate of apparent horizons during a binary black hole merger. I. Locating and understanding axisymmet- ric marginally outer trapped surfaces, Phys. Rev. D104, 084083 (2021), arXiv:2104.11343 [gr-qc]

  44. [44]

    Pook-Kolb, I

    D. Pook-Kolb, I. Booth, and R. A. Hennigar, Ultimate fate of apparent horizons during a binary black hole merger. II. The vanishing of apparent horizons, Phys. Rev. D104, 084084 (2021), arXiv:2104.11344 [gr-qc]

  45. [45]

    Pook-Kolb, R

    D. Pook-Kolb, R. A. Hennigar, and I. Booth, What Happens to Apparent Horizons in a Binary Black Hole Merger?, Phys. Rev. Lett.127, 181101 (2021), arXiv:2104.10265 [gr-qc]

  46. [46]

    Dynamics of marginally trapped surfaces in a binary black hole merger: Growth and approach to equilibrium

    A. Gupta, B. Krishnan, A. Nielsen, and E. Schnetter, Dy- namics of marginally trapped surfaces in a binary black hole merger: Growth and approach to equilibrium, Phys. Rev.D97, 084028 (2018), arXiv:1801.07048 [gr-qc]

  47. [47]

    Mourier, X

    P. Mourier, X. Jim´ enez Forteza, D. Pook-Kolb, B. Kr- ishnan, and E. Schnetter, Quasinormal modes and their overtones at the common horizon in a binary black hole merger, Phys. Rev. D103, 044054 (2021), arXiv:2010.15186 [gr-qc]

  48. [48]

    Khera, A

    N. Khera, A. Ribes Metidieri, B. Bonga, X. Jim´ enez Forteza, B. Krishnan, E. Poisson, D. Pook- Kolb, E. Schnetter, and H. Yang, Nonlinear Ringdown at the Black Hole Horizon, Phys. Rev. Lett.131, 231401 (2023), arXiv:2306.11142 [gr-qc]

  49. [49]

    Understanding the "anti-kick" in the merger of binary black holes

    L. Rezzolla, R. P. Macedo, and J. L. Jaramillo, Un- derstanding the ’anti-kick’ in the merger of binary black holes, Phys. Rev. Lett.104, 221101 (2010), arXiv:1003.0873 [gr-qc]

  50. [50]

    J. L. Jaramillo, R. P. Macedo, P. M¨ osta, and L. Rezzolla, Towards a cross-correlation approach to strong-field dy- namics in Black Hole spacetimes,Proceedings, Span- ish Relativity Meeting : Towards new paradigms. (ERE 2011): Madrid, Spain, August 29-September 2, 2011, AIP Conf. Proc.1458, 158 (2011), arXiv:1205.3902 [gr-qc]

  51. [51]

    J. L. Jaramillo, R. P. Macedo, P. M¨ osta, and L. Rezzolla, Black-hole horizons as probes of black-hole dynamics II: geometrical insights, Phys. Rev.D85, 084031 (2012), arXiv:1108.0061 [gr-qc]

  52. [52]

    Y. Chen, P. Kumar, N. Khera, N. Deppe, A. Dhani, M. Boyle, M. Giesler, L. E. Kidder, H. P. Pfeiffer, M. A. Scheel, and S. A. Teukolsky, Multipole moments on the common horizon in a binary-black-hole simulation, Phys- ical Review D106, 10.1103/physrevd.106.124045 (2022)

  53. [53]

    A generalized Damour-Navier-Stokes equation applied to trapping horizons

    E. Gourgoulhon, A Generalized Damour-Navier-Stokes equation applied to trapping horizons, Phys. Rev. D72, 104007 (2005), arXiv:gr-qc/0508003

  54. [54]

    Ashtekar, N

    A. Ashtekar, N. Khera, M. Kolanowski, and J. Lewandowski, Charges and fluxes on (perturbed) non- expanding horizons, JHEP02, 066, arXiv:2112.05608 [gr-qc]

  55. [55]

    Dynamical Black Holes: Approach to the Final State

    A. Ashtekar, M. Campiglia, and S. Shah, Dynamical Black Holes: Approach to the Final State, Phys. Rev. D88, 064045 (2013), arXiv:1306.5697 [gr-qc]

  56. [56]

    Local existence of dynamical and trapping horizons

    L. Andersson, M. Mars, and W. Simon, Local existence of dynamical and trapping horizons, Phys.Rev.Lett.95, 111102 (2005), arXiv:gr-qc/0506013 [gr-qc]

  57. [57]

    Stability of marginally outer trapped surfaces and existence of marginally outer trapped tubes

    L. Andersson, M. Mars, and W. Simon, Stability of marginally outer trapped surfaces and existence of marginally outer trapped tubes, Adv.Theor.Math.Phys. 12(2008), arXiv:0704.2889 [gr-qc]

  58. [58]

    Isolated, slowly evolving, and dynamical trapping horizons: geometry and mechanics from surface deformations

    I. Booth and S. Fairhurst, Isolated, slowly evolving, and dynamical trapping horizons: geometry and mechan- ics from surface deformations, Phys. Rev.D75, 084019 (2007), arXiv:gr-qc/0610032

  59. [59]

    A Fast Apparent-Horizon Finder for 3-Dimensional Cartesian Grids in Numerical Relativity

    J. Thornburg, A Fast apparent horizon finder for three- dimensional Cartesian grids in numerical relativity, Class. Quant. Grav.21, 743 (2004), arXiv:gr-qc/0306056

  60. [60]

    Event and Apparent Horizon Finders for 3+1 Numerical Relativity

    J. Thornburg, Event and Apparent Horizon Finders for 3+1 Numerical Relativity, Living Rev. Rel.10, 3 (2007), arXiv:gr-qc/0512169

  61. [61]

    The Merger of Small and Large Black Holes

    P. M¨ osta, L. Andersson, J. Metzger, B. Szil´ agyi, and J. Winicour, The Merger of Small and Large Black Holes, Class. Quant. Grav.32, 235003 (2015), arXiv:1501.05358 [gr-qc]

  62. [62]

    The time evolution of marginally trapped surfaces

    L. Andersson, M. Mars, J. Metzger, and W. Simon, The Time evolution of marginally trapped surfaces, Class. Quant. Grav.26, 085018 (2009), arXiv:0811.4721 [gr-qc]. 20

  63. [63]

    Ashtekar, J

    A. Ashtekar, J. Engle, T. Paw lowski, and C. Van Den Broeck, Multipole moments of isolated hori- zons, Class. Quant. Grav.21, 2549 (2004), arXiv:gr- qc/0401114

  64. [64]

    Gourgoulhon, A

    E. Gourgoulhon, A. L. Tiec, and M. Casals, Horizon mul- tipole moments of a kerr black hole, Phys. Rev. D (2026)

  65. [65]

    D. R. Brill and R. W. Lindquist, Interaction energy in geometrostatics, Phys. Rev.131, 471 (1963)

  66. [66]

    Pook-Kolb, O

    D. Pook-Kolb, O. Birnholtz, B. Krishnan, E. Schnetter, and V. Zhang, MOTS Finder version 1.3 (2019)

  67. [67]

    The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics

    F. L¨ offler, J. Faber, E. Bentivegna, T. Bode, P. Diener, R. Haas, I. Hinder, B. C. Mundim, C. D. Ott, E. Schnet- ter, G. Allen, M. Campanelli, and P. Laguna, The Ein- stein Toolkit: A Community Computational Infrastruc- ture for Relativistic Astrophysics, Class. Quantum Grav. 29, 115001 (2012), arXiv:1111.3344 [gr-qc]

  68. [68]

    EinsteinToolkit, Einstein Toolkit: Open software for rel- ativistic astrophysics,http://einsteintoolkit.org/

  69. [69]

    A single-domain spectral method for black hole puncture data

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

  70. [70]

    J. D. Brown, P. Diener, O. Sarbach, E. Schnetter, and M. Tiglio, Turduckening black holes: An Analytical and computational study, Phys. Rev. D79, 044023 (2009), arXiv:0809.3533 [gr-qc]

  71. [71]

    S. Husa, I. Hinder, and C. Lechner, Kranc: a Mathemat- ica application to generate numerical codes for tensorial evolution equations, Comput. Phys. Commun.174, 983 (2006), arXiv:gr-qc/0404023

  72. [72]

    Kranc, Kranc: Kranc assembles numerical code

  73. [73]

    S. W. Hawking, Gravitational radiation from colliding black holes, Phys. Rev. Lett.26, 1344 (1971)

  74. [74]

    Bhagavan, B

    S. Bhagavan, B. de Koning, S. Maddhashiya, and C. Rackauckas, Datainterpolations.jl: Fast interpolations of 1d data, Journal of Open Source Software9, 6917 (2024)

  75. [75]

    Borhanian, K

    S. Borhanian, K. G. Arun, H. P. Pfeiffer, and B. S. Sathyaprakash, Comparison of post-Newtonian mode amplitudes with numerical relativity simulations of bi- nary black holes, Class. Quant. Grav.37, 065006 (2020), arXiv:1901.08516 [gr-qc]

  76. [76]

    X. J. Forteza and P. Mourier, High-overtone fits to numerical relativity ringdowns: Beyond the dismissed n=8 special tone, Phys. Rev. D104, 124072 (2021), arXiv:2107.11829 [gr-qc]

  77. [77]

    Black hole spectroscopy: from theory to experiment

    E. Bertiet al., Black hole spectroscopy: from theory to experiment (2025), arXiv:2505.23895 [gr-qc]