Skip to content

Extrinsic Curvature and Metric Evolution

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.

Extrinsic curvature is the normal velocity of the spatial metric, with a sign fixed by convention. For the convention used here, expansion toward the future gives negative KijK_{ij}. Once lapse and shift are restored, the coordinate-time change of γij\gamma_{ij} separates cleanly into normal geometric evolution and a spatial Lie derivative caused only by relabeling.

Required background. Hypersurface Normals and Tangential Projectors supplies the induced connection and projections, while The Einstein Equation supplies the dynamics for which KijK_{ij} becomes initial and evolved data.

Helpful background. The Gibbons-Hawking-York Boundary Term shows why normal metric derivatives and extrinsic curvature enter the gravitational variational problem.

Extrinsic curvature as an embedding tensor

Section titled “Extrinsic curvature as an embedding tensor”

Let nμn^\mu be the future unit normal and γμν=δμν+nμnν\gamma^\mu{}_\nu=\delta^\mu{}_\nu+n^\mu n_\nu the spatial projector. This volume fixes

Kμν=−γμργνσ∇ρnσ.K_{\mu\nu} =-\gamma_\mu{}^\rho\gamma_\nu{}^\sigma\nabla_\rho n_\sigma.

Both indices are spatial. Hypersurface orthogonality makes KμνK_{\mu\nu} symmetric. For tangent vectors XμX^\mu and YμY^\mu, metric compatibility gives

(Lng)(X,Y)=g(∇Xn,Y)+g(X,∇Yn).(\mathcal L_n g)(X,Y) =g(\nabla_Xn,Y)+g(X,\nabla_Yn).

The normal terms in γμν=gμν+nμnν\gamma_{\mu\nu}=g_{\mu\nu}+n_\mu n_\nu vanish when evaluated on XX and YY. Symmetry of KK then yields

(Lnγ)(X,Y)=−2K(X,Y),(\mathcal L_n\gamma)(X,Y)=-2K(X,Y),

or equivalently

Kij=−12Lnγij.K_{ij}=-\frac12\mathcal L_n\gamma_{ij}.

This equation fixes the sign before any coordinate calculation. Sources that define Kij=+12LnγijK_{ij}=+\tfrac12\mathcal L_n\gamma_{ij} describe the same embedding with the opposite tensor; every term odd in KijK_{ij} must then be translated, not just the definition.

The trace and trace-free parts are

K=γijKij,Aij=Kij−13Kγij.K=\gamma^{ij}K_{ij}, \qquad A_{ij}=K_{ij}-\frac13K\gamma_{ij}.

The trace measures local volume expansion of the normal congruence, while AijA_{ij} measures shear of the slice geometry.

The coordinate evolution vector is fixed as

(∂t)μ=αnμ+βμ.(\partial_t)^\mu=\alpha n^\mu+\beta^\mu.

After projecting all free indices onto a slice, normal evolution is represented by ∂t−Lβ=αproj⁡ΣLn\partial_t-\mathcal L_\beta=\alpha\operatorname{proj}_{\Sigma}\mathcal L_n. For the spatial metric, with its indices pulled back to the slice, this gives the kinematic evolution equation

(∂t−Lβ)γij=−2αKij.(\partial_t-\mathcal L_\beta)\gamma_{ij} =-2\alpha K_{ij}.

Because Dkγij=0D_k\gamma_{ij}=0,

(Lβγ)ij=Diβj+Djβi.(\mathcal L_\beta\gamma)_{ij} =D_i\beta_j+D_j\beta_i.

Therefore the coordinate formula for extrinsic curvature is

Kij=−12α(∂tγij−Diβj−Djβi).K_{ij} =-\frac1{2\alpha} \left( \partial_t\gamma_{ij} -D_i\beta_j-D_j\beta_i \right).

The shift-dependent terms remove coordinate deformation within a slice. It is incorrect to infer nonzero extrinsic curvature from ∂tγij\partial_t\gamma_{ij} alone when the shift is nonzero.

Taking the trace gives the spatial-volume law. For the volume form ϵΣ=γ d3x\boldsymbol\epsilon_\Sigma=\sqrt\gamma\,d^3x,

(∂t−Lβ)ϵΣ=−αKϵΣ.(\partial_t-\mathcal L_\beta)\boldsymbol\epsilon_\Sigma =-\alpha K\boldsymbol\epsilon_\Sigma.

