Darmois-Israel Junction Conditions

The Darmois–Israel conditions identify the weakest nonnull matching admitted by ordinary distributional Einstein gravity. If the induced metric is continuous, a finite jump in extrinsic curvature produces a Dirac-delta Einstein tensor rather than an ill-defined product of distributions. Its coefficient is exactly the surface stress-energy tensor. A vanishing coefficient recovers a regular Darmois join.

Required background. We use Hypersurface Normals and Tangential Projectors, The Einstein Equation, and Flat-Background Hilbert Stress-Energy from Metric Variation.

Helpful background. The Gibbons-Hawking-York Boundary Term shows why fixing the first fundamental form removes normal metric variations on each smooth nonnull boundary segment.

Gaussian normal form and singular curvature

Let the nonnull interface $\Sigma$ divide $M_-$ from $M_+$. Choose a signed Gaussian normal coordinate $\ell$ that increases from the minus side to the plus side. Its unit normal is $n^\mu=(\partial_\ell)^\mu$, with $n^\mu n_\mu=\epsilon$, where $\epsilon=+1$ for a timelike interface and $\epsilon=-1$ for a spacelike one. Near $\Sigma$,

$$ds^2=\epsilon\,d\ell^2+h_{ab}(\ell,y)dy^a dy^b.$$

All side limits are pulled back to one abstract copy of $\Sigma$, and the jump convention is

$$[Q]\equiv Q_+-Q_-.$$

The normal is one continuous coorientation, not the two outward normals of the regions. With the convention used throughout this volume,

$$K_{ab}=-e^\mu_a e^\nu_b\nabla_\mu n_\nu =-\frac12\partial_\ell h_{ab}.$$

Assume $[h_{ab}]=0$ while $K_{ab}$ has finite one-sided limits. The connection is then piecewise bounded. Differentiating its step discontinuity once produces a delta function, while the quadratic connection terms remain locally integrable. Keeping only the normal derivatives gives

$$\left(G_{ab}\right)_{\mathrm{sing}} =\epsilon\bigl([K_{ab}]-h_{ab}[K]\bigr)\delta(\ell), \qquad \left(G_{\ell\mu}\right)_{\mathrm{sing}}=0,$$

where $K=h^{ab}K_{ab}$. This short Gaussian-normal calculation is the distributional core of the junction law. Israel’s invariant construction and its Gaussian-coordinate translation appear in Israel 1966, §§ 3, 5, and 7, pp. 6–10. His paper defines extrinsic curvature with the opposite sign, so its Lanczos equation must be sign-translated before comparison with the equation above.

The junction equation

Write the singular matter tensor as

$$T^{\mu\nu}_{\mathrm{sing}} =S^{ab}e^\mu_a e^\nu_b\,\delta(\ell), \qquad S^{\mu\nu}n_\nu=0.$$

Equating delta coefficients in $G_{\mu\nu}=8\pi G T_{\mu\nu}$ yields

$$8\pi G\,S_{ab} =\epsilon\bigl([K_{ab}]-h_{ab}[K]\bigr).$$

This form includes the causal sign $\epsilon$. For the timelike world tube of an ordinary material shell, $\epsilon=+1$. Taking the three-dimensional trace and substituting it back gives the equivalent inverse relation

$$[K]=-4\pi G\,\epsilon S, \qquad [K_{ab}] =8\pi G\,\epsilon \left(S_{ab}-\frac12h_{ab}S\right).$$

The Darmois conditions are the zero-layer specialization:

$$[h_{ab}]=0, \qquad [K_{ab}]=0.$$

They do not require the bulk stress tensor itself to be continuous. The Gauss–Codazzi constraints show instead that a regular join makes its normal-normal and normal-tangential projections compatible. Tangential stress or density may jump without creating a delta layer, as happens at an idealized stellar surface.

The distributional Bianchi identity supplies an additional consistency condition. For a timelike shell with the stated continuous normal,

$$D_a S^a{}_b =-[T_{\mu\nu}n^\mu e^\nu_b].$$

Thus surface stress-energy is conserved only when no net bulk flux crosses the layer. This relation is not an independent evolution law for the embedding; it is a consequence of the junction equation together with bulk Einstein equations.

Checks and limits of the formalism

There are three immediate checks. First, if the metric is actually $C^1$ across $\Sigma$, then $[K_{ab}]=0$ and the surface tensor vanishes. Second, reversing the single coorientation exchanges the labels $+\leftrightarrow-$ and sends both $n$ and every $K^\pm$ through the corresponding sign changes; the geometric shell tensor is unchanged when the entire convention is transformed consistently. Third, the dimensions agree: $K$ has inverse-length units, so $G S$ has inverse-length units, or equivalently $S$ has units $1/(G\,\text{length})$, as a surface energy density should.

The law fails if one inserts two outward-normal extrinsic curvatures into the continuous-normal jump. One side must first be reversed. It also fails when $[h_{ab}]\ne0$: then the connection itself contains a delta and curvature formally contains products such as $\delta^2$. Finally, the unit-normal derivation does not have a smooth null limit. A null interface has a degenerate induced metric and requires transverse curvature instead.

Exercises

Consider a timelike interface at $\ell=0$ with

$$ds^2=d\ell^2-d\tau^2+a(\ell)^2(dx^2+dy^2), \qquad a(0)=1.$$

Let the one-sided slopes be $a'_-=A$ and $a'_+=B$. Compute the surface density and the two equal surface pressures in the orthonormal coordinates $(\tau,x,y)$.

Solution

Set $\Delta=B-A$. The site convention gives

$$[K_{\tau\tau}]=0, \qquad [K_{xx}]=[K_{yy}]=-\Delta, \qquad [K]=-2\Delta.$$

Because the interface is timelike, $\epsilon=+1$. The junction equation therefore gives

$$S_{\tau\tau} =\frac{0-(-1)(-2\Delta)}{8\pi G} =-\frac{\Delta}{4\pi G},$$

and

$$S_{xx}=S_{yy} =\frac{-\Delta-1(-2\Delta)}{8\pi G} =\frac{\Delta}{8\pi G}.$$

Hence the surface energy density measured by $\partial_\tau$ is $\sigma=S_{\tau\tau}=-\Delta/(4\pi G)$ and the pressures are $p_x=p_y=\Delta/(8\pi G)$. All three vanish when $A=B$, as the Darmois conditions require.

References

• Israel, Werner. “Singular Hypersurfaces and Thin Shells in General Relativity.” Il Nuovo Cimento B 44, no. 1 (1966): 1–14; erratum 48, no. 2 (1967): 463. doi:10.1007/BF02710419. Erratum doi:10.1007/BF02712210.