research32 min studyLesson 847 in Physics

Black-hole punctures and excision

Compare two singularity-handling strategies and interpret horizon mass, spin, and characteristic diagnostics.

On this page

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.

ψ=1+m2r,Rareal=ψ2r,rthroat=m/2 ⟹ Rareal=2m\psi=1+\frac{m}{2r},\qquad R_{\rm areal}=\psi^2r,\qquad r_{\rm throat}=m/2\ \Longrightarrow\ R_{\rm areal}=2m
These formulas apply to the time-symmetric isotropic Schwarzschild slice away from its puncture. The parameter m and both radii have length units; this is an exact special initial geometry.↔ Scroll wide equations horizontally.

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.

v0=−βs,v±=−βs±αfor the declared unit-normal toy characteristic familyv_0=-\beta^s,\qquad v_\pm=-\beta^s\pm\alpha\qquad\text{for the declared unit-normal toy characteristic family}
The outward normal here points from retained grid into excision. Positive speeds are outflow for a fixed boundary; this three-speed family is an explicit toy, not a complete formulation spectrum.↔ Scroll wide equations horizontally.

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.

Mirr=A16π,MH2=Mirr2+JH24Mirr2,χH=JHMH2M_{\rm irr}=\sqrt{\frac{A}{16\pi}},\qquad M_H^2=M_{\rm irr}^2+\frac{J_H^2}{4M_{\rm irr}^2},\qquad \chi_H=\frac{J_H}{M_H^2}
Area and angular momentum have length-squared units, masses have length units, and chi is dimensionless. The mass expression needs the stated uncharged Kerr or compatible isolated-horizon prescription.↔ Scroll wide equations horizontally.

Area and spin in a declared horizon prescription

For an uncharged axisymmetric horizon using the displayed mass prescription, let area be 64πL264\pi L^2 and angular momentum 4L24L^2. Find irreducible mass, horizon mass and dimensionless spin. Explain which quantities need numerical uncertainty estimates.

  1. Dividing area by sixteen pi gives 4L24L^2, so irreducible mass is 2L2L. This is an area-derived mass, not a coordinate radius.

  2. The rotational contribution to mass squared is 16L4/(4×4L2)=L216L^4/(4\times4L^2)=L^2. Thus horizon mass squared is 5L25L^2 and horizon mass is 5L≈2.236L\sqrt5L\approx2.236L.

  3. The dimensionless spin is 4/5=0.84/5=0.8. 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 result

The irreducible mass is 2 L, the prescribed horizon mass is square-root-five L, and the dimensionless spin is 0.8.

Explore the model

Interactive lab

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.

A root is not automatically a horizonBisection 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. Current readings: Midpoint after 0 refinements: 1.5; Bracket error bound: 1.5; Computed residual · x² − c: 0.25. All parameters are adjustable with the labeled controls below.Bisection on x squared minus 2: bracket from 0 to 3 after 0 refinements.00.751.52.253−2.360.192.745.297.84x · bracket stays nonnegativef(x) = x² − cn = 0 · [ 0, 3 ]
Quadratic / midpointCurrent bracketPositive root
Midpoint after 0 refinements
1.5
Bracket error bound
1.5
Computed residual · x² − c
0.25
t = 0 s

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.

Check your understanding

Why is placing an excision surface inside an apparent horizon insufficient by itself to justify omitting all boundary data?

Try it yourself

For the stated fixed-boundary toy family use α=0.5\alpha=0.5 and βs=−1.5\beta^s=-1.5. Find the three speeds and decide whether they are all outflow. Repeat with βs=0.2\beta^s=0.2. 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.

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.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.