Skip to content

Impulsive Waves and Thin Sources

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.

Thin shells and impulsive waves concentrate curvature on a hypersurface while leaving ordinary smooth Einstein geometry on either side. They are controlled singularities only when the matching variables and distributional products are well defined. A non-null shell is governed by a jump in extrinsic curvature; a null shell needs a different transverse geometry; an impulsive wave can also appear through a distributional Brinkmann profile whose interpretation must be tied to a regularization or continuous coordinate form.

Required background. Analysis for Mathematical Relativity supplies distributions and traces; The Einstein Equation supplies the singular field equation; PDEs and Green Functions for Gravitation supplies weak solutions.

Helpful background. Causality and Global Structure helps distinguish a wavefront from a spacetime boundary or singular endpoint.

Let a timelike hypersurface Σ\Sigma divide spacetime into regions M−M^- and M+M^+. Choose a spacelike unit normal ss pointing from M−M^- to M+M^+, so g(s,s)=+1g(s,s)=+1, and require the induced metric habh_{ab} to agree from both sides. For this page define

Kab±=haμhbν∇μ±sν,[Kab]=Kab+−Kab−.\mathcal K^{\pm}_{ab} =h_a{}^\mu h_b{}^\nu\nabla^{\pm}_\mu s_\nu, \qquad [\mathcal K_{ab}]=\mathcal K^+_{ab}-\mathcal K^-_{ab}.

Here, K\mathcal K uses the displayed plus sign and is not the spacelike-slice KijK_{ij} of the volume convention.

Continuity of hh keeps the metric continuous. A jump in its first normal derivative makes the connection discontinuous and produces a delta term in curvature, but no delta squared. Einstein’s equation gives the Israel junction condition

8πGSab=−([Kab]−hab[K]),8\pi G S_{ab} =-\left( [\mathcal K_{ab}]-h_{ab}[\mathcal K] \right),

where SabS_{ab} is the surface stress tensor and [K]=hab[Kab][\mathcal K]=h^{ab}[\mathcal K_{ab}]. The labels and normal are linked: if the normal is reversed, interchange M+M^+ and M−M^- as well. Then [K]new=(−K−)−(−K+)=[K][\mathcal K]_{\mathrm{new}}=(-\mathcal K^-)-(-\mathcal K^+) =[\mathcal K], and the physical SabS_{ab} is unchanged. Reversing only the normal while retaining the jump labels would incorrectly reverse the source.

The formula has three hypotheses that are often hidden: the hypersurface is non-null, the first fundamental form matches, and the metric is sufficiently regular on each side to possess traces of K±\mathcal K^\pm. A discontinuous induced metric is not an ordinary thin shell.

On a null hypersurface, the normal is tangent and the induced metric is degenerate. There is no unit normal and no inverse three-metric with which to take the trace used above. One introduces a transverse null or non-null vector, a two-dimensional screen metric, and a transverse curvature. Its jump determines the surface stress and any impulsive gravitational part.

The split depends on the normalization of the generator and transverse vector, while the final distributional Einstein tensor is geometric. Barrabès and Israel construct this null-shell formalism in Barrabès–Israel 1991, §§II–IV. Substituting a null normal directly into the non-null Israel formula is therefore not a limiting derivation.

In Brinkmann coordinates, an idealized impulsive pp-wave can be written

ds2=−2 du dv+dx2+dy2+H(x,y)δ(u) du2.ds^2=-2\,du\,dv+dx^2+dy^2 +H(x,y)\delta(u)\,du^2.

The wavefront is u=0u=0. Formally, the nonzero distributional curvature components are

Ruiuj=−12 ∂i∂jH δ(u),i,j∈{x,y},R_{uiuj} =-\frac12\,\partial_i\partial_jH\,\delta(u), \qquad i,j\in\{x,y\},

and hence

Ruu=−12 Δ⊥H δ(u).R_{uu} =-\frac12\,\Delta_\perp H\,\delta(u).

If HH is harmonic in the transverse plane on the region considered, the Ricci distribution vanishes there while the trace-free Hessian carries an impulsive gravitational wave. If Δ⊥H\Delta_\perp H is nonzero, the wavefront also carries null matter.

The displayed metric coefficient itself contains a delta, so ordinary distribution theory cannot justify every intermediate inverse and product. One rigorous route regularizes the impulse and proves a distributional or generalized-function limit. Another uses a continuous coordinate form in which the metric is locally Lipschitz and geodesics can be matched across the wavefront. Podolský compares distributional and continuous forms and their geodesics in Podolský 2002, §§2–5.

For the displayed pp-wave, the regions u≠0u\ne0 are flat. Choose uu as an affine parameter for a geodesic crossing the wavefront, and replace δ\delta by a smooth unit-integral pulse before taking the limit. The transverse equations are

d2xidu2=12∂iH δ(u),[dxidu]=12∂iH(x0,y0).\frac{d^2x^i}{du^2}=\frac12\partial_iH\,\delta(u), \qquad \left[\frac{dx^i}{du}\right]=\frac12\partial_iH(x_0,y_0).

