Frame fields in general relativity

In general relativity, a frame field (also called a tetrad or vierbein) is a set of four orthonormal vector fields, one timelike and three spacelike, defined on a Lorentzian manifold that is physically interpreted as a model of spacetime. The timelike unit vector field is often denoted by $\vec{e}_0$ and the three spacelike unit vector fields by $\vec{e}_1, \vec{e}_2, \, \vec{e}_3$. All tensorial quantities defined on the manifold can be expressed using the frame field and its dual coframe field.

Frames were introduced into general relativity by Hermann Weyl in 1929^[1].

The general theory of tetrads (and analogs in dimensions other than 4) is described in the article on Cartan formalism; the index notation for tetrads is explained in tetrad (index notation).

Physical interpretation

Frame fields always correspond to a family of ideal observers immersed in the given spacetime; the integral curves of the timelike unit vector field are the worldlines of these observers, and at each event along a given worldline, the three spacelike unit vector fields specify the spatial triad carried by the observer. The triad may be thought of as defining the spatial coordinate axes of a local laboratory frame, which is valid very near the observer's worldline.

In general, the worldlines of these observers need not be timelike geodesics. If any of the worldlines bends away from a geodesic path in some region, we can think of the observers as test particles that accelerate by using ideal rocket engines with a thrust equal to the magnitude of their acceleration vector. Alternatively, if our observer is attached to a bit of matter in a ball of fluid in hydrostatic equilibrium, this bit of matter will in general be accelerated outward by the net effect of pressure holding up the fluid ball against the attraction of its own gravity. Other possibilities include an observer attached to a free charged test particle in an electrovacuum solution, which will of course be accelerated by the Lorentz force, or an observer attached to a spinning test particle, which may be accelerated by a spin-spin force.

It is important to recognize that frames are geometric objects. That is, vector fields make sense (in a smooth manifold) independently of choice of a coordinate chart, and (in a Lorentzian manifold), so do the notions of orthogonality and length. Thus, just like vector fields and other geometric quantities, frame fields can be represented in various coordinate charts. But computations of the components of tensorial quantities, with respect to a given frame, will always yield the same result, whichever coordinate chart is used to represent the frame.

These fields are required to write the Dirac equation in curved spacetime.

Specifying a frame

To write down a frame, a coordinate chart on the Lorentzian manifold needs to be chosen. Then, every vector field on the manifold can be written down as a linear combination of the four coordinate basis vector fields:

$\vec{X} = X^j \, \partial_{x^j}$

(Here, the Einstein summation convention is used, and the vector fields are thought of as first order linear differential operators, and the components X^j are often called contravariant components.) In particular, the vector fields in the frame can be expressed this way:

$\vec{e}_a = {e_a}^j \, \partial_{x^j}$

In "designing" a frame, one naturally needs to ensure, using the given metric, that the four vector fields are everywhere orthonormal.

Once a signature is adopted (in the case of a four-dimensional Lorentzian manifold, the signature is -1 + 3), by duality every vector has a dual covector and conversely. Thus, every frame field is associated with a unique coframe field, and vice versa.

Specifying the metric using a coframe

Alternatively, the metric tensor can be specified by writing down a coframe in terms of a coordinate basis and stipulating that the metric tensor is given by

$g = -\sigma^0 \otimes \sigma^0 + \sigma^1 \otimes \sigma^1 + \sigma^2 \otimes \sigma^2 + \sigma^3 \otimes \sigma^3$

This is just a fancy way of saying that the coframe is orthonormal. Whether this is used to obtain the metric tensor after writing down the frame (and passing to the dual coframe), or starting with the metric tensor and using it to verify that a frame has been obtained by other means, it must always hold true.

Relationship with metric tensor, in a coordinate basis

The vierbein field, $e^{\mu}_{\ a}$, has two indices: $\mu \,$ labels the general spacetime coordinate and $a \,$ labels the local lorentz spacetime or local laboratory coordinates.

The vierbein field or frame fields can be regarded as the square root of the metric tensor, $g^{\mu \nu} \,$, since in a coordinate basis,

$g^{\mu \nu}= e^{\mu}_{\ a} e^{\nu}_{\ b} \eta^{ab} \,$

where $\eta^{ab} \,$ is the Lorentz metric.

Local lorentz indices are raised and lowered with the lorentz metric in the same way as general spacetime coordinates are raised and lowered with the metric tensor. For example:

$T^a = \eta^{ab} T_b \,$

The vierbein fields enables conversion between spacetime and local lorentz indices. For example:

$T_a = e^\mu_a T_\mu \,$

The vierbein field itself can be manipulated in the same fashion:

