Affine-Null Characteristic Formulation

• Affine-null formulation is a coordinate method where null hypersurfaces are parametrized by an affine parameter rather than the areal radius, ensuring regularity at vanishing null expansion.
• The method reorganizes Einstein's equations by introducing auxiliary variables to restore a sequential radial integration hierarchy when the standard Bondi-Sachs hierarchy breaks down.
• Tetrad-based and first-order reductions within this framework yield robust hyperbolic formulations, supporting applications in black hole, scalar collapse, and cosmological models.

Affine-null characteristic formulation is a characteristic formulation of Einstein’s equations in which null hypersurfaces are labeled by a null coordinate and the generators of those hypersurfaces are parametrized by an affine parameter rather than by the Bondi-Sachs areal radius. In the general-relativistic literature, this choice appears in several closely related settings: outgoing or ingoing single-null evolution for gravitational radiation and black-hole spacetimes (Winicour, 2013), characteristic initial-boundary value problems with two intersecting null hypersurfaces (Mädler, 2018), spherically symmetric collapse and black-hole interior evolution (Crespo et al., 2019), hyperbolicity analyses of Bondi-like systems including affine gauge (Gundlach, 2024), and conformal or compactified formulations reaching null infinity (Mädler et al., 26 Apr 2025). The central geometric point is that the affine parameter remains regular where the null expansion vanishes, so the areal radius becomes a metric function rather than a coordinate; the central analytic point is that the Einstein equations typically lose the naive Bondi-Sachs hierarchy and must be reorganized through auxiliary variables to recover a sequential radial integration scheme (Winicour, 2013).

1. Geometric definition and coordinate gauges

In the outgoing single-null affine-null gauge, the metric form used in characteristic relativity is

$$ds^2 = V\,du^2 + 2\,du\,dr - h_{ab}\left(d\theta^a+W^a du\right)\left(d\theta^b+W^b du\right),$$

with coordinate conditions

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$

This makes $u=\mathrm{const}$ a null hypersurface and $r$ an affine parameter along the outgoing null geodesics generated by $\partial_r$ (Ripley, 2021). In the metric formulation designed for characteristic numerical evolution, the affine-null coordinates are instead written as $x^\alpha=(u,\lambda,x^A)$, with metric

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$

and angular decomposition

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$

so that $r$ is no longer a coordinate but a metric function $r(u,\lambda,x^A)$ (Winicour, 2013).

A closely related double-affine characteristic gauge is used for two intersecting null hypersurfaces. With coordinates

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$0

the affine-null metric ansatz is

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$1

with

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$2

and

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$3

(Mädler, 2018). In this setting, $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$4 labels the family of null hypersurfaces $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$5, $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$6 is an affine parameter along the null generators of $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$7, and the null boundary $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$8 sits at $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$9 (Mädler, 2018).

In spherical symmetry, the affine-null metric takes the particularly simple form

$u=\mathrm{const}$0

(Crespo et al., 2019), or equivalently

$u=\mathrm{const}$1

(Mädler et al., 2024). In $u=\mathrm{const}$2-dimensional spherical symmetry, the corresponding Einstein-Maxwell-scalar formulation uses

$u=\mathrm{const}$3

with $u=\mathrm{const}$4 (Bridera et al., 11 Jun 2026).

The geometrical motivation is stated consistently across these formulations: $u=\mathrm{const}$5 is tied to surface area and can fail at vanishing expansion, whereas the affine parameter is tied to the geodesic parametrization of the null generators and remains regular until caustics form (Winicour, 2013). In cosmological past-null-cone applications, the diameter distance reaches a maximum at the apparent horizon and then decreases, so the same value of $u=\mathrm{const}$6 can correspond to multiple points on the cone, while the affine parameter remains a unique label along each null geodesic all the way through the refocussed region (Walt et al., 2011).

2. Relation to Bondi-Sachs and the loss of the naive hierarchy

The Bondi-Sachs formulation uses outgoing null hypersurfaces $u=\mathrm{const}$7 and areal radius $u=\mathrm{const}$8, with metric

$u=\mathrm{const}$9

and determinant condition

