Black-hole punctures and excision
Compare two singularity-handling strategies and interpret horizon mass, spin, and characteristic diagnostics.
On this page
The question for today
What you’ll understand
- Distinguish a puncture coordinate singularity from an excision boundary.
- Evaluate a horizon-area and spin diagnostic under its stated assumptions.
- Check outgoing characteristic signs at a toy excision surface.
A black-hole evolution must avoid letting the physical interior singularity destroy the calculation. Puncture and excision methods address that challenge through different coordinate and domain constructions. Neither removes the black hole from the physics. Their success depends on gauge, evolved variables, causal characteristics, and resolution. Horizon diagnostics then measure geometric properties of the resulting solution, with uncertainties and assumptions that are different from the coordinate position of the hole.
Understand what a puncture represents
In time-symmetric conformally flat Schwarzschild initial data, the conformal factor has a simple inverse-radius singular term. The zero-coordinate-radius puncture represents another asymptotic end of the initial spatial geometry rather than an ordinary grid point at the physical curvature singularity. Isotropic radius and areal radius are different: the latter includes the square of the conformal factor. Treating the coordinate puncture distance as a physical separation therefore gives a misleading interpretation.
Moving-puncture evolution uses suitably chosen variables, lapse and shift so that puncture coordinates can move while the numerically evolved fields remain manageable. Inverse conformal-factor variables help avoid directly representing a divergent logarithmic conformal factor. Initial regularity, stencil behavior and gauge choices still matter. This method does not simply assume that every field at the puncture is smooth to arbitrary order or that errors inside a coordinate horizon are automatically irrelevant to all evolved modes.
Make excision an actual characteristic outflow problem
Excision removes an interior region from the numerical domain and places a boundary within the black-hole region. Physical causal intuition motivates this strategy, but a successful PDE boundary also needs the characteristics of the chosen formulation and gauge. With the normal pointing out of the retained domain into the excised region, positive normal speeds flow out of the computational domain. Incoming fields would require appropriate boundary information and invalidate an assumption of pure outflow.
A coordinate apparent horizon is not by itself a guarantee that every gauge or constraint characteristic has the needed sign at the excision surface. The surface can move, so characteristic speeds must be measured relative to its motion. Coordinate control systems and boundary placement are consequently part of a practical method. The toy speeds below illustrate the sign test only; a real first-order reduction may have additional modes and parameters that must all be checked.
Measure horizon geometry and state the mass prescription
An apparent horizon is an outermost closed marginally outer trapped surface on a chosen slice; locating it is an elliptic surface problem and depends on the slicing. An event horizon instead uses global causal information. The area of a located surface and an angular momentum derived with an axial or appropriately approximated symmetry vector can support a quasi-local diagnostic. Finding a surface is therefore different from assigning a coordinate center and radius.
For an uncharged Kerr horizon, or a suitable canonical axisymmetric isolated-horizon mass prescription, area defines an irreducible mass and area together with angular momentum gives a Christodoulou-type mass. The dimensionless spin is angular momentum divided by mass squared. During a distorted dynamical phase, approximate symmetry, horizon finding, resolution and the adopted quasi-local definitions contribute uncertainty. A Kerr formula used diagnostically should not silently become an exact invariant characterization of every evolving surface.
Work it through
Area and spin in a declared horizon prescription
For an uncharged axisymmetric horizon using the displayed mass prescription, let area be and angular momentum . Find irreducible mass, horizon mass and dimensionless spin. Explain which quantities need numerical uncertainty estimates.
Dividing area by sixteen pi gives , so irreducible mass is . This is an area-derived mass, not a coordinate radius.
The rotational contribution to mass squared is . Thus horizon mass squared is and horizon mass is .
The dimensionless spin is . An actual measurement would propagate horizon-location, area quadrature, angular-momentum and approximate-symmetry errors, along with the assumptions behind assigning this horizon mass.
The irreducible mass is 2 L, the prescribed horizon mass is square-root-five L, and the dimensionless spin is 0.8.
The living notebook
Explore the model
Interactive lab · root-finding
A root is not automatically a horizon
Reset and set target one, upper bound two, then Step once; predict root one and zero bracket bound. Compare a fresh Reset run at target 2.25 and upper bound three; predict midpoint root 1.5 after one Step. This monotone scalar bisection has no null expansion, metric, surface shape, or horizon equation.
- Midpoint after 0 refinements
- 1.5
- Bracket error bound
- 1.5
- Computed residual · x² − c
- 0.25
Bisection of x² − c on [0,U], capped at 50 iterations for floating-point resolution. The error bound is the maximum distance from the displayed midpoint to either bracket endpoint (half the width in exact arithmetic). Compensated square evaluation preserves the bracket when ordinary multiplication rounds a residual to zero. A truly exact endpoint or midpoint collapses the bracket. This is one monotone quadratic, not a general root solver. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forA numerical root is meaningful only for its specified equation; this instrument does not find black-hole horizons.
Pause & reason
Check your understanding
Your turn
Try it yourself
For the stated fixed-boundary toy family use and . Find the three speeds and decide whether they are all outflow. Repeat with . The normal points from retained domain into excision.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Apply the same negative shift offset to the advective mode and the two lapse-offset modes. Negative speeds enter the retained computational domain.
Compare with the solution
With shift minus 1.5 the speeds are 1.5, 2 and 1, all positive and hence outflow for this toy family. With shift 0.2 they are minus 0.2, 0.3 and minus 0.7, so two modes are incoming. The latter case cannot use a pure-outflow assumption for this family. A real excision calculation must check all its additional characteristic fields and account for the surface's motion rather than transferring this toy result directly.
Keep the idea
Take it with you
- Punctures and excision handle the numerical domain in different ways.
- Excision needs characteristic outflow, including gauge and constraint fields.
- Horizon area and spin diagnostics require a declared geometric prescription and uncertainty.
Read further
Original sources and deeper explanations to follow at your own pace.
- Gourgoulhon, section 8.2.5, Schwarzschild puncture initial geometry
- Campanelli et al., Accurate Evolutions of Orbiting Black-Hole Binaries Without Excision
- Thornburg, Event and Apparent Horizon Finders for 3+1 Numerical Relativity, apparent-horizon and excision discussion
- Ashtekar and Krishnan, Isolated and Dynamical Horizons and Their Applications, section 4.1, canonical horizon mass and angular momentum
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.