$e^\nu_a = e^\mu_a e^\nu_\mu \,$, since $e^\nu_\mu = \delta^\nu_\mu \,$

And these can combine.

$T^a = e_\mu^a T^\mu \,$

A few more examples: Spacetime and local lorentz coordinates can be mixed together:

$T^{\mu a}=e_\nu^a T^{\mu \nu}$

The local lorentz coordinates transform differently from the general spacetime coordinates. Under a general coordinate transformation we have:

$T'^{\mu a} = \frac{\partial x'^\mu}{\partial x^\nu}T^{\nu a}$

whilst under a local lorentz transformation we have:

$T'^{\mu a} = \Lambda(x)^a_{\ b} T^{\mu b}$.

Comparison with coordinate basis

Coordinate basis vectors have the special property that their Lie brackets pairwise vanish. Except in locally flat regions, at least some Lie brackets of vector fields from a frame will not vanish. The resulting baggage needed to compute with them is acceptable, as components of tensorial objects with respect to a frame (but not with respect to a coordinate basis) have a direct interpretation in terms of measurements made by the family of ideal observers corresponding the frame.

Coordinate basis vectors can very well be null, which, by definition, cannot happen for frame vectors.

Nonspinning and inertial frames

Some frames are nicer than others. Particularly in vacuum or electrovacuum solutions, the physical experience of inertial observers (who feel no forces) may be of particular interest. The mathematical characterization of an inertial frame is very simple: the integral curves of the timelike unit vector field must define a geodesic congruence, or in other words, its acceleration vector must vanish:

$\nabla_{\vec{e}_{0}} \, \vec{e}_0 = 0$

It is also often desirable to ensure that the spatial triad carried by each observer does not rotate. In this case, the triad can be viewed as being gyrostabilized. The criterion for a nonspinning inertial (NSI) frame is again very simple:

$\nabla_{\vec{e}_0} \, \vec{e}_j = 0, \; \; j = 0 \dots 3$

This says that as we move along the worldline of each observer, their spatial triad is parallel-transported. Nonspinning inertial frames hold a special place in general relativity, because they are as close as we can get in a curved Lorentzian manifold to the Lorentz frames used in special relativity (these are special nonspinning inertial frames in the Minkowski vacuum).

More generally, if the acceleration of our observers is nonzero, $\nabla_{\vec{e}_0}\,\vec{e}_0 \neq 0$, we can replace the covariant derivatives

$\nabla_{\vec{e}_0} \, \vec{e}_j, \; j = 1 \dots 3$

with the (spatially projected) Fermi-Walker derivatives to define a nonspinning frame.

Given a Lorentzian manifold, we can find infinitely many frame fields, even if we require additional properties such as inertial motion. However, a given frame field might very well be defined on only part of the manifold.

Example: Static observers in Schwarzschild vacuum

It will be instructive to consider in some detail a few simple examples. Consider the famous Schwarzschild vacuum that models spacetime outside an isolated nonspinning spherically symmetric massive object, such as a star. In most textbooks one finds the metric tensor written in terms of a static polar spherical chart, as follows:

$ds^2 = -(1-2m/r) \, dt^2 + \frac{dr^2}{1-2m/r} + r^2 \, \left( d\theta^2 + \sin(\theta)^2 \, d\phi^2 \right)$

$-\infty < t < \infty, \; 2 m < r < \infty, \; 0 < \theta < \pi, \; -\pi < \phi < \pi$

More formally, the metric tensor can be expanded with respect to the coordinate cobasis as

$g = -(1-2m/r) \, dt \otimes dt + \frac{1}{1-2m/r} \, dr \otimes dr + r^2 \, d\theta \otimes d\theta + r^2 \sin(\theta)^2 \, d\phi \otimes d\phi$

A coframe can be read off from this expression:

$\sigma^0 = -\sqrt{1-2m/r} \, dt, \; \sigma^1 = \frac{dr}{\sqrt{1-2m/r}}, \; \sigma^2 = r d\theta, \; \sigma^3 = r \sin(\theta) d\phi$

To see that this coframe really does correspond to the Schwarzschild metric tensor, just plug this coframe into

$g = -\sigma^0 \otimes \sigma^0 + \sigma^1 \otimes \sigma^1 + \sigma^2 \otimes \sigma^2 + \sigma^3 \otimes \sigma^3$

The frame dual to the coframe is

$\vec{e}_0 = \frac{1}{\sqrt{1-2m/r}} \partial_t, \; \vec{e}_1 = \sqrt{1-2m/r} \partial_r, \; \vec{e}_2 = \frac{1}{r} \partial_\theta, \; \vec{e}_3 = \frac{1}{r \sin(\theta)} \partial_\phi$

