CongruenceProperties - Maple Help

Tensor[CongruenceProperties] - calculate properties of a congruence of curves

Calling Sequences

CongruenceProperties(${\mathbit{g}}$, U)

CongruenceProperties(${\mathbit{g}}$, K, L)

CongruenceProperties(${\mathbit{g}}$, K)

CongruenceProperties(${\mathbit{g}}$, NT)

Parameters

g     - a metric tensor

U     - a unit vector

K,L   - normalized null vectors, the vector defines an affinely parameterized, geodesic null congruence.

NT    - a list of 4 vectors, defining a null tetrad, the first vector in the tetrad defines the geodesic null congruence.

Description

 • The command CongruenceProperties returns a table of properties associated to a line congruence defined by a unit (time-like or space-like) vector field $U$ or a null vector field $K$.
 • Let , set . The following scalar and tensor fields are calculated by the first calling sequence.

- Acceleration:

- Expansion: Θ = ${\nabla }_{a}{U}^{a}$ .

- Rotation Tensor : 1/2 (

- Shear Tensor: 1/2 (

 • The left-hand side of the Raychaudhuri equation ${U}^{a}{\nabla }_{a}$valid when the congruence is geodesic (${A}_{a}=$0), where is the Ricci tensor and is also calculated.
 • The first calling sequence returns a table with indices "Acceleration", "Expansion", "RotationTensor", "ShearTensor", "Raychaudhuri".
 • The remaining three calling sequences apply only to an affinely parameterized, geodesic null congruence , that is,  and
 • The second calling sequence requires where Setand

- Expansion: Θ = ${\nabla }_{a}{K}^{a}$ .

- Rotation Tensor:

- Rotation Scalar:

- Complex expansion: .

- Shear Tensor:

The Raychaudhuri equation is as above but using these definitions ofand and with

 • The second calling sequence returns a table with 8 indices "Expansion", "RotationNormSquared" "ShearNormSquared", "RotationTensor", "RotationScalar", "ShearTensor" , "ComplexExpansion" and "Raychaudhuri".
 • The third calling sequence calculates: Expansion: Θ = Rotation norm squared = and Shear norm squared = The definitions are as in the second calling sequence but, as these scalars do not in fact depend upon the choice of L, only the vector K is needed as input. The third calling sequence returns a table with indices "Expansion", "RotationNormSquared", "ShearNormSquared" and "Raychaudhuri".
 • Finally, from the 4th calling sequence we set and and calculate, in addition to the 8 quantities calculated for the second calling sequence , Newman-Penrose Spin Coefficients.

Examples

 > $\mathrm{with}\left(\mathrm{DifferentialGeometry}\right):$$\mathrm{with}\left(\mathrm{Tensor}\right):$

Example 1.

For our first example we use the standard metric on the sphere.

 > $\mathrm{DGsetup}\left(\left[\mathrm{θ},\mathrm{φ}\right],M\right)$
 ${\mathrm{frame name: M}}$ (2.1)
 > $g≔\mathrm{evalDG}\left({R}^{2}\left(\mathrm{dtheta}&t\mathrm{dtheta}+{\mathrm{sin}\left(\mathrm{θ}\right)}^{2}\mathrm{dphi}&t\mathrm{dphi}\right)\right)$
 ${\mathrm{_DG}}{}\left(\left[\left[{"tensor"}{,}{M}{,}\left[\left[{"cov_bas"}{,}{"cov_bas"}\right]{,}\left[{}\right]\right]\right]{,}\left[\left[\left[{1}{,}{1}\right]{,}{{R}}^{{2}}\right]{,}\left[\left[{2}{,}{2}\right]{,}{{R}}^{{2}}{}{{\mathrm{sin}}{}\left({\mathrm{θ}}\right)}^{{2}}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"tensor"}{,}{M}{,}\left[\left[{"cov_bas"}{,}{"cov_bas"}\right]{,}\left[{}\right]\right]\right]{,}\left[\left[\left[{1}{,}{1}\right]{,}{{R}}^{{2}}\right]{,}\left[\left[{2}{,}{2}\right]{,}{{R}}^{{2}}{}{{\mathrm{sin}}{}\left({\mathrm{θ}}\right)}^{{2}}\right]\right]\right]\right)$ (2.2)

Define a unit vector field $U$.

 M > $U≔\mathrm{evalDG}\left(\frac{1\mathrm{D_phi}}{R\mathrm{sin}\left(\mathrm{θ}\right)}\right)$
 ${\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}\frac{{1}}{{R}{}{\mathrm{sin}}{}\left({\mathrm{θ}}\right)}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}\frac{{1}}{{R}{}{\mathrm{sin}}{}\left({\mathrm{θ}}\right)}\right]\right]\right]\right)$ (2.3)

