Skip to content

Foliations, Lapse, and Shift

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

A spacelike foliation writes a spacetime region as a one-parameter family of spacelike hypersurfaces. The lapse tells how much proper time separates neighboring slices along their future unit normals; the shift tells how the spatial coordinates slide tangentially from one slice to the next. Together with the induced metric they reconstruct the spacetime metric kinematically.

Because the construction is geometric, one can change adapted coordinates and obtain different lapse and shift while preserving the same foliation and spacetime metric; their transformation must be followed together with the coordinate labels.

The split separates three choices that are often blended together: which hypersurfaces count as equal time, how quickly the time label advances between them, and how spatial points are identified from one slice to the next. The first fixes the normal direction, the lapse normalizes the label rate, and the shift fixes the tangential identification. None of these choices should be mistaken for an observable flow without an additional physical specification.

Required background. Affine Connections and Covariant Derivatives supplies covariant differentiation, while Flows and Lie Derivatives supplies the flow-based derivative used for evolution; The Riemann Curvature Tensor supplies later curvature projections; Pseudo-Riemannian Metrics and Signature supplies spacelike slices and timelike normals. Helpful background. Calculus and Linear Algebra for Relativity supplies block decompositions; Exterior Algebra and Differential Forms supplies the form and pullback language used for level-set time functions.

Let tt be a smooth time function whose level sets Σt\Sigma_t are spacelike. Their future-directed unit normal can be written

nμ=−N∇μt,N>0,n_\mu=-N\nabla_\mu t, \qquad N>0,

where normalization gives

N=(−gμν∇μt∇νt)−1/2.N=\left(-g^{\mu\nu}\nabla_\mu t\nabla_\nu t\right)^{-1/2}.

Choose coordinates xix^i on one slice and propagate their labels so that the coordinate evolution vector is tμ=(∂/∂t)μt^\mu=(\partial/\partial t)^\mu. Since tμ∇μt=1t^\mu\nabla_\mu t=1 and nμ∇μt=1/Nn^\mu\nabla_\mu t=1/N, the unique normal–tangential split is

tμ=Nnμ+βμ,nμβμ=0.t^\mu=Nn^\mu+\beta^\mu, \qquad n_\mu\beta^\mu=0.

The lapse NN is a scalar on each slice. The shift βμ\beta^\mu is a tangential vector. Neither is determined by the intrinsic metric alone; both describe how the foliation and coordinate labels are threaded through spacetime.

The vector triangle below makes the two roles geometrically distinct. Follow the normal leg from the lower leaf, then the tangential leg within the upper leaf, and compare their sum with the coordinate-time displacement.

Between two spacelike slices, a diagonal coordinate-time arrow is resolved into a normal arrow between slices followed by a tangent arrow along the upper slice.

A coordinate-time displacement is the sum of normal proper advance N n dt and tangential coordinate shift beta dt. A coordinate line is normal only when the shift vanishes. Schematic and not to scale.

In coordinates adapted to the foliation,

ds2=−N2dt2+hij(dxi+βidt)(dxj+βjdt).ds^2 =-N^2dt^2 +h_{ij}(dx^i+\beta^i dt)(dx^j+\beta^j dt).

Thus

gtt=−N2+hijβiβj,gti=hijβj,gij=hij.g_{tt}=-N^2+h_{ij}\beta^i\beta^j, \qquad g_{ti}=h_{ij}\beta^j, \qquad g_{ij}=h_{ij}.

The inverse components are

gtt=−1N2,gti=βiN2,gij=hij−βiβjN2.g^{tt}=-\frac1{N^2}, \qquad g^{ti}=\frac{\beta^i}{N^2}, \qquad g^{ij}=h^{ij}-\frac{\beta^i\beta^j}{N^2}.

Multiplying the two block matrices verifies these expressions. The determinant satisfies −g=Nh\sqrt{-g}=N\sqrt h. These algebraic checks catch the common mistake of omitting the shift-squared term in gttg_{tt}.

The spacetime volume element consequently factors as normal proper separation times spatial volume:

dVg=N dt dVh.dV_g=N\,dt\,dV_h.

