In general relativity, a congruence (more properly, a congruence of curves) is the set of integral curves of a (nowhere vanishing) vector field in a four-dimensional Lorentzian manifold which is interpreted physically as a model of spacetime. Often this manifold will be taken to be an exact or approximate solution to the Einstein field equation.
Congruences generated by nowhere vanishing timelike, null, or spacelike vector fields are called timelike, null, or spacelike respectively.
\vec{X}
\nabla\vec{X
The integral curves of the vector field are a family of non-intersecting parameterized curves which fill up the spacetime. The congruence consists of the curves themselves, without reference to a particular parameterization. Many distinct vector fields can give rise to the same congruence of curves, since if
f
\vec{X}
\vec{Y}=f\vec{X}
However, in a Lorentzian manifold, we have a metric tensor, which picks out a preferred vector field among the vector fields which are everywhere parallel to a given timelike or spacelike vector field, namely the field of tangent vectors to the curves. These are respectively timelike or spacelike unit vector fields.
In general relativity, a timelike congruence in a four-dimensional Lorentzian manifold can be interpreted as a family of world lines of certain ideal observers in our spacetime. In particular, a timelike geodesic congruence can be interpreted as a family of free-falling test particles.
Null congruences are also important, particularly null geodesic congruences, which can be interpreted as a family of freely propagating light rays.
Warning: the world line of a pulse of light moving in a fiber optic cable would not in general be a null geodesic, and light in the very early universe (the radiation-dominated epoch) was not freely propagating. The world line of a radar pulse sent from Earth past the Sun to Venus would however be modeled as a null geodesic arc. In dimensions other than four, the relationship between null geodesics and "light" no longer holds: If "light" is defined as the solution to the Laplacian wave equation, then the propagator has both null and time-like components in odd space-time dimensions and is no longer a pure Dirac delta function in even space-time dimensions greater than four.
Describing the mutual motion of the test particles in a null geodesic congruence in a spacetime such as the Schwarzschild vacuum or FRW dust is a very important problem in general relativity. It is solved by defining certain kinematical quantities which completely describe how the integral curves in a congruence may converge (diverge) or twist about one another.
It should be stressed that the kinematical decomposition we are about to describe is pure mathematics valid for any Lorentzian manifold. However, the physical interpretation in terms of test particles and tidal accelerations (for timelike geodesic congruences) or pencils of light rays (for null geodesic congruences) is valid only for general relativity (similar interpretations may be valid in closely related theories).
Consider the timelike congruence generated by some timelike unit vector field X, which we should think of as a first order linear partial differential operator. Then the components of our vector field are now scalar functions given in tensor notation by writing
\vec{X}f=f,aXa
\nabla\vec{X
X |
a=
a} | |
{X | |
;b |
Xb
\left(
X |
aXb+
a} | |
{X | |
;b |
\right)Xb=
a} | |
{X | |
;b |
Xb-
X |
a=0
a} | |
{X | |
;b |
hab=gab+XaXb
\vec{X}
X |
aXb+Xa;b=
m} | |
{h | |
a |
n} | |
{h | |
b |
Xm;
X |
aXb+Xa;b=\thetaab+\omegaab
\thetaab=
m} | |
{h | |
a |
n} | |
{h | |
b |
X(m;n)
\omegaab=
m} | |
{h | |
a |
n} | |
{h | |
b |
X[m;n]
Because these tensors live in the spatial hyperplane elements orthogonal to
\vec{X}
\theta
\thetaab=\sigmaab+
1 | |
3 |
\thetahab
Xa;b=\sigmaab+\omegaab+
1 | |
3 |
\thetahab-
X |
aXb
The expansion scalar, shear tensor (
\sigmaab
See the citations and links below for justification of these claims.
By the Ricci identity (which is often used as the definition of the Riemann tensor), we can write:
Xa;bn-Xa;nb=RambnXm
In the famous slogan of John Archibald Wheeler:
Spacetime tells matter how to move; matter tells spacetime how to curve.We now see how to precisely quantify the first part of this assertion; the Einstein field equation quantifies the second part.
In particular, according to the Bel decomposition of the Riemann tensor, taken with respect to our timelike unit vector field, the electrogravitic tensor (or tidal tensor) is defined by:
E[\vec{X}]ab=RambnXmXn
\left(Xa:bn-Xa:nb\right)Xn=E[\vec{X}]ab
\begin{align} E[\vec{X}]ab&=
2 | |
3 |
\theta\sigmaab-\sigmaam
m} | |
{\sigma | |
b |
-\omegaam
m} | |
{\omega | |
b |
\\ &-
1 | |
3 |
\left(
\theta |
+
\theta2 | |
3 |
\right)hab-
m} | |
{h | |
a |
n} | |
{h | |
b |
\left(
\sigma |
mn-
X |
(m;n)\right)-
X |
a
X |
b\\ \end{align}
In this section, we turn to the problem of obtaining evolution equations (also called propagation equations or propagation formulae).
It will be convenient to write the acceleration vector as
X |
a=Wa
Jab=Xa:b=
\theta | |
3 |
hab+\sigmaab+\omegaab-
X |
aXb
J |
ab=Jan;bXn-E[\vec{X}]ab
\left(JanXn\right);b=Jan;bXn+Jan
n} | |
{X | |
;b |
=Jan;bXn+Jam
m} | |
{J | |
b |
J |
ab=-Jam
m} | |
{J | |
b |
-{E[\vec{X}]}ab+Wa;b
Jab
Consider first the easier case when the acceleration vector vanishes. Then (observing that the projection tensor can be used to lower indices of purely spatial quantities), we have:
Jam
m} | |
{J | |
b |
=
\theta2 | |
9 |
hab+
2\theta | |
3 |
\left(\sigmaab+\omegaab\right)+\left(\sigmaam
m} | |
{\sigma | |
b |
+\omegaam
m} | |
{\omega | |
b |
\right)+\left(\sigmaam
m} | |
{\omega | |
b |
+\omegaam
m} | |
{\sigma | |
b |
\right)
J |
ab=-
\theta2 | |
9 |
hab-
2\theta | |
3 |
\left(\sigmaab+\omegaab\right)-\left(\sigmaam
m} | |
{\sigma | |
b |
+\omegaam
m} | |
{\omega | |
b |
\right)-\left(\sigmaam
m} | |
{\omega | |
b |
+\omegaam
m} | |
{\sigma | |
b |
\right)-{E[\vec{X}]}ab
\Sigma,\Omega
\Sigma2+\Omega2
\Sigma\Omega+\Omega\Sigma
\theta |
=\omega2-\sigma2-
\theta2 | |
3 |
-
m} | |
{E[\vec{X}] | |
m |
\sigma |
ab=-
2\theta | |
3 |
\sigmaab-\left(\sigmaam
m} | |
{\sigma | |
b |
+\omegaam
m} | |
{\omega | |
b |
\right)-{E[\vec{X}]}ab+
\sigma2-\omega2+{E[\vec{X | |
] |
m} | |
m}{3} |
hab
\omega |
ab=-
2\theta | |
3 |
\omegaab-\left(\sigmaam
m} | |
{\omega | |
b |
+\omegaam
m} | |
{\sigma | |
b |
\right)
\sigma2=\sigmamn\sigmamn, \omega2=\omegamn\omegamn
\sigma,\omega
a} | |
{E[\vec{X}] | |
a |
=RmnXmXn