We see that the congruence is geodesic on the equator () but is accelerating elsewhere. It is shearing, rotating and non-expanding.

 M > $\mathrm{CongruenceProperties}\left(g,U\right)$
 ${\mathrm{table\left( \left[\left( "Raychaudhuri" \right) = /\left(1, *\left(^\left(R, 2\right)\right)\right), \left( "Acceleration" \right) = _DG\left(\left[\left["tensor", M, \left[\left["con_bas"\right], \left[\right]\right]\right], \left[\left[\left[1\right], +\left(-\left(/\left(*\left(cos\left(theta\right)\right), *\left(^\left(R, 2\right), *\left(sin\left(theta\right)\right)\right)\right)\right)\right)\right]\right]\right]\right), \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right), \left( "Expansion" \right) = 0, \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}{,}{\mathrm{table\left( \left[\left( "Raychaudhuri" \right) = /\left(1, *\left(^\left(R, 2\right)\right)\right), \left( "Acceleration" \right) = _DG\left(\left[\left["tensor", M, \left[\left["con_bas"\right], \left[\right]\right]\right], \left[\left[\left[1\right], +\left(-\left(/\left(*\left(cos\left(theta\right)\right), *\left(^\left(R, 2\right), *\left(sin\left(theta\right)\right)\right)\right)\right)\right)\right]\right]\right]\right), \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right), \left( "Expansion" \right) = 0, \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}$ (2.4)

Example 2.

For the next example we consider a class of Robinson-Trautman metrics. These are of Petrov type II and admit a null congruence which is shear-free.

 M > $\mathrm{DGsetup}\left(\left[u,r,\mathrm{zeta},\mathrm{zetab}\right],\mathrm{RT}\right)$
 ${\mathrm{frame name: RT}}$ (2.5)
 RT > $g≔\mathrm{evalDG}\left(2{r}^{2}{P\left(\mathrm{zeta},\mathrm{zetab},u\right)}^{-2}\mathrm{dzeta}&s\mathrm{dzetab}-2\mathrm{du}&s\mathrm{dr}-2H\left(\mathrm{zeta},\mathrm{zetab},r,u\right)\mathrm{du}&t\mathrm{du}\right)$
 ${\mathrm{_DG}}{}\left(\left[\left[{"tensor"}{,}{\mathrm{RT}}{,}\left[\left[{"cov_bas"}{,}{"cov_bas"}\right]{,}\left[{}\right]\right]\right]{,}\left[\left[\left[{1}{,}{1}\right]{,}{-}{2}{}{H}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{r}{,}{u}\right)\right]{,}\left[\left[{1}{,}{2}\right]{,}{-}{1}\right]{,}\left[\left[{2}{,}{1}\right]{,}{-}{1}\right]{,}\left[\left[{3}{,}{4}\right]{,}\frac{{{r}}^{{2}}}{{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}^{{2}}}\right]{,}\left[\left[{4}{,}{3}\right]{,}\frac{{{r}}^{{2}}}{{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}^{{2}}}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"tensor"}{,}{\mathrm{RT}}{,}\left[\left[{"cov_bas"}{,}{"cov_bas"}\right]{,}\left[{}\right]\right]\right]{,}\left[\left[\left[{1}{,}{1}\right]{,}{-}{2}{}{H}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{r}{,}{u}\right)\right]{,}\left[\left[{1}{,}{2}\right]{,}{-}{1}\right]{,}\left[\left[{2}{,}{1}\right]{,}{-}{1}\right]{,}\left[\left[{3}{,}{4}\right]{,}\frac{{{r}}^{{2}}}{{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}^{{2}}}\right]{,}\left[\left[{4}{,}{3}\right]{,}\frac{{{r}}^{{2}}}{{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}^{{2}}}\right]\right]\right]\right)$ (2.6)

Here is a null tetrad for this metric.

 RT > $\mathrm{NT}≔\mathrm{evalDG}\left(\left[\mathrm{D_r},\mathrm{D_u}-H\left(\mathrm{zeta},\mathrm{zetab},r,u\right)\mathrm{D_r},\frac{P\left(\mathrm{zeta},\mathrm{zetab},u\right)\mathrm{D_zeta}}{r},\frac{P\left(\mathrm{zeta},\mathrm{zetab},u\right)\mathrm{D_zetab}}{r}\right]\right)$
 $\left[{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{1}\right]{,}{1}\right]{,}\left[\left[{2}\right]{,}{-}{H}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{r}{,}{u}\right)\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{3}\right]{,}\frac{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}{{r}}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{4}\right]{,}\frac{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}{{r}}\right]\right]\right]\right)\right]{,}\left[{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{1}\right]{,}{1}\right]{,}\left[\left[{2}\right]{,}{-}{H}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{r}{,}{u}\right)\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{3}\right]{,}\frac{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}{{r}}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{4}\right]{,}\frac{{P}{}\left({\mathrm{ζ}}{,}{\mathrm{zetab}}{,}{u}\right)}{{r}}\right]\right]\right]\right)\right]$ (2.7)