(The minus sign on σ⁰ ensures that $\vec{e}_0$ is future pointing.) This is the frame that models the experience of static observers who use rocket engines to "hover" over the massive object. The thrust they require to maintain their position is given by the magnitude of the acceleration vector

$\nabla_{\vec{e}_0} \vec{e}_0 = \frac{m/r^2}{\sqrt{1-2m/r}} \, \vec{e}_1$

This is radially outward pointing, since the observers need to accelerate away from the object to avoid falling toward it. On the other hand, the spatially projected Fermi derivatives of the spatial basis vectors (with respect to $\vec{e}_0$) vanish, so this is a nonspinning frame.

The components of various tensorial quantities with respect to our frame and its dual coframe can now be computed.

For example, the tidal tensor for our static observers is defined using tensor notation (for a coordinate basis) as

$E[X]_{ab} = R_{ambn} \, X^m \, X^n$

where we write $\vec{X} = \vec{e}_0$ to avoid cluttering the notation. Its only non-zero components with respect to our coframe turn out to be

$E[X]_{11} = -2m/r^3, \; E[X]_{22} = E[X]_{33} = m/r^3$

The corresponding coordinate basis components are

$E[X]_{rr} = -2m/r^3/(1-2m/r), \; E[X]_{\theta \theta} = m/r, \; E[X]_{\phi \phi} = m \sin(\theta)^2/r$

