Pith. Sign in

REVIEW

Introduction to dynamical horizons in numerical relativity

General Relativity and Quantum Cosmology · Submitted 4 Apr 2006 · 95 references

Pith's verdict: A genuinely useful numerical implementation of dynamical horizon diagnostics, but the headline claim of coordinate-independent mass/spin/multipoles is only demonstrated in axisymmetry; the non-axisymmetric example explicitly disclaims those quantities.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDgr-qc/0604015 (version 2)
Pith IDARIXMIR2, Pith's permanent record of this paper and its review
CitationBadri Krishnan, Erik Schnetter, Florian Beyer. "Introduction to dynamical horizons in numerical relativity." arXiv:gr-qc/0604015 (2006). https://pith.science/paper/ARIXMIR2
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Dynamical horizons turn numerical simulations into gauge-free readouts of black hole mass, spin, and growth.

The reading

The paper argues that the dynamical horizon framework, which treats a black hole boundary as a stack of marginally trapped surfaces, can be turned into a practical diagnostic for numerical relativity. It claims that mass, angular momentum, source multipole moments, and energy fluxes can all be computed in a coordinate-independent way once the horizon world tube is located, and it implements the calculations in three simulations: head-on black hole collision, non-axisymmetric black hole inspiral, and neutron star collapse. In each case the method resolves how the hole forms, grows, and approaches a stationary Kerr state, and it exposes an inner horizon that behaves very differently from the outer one. A sympathetic reader would take the paper's central message to be that quasi-local horizons, not event horizons, are the right tool for reading off real-time black hole physics from a simulation.

What carries the argument

The central object is the marginally trapped tube (MTT): a smooth three-surface foliated by marginally outer trapped surfaces, each a closed spacelike 2-surface on which the outgoing null expansion vanishes. The signature of the tube is determined by the sign of $\det q_{ab}$, or equivalently by the sign of $\alpha = [\sigma_{(\ell)ab}\sigma^{ab}_{(\ell)} + R_{ab}\ell^a\ell^b]/\mathcal{L}_n\Theta_{(\ell)}$: spacelike for a dynamical horizon, null for an isolated horizon, timelike for a timelike membrane. This signature controls the area law, because spacelike tubes have increasing area along the outward direction while timelike tubes with $\Theta_{(n)}<0$ have decreasing area, and it feeds into the balance laws for energy and angular momentum. The paper computes the signature numerically from the induced metric of the world tube and uses a rotational vector field on each cross-section to define angular momentum, mass, and multipole moments.

Extended reading notes

Core claim

The central claim is that an evolving black hole's world tube of marginally outer trapped surfaces, called a marginally trapped tube, carries enough gauge-invariant information to define the hole's mass, angular momentum, multipole moments, and energy flux, and that the tube's causal character can be read off from the sign of the determinant of its induced three-metric. Applied to three simulations, the tube is null for initially isolated individual horizons, spacelike for the newly formed outer common horizon, and becomes timelike for the inner common horizon, whose area decreases monotonically. The outer horizon's mass multipoles decay toward Schwarzschild values in the head-on case, while in the neutron-star collapse the mass quadrupole and angular-momentum octupole approach their Kerr values, indicating that the hole settles down to Kerr. The paper treats these results as evidence that the dynamical horizon formalism provides a coherent, coordinate-invariant picture of black hole growth and relaxation.

Load-bearing premise

The argument assumes, without proof, that the inner marginally trapped surface stays smooth as it evolves; only the outer surface is covered by an existence theorem.

Editorial extensions