The null congruence is very simple:

 RT > $U≔{\mathrm{NT}}_{1}$
 ${\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{\mathrm{RT}}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)$ (2.8)

First calling sequence:

 RT > $\mathrm{CongruenceProperties}\left(g,\mathrm{D_r}\right)$
 ${\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = 0, \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "Expansion" \right) = +\left(/\left(*\left(2\right), *\left(r\right)\right)\right) \right] \right)}}{,}{\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = 0, \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "Expansion" \right) = +\left(/\left(*\left(2\right), *\left(r\right)\right)\right) \right] \right)}}$ (2.9)

Third calling sequence:

 RT > $\mathrm{CongruenceProperties}\left(g,{\mathrm{NT}}_{1},{\mathrm{NT}}_{2}\right)$
 ${\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = 0, \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right), \left( "Expansion" \right) = +\left(/\left(*\left(2\right), *\left(r\right)\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}{,}{\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = 0, \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right), \left( "Expansion" \right) = +\left(/\left(*\left(2\right), *\left(r\right)\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}$ (2.10)

Fourth calling sequence

 RT > $\mathrm{CongruenceProperties}\left(g,\mathrm{NT}\right)$
 ${\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = 0, \left( "RotationNormSquared" \right) = 0, \left( "sigma" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right), \left( "rho" \right) = +\left(-\left(/\left(1, *\left(r\right)\right)\right)\right), \left( "Expansion" \right) = +\left(/\left(*\left(2\right), *\left(r\right)\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}{,}{\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = 0, \left( "RotationNormSquared" \right) = 0, \left( "sigma" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right), \left( "rho" \right) = +\left(-\left(/\left(1, *\left(r\right)\right)\right)\right), \left( "Expansion" \right) = +\left(/\left(*\left(2\right), *\left(r\right)\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", RT, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}$ (2.11)

Example 3.

Here is an example of a Newman-Tamburino metric of Petrov type I and which admits a null geodesic congruence with non-vanishing shear.

 RT > $\mathrm{DGsetup}\left(\left[u,r,x,y\right],M\right)$
 ${\mathrm{frame name: M}}$ (2.12)
 M > $g≔\mathrm{evalDG}\left({r}^{2}\mathrm{dx}&t\mathrm{dx}+{x}^{2}\mathrm{dy}&t\mathrm{dy}-\frac{2r\mathrm{du}&s\mathrm{dx}}{x}-2\mathrm{du}&s\mathrm{dr}+\frac{1\left(c+\mathrm{ln}\left({r}^{2}{x}^{4}\right)\right)\mathrm{du}&t\mathrm{du}}{{x}^{2}}\right)$
 ${\mathrm{_DG}}{}\left(\left[\left[{"tensor"}{,}{M}{,}\left[\left[{"cov_bas"}{,}{"cov_bas"}\right]{,}\left[{}\right]\right]\right]{,}\left[\left[\left[{1}{,}{1}\right]{,}\frac{{c}{+}{\mathrm{ln}}{}\left({{r}}^{{2}}{}{{x}}^{{4}}\right)}{{{x}}^{{2}}}\right]{,}\left[\left[{1}{,}{2}\right]{,}{-}{1}\right]{,}\left[\left[{1}{,}{3}\right]{,}{-}\frac{{r}}{{x}}\right]{,}\left[\left[{2}{,}{1}\right]{,}{-}{1}\right]{,}\left[\left[{3}{,}{1}\right]{,}{-}\frac{{r}}{{x}}\right]{,}\left[\left[{3}{,}{3}\right]{,}{{r}}^{{2}}\right]{,}\left[\left[{4}{,}{4}\right]{,}{{x}}^{{2}}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"tensor"}{,}{M}{,}\left[\left[{"cov_bas"}{,}{"cov_bas"}\right]{,}\left[{}\right]\right]\right]{,}\left[\left[\left[{1}{,}{1}\right]{,}\frac{{c}{+}{\mathrm{ln}}{}\left({{r}}^{{2}}{}{{x}}^{{4}}\right)}{{{x}}^{{2}}}\right]{,}\left[\left[{1}{,}{2}\right]{,}{-}{1}\right]{,}\left[\left[{1}{,}{3}\right]{,}{-}\frac{{r}}{{x}}\right]{,}\left[\left[{2}{,}{1}\right]{,}{-}{1}\right]{,}\left[\left[{3}{,}{1}\right]{,}{-}\frac{{r}}{{x}}\right]{,}\left[\left[{3}{,}{3}\right]{,}{{r}}^{{2}}\right]{,}\left[\left[{4}{,}{4}\right]{,}{{x}}^{{2}}\right]\right]\right]\right)$ (2.13)

Here is a null tetrad for this metric.

 M > $\mathrm{NT}≔\left[\mathrm{D_r},\mathrm{D_u}+\frac{\left(c+\mathrm{ln}\left({r}^{2}{x}^{4}\right)\right)\mathrm{D_r}}{2{x}^{2}},-\frac{\sqrt{2}\mathrm{D_r}}{x}+\frac{\sqrt{2}\mathrm{D_x}}{2r}+\frac{I\frac{1}{2}\sqrt{2}\mathrm{D_y}}{x},-\frac{\sqrt{2}\mathrm{D_r}}{x}+\frac{\sqrt{2}\mathrm{D_x}}{2r}-\frac{I\frac{1}{2}\sqrt{2}\mathrm{D_y}}{x}\right]$
 $\left[{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{1}\right]{,}{1}\right]\right]\right]\right){+}\frac{{1}}{{2}}{}\frac{\left({c}{+}{\mathrm{ln}}{}\left({{r}}^{{2}}{}{{x}}^{{4}}\right)\right){}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)}{{{x}}^{{2}}}{,}{-}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)}{{x}}{+}\frac{{1}}{{2}}{}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{3}\right]{,}{1}\right]\right]\right]\right)}{{r}}{+}\frac{\frac{{1}}{{2}}{}{I}{}\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{4}\right]{,}{1}\right]\right]\right]\right)}{{x}}{,}{-}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)}{{x}}{+}\frac{{1}}{{2}}{}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{3}\right]{,}{1}\right]\right]\right]\right)}{{r}}{-}\frac{\frac{{1}}{{2}}{}{I}{}\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{4}\right]{,}{1}\right]\right]\right]\right)}{{x}}\right]{,}\left[{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right){,}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{1}\right]{,}{1}\right]\right]\right]\right){+}\frac{{1}}{{2}}{}\frac{\left({c}{+}{\mathrm{ln}}{}\left({{r}}^{{2}}{}{{x}}^{{4}}\right)\right){}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)}{{{x}}^{{2}}}{,}{-}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)}{{x}}{+}\frac{{1}}{{2}}{}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{3}\right]{,}{1}\right]\right]\right]\right)}{{r}}{+}\frac{\frac{{1}}{{2}}{}{I}{}\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{4}\right]{,}{1}\right]\right]\right]\right)}{{x}}{,}{-}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)}{{x}}{+}\frac{{1}}{{2}}{}\frac{\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{3}\right]{,}{1}\right]\right]\right]\right)}{{r}}{-}\frac{\frac{{1}}{{2}}{}{I}{}\sqrt{{2}}{}{\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{4}\right]{,}{1}\right]\right]\right]\right)}{{x}}\right]$ (2.14)