(A quick note concerning notation: many authors put carets over abstract indices referring to a frame. When writing down specific components, it is convenient to denote frame components by 0,1,2,3 and coordinate components by t,r,θ,φ. Since an expression like S_ab = 36m / r doesn't make sense as a tensor equation, there should be no possibility of confusion.)

Compare the tidal tensor Φ of Newtonian gravity, which is the traceless part of the Hessian of the gravitational potential U. Using tensor notation for a tensor field defined on three-dimensional euclidean space, this can be written

$\Phi_{ij} = U_{,i j} - \frac{1}{3} {U^{,k}}_{,k} \, \eta_{ij}$

The reader may wish to crank this through (notice that the trace term actually vanishes identically when U is harmonic) and compare results with the following elementary approach: we can compare the gravitational forces on two nearby observers lying on the same radial line:

$m/(r+h)^2 - m/r^2 = -2m/r^3 \, h + 3m/r^4 \, h^2 + O(h^3)$

Because in discussing tensors we are dealing with multilinear algebra, we retain only first order terms, so Φ₁₁ = − 2m / r³. Similarly, we can compare the gravitational force on two nearby observers lying on the same sphere r = r₀. Using some elementary trigonometry and the small angle approximation, we find that the force vectors differ by a vector tangent to the sphere which has magnitude

$\frac{m}{r_0^2} \, \sin(\theta) \approx \frac{m}{r_0^2} \, \frac{h}{r_0} = \frac{m}{r_0^3} \, h$

By using the small angle approximation, we have ignored all terms of order O(h²), so the tangential components are Φ₂₂ = Φ₃₃ = m / r³. Here, we are referring to the obvious frame obtained from the polar spherical chart for our three-dimensional euclidean space:

$\vec{\epsilon}_1 = \partial_r, \; \vec{\epsilon}_2 = \frac{1}{r} \, \partial_\theta, \; \vec{\epsilon}_3 = \frac{1}{r \sin \theta} \, \partial_\phi$

Plainly, the coordinate components $E[X]_{\theta \theta}, \, E[X]_{\phi \phi}$ computed above don't even scale the right way, so they clearly cannot correspond to what an observer will measure even approximately. (By coincidence, the Newtonian tidal tensor components agree exactly with the relativistic tidal tensor components we wrote out above.)

Example: Lemaître observers in the Schwarzschild vacuum

To find an inertial frame, we can boost our static frame in the $\vec{e}_1$ direction by an undetermined boost parameter (depending on the radial coordinate), compute the acceleration vector of the new undetermined frame, set this equal to zero, and solve for the unknown boost parameter. The result will be a frame which we can use to study the physical experience of observers who fall freely and radially toward the massive object. By appropriately choosing an integration constant, we obtain the frame of Lemaître observers, who fall in from rest at spatial infinity. (This phrase doesn't make sense, but the reader will no doubt have no difficulty in understanding our meaning.) In the static polar spherical chart, this frame can be written

$\vec{f}_0 = \frac{1}{1-2m/r} \, \partial_t - \sqrt{2m/r} \, \partial_r$

$\vec{f}_1 = \partial_r - \frac{\sqrt{2m/r}}{1-2m/r} \, \partial_t$

$\vec{f}_3 = \frac{1}{r} \, \partial_\theta$

$\vec{f}_3 = \frac{1}{r \sin(\theta)} \, \partial_\phi$

Note that $\vec{e}_0 \neq \vec{f}_0, \; \vec{e}_1 \neq \vec{f}_1$, and that $\vec{e}_0$ "leans inwards", as it should, since its integral curves are timelike geodesics representing the world lines of infalling observers. Indeed, since the covariant derivatives of all four basis vectors (taken with respect to $\vec{e}_0$) vanish identically, our new frame is a nonspinning inertial frame.

If our massive object is in fact a (nonrotating) black hole, we probably wish to follow the experience of the Lemaître observers as they fall through the event horizon at r = 2m. Since the static polar spherical coordinates have a coordinate singularity at the horizon, we'll need to switch to a more appropriate coordinate chart. The simplest possible choice is to define a new time coordinate by

$T(t,r) = t - \int \frac{\sqrt{2m/r}}{1-2m/r} \, dr = t + 2 \sqrt{2mr} + 2m \log \left( \frac{\sqrt{r}-\sqrt{2m}}{\sqrt{r}+\sqrt{2m}} \right)$

This gives the Painlevé chart. The new line element is

$ds^2 = -dT^2 + \left( dr + \sqrt{2m/r} \, dT \right)^2 + r^2 \left( d\theta^2 + \sin(\theta)^2 \, d\phi^2 \right)$

$-\infty < T < \infty, \; 0 < r < \infty, \; 0 < \theta < \pi, \; -\pi < \phi < \pi$

With respect to the Painlevé chart, the Lemaître frame is

$\vec{f}_0 = \partial_T - \sqrt{2m/r} \, \partial_r$

$\vec{f}_1 = \partial_r$

$\vec{f}_2 = \frac{1}{r} \, \partial_\theta$

$\vec{f}_3 = \frac{1}{r \sin(\theta)} \, \partial_\phi$

Notice that their spatial triad looks exactly like the frame for three-dimensional euclidean space which we mentioned above (when we computed the Newtonian tidal tensor). Indeed, the spatial hyperslices T = T₀ turn out to be locally isometric to flat three-dimensional euclidean space! (This is a remarkable and rather special property of the Schwarzschild vacuum; most spacetimes do not admit a slicing into flat spatial sections.)

The tidal tensor taken with respect to the Lemaître observers is

$E[Y]_{ab} = R_{ambn} \, Y^m \, Y^n$

where we write $Y = \vec{f}_0$ to avoid cluttering the notation. This is a different tensor from the one we obtained above, because it is defined using a different family of observers. Nonetheless, its nonvanishing components look familiar: $E[Y]_{11} = -2m/r^3, \, E[Y]_{22} = E[Y]_{33} = m/r^3$. (This is again a rather special property of the Schwarzschild vacuum.)

Notice that there is simply no way of defining static observers on or inside the event horizon. On the other hand, the Lemaître observers are not defined on the entire exterior region covered by the static polar spherical chart either, so in these examples, neither the Lemaître frame nor the static frame are defined on the entire manifold.

Example: Hagihara observers in the Schwarzschild vacuum

In the same way that we found the Lemaître observers, we can boost our static frame in the $\vec{e}_3$ direction by an undetermined parameter (depending on the radial coordinate), compute the acceleration vector, and require that this vanish in the equatorial plane θ = π / 2. The new Hagihara frame describes the physical experience of observers in stable circular orbits around our massive object. It was apparently first discussed by the distinguished (and mathematically gifted) astronomer Yusuke Hagihara.

In the static polar spherical chart, the Hagihara frame is

$\vec{h}_0 = \frac{1}{\sqrt{1-3m/r}} \, \partial_t + \frac{\sqrt{m/r^3}}{\sqrt{1-3m/r} \, \sin(\theta)} \, \partial_\phi$

$\vec{h}_1 = \sqrt{1-2m/r} \, \partial_r$

$\vec{h}_2 = \frac{1}{r} \, \partial_\theta$

$\vec{h}_3 = \frac{\sqrt{1-2m/r}}{\sqrt{1-3m/r} \,\sin(\theta)} \, \partial_\phi - \frac{\sqrt{m/r^3}}{\sqrt{1-2m/r} \, \sqrt{1-3m/r}} \, \partial_t$

which in the equatorial plane becomes

$\vec{h}_0 = \frac{1}{\sqrt{1-3m/r}} \, \partial_t + \frac{\sqrt{m/r^3}}{\sqrt{1-3m/r}} \, \partial_\phi$

$\vec{h}_1 = \sqrt{1-2m/r} \, \partial_r$

$\vec{h}_2 = \frac{1}{r} \, \partial_\theta$

$\vec{h}_3 = \frac{\sqrt{1-2m/r}}{\sqrt{1-3m/r}} \, \partial_\phi - \frac{\sqrt{m/r^3}}{\sqrt{1-2m/r} \, \sqrt{1-3m/r}} \, \partial_t$

The tidal tensor E[Z]_ab where $\vec{Z} = \vec{h}_0$ turns out to be given (in the equatorial plane) by

$E[Z]_{11} = -\frac{m}{r^3} \, \frac{2-3m/r}{1-2m/r} = -\frac{2m}{r^3} - \frac{m^2}{r^4} + O(1/r^5)$

$E[Z]_{22} = \frac{m}{r^3} \, \frac{1}{1-3m/r} = -\frac{m}{r^3} + \frac{3m^2}{r^4} + O(1/r^5)$

$E[Z]_{33} = \frac{m}{r^3}$

Thus, compared to a static observer hovering at a given coordinate radius, a Hagihara observer in a stable circular orbit with the same coordinate radius will measure radial tidal forces which are slightly larger in magnitude, and transverse tidal forces which are no longer isotropic (but slightly larger orthogonal to the direction of motion).

Note that the Hagihara frame is only defined on the region r > 3m. Indeed, stable circular orbits only exist on r > 6m, so the frame should not be used inside this locus.

Computing Fermi derivatives shows that the frame field just given is in fact spinning with respect to a gyrostabilized frame. The principal reason why is easy to spot: in this frame, each Hagihara observer keeps his spatial vectors radially aligned, so $\vec{h}_1, \; \vec{h}_3$ rotate about $\vec{h}_2$ as the observer orbits around the central massive object. However, after correcting for this observation, a small precession of the spin axis of a gyroscope carried by a Hagihara observer still remains; this is the de Sitter precession effect (also called the geodetic precession effect).

Generalizations

This article has focused on the application of frames to general relativity, and particularly on their physical interpretation. Here we very briefly outline the general concept. In an n-dimensional Riemannian manifold or pseudo-Riemannian manifold, a frame field is a set of orthonormal vector fields which forms a basis for the tangent space at each point in the manifold. As before, frames can be specified in terms of a given coordinate basis, and in a non-flat region, some of their pairwise Lie brackets will fail to vanish.

In fact, given any inner-product space V, we can define a new space consisting of all tuples of orthonormal bases for V. Applying this construction to each tangent space yields the orthonormal frame bundle of a (pseudo-)Riemannian manifold and a frame field is a section of this bundle. More generally still, we can consider frame bundles associated to any vector bundle, or even arbitrary principal fiber bundles. The notation becomes a bit more involved because it is harder to avoid distinguishing between indices referring to the base, and indices referring to the fiber. Many authors speak of internal components when referring to components indexed by the fiber.

See also

• Exact solutions in general relativity
• Georges Lemaître
• Karl Schwarzschild
• Method of moving frames
• Paul Painlevé
• Vierbein
• Yusuke Hagihara

References

1. ^ Hermann Weyl "Elektron und Gravitation I", Zeitshrift Physik, 56, p330-352, 1929.

• Flanders, Harley (1989). Differential Forms with Applications to the Physical Sciences. New York: Dover. ISBN 0-486-66169-5.  See Chapter IV for frames in E³, then see Chapter VIII for frame fields in Riemannian manifolds. This book doesn't really cover Lorentzian manifolds, but with this background in hand the reader is well prepared for the next citation.
• Misner, Charles; Thorne, Kip S. & Wheeler, John Archibald (1973). Gravitation. San Francisco: W. H. Freeman. ISBN 0-7167-0344-0.  In this book, a frame field (coframe field) is called an anholonomic basis of vectors (covectors). Essential information is widely scattered about, but can be easily found using the extensive index.
• Landau, L. D. & Lifschitz, E. F. (1980). The Classical Theory of Fields (4th ed.). London: Butterworth-Heinemann. ISBN 0-7506-2768-9.  In this book, a frame field is called a tetrad (not to be confused with the now standard term NP tetrad used in the Newman-Penrose formalism). See Section 98.
• De Felice, F.; and Clarke, C. J. (1992). Relativity on Curved Manifolds. Cambridge: Cambridge University Press. ISBN 0-521-42908-0.  See Chapter 4 for frames and coframes. If you ever need more information about frame fields, this might be a good place to look!