$r$0

(Winicour, 2013). Its main equations take the hierarchical radial form

$r$1

$r$2

$r$3

$r$4

so one can integrate radially in the sequence

$r$5

(Winicour, 2013).

The affine-null system is related to Bondi-Sachs by

$r$6

but the paper emphasizes that this is misleading at the PDE level because the hypersurfaces $r$7 const and $r$8 const are not the same (Winicour, 2013). In affine-null coordinates, the corresponding equations initially look like

$r$9

$\partial_r$0

$\partial_r$1

$\partial_r$2

The key difficulty is the third equation: it contains $\partial_r$3 on the left, so the system is not a pure radial hypersurface hierarchy in the same clean sense as Bondi-Sachs (Winicour, 2013).

This same obstruction reappears in later spherical and conformal formulations. In the affine-null spherical Einstein-scalar system, the raw field equations are

$\partial_r$4

$\partial_r$5

$\partial_r$6

(Crespo et al., 2019). In the conformal affine-null system, the specialized PDEs for $\partial_r$7 still contain a derivative $\partial_r$8, which destroys the straightforward Bondi-Sachs-like hierarchy (Mädler et al., 26 Apr 2025).

A common misconception is that replacing $\partial_r$9 by an affine parameter is merely a coordinate relabeling. The cited formulations argue otherwise: the change is geometrically natural, but analytically it alters the structure of the reduced Einstein system and requires a different hierarchy (Winicour, 2013).

3. Restoration of the characteristic hierarchy by auxiliary variables

The standard method for restoring a sequential radial integration scheme is the introduction of auxiliary variables. In the general metric formulation, the auxiliary variables are

$x^\alpha=(u,\lambda,x^A)$0

$x^\alpha=(u,\lambda,x^A)$1

$x^\alpha=(u,\lambda,x^A)$2

together with

$x^\alpha=(u,\lambda,x^A)$3

$x^\alpha=(u,\lambda,x^A)$4

$x^\alpha=(u,\lambda,x^A)$5

(Winicour, 2013). The main equations then become

$x^\alpha=(u,\lambda,x^A)$6

$x^\alpha=(u,\lambda,x^A)$7

$x^\alpha=(u,\lambda,x^A)$8

$x^\alpha=(u,\lambda,x^A)$9

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$0

so that once $$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$1 is known on a null hypersurface, one can integrate these equations in order to obtain

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$2

and then reconstruct $$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$3 from $$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$4 (Winicour, 2013).

In the two-null-boundary formulation, the reorganization uses

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$5

together with

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$6

(Mädler, 2018). The hypersurface hierarchy then becomes

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$7

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$8

$$ds^2 = -(V - g_{AB} W^A W^B)\,du^2 - 2\,du\,d\lambda - 2 g_{AB} W^B\,du\,dx^A + g_{AB}\,dx^A dx^B,$$9

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$0

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$1

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$2

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$3

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$4

with

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$5

(Mädler, 2018).

In spherical symmetry, the auxiliary variables simplify further. One set is

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$6

leading to

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$7

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$8

$$g_{AB} = r^2 h_{AB}, \qquad \det(h_{AB}) = \det(q_{AB}),$$9

(Crespo et al., 2019). For horizon-penetrating evolution, the regularized variables

$r$0

lead to the final regular system

$r$1

$r$2

$r$3

$r$4

(Crespo et al., 2019). The closely related apparent-horizon tracking formulation uses

$r$5

$r$6

$r$7

and yields

$r$8

$r$9

$r(u,\lambda,x^A)$0

$r(u,\lambda,x^A)$1

(Mädler et al., 2024).

In the conformal compactified Einstein-scalar system, the auxiliary fields are $r(u,\lambda,x^A)$2 and $r(u,\lambda,x^A)$3, with

$r(u,\lambda,x^A)$4

$r(u,\lambda,x^A)$5

$r(u,\lambda,x^A)$6

and the final hierarchy

$r(u,\lambda,x^A)$7

$r(u,\lambda,x^A)$8

$r(u,\lambda,x^A)$9

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$00

(Mädler et al., 26 Apr 2025).