Again we consider the first leg of this tetrad.

 M > $U≔\mathrm{D_r}$
 ${\mathrm{_DG}}{}\left(\left[\left[{"vector"}{,}{M}{,}\left[{}\right]\right]{,}\left[\left[\left[{2}\right]{,}{1}\right]\right]\right]\right)$ (2.15)

First calling sequence:

 RT > $\mathrm{CongruenceProperties}\left(g,U\right)$
 ${\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = +\left(/\left(*\left(/\left(1, 2\right)\right), *\left(^\left(r, 2\right)\right)\right)\right), \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "Expansion" \right) = /\left(1, *\left(r\right)\right) \right] \right)}}{,}{\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = +\left(/\left(*\left(/\left(1, 2\right)\right), *\left(^\left(r, 2\right)\right)\right)\right), \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "Expansion" \right) = /\left(1, *\left(r\right)\right) \right] \right)}}$ (2.16)

Third calling sequence:

 RT > $\mathrm{CongruenceProperties}\left(g,{\mathrm{NT}}_{1},{\mathrm{NT}}_{2}\right)$
 ${\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = +\left(/\left(*\left(/\left(1, 2\right)\right), *\left(^\left(r, 2\right)\right)\right)\right), \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[3, 3\right], +\left(*\left(/\left(1, 2\right), *\left(r\right)\right)\right)\right], \left[\left[4, 4\right], +\left(-\left(/\left(*\left(/\left(1, 2\right), *\left(^\left(x, 2\right)\right)\right), *\left(r\right)\right)\right)\right)\right]\right]\right]\right), \left( "Expansion" \right) = /\left(1, *\left(r\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}{,}{\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = +\left(/\left(*\left(/\left(1, 2\right)\right), *\left(^\left(r, 2\right)\right)\right)\right), \left( "RotationNormSquared" \right) = 0, \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[3, 3\right], +\left(*\left(/\left(1, 2\right), *\left(r\right)\right)\right)\right], \left[\left[4, 4\right], +\left(-\left(/\left(*\left(/\left(1, 2\right), *\left(^\left(x, 2\right)\right)\right), *\left(r\right)\right)\right)\right)\right]\right]\right]\right), \left( "Expansion" \right) = /\left(1, *\left(r\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}$ (2.17)

Fourth calling sequence:

 RT > $\mathrm{CongruenceProperties}\left(g,\mathrm{NT}\right)$
 ${\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = +\left(/\left(*\left(/\left(1, 2\right)\right), *\left(^\left(r, 2\right)\right)\right)\right), \left( "RotationNormSquared" \right) = 0, \left( "sigma" \right) = +\left(-\left(/\left(*\left(/\left(1, 2\right)\right), *\left(r\right)\right)\right)\right), \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[3, 3\right], +\left(*\left(/\left(1, 2\right), *\left(r\right)\right)\right)\right], \left[\left[4, 4\right], +\left(-\left(/\left(*\left(/\left(1, 2\right), *\left(^\left(x, 2\right)\right)\right), *\left(r\right)\right)\right)\right)\right]\right]\right]\right), \left( "rho" \right) = +\left(-\left(/\left(*\left(/\left(1, 2\right)\right), *\left(r\right)\right)\right)\right), \left( "Expansion" \right) = /\left(1, *\left(r\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}{,}{\mathrm{table\left( \left[\left( "ShearNormSquared" \right) = +\left(/\left(*\left(/\left(1, 2\right)\right), *\left(^\left(r, 2\right)\right)\right)\right), \left( "RotationNormSquared" \right) = 0, \left( "sigma" \right) = +\left(-\left(/\left(*\left(/\left(1, 2\right)\right), *\left(r\right)\right)\right)\right), \left( "Raychaudhuri" \right) = 0, \left( "RotationScalar" \right) = 0, \left( "ShearTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[3, 3\right], +\left(*\left(/\left(1, 2\right), *\left(r\right)\right)\right)\right], \left[\left[4, 4\right], +\left(-\left(/\left(*\left(/\left(1, 2\right), *\left(^\left(x, 2\right)\right)\right), *\left(r\right)\right)\right)\right)\right]\right]\right]\right), \left( "rho" \right) = +\left(-\left(/\left(*\left(/\left(1, 2\right)\right), *\left(r\right)\right)\right)\right), \left( "Expansion" \right) = /\left(1, *\left(r\right)\right), \left( "RotationTensor" \right) = _DG\left(\left[\left["tensor", M, \left[\left["cov_bas", "cov_bas"\right], \left[\right]\right]\right], \left[\left[\left[1, 1\right], 0\right]\right]\right]\right) \right] \right)}}$ (2.18)
 M >