This factorization is local and orientation sensitive. It is useful for integrals over a foliated region, but it does not assert that the coordinates cover the entire spacetime or that all slices are complete.

Using the convention

Kij=−hiμhjν∇μnν,K_{ij} =-h_i{}^\mu h_j{}^\nu\nabla_\mu n_\nu,

the normal derivative of the spatial metric is Lnhij=−2Kij\mathcal L_nh_{ij}=-2K_{ij}. Because tμ=Nnμ+βμt^\mu=Nn^\mu+\beta^\mu,

∂thij=−2NKij+Lβhij,\partial_t h_{ij} =-2NK_{ij}+\mathcal L_\beta h_{ij},

or

Kij=−12N(∂thij−Lβhij).K_{ij} =-\frac1{2N} \left(\partial_t h_{ij}-\mathcal L_\beta h_{ij}\right).

This identity says how the embedding changes under the chosen threading. It does not prescribe dynamics for hijh_{ij}.

The normal congruence generally accelerates. Projection of its derivative gives

ai=hiμnν∇νnμ=Diln⁡N.a_i=h_i{}^\mu n^\nu\nabla_\nu n_\mu=D_i\ln N.

A lapse constant on each slice makes the normal curves geodesic, although it may still vary from slice to slice. The shift does not enter this invariant acceleration; it changes coordinate threading rather than the normal itself. These relations are derived in Gourgoulhon 2007, §3.3, pp. 41–46 and §4.2, pp. 54–59.

Cartesian Minkowski slicing has N=1N=1, βi=0\beta^i=0, time-independent hijh_{ij}, and hence Kij=0K_{ij}=0. The same spacetime in uniformly translated spatial coordinates has a nonzero constant shift, but Lβh=0\mathcal L_\beta h=0 and KK remains zero.

Milne coordinates inside the future cone have

ds2=−dτ2+τ2γijdxidxj.ds^2=-d\tau^2+\tau^2\gamma_{ij}dx^idx^j.

Here N=1N=1, β=0\beta=0, and

Kij=−12∂τhij=−1τhij.K_{ij}=-\frac12\partial_\tau h_{ij} =-\frac1\tau h_{ij}.

This is an essential check: extrinsic curvature depends on the slicing even when the ambient Riemann tensor vanishes.

Rindler slicing,

ds2=−ρ2dη2+dρ2+dy2+dz2,ds^2=-\rho^2d\eta^2+d\rho^2+dy^2+dz^2,

has N=ρN=\rho, zero shift, and time-independent spatial metric, so Kij=0K_{ij}=0 while aρ=Dρln⁡ρ=1/ρa_\rho=D_\rho\ln\rho=1/\rho. This distinguishes acceleration of the normals from bending of the slices.

The decomposition may fail globally even when it works locally. A time function can cease to have spacelike level sets, a lapse can approach zero at a coordinate boundary, or normal curves can intersect. These are failures of a selected foliation or adapted chart unless an independent invariant diagnosis establishes something stronger.

Lapse is not generally proper time per coordinate time for every worldline. The statement applies along the normal displacement between slices. A coordinate line with nonzero shift is not normal.

Shift is not a spacetime vector independent of the foliation. It is the tangential part of the selected evolution field. Changing spatial coordinates changes its components and may change the field itself.

The metric-change equation is kinematic. Constraints, gauge choices for evolution, and numerical formulations require gravitational dynamics and lie outside this page.

Starting from the reconstructed line element, derive the future unit normal components in adapted coordinates.

Solution

Because nμ=−N∇μtn_\mu=-N\nabla_\mu t,

nμ=(−N,0,0,0).n_\mu=(-N,0,0,0).

Raising with the inverse metric gives

nt=1N,ni=−βiN.n^t=\frac1N, \qquad n^i=-\frac{\beta^i}{N}.

Then nμnμ=−1n^\mu n_\mu=-1, and

∂t=Nn+β\partial_t=Nn+\beta

follows component by component. This verifies both the shift sign in the line element and the normalization of the lapse.

  • Gourgoulhon, Éric. 3+1 Formalism and Bases of Numerical Relativity. Lecture notes, 2007. arXiv:gr-qc/0703035.