The transverse positions (x,y)(x,y) remain continuous. The Brinkmann coordinate vv, however, has the jump [v]=H(x0,y0)/2[v]=H(x_0,y_0)/2; it is not correct to say that every coordinate position remains continuous in this distributional chart. In a continuous chart the same geometry is represented by continuous geodesics with the appropriate regularity. The coordinate relation is described in Griffiths–Podolský 2009, §20.2, Eqs. (20.3)–(20.6), pp. 395–396. Their profile enters as −2Hδ(u)du2-2\mathcal H\delta(u)du^2, so our profile is H=−2HH=-2\mathcal H.

Multiplying the geodesic equation by distributions without regularization can create ambiguous products such as δ(u)u˙\delta(u)\dot u evaluated on a nonsmooth trajectory. A valid result states the coordinate form, regularization or generalized solution notion, and which limiting quantities are invariant.

The transverse-curvature jump in continuous coordinates

Section titled “The transverse-curvature jump in continuous coordinates”

For H=A(x2−y2)H=A(x^2-y^2), a continuous form of the metric near the impulse is

ds2=−2 du dV+(1+Au+)2dX2+(1−Au+)2dY2,u+=max⁡(u,0).ds^2=-2\,du\,dV +(1+A u_+)^2dX^2+(1-A u_+)^2dY^2, \qquad u_+=\max(u,0).

Work where ∣Au∣<1|Au|\lt1 to exclude the coordinate caustic. If the coordinates have length units, AA has units of inverse length. The null generator is k=∂Vk=\partial_V and a transverse vector is N=∂uN=\partial_u, with g(k,N)=−1g(k,N)=-1. With screen metric qABq_{AB} and the convention CAB=g(∇∂AN,∂B)=12∂uqABC_{AB}=g(\nabla_{\partial_A}N,\partial_B)=\tfrac12\partial_uq_{AB},

[CAB]=(A00−A),qAB[CAB]=0.[C_{AB}]=\begin{pmatrix}A&0\\0&-A\end{pmatrix}, \qquad q^{AB}[C_{AB}]=0.

The trace-free jump produces impulsive Weyl curvature and no null Ricci source. This agrees with Δ⊥H=0\Delta_\perp H=0 in the distributional chart. The corresponding kicks are [x′]=Ax0[x']=Ax_0 and [y′]=−Ay0[y']=-Ay_0.

A constant reparametrization u′=buu'=bu, V′=V/bV'=V/b, b>0b\gt0, changes k′=bkk'=bk, N′=N/bN'=N/b, [CAB′]=[CAB]/b[C'_{AB}]=[C_{AB}]/b, and the Brinkmann profile to H′=H/bH'=H/b. Nevertheless the singular curvature tensor is unchanged, since δ(u′)du′⊗du′=bδ(u)du⊗du\delta(u')du'\otimes du'=b\delta(u)du\otimes du. This explicit cancellation illustrates why a null-shell density or transverse-curvature component alone depends on normalization, while the complete tensor distribution does not.

Solved diagnostic: vacuum impulse or null source?

Section titled “Solved diagnostic: vacuum impulse or null source?”

For each of the two transverse profiles below, compute Δ⊥H\Delta_\perp H and decide whether the impulse is Ricci-flat or carries a null surface source.

Solution

Take

H(x,y)=A(x2−y2).H(x,y)=A(x^2-y^2).

Then Δ⊥H=2A−2A=0\Delta_\perp H=2A-2A=0, so Ruu=0R_{uu}=0. Yet

∂x2H=2A,∂y2H=−2A,\partial_x^2H=2A, \qquad \partial_y^2H=-2A,

and the trace-free curvature components are nonzero on u=0u=0. This is a vacuum gravitational impulse.

For H=A(x2+y2)H=A(x^2+y^2), Δ⊥H=4A\Delta_\perp H=4A, so

Ruu=−2A δ(u).R_{uu}=-2A\,\delta(u).

Here the scalar curvature vanishes, so Einstein’s equation requires Tuu=−Aδ(u)/(4πG)T_{uu}=-A\delta(u)/(4\pi G). Nonnegative null energy therefore requires A≤0A\le0 for this particular sign convention and profile. The Laplacian, not the mere presence of a delta-profile metric term, distinguishes the Ricci source from the free gravitational impulse.

The final page asks what happens when singular or rapidly oscillating solutions are taken in sequences rather than prescribed on one fixed hypersurface: Compactness and Limits of Spacetimes.

  • Barrabès, Claude, and Werner Israel. “Thin Shells in General Relativity and Cosmology: The Lightlike Limit.” Physical Review D 43 (1991): 1129–1142. DOI.
  • Griffiths, Jerry B., and Jiří Podolský. Exact Space-Times in Einstein’s General Relativity. Cambridge Monographs on Mathematical Physics. Cambridge University Press, 2009. DOI.
  • Podolský, Jiří. “Exact Impulsive Gravitational Waves in Spacetimes of Constant Curvature.” In Gravitation: Following the Prague Inspiration, 205–246. World Scientific, 2002. arXiv:gr-qc/0201029.
  • Geroch, Robert, and Jennie Traschen. “Strings and Other Distributional Sources in General Relativity.” Physical Review D 36 (1987): 1017–1031. DOI.