A plausible implication is that the “affine-null characteristic formulation” is best understood not as a single fixed reduced system but as a family of null-coordinate reductions that share one invariant feature: the affine gauge is geometrically simple, whereas the PDE hierarchy must be reconstructed.

4. Hyperbolicity, well-posedness, and tetrad-based reformulations

The hyperbolicity question is a central technical issue. One line of work shows that the second-order Einstein equations in Bondi-like coordinates, including affine gauge, are not strongly hyperbolic (Gundlach, 2024). In the Bondi-like metric

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$01

affine gauge corresponds to prescribing $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$02, so the null generator

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$03

is affinely parametrized (Gundlach, 2024). The paper states that in the frozen-coefficient Minkowski limit, the problematic mode is

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$04

whose multiplicity exceeds the dimension of its eigenspace, so the second-order Bondi-like formulations are only weakly hyperbolic, not strongly hyperbolic (Gundlach, 2024). Giannakopoulos, Hilditch, and Zilhão are likewise cited as having shown that the standard metric-based Bondi / affine-null evolution systems are only weakly hyperbolic, not strongly hyperbolic (Ripley, 2021).

A distinct response is the tetrad-based Newman-Penrose reformulation in affine-null coordinates. For the outgoing case, the tetrad is chosen as

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$05

with gauge fixing

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$06

(Ripley, 2021). The evolved variables are

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$07

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$08

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$09

(Ripley, 2021). The principal part is block diagonal, with

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$10

and a $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$11 Weyl block $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$12 explicitly given in the paper (Ripley, 2021). The principal symbol is Hermitian, and for

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$13

the contraction with the principal matrix is positive definite, so the system is symmetric hyperbolic (Ripley, 2021).

The contrast is sharp. The second-order metric formulations in Bondi-like coordinates, including affine gauge, fail strong hyperbolicity (Gundlach, 2024). The tetrad-based affine-null Newman-Penrose/GHP system is symmetric hyperbolic (Ripley, 2021). This suggests that the hyperbolicity problem is not a generic obstruction to affine-null coordinates as such, but to particular variable choices and reductions. That suggestion is reinforced by the first-order symmetric hyperbolic formulation of the linearized vacuum Einstein equations in Bondi gauge about Schwarzschild, which establishes a well-posed characteristic initial-boundary value problem in an $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$14-type energy norm for data on an outgoing null cone and a timelike cylinder or an ingoing null cone (Gundlach, 2024).

5. Spherical symmetry, black holes, and matter-coupled extensions

Spherical symmetry has provided the main testing ground for the affine-null characteristic formulation. In vacuum, the two-intersecting-null-hypersurface framework yields a black-hole branch

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$15

identified as the double-null Israel black hole solution (Mädler, 2018). The horizon is at

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$16

and the singularity occurs where

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$17

(Mädler, 2018).

The same framework is used in spherical symmetry to recover Reissner-Nordström from the Einstein-Maxwell hierarchy with affine-null metric

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$18

and Maxwell gauge

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$19

(Gallo et al., 2021). In the non-extremal case, the direct coordinate-construction approach gives

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$20

with

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$21

and

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$22

(Gallo et al., 2021).

For dynamical scalar collapse, the affine-null system is especially well adapted to horizon penetration. The outgoing expansion is

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$23

so an apparent horizon occurs where

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$24

(Mädler et al., 2024). The paper derives analytically that the apparent horizon is spacelike: $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$25 (Mädler et al., 2024). It also shows that the final singularity traced by the caustics is spacelike (Mädler et al., 2024). Earlier numerical work had already emphasized that the affine parameter remains well behaved across the horizon and inside black-hole regions up to the final singularity (Crespo et al., 2019).

The spherically symmetric Einstein-massless-scalar system also yields asymptotic quantities in affine-null form. The Misner-Sharp mass is

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$26

and the Bondi mass is

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$27

(Crespo et al., 2019). The mass-loss law becomes

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$28

(Crespo et al., 2019).

Matter-coupled generalizations preserve the same hierarchical logic. In $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$29-dimensional spherical symmetry for the Einstein-Maxwell system with a charged massless complex scalar field, the Maxwell potential is chosen in radial gauge,

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$30