If this is right

  • The outer common MTT is spacelike and its area increases while it accretes matter and radiation, giving a local second law for black holes that event horizons cannot provide.
  • The inner MTT is generically not strictly-stably-outermost; it can be timelike and its area decreases, so the trapped region's world tube has richer structure than a single one-way membrane.
  • The multipole moments computed on the horizon can be compared directly with Kerr values, providing a gauge-invariant test of whether and how fast a black hole settles down.
  • The energy-flux balance law $(R_2 - R_1)/2 = F^{(m)} + F^{(g)}$ gives an exact, quasi-local analogue of the Bondi mass-loss formula at the horizon.
  • Because the formalism is slicing-dependent in which MTT it picks out, comparisons between simulations must specify the foliation, but once that is fixed the extracted quantities are coordinate-invariant.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper's tracking of inner MTTs is limited by surface finders that require star-shaped surfaces; a tracker allowing non-star-shaped or higher-genus surfaces could directly test whether inner tubes remain smooth or pinch off.
  • If the inner MTT's timelike signature and area decrease are generic, they may provide a quasi-local probe of the region near the singularity, analogous to the classical inner-horizon instability, without resolving the singularity itself.
  • The horizon multipole moments could be used to quantify the approach to Kerr by a single decay rate, such as the falloff of $M_2$, which might be compared with quasi-normal-mode frequencies extracted from gravitational waveforms.
Original abstract

This paper presents a quasi-local method of studying the physics of dynamical black holes in numerical simulations. This is done within the dynamical horizon framework, which extends the earlier work on isolated horizons to time-dependent situations. In particular: (i) We locate various kinds of marginal surfaces and study their time evolution. An important ingredient is the calculation of the signature of the horizon, which can be either spacelike, timelike, or null. (ii) We generalize the calculation of the black hole mass and angular momentum, which were previously defined for axisymmetric isolated horizons to dynamical situations. (iii) We calculate the source multipole moments of the black hole which can be used to verify that the black hole settles down to a Kerr solution. (iv) We also study the fluxes of energy crossing the horizon, which describes how a black hole grows as it accretes matter and/or radiation. We describe our numerical implementation of these concepts and apply them to three specific test cases, namely, the axisymmetric head-on collision of two black holes, the axisymmetric collapse of a neutron star, and a non-axisymmetric black hole collision with non-zero initial orbital angular momentum.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,290 input tokens, $0.0003 reported.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.09 · human participants 0.01 · controlled comparison 0.02 · randomized assignment 0.01 · states quantitative result 0.04 · compares with prior results 0.32 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

This paper presents a quasi-local method of studying the physics of dynamical black holes in numerical simulations. This is done within the dynamical horizon framework, which extends the earlier work on isolated horizons to time-dependent situations. In particular: (i) We locate various kinds of marginal surfaces and study their time evolution. An important ingredient is the calculation of the signature of the horizon, which can be either spacelike, timelike, or null. (ii) We generalize the calculation of the black hole mass and angular momentum, which were previously defined for axisymmetric isolated horizons to dynamical situations. (iii) We calculate the source multipole moments of the black hole which can be used to verify that the black hole settles down to a Kerr solution. (iv) We also study the fluxes of energy crossing the horizon, which describes how a black hole grows as it accretes matter and/or radiation. We describe our numerical implementation of these concepts and apply them to three specific test cases, namely, the axisymmetric head-on collision of two black holes, the axisymmetric collapse of a neutron star, and a non-axisymmetric black hole collision with non-zero initial orbital angular momentum.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 This paper presents a quasi-local method of studying the physics of dynamical black holes in numerical simulations. contribution 0.99 0.89 0.16 asserted 1.00
s1 This is done within the dynamical horizon framework, which extends the earlier work on isolated horizons to time-dependent situations. background 0.51 0.15 0.05 asserted 0.99
s2 In particular: (i) We locate various kinds of marginal surfaces and study their time evolution. contribution 0.97 0.94 0.05 asserted 0.98
s3 An important ingredient is the calculation of the signature of the horizon, which can be either spacelike, timelike, or null. method 0.77 0.47 0.05 asserted 0.98
s4 (ii) We generalize the calculation of the black hole mass and angular momentum, which were previously defined for axisymmetric isolated horizons to dynamical situations. contribution 1.00 0.92 0.07 asserted 1.00
s5 (iii) We calculate the source multipole moments of the black hole which can be used to verify that the black hole settles down to a Kerr solution. contribution 0.97 0.95 0.07 asserted 0.84
s6 (iv) We also study the fluxes of energy crossing the horizon, which describes how a black hole grows as it accretes matter and/or radiation. contribution 0.99 0.86 0.08 asserted 0.99
s7 We describe our numerical implementation of these concepts and apply them to three specific test cases, namely, the axisymmetric head-on collision of two black holes, the axisymmetric collapse of a neutron star, and a non-axisymmetric black hole collision with non-zero initial orbital angular momentum. method 0.61 0.96 0.09 asserted 1.00