In a fixed coordinate chart this is

∂tln⁡γ=−αK+Diβi.\partial_t\ln\sqrt\gamma =-\alpha K+D_i\beta^i.

Thus K<0K<0 means that normal observers see the local spatial volume expand toward the future under the convention selected here. Gourgoulhon derives the Lie-derivative identity and its adapted-coordinate form in Gourgoulhon 2007, §3.3.5, pp. 45–46 and §4.3.1, pp. 59–60, PDF.

For inertial Minkowski slices,

α=1,βi=0,γij=δij,\alpha=1, \qquad \beta^i=0, \qquad \gamma_{ij}=\delta_{ij},

so ∂tγij=0\partial_t\gamma_{ij}=0 and hence Kij=0K_{ij}=0. The same conclusion holds for a rigidly rotating spatial grid: ∂tγij=0\partial_t\gamma_{ij}=0 and Lβγij=0\mathcal L_\beta\gamma_{ij}=0 because the rotation is a spatial Killing field.

By contrast, the Milne foliation of the same flat spacetime has

γij=τ2qij,α=1,βi=0,\gamma_{ij}=\tau^2q_{ij}, \qquad \alpha=1, \qquad \beta^i=0,

where qijq_{ij} is the unit hyperbolic metric. Direct differentiation gives

∂τγij=2τqij=2τγij,\partial_\tau\gamma_{ij}=2\tau q_{ij} =\frac2\tau\gamma_{ij},

and therefore

Kij=−1τγij,K=−3τ.K_{ij}=-\frac1\tau\gamma_{ij}, \qquad K=-\frac3\tau.

The volume element satisfies γ=τ3q\sqrt\gamma=\tau^3\sqrt q and ∂τln⁡γ=3/τ=−K\partial_\tau\ln\sqrt\gamma=3/\tau=-K, exactly as the trace equation requires. Nonzero KijK_{ij} here records how hyperboloids sit inside Minkowski spacetime; it is not four-dimensional curvature.

A complementary warning comes from static Schwarzschild-time slices. Their spatial metric is intrinsically curved, but with zero shift and no time dependence they have Kij=0K_{ij}=0. Intrinsic curvature and extrinsic curvature answer different questions.

Do not import the opposite sign selectively. Changing the sign of KijK_{ij} also changes the visible signs of momentum constraints, canonical momenta, junction formulas, and every term linear in KijK_{ij}.

Coordinate-time change is not normal change. The spatial Lie derivative Lβγij\mathcal L_\beta\gamma_{ij} must be removed before reading off the embedding velocity.

Nonzero extrinsic curvature does not imply curved spacetime. Milne hyperboloids have Kij≠0K_{ij}\ne0 in a flat ambient geometry.

Exercise: an expanding flat spatial metric

Section titled “Exercise: an expanding flat spatial metric”

Consider

ds2=−α(t)2dt2+a(t)2δijdxidxjds^2=-\alpha(t)^2dt^2+a(t)^2\delta_{ij}dx^idx^j

with zero shift. Derive KijK_{ij} and KK, and verify the volume-evolution equation. Do not set α=1\alpha=1 during the derivation.

Solution

The spatial metric is γij=a2δij\gamma_{ij}=a^2\delta_{ij}, so

∂tγij=2aa˙ δij=2a˙aγij.\partial_t\gamma_{ij} =2a\dot a\,\delta_{ij} =2\frac{\dot a}{a}\gamma_{ij}.

With βi=0\beta^i=0,

Kij=−12α∂tγij=−a˙αaγij,K=−3a˙αa.K_{ij} =-\frac1{2\alpha}\partial_t\gamma_{ij} =-\frac{\dot a}{\alpha a}\gamma_{ij}, \qquad K=-3\frac{\dot a}{\alpha a}.

Since γ=a3\sqrt\gamma=a^3,

∂tln⁡γ=3a˙a=−αK.\partial_t\ln\sqrt\gamma=3\frac{\dot a}{a} =-\alpha K.

The lapse converts coordinate-time expansion into expansion per unit normal proper time. For a(t)=ta(t)=t and α=1\alpha=1 with hyperbolic rather than Euclidean unit spatial metric, this sign agrees with the Milne result.

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