with enclosed charge

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$31

(Bridera et al., 11 Jun 2026). The auxiliary variables are

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$32

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$33

(Bridera et al., 11 Jun 2026). The hierarchical equations include

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$34

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$35

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$36

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$37

and scalar transport

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$38

(Bridera et al., 11 Jun 2026). As a consistency check, the scalar-free limit recovers both non-extremal and extremal Reissner-Nordström-Tangherlini branches directly from the hierarchy (Bridera et al., 11 Jun 2026).

6. Asymptotics, conformal compactification, and applications

Affine-null formulations are closely tied to null infinity because outgoing affine-null coordinates are natural for radiation reaching future null infinity and are useful in Cauchy-characteristic extraction (Ripley, 2021). The original motivation for the metric formulation was precisely that, in Cauchy-characteristic extraction, the inner worldtube supplied by a Cauchy code is generally not located on the $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$39 const grid of a Bondi-Sachs code, so one must first interpolate Cauchy data into an auxiliary affine-null system anyway (Winicour, 2013). Using affine-null directly can eliminate that extra step and reduce interpolation error (Winicour, 2013).

In the conformal compactification of the affine-null Einstein-scalar system, the physical coordinates are

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$40

with line element

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$41

(Mädler et al., 26 Apr 2025). The conformal factor and compactified coordinate are chosen as

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$42

so that $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$43 sits at $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$44 (Mädler et al., 26 Apr 2025). With

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$45

the unphysical line element becomes

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$46

(Mädler et al., 26 Apr 2025). The asymptotic expansion yields the Bondi mass

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$47

and the balance law

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$48

(Mädler et al., 26 Apr 2025). In an inertial conformal frame on $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$49,

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$50

(Mädler et al., 26 Apr 2025).

The same paper shows that the conformal hierarchy is exactly equivalent to the system obtained by compactifying the affine coordinate in the physical spacetime and introducing regularized variables

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$51

(Mädler et al., 26 Apr 2025). This equivalence is presented as an identity of PDE systems, differing only in interpretation: in the conformal picture the boundary is geometrically built in, whereas in the physical compactified picture it is represented by a finite grid point after renormalization (Mädler et al., 26 Apr 2025).

Cosmological applications exploit the same geometric advantage of affine parametrization. The affine-null cosmology metric is

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$52

with regularity conditions at the vertex

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$53

(Walt et al., 2011). The radial geodesic relation

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$54

shows how the affine parameter is tied to the null generator (Walt et al., 2011). The numerical calculations are reported to have the same stability and convergence properties as the standard characteristic formalism (Walt et al., 2011).

Finally, affine-null coordinates have also been used to derive stationary rotating solutions directly from the characteristic hierarchy. In the slowly rotating Kerr construction, the affine-null metric ansatz is

$g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$55

with $g_{rr}=0,\qquad g_{ra}=0,\qquad g_{ru}=1.$56 (Mädler et al., 2023). By solving the vacuum Einstein equations order by order in slow rotation, the final metric is shown to be equivalent to Kerr in the slowly rotating approximation after an explicit transformation from Boyer-Lindquist coordinates (Mädler et al., 2023).

Taken together, these results support a broad characterization. The affine-null characteristic formulation is a null-coordinate framework in which the affine parameter regularizes horizon and refocussing behavior, restores direct contact with null infinity and null boundaries, and admits hierarchical evolution only after suitable auxiliary-field reformulation (Winicour, 2013). The main controversy is not over its geometric value, which is used repeatedly across vacuum, scalar, Maxwell, cosmological, and black-hole problems [(Walt et al., 2011); (Mädler, 2018); (Crespo et al., 2019); (Bridera et al., 11 Jun 2026)], but over the PDE structure of specific reductions: raw second-order metric systems in Bondi-like gauges are only weakly hyperbolic (Gundlach, 2024), whereas tetrad-based or more elaborate first-order reductions can recover symmetric hyperbolicity or a controlled characteristic initial-boundary value problem (Ripley, 2021, Gundlach, 2024).