What kind of paper the classifier reports this is: method or tool 0.93 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. shows also a Cauchy surface Σ intersecting a dynamical horizon H. This intersection S will always be assumed to be one of the given cross-sections of H. The unit timelike normal to the horizon is T a and the unit outward pointing spacelike normal to S within Σ is R a . The three metric and extrinsic curvature of Σ are denoted by q¯ab and K¯ ab respectively. The fiducial set of null normals to S arising naturally from Σ are ¯ℓ a = 1 √ 2 (T a + R a ), n¯ a = 1 √ 2 (T a − R a ). (2.8) A boost transformation of the form of equation (2.2) con￾nects (ℓ a , n a ) and
Figure 2. FIG. 2: Two MOTSs
Figure 3. FIG. 3: Coordinate shapes of the horizons at
Figure 4. FIG. 4: Determinant of the horizon world tube’s three-metri
Figure 5. FIG. 5: Irreducible mass vs. time for the individual and the
Figure 6. demonstrates how the common outer ap￾parent horizon grows. The energy flux vanishes at the poles, and the shear (but not the total flux) is maximum at the equator. The horizon is spacelike all the time, but it becomes increasingly isolated at late times as it ap￾proaches equilibrium. Thus the rate of area increase be￾comes smaller and the fluxes also becomes correspond￾ingly smaller. Let us now consider the higher mass multipoles Mn (all the Jns vanish identically). Here, since all quanti￾ties are symmetric with respect to a reflection about the equatorial plane, Mn = 0 for odd n
Figure 7. FIG. 7: Some mass multipole moments vs. time for the inner
Figure 8. FIG. 8: Coordinate shapes of the MOTSs at
Figure 9. FIG. 9: Determinant of the MTT three-metric at
Figure 11. FIG. 11: Energy flux through the common horizon and shear
Figure 12. FIG. 12: Lie derivative of the two-metric
Figure 13. FIG. 13: The average coordinate radius and the area radius as
Figure 14. FIG. 14: Average of the determinant of the horizon world
Figure 15. FIG. 15: The total mass

Showing 14 of the paper's 15 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:gr-qc/0604015 by the authors.

References
How the literature uses it

32 arXiv papers cite this in 53 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 29 sentences without a confident label

Uses its method or tool 15 papers
Cites it as background 8 papers
  • The problem of defining the analog of source multipole moments for black holes, which correspond to vacuum solutions ($T_ab=0$) of the Einstein equation, has only been tackled much more recently [ref] .

    Horizon Multipole Moments of a Kerr Black Hole · Motivation · 0.86

  • They are regularly used in numerical relativity simulations to compute physical quantities [ref] , and this formalism leads naturally to various versions of quasi-local black hole horizons.

    Self-intersecting marginally outer trapped surfaces · Marginally outer trapped surfaces · 0.99

  • We will mainly be interested in 2-surfaces $𝒮$ on which $θ_(l)=0$, called $marginally outer trapped surfaces$ (MOTSs) following the terminology in [ref] .

    Horizon dynamics of distorted rotating black holes · Basic Definitions and Concepts · 0.95

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.